Dans une première partie, on étudie la croissance d'une couche de glace lorsque sa température de surface est controlée, en utilisant pour cela l'approximation quasi stationnaire. Puis l'on examine les effets d'un manteau de neige sur la formation de la couche de glace à la surface d'un lac. On évoque enfin dans la dernière partie l'évolution saisonnière de la glace arctique.

Toutes les parties sont largement indépendantes les unes des autres. Sauf indication contraire, la pression

est constante et égale à la pression atmosphérique moyenne, soit

Une grande attention devra être apportée aux applications numériques.

Données numériques :

Capacité thermique massique de la glace
Capacité thermique massique de l'eau
Conductibilité thermique de la glace
Conductibilité thermique de la neige
Masse volumique de l'eau
Masse volumique de la glace
Masse volumique de la neige
Enthalpie de fusion de la glace

les données précédentes sont supposées indépendantes de la température.

é

Bilan radiatif à la surface de la banquise arctique :

Coefficients

I. Le problème de Stefan

La figure 1 illustre le problème de la formation d'une couche de glace tel qu'il fut formulé dans le travail pionnier de Stefan (1891). La surface d'un volume d'eau initialement à la température de fusion

est mis en contact à l'instant

avec une paroi plane, maintenue en position fixe et à température

. Une couche de glace apparaît et se développe progressivement au sein du fluide. On note

la position de l'interface entre l'eau et la glace; la glace occupe l'espace

. Soit

le champ de température dans la glace, supposé unidimensionnel. On suppose que

Figure 1

La diffusion thermique.
a) Exprimer la loi de Fourier reliant au sein de la glace la densité de courant d'énergie au gradient de température.
b) Effectuer un bilan énergétique sur un volume élémentaire de glace pour obtenir l'équation de la diffusion thermique, dite «de la chaleur ».
c) Quelles sont les conditions aux limites pour le champ de température de la glace? Permettent-elles de déterminer ?
d) Que peut-on dire de la température au sein de l'eau?
e) Pourquoi l'eau est-elle mise en mouvement par l'avancée de l'interface?
Soient l'enthalpie massique de la glace et celle de l'eau que l'on suppose indépendantes de la température. On désigne par la vitesse de l'interface et par la vitesse verticale de l'eau.
a) En raisonnant sur un cylindre vertical de section , exprimer à l'aide de la masse d'eau qui s'est transformée en glace entre les instants et .
b) Effectuer le bilan enthalpique de cette masse entre ces deux instants (on négligera la variation d'énergie cinétique de l'eau qui gèle).
c) En déduire la relation suivante :

On suppose que est suffisamment faible pour admettre que la distribution de température dans la glace est à tout instant celle de l'état stationnaire pour l'épaisseur de glace formée à cet instant (approximation quasi stationnaire).
a) Pourquoi n'a-t-on jamais rigoureusement de régime permanent?
b) Que devient l'équation de la chaleur dans l'approximation quasi stationnaire? En déduire le profil puis le gradient de température au sein de la glace.
c) Déduire alors de l'équation (1) une équation différentielle portant sur . Montrer que où est une constante que l'on explicitera.
d) Application numérique : calculer pour . Calculer l'épaisseur de glace après un jour, une semaine, un mois, six mois.

II. Effet d'une couche de neige

On souhaite étudier l'effet d'un couvert de neige sur la croissance de la glace. On suppose qu'il existe une couche de neige d'épaisseur

constante, présente dès l'instant initial sur une très mince couche de glace (figure 2). On note

la température à l'interface neige/glace.

Figure 2

Quelle est la forme des profils de température au sein de la neige et de la glace en régime quasi stationnaire? Quelle condition doit être vérifiée à l'interface neige/glace?
Soit la composante verticale de la densité de courant d'énergie . Exprimer en fonction de , puis de . Exprimer alors en fonction de et de .
En déduire la nouvelle équation différentielle portant sur . Montrer que la solution satis-
faisant aux conditions initiales est :

où

est une longueur caractéristique que l'on explicitera.
4. Application numérique : Calculer l'épaisseur de glace obtenue après un jour, une semaine, un mois et six mois pour

.
5. La neige joue-t-elle un rôle dans la croissance des couverts de glace?

III. Variation saisonnière de la glace arctique

Dans cette partie, toutes les températures sont exprimées en degré Celsius et toutes les durées en jour. On admet que l'évolution de la banquise au-delà du cercle polaire est entièrement contrôlée par l'équilibre qu'elle entretient avec l'atmosphère. Un modèle simple prédit que le flux surfacique d'énergie reçue par la banquise est de la forme

où

est la température de surface de la banquise.

Les paramètres

peuvent prendre deux valeurs suivant la saison. On admettra qu'il n'existe que deux saisons appelées saison chaude

et saison froide

. Chacune dure six mois. On ne prend pas en compte la salinité de l'eau de mer; on considère que la banquise gèle à

et qu'elle est entièrement caractérisée par sa température de surface

(en contact avec l'air) et son épaisseur

(figure 3). On se place dans l'approximation quasi-stationnaire (cf. I.3.) et l'on note

la durée d'une demie année.

Figure 3

La saison froide

Au début de la saison froide, la banquise a une épaisseur

et une température uniforme de

Première phase. On admet que la banquise ne fait que se refroidir et son épaisseur reste égale à . On modélise à tout instant la distribution de température de à dans la banquise par une loi linéaire.
a) Quelle quantité d'énergie (par unité de surface) doit-on fournir à la banquise quand sa température de surface change de ?
b) En déduire l'équation d'évolution de .
c) On admet que cette phase dure tant que le flux thermique dans la banquise calculé dans ce modèle reste inférieur au flux surfacique. À quelle température a-t-on égalité des densités de courant thermique en surface?
d) Au bout de quelle durée la température de surface a-t-elle atteint ?

On introduira

.
e) Application numérique : calculer

puis

pour une épaisseur initiale

et 5 m .
2. Deuxième phase : la couche de glace se met à croître ; on suppose toujours la distribution de température linéaire mais avec un flux thermique égal à celui imposé à la surface.
a) Exprimer la densité de courant thermique au sein de la banquise en fonction de

et de

. Montrer que la condition d'égalité des flux détermine

en fonction de

; exprimer

en fonction de

à l'aide de

et de

.
b) Reprendre dans le cadre de ce modèle l'équation (1).

On pose

. Montrer que, pour

c) Application numérique : calculer l'épaisseur de la banquise

et sa température de surface

à la fin de la saison froide pour une épaisseur initiale

et 5 m .

La saison chaude

L'apparition du soleil change le bilan thermique au niveau de la surface de la banquise. Il devient positif et la banquise va se réchauffer avant de fondre en surface.
3. Troisième phase. La banquise commence d'abord par se réchauffer jusqu'à ce que toute sa masse atteigne

. On adopte le même modèle qu'en III.1.
a) On pose

. Montrer que la durée

de ce réchauffement s'écrit

.
b) Application numérique : calculer

pour une épaisseur initiale

et 5 m .
4. Quatrième phase. La banquise fond par sa surface au contact de l'air.
a) Montrer que, dans cette étape, l'épaisseur de la banquise décroît linéairement avec le temps.
b) Application numérique : calculer les épaisseurs obtenues à la fin de la saison chaude pour un couvert en début de saison froide

et 5 m .

À partir des résultats numériques, montrer que ce modèle rend plausible pour

l'existence d'une solution périodique de période un an.