ECOLE POLYTECHNIQUE

ECOLES NORMALES SUPERIEURES

CONCOURS D’ADMISSION 2021

MERCREDI 14 AVRIL 2021

08h00 - 12h00
FILIERE MP - Epreuve n° 5

PHYSIQUE (XULCR)

Durée : 4 heures
L’utilisation des calculatrices n’est pas
autorisée pour cette épreuve
 Thermodynamique de l’atmosphère
Ce sujet aborde la structure de l’atmosphère sous l’angle de la thermodynamique. Il comporte
trois parties. Bien que le contexte général soit établi dans la première partie, les trois parties du
sujet peuvent être abordées indépendamment. Certaines définitions et relations de thermody-
namique sont rappelées en annexe à la fin du texte, ainsi que les valeurs de certaines grandeurs
apparaissant dans les applications numériques. Les résultats des applications numériques seront
donnés avec un seul chiffre significatif, à l’exception des trois premières questions pour lesquelles
on donnera deux chiffres significatifs.

1 Structure verticale de l’atmosphère

On considère dans cette section des modèles d’atmosphères de complexité croissante, afin
d’isoler les processus physiques qui gouvernent le profil de température vertical de l’atmosphère.
On considère tout d’abord les processus radiatifs : on admet qu’un système plan infini à
température T émet par rayonnement un flux surfacique d’énergie φ = sT 4 (en watts par
mètre carré), le flux étant dirigé vers l’extérieur. Dans les situations considérées ci-après, ce
flux correspond à du rayonnement infrarouge. Le coefficient s = 5, 67 × 10−8 W.m−2 .K−4 est
appelé constante de Stefan. On néglige la courbure de la Terre dans cette première partie, et
on considère donc la surface de la Terre comme localement plate.

1) Dans un premier modèle, on considère que la surface terrestre absorbe un flux solaire incident
I = 220 W.m−2 , considéré uniforme et indépendant du temps. A l’équilibre, la Terre ré-
émet l’intégralité de l’énergie reçue par rayonnement vers le haut. En déduire la température
d’équilibre de surface T0 de la Terre dans ce modèle. Faites une application numérique en vous
aidant de l’annexe. Commentez le résultat.

2) Un raffinement du modèle consiste à inclure une atmosphère infiniment fine située à une
hauteur h au-dessus de la surface de la Terre. Cette atmosphère est transparente au rayon-
nement solaire incident, mais elle est entièrement absorbante pour le rayonnement infrarouge
émis par la Terre. Elle atteint une température d’équilibre Ta , si bien qu’elle rayonne elle aussi
un flux d’énergie sTa4 vers le haut, mais aussi un flux sTa4 vers le bas. Calculez la température
Ta d’équilibre de l’atmosphère, ainsi que la nouvelle température d’équilibre T0 de surface de
la Terre. Faites une application numérique en vous aidant de l’annexe. Commentez le résultat
obtenu. Comment s’appelle l’effet mis en évidence par ce modèle simple ?
On oublie un instant les modèles radiatifs précédents, pour se concentrer sur la thermody-
namique du fluide atmosphérique. On suppose que l’atmosphère est un gaz parfait diatomique.

3) Rappelez quels sont les deux constituants principaux de l’atmosphère ainsi que leurs pro-
portions. En déduire une valeur approchée de la masse molaire moyenne M de l’air.

1
 On donne la différentielle de l’entropie du gaz parfait:
h  1
i
dS = nCp d ln T P γ −1 , (1)

où S, T , et P désignent l’entropie, la température et la pression du gaz. On désigne de plus

par n la quantité de matière (nombre de moles), et on définit γ = Cp /Cv , où Cp est la capacité
calorifique molaire à pression constante du gaz, et Cv est sa capacité calorifique molaire à
volume constant. Ces deux grandeurs sont également reliées par la relation de Mayer, rappelée
en annexe.

4) On s’intéresse au modèle isentropique, pour lequel l’entropie du gaz est supposée indépendante
de la coordonnée verticale z. Ecrire tout d’abord l’équilibre hydrostatique d’une parcelle
d’atmosphère, puis calculer le profil de température T (z) associé à ce modèle, en fonction
de la température T0 à la surface z = 0 de la Terre, de l’accélération g de la pesanteur (sup-
posée constante), de la capacité calorifique molaire Cp de l’air (supposée constante également),
et de la masse molaire moyenne M de l’air.

