Accueil/Leçons/Leçon 4

Le premier principe de la thermodynamique

Travail, chaleur, énergie interne et premier principe.

Énergie internePremier principeTravail des forces de pressionTransfert thermiqueFonction d'étatDifférentielle exacteTransformation quasi-statiqueTransformation adiabatiqueExtensivité

La leçon précédente a introduit les systèmes thermodynamiques, les états d'équilibre et les fonctions d'état. Nous allons maintenant répondre à une question centrale : comment intégrer l'énergie du système à cette représentation macroscopique ?

Dans certains cours, cette question est traitée trop rapidement : le premier principe y est présenté comme une simple conservation de l'énergie, puis vient directement la formule du bilan énergétique sans justifications supplémentaires. Nous verrons ces formules aux sections 3 et 4.3, et le lecteur pressé pourra en effet s'y rendre directement (voir aussi la synthèse rapide en Section 7).

Mais la construction qui aboutit à ces formules est en réalité subtile. D'abord, le premier principe affirme davantage que la conservation de l'énergie, qui semble évidente (de nos jours) d'un point de vue microscopique. Nous allons montrer que ce principe comble en réalité le vide logique entre l'existence d'une énergie UmicroU_{\mathrm{micro}} définie mécaniquement sur un espace comportant de l'ordre de 6N6N coordonnées, et celle d'une fonction d'état énergie UthU_{\mathrm{th}} ne dépendant que de quelques variables macroscopiques.

On montrera ensuite que l'additivité de l'énergie est une hypothèse supplémentaire nécessaire pour passer de la conservation de l'énergie d'un système isolé à un bilan d'énergie aux frontières d'un système fermé. C'est cette dernière étape qui permet d'enfin préciser le statut épistémique et logique de la chaleur : celle-ci est définie comme tout transfert d'énergie restant dans ce bilan, une fois pris en compte le travail macroscopique effectué sur un système fermé : Q  =def  ΔUWQ \equiv \Delta U - W.

Dans cette leçon, nous considérerons d'abord les systèmes fermés : aucune matière ne traverse leur frontière. La raison en est simple : dans un système ouvert, la matière transporte elle-même de l'énergie à travers la frontière et complique le bilan énergétique. Le cas des systèmes ouverts sera abordé plus loin dans le cours.

1. L'énergie interne microscopique

Considérons un système classique, non quantique et non relativiste dans un référentiel galiléen1. Considérons NN particules ponctuelles, de masses mim_i, de positions ri\vec r_i et de vitesses vi\vec v_i. À tout instant, son état microscopique est complètement déterminé par

Note 1 : Ces cas nécessitent un traitement adapté, hors-programme ici.
(r1,...,rN,p1,...,pN)Γ,avecpi=mivi,\boxed{ (\vec r_1,...,\vec r_N, \vec p_1,...,\vec p_N) \in\Gamma, \, \mathrm{avec} \,\,\, \vec p_i=m_i\vec v_i, }

Γ\Gamma désigne l'espace des phases classique du système, de dimension dimΓ=6N\dim\Gamma=6N dans ce modèle simple. On suppose que toutes les forces internes, par exemple les forces d'attraction ou de répulsion entre molécules, sont conservatives. On peut leur associer une énergie potentielle totale Epint(r1,...,rN),avecFiint=riEpintE_p^{\mathrm{int}}(\vec r_1,...,\vec r_N), \, \mathrm{avec} \,\,\, \vec F_i^{\,\mathrm{int}} =-\nabla_{\vec r_i}E_p^{\mathrm{int}}. Le système peut également être soumis à des forces conservatives extérieures, par exemple son poids. On leur associe de même une énergie potentielle extérieure Epext(r1,...,rN)E_p^{\mathrm{ext}}(\vec r_1,...,\vec r_N).

Le système n'est pas nécessairement au repos dans le référentiel choisi. On peut à tout instant définir son centre de masse et sa vitesse V\vec V, éventuellement dépendante du temps. En introduisant alors la masse totale du système et les vitesses relatives au centre de masse

M=i=1Nmi,V=1Mi=1Nmivi,ui=viV,M=\sum_{i=1}^{N} m_i, \qquad \vec V=\frac{1}{M}\sum_{i=1}^{N} m_i\vec v_i, \qquad \vec u_i=\vec v_i-\vec V,

on peut alors écrire l'énergie cinétique totale comme

Ec=12i=1Nmivi2=12i=1Nmi(V+ui)2.E_c=\frac12 \sum_{i=1}^{N}m_i \vec {v_i}^2 = \frac12 \sum_{i=1}^{N}m_i (\vec V + \vec u_i)^2.

En développant le carré, on constate que les termes croisés s'annulent, puisque imiui=0\sum_i m_i\vec u_i=\vec 0. Ce résultat est connu sous le nom de théorème de König en mécanique classique :

Ec=12MV2+12i=1Nmiui2+Vi=1NmiuiimiviMV=0E_c =\frac12MV^2 +\frac12 \sum_{i=1}^{N}m_i \vec{u_i}^2 +\vec V\mathbin{\cdot} \underbrace{\sum_{i=1}^{N} m_i\vec u_i}_{\sum_i m_i\vec v_i-M\vec V=\vec0}

Il vient :

Ec=Ecmacro+Ecmicro\boxed{E_c = E_c^{\mathrm{macro}} + E_c^{\mathrm{micro}}}

(1)

L'énergie cinétique totale se décompose donc toujours2 en l'énergie cinétique macroscopique du centre de masse, et celle des mouvements relatifs résiduels, « désordonnés », qui constituent ce que l'on appelle l'agitation thermique.

