Book Phys 15-05 Eleve
Book Phys 15-05 Eleve
Marc Beutier
obelix56@[Link]
[Link]
cpge TSI
Établissement St Joseph - LaSalle
Lorient
P0 Rappels de thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
0.1 Système thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
0.2 Équilibre thermodynamique et évolution quasi-statique . . . . . . . . . . . . . . . . . . . 29
0.2.1 Équilibre thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
0.2.2 Principe zéro . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
0.3 Grandeurs intensives et extensives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
0.4 Grandeur et fonction d’état . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
0.5 Premier principe de la thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . 30
0.5.1 Énergie totale d’un système . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
0.5.2 Travail en thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
0.5.3 Transfert thermique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
0.5.4 Types de transformations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
[Link] Transformation quelconque . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
[Link] Transformation quasi-statique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
0.5.5 Transformation mécaniquement réversible . . . . . . . . . . . . . . . . . . . . . . . . 32
0.5.6 Transformation réversible . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
0.5.7 Récapitulatif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
0.6 Différents types de travaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
0.6.1 Évolution adiabatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
0.6.2 Évolution isochore . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
0.6.3 Évolution isoénergétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
0.6.4 Évolution monobare . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
0.6.5 Évolution isobare . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
0.6.6 Évolution quasi-statique et isotherme . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
0.6.7 Évolution polytropique d’indice k . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
0.7 Enthalpie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
0.8 Gaz parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
0.8.1 Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
0.8.2 Équation d’état des gaz parfaits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
0.8.3 Lois de Joule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
[Link] Première loi de Joule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
[Link] Deuxième loi de Joule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
0.8.4 Rapport ou cœfficient γ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
0.8.5 Relation de Mayer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
0.9 Phases condensées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
3
P2 Principes de thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.1 Sous-systèmes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.2 Premier principe de la thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.2.1 Énergie totale d’un système . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.2.2 Premier principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
2.2.3 Capacités thermiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
[Link] Capacité thermique à pression constante . . . . . . . . . . . . . . . . . . . . . . . . 69
[Link] Capacité thermique à volume constant . . . . . . . . . . . . . . . . . . . . . . . . . 69
2.3 gaz parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
2.3.1 Lois de Joule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
[Link] Première loi de Joule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
[Link] Deuxième loi de Joule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
2.3.2 Rapport ou cœfficient γ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
2.3.3 Relation de Mayer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
2.4 Phases condensées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.5 Deuxième principe de la thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.5.1 Entropie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.5.2 Énoncés du deuxième principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.6 Identités thermodynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
2.6.1 Première identité thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
2.6.2 Deuxième identité thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
2.6.3 Définition de la température . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
2.7 Entropie d’un corps pur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
2.7.1 Variation de l’entropie pour un Gaz Parfait (G.P.) . . . . . . . . . . . . . . . . . . . . 73
2.7.2 Rappels : Loi de Laplace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
2.7.3 Variation de l’entropie pour une phase condensée . . . . . . . . . . . . . . . . . . . . . 75
2.8 Machines thermiques : rappels et compléments . . . . . . . . . . . . . . . . . . . . . . . 75
2.8.1 Moteurs et machines frigorifiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
2.8.2 Coefficient d’efficacité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
[Link] Cycle moteur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
[Link] Cycle frigorifique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
[Link] Cycle de pompe à chaleur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
[Link] Les différents cycles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
2.8.3 Cycle de Carnot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
[Link] Dans le diagramme de Watt (p, V ) . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
[Link] Dans le diagramme entropique (T, s) . . . . . . . . . . . . . . . . . . . . . . . . . . 80
2.9 Exercices : principes de la thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . 81
2.9.1 Énoncés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
Électricité 161
P8 Montages de base avec l’ALI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167
8.1 Montage suiveur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167
8.2 Montage amplificateur inverseur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
8.3 Montage amplificateur non inverseur . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
8.4 Montage sommateur inverseur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 169
8.5 Montage dérivateur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 171
8.6 Montage intégrateur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172
8.7 Montage pseudo-intégrateur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Électromagnétisme 238
Optique 353
Chimie 414
P4 Calculs de pH . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 457
4.1 Eau pure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 457
4.2 Acide fort dans l’eau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 458
4.3 Base forte dans l’eau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 459
4.4 Acide faible dans l’eau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 462
4.5 Base faible dans l’eau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 464
Annexes 497
P1 Analyse Vectorielle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 503
1.1 Produit scalaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 503
1.2 Produit vectoriel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 503
1.3 Dérivation d’un vecteur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 503
1.4 Les différents systèmes de coordonnées . . . . . . . . . . . . . . . . . . . . . . . . . . . 504
1.4.1 Coordonnées cartésiennes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 504
[Link] Définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 504
[Link] Schéma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 504
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
PIndex général, des noms propres, des symboles, des unités . . . . . . . . . . . . . . . . . . . 580
P2 Principes de thermodynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.1 Sous-systèmes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
[Link] Enthalpie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
[Link] Entropie de changement d’état . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
3.1.8 Équilibre liquide-vapeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
[Link] Règle des moments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
[Link] Variations élémentaires d’enthalpie et d’entropie de changement d’état . . . . . . . . 89
3.1.9 Autres propriétés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
3.2 Tables et diagrammes thermodynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . 90
3.2.1 Tables thermodynamiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
3.2.2 Diagrammes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
[Link] Diagramme (p, h) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
[Link] Diagramme (h, s) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
3.3 Allure des isoparamétriques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
3.3.1 Diagramme (p, V ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
[Link] Dans le domaine vapeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
[Link] Dans le domaine liquide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
[Link] Dans le domaine liquide-vapeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
3.3.2 Diagramme (T, s) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
[Link] Dans le domaine vapeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
[Link] Dans le domaine liquide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
[Link] Dans le domaine liquide-vapeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
3.4 Exercices : fluides réels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
3.4.1 Énoncés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
P0
— Système thermodynamique —
Un système, thermodynamique ou non, est une surface fermée qui délimite un milieu intérieur d’un
milieu extérieur. Il en existe trois types :
• Système ouvert : Les échanges de matière et d’énergie sont possibles.
• Système fermé : L’échange d’énergie est possible uniquement.
• Système isolé : Aucun échange n’est possible avec le milieu extérieur.
— É quilibre thermodynamique —
Dans ce cas, les grandeurs moyennées telles que T , p, U sont définies, ce qui n’est pas le cas hors équilibre.
— Prin ipe zé ro —
Ce principe est celui de l’existence d’état d’équilibre. Dans ces états d’équilibre, les grandeurs thermody-
namiques, grandeurs macroscopiques, sont parfaitement définies.
Une grandeur intensive est indépendante de la taille du système. Une grandeur extensive dépend de
la taille du système. Elle peut être conservative ou non conservative.
Autrement dit, une grandeur est extensive si l’addition a un sens pour elle. Si ce n’est pas le cas, la grandeur
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Une grandeur est dite d’état si elle caractérise un système. Elle est définie par l’expérience. Pour
un gaz, on a par exemple : p, V , T .
Une fonction d’état est un lien entre les grandeurs d’état et sa variation entre 2 états ne dépend pas
du chemin suivi.
Le premier principe est un principe de conservation. Il existe une grandeur, appelé énergie totale,
extensive et conservative, que l’on peut définir dans tout système fermé. On appelle énergie totale
d’un système, toute l’énergie présente, peu importe sa forme.
Etotale = Ec + U
∆Etotale = W + Q
où W est le travail reçu par le système, et Q la quantité de chaleur (ou transfert thermique) reçue.
On obtient alors :
∆U + ∆Ec = W + Q
— Remarque —
Dans le cas de la détente de Joule Gay-Lussac, on fait tendre W et Q vers 0, ce qui en fait une
transformation à énergie interne constante.
Les unité s :
δ 3 W = −pext(M ) d3 V
ou plus simplement :
δW = −pext dV
On en déduit donc que, si le système se dilate (dV > 0), on a δW < 0 : le milieu extérieur est
résistant. Si le système se comprime (dV < 0), on a δW > 0 : le milieu extérieur est moteur. Si
dV = 0, alors δW = 0 et l’évolution est isochore.
Les unité s :
Le transfert thermique, également appelé quantité de chaleur, peut être positif, négatif ou nul (δQ > 0,
δQ < 0, δQ = 0).
— Proprié té —
Pour ce dernier cas (δQ = 0), on dit que l’évolution est adiabatique.
— Remarque —
Lors d’une transformation quelconque, il peut ne pas y avoir d’équilibres intermédiaires, et on ne pourra
alors pas définir les variables d’état (U , p, T , . . .). On ne peut pas écrire d’équations pour décrire l’évolution
du système, et on ne peut donc pas l’étudier. On a besoin d’hypothèses supplémentaires . . .
— É volution quasi-statique —
Une transformation ou évolution est quasi-statique ou infiniment lente si elle amène le système de
l’état initial à l’état final par une succession d’états d’équilibre interne très proches.
La constante de temps de l’évolution du système doit être supérieure aux constantes de temps de relaxation
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Conséquence :
p et T peuvent être définies à chaque instant dans le système (mais elles peuvent différer des pression et
température extérieures pext et Text )
Une transformation est réversible s’il n’y a pas de phénomènes non conservatifs (frottements, . . .), ce qui
est quasiment impossible avec un système thermodynamique réel. De manière idéale, on doit supposer que
l’évolution est quasi-statique, et que le système est à tout instant à l’équilibre avec l’extérieur. L’équilibre
mécanique est en général le plus rapide à s’établir.
Conséquence :
p et T peuvent être définies à chaque instant dans le système ET la pression interne p est égale à la pression
externe pext (Mais T peut être différente de Text )
On suppose que la transformation est quasi-statique et que les deux équilibres mécanique et thermique sont
atteints à chaque instant.
p et T sont définies et égales à Pext et Text à tout instant.
0.5.7 Récapitulatif
En résumé :
Type de transformation T p
quelconque T non définie p non définie
quasi-statique T définie : peut être 6= Text p définie : peut être 6= pext
mécaniquement réversible T définie : peut être 6= Text p définie : à tout instant, p = pext
(équilibre mécanique à tout instant)
réversible T définie : à tout instant, T = Text p définie : à tout instant, p = pext
(sous entendu mécanique et thermique) (équilibre thermique à tout instant) (équilibre mécanique à tout instant)
∆U = W + Q
Une évolution isochore est caractérisée par un volume constant, donc par un travail nul (d3 V = 0 ou
dV = 0). Dans le cas d’une évolution isochore, on a donc :
W =0
et on obtient alors :
∆U = QV = Q
W = −Q
Une évolution est monobare si la pression extérieure est constante. On la note p2 car on sait que dans l’état
final, l’équilibre mécanique est réalisé et que la pression du système est alors égale à la pression extérieure.
On obtient alors l’expression du travail en fonction des volumes initial V1 et final V2 :
W = −p2 (V2 − V1 )
L’évolution quasi-statique de n moles de gaz parfait (G.P.) est caractérisée par une pression inté-
rieure qui est égale à tout moment à la pression extérieure. On a donc : p(t) = pext = p1 = p2
L’expression du travail est :
W = −p1 (V2 − V1 )
Pour un gaz parfait qui vérifie l’équation d’état des gaz parfaits, on obtient :
Z Z
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
V2
dV
W =− pdV = −nRT1
V1 V
Soit :
Å ã
V1
W = nRT1 ln
V2
On dit d’une évolution qu’elle est polytropique d’indice k si elle est quasi-statique et que :
Soit :
Ç å Ç å
V21−k − V11−k cte × V21−k − cte × V11−k
W = −cte =−
1−k 1−k
Ç å
p2 V2k × V21−k − p1 V1k × V11−k p2 V2 − p1 V1
W = =
k−1 k−1
nR(T2 − T1 )
W =
k−1
— Remarque —
0.7 Enthalpie
H = U + pV
— Transformation monobare —
Cette fonction enthalpie est notamment utilisée dans une transformation monobare, c’est-à-dire à
pression extérieure constante. On montre que dans le cas d’une transformation monobare (pression
extérieure constante), on a :
Qp = ∆H
0.8.1 Définition
Sur les plans macroscopiques, on appelle gaz parfait,noté G.P., tout gaz vérifiant simultanément :
• La loi de Boyle-Mariotte :
— Loi de Boyle-Mariotte —
• La loi d’Avogadro :
— Loi d'Avogadro —
Tous les gaz ont le même volume molaire dans les mêmes conditions de pression et de tem-
pérature.
• La loi de Charles :
— Loi de Charles —
• La loi de Gay-Lussac :
— Loi de Gay-Lussa —
• La loi de Dalton :
— Loi de Dalton —
La pression totale exercée par un mélange de gaz est égal à la somme des pressions que
chaque gaz exercerait s’il était seul.
pV = nRT
Avec :
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Les unité s :
Cette équation est éloignée des résultats expérimentaux. On utilise parfois un autre modèle que celui du gaz
parfait, celui de Van der Waals.
Le première loi de Joule postule que l’énergie interne U d’un gaz parfait est une fonction qui
dépend uniquement de la température : U = f1 (T ).
∆U = CV ∆T = nCV,m ∆T = mcV ∆T
Avec Cv la capacité calorifique à volume constant (CV,m et cV sont les capacités calorifiques respectivement
molaires et massiques à volume constant).
On utilise donc U dans un système isochore.
De même, la deuxième loi de Joule postule que l’enthalpie interne H d’un gaz parfait est également
une fonction qui dépend uniquement de la température : H = f2 (T ).
Avec Cp la capacité calorifique à pression constante. (Cp,m et cp sont les capacités calorifiques molaires et
massiques à pression constante).
On utilise donc H dans une transformation monobare.
Cp,m Cp cp
On définit la grandeur γ = = = >1
CV,m CV cV
— Relation de Mayer —
Cp − CV = nR
donc
Cp,m − CV,m = R
ou encore
n R
cp − cV = R= =r
m M
— Phases ondensé es —
Pour les phases condensées (solides ou liquides), on peut faire l’hypothèse que U (T ) et H(T ) sont
équivalentes.
On obtient alors :
Cp = CV = C
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Et :
0.11.1 Entropie
— Entropie —
Pour tout système fermé, on peut définir une fonction d’état, notée S, appelée entropie, qui est
une grandeur extensive, mais non conservative. Cette grandeur peut être créée ou non, mais jamais
détruite.
— Bilan entropique —
Les unité s :
— Remarque —
— Rappel —
Une transformation (ou évolution) est quasi-statique ou infiniment lente si elle amène le système
de l’état initial à l’état final par une succession d’états d’équilibre très proches.
Une transformation ou évolution est renversable si une modification infinitésimale des conditions
permet à un système qui a évolué sous l’influence de ce phénomène de retrouver son état immé-
diatement antérieur. Cette transformation ou évolution devient réversible si, en plus de la condition
nécessaire précédente, il en est de même du milieu extérieur.
— É volution ré versible —
Cela signifie que l’entropie créée lors de cette transformation est nulle : Scréée = 0. Le bilan entropique se
traduit alors par :
Z Z
δQ δQ
∆S = =
Text T
Le troisième principe de la thermodynamique a été formulé par Nernst en 1906 et s’énonce de la façon
suivante :
L’entropie d’un corps pur en équilibre thermodynamique tend vers zéro lorsque la température tend
vers zéro.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Les unité s :
Tf Vf
∆S = nCV,m ln + nR ln
Ti Vi
ou encore :
Tf pf
∆S = nCp,m ln − nR ln
Ti pi
Tf Tf Tf
∆S = C ln = nCm ln = mc ln
Ti Ti Ti
— Transformation isentropique —
∆S = 0
— loi de Lapla e —
pV γ = cte
Cp Cp,m cp
γ= = =
CV CV,m cV
T γ p1−γ = cte′
ou bien :
T V γ−1 = cte′′
— Remarque —
Dans le cadre d’une évolution isentropique, la pente de cette évolution dans un diagramme de Watt
est γ fois celle d’une isotherme (ceci est démonté page 78).
— É volution y lique —
On dit d’un système qu’il subit une évolution cyclique si, périodiquement, le système passe par un
même état thermodynamique. Le bilan de toute fonction d’état, est donc, sur un cycle, nul.
Une machine ne peut pas effectuer un travail si elle n’est en contact qu’avec une seule source de
chaleur.
Autrement dit, dans un système cyclique, il est impossible de convertir totalement en travail de la chaleur
extraite d’une seule source à température constante.
— É non é de Clausius —
Il ne peut pas y avoir de transformation dont le seul effet serait de transporter de la chaleur d’une
source froide vers une source chaude.
Autrement dit, dans un système cyclique, il est impossible de transférer de la chaleur d’une source froide à
une source chaude sans un apport extérieur.
• isobare,
• isochore,
• isotherme,
• isentropique.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Dans un diagramme de
I Watt (pV ), un cycle est dit moteur s’il décrit une évolution dans le sens horaire.
Dans ce diagramme, pdV est égale à l’aire du cycle. Il en découle que | W | est l’aire du cycle dans ce
diagramme.
Pour un tel cycle :
∆Ucycle = W + Q = 0
Dans un diagramme de Clapeyron (pv), un cycle est dit moteur s’il décrit également une évolution dans le
sens horaire. I
|W |
Dans ce diagramme, pdv est égale à l’aire du cycle. Il en découle que | w |= est l’aire du cycle
m
dans ce diagramme.
Pour un cycle :
∆ucycle = w + q = 0
• les thermostats ou souces idéales : leur température est une constante (océan, lac, . . .),
• les pseudo-sources ou sources réelles : leur température peut évoluer sous l’influence de la machine
thermique (cuve, . . .).
Soit Qf le transfert thermique de la source froide vers le système, en contact avec un thermostat de tempé-
rature Tf et Qc le transfert thermique de la source chaude vers le système en contact avec un thermostat de
température Tc .
◦ W >0
Attention aux conventions : Si par exemple, la source chaude est constituée d’eau de capacité
calorifique massique ce dont la température est égale à Te et variant de Te,I à Te,F , on a alors la
relation suivante :
Z Te,F
car Qc est la chaleur reçue par le système de la part de la source chaude. C’est donc la chaleur cédée
par l’eau chaude, soit l’opposé de la chaleur reçue par celle-ci !
Dans cette étude, sachant que le système est en contact avec deux thermostats, on a les formules suivantes :
Premier principe :
∆Ucycle = W + Qf + Qc = 0
Deuxième principe :
Qc Qf
∆Scycle = + + Scréation = 0
Tc Tf
Comme l’entropie créée Scréation est positive ou nulle, on en déduit l’inégalité de Clausius :
Qc Qf
+ 60
Tc Tf
— Remarques —
Comme précédemment, l’entropie créée Scréation est positive ou nulle et on en déduit l’inégalité de
Clausius :
Z Z
δQc δQf
+ 60
Tc Tf
— Effi a ité —
Énergie utile
η=
Énergie dépensée
Les unité s :
• η=− .
Qc
Dans une machine frigorifique, un fluide caloporteur (CF C, N H3 , ...) circule. Sous l’action du compres-
seur, le fluide caloporteur passe de sa phase gazeuse à sa phase liquide. Ce faisant, il restitue de l’énergie
à la source chaude. Sous l’action du détendeur, le fluide caloporteur passe de sa phase liquide à sa phase
gazeuse. Ce faisant, il absorbe de l’énergie au contact de la source froide.
La machine frigorifique fonctionne donc à l’aide des transitions de phases, qui réclament des transferts
thermiques.
Le fonctionnement en machine frigorifique est caractérisé par :
• W > 0,
• Qc < 0,
• Qf > 0,
Qc
• η=− .
W
La pente de la tangente à la courbe représentative d’une adiabatique réversible (donc isentropique) est plus
grande que celle représentant une isotherme.
Les isothermes du gaz parfait sont des hyperboles, l’isotherme Tc étant située au-dessus de l’isotherme Tf .
Les deux isothermes sont reliées par deux adiabatiques, de pentes plus importantes, d’où l’allure du cycle :
p(P a)
2
P2
2′
1′ Tc
P1 1
Tf
V (m3 )
Le rendement de ce moteur correspond au rapport du travail fourni sur un cycle par le moteur à la chaleur
qu’il absorbe depuis la source chaude pour pouvoir fonctionner :
W Wfourni
η=− =
Qc Qc
∆Ucycle = 0 = Qc + Qf + W
Qc Qf
∆Scycle = 0 = +
Tc Tf
On en déduit :
Qf Tf
η =1+ =1−
Qc Tc
Ce rendement est indépendant de la nature du fluide, puisque seules interviennent les températures des
sources. Pour un moteur fonctionnant entre les deux mêmes sources, le second principe s’écrit, dans le cas
général, avec Scréée :
Qc Qf
∆Scycle = 0 = + + Scréée
Tc Tf
Soit :
Qc Qf
+ = −Scréée 6 0
Tc Tf
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
On en déduit :
Qf Tf
6−
Qc Tc
et donc :
W Qf Tf
η=− =1− 61− = ηmax
Qc Qc Tc
Le rendement du cycle de Carnot est maximal pour des sources dont la température est fixée.
— Diagramme entropique —
Dans le diagramme entropique, les isothermes 1 → 1′ et 2 → 2′ sont horizontales, tandis que les transfor-
mations adiabatiques réversibles, donc isentropiques 1′ → 2 et 2′ → 1 sont des verticales :
T (K)
Tc
2 2′
Tf
1′ 1
s( J.K −1 .kg −1 )
smin smax
De même, comme le fluide cède réversiblement de la chaleur à la source froide, son entropie décroît et :
Qf
∆S2→2′ = Smin − Smax =
Tf
Le cycle de Carnot traduisant le meilleur fonctionnement possible, il est par conséquent évident que
les coefficients d’efficacité réels sont en réalité inférieurs à ceux qui sont donnés sur la dernière ligne
du tableau.
Ce cycle est constitué de deux isochores (BC et DA) et de deux isentropiques (AB et CD) La formule de
Séchangée développée dans l’étude n’est plus valable ici, les sources thermiques n’étant plus des thermostats,
mais des pseudo-sources.
On assimile un système à un corps pur quand on peut considérer que ce système est constitué
d’entités chimiques de même nature. Le corps pur est un modèle, dans lequel tout corps étranger est
supposé infiniment dilué.
On appelle corps pur simple un corps composé d’une unique entité comme H2 ou O2 .
— Système diphasé —
On dit qu’un système est diphasé s’il contient un corps pur sous deux phases (eau liquide et gaz
par exemple).
[Link] Transition
En voici un résumé :
Gaz
Co blim
risa tion
n
Su
nd atio
tio
c
Va uéfa
en s
atio n
Liq
po
n
Liquide Solide
Solidification
Fusion
S L
S/L Fluide
C
b
hypercritique
L/G
S/G
T
G
p(P a)
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
LIQUIDE
N M
p2
p1
SOLIDE
T
VAPEUR
T (K)
0 T1 T2
F IGURE 5 – Diagramme p = f (T ) : cas général
p(P a)
LIQUIDE
N M
p2
p1
SOLIDE
T
VAPEUR
T (K)
0 T1 T2
F IGURE 6 – Diagramme p = f (T ) : cas de l’eau
La pente des courbes est toujours croissante, c’est à dire qu’une augmentation de T impose une augmenta-
tion de pression dans un corps pur.
Il existe une exception : pour l’eau, la pente de la courbe de fusion (frontière Solide → Liquide) est négative.
Il existe deux points caractéristiques sur ce diagramme :
— Point triple —
T , appelé point triple, correspond au domaine de coexistence des trois phases du corps pur consi-
déré. Ce point est une propriété intrinsèque du corps.
— Point ritique —
Le point critique C est la limite de coexistence des phases liquide et gazeuse d’un corps pur. Au
dessus de ce point, on passe de manière continue de l’état gazeux à l’état liquide sans observer
d’état diphasé. On observe ce qu’on appelle un état fluide.
p(P a)
N M VAPEUR
ps (T2 )
LIQUIDE L LIQ ⇄ VAP V
ps (T1 ) isotherme T2
isotherme T1
[Link] Commentaires
Parcourons une isotherme (appelée isotherme d’Andrews) placée sous le point critique en augmentant le
volume :
La transformation qui décrit la courbe LV s’effectue à pression constante. Cette pression, qui est notée
ps (T ), est la pression de vapeur saturante. Elle augmente avec la température.
• Entre ces deux points, le système est un composé d’un mélange liquide - gaz. Le gaz est appelé "va-
peur saturante".
• On considère, en général, que la vapeur sèche ainsi que la vapeur saturante peuvent être assimilées à
des gaz parfaits (ce n’est pas le cas lorsque la pression est trop élevée).
Fluide supercritique : Lorsque T > TC et que p > pC , on dit alors que le système est à l’état de fluide
supercritique, on passe alors de l’état liquide à l’état gazeux de manière continue.
Opalescence optique : Au point critique, pour T = TC , la densité du milieu connaît de très fortes
fluctuations. On montre que l’indice d’un milieu est lié à la densité de celui-ci. Dans ces conditions, le
corps pur diffuse de manière particulière la lumière, on appelle ce phénomène opalescence optique.
Stockage des fluides D’après l’équation d’état de l’eau liquide, on à montré qu’une faible augmentation
de température, lors d’une évolution isochore, entrainait une très forte augmentation de pression. Pour cette
raison, on stocke les fluides dans un état diphasé. On les stocke dans un volume V supérieur au volume
VC du point critique, car une augmentation de température entraîne alors une augmentation de pression
relativement faible.
— Proprié té —
∆H = mℓ = L = Qp
On appelle chaleur latente (ou enthalpie) de changement d’état d’une masse ou quantité de matière
la quantité d’énergie qu’il faut lui fournir à pression constante pour qu’elle passe de l’état initial
(solide, liquide ou gazeux) à un autre état.
La variation d’enthalpie d’un corps pur au cours d’une transition de phase est également appelée chaleur
latente de transition de phase, notée L.
De plus, nous savons que Q est positif pour une transition de phase d’un état ordonné vers un état moins
ordonné. Ceci implique que dans ce cas, L est positif. Dans la situation inverse, L est négatif, tout comme
Q.
Dans le cadre de cette évolution isotherme et isobare, la variation d’entropie d’une masse m de corps pur
est donnée par :
mℓ
∆S =
T
[Link] Hygrométrie
Considérons de l’eau. Cette eau se vaporise, à une température T0 , jusqu’à atteindre la pression de vapeur
saturante, ps , qui dépend de T . Dans l’atmosphère, l’eau se vaporise jusqu’à que la pression partielle de
l’eau soit égale à ps . On définit donc τ , le taux d’hygrométrie, qui est le rapport de la pression partielle de
l’eau sur la pression de vapeur saturante :
pH2 O
τ=
ps (T0 )
— É tat mé tastable —
Très souvent, dans les transitions de phases, il existe un retard. On parle durant ce retard, d’état
métastable.
Dans ces deux derniers cas, il y retard à la solidification appelé phénomène de surfusion.
P1
Éléments de statique des fluides
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Dans ce chapitre, on étudie les propriétés d’un fluide au repos dans un référentiel galiléen en tenant compte
du champ de pesanteur.
— Fluide —
Outre le fait qu’un fluide est déformable, il est également sans rigidité : il peut s’écouler et épouser la forme
du récipient qui le contient.
On considère dans cette définition les liquides et les gaz.
Les liquides possèdent un volume propre : ils sont peu compressibles et les actions à courte distance (les
forces de Van Der Waals) interviennent de façon importante dans les propriétés des liquides.
Les gaz occupent tout le volume qui leur sont offert et ils sont compressibles.
Comme nous venons de le voir, les liquides peuvent être, avec une bonne précision, considérés
comme incompressibles, à l’inverse des gaz.
Cependant, la masse volumique des gaz peut dans certains cas être considérée comme uniforme et
dans ce cas, ces gaz pourront être assimilés à des fluides incompressibles : cf. page 99.
⇒ Activité 1. 1
Données :
Masse molaire de l’eau : M = 18 [Link]−1,
Masse volumique de l’eau : ρ = 1,00 kg.L−1,
Nombre d’Avogadro : N = 6,02 1023 mol−1 ,
Volume molaire de l’air à 25o C : Vm = 25 [Link]−1.
L’approximation des milieux continus consiste alors à définir des champs macroscopiques en faisant
des moyennes sur les éléments de volume d3 V(M) .
tube à obturateur
— Pression —
On appelle pression p(M) exercée par le fluide en M , la grandeur scalaire définie par :
−
→
d2 Fp = p(M) d2 S(M) −
→
n
−
→
d2 Fp est appelée force pressante élémentaire.
Les unité s :
Pour mesurer une pression, on peut se servir d’une capsule manométrique, raccordée à un manomètre en
U:
h
tube de raccordement
avec −
n→
b
b
b
b
b b
−−→
On peut également poser d2 S (M) = d2 S(M) −
n→ 2 −→
ext = d S(M) next .
Considérons un fluide quelconque. Soit M un point quelconque dans le fluide. Soit ρ(M) la masse volu-
mique dans le fluide, au voisinage de M , d3 m(M) la masse d’un volume élémentaire d3 V(M) au voisinage
de M . On a :
Parmi les actions extérieures subies par une partie quelconque d’un fluide, il faut distinguer :
• les forces volumiques : la force dF~v subie par une particule de fluide est proportionnelle à dV(M) .
Cela permet de définir une densité volumique de force :
dF~v
f~v =
dV(M)
• les forces surfaciques : la force dF~s subie par une particule de fluide est proportionnelle à dS. Cela
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
dF~s
f~s =
dS
Ainsi, pour les forces de pression, on obtient :
dF~p = −p(M ) dS ~n
fluide
Soit une particule de fluide de volume dV(M) en équilibre, dé-
coupé en coordonnées cartésiennes (dV(M) = dx dy dz), sou-
mise au champ de pesanteur et aux forces de pression. On a alors dz
−
→
g
à l’équilibre la relation :
dF~p + dF~v = ~0
dx
Avec un axe z vertical ascendant, on a ~g = −g →
−
ez et :
−p(x) dydz(− −
→
ex ) − p(x + dx) dy dz(+ −
→
ex ) − (p(x + dx) − p(x)) dy dz −
→
ex 0−→
ex
−
→ →
− →
− →
− →
−
dFp = −p(y) dxdz(− ey ) − p(y + dy) dx dz(+ ey ) = − (p(y + dy) − p(y)) dx dz ey = 0 ey
−p(z) dxdy(− −
→
ez ) − p(z + dz) dx dy(+ −
→
ez ) − (p(z + dz) − p(z)) dx dy −
→
ez ρ g dV(M) −
→
ez
¶
− (p(x + dx) − p(x)) dy dz = 0 − (p(y + dy) − p(y)) dx dz = 0 − (p(z + dz) − p(z)) dx dy = ρ g dx dy dz
La première équation indique que p(x + dx) = p(x), donc que la pression est indépendante de la variable
d’espace x.
Même raisonnement pour p(y) avec la deuxième équation : la pression est donc également indépendante de
la variable y.
— Remarque —
Les 2 résultats précédents pouvaient être obtenus directement en remarquant que les forces de pres-
sion sur les surfaces latérales s’annulaient 2 à 2.
On en déduit :
dp
= −ρ g
dz
C’est sous cette forme que nous retiendrons ce principe fondamental de la statique des fluides.
dp
= +ρ g
dz
— Fluide in ompressible —
ρM = cte = ρ0
p(M) = p0 − ρ0 gz(M)
avec ici, z axe vertical ascendant.
Toute variation de pression en un point d’un liquide au repos se transmet intégralement à tous les
autres points du liquide.
• Principe des vases communicants : la surface libre d’un fluide est contenue dans un plan horizontal.
• L’interface entre deux fluides de densités différentes est contenue dans un plan horizontal.
• Les isobares sont perpendiculaires à l’accélération de la pesanteur ~g .
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
p∗ = p + ρgz
⇒ Activité 1. 2
∆p
Calculer la variation relative de pression entre les altitudes 0 et 100 m.
p0
L’atmosphère terrestre s’étend sur plusieurs dizaines de kilomètres et apparaît donc en revanche comme un
système suffisamment étendu pour que l’influence de la pesanteur s’y fasse sentir.
Néanmoins, le modèle de l’atmosphère isotherme ne s’applique qu’à la haute atmosphère, pour des couches
d’air dont l’altitude est comprise entre 11 et 30 km avec une température de l’ordre de 223 K. En effet,
l’uniformité de la température suppose un brassage suffisant des couches atmosphériques ce qui n’est pas
le cas à basse altitude où le modèle du gradient de température est mieux adapté.
Considérons un gaz parfait isotherme de température T0 .
pV = nRT
Les unité s :
m Mp
ρ= =
En considérant un axe z vertical ascendant, la loi fondamentale de la statique des fluides s’écrit :
dp Mg
= −ρg = − p
dz RT0
dp Mg
=− dz
p RT0
On obtient alors :
M gz
−
p(z) = p0 e RT0
RT
En posant H = , on peut aussi exprimer la pression sous la forme suivante :
Mg
z
−
p(z) = p0 e H
La masse volumique ρm de l’eau de mer est environ égale à 1,025 kg.L−1, donc un peu plus grande que
l’eau douce.
Si on descend dans l’eau de mer d’une profondeur ∆z, la pression varie de ∆p = ρm g∆z.
Ainsi, pour une variation de profondeur de 10 m, la pression varie d’environ 103 × 10 × 10 = 105 P a, soit
1 bar.
La pression à la surface de la mer étant d’environ égale à 1 bar, elle vaut par exemple 11 bar à 100 m.
— Remarques —
• Le calcul exact en prenant g = 9,81 m.s−2 et ρm = 1,025 kg.L−1 modifie très peu le résultat
précédent.
• La masse volumique de l’eau de mer varie en fonction de sa salinité. Cette dernière varie à la
surface du globe et en un endroit donné, elle varie en fonction de la profondeur.
• L’accélération de la pesanteur g varie également en fonction de la profondeur . . .
Tous corps immergé dans un ou plusieurs fluides subit une force pressante résultante opposée au
poids et égale en module au poids du fluide déplacé.
→
−
Cette force est notée Π et est appelée poussée d’Archimède. Elle est appliquée au centre d’inertie
du fluide déplacé, point appelé centre de poussée. On le note souvent C.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Les unité s :
Le modèle de l’atmosphère isotherme n’étant pas très réaliste, on peut admettre que jusqu’à une altitude de
l’ordre de 10 km, la température vérifie la loi T (z) = T0 (1 − az). C’est ce qu’on appelle un gradient de
température.
dp Mp
La relation = −ρg avec ρ = donne :
dz RT
dp pM g M gdz
=− dz = −
dz RT pRT0 (1 − az)
Mg
En posant β = , et en séparant les variables, on obtient l’expression de la pression p :
aRT0
β
p(z) = p0 (1 − az)
Un autre modèle est encore possible, c’est celui de l’atmosphère isentropique (avec pV γ = cte).
p p p0
On a alors γ = cte, d’où γ = γ et :
ρ ρ ρ0
Å ã1/γ
p
ρ = ρ0
p0
— Remarque —
ï ò γ
γ − 1 Mg γ − 1
p(z) = p0 1 −
γ RT0
De plus, on a :
γ − 1 Mg
T (z) = T0 − z
γ R
Considérons une particule M de masse dm(M) au repos dans un référentiel non galiléen. On ajoute deux
forces dans la somme des forces du principe fondamental de la dynamique (P.F.D.) :
−→
• la force d’inertie d’entraînement d3 Fie = −d3 m(M) −
→
ae ,
−
→ →
− −
→
• la force d’inertie de Coriolis d Fic = −d m(M) ac = 0 car −
3 3 →
− →
vr = 0 .
• rotation uniforme autour d’un axe fixe d’un référentiel galiléen, alors :
a = −ω 2 −
→
− →r
e
avec k →
−
r k la distance entre l’axe de rotation et le point M.
1.10.1 Énoncés
Dans tous les exercices, sauf précision contraire, la pression atmosphérique sera prise égale à p0 =
patm = 105 P a, l’accélération de la pesanteur à g = 10 m.s−2 et la masse volumique de l’eau à
ρeau = 1,00 kg.L−1.
◮ Exercice 1. 3 : Tube en U
Pour mesurer la pression dans un gaz, on utilise le tube en U ci-dessous (cf. figure 1.5), rempli de liquide,
de masse volumique ρℓ .
Lorsque la pression dans le gaz est différente de la pression atmosphérique, la dénivellation dans le liquide
est différente de 0.
On considère un barrage de largeur ℓ, dont la surface de contact avec l’eau est supposée plane. La hauteur
d’eau est h.
z
1. Déterminer la résultante des forces exercées par
l’eau sur le barrage. h
2. Déterminer le moment résultant de ces forces
par rapport à la base du barrage.
3. En déduire la position du centre de poussée,
c’est-à-dire du point où serait exercée la
résultante des forces pour avoir le moment x
résultant. O
Un récipient contient de l’eau, de masse volumique 1,00 kg.L−1 sur une hauteur de 5,50 m et au-dessus,
de l’huile, de masse volumique 800 kg.m−3 sur une hauteur de 1,50 m.
Déterminer la pression exercée sur le fond du récipient par la colonne de liquide ainsi formée.
Données : accélération de la pesanteur g = 9,81 m.s−2 .
On adopte le modèle dit de l’atmosphère isotherme : la température y est uniforme, on la note T . L’air est
assimilé à un gaz parfait de masse molaire M = 29 [Link]−1 .
1. Établir la loi donnant la pression p(z) en fonction de l’altitude z. On notera p0 la pression au niveau
du sol et on exprimera P (z) en fonction de M , R, T0 ,p0 , g et z.
2. On note m(z) la masse de la colonne d’air, de section S, comprise entre les altitudes z = 0 et z.
(a) En utilisant la loi fondamentale de l’hydrostatique, établir la relation :
S[p0 − p(z)]
m(z) =
g
(b) En déduire l’expression de la masse totale mtot de l’atmosphère reposant sur la surface S.
(c) On définit la hauteur de l’atmosphère H comme étant l’altitude telle que m(H) = 0,99 mtot .
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
A p0
Soit un récipient ouvert sur l’air atmosphérique, de
z forme cylindrique, de masse m0 = 100 g, de hauteur
πA2
b = 15,0 cm et de section S = = 10,0 cm2.
4
b Ce récipient flotte à la surface de l’eau, de masse
volumique ρ. On repère le haut du récipient, par
rapport à la surface libre du liquide, par la côte
algébrique z.
1. Déterminer la valeur z0 que prend z à l’équilibre. Montrer que cette position d’équilibre est stable.
Donner la pulsation des petites oscillations du récipient autour de sa position d’équilibre. Réaliser les
applications numériques correspondantes.
2. On considère maintenant que le récipient est partiellement rempli d’eau. Quelle est la hauteur maxi-
male d’eau que peut contenir le récipient afin qu’il soit toujours en équilibre ?
Une bulle d’air de 5,00 mm de diamètre est émise au fond d’un étang. Arrivée à la surface, son diamètre
vaut 10,0 mm.
Déterminer la profondeur de l’étang en considérant que la température de l’eau est constante.
P2
2.1 Sous-systèmes
— Rappel —
Comme déjà vu page 55, on distingue trois échelles de longueur dans la description d’un fluide :
• l’échelle microscopique, définie au niveau moléculaire : c’est le libre parcours moyen ℓmoyen
caractérisant le mouvement des molécules du fluides, de l’ordre de 0,1 µm pour un gaz et de
0,1 nm pour un liquide.
• l’échelle macroscopique L, caractéristique de l’écoulement (diamètre d’une canalisation, lar-
geur d’une rivière, . . .).
• l’échelle mésoscopique a, intermédiaire entre les deux échelles précédentes : ℓmoyen <<
a << L.
Un volume de fluide de dimension mésoscopique a3 est suffisamment petit à l’échelle macro-
scopique pour être considéré comme ponctuel, et suffisamment grand à l’échelle microsco-
pique pour contenir un grand nombre de molécules. Typiquement, ce volume est de l’ordre
de la dizaine (ou de quelques dizaines) de µm.
L’échelle mésoscopique permet de définir des grandeurs locales, définies statistiquement sur
le grand nombre de molécules du volume a3 situé au point considéré.
Une particule de fluide constitue un système fermé, de volume mésoscopique dV(M) = a3 . Elle
permet de définir statistiquement des grandeurs locales intensives (température, pression, masse
volumique, vitesse).
L’approche mésoscopique définit une échelle intermédiaire entre le macroscopique et le micro-
scopique. À cette échelle, un petit volume est assez grand pour contenir un très grand nombre de
particules et on peut considérer le milieu comme continu mais suffisamment petit pour considérer
la densité de molécules comme uniforme.
Le premier principe est un principe de conservation. Il existe une grandeur, appelée énergie totale, extensive
et conservative, que l’on peut définir dans tout système fermé. On appelle énergie totale d’un système, toute
l’énergie présente, sous toutes ses formes.
Etotale = Ep + Ec + U + Enucl + Eautre
Etotale = Ec + U
Les autres énergies seront considérées comme constantes et leurs variations seront donc nulles.
Les unité s :
L’énergie est une fonction d’état. Elle ne dépend pas du chemin suivi, mais uniquement de l’état initial et
de l’état final.
La variation de cette énergie est donnée par :
∆Etotale = W + Q
On obtient alors :
∆U + ∆Ec = W + Q
dU + dEc = δW + δQ
Les unité s :
Le première loi de Joule postule que l’énergie interne U d’un gaz parfait est une fonction qui
dépend uniquement de la température : U = f1 (T ).
dU = CV dT = nCV,m dT = mcV dT
Avec Cv la capacité thermique à volume constant (CV,m et cV sont les capacités thermiques respectivement
molaires et massiques à volume constant).
On utilise donc U dans un système isochore.
De même, la deuxième loi de Joule postule que l’enthalpie interne H est également une fonction
qui dépend uniquement de la température : H = f2 (T ).
dH = Cp dT = nCp,m dT = mcp dT
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Avec Cp la capacité thermique à pression constante. (Cp,m et cp sont les capacités thermiques molaires et
massiques à pression constante).
On utilise donc H dans une transformation monobare.
Cp,m Cp cp
On définit la grandeur γ = = = >1
CV,m CV cV
— Relation de Mayer —
Cp − CV = nR
donc
Cp,m − CV,m = R
ou encore
n R
cp − cV = R= =r
m M
Pour les phases condensées (solides ou liquides), on peut faire l’hypothèse que U (T ) et H(T ) sont
équivalentes.
On obtient alors :
Cp = CV = C
Et :
2.5.1 Entropie
— Entropie —
Pour tout système fermé, on peut définir une fonction d’état, notée S, appelée entropie, qui est
une grandeur extensive, mais non conservative. Cette grandeur peut être créée ou non, mais jamais
détruite.
Ainsi, en pratique, toute transformation d’un système thermodynamique s’effectue avec augmen-
tation de l’entropie totale incluant l’entropie du système ainsi que celle du milieu extérieur. On dit
alors qu’il y a création d’entropie notée Scréation ou Scréée .
— Bilan entropique —
— Remarque —
Dans le cas d’une transformation réversible, l’entropie créée Scréée = Scréation est nulle.
dU = T dS − pdV
— Remarque —
dU = δW + δQ
δW = −pdV etδQ = T dS
La définition de l’enthalpie :
H = U + pV
permet d’écrire :
dH = dU + pdV + V dp
Soit :
dH = T dS + V dp
— Tempé rature —
À l’aide de l’expression de dS dans le cas d’un système caractérisé par T , V , en identifiant les
expressions suivantes :
— Pression —
dU = T dS − pdV
Å ã Å ã
∂U ∂U
=− dS + dV
∂S V ∂V S
Scréée = ∆S − Séchangée
dU = nCV,m dT
et la variation d’enthalpie :
dH = nCp,m dT
Par conséquent, pour un gaz parfait G.P. la variation d’entropie élémentaire dS vaut :
dT dV
dS = nCV,m + nR
T V
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
ou encore :
dT dp
dS = nCp,m − nR
T p
En intégrant entre les valeurs initiales (pi , Ti , Vi ) et les valeurs finales (pf , Tf , Vf ), on retrouve les formules
vues en première année :
Z f
δQ Tf Vf
∆S = = nCV,m ln + nR ln
i T Ti Vi
ou encore :
Z f
δQ Tf pf
∆S = = nCp,m ln − nR ln
i T Ti pi
En utilisant les lois de Joule ainsi que les identités thermodynamiques, on démontre la loi de Laplace.
— Loi de Lapla e —
pV γ = cte
T γ p1−γ = cte′
ou bien :
T V γ−1 = cte′′
⇒ Activité 2. 1
dT dT dT
dS = C = nCm = mc
T T T
En intégrant entre les températures initiale Ti et finale Tf , on retrouve les formules vues en première année :
Z Tf
dT Tf Tf Tf
∆S = C = C ln = nCm ln = mc ln
Ti T Ti Ti Ti
Soit Qf le transfert thermique de la source froide vers le système, en contact avec un thermostat de tempé-
rature Tf et Qc le transfert thermique de la source chaude vers le système en contact avec un thermostat de
température Tc .
◦ Qf < 0
◦ Qc > 0
◦ W <0
◦ Qf > 0
◦ Qc < 0
◦ W >0
Attention aux conventions : Si par exemple, la source chaude est constituée d’eau de capacité
calorifique massique ce dont la température est égale à Te et variant de Te,I à Te,F , on a alors la
relation suivante :
Z Te,F
Qc = − ce dTe
Te,I
car Qc est la chaleur reçue par le système de la part de la source chaude. C’est donc la chaleur cédée
par l’eau chaude, soit l’opposé de la chaleur reçue par celle-ci !
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Dans cette étude, sachant que le système est en contact avec deux thermostats, on a les formules suivantes :
Premier principe :
∆Ucycle = W + Qf + Qc = 0
Deuxième principe :
Qc Qf
∆Scycle = + + Scréation = 0
Tc Tf
Comme l’entropie créée Scréation est positive ou nulle, on en déduit l’inégalité de Clausius :
Qc Qf
+ 60
Tc Tf
— Remarques —
Comme précédemment, l’entropie créée Scréation est positive ou nulle et on en déduit l’inégalité de
Clausius :
Z Z
δQc δQf
+ 60
Tc Tf
— Effi a ité —
Énergie utile
η=
Énergie dépensée
Les unité s :
Dans une machine frigorifique, un fluide caloporteur (CF C, N H3 , ...) circule. Sous l’action du compres-
seur, le fluide caloporteur passe de sa phase gazeuse à sa phase liquide. Ce faisant, il restitue de l’énergie
à la source chaude. Sous l’action du détendeur, le fluide caloporteur passe de sa phase liquide à sa phase
gazeuse. Ce faisant, il absorbe de l’énergie au contact de la source froide.
La machine frigorifique fonctionne donc à l’aide des transitions de phases, qui réclament des transferts
thermiques.
Le fonctionnement en machine frigorifique est caractérisé par :
• W > 0,
• Qc < 0,
• Qf > 0,
Qc
• η=− .
W
La pente de la tangente à la courbe représentative d’une adiabatique réversible (donc isentropique) est plus
grande que celle représentant une isotherme.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Pour un gaz parfait, une transformation isotherme est décrite par la loi du gaz parfait :
pV = nRT = cste
dp dV
+ =0
p V
On obtient la pente :
dp p
=−
dV V
Pour une adiabatique réversible, donc isentropique, la première identité thermodynamique donne :
dU = T dS − pdV = −pdV
nR
dU = CV dT = dT
γ−1
On en déduit :
pdV + V dp
dU = = −pdV
γ−1
Soit :
pdV (1 + γ − 1) = −V dp
Et :
dp p
= −γ
dV V
Å ã Å ã
∂p ∂p
=γ
∂V S ∂V T
Les isothermes du gaz parfait sont des hyperboles, l’isotherme Tc étant située au-dessus de l’isotherme Tf .
Les deux isothermes sont reliées par deux adiabatiques, de pentes plus importantes, d’où l’allure du cycle :
p(P a)
2
P2
1′ Tc
P1 1
Tf
V (m3 )
Le rendement de ce moteur correspond au rapport du travail fourni sur un cycle par le moteur à la chaleur
qu’il absorbe depuis la source chaude pour pouvoir fonctionner :
W Wfourni
η=− =
Qc Qc
∆Ucycle = 0 = Qc + Qf + W
Qf Tf
η =1+ =1−
Qc Tc
Ce rendement est indépendant de la nature du fluide, puisque seules interviennent les températures des
sources. Pour un moteur fonctionnant entre les deux mêmes sources, le second principe s’écrit, dans le cas
général, avec Scréée :
Qc Qf
∆Scycle = 0 = + + Scréée
Tc Tf
Soit :
Qc Qf
+ = −Scréée 6 0
Tc Tf
On en déduit :
Qf Tf
6−
Qc Tc
et donc :
W Qf Tf
η=− =1− 61− = ηmax
Qc Qc Tc
Le rendement du cycle de Carnot est maximal pour des sources dont la température est fixée.
— Diagramme entropique —
δQ
∆S = + Scréée
T
Si la transformation est adiabatique , Q = 0 et si elle est réversible, Scréée = 0.
On en déduit ∆S = O : la transformation est alors isentropique.
Dans le diagramme entropique, les isothermes 1 → 1′ et 2 → 2′ sont horizontales, tandis que les transfor-
mations adiabatiques réversibles, donc isentropiques 1′ → 2 et 2′ → 1 sont des verticales :
T (K)
Tc
2 2′
Tf
1′ 1
s( J.K −1 .kg −1 )
smin smax
2.9.1 Énoncés
• une compression pendant laquelle un gaz parfait diatomique passe de l’état (p1 = 1,0 bar, T1 =
293 K) à l’état (p2 = 9,0 bar, T2 = 609 K),
• une détente au cours de laquelle ce même gaz passe de l’état (p3 = 8,5 bar, T3 = 1300 K) à l’état
(p4 = 1,0 bar, T4 = 793 K).
2. Que vaudrait k si les transformations étaient isentropiques (adiabatiques réversibles) ? Quelles consé-
quences peut-on en tirer ?
3. Rappeler la première loi de joule pour un gaz parfait. À l’aide de la première identité thermodyna-
mique, établir l’expression de la différentielle de l’entropie massique ds.
dT
4. En fait, on peut poser ds = c . Établir une relation liant c, cp , γ et k.
T
5. Que vaut k pour une transformation :
(a) isentropique ?
(b) isobare ?
6. Pour chacune des transformations, exprimer l’entropie créée et effectuer les applications numériques.
Un briquet à air permet de réduire un volume V0 d’air, initialement à T0 = 12,0 oC et p0 = 1,00 bar, au
centième de sa valeur.
1. En prenant un modèle adiabatique réversible pour cette transformation, déterminer la pression et la
température finales si γ = 1,4.
2. Calculer ensuite la variation d’énergie interne et d’enthalpie de cette masse d’air si V0 = 2,00 cm3 .
On étudie le chauffage d’une école pendant une journée d’hiver. On appelle Text la température de l’air
extérieur. On suppose qu’à chaque instant, toute l’école est à la même température T . Quand l’école reçoit
une quantité de chaleur élémentaire δQ, sa température T varie de dT . On note C la capacité thermique de
toute l’école. On suppose que la chaleur perdue par l’école (à cause des déperditions thermiques à travers
les murs, le toit, . . .) pendant la durée dt est égale à δQperdue = aC (T − Text ) dt, a étant une constante.
On donne Text = 263 K, C = 7,60.107 J.K −1 , a = 7,90.10−5 s−1 .
Dans le cycle de Diesel, le premier temps (représenté par AB) est une compression adiabatique de l’air seul,
avec un rapport volumétrique assez élevé. Le carburant n’est injecté dans le cylindre qu’à partir de B ; la
température est alors suffisante pour que le mélange s’enflamme spontanément (sans l’aide de bougies). Le
taux d’injection est réglé de telle façon que la pression reste constante pendant la détente BC ; on néglige
alors la quantité de carburant injectée et on raisonne comme avec un gaz unique recevant de l’extérieur une
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
B C
En considérant les gaz comme parfaits, exprimer le
Cp
rendement du cycle en fonction de γ = (supposé
Cv
indépendant de la composition du mélange),
Å du ã
rap-
VA
port volumétrique de la compression a = et D
Å VBã
VA
du rapport volumétrique de la détente b = . A′ A
VC
V
Écrire la relation fondamentale de la dynamique pour la bille et déterminer la pulsation ω0 ainsi que la
période T0 des oscillations en fonction de pe , V0 , s, m et γ.
P3
On assimile un système à un corps pur quand on peut considérer que ce système est constitué
d’entités chimiques de même nature. Le corps pur est un modèle, dans lequel tout corps étranger est
supposé infiniment dilué.
On appelle corps pur simple un corps composé d’une unique entité comme H2 ou O2
— Système diphasé —
On dit qu’un système est diphasé s’il contient un corps pur sous deux phases (Eau liquide et gaz
par exemple).
[Link] Transition
• solide,
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
• liquide,
• gaz.
En voici un résumé :
Gaz
Co blim
risa tion
n
Su
nd atio
tio
c
Va uéfa
en s
atio n
Liq
po
Liquide Solide
Solidification
Fusion
L/G
S/G
T
G
p(bar)
LIQUIDE
N M
p2
p1
SOLIDE
T
VAPEUR
T (K)
0 T1 T2
F IGURE 3.3 – Diagramme p = f (T )
p(bar)
LIQUIDE
N M
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
p2
p1
SOLIDE
T
VAPEUR
T (K)
0 T1 T2
F IGURE 3.4 – Diagramme p = f (T )
Pour l’eau, la pente de la courbe de fusion est négative, contrairement aux autres corps purs.
La pente de la courbe est toujours croissante, c’est à dire qu’une augmentation de T impose une augmenta-
tion de pression dans un corps pur.
Il existe une exception. Pour l’eau, la pente de la frontière Solide → Liquide est négative.
— Point triple —
T , appelé point triple, correspond au domaine de coexistence des trois phases du corps pur consi-
déré. Ce point est une propriété intrinsèque du corps. C’est le point triple de l’eau pur qui définit
l’échelle de température. On le fixe à 273,15 K.
— Point ritique —
Le point critique C est la limite de coexistence des phases liquide et gazeuse d’un corps pur. Au
dessus de ce point, on passe de manière continue de l’état gazeux à l’état liquide sans observer
d’état diphasé. On observe ce qu’on appelle un état fluide.
N M VAPEUR
ps (T2 )
LIQUIDE L LIQ ⇄ VAP V
ps (T1 ) isotherme T2
isotherme T1
[Link] Commentaires
Parcourons une isotherme (appelée isotherme d’Andrews) placée sous le point critique en augmentant le
volume :
La transformation qui décrit la courbe LV s’effectue à pression constante. Cette pression, qui est notée
ps (T ), est la pression de vapeur saturante. Elle augmente avec la température.
Fluide supercritique Lorsque T > TC et que p > pC , on dit alors que le système est à l’état de fluide
supercritique : on passe alors de l’état liquide à l’état gazeux de manière continue.
Opalescence Optique Au point critique, pour T = Tc , χt tend vers ∞. Ceci implique donc que la densité
du milieu connaît de très fortes fluctuations. On montre que l’indice d’un milieu est lié à la densité de celui-
ci. Dans ces conditions, le corps pur diffuse de manière particulière la lumière, on appelle ce phénomène
opalescence optique.
Stockage des fluides D’après l’équation d’état de l’eau liquide, on à montré qu’une faible augmentation
de température, lors d’une évolution isochore, entrainait une très forte augmentation de pression. Pour cette
raison, on stocke les fluides dans un état diphasé. On les stocke dans un volume V supérieur au volume Vc
du point critique, car une augmentation de température entraîne un augmentation de pression relativement
faible.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
[Link] Enthalpie
— Proprié té —
∆H = mℓ = L = Qp
— haleur latente —
On appelle chaleur latente (ou enthalpie) de changement d’état d’une masse ou quantité de matière
la quantité d’énergie qu’il faut lui fournir à pression constante pour qu’elle passe de l’état initial
(solide, liquide ou gazeux) à un autre état.
La variation d’enthalpie d’un corps pur au court d’une transition de phase est appelée chaleur latente de
transition de phase, notée L.
De plus, nous savons que Q est positif pour une transition de phase d’un état ordonné vers un état moins
ordonné. Ce qui implique que dans ce cas, L est positif. Dans la situation inverse, L est négatif, tout comme
Q.
— Exemples —
Chaleurs latentes :
• de vaporisation du fréon à 263 K : 159 [Link] −1
• de vaporisation de l’eau à 373 K : 2258 [Link] −1
• de vaporisation de l’eau à 485 K : 1892 [Link] −1
• de fusion de la glace à 273 K : 335 [Link] −1
Dans le cadre de cette évolution isotherme et isobare, la variation d’entropie d’une masse m de corps pur
est donnée par :
mℓ
∆S =
T
Si on exprime
LM
xv =
LV
MV
xℓ = 1 − xv =
LV
— Remarque —
Cette règle des moments peut également s’exprimer en fonction d’autres grandeurs, par exemple :
H = mℓ hℓ + mv hv = mh, soit h = xv hv + (1 − xv )hℓ .
On en déduit :
h − hℓ
x=
hv − hℓ
• Lorsqu’un système est diphasé (liquide-vapeur), le volume du liquide est en principe négligeable
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Les systèmes industriels utilisent des fluides réels pour lesquels il est impossible de trouver des équations
d’état. Par conséquent, on utilise :
⇒ Activité 3. 2
Considérons la phase motrice d’un cycle de machine à vapeur pendant laquelle de la vapeur d’eau se détend
dans un cylindre fermé par un piston mobile.
On pourra considérer cette détente suffisamment rapide pour que les transferts thermiques n’aient pas le
temps d’être quantitatifs. Ainsi, la détente pourra être considérée comme adiabatique. Pour simplifier, sup-
posons qu’elle est également réversible, donc isentropique. L’état initial est constitué de vapeur juste satu-
rante, à à la pression p1 = 20 bar et à la température T1 = 485 K. Dans l’état final, l’eau est sous forme de
mélange liquide-vapeur et se trouve à la pression p2 = 1,0 bar et à la température T2 = 373 K.
1. Calculons le titre massique en vapeur xv,2 dans l’état final à partir des données suivantes :
−3 −2 3
485 20 1,18.10 909 2,45 9,98.10 2,801.10 6,35
−3 3
373 1,0 1,04.10 418 1,30 1,70 2,676.10 7,36
ln p
C
b
ln p2
liquide vapeur
ln p1
liquide + vapeur
h
Saint Joseph - LaSalle CPGE TSI
92 D IAGRAMMES D ’ ÉTAT DES FLUIDES RÉELS PURS
R
Dans le domaine vapeur, on considère que celle-ci est un gaz parfait et dans ce cas, en notant r = , on
M
peut écrire :
dh = cp dT
dT dp
ds = cp −r
T p
dT dV
ds = cV +r
T V
dp p nR
=− =−
dV V T
Ce sont des courbes à allure exponentielle décroissante,
• les isenthalpes sont isothermes,
• les isentropiques, comme démontré dans le chapitre précédent, sont définies par :
dp p nR
= −γ = −γ
dV V T
Ce sont des courbes à allure exponentielle décroissante, de pente plus grande que les isothermes.
Dans le domaine liquide, on considère que celui-ci est une phase condensée et on peut écrire :
dh = cdT
dT
ds = c
T
dh = ℓv dx
ℓv dx
ds =
Tv
R
Dans le domaine vapeur, on considère que celle-ci est un gaz parfait et dans ce cas, en notant r = , on
M
peut écrire :
dh = cp dT
dT dp
ds = cp −r
T p
dT dV
ds = cV +r
T V
• les isothermes sont horizontales,
• les isenthalpes sont également horizontales,
• les isentropiques sont verticales,
dT dT T
• les isobares vérifient ds = cp et ont donc une pente = : ce sont donc des branches
T ds cp
d’exponentielles croissantes,
dT dT T
• les isochores vérifient ds = cV et ont donc une pente = : ce sont donc des branches
T ds cV
d’exponentielles croissantes.
Dans le domaine liquide, on considère que celui-ci est une phase condensée et on peut écrire :
dh = cdT
dT
ds = c
T
• les isothermes sont horizontales,
• les isenthalpes sont également horizontales,
• les isentropiques sont verticales,
dT dT T
• les isobares vérifient ds = c et ont donc une pente = : ce sont donc des branches d’expo-
T ds c
nentielles croissantes,
dT dT T
• les isochores vérifient ds = c et ont donc une pente = : ce sont donc des branches
T ds c
d’exponentielles,
dh = ℓv dx
ℓv dx
ds =
Tv
• les isothermes sont horizontales,
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
3.4.1 Énoncés
Données à po = 1,00.105 P a :
• ℓvap = 2,660.103 [Link] −1 à Tvap = 373 K,
• ℓfus = 334 [Link] −1 à Tfus = 273 K,
• c = 4,20 [Link] −1.K −1 .
P4
Fluides en écoulement stationnaire
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
— Rappel —
En mécanique des fluides, on sera amené à étudier des systèmes macroscopiques, c’est-à-dire
constitué d’un grand nombre de molécules ou d’atomes.
Comme déjà vu page 67, trois échelles de description sont possibles :
En effet, parler de molécules et d’atomes consiste à décrire la matière à l’échelle microscopique.
C’est une modélisation discrète de la matière (par opposition à modélisation continue), à l’échelle
du nm. Or, l’étude d’un fluide ne peut pas consister à déterminer le mouvement de chacune de
ses molécules. Non seulement ce serait impossible techniquement, mais on ne saurait que faire de
toute cette information. C’est pourquoi en mécanique des fluides, on ne va s’intéresser qu’aux pro-
priétés macroscopiques des systèmes étudiés. L’échelle de description que l’on va adopter est donc
l’échelle macroscopique. À cette échelle, on va adopter une modélisation continue de la répartition
de matière, comme si la matière n’était pas constituée d’atomes, mais emplissait tout l’espace.
Souvent, les grandeurs définies localement (vitesse du fluide, masse volumique) n’ont pas néces-
sairement la même valeur partout dans le fluide.
Ces grandeurs sont définies à l’échelle mésoscopique.
L’échelle mésoscopique, de l’ordre de la dizaine de µm, est assez grande devant l’échelle mi-
croscopique pour adopter une modélisation continue de la matière et assez petite devant l’échelle
macroscopique pour considérer ces grandeurs comme localement uniformes.
On considère, autour d’un point M , un petit volume dV(M) qui contient donc un grand nombre de
molécules.
Ne pas confondre molécule et particule de fluide : une particule de fluide contient un grand nombre
de molécules.
— Exemples —
• On photographie en pause longue un objet en mouvement, par exemple les phares d’une
voiture la nuit.
• Un observateur suit des yeux une feuille à la surface de l’eau s’écoulant dans le lit d’une
rivière.
— Exemples —
• Un observateur regarde fixement une zone de la rivière et voit passer des objets à la surface
de l’eau quand ils traversent son champ de vision.
• On photographie à un instant t un objet en mouvement. Typiquement, une carte des vents est
une description eulérienne :
— Tube de ourant —
En description lagrangienne, le vecteur vitesse d’un point M du fluide représente le vecteur vitesse de la
particule de fluide qui l’entoure alors qu’en description eulérienne, c’est le vecteur vitesse de la particule
située en M à un instant t.
La trajectoire retrace l’histoire d’une particule de fluide alors que la ligne de courant représente un instan-
tané du champ des vitesses.
Comme ce sera vu plus tard, un écoulement stationnaire est tel que les champs définis dans le fluide
sont indépendants du temps, en particulier celui des vitesses.
Dans un tel écoulement, le champ des vitesses eulérien ne dépend pas explicitement du temps et
dans ce cas, lignes de courant et trajectoires coïncident. En effet, dans le cas d’un écoulement
~v (M
,t)
′
,t)
~v(M
ant
ur
~v (r~0 ,t)
co
de
e
gn
,t 0) Li
~v (r~0
Variables r~0 (position initiale) et t M et t
L’observateur suit la particule fixe le point M
Idéal pour la mécanique du point l’étude des écoulements.
Les unité s :
— Fluide in ompressible —
• les liquides,
• les gaz, lorsque leur vitesse d’écoulement est suffisamment faible, c’est-à-dire lorsque leur vitesse v
est très inférieure à la vitesse du son : v ≪ vson .
— Dé bit volumique —
On appelle débit volumique DV à travers une surface S ~ orientée, le volume δV de fluide qui
traverse cette surface (S) par unité de temps. Le débit volumique est donc une grandeur algébrique,
~ et négative dans le cas contraire.
positive si la vitesse ~v du fluide est de même sens que S,
ZZ
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
~ = δV
DV = ~v .dS
(S) dt
Les unité s :
— Dé bit massique —
On appelle débit massique Dm à travers une surface S ~ orientée, la masse δm de fluide qui traverse
cette surface (S) par unité de temps. Le débit massique est donc une grandeur algébrique, positive
~ et négative dans le cas contraire.
si la vitesse ~v du fluide est de même sens que S
ZZ
Dm = ~ = δm
ρ~v · dS
(S) dt
Les unité s :
Les unité s :
Le flux du vecteur densité de courant massique à travers une surface est égal au débit massique :
Σ
−
→ −
→
jm (x) V0 jm (x + dx)
Et
ZZ
−−−−−−−→ −−→
Dm,s = − jm (x + dx) · dSs
Soit :
ZZZ Å ã
−
→ ∂ρ
div jm + dV = 0
(V0 ) ∂t
Ceci étant vrai quelque soit le volume V0 , cela conduit à l’équation de conservation de la masse :
—É oulement stationnaire —
Un écoulement est stationnaire (ou permanent) si en un point donné, son champ des vitesses ainsi
que sa masse volumique ne varient pas :
∂−
→v ∂ρ
= 0 et =0
∂t ∂t
En régime stationnaire, la trajectoire des particules et les lignes de courant sont confondues (mais ce n’est
pas le cas en général).
−
→
De l’équation de conservation de la masse, on déduit que div jm = 0 . Avec :
ZZZ {−
−
→ → − →
div jm dV = jm · dS
(V0 ) Σ
{−
→ − → −
→
On obtient jm · dS = 0 . Cela signifie que le vecteur jm est à flux conservatif : le débit massique est
(Σ)
identique à travers tout tube de courant.
Le débit massique étant constant, on peut écrire, pour deux sections S1 et S2 d’une conduite :
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
ρ1 S1 v1 = ρ2 S2 v2
ou encore :
ρSv = cte
En différenciant, on obtient :
dρ S v + ρ dS v + ρ S dv
et en divisant par ρSv :
dρ dS dv
+ + =0
ρ S v
⇒ Activité 4. 1
Retrouver l’expression ci-dessus en utilisant les dérivées logarithmiques :
Soit :
DV1 = DV2
DV
On a alors, en définissant la vitesse moyenne par v = :
S
v1 S1 = v2 S2
Le resserrement des lignes de courant provoque une augmentation de la vitesse moyenne.
— Rappel —
Comme déjà signalé page 99, la frontière "compressible - incompressible" est liée à la vitesse de
l’écoulement :
• si le gaz possède une vitesse inférieure à 0,3 vson , le gaz peut être considéré comme incom-
pressible,
• sinon, il est considéré comme compressible.
—É oulement irrotationnel —
Un écoulement irrotationnel est un écoulement dans lequel il n’y a pas d’effets tourbillonnaires :
∆φ = 0
—É oulement tourbillonnaire —
~ tel que :
Un écoulement tourbillonnaire est caractérisé par le vecteur tourbillon non nul Ω
Imaginez que vous êtes devant une boulangerie en train de faire la queue, vous voyez des gens entrer les un
après les autres et d’autre sortir avec leur baguette . . .
Vous vous demandez par exemple combien de personnes sont à l’intérieur de la boulangerie et si le nombre
de ces personnes varie au cour du temps.
La réponse à cette question est simple : il suffit de compter le nombre de personnes qui entrent pendant une
durée ∆t (par exemple 5 min) et le nombre de personnes qui sortent au cours de la même durée.
Si vous comptez 5 personnes qui entre mais seulement 4 personnes qui sortent pendant la même durée, vous
en déduirez que le nombre de personnes à l’intérieur de la boulangerie à augmenté de 1 personne pendant
cette durée.
Imaginons maintenant qu’il y ait 3 portes d’entrée dans la boulangerie ainsi que 3 portes de sortie.
Les gens peuvent entrer par n’importe laquelle de ces 3 portes et sortir par n’importe laquelle des 3 portes
de sortie.
Appelons x les portes d’entrée/sortie situées sur la gauche de la boulangerie, y les portes d’entrée/sortie
situées au milieu et z, celles de droite.
Pour faire le calcul précédent vous devrez faire votre raisonnement sur chaque couple de portes : combien
de gens entrent par la porte d’entrée x et sortent par la porte de sortie x ? De même pour y et z, toujours
pendant la durée ∆t.
Admettons que 3 personnes entrent en x et 1 personne sorte en x. Le bilan vaut ∆px = +2 pour la
boulangerie : c’est la variation pour la porte x.
Vous constatez que la porte y laisse entrer 4 personnes mais en fait sortir 7 : la boulangerie a perdu 3
personnes par les portes y. ∆py = −3 .
Les portes z permettent à 5 personnes d’entrer et à 6 personnes de sortir. La variation vaut ∆pz = +1 .
Dans un écoulement divergent, le débit volumique ne se conserve pas le long de tube de courant.
Ainsi, le volume d’une particule varie, ce qui signifie que la masse volumique change également.
Mathématiquement, si div − →v > 0(< 0), les lignes de champ convergent (divergent).
La divergence caractérise comment un champ évolue dans sa propre direction car cet opérateur fait
intervenir des dérivées partielles non-croisées.
Si div −
→v 6= 0, on dit que le champ possède une source ou un puis de champ : il est à flux non
conservatif.
Si div −
→v = 0, le flux est dit à champ conservatif.
— Exemple —
Si en un point donne, le rotationnel du champ de vitesse est non nul, c’est que localement, les particules de
fluides subissent une rotation sur elles-mêmes. Attention à ne pas tout confondre avec le tourbillon.
Tout champ de vitesse tourbillonnaire n’est pas nécessairement rotationnel, qui peut être nul, alors qu’un
écoulement plan et rectiligne peut, lui, être à rotationnel non nul.
En effet, faites un tourbillon dans un seau d’eau : la théorie prévoit que loin du centre, l’écoulement est
irrotationnel. En fait, les particules d’eau subissent juste une translation circulaire. Une allumette ne tourne
pas sur elle-même : elle suit le mouvement de translation circulaire.
En revanche, près du centre, l’écoulement devient rotationnel, les particules tournent sur elles-mêmes et
l’allumette précédente tourne sur elle-même.
Prenons un fluide visqueux (de l’huile par exemple) dans une bouteille et inclinons la bouteille.
L’écoulement est plan et rectiligne. Or le rotationnel du champ de vitesse est non nul. En effet, Sur la paroi
de la bouteille, il y a continuité du champ de vitesse, donc celui-ci est nul. Or en surface il est non nul ;
il y a donc un gradient de vitesse selon la profondeur. Si x est la composante parallèle a la paroi et z la
profondeur, la dérivée croisée est non nulle. Le rotationnel est donc non nul.
Physiquement, considérons une particule de fluide. Sur la face de hauteur z + dz, elle possède une vitesse
différente de celle sur la face d’altitude z : elle va se mettre à tourner sur elle-même.
Dans un écoulement rotationnel, les lignes de courant se referment sur elles-mêmes (mais ce n’est
pas une condition suffisante).
−→ → − → −→ → −→ − →
−
Mathématiquement, si rot − v 6= 0 , les lignes de champ s’enroulent autour de rot −
v . Si rot →
v = 0,
les lignes de champ n’ont pas de courbure.
— Exemple —
4.8.3 Applications
O (∆)
×
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Cte Cte
Cas a : v = Cas b : v =
r2 r
(∆)
×
Cas c Cas d : v = rω
Cas e Cas f
⇒ Activité 4. 2
Pour les exemples précédents, indiquer dans le tableau ci-dessous si les écoulements sont divergents ou non,
rotationnels ou non.
−→ −
Cas div −
→
v rot →
v
a
b
c
d
e
f
4.9.1 Énoncés
+M R
r
−−→
v(M)
Données :
En coordonnées cylindriques :
−
→ 1 ∂(rDr ) 1 ∂Dθ ∂Dz
div D = + +
r ∂r r ∂θ ∂z
Å ã Å ã Å ã
−→ −→ 1 ∂Dz ∂Dθ − → ∂Dr ∂Dz → − 1 ∂(rDθ ) ∂Dr − →
rot D = − er + − eθ + − ez
r ∂θ ∂z ∂z ∂r r ∂r ∂θ
Un verre contenant de l’eau (fluide newtonien) est posé sur un plateau tournant.
P5
Fluides en écoulement laminaire
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
stationnaire
—É oulement laminaire —
Dans le chapitre 1, nous avons vu que dans un fluide au repos, les forces de pression sont normales à
l’élément de surface sur lequel elles s’exercent.
Pour un fluide en mouvement, il faut envisager d’autres forces : celles relatives à la viscosité. Ainsi, dans
les fluides "pâteux" (miel, verre en fusion, . . .), ces phénomènes de viscosité deviennent prépondérants.
Le cas limite idéal dans lequel les phénomènes de viscosité sont négligeables correspondent au fluide par-
fait :
— Fluide parfait —
ou encore :
v(M),⊥ = v(Mparoi ),⊥
5.1 Viscosité
Dans un fluide réel en écoulement, les forces de contact ne sont pas perpendiculaires aux éléments de
surface sur lesquels elles s’exercent.
Considérons un fluide réel, en écoulement unidimensionnel laminaire dont le champ des vitesses est de la
forme : −
→v = v(z) −→
ex .
v~0
Les actions de contact exercées par le fluide de la couche (2) sur le fluide de la couche (1) à travers la surface
dS se décomposent en une composante normale et une composante tangentielle :
−−−−→ −−−−−→ −−−−−→
dF2→1 = dFn,2→1 + dFt,2→1
La composante normale définit la pression p(M, t) :
−−−−−→ −
→
dFn,2→1 = p(M, t)dS
La composante tangentielle est appelée force de viscosité ou force de cisaillement. Dans le cas d’un fluide
dit newtonien, elle s’écrit :
Å ã
−−−−−→ ∂v
dFt,2→1 = −η dS −
→
ex
∂z
dz
Ä ä
→
−
v z+
dS
Ä dz
ä 2
→
−
v z−
~
dF 2
x
F IGURE 5.3 – Viscosité : champ des vitesses
— Fluide newtonien —
Un fluide newtonien est un fluide dans lequel les forces visqueuses sont proportionnelles à la diffé-
rence de vitesse, c’est-à-dire au gradient des vitesses.
La contrainte (force surfacique) σ exercée par une couche supérieure sur un élément de surface inférieure,
s’écrit :
−→
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
−
σ−−→ = dF = σ ~n + σ ~t
1→2 n t
dS
— Contrainte de isaillement —
∂v ∂v
où désigne le gradient des vitesses dans la direction normale à la surface, par exemple .
∂n ∂y
Pour un fluide newtonien, les conditions aux limites sont nettement différentes de celles d’un fluide parfait.
En effet, pour un point voisin de la surface d’une conduite :
−−→ = −
v(M)
→
0
−−→ − −
v(M) v−−−→ − →
(Mparoi ) = 0
ou encore :
v 6= 0 v=0
u
u
v maximale
u
— Remarque —
Évidemment, le fluide de la couche (1) exerce sur le fluide de la couche (2) une force opposée et
égale en module :
−−−−−→ −−−−−→
dFt,1→2 = −dFt,2→1
— Vis —
La viscosité d’un fluide peut se définir comme la propriété de résistance d’un fluide à l’écoulement.
Bien évidemment, l’écoulement d’un fluide visqueux ne peut être qu’irréversible, en raison des frottements
qui en découlent.
Les unité s :
La viscosité cinématique d’un fluide est le rapport entre son cœfficient de viscosité et sa masse
volumique :
η
ν=
ρ
La viscosité cinématique tient compte non seulement de la viscosité, mais également de l’inertie du fluide.
Un ordre de comparaison : νair ≃ 15νeau .
— Exemple —
Un cas très utilisé est celui d’une sphère de rayon r animé d’une vitesse v faible dans un fluide
(vitesse par rapport à celui-ci). Celle-ci subit de la part de ce fluide une force de frottement :
−
→
Ff = −6 π η r −
→
v
L’écoulement de l’eau à la sortie d’un robinet permet de mettre en évidence des phénomènes différents :
• à des débits faibles, mais continus, l’écoulement est stable : les lignes de courant sont disposées de
façon régulière les unes par rapport aux autres. On est alors en présence d’écoulement laminaire.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
• si on augmente le débit, des gouttes sont émises en périphérie de l’écoulement principal, de façon
périodique : l’écoulement est alors plus instable.
• si on augmente encore le débit, il est possible d’observer une alternance régulière de grosses et de
petites gouttes. Il peut également y avoir des zones dans lesquelles l’écoulement est irrégulier. L’écou-
lement est alors turbulent. Ce phénomène de turbulence peut augmenter : l’écoulement peut en effet
devenir très chaotique.
Considérons une cuve assez haute au fond de laquelle est percée un orifice de petite taille.
Suivant la nature du liquide, par exemple de l’huile ou de l’eau, la forme du jet est différente. Ainsi, l’écou-
lement de l’huile, lent, donne un jet vertical, stable et cylindrique. À l’inverse, de l’eau donne un jet plus
rapide, turbulent, dans lequel des gouttelettes se détachent de l’écoulement principal.
On voit ici que la viscosité influe sur la nature de l’écoulement.
— Nombre de Reynolds —
C’est un nombre sans dimension, caractérisant l’écoulement d’un fluide dans une conduite ou au-
tour d’un obstacle.
Son expression est :
avec :
• ρ la masse volumique du fluide,
• v sa vitesse,
• η sa viscosité,
• d une longueur caractéristique du système : le diamètre de la canalisation ou de l’obstacle.
Lorsque le nombre de Reynolds est faible, l’écoulement est laminaire. Lorsqu’il est élevé, il est turbulent.
• Aux valeurs intermédiaires du Reynolds (entre 1 et 2000 environ), les forces d’inertie sont prépon-
dérantes, mais l’écoulement reste laminaire. Cependant, il n’est plus réversible : si l’on stoppe les
Il faut cependant comprendre que ces valeurs sont indicatives car elles dépendent de la géométrie de l’écou-
lement.
— Exemples —
• Dans une conduite, l’écoulement est laminaire lorsque le nombre de Reynolds est inférieur à
une valeur critique pour laquelle se produit une transition assez brutale vers le turbulent. 2300
est la valeur généralement retenue pour cette transition mais, dans des conditions soignées
(paroi particulièrement lisse, stabilité de la vitesse), la transition peut se produire pour une
valeur plus élevée. On considère souvent que la transition peut se produire entre 2000 et
3000.
• Sur un cylindre à section circulaire placé dans un écoulement, on obtient un écoulement
proprement laminaire qui s’ajuste parfaitement à l’obstacle jusqu’à un nombre de Reynolds
de l’ordre de 1. Un sillage turbulent apparaît à l’aval aux environs de 105. Entre les deux, la
transition se fait à travers diverses formes de sillages tourbillonnaires.
• Avec une plaque plane située dans le lit de l’écoulement, la dimension caractéristique n’est
plus l’épaisseur de celle-ci mais la distance d’un point au bord d’attaque. En effet une couche
limite, dans laquelle interviennent la viscosité ou la turbulence, se développe à partir du bord
d’attaque. Si celui-ci présente une arête émoussée, la couche limite est turbulente dès le dé-
but. Dans le cas d’un bord effilé, la couche limite est laminaire sur une certaine longueur, plus
turbulente ensuite. Ce caractère laminaire se maintient jusqu’à la distance qui correspond au
Reynolds critique de l’ordre de 5.105 , la zone située au delà développant une couche limite
turbulente.
• Pour un profil d’aile, la distribution d’épaisseur le long de la corde (et le gradient de pression
négative associé) de certains profils dits "laminaires" stabilise la laminarité et permet de
reculer le point de transition bien au delà de 5.105 : des valeurs de 7.106 sont possibles dans
des conditions aérologiques non turbulentes (difficiles à obtenir en soufflerie) sur une surface
parfaitement lisse (ailes de planeurs).
• Un corps profilé comme un fuselage peut avoir une transition reculée jusqu’à 5.107 , dans des
conditions idéales également.
v=0 v=0
u
u
u
u
u
u
u v maximale v maximale
u
u
u
u
u
u
5.3.1 Le modèle
—É oulement parfait —
Un écoulement est dit parfait si tous les phénomènes diffusifs peuvent être négligés.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Pour une particule de fluide, cela correspond à négliger les phénomènes diffusifs suivants :
• de quantité de mouvement : absence de viscosité
• de transfert thermique : évolution adiabatique
— Proprié té —
Pour un fluide parfait, la viscosité η est nulle et le nombre de Reynolds Re est infini. Un écoulement de
fluide parfait devrait donc être turbulent. Ce n’est pas le cas car le modèle du fluide parfait n’est pas une
extrapolation du fluide réel de faible viscosité.
En effet, la vitesse d’un fluide réel, même de très faible viscosité, est nulle au contact du solide (conduite
ou obstacle) alors que celle d’un fluide parfait est non nulle et tangente à la surface du solide.
Un fluide transporte de la quantité de mouvement : ce type de transport dû au mouvement est appelé trans-
port convectif.
Le mouvement est également lié à la mise en mouvement de proche en proche des couches de fluide : ce
type de mouvement est diffusif. Le nombre de Reynolds peut également s’écrire :
vd
Re =
ν
Le numérateur correspond au terme convectif (de convection) et le dénominateur correspond au terme
diffusif (de diffusion).
Près d’un obstacle, le terme diffusif est toujours prépondérant sur le terme convectif. Il existe donc une
couche, dite couche limite, où le modèle de l’écoulement parfait n’est pas vérifié.
— Cou he limite —
Sur une petite épaisseur δ autour de l’obstacle se concentrent les phénomènes de viscosité. Au delà,
l’écoulement pourra être considéré comme parfait
• fluide parfait
• écoulement permanent (les trajectoires s’identifient aux lignes de courant)
B B′
Système
de conduites
A A′
• les forces de pesanteur, qui fournissent, avec un axe z ascendant, un travail élémentaire :
δWp = pA DV dt − pB DV dt
Dm
On obtient, avec DV = :
ρ
1 2 2
pA pB
Dm dt vB − vA = −Dm g dt (zB − zA ) + Dm dt − Dm dt
2 ρ ρ
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
On en déduit finalement :
v2
Pour un écoulement incompressible et stationnaire d’un fluide parfait, la quantité p + ρ + ep
2
se conserve le long d’une ligne de courant.
Ainsi, dans le champ de pesanteur, et avec z ascendant, l’énergie potentielle vaut ρgz + Cte et le long d’une
ligne de courant, on a :
p v2
+ + g z = cte
ρ 2
p
• est une énergie massique associée aux forces de pression,
ρ
• gz est l’énergie potentielle massique associée aux forces de pesanteur (dont le poids),
1 2
• v représente l’énergie cinétique massique.
2
• la pression statique p∗ = p + ρ g z,
1
• la pression dynamique v 2 .
2
— Remarques —
p∗ p
• =z+ est appelée hauteur piézométrique,
ρg ρg
1 v2 p 1 v2
• p∗ + =z+ + est appelée charge totale.
2 g ρg 2 g
5.4.3 Applications
La sonde de Pitot permet de mesurer la vitesse d’un écoulement de fluide autour d’une structure, par
exemple de l’air pour un avion. Cette sonde est composée de 2 prises de pression.
À l’entrée du tube en A, la vitesse du fluide est nulle : c’est un point d’arrêt.
M A
b b
p A − p B = ρℓ g h vA ≃ O
B
La vitesse du fluide peut alors en être dé- b
duite : v vB
2 (pA − pB )
v=
ρair h
ou encore :
…
ρℓ
v = 2gh
ρair
L’équation de Bernoulli permet de comprendre le principe de divers instruments et appareils tels que trompe
à eau, bec Bunsen, pistolet à peinture, carburateur, . . .
En revanche, elle ne permet pas d’expliquer les pertes de charge dans les tuyaux, car le fluide est supposé
non visqueux.
Dans cette expérience, un réservoir, qui garde son niveau constant grâce à un système de trop plein se
déverse dans une canalisation sur laquelle sont disposées un certain nombre de prises de pression statique.
[Link] Observations
On constate que :
• le niveau affiché dans les prises de pression est différent,
• ce niveau diminue avec l’éloignement du réservoir,
• cette diminution est proportionnelle à cet éloignement.
⇒ Activité 5. 1
Ces observations sont-elles compatibles avec le théorème de Bernoulli ?
Ceci vient dans la réalité de dissipations par frottements dans l’écoulement : l’égalité de Bernoulli n’est pas
vérifiée. En effet, entre l’amont (noté e pour entrée) et l’aval (noté s pour sortie), on a :
1 1
ps + ρ g zs + ρ vs2 < pe + ρ g ze + ρ ve2
2 2
5.5.2 Définition
— Perte de harge —
Les unité s :
5.5.3 Conséquences
Notons Pi la puissance échangée avec le fluide (Pi > 0 pour les pompes et Pi < 0 pour les turbines).
Bilan d’énergie cinétique :
Considérons comme système, le fluide situé à entre A et B à l’instant t et entre A′ et B ′ à l’instant t + dt :
B B′
Système
de conduites
A A′
Le régime étant permanent, l’énergie cinétique de la portion de fluide située entre A′ et B se conserve. On
obtient alors :
′ ′ 1
dEc = EcBB − EcAA = 2
Dm dt vB 2
− vA
2
• les forces de pesanteur, qui fournissent, avec un axe z ascendant, un travail élémentaire :
δWp = pA DV dt − pB DV dt
δWi = Pi dt
• le travail des forces visqueuses, résistantes et qui fournissent au fluide le travail élémentaire négatif
(puisque résistant) :
Dm
On obtient alors, avec DV = :
ρ
1 2 2
Pi
pB − pA + ρ vB − vA + ρ g (zB − zA ) = − ∆pη
2 DV
∆pη
ou encore, avec Dm = ρDV et JAB = − <0:
ρ
pB − pA 1 2 2
Pi
+ vB − vA + g (zB − zA ) = + JAB
ρ 2 Dm
∆pη
ou même, en notant hη = :
ρg
2 2
pB − pA 1 vB − vA Pi
+ + (zB − zA ) = − hη
ρg 2 g g Dm
Dans les 3 expressions encadrées ci-dessus, ∆pη , JAB et hη représentent les pertes de charge.
La perte de charge JAB peut être due à une perte de charge régulière (ou linéaire), notée JR ou JL ou à une
perte de charge singulière, notée JS .
Quand la conduite subit de brusques variations de section ou de direction, il se produit des pertes de charge
dites singulières.
Celles-ci peuvent être mesurées et calculées à l’aide de l’expression suivante :
v2
JS = −KS
Les pertes de charge régulières ou linéaires sont réparties de façon régulière le long de la conduite. Elles
sont proportionnelles à la longueur L de la conduite, inversement proportionnelles à son diamètre intérieur
d et proportionnelles à la vitesse v du fluide.
En notant λ le cœfficient de perte de charge linéaire (ou régulière) :
v2 L
JL = JR = −λ
2 d
— Proprié té —
En complément des valeurs indiquées page 114, on pourra pour la suite admettre les résultats
empiriques à titre indicatif :
• Si Re < 2000, l’écoulement est laminaire.
• Si Re > 2000, l’écoulement est turbulent :
◦ lisse si 2000 < Re < 105 ,
◦ rugueux si Re > 105 .
On peut calculer de façon approximative les pertes de charge régulières en utilisant les formules suivantes :
64
λ= formule de Poiseuille
Re
…
ε
λ = 0,79 formule de Blench
d
— Remarque —
5.6.1 Énoncés
Dans tous les exercices, la pression atmosphérique sera prise égale à 1,00.105 P a et la masse volu-
mique de l’eau à 1,00 kg.L−1.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Un fluide newtonien est réparti sur une hauteur e entre 2 plaques horizontales de grande dimension. La
plaque inférieure est immobile par rapport à un référentiel galiléen et la plaque supérieure est en mouvement
de translation rectiligne uniforme et horizontale −
→
v = v0 − →
ex .
v0 −
→
ex
1. Principe
(a) Exprimer les conditions aux limites vérifiées par l’écoulement du fluide.
(b) Proposer une formule (la plus simple possible) pour le champ des vitesses vérifiant les condi-
tions précédentes.
(c) Exprimer la composante horizontale de la force par unité de surface exercée par le fluide, sur la
plaque inférieure.
2. Application
Un bloc métallique parallélépipédique, de surface carrée de côté a = 10 cm et de masse m = 1,0 kg,
est posé sur un plan incliné, faisant un angle α = 45o par rapport à l’horizontale. Ce plan incliné est
lubrifié, c’est-à-dire enduit d’une huile de viscosité (dynamique) η. La plaque se met en mouvement.
On suppose que l’écoulement de l’huile peut être modélisé de la même manière qu’en première
partie d’exercice. L’huile possède une épaisseur e = 1,0 mm. Le champ de pesanteur est noté g et
vaut 9,8 m.s−2 .
(a) Déterminer l’équation du mouvement du bloc métallique.
(b) Au bout d’un certain temps, la vitesse du bloc atteint la vitesse limite vf = 5,0 cm.s−1 . En
déduire la viscosité η de l’huile.
(c) Déterminer la durée tt du régime transitoire.
Un tuyau horizontal de section circulaire de 3,0 cm de rayon se rétrécit progressivement jusqu’à 2,0 cm.
Lorsque l’eau s’écoule dans ce tuyau à une certaine vitesse, la pression manométrique au niveau des deux
sections est respectivement 32 kP a et 24 kP a.
Déterminer le débit massique dans le tuyau.
◮ Exercice 5. 5 : Écluse
Dans une écluse, l’eau du compartiment central s’élève à une hauteur de 3,0 m au-dessus du niveau du
passent de l’autre.
Soit A le point d’intersection, supposé fixe et
d’ordonnée zA , entre cette ligne de séparation
et l’obstacle.
À la sortie du cœur humain, l’aorte peut être considérée comme une conduite cylindrique de rayon a0 =
1,0 cm.
Le débit volumique de sang dans l’aorte est DV,0 = 6,0 [Link]−1 et on suppose l’écoulement stationnaire.
La viscosité η du sang est prise égale à 4,0.10−3 P ℓ, sa masse volumique ρ à 1,0 kg.L−1.
À la sortie de l’aorte, se trouvent N artères de rayon a qui se divisent elles-mêmes chacune en N ′ artérioles
de rayon a′ = 20 µm.
Le débit volumique au niveau d’une artère vaut DV,a = 2,0.10−6 m3 .s−1 .
1. En supposant que le champ des vitesses est uniforme sur une section droite, déterminer la vitesse v0
du sang dans l’aorte.
2. Calculer le nombre N d’artères.
Une pompe immergée à 20 m sous terre permet de puiser de l’eau jusqu’à une installation située en surface.
La pression de captage de l’eau est voisine de la pression atmosphérique et le débit nécessaire en surface
est d’environ 7,0 m3 .h−1 . La pression dans l’installation doit valoir environ 4,0 bar.
La section des conduites est constante et l’écoulement est considéré comme stationnaire.
ℓ v2
∆pc = 0,316Re−1/4 ρ
d 2
où ℓ est la longueur de canalisation.
∆pc
Déterminer la perte de charge par unité de longueur de canalisation dans ce cas et vérifier l’ordre
ℓ
∆zc
de grandeur proposé par la règle simple : = 0,10.
ℓ
3. La règle d’usage convient-elle pour un débit compris entre 2,0 m3 .h−1 et 7,0 m3 .h−1 ?
1. En supposant le fluide parfait et en considérant que la relation de Bernoulli est applicable, exprimer
la différence de pression p1 − p2 entre l’amont et l’aval en fonction de l’énergie cinétique volumique
1
du fluide en amont ec,V = ρv 2 .
2
Commenter son signe.
2. Un étude plus précise de cet écoulement permet de mettre en évidence une zone "morte" dans laquelle
le fluide est immobile. Cette zone est située juste derrière l’élargissement.
On montre alors dans ce cas que la différence de pression entre l’amont et l’aval (derrière cette zone
"morte") vaut p1 − p2 = ρv2 (v2 − v1 ) (formule de Bélanger).
En déduire la nouvelle expression de p1 − p2 en fonction de ec,V .
3. Déterminer l’écart relatif entre les résultats des deux questions précédentes.
Commenter cet écart.
4. Proposer une expression pour la perte de charge singulière ∆PC en fonction de ec,V dans le cadre de
cet écoulement.
Commenter le signe de cette perte de charge.
S0
Un capteur destiné à mesurer le débit d’un
fluide est constitué d’une conduite avec un
S1
étranglement. On note S0 et S1 les sections
au niveau des points A0 et A1 situés sur →
− →
−
A0 v1 A1 v2
une ligne de courant moyenne. b b
1. Appliquer la relation de Bernoulli entre les points A0 et A1 de façon à trouver une relation entre les
vitesses v0 et v1 , et les pressions p0 et p1 en ces points.
2. On admet que la formule utilisable dans le cadre de la statique des fluides, s’applique entre les points
A0 et B0 , ainsi qu’entre A1 et B1 .
En déduire une relation entre p0 , p1 et ∆h.
3. Que peut-on dire de la répartition des vitesses sur chacune des sections S0 et S1 ? En déduire une
relation entre le débit volumique et le dénivelé ∆h.
4. En pratique, cette loi n’est pas très bien suivie. Pour quelles raisons ?
5. On retient néanmoins la proportionnalité du débit à ∆hα . Quelle valeur de l’exposant α suggère
l’étude idéale ? Comment obtenir en pratique un capteur exploitable ?
Remarque : On pourra utiliser les formules du cours pour répondre aux questions.
• AB de longueur 6,00 m,
• CD de longueur 10,0 m,
• EF de longueur 5,00 m,
• GH de longueur 4,00 m,
• IJ de longueur 7,00 m,
• KL de longueur 8,00 m.
A B
C
K L
J
E F
D
G
H I
Le débit volumique dans la conduite est DV = 2,50 L.s−1 et la pression en A vaut pA = 3,00 bar.
Remarque : On pourra utiliser les formules du cours pour répondre aux questions.
P6
6.1.1 Énoncé
entrée Σ sortie
δme M (t) δms
x x + dx
À l’instant t, le système est composé du système fermé {entrée + Σ} de masse δme + M (t) et à l’instant
t + dt, il est composé du système fermé {sortie + Σ} de masse δms + M (t + dt).
En régime permanent, δme = δms = δm = Dm dt et MΣ (t) = MΣ (t + dt).
On a, en appelant uΣ l’énergie interne massique du système Σ, ue celle du système à l’entrée et us celle du
système à la sortie :
Soit :
dU
= (us − ue )Dm
dt
On obtient alors :
dEc dEp
Dm [(uΣ + us ) − (uΣ + ue )] + + = Pth + Pi + Ppression
dt dt
En négligeant les variations d’énergies potentielle et cinétique, l’expression précédente se simplifie pour
donner :
Si l’écoulement s’opère de l’entrée e vers la sortie s, en notant ve et vs les volumes massiques respective-
ment à l’entrée et la sortie, le travail élémentaire des forces de pressions est :
Ppression = Dm (pe ve − ps vs )
Soit :
Dm ∆h = Pth + Pi
∆h = wi + qe
— Remarque —
Dans le cas où les variations d’énergies potentielle et cinétique ne sont pas négligées, on a :
6.1.2 Exemples
Considérons un gaz parfait en écoulement stationnaire dans une conduite horizontale, rigide et calorifugée.
Le gaz traverse une zone, par exemple une bourre de coton, destinée à abaisser sa pression.
On en déduit :
hs = he
L’enthalpie d’un gaz parfait ne dépendant que de sa température (deuxième loi de Joule), on en déduit que
la température ne varie pas au cours d’une détente de Joule-Thomson.
Considérons un gaz parfait en écoulement stationnaire dans une tuyère horizontale, rigide, calorifugée et de
symétrie de révolution Oz.
−
→
w −
→
w
e s
O z
1 1
hs + ωs2 = he + ωe2
2 2
qui donne, en utilisant le modèle du gaz parfait :
ωs2 − ωe2
= −(hs − he ) = −cp (Ts − Te )
2
nRγ Rγ γr
Soit, avec cp = = = :
m (γ − 1) M (γ − 1) γ−1
Application numérique : ωe = 300 m.s−1 , ωs = 500 m.s−1, M = 29 [Link]−1, γ = 1,4, R = 8,314 J.K −[Link]−1
Ts − Te ≃ −80 K.
Une turbine comporte une partie tournante formée d’aubes mobiles reliées à un arbre. Les chocs des parti-
cules de fluide sur les aubes mettent ces dernières en mouvement et la vitesse du fluide diminue donc lors
de ce choc.
On cherche à récupérer une puissance mécanique sur l’arbre. Celle-ci est l’opposée de celle reçue par le
fluide :
Parbre = −Pi
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
En négligeant les transferts thermiques ainsi que la variation d’énergie potentielle, on obtient :
ωs2 ω2
hs + − he − e = wi
2 2
et :
ïÅ ã Å ãò
ω2 ω2
Parbre = −Pi = −Dm wi = Dm he + e − hs + s
2 2
Considérons un échangeur thermique calorifugé dans lequel circulent, sans se mélanger et en sens inverses,
2 fluides, l’un entre l’entrée e et la sortie s, l’autre entre l’entrée e′ et la sortie s′ .
s′ e′
e s
Si les fluides sont des liquides, on peut poser, pour ces phases condensées :
Dans les systèmes industriels, sont également utilisés d’autres éléments tels que :
• les compresseurs,
• les séparateurs,
• les mélangeurs,
• ...
entrée Σ sortie
δme M (t) δms
x x + dx
À l’instant t, le système est toujours composé du système fermé {entrée + Σ} de masse δme + M (t) et à
l’instant t + dt, il est composé du système fermé {sortie + Σ} de masse δms + M (t + dt).
On obtient alors, en appelant sΣ l’entropie massique du système Σ, se celle du système à l’entrée et ss celle
du système à la sortie, la variation élémentaire d’entropie :
Le deuxième principe donne pour cette variation d’entropie du système fermé considéré :
δQe
dS = + δScréée
Te
où δQe représente le transfert thermique échangé (reçu, en fait), Te la température extérieure, supposée
uniforme à la frontière du système, et δScréée l’entropie élémentaire créée à l’intérieur du système.
On en déduit :
s δQe
dSΣ (t) + [s δm]e = + δScréée
Te
En régime stationnaire (ou permanent), on a δme = δms = δm = Dm dt et MΣ (t) = MΣ (t + dt).
De plus, en régime stationnaire, l’entropie ne dépend pas du temps, ce qui donne dS = 0.
On obtient alors :
δQe s
+ δScréée = δm (s)e
Te
— Remarque —
L’entropie créée par unité de temps est obtenue en divisant l’expression précédente par dt :
δScréée s Pth
= Dm [s]e −
dt Te
Le terme de droite est calculable par les méthodes habituelles, le terme de gauche en est déduit.
6.2.2 Exemples
Considérons un gaz parfait en écoulement stationnaire dans une conduite horizontale, rigide et calorifugée.
Le gaz traverse une zone, par exemple une bourre de coton, destinée à abaisser sa pression.
On a alors :
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
s
δm [s]e = δScréée
s
Pour calculer la variation d’entropie massique [s]e , choisissons un chemin réversible (en plusieurs étapes
éventuellement) entre les états (pe , Te ) et (ps , Ts ). On peut écrire :
dT dp
dS = n Cp,m −V
T T
dp
= −n R
p
On en déduit :
Å ã
pe
Ss − Se = n R ln
ps
R
Et en posant r = , on obtient :
M
Å ã
R pe
ss − se = ln
M ps
Å ã
pe
= r ln
ps
— Remarques —
dp
ds = −v >0
T
Or :
δqe
ds = + δscréée = δscréée > 0
Te
On vérifie qu’il y a création d’entropie et que la détente est irréversible (c’est une conséquence de la
viscosité).
Å ã Å ã
6.3.1 Énoncés
Déterminer la vitesse maximale d’éjection d’un gaz parfait diatomique entrant dans une tuyère à la pression
p= 1,50 bar et à la température Te = 900 K. Le gaz sort à la pression ps = 1,00 bar. L’écoulement est
considéré adiabatique et réversible.
Donnée :
cp = 1,00.103 J.K −1 .kg −1
Un lac de retenue d’eau est connecté à une centrale hydroélectrique par l’intermédiaire d’une galerie sou-
terraine.
On note H la différence d’altitude entre le point de captage C de l’eau dans le lac et la station où se trouve
la turbine entraînée par la chute de l’eau.
On note c la capacité thermique massique de l’eau liquide, que l’on suppose indépendante de la température.
La masse volumique ρ de l’eau est égale à 1,00 kg.L−1 et l’accélération de la pesanteur g est prise égale à
9,80 m.s−2.
1. Rappeler l’expression de l’enthalpie massique et de l’entropie massique pour une phase condensée.
2. Qu’implique, pour le débit volumique, le fait que l’écoulement soit stationnaire ?
3. Entre le point de captage C et l’entrée E de la turbine, écrire le bilan enthalpique.
4. Entre le point de captage C et la sortie S de la turbine, écrire le bilan enthalpique. On considérera
que les évolutions sont suffisamment rapides pour pouvoir négliger les échanges thermiques avec les
parois.
5. Écrire le bilan entropique entre les points C et S et en déduire une inégalité entre les températures en
ces points.
En déduire une inégalité faisant apparaître le travail massique électrique wélec produit par la centrale
et la différence d’altitude H.
6. Exprimer finalement la puissance maximale par unité de longueur de dénivelé obtenue pour un débit
volumique DV égal à 1,00 m3.s−1 .
Un fluide (de l’eau, de masse volumique ρ) circule dans une turbomachine et décrit les cycles suivants :
• Une pompe A amène le liquide saturant, pris à la sortie du condenseur (état F ), jusqu’à la pression
p1 de la chaudière. Cette opération est pratiquement adiabatique et on peut considérer qu’à la sortie
de la pompe, le fluide est liquide (état G) à la température T2 du condenseur.
• L’eau est alors injectée dans une chaudière C où elle se vaporise de façon isobare à la pression p1 . À
la sortie de la chaudière, la vapeur est saturante sèche à la température T1 (état D).
• Elle subit ensuite une détente adiabatique et réversible dans une turbine T . À la sortie de la turbine, le
fluide est à la température T2 et à la pression p2 d’un condenseur C ′ (point E) dans lequel il achève
de se liquéfier de façon isobare (point F ).
Chaudière
Condenseur
C′
Données :
• T1 = 523 K, T2 = 293 K,
• enthalpie de vaporisation à 523 K : ℓ1 = 1714 [Link] −1,
• pression de vapeur saturante à 523 K : p1 = 39,7.105 P a,
• pression de vapeur saturante à 293 K : p2 = 2300 P a,
• enthalpie massique du liquide saturant à 293 K : hℓ = 84,0 [Link] −1,
• enthalpie massique de la vapeur saturante sèche à 293 K : hv = 2538 [Link] −1,
• capacité thermique massique du liquide : cliq = 4180 J.K −[Link] −1 ,
• volume massique du liquide : vliq = 1,00.10−3 m3 .kg −1 .
P7
Conduction thermique
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
La présence, dans un milieu matériel sans mouvement macroscopique, d’une inhomogénéité de température,
fait apparaître un transfert thermique par conduction.
— Proprié té s —
• Ce transfert a lieu des zones les plus chaudes vers les plus froides.
• Il est proportionnel à la surface à travers laquelle on évalue la puissance diffusée.
• Il croît linéairement avec la pente de la variation spatiale de température.
Les unité s :
−→
jth s’exprime en W.m−2 .
où λ est une constante caractéristique du milieu dans lequel a lieu la conduction thermique. Cette constante
est appelée conductivité thermique.
λ s’exprime en W.m−1 .K −1 .
Le signe négatif signifie que les transferts thermiques ont lieu dans le sens des températures décroissantes.
−→ − →
jth 6= 0 dès que T n’est pas uniforme.
Électrique Thermique
−
→ −−→ −→ −−→
je = −γ grad V jth = −λ gradT
ZZ ZZ
→ −−
− → −→ −−→
I= je · d2 S φth = Pth = jth · d2 S
Σ Σ
V T
γ λ
I Pth
T0 T1
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
0 ℓ
Après ce régime transitoire apparaît un régime stationnaire dans lequel la température dépend de x :
T = T (x).
En régime stationnaire, la quantité de chaleur qui entre est égale à la quantité de chaleur qui sort, autrement
dit, la quantité de chaleur reçue est nulle. Les transferts thermiques se faisant de manière isobare, on peut
écrire :
Qp = ∆H = 0
Isolons une tranche élémentaire de tige de section S située entre les abscisses x et x + dx :
x
x x + dx
Å ã
dT
Le développement de Taylor-Young de donne :
dx x+dx
Å ã Å ã Å ã
dT dT ∂2T dT d2 T
= + dx = + dx
dx x+dx dx x ∂x2 dx x dx2
De plus :
−→ −−→ dT −→
jt h = −λ gradT = −λ ex
dx
−−−−→ −−−−−−→
Avec d2 S(x) = d2 S −
→
ex et d2 S(x+dx) = −d2 S →
−
ex , la puissance thermique totale reçue est alors égale à la
somme des flux thermiques sortants :
Pth = φth
ZZ ZZ
−
→ −−−−→ −
→ −−−−−−→
= jth (x) · d2 S(x) + jth (x + dx) · d2 S(x+dx)
Σ(x) Σ(x+dx)
Å ã Å ã
dT dT
= −λ S + λS
dx x dx x+dx
d2 T
= λ S 2 dx
dx
Cette puissance totale reçue Pth étant nulle en régime stationnaire, on a nécessairement :
T0
T1
x
O ℓ
Le flux thermique conductif traversant toute section de la tige peut alors être exprimé par :
Å ã
dT λS
Pth = Φ = −λS = (T0 − T1 )
dx ℓ
On peut alors écrire une relation analogue à la loi d’Ohm (VP − VQ = Rélec I) :
en posant :
ℓ
Rth =
λS
Φth est du signe de T0 − T1 : cela signifie que si l’extrémité gauche de la barre possède une température
plus élevée que l’extrémité droite, alors le transfert thermique est dirigé de la gauche vers la droite.
Par analogie avec l’électrocinétique, la grandeur Rth est appelée résistance thermique. Cette dernière dé-
pend de la géométrie et de la nature du matériau.
Les unité s :
⇒ Activité 7. 1
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Rth,2
T0 T1
T0 T1 T2
F IGURE 7.2 – Résistances thermiques en série F IGURE 7.3 – Résistances thermiques en dérivation
— Remarque —
On peut représenter les résistances thermiques de la même façon que les résistances électriques :
Rth,2
Φth,1
Rth,1 Rth,2
Φth b
Φth b
b b b
T0 T1 T2 T0 Rth,1 T1
Φth,2
x
x x + dx
De même, le flux thermique sortant par la face droite, située à l’abscisse x + dx, vaut :
ZZ
−
→ −−→
Φth (x + dx,t) = jth (x + dx) · d2 S = jth (x + dx, t) S
Σ(x+dx)
Le développement de Taylor-Young pour une fonction f (x, t) à deux variables, s’écrit, au premier
ordre :
∂f (a, b) ∂f (a,b)
f (a + h, b + k) = f (a, b) + h+ k
∂x ∂y
∂T ∂T
D’autre part, avec dT = dx + dt, le transfert thermique reçu par le tronçon de barre peut s’écrire,
∂x ∂t Å ã
∂T
la température étant uniforme sur le tronçon de barre =0 :
∂x
∂T
δ 2 Q = dm cp dT = ρ S dx cp dt
∂t
δ2Q
On en déduit, en égalisant Φ(x,t) − Φ(x + dx,t) et :
dt
∂jth (x) ∂T
− S dx = ρ S dx cp
∂x ∂t
∂jth (x) ∂T
+ ρ cp =0
∂x ∂t
— Remarque —
−
→ −−→
La température ne dépendant spatialement que de l’abscisse x, la loi de Fourier jth (x) = −λ grad T s’écrit :
∂T
jth (x) = −λ
∂x
λ ∂2T ∂T
− =0
ρ cp ∂x2 ∂t
Ce n’est pas tout-à-fait une équation de propagation, comme on a pu rencontrer en première année, car la
dérivée temporelle est d’ordre 1 uniquement.
Les unité s :
Pour résoudre cette équation de diffusion, il est nécessaire de prendre en compte des valeurs particulières.
Ainsi, les conditions aux limites, par exemple en x = 0 et x = ℓ.
On pourra également s’intéresser aux conditions initiales, à l’instant t = 0.
Soit une paroi solide en contact avec un fluide. Celui-ci est le siège de phénomènes de convection : les
particules de fluide sont en mouvement et le fluide est donc animé de mouvements de convection provoquant
une homogénéisation de sa température.
Paroi Fluide
x
O
Considèrerons que la température du fluide est constante égale à TF et que le phénomène est unidimension-
nel (suivant x).
Dans la paroi, le transfert conductif est donc dirigé suivant l’axe des abscisses x. La température de la paroi
en x = 0− vaut TP .
Comme la paroi et le fluide ne sont pas à la même température, on peut proposer le profil de température
suivant :
T (x)
TP
TF
x
O
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
T (x)
TP
TF
x
O e
Pour une section d’aire S, on peut définir un flux thermique conducto-convectif, noté Φcc correspondant au
transfert thermique entre la paroi et le fluide :
Å ã
dT
Φcc = −λF S
dx 06x6e
En considérant que la température est une fonction affine de l’abscisse sur l’intervalle [0,e], on peut écrire :
Å ã
TF − TP
Φcc = −λF S
e
λF
Plus simplement, en posant h = , on obtient la loi de Newton :
e
ϕcc = h (TP − TF )
⇒ Activité 7. 2
Une ailette de refroidissement cylindrique de rayon r très
petit devant sa longueur ℓ, de conductivité thermique λ est en
contact par une de ses bases avec un thermostat, ici de forme
parallélépipédique, à la température T0 .
On considère le reste de l’ailette de refroidissement en contact
avec l’air ambiant à la température Tf .
On note h le coefficient de conducto-convection de l’ailette et
on rappelle la loi de Newton ϕcc = h (T (x) − Tf ) où ϕcc est le
x
x x + dx
ZZ
−
→ −−→
Φth (x, t) = jth (x) · d2 S = jth (x, t) S
(Σ(x))
De même, le flux thermique sortant par la face droite, située à l’abscisse x + dx, vaut :
ZZ
→
− −−→
Φth (x + dx, t) = jth (x + dx) · d2 S = jth (x + dx, t) S
(Σ(x+dx))
Le développement de Taylor-Young pour une fonction f (x, t) à deux variables, s’écrit, au premier
ordre :
∂f (a, b) ∂f (a, b)
f (a + h, b + k) = f (a, b) + h+ k
∂x ∂y
∂H
dH = H(t + dt) − H(t) = dt
∂t
∂T
δ 2 Qp = d2 H = dm cp dT = ρ S dx cp dt
∂t
δ 2 Qp
On en déduit, en égalisant Φ(x, t) − Φ(x + dx, t) et :
dt
∂jth (x) ∂T
− S dx = ρ S dx cp
∂x ∂t
∂jth (x) ∂T
+ ρ cp =0
∂x ∂t
— Remarque —
En régime stationnaire, on a :
−
→ −−→
La température ne dépendant spatialement que de l’abscisse x, la loi de Fourier jth (x) = −λ grad T s’écrit :
∂T
jth (x) = −λ
∂x
λ ∂2T ∂T
2
− =0
ρ cp ∂x ∂t
Ce n’est pas tout-à-fait une équation de propagation, comme on a pu rencontrer en première année, car la
dérivée temporelle est d’ordre 1 uniquement.
λ
On peut introduire un nouveau paramètre κ = , appelé diffusivité thermique ou cœfficient de
ρ cp
diffusion.
∂2T 1 ∂T
− =0
∂x2 κ ∂t
Les unité s :
Pour résoudre cette équation de diffusion, il est nécessaire de prendre en compte des valeurs particulières.
7.5.1 Dispositif
Résoudre l’équation de la diffusion thermique dans un cas non stationnaire consiste à déterminer l’évolution
de la température en fonction de l’espace et du temps. Si on se place toujours dans un cas unidimensionnel,
par exemple suivant la direction x, la température est alors de la forme T (x, t).
Nous considérons ici une barre cylindrique homogène de longueur ℓ et de section S, entièrement calorifu-
gée :
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
x x
0 ℓ x x + dx
−
→ −−→
• Les conditions aux limites en x = 0 et x = ℓ donnent, grâce à la loi de Fourier jth = −λ grad T :
∂T
(0, t) = 0
∂x
∂T (ℓ, t) = 0
∂x
• Les conditions initiales sont supposées connues. On notera pour la suite T0 (x) le profil de température
à l’instant initial.
En supposant par exemple que la barre était auparavant en régime stationnaire, on peut poser :
x
T (x,0) = T0 + (T1 − T0 )
ℓ
avec le profil de température déjà rencontré auparavant :
T0
T1
x
0 ℓ
Les échanges thermiques au sein de la barre vont mener à une homogénéisation de la température au sein
de celle-ci. Notons la température finale T∞ .
Un bilan enthalpique permet d’écrire :
Z x=ℓ Z x=ℓ
d2 Qp = d2 H
x=0 x=0
Z x=ℓ Z T∞
= dm cp dT
x=0 T0
Z x=ℓ Z T∞
= ρ S cp dT dx
Comme la barre n’échange pas de chaleur avec l’extérieur, cette grandeur est nulle et on en déduit alors :
Z ℓ
1
T∞ = T0 (x) dx
ℓ 0
Cette température correspond à la valeur moyenne de la température dans la barre à l’état initial.
Le régime n’est pas permanent et nous chercherons une solution (dite à variables séparées) de la forme :
∂2T 1 ∂T
− =0
∂x2 κ ∂t
donne alors :
∂2Θ 1 ∂Θ
=
∂x2 κ ∂t
1
f ′′ (x) g(t) = f (x) g ′ (t)
κ
ou encore :
f ′′ (x) 1 g ′ (t)
=
f (x) κ g(t)
La seule solution possible pour cette équation, si on veut qu’elle soit vérifiée à tout instant et n’importe où,
est :
f ′′ (x) 1 g ′ (t)
= = Cte = K
f (x) κ g(t)
( √ √
A −K − B −K = 0
√ Ä√ ä √ Ä √ ä
A −K exp −K ℓ − B −K exp − −K ℓ = 0
( √
−B K=0
√ Ä√ ä √ Ä√ ä
−A K sin K ℓ + B K cos Kℓ =0
Soit :
(
B=0
√ Ä√ ä
A K sin Kℓ =0
et donc que :
√ √ nπ n π 2
K ℓ = n π, soit K = et K = .
ℓ ℓ
Comme dans le cas de la cavité électromagnétique, l’expression de f (x) est quantifiée et vaut cette
fois :
x
fn (x) = An cos n π
ℓ
f ′′ (x) n π 2
Le rapport est alors égal à − et :
f (x) ℓ
g ′ (t) n π 2
= −κ
g(t) ℓ
ℓ2 τ
En posant τ = et τn = 2 , on a :
κ π2 n
g ′ (t) n2 1
=− =−
g(t) τ τn
⇒ Activité 7. 3
Donner la solution gn (t) de cette équation différentielle.
La solution générale de l’équation de la chaleur est alors une combinaison linéaire de tous les modes, et en
posant Θn = An Bn , on obtient :
Å ã
X x n π 2
Θ(x, t) = Θn cos n π exp − κt
n
ℓ ℓ
Θ
p=1 p=3
x
0
ℓ
p=2
Il nous reste ensuite à exploiter les conditions initiales, c’est-à-dire Θ0 (x) = T0 (x) − T∞ .
À l’instant initial t = 0, on a :
X x
Θ0 (x) = Θn cos n π
n
ℓ
Z
2 ℓ
x
Θn = Θ0 (x) cos n π dx
ℓ 0 ℓ
• Dans le temps, la fonction exponentielle montre bien une uniformisation.
• Dans l’espace, le mode fondamental va dominer dès que les harmoniques seront suffisamment atté-
nués : l’évolution devient rapidement sinusoïdale.
T0 T1
x
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
0 ℓ
x
x x + dx
d2 S = δ 2 S échangée + δ 2 S créée
L’entropie créée δ 2 S créée est positive ou nulle.
Avec δQ = Φth dt, l’entropie échangée vaut :
Φth (x) Φth (x + dx)
δ 2 S échangée = dt + dt
T (x) T (x + dx)
Or, on a :
Φth (x) = jth (x) S et Φth (x + dx) = −jth (x + dx) S
On a alors :
Å ã
jth (x) jth (x + dx) ∂ jth (x)
δ 2 S échangée = S dt − S dt = −S dt dx
T (x) T (x + dx) ∂x T (x)
En régime stationnaire, nous avons déjà établi que :
x
T (x) = T0 + (T1 − T0 )
ℓ
et que :
dT
jth (x) = −λ = Cte
dx
d2 T
On en déduit, avec = 0 ou avec jth = Cte :
dx2
Å ãÅ ã
dT 1 dT
δ 2 S échangée = −S dt −λ − 2 dx
dx T (x) dx
Å ã2 Å ã
dT 1
= −λ S dt dx
dx T 2 (x)
Å ã
1 dT 2
= −λ S dx dt
T (x) dx
Un carré étant toujours positif, cela signifie que l’entropie volumique créée est nécessairement po-
sitive. Le phénomène de conduction thermique s’accompagne donc toujours d’une création d’entro-
pie : c’est un phénomène irréversible.
7.7.1 Énoncés
1. Calculer le flux thermique à travers un simple vitrage dont la surface vaut S = 2,00 m2, d’épaisseur
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Données :
Conductivités thermiques (en W.m−1 .K −1 ) du verre et du béton :
λverre = 1,20, λbéton = 1,40.
1. Calculer la résistance thermique d’un mur de briques dont la surface vaut S = 15,0 m2 et d’épaisseur
e = 10,0 cm.
2. Calculer la résistance thermique d’un mur de même aire, composé de 10,0 cm de brique, 30,0 cm de
parpaing et 2,00 cm de plâtre.
Données
Conductivités thermiques (en W.m−1 .K −1 ) :
λbrique = 0,67, λplâtre = 0,80, λparpaing = 1,15.
1. Si la surface du mur est de 30,00 m2, calculer la perte calorifique par conduction pendant une journée
entière.
2. Déterminer la température Tm au milieu du mur.
1. Calculer les pertes thermiques par mètre linéaire sachant que la température T3 à la surface de l’isolant
est de 20,0 o C.
2. On utilise cette conduite, d’une longueur L = 100 m, pour transporter de l’eau chaude (de masse vo-
lumique ρ = 1,00 kg.m−3) dont le débit volumique DV est de 1000 L.h−1. La température d’entrée
Te de l’eau est de 86,0 oC et on désire qu’elle sorte à Ts = 84,0 oC.
La capacité thermique massique à pression constante de l’eau vaut cp = 4,185 [Link] −1.K −1 .
Données :
a = 0,40 mm,
λ = 200 W.m−1.K −1 (aluminium),
h = 100 W.K −1.m−2 .
On assimile localement la croute terrestre à un demi-espace (z > 0) constitué d’un milieu homogène qui
possède, comme caractéristiques :
T (0, t) = T0 + θ0 cos(ωt)
On considère le régime forcé sinusoïdal et on cherche la température, pour z > 0, sous la forme :
T (z, t) = T0 + α(z, t)
On associe alors à α(z, t) la représentation complexe α(z, t) que l’on propose d’écrire ainsi :
x x
0 ℓ x x + dx
On considère un barreau cylindrique de rayon a et de longueur ℓ calorifugé sur sa paroi latérale. On met en
contact une de ses extrémités (placée en x = 0) avec un thermostat de température T0 et l’autre extrémité
(placée en x = ℓ) avec un thermostat de température T1 . On suppose le régime permanent établi.
Le barreau possède une conductivité thermique λ.
P 10 Oscillateurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 199
10.1 Oscillateur quasi-sinusoïdal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 199
10.1.1 Le principe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 199
10.1.2 Démarrage et entretien des oscillations . . . . . . . . . . . . . . . . . . . . . . . . . . 200
10.2 Oscillateur à pont de Wien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 200
10.2.1 Le filtre de Wien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 200
10.2.2 Système bouclé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 201
10.3 Oscillateur à résistance négative . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
10.3.1 Pourquoi une résistance négative ? . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
10.3.2 Réalisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 203
10.3.3 Étude théorique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
[Link] Régime linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
[Link] Régime saturé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 205
[Link] Caractéristique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 205
[Link] Équation différentielle du circuit . . . . . . . . . . . . . . . . . . . . . . . . . . . . 205
10.3.4 Conditions d’oscillations et différents régimes . . . . . . . . . . . . . . . . . . . . . . 206
[Link] Démarrage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 206
[Link] Entretien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 208
[Link] Extinction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 208
10.4 Portrait de phase . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 209
10.4.1 Oscillateur linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 209
10.4.2 Oscillateur à pont de Wien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 210
10.4.3 Oscillateur à résistance négative . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 210
10.5 Oscillateur à relaxation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 210
10.5.1 Réalisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 210
10.5.2 Analyse qualitative . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 211
10.5.3 Calcul de la période . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 212
10.5.4 Allure des tensions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 212
10.5.5 Rapport cyclique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 213
10.6 Analyse spectrale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 213
10.7 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 214
10.8 Les programmes python . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 215
10.9 Exercices : oscillateurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
10.9.1 Énoncés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 220
CPGE TSI
Saint Joseph - LaSalle
166
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
167
P8
L’amplificateur linéaire intégré (ALI) est aussi appelé amplificateur opérationnel ou ampli-op (AOP).
Le modèle idéal a déjà été étudié en Sciences Industrielles en première année. Voici quelques rappels.
+ −
∞ ∞
− +
ve vs vs
ve
⇒ Activité 8. 1
Exprimer vs en fonction de ve .
R2
R1 −
∞
+
ve vs
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
⇒ Activité 8. 2
vs
Exprimer la fonction de transfert H = .
ve
+
∞
−
ve R2
vs
R1
⇒ Activité 8. 3
Exprimer vs en fonction de ve .
R1 −
∞
+
ve1 R2
vs
ve2
⇒ Activité 8. 4
Exprimer vs en fonction de ve .
R
C
ve vs
⇒ Activité 8. 5
vs
Exprimer la fonction de transfert H = .
ve
R −
∞
+
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
ve vs
⇒ Activité 8. 6
Exprimer vs en fonction de ve .
R2
R1
− ✄∞
+
e(t)
s(t)
⇒ Activité 8. 7
vs
Exprimer la fonction de transfert H = du montage et en déduire une équation différentielle entre vs et
ve
ve .
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
P9
9.1 Modèle
L’amplificateur linéaire intégré,noté ALI, est aussi appelé amplificateur opérationnel (ou "ampli-op"), noté
AOP.
−Vcc
v− − ✄∞
vs
v+ + v− +Vcc
+Vcc v+ vs
−Vcc
Par la suite, bien que l’alimentation (+Vcc , −Vcc ) de l’amplificateur linéaire intégré soit indispensable à son
bon fonctionnement, elle ne sera pas indiquée dans le montage.
— Remarque —
Bien que ce ne soit pas indispensable, l’alimentation de l’ALI est en principe symétrique.
−
+
v− vs
v+
À l’entrée inverseuse −, la tension est notée v− et à l’entrée non inverseuse +, elle est notée v+ .
Ces tensions sont bien sûr les différences de potentiel entre les entrées inverseuse/non inverseuse et la masse.
ε = v+ − v−
9.2.1 Propriétés
La différence de tensions ε est très faible et le gain différentiel Ad est donc très grand.
Dans le modèle de l’amplificateur idéal, ce gain différentiel Ad est considéré comme infini.
− vs
∞
Les intensités des courants entrant par les bornes inverseuse et non inverseuse sont négligeables devant les
autres courants. On considère qu’ils sont nuls :
i+ = i− = 0
Les courants de polarisation étant nuls, on dit que la résistance — ou mieux, l’impédance — d’entrée de
l’amplificateur est infinie.
La valeur de la tension de sortie Vs est indépendante du courant débité en sortie de l’amplificateur. L’étage de
sortie est alors assimilable à une source idéale de tension. On dit que la résistance — ou mieux, l’impédance
— de sortie est nulle.
La tension de sortie vs et la différence de potentiel ε sont proportionnelles dans une plage de linéarité :
vs = Ad ε
Vsat
Cette relation est valable tant que |ε| < .
Ad
En effet, la tension de sortie ne peut pas dépasser la plage définie par les alimentation +Vcc et −Vcc du
circuit.
Il y a alors saturation de la tension de sortie lorsque vs atteint la valeur +Vsat ≃ +Vcc .
Vsat
Ainsi, pour ε > , vs = Vsat ≃ Vcc : il y a saturation positive.
Ad
Vsat
Pour ε 6 − , vs = −Vsat ≃ −Vcc : il y a saturation négative.
Ad
Pour l’ALI parfait, on considère que +Vsat et −Vsat sont égales aux tensions d’alimentation. En réalité,
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
pour un ALI réel, la tension de saturation est égale à la tension d’alimentation Vcc moins une tension appelée
tension de déchet ∆V comprise entre quelques dixièmes de volts et quelques volts selon les ALI.
Le courant en sortie de l’ALI ne peut excéder une valeur limite dite de saturation. Celle-ci est de l’ordre de
20 mA.
Dans un amplificateur réel, il subsiste toujours des courants de polarisation i+ et i− de l’ordre de quelques
dizaines de picoampères. Ces courants de polarisation sont à l’origine d’une tension de décalage appelée
également tension d’offset qui subsiste alors entre les deux entrées inverseuse et non inverseuse.
− vs
+ +Vsat
v− vs saturation
v+
régime linéaire ε
saturation
Pas de signe ∞ −Vsat
[Link] Schéma
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
+
∞
−
ve R2
vs
R1
Le courant d’entrée i− étant nul, les courants circulant dans les résistances R1 et R2 sont identiques (on
reconnaît un pont diviseur de tension) :
v− − 0 vs − 0
=
R1 R1 + R2
v− = v+ = ve
et :
R1 + R2
vs = ve
R1
vs R1 + R2
G= =
ve R1
vs
R1 + R2
Soit, avec G = :
R1
Å ã
1 1
ve = vs +
Ad G
et :
vs G Ad Ad
= = 1
ve G + Ad Ad
1+
G
Cette stabilité du gain est due au bouclage opéré sur la borne négative v− .
Imaginons qu’un élément extérieur perturbe la tension de sortie vs (qui par exemple augmente sa valeur).
Comme la rétroaction est faite sur la borne inverseuse, le pont diviseur induira une perturbation de v− dans
le même sens (ici une augmentation).
De cette façon, la tension ε = v+ − v− variera de manière contraire (donc ici, une diminution) et corrigera,
par l’intermédiaire de la relation vs = Ad ε, la tension de sortie vs .
Ce principe est appelé rétroaction ("feed-back" en anglais).
Lorsque la rétroaction est faite sur la borne inverseuse, comme décrit plus haut, la rétroaction est
dite négative et constitue un bon indice de stabilité du montage.
Si elle est faite sur la borne non-inverseuse, elle est positive et peut créer une instabilité du montage,
comme ce sera vu plus tard.
9.4.4 Schéma-bloc
vs = Ad (v+ − v− ) = Aε
v+ = ve
R1 vs
v−
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
= vs = B vs =
R1 + R2 G
s’écrivent encore :
s = Aε
v+ = e
v− = B s = er
e = v+ ε s
+ A
−
er
B
Le comparateur-bloc A est constitué par l’amplificateur linéaire intégré (ALI) et le bloc de gain B représente
le diviseur de tension formé par les résistances R1 et R2 .
Les relations :
®
ε = v+ − v− = e−Bs
s = Aε
donnent :
s = A (e − B s) = A e − A B s
R2
R1 −
∞
+
ve vs R2
vs =H =−
ve R1
+
∞
−
ve R2
vs
R1 + R2
R1
vs = ve
R1
R1 −
∞
+
ve1 R2 R R
vs vs = − ve1 − ve
ve2 R1 R2 2
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
R
C
−
+
∞ dve
vs = −R C
dt
ve
ω
vs H = −j R Cω = −j
ω0
R −
∞
+ Z
1
vs = − ve dt
ve vs RC
L’amplificateur linéaire intégré possède, comme tout système physique, un temps de réponse intrinsèque.
On peut décrire son comportement dynamique en adoptant un modèle du premier ordre.
dvs (t)
τ + vs (t) = Ad ε(t)
dt
τ correspond au temps de réponse intrinsèque de l’ALI.
Typiquement, τ est de l’ordre de 5 ms.
On peut, en première approximation, considérer que les amplificateurs opérationnels réels se comportent
comme des systèmes du premier ordre ayant une fréquence de coupure inférieure voisine de 10 Hz et dont
le produit gain-bande passante en système bouclé est constant.
On peut alors poser le gain différentiel complexe :
vs A0
Ad (jω) = = ω
ε 1+j
ω0
• A0 est le gain différentiel statique, c’est-à-dire celui obtenu pour les signaux constants (noté Ad
précédemment).
1
• ω0 = désigne la pulsation de coupure.
τ
Dans le cas de l’amplificateur non-inverseur, la gain de la chaîne de retour est réel car c’est un rapport de
résistances. On note B(jw) = B sa fonction de transfert.
A0
Ad (jω) = ω
1+j
ω0
vs Ad (j ω)
H(jω) = = cf. 1
ve 1 + Ad (j ω) B
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
et donc :
A0
ω A0
1+j
A0
ω0 A0 1 + A0 B
H(jω) = = ω = ω = ω
A0 1+j + A0 B 1 + A0 B + j 1+j
1+ ω B ω0 ω0 (1 + A0 B) ω0
1+j
ω0
En posant :
A0 1
• un gain statique H0 = ≃
1 + A0 B B
• une pulsation de coupure ωc = ω0 (1 + A0 B)
On peut le mettre sous la forme :
— Gain en dé ibels —
ω H0
On peut poser x = . Dans ce cas, H(jω) = .
ωc 1+jx
Ici, avec G0 = 20 log H0 , le gain en décibels vaut :
H0
G = 20 log √ = G0 − 10 log 1 + x2
1+x 2
H0
H(jω) = ω
1+j
ωc
• Quand ω → 0, H ∼ H0 et GdB ∼ G0 : le gain est maximum.
H0 ωc
• Quand ω → ∞, H ∼ ω = ω H0 et GdB ∼ 20 log H0 + 20 log ωc − 20 log ω, soit une pente à
j
ωc
−20 dB.
H0 √
• Quand ω = ωc , H = √ et GdB = Gmax = 20 log H0 − 20 log 2 = G0 − 3 dB.
2
20 dB bande passante à −3 dB
G0 − 3 dB G0 = 20 log H0
10 dB
0 dB
-20 dB/décade
–10 dB
–20 dB
–30 dB
–40 dB
F IGURE 9.13 – Diagramme de Bode pour l’amplitude
— Dé phasage —
ϕ = arg H
ϕ = −arg [1 + j x] = − arctan x
H0
H(jω) = ω
1+j
ωc
Å ã
ω
ϕ = arg(H) = − arctan
ωc
• Quand ω → 0, ϕ → 0.
π
• Quand ω → ∞, ϕ → − .
2
π
• Quand ω = ωc , ϕ = − .
4
Le diagramme de phase est le suivant :
-π/4
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
-π/2
F IGURE 9.14 – Diagramme de Bode pour la phase
Une propriété fondamentale apparaît, lorsqu’on effectue le produit du gain statique et de la pulsation de
coupure :
H0 ωc = A0 ω0
Le produit du gain statique et de la pulsation de coupure est identique à celui du composant utilisé
pour la réalisation de l’amplificateur.
1
Comme H0 ≃ , on a :
B
G (dB)
20 log A0
1
20 log
B
log ω
log ω0 log ωc
1
En effectuant une analyse temporelle du système bouclé précédent, on peut écrire, en posant τ0 = :
ω0
dvs (t)
τ0 + vs (t) = A0 ε
dt
Cela mène à :
dvs (t)
τ0 + vs (t) (1 + A0 B) = A0 ve
dt
où encore, si 1 + A0 B 6= 0 :
τ0 dvs (t) A0
+ vs (t) = ve
1 + A0 B dt 1 + A0 B
— Remarque —
La constante de temps du système bouclé est égale au rapport de la constante de temps de l’actionneur par
le taux de rétroaction :
τ0
τ′ = ≪ τ0
1 + A0 B
— Exemple —
τ0 dvs (t)
+ vs (t) = 0
1 + A0 B dt
On constate que les cœfficients de l’équation différentielle du régime libre sont tous de même signe : cela
signifie que le système est stable.
τ0
En effet, en posant τ ′ = , la solution de cette équation différentielles est du type :
1 + A0 B
Si le bouclage avait été réalisé sur la borne non inverseuse, l’équation différentielle aurait été :
Avec un produit A0 B très supérieur à 1, les cœfficients de l’équation différentielle auraient été de signes
différents et le système aurait été instable.
τ0
En effet, en posant τ ′′ = , la solution de cette équation différentielles est du type :
1 − A0 B
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Pour obtenir un opérateur stable, il faut effectuer la rétroaction sur l’entrée inverseuse d’un amplifi-
cateur linéaire intégré.
− +
∞ ∞
+ −
E vs E vs
ve ve
Un comparateur à un seuil est parfois insuffisant. Il est alors plus judicieux d’utiliser 2 seuils.
Ainsi, un détecteur de lumière placé dans une zone où la luminosité varie autour de sa valeur seuil (par
exemple en présence de nuages ou d’ombres), passera son temps à déclencher dans un sens et dans l’autre.
Il sera alors plus judicieux d’utiliser un comparateur à 2 seuils, également appelé comparateur à Hystérésis.
En appelant E− et E+ les deux seuils du comparateur, on pourra alors donner son fonctionnement :
⇒ Activité 9. 2
À partir de la description du fonctionnement du comparateur simple non inverseur, représenter la caracté-
ristique entrée - sortie du comparateur à Hystérésis non inverseur.
+Vsat
E− E+ ve
−Vsat
9.7.3 Réalisation
−
∞
+
R2
ve R2
R1 + vs
∞
−
ve R1
vs
⇒ Activité 9. 3
Dans le cas où ve varie entre E− et E+ exclus, on constate que vs ne change pas de valeur. Cela signifie
par exemple que si ve = 0, l’état de vs dépend de la valeur précédente de ve . Le dispositif conserve donc
une mémoire de ses états antérieurs : on parle de fonction mémoire.
Ainsi, en notant l’état 1 haut de vs (+Vsat ) et 0 son état bas (vs = −Vsat ), cet état n’est pas modifié tant
que le seuil de basculement n’est pas franchi : on parle de transition commandée. Cette commande se fait
ici par l’intermédiaire de ve .
Le comportement de filtre du premier ordre de l’amplificateur linéaire intégré est lié à la présence d’un
condensateur au sein du circuit intégré.
En effet, l’étage différentiel fournit à l’étage d’amplification de tension un courant proportionnel à la diffé-
rence de tension entre les deux entrées. Ce courant sert majoritairement à charger la capacité présente dans
dVc
l’étage d’amplification en tension. La relation courant/tension est alors celle d’un condensateur : i = C .
dt
Ainsi, la vitesse de balayage (ou "slew rate") représente la vitesse de variation maximale de tension que
peut produire un amplificateur. Lorsque la vitesse de variation du signal de sortie d’un amplificateur est
supérieure à sa vitesse de balayage, sa tension de sortie est une droite de pente SR .
Å ã
dvs (t)
SR = max
dt
Les unité s :
La vitesse de balayage est de l’ordre de 0,67 V /µs pour un ALI de type A741, et de 15 V /µs pour un
T L081.
Imaginons que la tension de sortie soit en théorie un échelon de tension vs,th . La tension de sortie aura en
réalité l’allure de vs,prat :
vs,th
vs,prat
pente SR
De même, pour un signal sinusoïdal de fréquence très élevée en entrée, le signal de sortie sera limité par la
vitesse de balayage de l’amplificateur linéaire intégré.
Si de plus, l’ALI est utilisé en mode saturé, le signal de sortie ressemblera à la figure suivante :
vs
Vsat
t
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
−Vsat
9.8.1 Énoncés
+
∞
−
ve R2
vs
R1
Le gain statique Ad = A de l’ALI vaut 1,00.105, sa vitesse de balayage SR vaut 1,00.107 V /µs et les
résistances sont R1 = 10,0 kΩ et R2 = 90,0 kΩ. L’ALI est alimenté en ± 15,0 V .
La tension d’entrée ve s’exprime par :
ve = E cos (ω1 t + ϕ1 )
avec E = 0,500 V .
R2
R1 −
∞
+
ve vs
On donne R1 = 33 kΩ et R2 = 69 kΩ.
1. Observations expérimentales
En appliquant un signal sinusoïdal d’amplitude 5,00 V et de fréquence égale à 100 Hz à l’entrée du
montage, on observe l’oscillogamme suivant :
(d) Avec les mêmes signaux, mais en sélectionnant le mode X − Y , on obtient l’oscillogramme
suivant :
Expliquer.
Quelle propriété du montage peut-on en déduire de l’indication des curseurs ∆X = 12,4 V et
∆Y = 27,2 V ?
2. Modélisation par la rétroaction
On souhaite tenir compte du gain différentiel fini, noté A, de l’amplificateur linéaire intégré. On
suppose que l’amplitude du signal d’entrée est suffisamment faible pour qu’aucune saturation n’ap-
paraisse.
(a) On pose ε = v+ − v− . Proposer un schéma fonctionnel composé de 3 blocs linéaires, comme
indiqué ci-dessous.
ve ε vs
C + A
−
vs
(b) Exprimer en régime permanent le rapport d’amplification G = de l’opérateur et examiner le
ve
comportement lorsque A → ∞. Conclure.
3. Étude de la bande passante
(a) Justifier que la bande passante de l’opérateur complet est fixée par celle du réseau bouclé formé
par le soustracteur et les blocs A et B.
(b) En déduire l’expression de la bande passante de l’opérateur dans l’hypothèse d’un amplificateur
linéaire intégré du premier ordre de gain statique A0 et de fréquence de coupure f0 .
(c) Proposer une valeur numérique pour cette largeur de bande passante.
R2
R1 +
∞
−
ve vs
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
P 10
— Os illateur —
La non-linéarité est une idée essentielle pour les oscillateurs : elle engendre des harmoniques et stabilise les
oscillateurs.
— Exemple —
m ẍ + f ẋ + k x = 0
L’équation est non linéaire dès lors qu’un des facteurs m, f ou k est fonction de x.
Il est parfois intéressant d’avoir un système instable, c’est le cas dans les dispositifs oscillateurs.
10.1.1 Le principe
ve vs
Amplificateur
Quadripôle
• le démarrage des oscillations qui constitue un régime transitoire pendant lequel l’amplitude des os-
cillations croît,
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
• un régime permanent pendant lequel le signal de sortie est périodique. Durant cette phase, le signal
n’est pas parfaitement sinusoïdal car des phénomènes de non-linéarité interviennent, par exemple la
saturation d’un ALI.
Considérons le filtre suivant dans lequel on applique un signal vs à l’entrée et pour lequel on récupère
ve à la sortie. Cette configuration peu ordinaire tient au fait que la sortie du filtre est reliée à l’entrée de
l’amplificateur et que la sortie de l’amplificateur est injectée à l’entrée du filtre.
R C
i
vs R C ve
⇒ Activité 10. 1
R2
− ✄∞
S
R1
+
C R
E
vs (t)
ve (t)
C R
R1 + R2
vs = A ve = ve
R1
vs
Soit ve = et :
A
d2 vs 3 − A dvs vs
+ + 2 =0
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
dt2 τ dt τ
Une valeur de 3 pour A permettrait d’obtenir un oscillateur harmonique. Une telle valeur est cependant
risquée car la moindre dispersion énergétique dans le circuit entraînerait un amortissement du signal et
donc, sa disparition.
• Lorsqu’un système obéit à une équation différentielle du deuxième ordre, il est stable si les
3 cœfficients sont de même signe. Ainsi, sa solution ne diverge pas.
• Si un des cœfficients est de signe différent, alors, la solution diverge : choisissons le terme
d’ordre 2 est positif.
◦ Si le terme d’ordre 1 est négatif, cela correspond à un frottement "négatif", donc une
amplification du signal, ce que l’on souhaite ici.
◦ Si le terme d’ordre 0 est négatif, la solution n’est pas oscillatoire (en sin et cos ou
exp(i . . . ) et exp(−i . . . )) mais divergente (en ch et sh ou exp(+ . . . ) et exp(− . . . )).
◦ Si les termes d’ordres 0 et 1 sont négatifs, les deux effets se combinent.
1
La fréquence des oscillations est donnée par ω02 = , soit :
τ2
ω0
f0 =
2π
On voit ainsi que la fréquence des oscillations est fixée par le filtre de Wien :
1
f0 =
2πRC
⇒ Activité 10. 2
Quelles seront les conséquences sur le signal de sortie vs :
1. d’une augmentation de R2 ?
2. d’une diminution de R ?
Si on monte en série une bobine et un condensateur, on doit obtenir en principe en régime libre, pour
l’intensité, une équation du type :
La réponse à un échelon de tension devrait donc être des oscillations harmoniques (sinusoïdales).
Or, en pratique, il y a un amortissement, lié à la résistance du circuit, notamment dans la bobine : le circuit
LC est en réalité un circuit RLC.
C L
Pour obtenir des oscillations quasi-sinusoïdales, il faudrait annuler la résistance du circuit. Ceci est possible
en insérant dans le circuit un dipôle simulant une résistance négative.
10.3.2 Réalisation
R1
ie
−
∞
+
R2
ve vs
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
R0
ε = 0 = v+ − v−
Soit :
v+ = v− = ve
R0
v+ = vs
R0 + R2
De plus :
ve − vs
ie =
R1
Soit :
ve = R1 ie + vs
R0 + R2
ve = R1 ie + ve
R0
R1
et finalement, en posant RN = R0 :
R2
R1
ve = −R0 ie = −RN ie
R2
Cette relation est satisfaite tant que |vs | < Vsat , c’est-à-dire tant que :
ve = R1 ie + vs = R1 ie ± Vsat
Comme on peut le voir avec cette dernière expression, l’ordonnée à l’origine de ve vaut ±Vsat .
[Link] Caractéristique
ie
1
pente
R1 ve
+Ve,max
1
pente −
RN
pente −RN pente R1
+Ve,max ve ie
−Ve,max
pente R1
−Ve,max
1
pente
R1
(a) courant - tension (b) tension - courant
On insère maintenant ce dipôle simulant une résistance négative dans le circuit RLC précédent.
Suivant la tension que l’on désire observer, les dipôles ne seront pas placés dans le même ordre.
R1
R −
∞
+
L
R2
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
uC C R0
u C + u L + u R + u RN = 0
Avec :
q
uC =
C
di
uL = L
dt
uR = Ri
u RN = −RN i
dq duC
Avec i = , on a i = C et :
dt dt
di
uC + L + (R − RN ) i = 0
dt
soit :
q di
+L + (R − RN ) i = 00
C dt
dq
ou encore, avec i = , et en dérivant l’équation précédente :
dt
d2 i R − RN di 1
+ + i=0
dt2 L dt L C
1 R − RN
En posant ω02 = et 2 λ = , on peut également l’écrire :
LC L
d2 i di
2
+2λ + ω02 i = 0
dt dt
[Link] Démarrage
Le système est stable si les 3 termes de l’équation différentielle sont de même signe.
Pour qu’il y ait des oscillations, il faut donc que le terme du premier ordre soit négatif, soit RN > R (donc
Evolution temporelle
0.04
0.03
0.02
0.01
Intensité i (A)
0.00
−0.01
−0.02
−0.03
−0.04
0 5 10 15 20 25 30 35
Temps t (s)
En régime établi, l’amplificateur est saturé et le régime est gouverné par une autre équation différentielle
qui tient compte du fait que l’amplificateur linéaire intégré fonctionne en régime saturé :
uRN ∓ Vsat = R1 i
et :
di
uC + L + (R + R1 ) i = ∓ Vsat
dt
soit :
d2 i (R + R1 ) di 1
2
+ + i=0
dt L dt L C
1 R + R1
ou encore, avec ω02 = et 2 λ′ = :
LC L
d2 i di
+ 2λ′ + ω02 i = 0
dt2 dt
On voit ainsi que dans les deux régimes, transitoire et établi, l’équation différentielle régissant l’évo-
lution de l’intensité i est linéaire à cœfficients constants.
Le caractère non linéaire du régime établi provient uniquement de la variation du signe d’amortis-
sement (2 λ ou 2 λ′ ) suivant que l’amplificateur linéaire intégré fonctionne en régime linéaire ou
saturé.
[Link] Entretien
L’entretien des oscillations ne nécessite pas d’alimentation, comme on peut le voir sur la représentation du
circuit.
Il ne faut cependant pas oublier que l’amplificateur linéaire intégré est alimenté par une source externe de
tension : c’est elle qui fournit l’énergie nécessaire à l’entretien de ces oscillations.
On peut représenter l’entretien des oscillations de la façon suivante :
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Evolution temporelle
0.04
0.03
0.02
0.01
Intensité i (A)
0.00
−0.01
−0.02
−0.03
−0.04
30 32 34 36 38 40
Temps t (s)
[Link] Extinction
Si la valeur de la résistance négative est insuffisante, la résistance du circuit crée un amortissement des
oscillations jusqu’à leur disparition :
Evolution temporelle
0.04
0.03
0.02
0.01
Intensité i (A)
0.00
−0.01
−0.02
−0.03
−0.04
30 35 40 45 50 55 60 65 70
Temps t (s)
• Pour le démarrage, les allures de l’amplitude et du portrait de phase sont les suivantes :
Evolution temporelle Portrait de phase
0.04 0.04
0.03 0.03
0.02 0.02
0.01 0.01
Intensité i (A)
−0.01 −0.01
−0.02 −0.02
−0.03 −0.03
−0.04 −0.04
0 5 10 15 20 25 30 35 −0.06 −0.04 −0.02 0.00 0.02 0.04 0.06
Temps t (s) Intensité i (A)
0.03 0.03
0.02 0.02
0.01 0.01
Intensité i (A)
0.00 0.00
−0.01 −0.01
−0.02 −0.02
−0.03 −0.03
−0.04 −0.04
30 32 34 36 38 40 −0.06 −0.04 −0.02 0.00 0.02 0.04 0.06
Temps t (s) Intensité i (A)
0.03 0.03
0.02 0.02
0.01 0.01
Intensité i (A)
0.00 0.00
−0.01 −0.01
−0.02 −0.02
−0.03 −0.03
−0.04 −0.04
30 35 40 45 50 55 60 65 70 −0.06 −0.04 −0.02 0.00 0.02 0.04 0.06
Temps t (s) Intensité i (A)
— Remarque —
Les graphiques précédents ont été tracés grâce au programme [Link] que vous trouverez en fin de
chapitre.
6
Demarrage des oscillations 6
Portrait de phase réduit
4 4
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
2 2
0 0
−2 −2
−4 −4
−6 −6
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 −5 0 5
temps t (s) Entrée ve (V)
5
Sortie vs (V)
0 0
−5 −5
−10
−10
−15
−15 −20
0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 −30 −20 −10 0 10 20 30
temps t (s) Sortie vs (V)
Mêmes remarques que pour l’oscillateur à pont de Wien. On peut assimiler l’oscillateur à résistance négative
à un oscillateur de Van Der Pol et construire le portrait de phase (obtenus grâce au programme [Link])
suivant :
15 15
Intensité i(t)
10 10
Dérivée de l'intensité di(t)/ω
5 5
Intensité i(t)
0 0
−5 −5
−10 −10
−15
−20 −10 0 10 20
−15
0 2 4 6 8 10 12 14 16
10.5.1 Réalisation
R2
R1
+ ✄∞
− b
B
uC
C R
s1 (t)
s2 (t)
b
− A
✄∞ +
v+ > v− = 0
R1
c’est-à-dire lorsque ve > Vsat = E+ = V0 .
R2
C’est le seuil haut.
R1
E− = − Vsat = −V0
R2
R1
Cela signifie que s2 (t) = +Vsat tant que s1 (t) > − Vsat = −V0 et que s2 (t) = −Vsat tant que
R2
R1
s1 (t) < Vsat = V0 .
R2
Supposons que le comparateur à Hystérésis soit initialement dans l’état s2 (t) = −Vsat .
Vsat
Dans ce cas, ds1 (t) = dt > 0 et s1 (t) croît linéairement jusqu’à atteindre la valeur V0 , valeur pour
RC
laquelle s2 (t) bascule à +Vsat .
Vsat
On a alors ds1 (t) = − dt < 0 : s1 (t) décroît linéairement jusqu’à atteindre la valeur −V0 , valeur pour
RC
laquelle s2 (t) bascule à −Vsat .
uC (0− ) = uC (0+ ) = V0
Vsat
uC (t) = E+ − t
RC
Cet état de saturation positive du comparateur à Hystérésis a lieu tant que uC (t) > −V0 , c’est-à-dire tant
que :
V0 R1
0 < t < t1 = R C = RC
Vsat R2
V0 R1
t1 < t < t2 = 2 R C = 2RC
Vsat R2
R1
T = 4RC
R2
s1 , s2
+Vsat
t1 t2 t1
V0
−V0
−Vsat
Une des améliorations classiques d’un oscillateur à relaxation consiste à rendre les deux états dissymé-
triques : on peut alors les noter t+ (saturation positive) et t− (saturation négative).
— Rapport y lique —
Pratiquement, on utilise fréquemment une association de deux diodes, tête-bêches, placées en dérivation.
⇒ Activité 10. 3
Déterminer le rapport cyclique du montage étudié précédemment.
Dans le cas idéal des deux premiers oscillateurs (à pont de Wien et à résistance négative), l’analyse spectrale
des intensités circulant dans les circuits devrait fournir l’harmonique de rang 1.
En effet, on devrait obtenir l’équation différentielle :
d2 i
+ ω02 i = 0
dt2
10.7 Conclusion
Le filtre de Wien est peu sélectif et on observe en fait, en plus du fondamental, les harmoniques de rang 3
et 5.
L’oscillateur à résistance négative permet d’obtenir des signaux plus purs mais il est coûteux et encombrant.
En pratique, on utilise de préférence des oscillateurs à quartz ou des montages à transistors.
59 [Link]()
60 [Link](y2[:,0],y2[:,1])
61 [Link]("Portrait de phase")
62 [Link](’Intensité i (A)’)
63 [Link](’Dérivée di/dt (A/s)’)
64 [Link]([-0.04,0.04,-0.12,0.12])
65 [Link]("[Link]")
66 [Link]()
67
68 [Link]()
69 [Link](’equal’)
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
70 [Link](y2[:,0],y2[:,1]/[Link](8))
71 [Link]("Portrait de phase")
72 [Link](’Intensité i (A)’)
73 [Link](’1/$\omega$ di/dt (A/s)’)
74 [Link]("[Link]")
75 [Link]()
76
77 [Link]()
78 t=[Link](30,70,1000)
79 y3=odeint(f,[0.035,0],t)
80 [Link](t,y3[:,0])
81 [Link]("Evolution temporelle")
82 [Link](’Temps t (s)’)
83 [Link](’Intensité i (A)’)
84 [Link]([30,70,-0.04,0.04])
85 [Link]("[Link]")
86 [Link]()
87
88 [Link]()
89 [Link](y3[:,0],y3[:,1])
90 [Link]([-0.04,0.04,-0.12,0.12])
91 [Link]("Portrait de phase")
92 [Link](’Intensité i (A)’)
93 [Link](’Dérivée di/dt (A/s)’)
94 [Link]("[Link]")
95 [Link]()
96
97 [Link]()
98 [Link](’equal’)
99 [Link](y3[:,0],y3[:,1]/[Link](8))
100 [Link]("Portrait de phase réduit")
101 [Link](’Intensité i (A)’)
102 [Link](’1/$\omega$ di/dt (A/s)’)
103 [Link]("[Link]")
104 [Link]()
105
106 [Link]()
107 t=[Link](0,70,1000)
108 y4=odeint(f,[0.01,0],t)
109 [Link](t,y4[:,0])
110 [Link]("Evolution temporelle")
111 [Link](’Temps t (s)’)
112 [Link](’Intensité i (A)’)
113 [Link](0,0.01,’ro’)
114 [Link]("[Link]")
115 [Link]()
116
117 [Link]()
118 [Link](y4[:,0],y4[:,1])
119 [Link]([-0.04,0.04,-0.12,0.12])
120 [Link]("Portrait de phase")
Programme [Link] :
1 # -*- coding: utf-8 -*-
2
3 import numpy as np
4 import [Link] as plt # Outils graphiques
5
6 R=1e3 # Resistance du Pont de Wien
7 C=100e-9 # Capacite du Pont de Wien
8 K=3.2 # Gain de l’etage d’amplification
9 tfin=15e-3 # instant final de la simulation
10 dt=0.01e-3 # pas
11 vsat=15 # tension de saturation de l’AO
12
13 # Conditions Initiales
14 t=[0.0]
15 e=[0.01]
16 s=[0.0]
17 ds=[0.0]
18 de=[0.01]
19
20
21 while t[-1]<tfin:
22 [Link]((1/(R*C)*ds[-1]-3/(R*C)*de[-1]-e[-1]/(R*R*C*C))*dt+de[-1])
23 [Link](de[-1]*dt+e[-1])
24 if e[-1]<-vsat/K:
25 [Link](-vsat)
26 elif e[-1]>vsat/K:
27 [Link](vsat)
28 else:
29 [Link](K*e[-1])
30 [Link]((s[-1]-s[-2])/dt)
31 [Link](t[-1]+dt)
32
33 [Link]()
34 [Link](t,e)
35 [Link]("temps t (s)")
36 [Link]("Entrée $v_e$ (V)")
37 [Link]("Demarrage des oscillations")
38 [Link](0,0.01,’ro’)
39 [Link]("[Link]")
40 [Link]()
41
42 [Link]()
43 [Link](e,de)
44 [Link]("Entrée $v_e$ (V)")
45 [Link]("Dérivée d$v_e$/dt (V/s)")
46 [Link]("Portrait de phase")
47 [Link](0.01,0.01,’ro’)
48 [Link]("[Link]")
49 [Link]()
50
51 de_ar=[Link](de,float)
52 de_red=de_ar*(R*C)
53 [Link](’equal’)
54 [Link](e,de_red)
55 [Link]("Entrée $v_e$ (V)")
56 [Link]("RC * d$v_e$/dt (V/s)")
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Programme [Link] :
1 # -*- coding: utf-8 -*-
2
3 import numpy as np
4 import [Link] as plt
5 from [Link] import odeint
6
7 tau = 0.5
8 Q = - 5
9 omega = 2*[Link]/tau
10
11 def f(y,t):
12 if abs(y[0])<5:
13 return [y[1],-omega**2*y[0]-(omega/Q)*y[1]]
14 if abs(y[0])>5:
15 return [y[1],-omega**2*y[0]-(1/tau)*y[1]]
16
17 t=[Link](0,15,10000)
18 # Calcul en partant de y=0.01, y’=0
19 y1=odeint(f,[0.01,0],t)
20
21 [Link](t,y1[:,0],label="Intensité i(t)")
22 [Link](loc=’upper left’)
23 [Link](’Temps (t)’)
10.9.1 Énoncés
− −
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
∞ ∞
+ +
ve R2 ve C R2
vs vs
R1 R1
R2
− ✄∞
S
R1
+
C R
E
vs (t)
ve (t)
C R
1. (a) Déterminer l’équation différentielle vérifiée par vs (t) si l’amplificateur linéaire intégré fonc-
tionne en régime linéaire.
(b) À quelle condition le montage est-il instable ?
2. Pour éviter une croissance trop importante de l’amplitude de vs (t) ainsi que sa saturation, on rend le
vs
rapport µ = dépendant de l’amplitude de ve .
ve
(a) µ doit-il être une fonction croissante ou RL (Ω)
décroissante de l’amplitude de ve ?
On considère un amplificateur linéaire intégré idéal dont la tension de sortie vaut Vsat > 0 :
vs
+Vsat +
saturation (positive) ∞
−
régime linéaire v+
ε v− vs
saturation (négative)
−Vsat
R′
i
B′ +
∞
−
R2
v vs
R1
1. (a) Dans le cas où l’ALI conserve un régime de fonctionnement linéaire, montrer que le dipôle B B ′
est équivalent à une résistance négative de valeur −RN (avec RN > 0) que l’on exprimera en
fonction de R1 , R2 et R′ .
(b) Le régime de fonctionnement de l’ALI n’étant plus supposé linéaire a priori, déterminer la
caractéristique i = f (v) du dipôle B B ′ , lorsque v varie entre −Vsat et +Vsat .
Représenter la courbe i = f (v).
(c) On admettra, dans la suite de l’exercice, que l’on peut remplacer l’équation i = f (v) de la
caractéristique trouvée précédemment par l’équation approchée :
i = g(v) = a v −1 + b v 2
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
B′
L R C
P 11
Ce chapitre fait suite et rappelle quelques notions vues en première année en Sciences Industrielles dans le
chapitre Traitement de l’information. Vous pourrez vous y reporter avec intérêt.
— Signal —
Exemples : tension électrique (en V ), température (en ◦ C), niveau sonore (en dB), pH (sans unité), . . .
— Signal analogique —
U (V)
t(s)
0
U (V)
0 1 0 0 0 1 0 1 0 1 1 0 0 1
0 t(s)
L’enregistrement numérique consiste à convertir le signal électrique en une suite de nombres dont
chacun représente l’amplitude instantanée du signal originel à un instant significatif donné, puis à
enregistrer ces nombres après un codage qui permet de détecter, à la lecture, un défaut éventuel.
• L’enregistrement numérique est beaucoup plus tolérant que son équivalent analogique au niveau du
support d’enregistrement.
• Il permet de rendre indépendant le signal de la distance : lorsqu’un signal analogique est transporté
sur un canal de transmission, il subit de nombreuses modifications, comme l’atténuation ou l’ajout de
bruit, qui affectent la qualité de cette transmission. À l’arrivée, après amplification, le signal originel
est mêlé à du bruit, ce qui dans certains cas, rend difficile la compréhension du message.
Les signaux numériques ne prenant que deux valeurs, "0" ou "1", le bruit occasionné par les canaux
de transmission peut être enlevé de manière simple et efficace. Le signal arrivant est une réplique
exacte du message d’origine, d’où une qualité sans équivalent.
• L’enregistrement numérique permet également de pouvoir facilement multiplexer sur un même canal
de transmission plusieurs signaux de voix, qui sont agrégés sur un même lien physique.
• Il permet de bénéficier des développements et progrès informatiques.
• Il est moins coûteux.
Horloge FE
Les opérations précédentes sont cadencées par une horloge de fréquence FE , où FE correspond à la fré-
quence d’échantillonnage.
11.2 Stroboscopie
11.2.1 Principe
rayons. b b
b b
Lorsque le petit disque est immobile, on est dans la situation b
suivante :
Un stroboscope est un appareil qui émet des flashs lumineux à une fréquence réglable. On supposera pour
la suite que les durées mises en jeu sont grandes devant la persistance rétinienne, de façon à ce que l’obser-
vateur n’ait pas un sentiment de continuité. Ainsi, le mouvement paraîtra saccadé.
b
b b
d’émission du stroboscope.
Considérons la rotation du petit disque dans le sens horaire et une b b
disque. b
b
b
En notant TD la période de rotation et TS l’intervalle de temps entre 2 flashs successifs, on verra alors :
b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b
Si la fréquence d’émission du stroboscope est égale à la fréquence de rotation du disque (ou plus généra-
lement si la fréquence de rotation du disque est un multiple p de celle du stroboscope), alors le disque a
le temps d’effectuer un tour (ou plusieurs) dans l’intervalle de temps entre 2 flashs successifs et celui-ci
b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b
1
F IGURE 11.2 – Observation du disque avec le stroboscope pour fS = fD , p ∈ N
p
De la même façon, si la fréquence du stroboscope est égale à 2 fois celle de rotation du disque, ce dernier
aura le temps d’effectuer un demi-tour entre 2 flashs successifs :
b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b
11 11
Si on prend fD = fS , on a alors TS = TD : entre 2 flashs successifs, le disque n’a pas tout-à-fait le
12 12
temps d’effectuer un tour. On est alors dans la configuration suivante :
b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b
12
F IGURE 11.4 – Observation du disque avec le stroboscope pour fS = fD
11
On voit dans cet exemple que les images observées ne correspondent plus du tout au sens de rotation réelle
du disque.
C’est ce que l’on peut voir pour les roues de voiture dans des films ou dans la rue à la lumière de l’éclairage
public.
11.3 Échantillonage
—É hantillonnage —
L’échantillonnage d’un signal analogique sa (t) consiste à multiplier le signal sa (t) par une fonction
peigne de Dirac constituée d’impulsions très courtes se répétant périodiquement. Cette répétition
est caractérisée par la fréquence d’échantillonnage notée Fe .
Comme dans le cas précédent, on échantillonne donc un signal analogique sa (t) en relevant sa valeur à
intervalles de temps réguliers. La durée entre 2 intervalles successifs est appelé période de l’échantillonnage
et est notée Te . On récupère alors des valeurs du signal sa (t) à chaque instant défini par tn = n Te , avec
n ∈ Z.
On a alors :
11.3.2 Intérêt.
U (V)
t(s)
0
On peut obtenir :
U (V)
t(s)
0
Si le signal est trop faible, on peut même ne plus rien reconnaître du tout :
U (V)
t(s)
0
U (V) U (V)
donne
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
0 t(s) 0 t(s)
1
La fréquence d’échantillonnage Fe = correspond au nombre d’échantillons prélevés par unité de temps.
Te
On peut ainsi faire le parallèle avec la fréquence du stroboscope.
On cherche bien sûr à diminuer autant que possible cette fréquence d’échantillonnage, tout en ne perdant
pas de vue que les échantillons, mis bout à bout, doivent restituer le signal le plus fidèlement possible.
Ainsi, on ne pourra pas descendre en-deçà d’une certaine limite :
U (V)
t(s)
0
Te
Ici, la courbe rouge ne représente pas du tout le signal original. Il faut augmenter la fréquence d’échan-
tillonnage fe et donc diminuer la période Te :
U (V)
t(s)
0
Te
ou mieux :
U (V)
Plus la fréquence d’échantillonnage sera grande, plus la mesure sera précise et plus le signal pourra être res-
titué fidèlement. Mais cela demande un espace mémoire conséquent en raison du grand nombre de données
à traiter.
Considérons un signal sinusoïdal sa (t) de fréquence f = 300 Hz, donc de période T = 3,33 ms, échan-
tillonné à l’aide d’un signal impulsionnel de fréquence Fe = 2,00 kHz.
U (V)
t(s)
0
Pour réaliser de façon pratique un échantillonnage, on multiplie le signal analogique par ce que l’on appelle
un train d’impulsion de Dirac :
h(t)
t
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Amplitude
Fréquence
f Fe 2Fe 3Fe 4Fe
Dans la pratique, d’autres petites raies apparaissent : elles sont produites par les imperfections de la numé-
risation, on les appelle des artefacts :
Amplitude
Fréquence
f Fe 2Fe 3Fe 4Fe
sa (t) Amplitude
A 8A
π2
1
Décroissance en
n2
t
T
8A
32 π 2 Fréquence
-A
f 3f 5f 7f
(a) Signal triangulaire (b) Spectre correspondant
En effet, la décomposition en série de Fourier de ce signal triangulaire (pair) s’exprime de la façon suivante :
∞
8A 1 1 8 A X cos[(2n + 1) ω t]
sa (t) = [cos(ω t) + cos(3 ω t) + cos(5 ω t) + . . . ] = ]
π2 32 52 π 2 n=0 [(2n + 1) π]2
spectre de sa (t)
f
0 fmax
Amplitude
Fréquence
−2Fe −Fe 0 fmax Fe 2Fe
Amplitude
Fréquence
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
On constate ici que les intervalles spectraux ne sont plus disjoints. Dans ce cas, il ne sera plus possible
d’analyser le signal échantillonné et de retrouver le signal sa (t).
On parle alors de repliement de spectre (ou "aliasing") : certaines raies spectrales répliquées autour de la
fréquence Fe empiètent sur l’intervalle de fréquence [0,fmax ] occupé normalement par le spectre du signal
sa (t) d’origine :
Amplitude
fmax
Fe
2 Fréquence
Fe − fmax Fe Fe + fmax
Avec l’exemple précédent, on voit que pour éviter le recouvrement des spectres, il faut que fmax soit
Fe
inférieure à .
2
La condition pour que les intervalles spectraux soient disjoints est appelée condition de Nyquist-
Shannon :
— Exemples —
• Les signaux téléphoniques ont leur spectre limité à environ 3,4 kHz et sont échantillonnées
à 8,0 kHz.
• Les signaux musicaux possèdent une fréquence maximale à 20 kHz. Leur fréquence
d’échantillonnage minimale d’échantillonnage est Fe ≃ 44 kHz.
b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b b b b b b
b b b b b
4
F IGURE 11.6 – Observation du disque avec le stroboscope pour fS = fD
3
Amplitude
f′ f
Fe 3Fe Fe
2 4
11.3.8 Conclusion
On voit que l’échantillonnage, comme toute opération non linéaire que constitue ici la multiplication du
signal d’origine par le peigne de Dirac, enrichit le spectre. En effet, des fréquences nouvelles apparaissent
dans le spectre se . Cette apparition n’est évidement pas souhaitée et il faudra les éliminer.
11.4 Quantification
— Quantifi ation —
La quantification va attribuer à chaque échantillon une valeur numérique selon un principe de pro-
portionnalité. Ainsi, tout convertisseur analogique-numérique (CAN) a besoin d’une valeur de ré-
férence (en général la masse) et d’une valeur seuil qui sera choisie en fonction de la précision
voulue.
Par exemple, si on a un convertisseur 16 bits, qu’on utilise la masse comme référence et 5 V comme valeur
5
du seuil, on aura une précision de ≃ 76 µV , ce qui signifie qu’une variation d’une unité numérique-
65536
ment parlant équivaudra à une variation analogique de 76 µV . Le nombre de valeurs possibles avec 16 bits
vaut donc 65536 = 216 .
Deux facteurs vont jouer sur la qualité du signal :
• La fréquence d’échantillonnage. En effet, plus elle sera élevée plus le signal numérisé sera fidèle au
signal analogique initial. Une fréquence trop élevée demandera cependant un espace mémoire très
important.
• Le nombre de bits sur lequel on code les valeurs qui va nous donner le nombre de valeurs que peut
prendre un échantillon et la précision en fonction de la valeur seuil.
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
5
3.75 3
2.5 2
1.25 1
0 0 t(s)
3 3 1 1 1 3 3 3 1 2 1 3 3 3 2
11 11 01 01 01 11 11 11 01 10 01 11 11 11 10
On appelle pas, la différence de tension entre deux valeurs successives.
Il dépend du nombre n de bits utilisés pour le codage ainsi que de la valeur de la tension maximale Umax :
Umax
p=
2n
Le signal rouge est la restitution en tension du signal numérique. Il n’est pas ici très ressemblant au signal
analogique . . .
Il faut donc effectuer un découpage plus fin.
Sur 3 bits, on obtient :
U (V)
5
4.375 7
3.75 6
3.125 5
2.5 4
1.875 3
1.25 2
0.625 1
0 0 t(s)
7 5 2 3 2 5 7 7 2 3 1 5 6 7 3
111 101 010 011 010 101 111 111 010 011 001 101 110 111 011
Sur 4 bits, on a :
U (V)
5
4.6875 15
4.375 14
Même si le signal rouge s’approche un peu de signal réel, la ressemblance n’est tout de même pas frappante :
il faut augmenter la fréquence d’échantillonnage et diminuer le pas.
U (V)
5
4.6875 15
4.375 14
4.0625 13
3.75 12
3.4375 11
3.125 10
2.8125 9
2.5 8
2.1875 7
1.875 6
1.5625 5
1.25 4
0.9375 3
0.625 2
0.3125 1
0 0 t(s)
U (V)
t(s)
Les ordinateurs actuels codent les informations sur 64 voire 128 bits, c’est à dire une suite de 64 ou 128
zéros ou uns par échantillon et à des fréquences avoisinant le M Hz.
Cela laisse imaginer le nombre d’informations à stocker !
11.5 Filtrage
Une fois le signal échantillonné, il faut se débarrasser des composantes inutiles. En effet les spectres créés
autour des différents multiples de la fréquence d’échantillonnage sont gênants. Pour les éliminer, on uti-
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Fe
lise en général un filtre passe-bas analogique dont la fréquence de coupure fc est inférieure à , qui va
2
supprimer les harmoniques supérieurs à cette fréquence.
11.6.1 Énoncés
2π
Un disque animé d’une vitesse de rotation angulaire ωD = est filmé par une caméra qui délivre 25
TD
images par seconde. Le sens de rotation du disque n’est pas connu a priori. Sur le disque est collé un repère
comme dans le cours. Nous noterons Ti la durée écoulée entre 2 prises de vue et α l’angle entre les 2 re-
pères. Voici ci-dessous deux images successives de la scène filmée :
b b b b
disque :
TD
b b b
α
b
1. si Ti < ,
2
b b b b
TD
b b b b 2. si < Ti < TD ,
b b
2
(a) t = 0 (b) t = Ti 3. si Ti > TD .
P 13 Électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 252
13.1 Champ électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 252
13.1.1 Champ électrostatique créé par une charge ponctuelle . . . . . . . . . . . . . . . . . . 252
13.1.2 Champ électrostatique créé par un ensemble de charges ponctuelles . . . . . . . . . . . 252
13.1.3 Champ créé par une distribution continue de charges . . . . . . . . . . . . . . . . . . . 253
[Link] Distribution volumique de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . 253
[Link] Distribution surfacique de charges . . . . . . . . . . . . . . . . . . . . . . . . . . . 253
[Link] Distribution linéique de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 253
13.2 Symétries et invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 254
[Link] Principe de Curie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 254
[Link] Éléments de symétrie du champ électrique . . . . . . . . . . . . . . . . . . . . . . . 254
13.3 Circulation du champ électrostatique, potentiel électrostatique . . . . . . . . . . . . . . . 255
13.3.1 Circulation d’un champ de vecteurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . 255
13.3.2 Circulation du champ électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . 256
13.3.3 Potentiel électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 257
13.3.4 Lignes de champ et surfaces équipotentielles . . . . . . . . . . . . . . . . . . . . . . . 257
13.3.5 Théorème de l’énergie cinétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 258
13.4 Flux du champ électrostatique : Théorème de Gauss . . . . . . . . . . . . . . . . . . . . 259
13.4.1 Flux d’un champ de vecteurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259
13.4.2 Relation de continuité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 259
13.4.3 Conditions d’application du théorème de Gauss . . . . . . . . . . . . . . . . . . . . . 260
13.5 Distributions à haut degré de symétrie . . . . . . . . . . . . . . . . . . . . . . . . . . . . 260
13.5.1 Méthode générale pour le théorème de Gauss . . . . . . . . . . . . . . . . . . . . . . . 260
13.5.2 Sphère uniformément chargée en volume . . . . . . . . . . . . . . . . . . . . . . . . . 260
[Link] Symétries et invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 261
[Link] Choix de la surface de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 261
240 TABLE DES MATIÈRES
P 14 Magnétostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 276
14.1 Théorème de superposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 276
14.2 Symétries et invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 276
14.3 Caractère conservatif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 277
14.4 Théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 279
14.4.1 Orientation d’une surface et du contour associé . . . . . . . . . . . . . . . . . . . . . . 279
14.4.2 Forme globale ou intégrale du théorème d’Ampère . . . . . . . . . . . . . . . . . . . . 279
14.4.3 Exemples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 281
14.4.4 Divergence et rotationnel du champ magnétique . . . . . . . . . . . . . . . . . . . . . 281
[Link] Divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 281
[Link] Rotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 281
14.5 Distributions à haut degré de symétrie . . . . . . . . . . . . . . . . . . . . . . . . . . . . 282
14.5.1 Méthode générale pour le théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . 282
14.5.2 Fil rectiligne infini . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 282
[Link] Symétries et invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 282
[Link] Choix du contour d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 283
[Link] Application du théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . 283
14.5.3 Câble cylindrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 284
[Link] Symétries et invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 284
[Link] Choix du contour d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 285
[Link] Application du théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . 285
14.5.4 Solénoïde infini . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 286
[Link] Symétries et invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 286
[Link] Choix du contour d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 286
[Link] Application du théorème d’Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . 287
14.6 Ordres de grandeur du champ magnétique . . . . . . . . . . . . . . . . . . . . . . . . . . 288
14.7 Exercices : magnétostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 289
14.7.1 Énoncés . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 289
CPGE TSI
Saint Joseph - LaSalle
244
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
245
P 12
En première année a été étudiée l’interaction gravitationnelle. C’est une force qui s’exerce à distance mais
ce n’est pas la seule. En effet, la force de Lorentz s’exerce également à distance. Celle-ci a été utilisée
aussi en première année pour décrire le mouvement de particules chargées dans un champ électrique ou
magnétique. En deuxième année, nous allons expliquer comment sont créés ces champs.
12.1 Généralités
Comme déjà vu en première année, une particule de charge q placée en M et animée d’une vitesse →
−
v,
→−
− →
placée dans un champ électromagnétique [ E , B ], subit une force :
−
→ →
− →
−
FL = q E (M,t) + q −
→
v ∧ B (M,t)
Une distribution D de charges et de courants crée en tout point M de l’espace un champ électromagnétique
→
− →
− →−
− →
[ E (M,t) , B (M,t) ]. Ce couple [ E , B ] constitue une entité indissociable.
→
− →
−
En régime variable, les champs électrique E et magnétique B dépendent à la fois de l’espace et du temps.
En régime stationnaire (ou permanent), ces champs sont indépendants du temps. Il apparaît alors un décou-
plage et on peut alors étudier séparément les deux champs électrique et magnétique. C’est ce que l’on fera
→
− →
−
dans les chapitres sur l’électrostatique (pour E ) et sur la magnétostatique (pour B ).
La matière est électriquement neutre. Cette neutralité résulte d’une fine compensation entre les charges
positives et les charges négatives dans la matière.
La charge électrique est une grandeur extensive, conservative (non créée et non détruite) et quantifiée.
— Charge é le trique —
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Une charge électrique est dite ponctuelle quand on peut concentrer la totalité de la charge en un
point. On la note q.
Si la charge électrique ne peut pas être concentrée en un seul point (ou en différents points), on envisage une
modélisation de même nature que pour les fluides : une modélisation continue (par opposition à discrète).
• distribution volumique,
• distribution surfacique,
• distribution linéique.
Soit une distribution volumique de charges, de charge totale QT , contenue dans un volume total τ .
Soit M un point appartenant au voisinage duquel on définit un volume élémentaire d3 τ(M) . Ce volume
élémentaire d3 τ(M) contient la charge charge élémentaire d3 q(M) .
d3 q(M)
ρ(M) =
d3 τ(M)
ou plus simplement :
dq(M)
ρ(M) =
dτ(M)
Les unité s :
ρ(M) = cte ∀M
Q T = ρ0 τ
On peut également définir la densité volumique de particules, c’est-à-dire le nombre de particules par unité
de volume (exprimée en m−3 ). La densité volumique de charges ρ vaut :
ρ = nq
Dans le cas où coexistent plusieurs porteurs de charge différents, la densité volumique de charges peut
s’exprimer par :
X
ρ= nk qk
k
Si l’épaisseur ou la section sont négligeables devant la dimension du corps considéré, on peut effectuer les
approximations suivantes :
ou plus simplement :
dq
σ=
dS
L’explication est assez simple et provient de l’intégration, par exemple suivant z, de l’expression pour
la distribution volumique :
Z Z z=h
d3 q(M) = ρ(M) d3 τ(M)
z=0
qui donne :
Z
d2 q(M) = ρ(M) h d2 S(M)
ou encore :
dq
λ=
dℓ
QT = σ0 S
QT = λ0 ℓ
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Soient q1 et q2 deux charges ponctuelles. La force d’interaction électrostatique exercée par q1 sur q2 est
donnée par la loi de Coulomb :
— Loi de Coulomb —
−−−→ q1 q2 − → −−−→
F1→2 = er = −F2→1
4 π ε0 r 2
1
Cette force est proportionnelle à . On dit qu’elle est newtonienne. C’est une force conservative.
r2
Les unité s :
1
Numériquement, ε0 = ≃ 8,842.10−12 ≃ 10−11 F.m−1
36 π 109
⇒ Activité 12. 1
Calculez le rapport entre la force électrique exercée par le noyau d’un atome d’hydrogène sur son électron
(r = 10−10 m) et la force de gravitation exercée entre ces deux mêmes éléments.
On donne :
G = 6,67.10−11 S.I. et :
N
−−→ X qi q ′ − →
Fres = ui
i
4 π ε0 ri2
avec :
−−−→′
−
→ Mi M
ui = −−−→
||Mi M ′ ||
Ceci constitue le principe de superposition. On postule qu’il y a linéarité entre la cause et les effets.
On peut généraliser ceci pour une distribution volumique, surfacique ou linéique.
Le débit de charge à travers une surface est appelé l’intensité du courant électrique :
dq
i=
dt
Les unité s :
−
→
S
dh = v dt
v , les particules
qui traversent la surface S pendant la durée dt sont contenues dans un cylindre de base S et de hauteur
dh = v dt. Le nombre de ces particules est égal à :
dN = n S dh = n S v dt
i = nqS v
Si la normale à la surface et la vitesse ne sont pas colinéaires, l’angle entre les 2 intervient :
i = n q v S cos θ
−
→
S
θ
On obtient alors :
→ −
− →
i= j ·S
Cette relation n’est bien sûr valable que si la densité de courant j est uniforme sur toute la surface.
Dans un cas plus général, on pourra exprimer le courant électrique par la relation :
ZZ
− −−
→ →
i= j · d2 S
S
L’intensité du courant électrique est égal au flux, à travers une section, du vecteur densité de courant.
Comme déjà vu, une particule de charge q placée en M et animée d’une vitesse −
→
v , placée dans un champ
→−
− →
électromagnétique [ E , B ], subit une force :
−
→ →
− →
−
FL = q E (M,t) + q −
→
v ∧ B (M,t)
ou plus simplement :
En considérant un élément de volume dτ(M) contenant des porteurs de charges de densité volumique nk ,
de charge qk et animés des vitesses −
→
vk , on peut alors exprimer la force élémentaire s’exerçant sur toutes les
particules de ce volume élémentaire :
−−→ X → X
− →
−
dFL = nk dτ(M) qk E + nk dτ(M) qk − →
vk ∧ B
k k
Soit :
−−→ X →
− X →
−
dFL = nk qk E dτ(M) + nk qk −
→
vk ∧ B dτ(M)
k k
→ X
−
En exprimant le vecteur densité volumique de courant j = nk qk −
→
vk et la densité volumique de charges
X k
ρ= nk qk , on peut ainsi exprimer la force de Lorentz volumique (ou densité volumique de la force de
k
Lorentz) :
−−→
−
→ dFL →
− −→ − →
fL = =ρE + j ∧B
dτ
P 13
Électrostatique
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
L’intérêt de définir un champ est qu’on se libère du sujet d’étude. Un champ est une propriété intrinsèque
de l’espace.
−−−−→
Soient deux charges ponctuelles, P (q) de charge q, et M (q ′ ) de charge q ′ . La force FP →M exercée par q
sur q ′ , distantes de r, est donnée par la loi de Coulomb.
−−−−→
FP →M →
− q →
−
′
= E (M) = er
q 4 π ε0 r 2
+q +q
• •
Soient une distribution de charges ponctuelles, Pi de charges qi .
Cette distribution est qualifiée de discrète. En procédant par su- M
perposition, on peut exprimer le champ électrique créé en un
−
→ •
point M de l’espace : dE ′
→
−
dE
X X qi −−→
−
→ qi −→= Pi M
E (M) = 2 e r,i
4 π ε0 ri 4 π ε0 (Pi M )3
i i −
→
E
⇒ Activité 13. 1
On peut étendre l’expression précédente à une distribution D continue. Il faut alors considérer que chaque
élément de charge élémentaire dq, situé en P , crée en un point M de l’espace le champ électrique élémen-
taire :
−−→
−→ dq PM
dE (M) =
4 π ε0 (P M )3
Soit P un point quelconque d’une distribution volumique de charge, au voisinage duquel on définit un
volume élémentaire d3 τ(P ) , contenant la charge élémentaire d3 q(P ) = ρ(P ) d3 τ(P ) .
Soit M un point quelconque de l’espace.
En partant de l’expression du champ créé par un charge ponctuelle (Méthode de calcul directe), on obtient :
ZZZ −−→
−
→ 1 ρ(P ) P M 3
E (M) = d τ(P )
4 π ε0 (τ ) (P M )3
Appliqué à l’électrostatique, ce principe de Curie implique que tous les éléments de symétrie pour la charge
sont également éléments de symétrie pour le champ. Donc, si la distribution de charges admet un plan de
→
−
symétrie, alors E est contenu dans ce plan.
−
→
E
Le champ électrique appartient donc aux plans de symétrie d’une distribution de charges.
Le champ électrique est normal aux plans d’antisymétrie d’une distribution de charges.
⇒ Activité 13. 2
Représenter sur un schéma l’illustration du principe d’antisymétrie avec 2 charges +q et −q placées de part
et d’autre d’un plan d’antisymétrie (P ′ ).
De plus, le champ électrique possède les mêmes propriétés d’invariance (rotation, translation)
que la distribution de charges qui en est la source.
— Cir ulation —
−−−→
Soit F (M ) un champ de vecteurs au point M .
−−−→
On appelle circulation de F (M ) entre deux points M1 et M2 , le scalaire défini par :
Z M2 Z M2
−−−→ − →
C= dC(M ) = F (M ) · dℓ
M1 M1
— Proprié té —
ellaSaL - hpesoJ tniaS tnemessilbatÉ - - IST EGPC - eimihC - euqisyhP
Soit P (q) une charge ponctuelle de charge q placée à l’origine d’un repère.
→
−
La circulation du champ E (M) entre deux points quelconques M1 (r1 ) et M2 (r2 ) est :
Z M2
−
→ →
−
C= E (M) · dℓ
M1
→
−
En explicitant dℓ en coordonnées cylindriques, on obtient :
Z M2 Z r2
q −
→ q dr q 1 1
C= er · (dr →
−
er + r dθ →
−
eθ + dz →
−
ez ) = = ( − )
M1 4 π ε0 r 2 r1 4 π ε0 r 2 4 π ε0 r1 r2
−
→
En prenant r1 = r2 , on voit donc que E (M) est un champ à circulation conservative.
où (Σ) est une surface orientée s’appuyant sur le contour fermé (Γ).
Ce résultat étant vraie quelle que soit la surface (Σ), cela signifie que :
−→ −→ − →
rot E = 0
→
−
C’est une relation locale pour le champ E .