5) Faites une application numérique pour le gradient de température Γ = dT /dz dans ce modèle.
On pourra considérer l’air comme un gaz parfait diatomique.
On souhaite justifier la pertinence de ce modèle d’atmosphère isentropique. Dans ce but,
on considère une atmosphère ayant un profil de température initial T (z) quelconque (les profils
de densité ρ(z) et de pression P (z) se déduisent de ce profil de température par l’équilibre
hydrostatique et la loi des gaz parfaits). On isole une parcelle de fluide à l’altitude initiale zi qui
constitue le système fermé étudié. Ce système est initialement à la température Ti = T (zi ), et à
la pression Pi = P (zi ). On suppose qu’une perturbation externe déplace la parcelle considérée
à l’altitude zi + dz, de façon adiabatique et quasi-statique (réversible), où dz représente un
déplacement infinitésimal. Le système s’équilibre donc rapidement à la pression Pf = P (zi +dz).

6) Calculez la température Tf du système à l’issue de cette transformation, en fonction de Pi ,

Pf , Ti , R et Cp .

7) Discutez le mouvement vertical ultérieur de la parcelle d’atmosphère, selon que Tf < T (zi +
dz) ou Tf > T (zi +dz). Quelle situation correspond à une atmosphère stable ? à une atmosphère
instable ? On appelle ce phénomène l’instabilité convective.

8) Le cas où Tf = T (zi + dz) s’appelle le cas “marginalement stable”. Quel est le gradient de
température vertical Γ = dT /dz de l’atmosphère dans cette situation marginalement stable ?
Commentez, en lien avec les questions 4) et 5).
On revient à un modèle purement radiatif, que l’on complexifie en prenant en compte
l’absorption et la ré-émission de rayonnement par chaque tranche élémentaire d’atmosphère.
Pour un profil donné du coefficient d’absorption de rayonnement infrarouge par l’atmosphère,
ce modèle conduit au “profil radiatif” de température représenté en figure 1.

9) Après avoir estimé le gradient de température dans les hautes et basses couches de ce modèle

2
 .*+/#-

!"#$%&'()&"*+,-

Figure 1: Profil vertical de température issu d’un modèle purement radiatif prenant en compte
l’absorption et la ré-émission de chaque couche infinitésimale d’atmosphère.

d’atmosphère, indiquez approximativement quelle partie de l’atmosphère est stable, et laquelle

est instable, suivant ce modèle.
Le phénomène d’instabilité convective fait que le modèle purement radiatif de la figure 1
n’est pas satisfaisant. A terme, l’instabilité convective conduit à des mouvements erratiques
de l’air qui affectent fortement le profil de température mesuré dans l’atmosphère réelle. On
prend en compte ce phénomène en admettant qu’à toute altitude z où le profil de température
est instable, ces mouvements erratiques se produisent jusqu’à ce que le profil de température
retrouve localement un gradient de température proche (mais légèrement plus faible, comme
nous l’expliquerons dans les questions suivantes) du gradient marginalement stable déterminé à
la question 8). Cette région de l’atmosphère où l’instabilité convective joue un rôle prépondérant
s’appelle la troposphère. La couche d’atmosphère qui n’est pas affectée par cette instabilité
s’appelle la stratosphère.

10) Commentez alors l’allure qualitative du profil de température du modèle simplifié d’atmosphère
standard utilisé en aéronautique et fourni en figure 2, en lien avec le profil du modèle purement
radiatif de la figure 1. Identifiez en particulier la troposphère et la stratosphère.
Le gradient de température mesuré dans la troposphère est moins raide aux latitudes
équatoriales. La présence d’humidité dans l’atmosphère permet d’expliquer ce phénomène. On
considère alors un système fermé constitué d’une parcelle d’atmosphère contenant un mélange
air-vapeur d’eau. On note nv la quantité de matière (nombre de moles) de vapeur et na le
nombre de moles d’air dans le système. Le rapport xv = nv /(na + nv ) désigne la fraction
molaire de vapeur d’eau dans la phase gazeuse. Cette fraction est supposée très faible devant
1, xv  1. On admet que lorsqu’il y a condensation, l’eau liquide apparaı̂t sous la forme de
fines gouttelettes qui restent en suspension dans la parcelle d’atmosphère, et qui ne tombent
pas sous l’effet de la gravité.