Note 2 : La translation du centre de masse n'épuise pas nécessairement l'énergie cinétique macroscopique. Le système peut posséder un mouvement ordonné de rotation globale. Dans ce cas, on peut écrire la décomposition vi=V+ω×(riR)+ui,\vec v_i = \vec V +\vec\omega\times(\vec r_i-\vec R) +\vec u_i',R\vec R est la position du centre de masse et où ui\vec u_i' désigne maintenant le mouvement résiduel après soustraction de la translation et de la rotation globales. En choisissant ω\vec\omega de telle sorte que ce mouvement résiduel possède un moment cinétique total nul par rapport au centre de masse, on obtient Ec=12MV2+12ωICMω+12i=1Nmiui2,E_c = \frac12 MV^2 + \frac12\vec\omega\cdot \mathbf I_{\mathrm{CM}}\vec\omega + \frac12\sum_{i=1}^{N}m_i {u_i'}^2,ICM\mathbf I_{\mathrm{CM}} est le tenseur d'inertie par rapport au centre de masse. La séparation entre énergie cinétique macroscopique et énergie cinétique microscopique reste donc possible.

En regroupant les différents termes, l'énergie mécanique totale s'écrit finalement

Etot=Ecmacro+Epext+[Ecmicro+Epint(r1,...,rN)+Eintautres]eˊnergie interne.E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +\underbrace{\left[ E_c^{\mathrm{micro}} +E_p^{\mathrm{int}}(\vec r_1,...,\vec r_N) +E_{\mathrm{int}}^{\mathrm{autres}} \right]}_{\text{énergie interne}}.

Remarque : le terme EintautresE_{\mathrm{int}}^{\mathrm{autres}} regroupe les contributions des degrés de liberté internes que ce modèle ponctuel minimal ne décrit pas, par exemple les degrés de liberté rotationnels ou vibrationnels des molécules. Si ceux-ci sont modélisés explicitement, l'espace des phases Γ\Gamma doit être élargi en conséquence.

Définition 1 (Énergie interne microscopique)
Nous noterons Umicro=Ecmicro+Epint+EintautresU_{\mathrm{micro}}= E_c^{\mathrm{micro}} +E_p^{\mathrm{int}} +E_{\mathrm{int}}^{\mathrm{autres}} l'énergie interne microscopique. C'est une fonction définie sur l'espace des phases, Umicro:ΓR,U_{\mathrm{micro}}:\Gamma\longrightarrow\mathbb R,

qui dépend, en général, d'un nombre considérable de variables, et qui satisfait :

Etot=Ecmacro+Epext+Umicro.E_{\mathrm{tot}} =E_c^{\mathrm{macro}}+E_p^{\mathrm{ext}}+U_{\mathrm{micro}}.

(2)

Dans la suite, nous considérerons des systèmes macroscopiquement au repos, dépourvus de rotation globale, et dont l'énergie potentielle extérieure ne varie pas. À un choix d'origine près, on pourra donc écrire

Etot=Umicro.E_{\mathrm{tot}}=U_{\mathrm{micro}}.

Nous allons maintenant introduire la fonction d'état énergie thermodynamique, provisoirement notée UthU_{\mathrm{th}}, et comprendre pourquoi elle ne doit pas être confondue a priori avec UmicroU_{\mathrm{micro}}.

2. L'énergie interne thermodynamique

Considérons comme système une casserole d'eau placée sur une plaque chauffante. Nous restons dans le cadre simple décrit plus haut : la casserole est macroscopiquement immobile et son énergie potentielle extérieure ne varie pas. On peut donc écrire Etot=UmicroE_{\mathrm{tot}}=U_{\mathrm{micro}} à une constante près.

Lorsque l'on chauffe l'eau, aucune énergie cinétique macroscopique n'apparaît et la casserole ne change pas de position. En admettant la conservation stricte de l'énergie, sous toutes ses formes, l'énergie qui lui est apportée se retrouve donc nécessairement dans son énergie interne microscopique :

ΔEtot=ΔUmicro.\Delta E_{\mathrm{tot}}=\Delta U_{\mathrm{micro}}.

La thermodynamique cherche à décrire exactement la même variation. On souhaiterait donc construire une grandeur UthU_{\mathrm{th}} telle que

ΔUth=ΔUmicro.\Delta U_{\mathrm{th}}=\Delta U_{\mathrm{micro}}.

La solution semble évidente. Ne suffit-il pas d'identifier l'énergie interne thermodynamique à l'énergie interne microscopique ?

Cette identification serait en fait mathématiquement incorrecte, puisque ces deux fonctions n'ont pas le même domaine de définition. UmicroU_{\mathrm{micro}} est une fonction définie sur l'espace des phases et dépend d'au moins 6N6N variables, tandis que UthU_{\mathrm{th}} ne dépend que d'un petit nombre dd de variables indépendantes décrivant l'état d'équilibre X=(x1,...,xd)EX=(x^1,...,x^d)\in\mathcal E (on rappelle que les xix^i sont des variables d'état macroscopiques, comme la température, la pression ou le volume).

Ces deux fonctions doivent néanmoins être reliées de telle sorte que Uth(X)U_{\mathrm{th}}(X) représente, à l'échelle macroscopique, l'énergie interne correspondant à l'état d'équilibre XX. Il faut expliquer comment fonctionne cette réduction dimensionnelle. Cette question est profonde et dépasse le cadre de ce cours d'introduction. Une description microscopique est nécessaire si l'on souhaite dériver ce passage entre les deux échelles. La sous-section 6.2, hors programme, donne les grandes lignes de la manière dont la physique statistique résout ce problème.

Mais d'un point de vue strictement thermodynamique, il n'y a pas de justification possible de cette identification microscopique, et l'existence de UthU_{\mathrm{th}} comme fonction de quelques variables macroscopiques seulement doit être postulée. Ce sera le rôle du premier principe de la thermodynamique, qui n'affirme donc pas seulement la conservation de l'énergie, mais aussi, et surtout, que cette énergie peut être représentée, à l'équilibre, par une fonction d'état définie sur l'espace macroscopique E\mathcal E.

Une fois que nous avons choisi les dd variables d'état indépendantes xix^i, par exemple (TT, VV, NN) pour un gaz parfait, la fonction U(T,V,N)U(T,V,N) définit un champ scalaire sur l'espace des états ; dans l'espace augmenté des coordonnées U,T,V,NU,T,V,N, son graphe est une hypersurface, dite surface d'équilibre, dont le rôle dans la suite de ce cours sera crucial, cf. Figure 1.

L'hypersurface énergie définie sur l'espace des états d'équilibre.
Figure 1. L'hypersurface énergie définie sur l'espace des états d'équilibre.

Remarquons pour finir un corollaire immédiat qui illustre bien la différence entre UthU_{\mathrm{th}} et UmicroU_{\mathrm{micro}} : dans la thermodynamique d'équilibre développée ici, l'énergie interne thermodynamique n'est en général pas définie en dehors des états d'équilibre. Lors d'une transformation brutale d'un gaz, par exemple, il est possible, et même courant, que l'on ne puisse plus lui attribuer une température ou une pression unique. Dans ce cas, on ne peut pas non plus lui associer directement l'énergie interne thermodynamique UthU_{\mathrm{th}}, puisqu'elle n'est définie que sur l'espace des états d'équilibre E\mathcal E. Le problème n'est pas que l'énergie n'existe plus, puisque son énergie microscopique reste parfaitement définie à chaque instant, mais que le système n'est simplement plus représenté par un point XEX\in\mathcal E.

Tout ceci ayant été dit, à partir de la section suivante, la fonction mécanique UmicroU_{\mathrm{micro}} n'interviendra plus directement : nous écrirons simplement UU pour l'énergie interne thermodynamique UthU_{\mathrm{th}}, sous-entendu définie uniquement sur l'espace des états d'équilibre.

Remarque 1 (Extensions proche-équilibre)
Signalons que pour des systèmes suffisamment proches de l'équilibre, on peut étendre la description thermodynamique en agrandissant l'espace des états. Par exemple, dans l'exemple déjà rencontré de la barre métallique située entre deux sources de température, on peut introduire un champ local stationnaire de température T(x)T(x), plutôt qu'une température unique, et définir une densité locale d'énergie interne. L'énergie totale devient alors une fonctionnelle de ces champs. Nous ne développerons pas ici cette thermodynamique dite proche-équilibre.

3. Le premier principe de la thermodynamique

La section précédente a largement préparé le terrain pour énoncer maintenant :

Premier principe de la thermodynamique pour un système fermé
Pour tout système thermodynamique fermé, il existe une fonction d'état U:ER,U:\mathcal E\longrightarrow\mathbb R,

dite énergie interne, seulement déterminée à une constante additive près, au moins deux fois différentiable, telle que, pour tout état d'équilibre XEX\in\mathcal E,

Etot=Ecmacro+Epext+U(X).E_{\mathrm{tot}} = E_c^{\mathrm{macro}} +E_p^{\mathrm{ext}} +U(X).
(3)

Le premier principe affirme également la conservation de l'énergie : pour un système isolé,

ΔEtot=0.\Delta E_{\mathrm{tot}}=0.

La régularité de la fonction UU est aussi une hypothèse non triviale. En pratique, on pourra considérer UU comme C\mathcal C^\infty. Comme UU n'est déterminée qu'à une constante près à cause des énergies potentielles impliquées dans l'égalité 3, seules ses variations portent une signification physique. Pour simplifier, nous nous limiterons dans la suite aux transformations pour lesquelles les énergies cinétique et potentielle macroscopiques ne varient pas. On aura alors ΔEtot=ΔU\Delta E_{\mathrm{tot}}=\Delta U.

Considérons maintenant un système fermé SS, qui peut échanger de l'énergie, mais pas de matière, avec son environnement. Le système global

S=S+ExtS'=S+\mathrm{Ext}

peut toujours être inclus dans un système isolé. Le premier principe impose alors

ΔES=0.\Delta E_{S'}=0.

Sous l'hypothèse supplémentaire que l'énergie est additive, c'est-à-dire que leur énergie d'interaction est négligeable, on a : ES=ES+EExtE_{S'}=E_S+E_{\mathrm{Ext}}, de sorte que

ΔES=ΔEExt.\Delta E_S=-\Delta E_{\mathrm{Ext}}.

L'hypothèse d'additivité permet ainsi de traduire la conservation de l'énergie sous la forme d'un bilan entre sous-systèmes : toute énergie gagnée par le système est perdue par son environnement, et réciproquement. Il reste maintenant à préciser comment cette énergie peut être transférée à travers la frontière du système. Nous distinguerons deux modes de transfert : le travail, noté WW, qui possède une définition indépendante issue de la mécanique, et la chaleur, notée QQ.

Le bilan énergétique prendra alors la forme

ΔES=ΔU=Q+W\boxed{ \Delta E_S = \Delta U = Q+W }
(4)

WW désigne le travail reçu par le système et QQ la chaleur reçue par celui-ci.

Remarquez que l'hypothèse d'additivité de l'énergie est non triviale. À un niveau élémentaire, on pourra simplement l'admettre pour les systèmes thermodynamiques usuels. Elle n'est cependant pas toujours valide. Les interactions à longue portée, en particulier, peuvent mettre simultanément en défaut l'additivité et l'extensivité de l'énergie, voir la Section 6.4.

Il nous reste maintenant à définir ces deux modes de transfert. Mais avant cela, il est impératif de noter une convention de signe, très standard, qu"on utilise dans toute la suite de ce cours.

À retenir (La convention de signe du banquier)
Dans tout ce livre, un transfert d'énergie est compté positivement lorsqu'il est reçu par le système et négativement lorsqu'il est fourni par celui-ci. Ainsi, Q>0, W>0Q>0,\ W>0

correspondent à de l'énergie reçue par le système, tandis que

Q<0, W<0Q<0,\ W<0

correspondent à de l'énergie cédée à l'extérieur.

4. Travail, chaleur et bilan énergétique

4.1. Le travail

L'un des avantages du travail est qu'il possède une définition indépendante de la thermodynamique, héritée de la mécanique. Lorsqu'une force extérieure Fext\vec F_{\mathrm{ext}} s'exerce sur un point de la frontière du système et que ce point se déplace de drd\vec r, le travail élémentaire fourni par cette force, et donc reçu par le système, vaut

δW=Fextdr.\delta W = \vec F_{\mathrm{ext}}\cdot d\vec r.

Pour un système étendu, il faut sommer les contributions de toutes les forces extérieures qui travaillent sur sa frontière.

Le cas qui nous servira le plus souvent en pratique est celui d'un fluide contenu dans une enceinte dont une paroi est mobile, ce qui conduit à évaluer le travail des forces de pression. Considérons un gaz enfermé dans un cylindre par un piston plan de section S\mathcal S, cf. Figure 2. L'axe xx est orienté vers l'extérieur du gaz. Lorsque le piston se déplace d'une quantité dxdx, le volume varie de

dV=Sdx.dV=\mathcal S\,dx.

Notons PextP_{\mathrm{ext}} la pression extérieure exercée sur la frontière du gaz. La force reçue par le gaz est dirigée vers l'intérieur et vaut

Fext=PextSex.\vec F_{\mathrm{ext}} = -P_{\mathrm{ext}}\mathcal S\,\vec e_x.

Le travail élémentaire reçu par le gaz est donc

δWpression=Fext(dxex)=PextSdx,\delta W_{\mathrm{pression}} = \vec F_{\mathrm{ext}}\cdot(dx\,\vec e_x) = -P_{\mathrm{ext}}\mathcal S\,dx,

soit

δWpression=PextdV\boxed{ \delta W_{\mathrm{pression}}=-P_{\mathrm{ext}} dV }
(5)

La règle du banquier permet immédiatement de vérifier le signe. Lors d'une compression, dV<0dV<0, donc δW>0\delta W>0 : le gaz reçoit du travail. Lors d'une détente, dV>0dV>0, donc δW<0\delta W<0 : le gaz fournit du travail à son environnement.

À noter qu'il s'agit bien dans cette formule de la pression extérieure au système, qui n'a pas de raisons, en général, d'être égale à la pression du gaz elle-même. C'est en particulier important de le garder en mémoire dans le cas des transformations brutales, hors équilibre, où la pression interne n'est même pas définie en général.

Travail reçu par un gaz lors du déplacement d'un piston de section S. Lors d'une détente, dV>0 tandis que la force exercée par la pression extérieure est opposée au déplacement : le travail reçu est donc négatif.
Figure 2. Travail reçu par un gaz lors du déplacement d'un piston de section S\mathcal S. Lors d'une détente, dV>0dV>0 tandis que la force exercée par la pression extérieure est opposée au déplacement : le travail reçu est donc négatif.

En intégrant, on obtient la formule suivante pour une transformation finie ABA \to B au cours de laquelle la pression extérieure est connue :

Wpression[AB]=ABPextdV\boxed{ W_{\mathrm{pression}}[A \to B] = -\int_{A}^B P_{\mathrm{ext}} dV }
(6)

Cette intégrale montre que, pour calculer le travail total, il faut connaître la valeur de PextP_{\mathrm{ext}} tout le long de la transformation. Cela a une conséquence majeure : le travail ne peut en général pas être déterminé à partir des seuls états initial et final AA et BB : son expression dépend du processus suivi.

Trois cas particuliers sont importants à signaler ici.

  1. Si PextP_{\mathrm{ext}} est constante, alors Wpression[AB]=PextΔVW_{\mathrm{pression}}[A\to B] = -P_{\mathrm{ext}}\Delta V.
  2. Le long d'une transformation isochore, dV=0dV=0, et le travail des forces de pression est nul : δWpression=0\delta W_{\mathrm{pression}}=0 et Wpression[AB]=0W_{\mathrm{pression}}[A\to B]=0.
  3. Lors d'une détente dans le vide, Pext=0P_{\mathrm{ext}}=0, donc le travail des forces de pression est également nul.

Remarque 2 (Autres formes de travail)
Le travail des forces de pression n'est qu'un exemple parmi d'autres. Pour un fil que l'on étire, si LL désigne sa longueur et FextF_{\mathrm{ext}} la force de traction reçue, δW=FextdL.\delta W=F_{\mathrm{ext}}\,dL.

Pour une interface dont on augmente l'aire AA, la tension superficielle γ\gamma conduit, dans les conditions usuelles, à un travail de la forme

δW=γdA.\delta W=\gamma\,dA.

On peut également transférer de l'énergie sans déplacement mécanique visible. Par exemple, le déplacement d'une charge dqdq à travers une différence de potentiel peut conduire à un travail électrique de la forme

δWeˊlec=Vextdq,\delta W_{\mathrm{élec}}=V_{\mathrm{ext}}\,dq,

VextV_{\mathrm{ext}} est le potentiel électrique extérieur, souvent noté, dans ce cadre, Φext\Phi_{\mathrm{ext}} pour ne pas le confondre avec le volume. Ici le signe de dqdq est aussi choisi conformément à la convention du banquier.

La sous-section 6.3 précisera la forme générale des travaux extérieurs.

4.2. La chaleur

Le travail ne peut pas être le seul mode de transfert d'énergie, puisque chauffer un gaz contenu dans un récipient rigide augmente sa température et son énergie sans qu'aucun travail extérieur ne soit réalisé. De l'énergie a donc traversé la frontière du système sans avoir été transférée sous forme de travail. Ce second mode de transfert d'énergie est appelé chaleur, ou transfert thermique.

Dans la construction adoptée ici (voir sous-section 6.1 pour une construction alternative équivalente), pour une transformation physique ABA \to B reliant deux états d'équilibre AA et BB, une fois tous les travaux reçus par le système identifiés, la chaleur reçue est définie par

Q[AB]  =def  U(B)U(A)W[AB].\boxed{ Q[A \to B] \equiv U(B)-U(A)-W[A \to B]. }
(7)

de sorte qu'on retrouve ainsi la relation annoncée plus haut,

ΔU=Q+W\boxed{ \Delta U=Q+W }

Puisque le travail dépend du chemin de la transformation ABA \to B, et puisque la variation d'énergie ΔU=U(B)U(A)\Delta U = U(B)-U(A) n'en dépend pas, c'est que nécessairement la quantité de chaleur échangée en dépend elle aussi.

Dans le langage thermodynamique, on dit que l'énergie interne UU est une fonction d'état (par postulat), tandis que le travail WW et la chaleur QQ n'en sont pas. On rappelle que cela signifie : le fait que WW et QQ ne peuvent pas être des fonctions de l'état du système à un instant donné signifie que le système ne possède pas « une quantité de travail ou de chaleur », mais plutôt que WW et QQ sont des transferts d'énergie à sa frontière pendant une transformation quelconque ABA \to B, exactement ce que l'on souhaitait obtenir, au vu des diverses expériences détaillées dans la leçon n°2.

Notons que lors d'une transformation cyclique AAA \to A, le système retrouve son état initial, donc

ΔUcycle=0.\Delta U_{\mathrm{cycle}}=0.

Le premier principe impose alors

Qcycle+Wcycle=0.\boxed{Q_{\mathrm{cycle}}+W_{\mathrm{cycle}}=0.}

Une machine cyclique ne peut donc fournir indéfiniment du travail sans recevoir une quantité égale d'énergie de son environnement. C'est, historiquement, ce que l'on appelle l'impossibilité du mouvement perpétuel de première espèce.

4.3. Forme différentielle du premier principe

La relation :

ΔU=Q+W\Delta U=Q+W

relie deux états d'équilibre AA et BB. Elle est toujours vraie, dans le sens où elle ne suppose pas que tous les états intermédiaires puissent eux-mêmes être représentés par des points de l'espace des états d'équilibre E\mathcal E. Si la transformation est brutale, le système peut quitter temporairement la surface d'équilibre définie sur E\mathcal E, alors même que U(A)U(A) et U(B)U(B) restent parfaitement définies.

Un cas spécial doit maintenant être détaillé. Supposons que la transformation soit quasi-statique. À chaque instant le système est alors, par définition, dans un état d'équilibre. Cette transformation peut alors être représentée par un chemin

γE\gamma\subset\mathcal E

constitué d'états d'équilibre infiniment voisins, qui ne quitte pas la surface d'équilibre. Puisque UU est une fonction d'état différentiable sur E\mathcal E, sa variation entre deux états infiniment voisins est une différentielle exacte, notée dUdU.

On note alors δQ\delta Q et δW\delta W les transferts élémentaires de chaleur et de travail correspondants. Le premier principe prend dans ce cas la forme

dU=δQ+δW\boxed{ dU=\delta Q+\delta W }
(8)

La différence de notation entre dUd U et δW\delta W ou δQ\delta Q n'est pas que cosmétique. Elle indique que la première est une différentielle exacte, et que les deux autres n'en sont pas, ce qui est la traduction mathématique du fait que la variation d'énergie ne dépend pas du chemin suivi alors que QQ et WW en dépendent. La section suivante rappellera au lecteur novice les mathématiques nécessaires pour bien comprendre ce point crucial.

Dans le cas particulier où le seul travail est celui des forces de pression et où l'on peut identifier PextP_{\mathrm{ext}} à la pression PP du système, cette relation devient

dU=δQPdV.dU=\delta Q-P\,dV.
À retenir (Deux formulations du premier principe)
Sa forme intégrée ΔU=Q+W\Delta U=Q+W

relie deux états d'équilibre et reste utilisable même lorsque le processus intermédiaire n'est pas quasi-statique. Sa forme différentielle

dU=δQ+δWdU=\delta Q+\delta W

n'est valable que pour une transformation quasi-statique qui peut être représentée par un chemin dans l'espace des états d'équilibre.

5. Formes différentielles exactes et inexactes

Cette section rassemble le minimum de mathématiques nécessaire pour donner un sens précis à la distinction entre dUdU d'une part, δQ\delta Q et δW\delta W d'autre part. Le lecteur qui la connaît déjà peut la passer.

5.1. La différentielle d'une fonction

Soit une fonction ff des variables x1,...,xdx^1,...,x^d, supposée dérivable. Sa différentielle est l'expression

df=i=1dfxidxi,df=\sum_{i=1}^{d}\frac{\partial f}{\partial x^i}\,dx^i,

qui mesure la variation de ff lorsque l'on passe du point XX au point voisin X+dXX+dX. Le point important est le suivant : si l'on suit un chemin γ\gamma allant d'un point AA à un point BB, on obtient en sommant toutes ces variations élémentaires

γdf=f(B)f(A).\int_\gamma df=f(B)-f(A).

En particulier, sur un chemin fermé, on a :

df=0.\oint df=0.

5.2. Les formes différentielles

Considérons maintenant une expression du même type, mais qu'il ne faut absolument pas confondre avec la différentielle d'une fonction. Soit des fonctions quelconques Ai(x1,...,xd)A_i(x^1,...,x^d). On définit alors ω\omega par :

ω=i=1dAidxi.\omega=\sum_{i=1}^{d}A_i\,dx^i.

On appelle un tel objet une forme différentielle. Si l'on se donne l'expression des fonctions AiA_i, on sait l'intégrer sans problèmes le long d'un chemin γ\gamma. Mais rien ne garantit qu'il existe une fonction ff dont ω\omega serait la différentielle, c'est-à-dire telle que Ai=f/xiA_i=\partial f/\partial x^i pour tout ii.

  • Si une telle fonction ff existe, la forme est dite exacte, et on note ω=df\omega = d f. L'intégrale γω=γdf=f(B)f(A)\int_\gamma\omega = \int_\gamma df = f(B)-f(A) ne dépend alors que des extrémités.
  • Sinon, la forme est dite inexacte, et on la note alors δω\delta\omega plutôt que ω\omega, pour rappeler qu'elle n'est la différentielle d'aucune fonction. (Les mathématiciens noteraient simplement la forme ω\omega.)

5.3. Le critère de Schwarz

Un critère simple pour déterminer si une forme différentielle est exacte est le critère de Schwarz. Plaçons-nous pour simplifier à deux variables (la généralisation est immédiate), avec

ω=A(x,y)dx+B(x,y)dy.\omega=A(x,y)\,dx+B(x,y)\,dy.

Si ω\omega était exacte, on aurait A=f/xA=\partial f/\partial x et B=f/yB=\partial f/\partial y, donc

Ay=2fyx=2fxy=Bx,\frac{\partial A}{\partial y} =\frac{\partial^2 f}{\partial y\,\partial x} =\frac{\partial^2 f}{\partial x\,\partial y} =\frac{\partial B}{\partial x},

puisque l'ordre des dérivations est indifférent pour une fonction deux fois continûment dérivable. On dispose donc d'un test commode :

AyBx    ω n’est pas exacte.\boxed{ \frac{\partial A}{\partial y}\neq\frac{\partial B}{\partial x} \;\Longrightarrow\; \omega\ \text{n'est pas exacte.} }

A des subtilités près qui n'arriveront pas en thermodynamique, la réciproque est vraie.

5.4. Exemple du travail des forces de pression

Lorsque les transferts élémentaires peuvent être exprimés en fonction des variables d'état le long du chemin quasi-statique, δQ\delta Q et δW\delta W deviennent alors des formes différentielles sur E\mathcal E, et on peut vérifier si elles sont exactes ou non.

En guise d'exemple, considérons nn moles de gaz parfait subissant une transformation quasi-statique, avec Pext=PP_{\mathrm{ext}} = P. On a alors δW=PdV\delta W=-P\,dV. Mais comme P=nRT/VP=nRT/V, le travail élémentaire s'écrit

δW=PdV=nRTVdV=0×dTnRTVdV,soitA=0,B=nRTV,\delta W=-P dV = -\frac{nRT}{V}\,dV = 0\times dT-\frac{nRT}{V}\,dV, \qquad\text{soit}\qquad A=0, \quad B=-\frac{nRT}{V},

en identifiant x=Tx=T et y=Vy=V. C'est bien une forme différentielle définie sur l'espace des états. Le critère de Schwarz donne

AV=0,BT=nRV0.\frac{\partial A}{\partial V}=0, \qquad \frac{\partial B}{\partial T}=-\frac{nR}{V}\neq0.

Les deux dérivées croisées ne sont pas égales : le travail élémentaire n'est effectivement pas une différentielle exacte. Il n'existe donc dans ce cas aucune fonction d'état W(T,V)W(T,V) dont δW\delta W serait la variation, et le travail reçu entre deux états dépend bel et bien du chemin suivi. Le même raisonnement vaut pour δQ=dUδW\delta Q=dU-\delta W : comme dUdU est exacte et δW\delta W ne l'est pas, leur différence ne peut pas être exacte non plus.

À retenir (pourquoi δ\delta et non dd)
dUdU est une différentielle exacte : sa variation ne dépend que des états initial et final, et son intégrale sur un cycle est nulle. δQ\delta Q et δW\delta W sont des formes inexactes : elles dépendent du chemin suivi, et il n'existe ni fonction QQ ni fonction WW de l'état du système telles que δW=dW\delta W = d W et δQ=dQ\delta Q = dQ.

6. Pour aller plus loin

Cette section rassemble des considérations plus avancées annoncées précédemment. Elle peut être passée en première lecture.

6.1. Deux constructions possibles du premier principe

Nous avons choisi ici l'ordre logique suivant : l'existence de la fonction d'état UU est postulée par le premier principe, le travail est défini indépendamment par la mécanique, puis la chaleur est définie par le bilan

Q=ΔUW.Q=\Delta U-W.

On en tire un fait remarquable. Si la transformation est adiabatique, alors ΔU=W\Delta U = W, ce qui impose que le travail adiabatique WadW_{\mathrm{ad}} est, lui, toujours indépendant du chemin suivi.

Cela pointe vers une autre construction possible de l'énergie interne : opérationnelle, donc indépendante de toute considération microscopique, et en l'occurrence plus proche du cheminement historique. On commence par caractériser les transformations adiabatiques sans introduire préalablement la grandeur QQ, afin d'éviter toute circularité de raisonnement.

L'idée est astucieuse : on ne définit pas le transfert lui-même, mais le dispositif. Une paroi est dite adiabatique lorsque l'état du système qu'elle enferme ne peut être modifié qu'en déplaçant les coordonnées mécaniques extérieures, c'est-à-dire le piston, un agitateur ou un courant électrique. Le test est alors direct : on maintient ces coordonnées fixes et on modifie arbitrairement l'extérieur, en plongeant par exemple l'enceinte dans un bain glacé ou en l'approchant d'une flamme. Si aucune variable d'état du système ne varie, c'est que la paroi est adiabatique.

Cette caractérisation ne fait intervenir que des états et des déplacements mécaniques, jamais un transfert d'énergie : elle précède donc toute notion de chaleur. En pratique, on s'en approche par une bonne isolation, ou en opérant vite devant le temps de relaxation thermique.

Le travail reçu, lui, se mesure de façon purement mécanique ou électrique. C'est tout le contenu empirique des expériences de Joule, menées entre 1843 et 1850 et reprises depuis avec une précision croissante : une masse mm qui tombe d'une hauteur hh fournit mghmgh, sans qu'aucun transfert thermique n'intervienne. Tous les dispositifs qu'il a imaginés et réalisés en enceinte adiabatique au sens ci-dessus, agitateur à palettes, résistance chauffante ou compression d'un gaz, l'ont conduit au même constat : un même travail fourni produit un même changement d'état, ce qui constituait pour lui « l'équivalent mécanique de la chaleur », comme nous l'avons expliqué en leçon n°2. Joule fut ainsi le premier à montrer (une indication expérimentale du fait) que le travail adiabatique ne dépend pas du chemin suivi.

C'est alors Carathéodory, en 1909 [2], puis surtout Born, en 1921, qui ont proposé d'ériger ce fait expérimental en postulat. Si on l'admet, il devient facile de construire une fonction d'état UU en posant

U(B)U(A)=Wad(AB),U(B)-U(A)=W_{\mathrm{ad}}(A\to B),

puis de définir la chaleur par différence pour les transformations générales, comme nous l'avons fait.

Cette construction est reprise dans de nombreux ouvrages, notamment ceux de Pippard [3] et de Callen [4]. Elle est équivalente à la nôtre : postuler UU, puis déduire l'indépendance du travail adiabatique est équivalent à postuler l'indépendance de ce travail, et de construire UU. On trouvera dans Rosenberg [5] le récit détaillé de cette évolution conceptuelle, de Joule à Carathéodory et Born.

6.2. Micro et macro-états

Si l'on ne choisit pas de suivre la voie opérationnelle décrite ci-dessus, alors la section 2 a mis au jour une question fondamentale : comment une énergie définie sur un espace de 6N6N coordonnées peut-elle se résumer à une fonction de quelques variables macroscopiques seulement ? La physique statistique construit ce passage entre les deux échelles. En voici les grandes lignes.

En physique statistique, on appellera un point XX de l'espace des états d'équilibre E\mathcal E un macro-état. À un même macro-état correspond en général un nombre gigantesque de micro-états γΓ\gamma\in\Gamma compatibles avec les mêmes contraintes macroscopiques.

On démontre alors que les énergies des ces micro-états se distribuent de façon extrêmement concentrée autour de sa valeur moyenne lorsque NN\to\infty, limite dite thermodynamique. Autrement dit, les fluctuations relatives de l'énergie tendent vers zéro :

σ(Umicro)Umicro0.\frac{\sigma(U_{\mathrm{micro}})} {\langle U_{\mathrm{micro}}\rangle}\longrightarrow0.

On peut alors poser Uth(X)=UmicroXU_{\mathrm{th}}(X)= \langle U_{\mathrm{micro}} \rangle_X, la moyenne portant sur tous les micro-états compatibles avec le macro-état. La concentration évoquée ci-dessus est ce qui donne son sens à cette définition. Le système occupe à chaque instant un micro-état unique, et non la moyenne : ce n'est que parce que presque tous portent la même énergie que l'on peut parler de l'énergie du macro-état, et que deux préparations identiques donnent la même mesure.

On peut, pour se représenter les choses, imaginer qu'à un macro-état correspond un très grand ensemble de micro-états microscopiquement différents mais macroscopiquement indiscernables. Cette image ressemble à celle d'une classe d'équivalence sur l'espace des phases. La construction précise sera cependant légèrement différente et fera intervenir une distribution de probabilité sur l'espace des phases.

La figure suivante résume ce procédé.

Un même macro-état X est compatible avec un nombre immense de micro-états. Leurs énergies sont, à notre échelle, extrêmement concentrées autour d'une même valeur macroscopique, ce qui permet de définir U_ th(X).
Figure 3. Un même macro-état XX est compatible avec un nombre immense de micro-états. Leurs énergies sont, à notre échelle, extrêmement concentrées autour d'une même valeur macroscopique, ce qui permet de définir Uth(X)U_{\mathrm{th}}(X).

6.3. Travaux généralisés

Les différentes formes de travail rencontrées plus haut possèdent une structure commune. Nous avons vu les formes

δW=PextdV\delta W=-P_{\mathrm{ext}}\,dV

pour le travail de pression,

δW=FextdL\delta W=F_{\mathrm{ext}}\,dL

pour l'étirement d'un fil, ou encore

δW=γextdA\delta W=\gamma_{\mathrm{ext}}\,dA

pour un travail de surface. On note qu'à chaque fois, le terme différentié est extensif alors que son préfactuer est intensif. Généralisant, on écrira tout terme de travail mécanique sous la forme :

δW=iνiextdxi,\boxed{ \delta W=\sum_i \nu_i^{\mathrm{ext}}\,dx^i, }

νiext\nu_i^{\mathrm{ext}} désigne la force généralisée extérieure (en général intensive) conjuguée à la coordonnée xix^i (en général extensive), avec le signe correspondant à notre convention.

Dans le cas où les forces généralisées extérieures s'identifient aux forces thermodynamiques correspondantes du système (exemple : Pext=PP_{\mathrm{ext}}=P), on peut alors écrire le premier principe sous la forme

dU=δQ+iνidxi.\boxed{ dU=\delta Q+\sum_i\nu_i\,dx^i.}

6.4. L'extensivité de l'énergie interne

Nous avons présenté dans la leçon précédente l'énergie interne comme une grandeur extensive : multiplier la taille d'un système homogène par un facteur λ\lambda, à variables intensives fixées, multiplie son énergie par λ\lambda,

UλU.U\longrightarrow\lambda U.

Cette propriété n'est pas postulée par le premier principe. C'est une hypothèse supplémentaire que l'on prendra souvent, mais qui ne tient pas toujours. On s'attend en effet à ce que l'existence de forces de longues portée (notamment gravitationnelles) entre les constituants du système ruine son extensivité.

Précisons cela. Considérons NN constituants distribués à densité constante dans un espace de dimension DD, et supposons que leur énergie potentielle d'interaction se comporte comme

Ep(r)1rα.E_p(r)\sim\frac{1}{r^\alpha}.

À densité ρ=N/V\rho=N/V fixée, la taille linéaire typique du système croît comme

LN1/D.L\sim N^{1/D}.

Évaluons alors l'énergie d'interaction totale en trois temps.

Compter les voisins.

Fixons un constituant et demandons-nous combien il en a d'autres à distance comprise entre rr et r+drr+dr. C'est la densité multipliée par le volume de la coquille correspondante :

dn(r)=ρSDrD1dr,dn(r)=\rho\,S_D\,r^{D-1}\,dr,

SDS_D désigne l'aire de la sphère unité en dimension DD, soit 4π4\pi en dimension trois.

Sommer sur les distances.

Chacun de ces voisins contribue pour Ep(r)rαE_p(r)\sim r^{-\alpha}. L'énergie d'interaction d'un seul constituant avec tous les autres vaut donc

uρaLrD11rαdr=ρaLrD1αdr.u\sim\rho\int_a^L r^{D-1}\,\frac{1}{r^{\alpha}}\,dr =\rho\int_a^L r^{D-1-\alpha}\,dr .

La borne inférieure aa est la distance minimale d'approche, en deçà de laquelle la loi en 1/rα1/r^\alpha cesse d'être valable et qui empêche l'intégrale de diverger en r0r\to0. La borne supérieure LL est la taille du système : il n'y a pas de voisin au-delà.

Sommer sur les constituants.

On multiplie par le nombre de constituants (en divisant par deux pour ne pas compter chaque paire deux fois) :

UintN2uNρaLrD1αdr,U_{\mathrm{int}}\sim\frac{N}{2}\,u \sim N\rho \int_a^L r^{D-1-\alpha}\,dr,

où l'on a absorbé le facteur géométrique SD/2S_D/2 dans l'ordre de grandeur. Regardons alors la valeur de l'intégrale :

aLrD1αdr=[rDαDα]aL(αD).\int_a^L r^{D-1-\alpha}\,dr =\left[\frac{r^{D-\alpha}}{D-\alpha}\right]_a^L \qquad(\alpha\neq D).

Si α>D\alpha>D, l'exposant est négatif : l'intégrale converge lorsque LL\to\infty et est dominée par sa borne inférieure, donc de l'ordre de la constante aDα/(αD)a^{D-\alpha}/(\alpha-D). Si α<D\alpha<D, elle est au contraire dominée par sa borne supérieure et vaut de l'ordre de LDα=N1α/DL^{D-\alpha}=N^{1-\alpha/D}. Enfin, dans le cas α=D\alpha=D, la primitive est un logarithme et l'intégrale vaut ln(L/a)=1DlnN\ln(L/a)=\tfrac1D\ln N. En reportant, on obtient ainsi

Uintcste×{N,α>D,NlnN,α=D,N2α/D,α<D.\begin{aligned} U_{\mathrm{int}} \sim \text{cste} \times \begin{cases} N, & \alpha>D,\\[1mm] N\ln N, & \alpha=D,\\[1mm] N^{\,2-\alpha/D}, & \alpha<D. \end{cases} \end{aligned}

Physiquement, l'interprétation est claire. Lorsque α>D\alpha>D, nous avons des forces de courte portée : la décroissance de EpE_p l'emporte sur la croissance du nombre de voisins. Chaque constituant ne sent que son voisinage immédiat, de sorte que son énergie ne dépend pas de la taille du système, et le total est proportionnel à NN, ce qui est extensif.

Lorsque αD\alpha\leq D, c'est l'inverse : les voisins lointains, bien plus nombreux, l'emportent, chaque constituant sent le système tout entier, et son énergie propre croît avec la taille de celui-ci. Le total croît alors plus vite que NN, et l'extensivité est perdue. Des interactions suffisamment courtes conduisent donc naturellement à une énergie extensive, tandis que des interactions à longue portée peuvent détruire cette propriété.

En gravitation newtonienne tridimensionnelle, par exemple,

Ep(r)1r,D=3,α=1,E_p(r)\sim-\frac1r, \qquad D=3, \qquad \alpha=1,

et l'estimation précédente donne, à densité fixée,

UgravN5/3.U_{\mathrm{grav}}\sim-N^{5/3}.

Le signe négatif, que l'argument d'ordre de grandeur ne fournit pas à lui seul, traduit le caractère attractif de la gravitation.

L'extensivité de l'énergie sera massivement utilisée dans la suite de cet ouvrage. La gravitation est bien sûr universelle, mais c'est une force très faible : si l'on étudie deux volumes de gaz mis en contact, on peut la négliger en pratique. D'un autre côté, les forces internes au gaz, de type van der Waals, décroissent très rapidement (en 1/r61/r^6), de sorte que l'extensivité est dans ce cas quasi exacte.

En revanche, pour décrire les systèmes dits auto-gravitants (étoiles, galaxies, ...), la gravitation est bien sûr l'ingrédient essentiel, et toute l'analyse thermodynamique doit être reprise depuis le début puisque l'on perd nécessairement l'extensivité de l'énergie. De ce fait, la thermodynamique de ces systèmes est une science presque à part, et montre des comportements inattendus : par exemple, une étoile qui rayonne de l'énergie se réchauffe au lieu de se refroidir. Nous y reviendrons dans la partie avancée de l'ouvrage.

7. Synthèse

Résumons ce que cette leçon a mis en place, et qui servira dans toute la suite.

  • Le premier principe postule l'existence d'une fonction d'état U:ERU:\mathcal E\longrightarrow\mathbb R, l'énergie interne, et affirme la conservation de l'énergie d'un système isolé.
  • Dans la suite de l'ouvrage, nous la supposerons aussi extensive, et suffisamment continue et différentiable dans toutes ses variables.
  • Complété par l'additivité de l'énergie, le premier principe permet de formuler la conservation comme un bilan énergétique entre un système fermé et son environnement. Le travail étant défini indépendamment par la mécanique, la chaleur est alors définie comme le transfert restant, ce qui conduit à la formule utile ΔU=Q+W\Delta U=Q+W.
  • Pour une transformation quasi-statique, ce bilan d'énergie prend la forme différentielle dU=δQ+δWdU=\delta Q+\delta W, dans laquelle seul dUdU est une différentielle exacte.

L'analyse des systèmes ouverts sera discutée ailleurs dans le cours.

8. Références

Sur l'évolution conceptuelle du premier principe, de Joule à Carathéodory et Born, voir Rosenberg [5]. Pour une présentation classique, voir Pippard [3] et Callen [4].

  1. J. P. Joule, « On the Mechanical Equivalent of Heat », Phil. Trans. R. Soc. Lond. 140, 61—82 (1850)
  2. C. Carathéodory, « Untersuchungen über die Grundlagen der Thermodynamik », Math. Ann. 67, 355—386 (1909)
  3. A. B. Pippard, Elements of Classical Thermodynamics for Advanced Students of Physics, Cambridge University Press (1957)
  4. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2e éd., Wiley (1985)
  5. R. M. Rosenberg, « From Joule to Caratheodory and Born: A Conceptual Evolution of the First Law of Thermodynamics », J. Chem. Educ. 87, 691—693 (2010)
  6. A. M. Steane, « First Law, internal energy », chap. 7 dans Thermodynamics: A Complete Undergraduate Course, Oxford University Press (2017)
  7. E. A. Gislason et N. C. Craig, « Cementing the foundations of thermodynamics: Comparison of system-based and surroundings-based definitions of work and heat », J. Chem. Thermodynamics 37, 954—966 (2005)