3
 20

15

z [km]
10

0
220 240 260 280 300
Température [K]

Figure 2: Profil vertical de température du modèle international d’atmosphère standard utilisé

en aéronautique.

11) On considère cette parcelle d’atmosphère à une altitude donnée dans le profil de température
marginalement stable calculé à la question 8). On suppose de plus que cette parcelle est initiale-
ment à la pression de vapeur saturante en vapeur d’eau. On imagine que cette parcelle subit un
léger déplacement vers le haut, de façon quasi-statique et adiabatique (réversible). Expliquez
qualitativement pourquoi, sous certaines conditions, la parcelle peut néanmoins continuer son
mouvement ascendant alors même que l’atmosphère environnante n’est pas instable vis-à-vis
du critère déterminé à la question 7). Que voit-on alors apparaitre dans le ciel ?

12) Proposez une expression approchée pour dxv en fonction de dnv et na uniquement, lorsque
la quantité de vapeur en phase gazeuse varie. On utilisera cette expression approchée dans ce
qui suit.

13) On étudie l’évolution de la parcelle au cours de son mouvement ascendant, toujours supposé
quasi-statique et sans transfert thermique avec l’atmosphère environnante, dont le profil de
température est maintenant considéré arbitraire. On considère toujours nv  na , si bien que
l’on peut continuer à décrire la phase gazeuse par la relation (1), dans laquelle on remplacera
n par na (on néglige donc la vapeur dans l’évaluation de la densité, de la quantité de matière
et de la masse molaire de la parcelle). Dans cette relation, la variation d’entropie est donnée
par dS = dHcond /T , où dHcond > 0 est l’enthalpie libérée par la condensation. En combinant
cette relation avec l’équilibre hydrostatique, montrez alors que :
dxv dT
− ∆vap H = Cp + Mg , (2)
dz dz
où T et xv sont respectivement la température et la fraction molaire de vapeur dans la phase
gazeuse, lorsque la parcelle passe à l’altitude z. ∆vap H est l’enthalpie molaire de vaporisation
de l’eau. Elle est considérée constante dans ce qui suit, de même que Cp et g.
Au voisinage immédiat de l’océan équatorial, l’air se charge en vapeur d’eau jusqu’à satu-

4
 z=h

z=0 T0 T1 T0
surface océanique
latitude
30°S 0° 30°N

Figure 3: La circulation de Hadley correspond aux mouvements atmosphériques dominants

dans un plan méridional (latitude - altitude). La température de surface au niveau de l’océan
est supérieure au voisinage de l’équateur (T1 ) qu’aux latitudes subtropicales (T0 < T1 ). Les
flèches illustrent la direction dominante du vent.

(i)
ration. La fraction molaire en vapeur d’eau initiale, xv , est alors de 4%, et la température
de l’ordre de 30 o C. Cette atmosphère humide consiste en de nombreuses parcelles isolées
d’atmosphère, identiques au système précédent, qui ont un mouvement ascendant depuis le
niveau de la mer z = 0 jusqu’au sommet de la troposphère, à une altitude h = 10 km. On
suppose que toute la vapeur est condensée lorsqu’une telle parcelle atteint l’altitude h.

14) Calculez alors le gradient de température moyen Γh = (T (h) − T (0))/h en fonction de M ,

(i)
g, Cp , h, ∆vap H, et xv . Γh est appelé le “gradient de température de l’atmosphère humide”.

15) Donnez la valeur numérique du gradient de température de l’atmosphère humide Γh , en

considérant l’atmosphère comme un gaz parfait diatomique. On donne l’enthalpie de vapori-
sation de l’eau : ∆vap H = 4, 1 × 104 J.mol−1 . Commentez le résultat, en lien avec le profil du
modèle international d’atmosphère standard fourni en figure 2.

16) On représente sur la figure 3 la circulation de Hadley, qui est le principal courant méridional
de l’atmosphère au voisinage de l’équateur. Pourquoi la température de surface est-elle supérieure
au voisinage de l’équateur (T1 ), par rapport aux latitudes subtropicales (T0 < T1 )? Com-
mentez alors la structure de cette circulation : direction des mouvements d’air, présence de
précipitations équatoriales, etc.
On souhaite établir la dépendance de la pression de vapeur saturante de vapeur d’eau
Psat (T ) avec la température. Pour cela, on considère un système fermé, formé uniquement
d’eau liquide et de vapeur d’eau en co-existence (sans air). On note Gl,m (T, P ) l’enthalpie
libre molaire de la phase liquide et Gv,m (T, P ) celle de la phase vapeur. De même, on note
Vl,m (T, P ) et Vv,m (T, P ) les volumes molaires de ces deux phases, et Sl,m (T, P ) et Sv,m (T, P )
leurs entropies molaires.

17) Exprimez la différentielle de Gl,m (T, P ) en fonction de Vl,m (T, P ), Sl,m (T, P ), et des diffé-

5
rentielles de T et P . Exprimez de même la différentielle de Gv,m (T, P ) en fonction de Vv,m (T, P ),
Sv,m (T, P ), et des différentielles de T et P . En utilisant l’extensivité de l’enthalpie libre G du
système complet, donnez son expression en fonction de Gl,m (T, P ), Gv,m (T, P ), de la quantité
de matière en phase vapeur nv et de la quantité de matière totale n.

18) On considère la situation où T et P sont imposées par un thermostat et un réservoir de

pression. On admet que G est un potentiel thermodynamique dans cette situation. Exprimez
alors la condition d’équilibre liquide-vapeur en fonction de Gv,m et Gl,m .

19) On considère deux points sur la courbe de coexistence liquide-vapeur, de coordonnées

(T, Psat (T )) et (T + dT, Psat (T + dT )). Ecrire la condition d’équilibre liquide-vapeur aux deux
points considérés en fonctions des enthalpies libres molaires des deux phases. Soustraire ces
deux conditions d’équilibre pour obtenir la relation de Clapeyron :
dPsat ∆vap H
= (3)
dT T (Vv,m − Vl,m )

En déduire le sens de variation de la pression de vapeur saturante avec la température.

20) On néglige le volume molaire de la phase liquide devant celui de la phase gazeuse, et on
assimile la vapeur à un gaz parfait à la pression Psat . On considère de plus que ∆vap H est
indépendante de la température. Montrez alors que la pression de vapeur saturante s’exprime
sous la forme :    
Psat (T ) 1 1
ln =C − , (4)
Psat (T0 ) T0 T
où T0 est une température de référence arbitraire, et où vous donnerez l’expression de C.

21) On considère que l’eau existe dans l’océan et l’atmosphère sous deux formes, H2 O16 et H2 O18 ,
formées respectivement avec les isotopes 16 et 18 de l’oxygène. Ces deux isotopes sont stables
du point de vue de la radioactivité. Les deux types de molécules d’eau se caractérisent par
des enthalpies molaires de vaporisation différentes, et donc des pressions de vapeur saturante
différentes. Expliquez comment la fraction de H2 O18 dans les calottes glaciaires fournit un
historique de la température de l’atmosphère terrestre.

2 La troposphère vue comme une machine thermique

On se propose de modéliser le courant atmosphérique représenté sur la figure 3 par une

machine thermique cyclique. La thermodynamique nous permet alors de déterminer une
borne supérieure sur la fraction du transfert thermique incident qui est transformée en travail
mécanique sous la forme de vents. Dans ce but, on revient dans un premier temps sur le fonc-
tionnement et le rendement des moteurs dithermes cycliques. Pour les applications numériques
demandées dans les deux dernières questions de cette partie, on considèrera la Terre comme
une sphère dont le rayon vaut approximativement 6400 km.

6
22) On considère tout d’abord un moteur ditherme cyclique en interaction avec un thermostat
chaud de température Tc et un thermostat froid de température Tf . Définissez le rendement
ξ de ce moteur. Etablissez la valeur théorique maximale ξc de ce rendement (le rendement de
Carnot), atteinte si les transformations associées au cycle sont réversibles.
Une machine thermique entièrement réversible n’est pas un bon modèle d’une machine
réelle : les processus mis en jeu doivent être infiniment lents, en particulier les transferts
thermiques. Si le travail fourni par la machine est non nul, ce travail est fourni en un temps
infini ! Autrement dit, la puissance engendrée (travail fourni par unité de temps) tend vers
zéro dans la limite de transformations parfaitement réversibles. Pour résoudre ce problème on
considère la situation représentée sur la figure 4. La machine comprend un fluide, le système
F à l’intérieur du cadre pointillé, qui subit de façon cyclique 4 transformations :

1. Une étape isotherme, à température T1 , pendant laquelle le fluide reçoit un transfert

thermique algébrique Q1 depuis le thermostat chaud, à température Tc . Le thermostat
est connecté au système par une résistance thermique R. Cette étape est supposée durer
un temps τ . On notera que la température au voisinage immédiat du système F (au bord
du domaine pointillé) est bien T1 : c’est la température effective de la source chaude pour
le système F.

2. Une détente adiabatique rapide du fluide. Si nécessaire, une telle détente peut être ef-
fectuée rapidement tout en restant réversible (il suffit que le temps soit long devant le
temps de propagation des ondes sonores dans le fluide). Sa durée est donc négligeable
par rapport au temps τ .

3. Une étape isotherme, à température T2 , pendant laquelle le fluide reçoit un transfert

thermique algébrique Q2 , à travers une résistance thermique R, depuis un thermostat
froid à température Tf . Cette étape dure un temps τ égal à celui de l’étape 1. On notera
encore une fois que la température au voisinage immédiat du système F est T2 : c’est la
température effective de la source froide pour le système F.

4. Une compression adiabatique rapide du fluide (qui peut néanmoins être réversible, si
nécessaire), dont la durée est négligeable par rapport au temps τ .

On note W le travail algébrique reçu par le système pendant un cycle. Les températures des
thermostats Tc et Tf sont fixées, mais on suppose que l’on peut ajuster les températures de
fonctionnement T1 et T2 du moteur pour obtenir un fonctionnement optimal.

23) Donnez les signes de Q1 , Q2 et W dans ce qui vous semble être le régime de fonctionnement
normal du moteur. Ordonnez les températures Tc , Tf , T1 et T2 dans ce régime de fonction-
nement, en justifiant très brièvement. On supposera cet ordre des températures respecté dans
la suite.

24) Exprimez les transferts thermiques algébriques Q1 et Q2 reçus par le système F pendant
les étapes 1 et 3, en fonction de R, T1 , T2 , Tc , Tf et τ .

7
 travail -W
Système F
2
résistance
résistance

thermique
thermique

source
Tc T1 T2 Tf source

chaude
R 1 3 R froide

Q1 Q2

Figure 4: Une machine cyclique est connectée aux sources chaude et froide par le biais de
résistances thermiques. Le système F subit alors de façon cyclique les 4 étapes décrites dans
le texte, représentées schématiquement par les chiffres entourés.

25) Définissez la puissance motrice P fournie par le moteur. Exprimez P en fonction de R, T1 ,

T2 , Tf et Tc .

26) A l’aide du second principe de la thermodynamique, établir une inégalité entre T1 , T2 , Tf

et Tc . Pour simplifier les calculs, on admet que le rendement maximal est obtenu lorsque les
transformations sont réversibles pour le système F. On remplacera alors cette inégalité par une
égalité dans ce qui suit.

27) A l’aide des réponses aux deux questions précédentes, éliminez T2 pour exprimer la puissance
P sous la forme :
Tc
P = × f (X, η) , (5)
2R
où l’on a introduit les paramètres X = T1 /Tc et η = Tf /Tc .

28) A l’aide des réponses aux questions précédentes (en particulier la question 23), montrez que
X > 1/2.

29) En pratique, les moteurs réels fonctionnent dans un régime où T1 (et donc T2 ) est choisie
pour maximiser la puissance motrice délivrée, à Tc et Tf fixées. Calculez la valeur X+ de X
associée à ce maximum, en fonction de η. Donnez la valeur du maximum de puissance en
fonction de R, Tc et η.

30) Donnez l’expression du rendement ξ du moteur en fonction de X et η. Montrez que le

rendement du moteur à puissance maximale vaut ξ(X+ ) = 1 − η 1/2 . Comment cette expression
du rendement se compare-t-elle au rendement de Carnot ?

31) La différence de valeur entre ces deux rendements signifie que certains processus apparais-
sant sur la figure 4 sont irréversibles. Quels sont-ils ? On définit le taux de création d’entropie
Ṡ comme le rapport de l’entropie créée par ces processus irréversibles pendant un cycle divisée

8
par la durée d’un cycle. En appliquant le second principe de la thermodynamique à un système
judicieusement choisi, calculez Ṡ en fonction de R et η dans le régime de fonctionnement à
puissance maximale du moteur.
On souhaite appliquer le modèle précédent aux courants atmosphériques. Le système est la
totalité de la troposphère, que l’on considère comme une couche de fluide de hauteur h = 10 km
à la surface de la Terre, et de densité ρ = 1 kg.m−3 considérée constante pour ces calculs d’ordre
de grandeur. Le transfert thermique Q1 reçu par la troposphère est dû au rayonnement solaire
incident qui atteint la surface terrestre. Ce rayonnement a une intensité I = 220 W.m−2 en
moyenne sur une journée et sur toute la surface de la Terre. Q2 correspond alors au transfert
thermique sous forme de rayonnement infrarouge émis au niveau du sommet de la troposphère
(on néglige la stratosphère), tandis que |W | est l’énergie qui est transformée en travail moteur
sous la forme de vents, ou courants atmosphériques. La température moyenne au niveau du
sol est Tf = 290 K, tandis que la température moyenne au sommet de la troposphère est
Tc = 220 K.

32) En supposant que le système fonctionne à puissance maximale, calculez approximativement

l’énergie cédée aux vents et courants atmosphériques par unité de temps.

33) On admet que les courants atmosphériques ainsi engendrés dissipent leur énergie cinétique
à un taux D (énergie dissipée par unité de temps), de l’ordre de D ' 0, 05 × MU 3 /h, où M
est la masse totale de la troposphère et U la vitesse typique des courants atmosphériques. En
déduire un ordre de grandeur de U . Commentez.

3 Nucléation et formation des nuages

On s’intéresse à la formation des nuages, et on considère donc un système constitué d’eau

liquide et d’un mélange gazeux air-vapeur d’eau. On note Pv la pression partielle de vapeur
dans le gaz, et Psat (T ) la pression de vapeur saturante.

34) Donnez la relation entre le potentiel chimique µl (T ) de l’eau liquide, supposé ne dépendre
que de la température, et le potentiel µv [T, Pv = Psat (T )] de la vapeur d’eau à la pression de
vapeur saturante ?
On considère maintenant une pression de vapeur quelconque, et on s’intéresse à l’évolution d’une
gouttelette d’eau liquide dans l’atmosphère. L’existence d’une interface de surface A entre le
liquide et le gaz conduit à un terme d’énergie de surface supplémentaire, σA, qu’il convient
d’ajouter à l’énergie du système (ou au potentiel thermodynamique considéré). Le coefficient
σ est une constante positive appelée tension de surface. Son unité est le kg.s−2 . Le système
considéré contient un nombre fixe de molécules d’air, et un nombre total fixe ntot = nl + nv de
moles d’eau, où nl est le nombre de moles d’eau liquide et nv le nombre de moles de vapeur
d’eau. Le système est maintenu à pression et température constante.

35) Exprimez la différentielle de l’enthalpie libre du système en fonction de σ, A, µl , µv , nl

9
 atmosphère

(air+vapeur d’eau)
particule
solide

r
θ

eau liquide

Figure 5: Une particule fine solide comporte une cavité conique ouverte sur l’extérieur. On
considère la possibilité que de l’eau liquide apparaisse sous la forme d’un domaine conique au
fond de la cavité.

et/ou leurs différentielles.

36) On suppose que la phase liquide consiste à tout instant en une gouttelette d’eau sphérique
de rayon r et de masse volumique constante ρ0 = 1000 kg.m−3 . Exprimez la différentielle de
l’enthalpie libre en fonction de µl , µv , ρ0 , r, σ, MH2 O et dr, où MH2 O est la masse molaire de
l’eau.

37) Intégrez cette relation pour obtenir l’enthalpie libre G(r) (à une constante près, qui ne joue
pas de rôle dans ce qui suit). Tracez l’allure de G(r) pour µl et µv fixés. On distinguera les cas
µl < µv et µl > µv .

38) Discutez l’évolution de la gouttelette suivant son rayon initial, en distinguant à nouveau
les cas µl < µv et µl > µv . Quand cela est pertinent, on identifiera une valeur remarquable
rc de r, dont on donnera l’expression en fonction des mêmes paramètres qu’à la question 36).
On définira également une “énergie d’activation” pour l’apparition de la phase liquide, lorsque
l’on part uniquement d’un mélange gazeux air-vapeur d’eau, dont on donnera l’expression en
fonction de ces mêmes paramètres.

39) On modélise la phase gazeuse comme un mélange de gaz parfaits. A l’aide de l’expression du
potentiel chimique donnée en annexe, que l’on appliquera à la vapeur d’eau, exprimez le rayon
critique au-delà duquel une gouttelette croı̂t spontanément en fonction de la “sursaturation”
χ = Pv /Psat (T ) > 1 et de MH2 O , σ, ρ0 , R et T .

40) A partir de vos connaissances sur la loi d’Arrhénius, ou sur le facteur de Boltzmann,
proposez une expression pour la dépendance de la vitesse d’apparition des gouttelettes dans
l’atmosphère en fonction de la sursaturation χ (cette expression sera obtenue à une constante
multiplicative près). Commentez cette dépendance en χ.
On considère la situation représentée sur la figure 5 : une particule solide en suspension dans
l’atmosphère présente une cavité conique ouverte sur l’extérieur, de demi-angle d’ouverture θ.

10
On étudie la possibilité que de l’eau liquide se forme à partir du fond de la cavité, selon la
géométrie représentée sur la figure 5. L’eau liquide occupe alors un domaine conique dont la
base est un disque de rayon r perpendiculaire à l’axe du cône. Le volume d’un tel cône est
πr2 `/3, où ` est la hauteur du cône.

41) On fait l’hypothèse que seule l’interface liquide-gaz intervient dans le terme de tension de
surface, et donc que l’interface solide-liquide n’engendre aucun terme de tension de surface.
Exprimez la différentielle de l’enthalpie libre en fonction de µl , µv , ρ0 , r, σ, MH2 O , θ et dr.

42) Comme à la question 38), identifiez une valeur remarquable de r. Exprimez l’énergie
d’activation associée en fonction de σ, MH2 O , ρ0 , µl , µv et θ.

43) Expliquez alors pourquoi la présence de particules en suspension dans l’atmosphère favorise
l’apparition de nuages.

Quelques rappels

Pour simplifier certaines applications numériques, on donne les racines quatrièmes suivantes :
(220/5, 67)1/4 ' 2, 5 et 21/4 ' 1, 2.

Thermodynamique générale

Pour un système formé de N constituants (ou phases) de potentiels chimiques µi , avec les
quantités de matière ni , la différentielle de l’énergie interne s’écrit :

dU = T dS − P dV + ΣN
i=1 µi dni , (6)

et l’enthalpie libre est définie par G = U + P V − T S. Le dernier terme de la différentielle

(6) n’apparaı̂t pas lorsque l’on considère une phase pure ou un système sans transformation de
phases (dni = 0).

Thermodynamique des gaz parfaits

On rappelle la relation de Mayer pour un gaz parfait, Cp − Cv = R, où la constante des gaz
parfaits vaut R = 8.314 J.K−1 .mol−1 . Dans un mélange de gaz parfaits, le potentiel chimique
µi d’un constituant i est donné par la relation :

(0) Pi
µi (Pi , T ) = µi (T ) + RT ln , (7)
P0
(0)
où µi est une fonction qui ne dépend que de la température, Pi désigne la pression partielle
du constituant i, et P0 = 1 bar est une pression de référence.

11
Masses molaires

On rappelle les masses molaires de l’hydrogène, du carbone, de l’azote et de l’oxygène : MH =

1, 0 g.mol−1 , MC = 12 g.mol−1 , MN = 14 g.mol−1 , MO = 16 g.mol−1 .

12