LMDZPedia - Contributions de l’utilisateur [fr]

Eau Condensee Fisrtilp

2021-09-14T11:20:16Z

Jyg :

$$\newcommand{\Condl}{Conditions aux limites }
\newcommand{\condl}{conditions aux limites }
\newcommand{\cad}{c'est-\`a-dire }
\newcommand{\Cad}{C'est-\`a-dire }
\newcommand{\tem}{temp\'erature }
\newcommand{\tems}{temp\'eratures }
\newcommand{\Tpni}{T_{\text{p,ni}}}
\newcommand{\Tp}{T_{\text{p}}}
\newcommand{\qpni}{q_{\text{p,ni}}}
\newcommand{\qp}{q_{\text{p}}}
\newcommand{\qsati}{q_{\text{sat,i}}}
\newcommand{\qsat}{q_{\text{sat}}}
\newcommand{\hpni}{h_{\text{p,ni}}}
\newcommand{\hp}{h_{\text{p}}}
%*KEEP,mathwake,T=LATEX.
\newcommand{\nab}{\vec{\nabla}}
\newcommand{\V}{\vec{V}}
\newcommand{\Vt}{\vec{\tilde{V}}}
\newcommand{\F}{\vec{F}}
\newcommand{\NF}{\|\vec{F}\|}
\newcommand{\CD}{C_D}
\newcommand{\pt}{\tilde{p}}
\newcommand{\Dt}{\tilde{D}}
\newcommand{\thetat}{\tilde{\theta}}
\newcommand{\alphat}{\tilde{\alpha}}
\newcommand{\Mt}{\tilde{M}}
\newcommand{\rhot}{\tilde{\rho}}
\newcommand{\qt}{\tilde{q}}
\newcommand{\Tt}{\tilde{T}}
\newcommand{\wt}{\tilde{w}}
\newcommand{\wbar}{\overline{w}}
\newcommand{\ubar}{\overline{u}}
\newcommand{\vbar}{\overline{v}}
\newcommand{\qbar}{\overline{q_v}}
\newcommand{\omgb}{\overline{\omega}}
\newcommand{\omgwb}{\overline{\omega_w}}
\newcommand{\omgub}{\overline{\omega_u}}
\newcommand{\omgxb}{\overline{\omega_x}}
\newcommand{\thetab}{\overline{\theta}}
\newcommand{\qvb}{\overline{q_v}}
\newcommand{\qvwb}{\overline{q_{vw}}}
\newcommand{\qvxb}{\overline{q_{vx}}}
\newcommand{\qvw}{q_{vw}}
\newcommand{\qvx}{q_{vx}}
\newcommand{\thetawb}{\overline{\theta_w}}
\newcommand{\thetaub}{\overline{\theta_u}}
\newcommand{\thetaxb}{\overline{\theta_x}}
\newcommand{\uwb}{\overline{u_w}}
\newcommand{\uub}{\overline{u_u}}
\newcommand{\uxb}{\overline{u_x}}
\newcommand{\vwb}{\overline{v_w}}
\newcommand{\vub}{\overline{v_u}}
\newcommand{\vxb}{\overline{v_x}}
\newcommand{\vwhat}{\hat{v_w}}
\newcommand{\vuhat}{\hat{v_u}}
\newcommand{\vxhat}{\hat{v_x}}
\newcommand{\Kb}{\overline{K}}
\newcommand{\alphab}{\overline{\alpha}}
\newcommand{\alphawb}{\overline{\alpha_w}}
\newcommand{\alphaub}{\overline{\alpha_u}}
\newcommand{\alphaxb}{\overline{\alpha_x}}
\newcommand{\rhob}{\overline{\rho}}
\newcommand{\sigw}{\ensuremath{\sigma_w}}
\newcommand{\dlt}{\delta t}
\newcommand{\Xrz}{X_{r}}
\newcommand{\Xr}[1]{X_{r(#1)}}
\newcommand{\Xrpz}{X_{rp}}
\newcommand{\Xrp}[1]{X_{rp(#1)}}
\newcommand{\Xrvz}{X_{rv}}
\newcommand{\Xrv}[1]{X_{rv(#1)}}
\newcommand{\Xrhz}{\hat{X}_{r}}
\newcommand{\Xrh}[1]{\hat{X}_{r(#1)}}
%%%\newcommand{\dlts2}{\frac{\dlt}{2}}
\newcommand{\dsdlt}{\frac{2}{\dlt}}
\newcommand{\vphi}{\vec{\varphi}}
\newcommand{\veta}{\vec{\eta}}
\newcommand{\vtheta}{\vec{\theta}}
\newcommand{\veps}{\vec{\varepsilon}}
\newcommand{\delphi}{\vec{\delta \varphi}}
\newcommand{\vG}{\vec{G}}
\newcommand{\vf}{\vec{f}}
\newcommand{\n}{\vec{n}}
\newcommand{\rl}{\rho_{l}}
\newcommand{\rv}{\rho_{v}}
\newcommand{\vv}{\vec{V}_{v}}
\newcommand{\vl}{\vec{V}_{l}}
%%%\newcommand{\v}{\vec{V}}
\newcommand{\g}{\vec{g}}
\newcommand{\dsdt}{\frac {\partial }{\partial t}}
\newcommand{\dst}{\partial_t}
\newcommand{\dsp}{\partial_p}
\newcommand{\dsz}{\partial_z}
\newcommand{\dsy}{\partial_y}
\newcommand{\dsx}{\partial_x}
\newcommand{\dsth}{\partial_{\theta}}
\newcommand{\dsq}{\partial_q}
\newcommand{\dsv}{\partial_V}
\newcommand{\dalpha}{\delta \alpha}
\newcommand{\dw}{\delta w}
\newcommand{\dqv}{\delta q_v}
\newcommand{\dtheta}{\delta \theta}
\newcommand{\domg}{\delta \omega}
\newcommand{\domgcv}{\delta \omega^{\text{cv}}}
%%\newcommand{\domgtop}{\delta \omega_{\text{top}}}
\newcommand{\domgtop}{\delta \omega_{w}}
%%\newcommand{\domgupper}{\delta \omega_{\text{upper}}}
\newcommand{\domgupper}{\delta \omega_{h}}
%%\newcommand{\ptop}{p_{\text{top}}}
\newcommand{\ptop}{p_{w}}
%%\newcommand{\pupper}{p_{\text{upper}}}
\newcommand{\pupper}{p_{h}}
%%\newcommand{\htop}{h_{\text{top}}}
\newcommand{\htop}{h_{w}}
%%\newcommand{\hupper}{h_{\text{upper}}}
\newcommand{\hupper}{h_{h}}
\newcommand{\intw}{\int_{\Sigma_w}}
\newcommand{\intu}{\int_{\Sigma_u}}
\newcommand{\intx}{\int_{\Sigma_x}}
\newcommand{\inte}{\int_{\Sigma_e}}
\newcommand{\intp}{\int_{\Sigma'}}
\newcommand{\intgw}{\int_{\Gamma_{w,in}}}
\newcommand{\intgwp}{\int_{\Gamma_{w,in}^+}}
\newcommand{\intgwm}{\int_{\Gamma_{w,in}^-}}
\newcommand{\intg}{\int_{\Gamma}}
\newcommand{\intgp}{\int_{\Gamma'}}
\newcommand{\intgpe}{\int_{\Gamma'_e}}
\newcommand{\intgpu}{\int_{\Gamma'_u}}
\newcommand{\intgpx}{\int_{\Gamma'_x}}
\newcommand{\intgpw}{\int_{\Gamma'_w}}
\newcommand{\kgw}{k_{\text{gw}}}
\newcommand{\taugw}{\tau_{\text{gw}}}
\newcommand{\klift}{k_{\text{lift}}}
%%\newcommand{\kliftw}{k_{\text{lift,wk}}}
\newcommand{\kliftw}{\mathcal{E}_{\text{lift}}^{\text{wk}}}
\newcommand{\Eliftw}{\mathcal{E}_{\text{lift}}^{\text{wk}}}
\newcommand{\Pb}{\ensuremath{P_{\text{buoy}}}}
\newcommand{\Pliftw}{\ensuremath{P_{\text{lift}}^{\text{wk}}}}
\newcommand{\Pliftwh}{\ensuremath{\widehat{P_{\text{lift}}^{\text{wk}}}}}
\newcommand{\Pliftwp}{\ensuremath{\widetilde{P_{\text{lift}}^{\text{wk}}}}}
\newcommand{\LG}{\ensuremath{L_{\Gamma}}}
\newcommand{\LGh}{\ensuremath{\widehat{L_{\Gamma}}}}
\newcommand{\LGp}{\ensuremath{\widetilde{L_{\Gamma}}}}
\newcommand{\LGt}{\ensuremath{\widetilde{L_{\Gamma}}}}
\newcommand{\Dw}{\ensuremath{D_{\text{wk}}}}
\newcommand{\QI}{\ensuremath{Q_1}}
\newcommand{\DQI}{\ensuremath{\delta Q_1}}
\newcommand{\QIsat}{\ensuremath{Q_{1}^{\text{sat}}}}
\newcommand{\QIcvsat}{\ensuremath{Q_{1}^{\text{cv,sat}}}}
\newcommand{\QIunsat}{\ensuremath{Q_{1}^{\text{unsat}}}}
\newcommand{\QIcvunsat}{\ensuremath{Q_{1}^{\text{cv,unsat}}}}
\newcommand{\QIIsat}{\ensuremath{Q_{2}^{\text{sat}}}}
\newcommand{\QIIcvsat}{\ensuremath{Q_{2}^{\text{cv,sat}}}}
\newcommand{\QIIunsat}{\ensuremath{Q_{2}^{\text{unsat}}}}
\newcommand{\QIIcvunsat}{\ensuremath{Q_{2}^{\text{cv,unsat}}}}
\newcommand{\QIw}{\ensuremath{Q_{1,w}}}
\newcommand{\QIcvw}{\ensuremath{Q_{1,w}^{\text{cv}}}}
\newcommand{\QIx}{\ensuremath{Q_{1,x}}}
\newcommand{\QIcvx}{\ensuremath{Q_{1,x}^{\text{cv}}}}
\newcommand{\QIcv}{\ensuremath{Q_{1}^{\text{cv}}}}
\newcommand{\DQIcv}{\ensuremath{\delta Q_{1}^{\text{cv}}}}
\newcommand{\QIbl}{\ensuremath{Q_{1}^{\text{bl}}}}
\newcommand{\QIcl}{\ensuremath{Q_{1}^{\text{cl}}}}
\newcommand{\QIwk}{\ensuremath{Q_{1}^{\text{wk}}}}
\newcommand{\DQIwk}{\ensuremath{\delta Q_{1}^{\text{wk}}}}
\newcommand{\QIp}{\ensuremath{Q_{1,p}}}
\newcommand{\QIS}{\ensuremath{\mathcal{Q}_{1}^S}}
\newcommand{\QID}{\ensuremath{\mathcal{Q}_{1}^D}}
\newcommand{\QIT}{\ensuremath{\mathcal{Q}_{1}^T}}
\newcommand{\QIE}{\ensuremath{\mathcal{Q}_{1}^E}}
\newcommand{\QII}{\ensuremath{Q_2}}
\newcommand{\DQII}{\ensuremath{\delta Q_2}}
\newcommand{\QIIw}{\ensuremath{Q_{2,w}}}
\newcommand{\QIIcvw}{\ensuremath{Q_{2,w}^{\text{cv}}}}
\newcommand{\QIIx}{\ensuremath{Q_{2,x}}}
\newcommand{\QIIcvx}{\ensuremath{Q_{2,x}^{\text{cv}}}}
\newcommand{\QIIcv}{\ensuremath{Q_{2}^{\text{cv}}}}
\newcommand{\DQIIcv}{\ensuremath{\delta Q_{2}^{\text{cv}}}}
\newcommand{\QIIbl}{\ensuremath{Q_{2}^{\text{bl}}}}
\newcommand{\QIIcl}{\ensuremath{Q_{2}^{\text{cl}}}}
\newcommand{\QIIwk}{\ensuremath{Q_{2}^{\text{wk}}}}
\newcommand{\DQIIwk}{\ensuremath{\delta Q_{2}^{\text{wk}}}}
\newcommand{\QIIp}{\ensuremath{Q_{2,p}}}
\newcommand{\QIIS}{\ensuremath{\mathcal{Q}_{2}^S}}
\newcommand{\QIID}{\ensuremath{\mathcal{Q}_{2}^D}}
\newcommand{\QIIT}{\ensuremath{\mathcal{Q}_{2}^T}}
\newcommand{\QIIE}{\ensuremath{\mathcal{Q}_{2}^E}}
\newcommand{\Sig}{\ensuremath{\Sigma\;}}
\newcommand{\Sigp}{\ensuremath{\Sigma'\;}}
\newcommand{\Sigw}{\ensuremath{\Sigma_w\;}}
\newcommand{\Sige}{\ensuremath{\Sigma_e}}
\newcommand{\Sigu}{\ensuremath{\Sigma_u}}
\newcommand{\Sigx}{\ensuremath{\Sigma_x}}
\newcommand{\Wb}{\ensuremath{W_{\text{buoy}}}}
\newcommand{\Wbb}{\bar{\bar{W}}}
\newcommand{\wa}{\langle w\rangle}
\newcommand{\wpa}{\langle w'\rangle}
\newcommand{\wwa}{\langle w^2\rangle}
\newcommand{\wwwa}{\langle w^3\rangle}
\newcommand{\wwpa}{\langle w^{\prime\prime\prime 2}\rangle}$$

Le but de cette note est d'établir les formules menant à la détermination des quantités
d'eau condensées sous forme d'eau liquide et de glace dans le sous-programme Fisrtilp. La
principale difficulté réside dans le fait que la condensation dégage de la chaleur latente
et que $$\qsat$$ dépend de la température. Dans tous ces calculs, on considère une couche
atmosphérique sans condensat de température $$T$$ et d'humidité spécifique $$q$$ et on se
donne la fraction $$f$$ de la maille où se produit la condensation. Le principe est d'écrire
que l'énergie statique humide et la quantité d'eau totale sont conservées dans le
processus de condensation.

Dans la première version créée par Laurent Li, on suppose que seule la partie
nuageuse de la maille voit sa température augmenter du fait de la condensation. Alors le
système que l'on décrit présente des températures différentes dans le nuage et hors du
nuage.

Laurent Li résoud :
\begin{equation}
C_p T' \;+\; L_v \qsat(T') \;=\; C_p T \;+\; L_v q
\label{eq.init}
\end{equation}
étant entendu que cela s'applique à la partie nuageuse de la maille.

Nous voulons, au contraire, écrire que les températures sont toujours uniformisées
horizontalement très rapidement (et même on va supposer instantanément). En outre nous
allons généraliser les calculs au cas où de la condensation solide est présente et nous
allons profiter de l'occasion pour écrire des formules complètes prenant en compte la
variation des chaleurs latentes avec la température (voir dans l'annexe les formules donnant
l'énergie statique humide $$h$$).

== Écriture des équations de conservation ==

Soit donc $$f$$ la fraction nuageuse. On écrit que sur $$f$$ il y a condensation et que sur
$$(1-f)$$ il n'y a pas de changement de phase. On suppose connu le contenu massique total
d'eau $$q_n$$ (i.e. vapeur + liquide + glace) dans le nuage. En désignant par $$T'$$ la
température (uniforme) après condensation et par $$q"$$ l'humidité spécifique hors des
nuages, les énergies statiques massiques initiale $$h$$, dans le nuage $$h_n$$ et hors du
nuage $$h"$$ s'écrivent :
\begin{equation}
\left \lbrace
\begin{array}{ll}
h_n &=\;[(1-q_n)C_{pd} + q_n C_l] T' + L_v \qsat - L_f q_i + gz \\
h" &=\;[(1-q")C_{pd} + q"C_l] T' + L_v q" + gz \\
h &=\;[(1-q)C_{pd} + q C_l] T + L_v q + gz \\
\end{array}
\right .
\end{equation}
où $$\qsat$$ désigne l'humidité saturante. Ce peut être l'humidité saturante à la
température $$T'$$ sur l'eau liquide ou sur la glace. A voir plus loin.

La conservation de l'énergie statique humide et de l'eau s'écrit :
\begin{equation}
\left \lbrace
\begin{array}{ll}
f h_n + (1-f)h" &=\;h \\
f q_n + (1-f)q" &=\;q \\
\end{array}
\right .
\end{equation}
Ce qui donne :
\begin{equation}
\left \lbrace
\begin{array}{l}
T'[C_{pd} + (f q_n + (1-f) q") (C_l-C_{pd})] + L_v [f \qsat + (1-f)q"] - L_f f q_i \\
\qquad \qquad = [C_{pd} + q (C_l-C_{pd})] T + L_v q \\
\\
q = f q_n + (1-f) q"\\
\end{array}
\right .
\end{equation}
Soit, en posant $$C_p = C_{pd} + (C_l-C_{pd}) q$$ :
\begin{equation}
\left \lbrace
\begin{array}{l}
C_p T' - f [L_v(q_n-\qsat) + L_f q_i] = C_p T \\
q = f q_n + (1-f) q"\\
\end{array}
\right .
\end{equation}

Maintenant nous supposons que la glace représente une fraction $$\gamma$$ du condensat :
$$q_i = \gamma (q_n - \qsat)$$. Finalement les équations de conservation s'écrivent :
\begin{equation}
\label{eq.conserv_fin}
\left \lbrace
\begin{array}{l}
C_p T' - f (L_v+\gamma L_f)(q_n-\qsat) = C_p T \\
q = f q_n + (1-f) q"\\
\end{array}
\right .
\end{equation}

== Mise en oeuvre ==
Nous supposons maintenant que $$\qsat$$ est égale à l'humidité saturante sur eau liquide
lorsqu'il n'y a pas de glace et à l'humidité saturante sur glace dès que la glace est
présente.

L'équation (\ref{eq.conserv_fin}) pour la température peut être écrite sous une forme
analogue à l'équation initiale (\ref{eq.init}) de Laurent Li :
\begin{equation}
\label{eq.conserv_fin_lli}
C_p T' \;+\; f (L_v+\gamma L_f) \qsat(T')] \;=\; C_p T \;+\; f (L_v+\gamma L_f) q_n
\end{equation}

On constate que l'on passe de l'équation (\ref{eq.init}) à l'équation
(\ref{eq.conserv_fin_lli}) en remplaçant $$L_v$$ par $$f (L_v+\gamma L_f)$$ et $$q$$ par $$q_n$$.
L'itération convergeant vers $$T'$$ s'obtient ainsi directement à partir des formules de
Laurent :
\begin{equation}
\displaystyle{
T'_{j+1}\;=\;T'_j\;+\;\frac{T-T'_j\;+\;\frac{f (L_v+\gamma L_f)}{C_p}[q_n\;-\;\qsat(T'_j)]}
{1\;+\;\frac{f (L_v+\gamma L_f)}{C_p}\partial_T\qsat(T'_j)} }
\label{eq.iter}
\end{equation}

L'expérience montre que cette itération converge assez vite, en 3 ou 4 passages. Elle va
fournir une valeur approchée $$T'_a$$ à partir de laquelle on peut calculer l'eau condensée :
$$
c\;=\;f[q_n\;-\;\qsat(T'_a)]
$$
et les estimations finales de la température $$T'$$ et de l'humidité $$q'$$ telles que
l'énergie et l'eau soient exactement conservées :
\begin{equation}
\left \lbrace
\begin{array}{l}
T'\;=\; T\;+\;\frac{L_v+\gamma L_f}{C_p}c \\
q'\;=\;q\;-\;c \\
\end{array}
\right .
\end{equation}

== Annexe: définitions des enthalpies et énergies statiques massiques ==
=== I - Enthalpies massiques: ===
Grandeurs conservées dans les transformations adiabatiques
isobares.
\begin{equation}
\begin{array}{ll}
\left \lbrace
\begin{array}{lll}
k &=& (1-q_t) k_d + q_v k_v + q_l k_l + q_i k_i \\
L_v(T) &=& k_v - k_l \\
L_s(T) &=& k_v - k_i\;, \quad L_f(T) = L_s(T)-L_v(T)\\
k_d &=& C_{pd} T \\
k_l &=& C_l T \\
k_i &=& C_i T \\
\end{array}
\right \rbrace & \left .
\begin{array}{l}
\rm{enthalpie} \\
k=[(1-q_t) C_{pd} +q_t C_l] T + L_v q_v - L_f q_i\\
\end{array} \right .
\end{array}
\end{equation}

Comme $$L_v-(C_{pv}-C_l)T = \rm{cste}$$,
l'enthalpie avec eau liquide $$k_w=k-q_t[L_v-(C_{pv}-C_l)T]$$
est aussi conservée dans les transformations adiabatiques isobares :
$$
k_w = [(1-q_t) C_{pd} + q_t C_{pv}] T - L_v q_l - L_s q_i
$$
On peut faire le même genre de chose pour obtenir une expression sans $$q_i$$,
ce qui donne : et comme $$L_f-(C_l-C_i)T = \rm{cste}$$, l'enthalpie (sans nom)
$$k_i=k-q_t[L_f-(C_l-C_i)T]$$ est aussi conservée dans les transformations
adiabatiques isobares :
$$
k_i = [(1-q_t) C_{pd} + q_t C_i] T+ L_s q_v + L_f q_l
$$

=== II - Energies statiques: ===
Grandeurs conservées dans les
transformations adiabatiques où le changement
de pression est purement hydrostatique.
\begin{equation}
\begin{array}{ll}
\left .
\begin{array}{lll}
\rm{Adiabatisme :} & dk &\displaystyle{=\frac{dp}{\rho} }\\
& & \\
\rm{Hydrostasie :} & \displaystyle{\frac{dp}{\rho} }&= -g dz
\end{array}
\right \rbrace & \left .
\begin{array}{l}
dk + g dz = 0 \\
\end{array} \right .
\end{array}
\end{equation}

{| border="1" style="border-collapse:collapse"
| Energie statique humide :
| $$h=[(1-q_t)C_{pd}+q_tC_l]T + L_v q_v - L_f q_i + gz$$
|-
| Energie statique avec eau liquide :
| $$h_w=[(1-q_t)C_{pd}+q_tC_{pv}]T - L_vq_l - L_s q_i + gz$$
|-
| Energie statique sans nom :
| $$h_i=[(1-q_t) C_{pd} + q_t C_i] T+ L_s q_v + L_f q_l + gz$$
|-
| Energie statique sèche :
| $$h_d=[(1-q)C_{pd}+qC_{pv}]T +gz$$ 
(= $$h_w$$ avec $$q_l\;=\;0$$ et $$q_i\;=\;0$$)
|}

[[Category:Physique]]
[[Category:Notes Techniques]]

Eau Condensee Fisrtilp

2021-09-13T15:10:23Z

Jyg :

$$\newcommand{\Condl}{Conditions aux limites }
\newcommand{\condl}{conditions aux limites }
\newcommand{\cad}{c'est-\`a-dire }
\newcommand{\Cad}{C'est-\`a-dire }
\newcommand{\tem}{temp\'erature }
\newcommand{\tems}{temp\'eratures }
\newcommand{\Tpni}{T_{\text{p,ni}}}
\newcommand{\Tp}{T_{\text{p}}}
\newcommand{\qpni}{q_{\text{p,ni}}}
\newcommand{\qp}{q_{\text{p}}}
\newcommand{\qsati}{q_{\text{sat,i}}}
\newcommand{\qsat}{q_{\text{sat}}}
\newcommand{\hpni}{h_{\text{p,ni}}}
\newcommand{\hp}{h_{\text{p}}}
%*KEEP,mathwake,T=LATEX.
\newcommand{\nab}{\vec{\nabla}}
\newcommand{\V}{\vec{V}}
\newcommand{\Vt}{\vec{\tilde{V}}}
\newcommand{\F}{\vec{F}}
\newcommand{\NF}{\|\vec{F}\|}
\newcommand{\CD}{C_D}
\newcommand{\pt}{\tilde{p}}
\newcommand{\Dt}{\tilde{D}}
\newcommand{\thetat}{\tilde{\theta}}
\newcommand{\alphat}{\tilde{\alpha}}
\newcommand{\Mt}{\tilde{M}}
\newcommand{\rhot}{\tilde{\rho}}
\newcommand{\qt}{\tilde{q}}
\newcommand{\Tt}{\tilde{T}}
\newcommand{\wt}{\tilde{w}}
\newcommand{\wbar}{\overline{w}}
\newcommand{\ubar}{\overline{u}}
\newcommand{\vbar}{\overline{v}}
\newcommand{\qbar}{\overline{q_v}}
\newcommand{\omgb}{\overline{\omega}}
\newcommand{\omgwb}{\overline{\omega_w}}
\newcommand{\omgub}{\overline{\omega_u}}
\newcommand{\omgxb}{\overline{\omega_x}}
\newcommand{\thetab}{\overline{\theta}}
\newcommand{\qvb}{\overline{q_v}}
\newcommand{\qvwb}{\overline{q_{vw}}}
\newcommand{\qvxb}{\overline{q_{vx}}}
\newcommand{\qvw}{q_{vw}}
\newcommand{\qvx}{q_{vx}}
\newcommand{\thetawb}{\overline{\theta_w}}
\newcommand{\thetaub}{\overline{\theta_u}}
\newcommand{\thetaxb}{\overline{\theta_x}}
\newcommand{\uwb}{\overline{u_w}}
\newcommand{\uub}{\overline{u_u}}
\newcommand{\uxb}{\overline{u_x}}
\newcommand{\vwb}{\overline{v_w}}
\newcommand{\vub}{\overline{v_u}}
\newcommand{\vxb}{\overline{v_x}}
\newcommand{\vwhat}{\hat{v_w}}
\newcommand{\vuhat}{\hat{v_u}}
\newcommand{\vxhat}{\hat{v_x}}
\newcommand{\Kb}{\overline{K}}
\newcommand{\alphab}{\overline{\alpha}}
\newcommand{\alphawb}{\overline{\alpha_w}}
\newcommand{\alphaub}{\overline{\alpha_u}}
\newcommand{\alphaxb}{\overline{\alpha_x}}
\newcommand{\rhob}{\overline{\rho}}
\newcommand{\sigw}{\ensuremath{\sigma_w}}
\newcommand{\dlt}{\delta t}
\newcommand{\Xrz}{X_{r}}
\newcommand{\Xr}[1]{X_{r(#1)}}
\newcommand{\Xrpz}{X_{rp}}
\newcommand{\Xrp}[1]{X_{rp(#1)}}
\newcommand{\Xrvz}{X_{rv}}
\newcommand{\Xrv}[1]{X_{rv(#1)}}
\newcommand{\Xrhz}{\hat{X}_{r}}
\newcommand{\Xrh}[1]{\hat{X}_{r(#1)}}
%%%\newcommand{\dlts2}{\frac{\dlt}{2}}
\newcommand{\dsdlt}{\frac{2}{\dlt}}
\newcommand{\vphi}{\vec{\varphi}}
\newcommand{\veta}{\vec{\eta}}
\newcommand{\vtheta}{\vec{\theta}}
\newcommand{\veps}{\vec{\varepsilon}}
\newcommand{\delphi}{\vec{\delta \varphi}}
\newcommand{\vG}{\vec{G}}
\newcommand{\vf}{\vec{f}}
\newcommand{\n}{\vec{n}}
\newcommand{\rl}{\rho_{l}}
\newcommand{\rv}{\rho_{v}}
\newcommand{\vv}{\vec{V}_{v}}
\newcommand{\vl}{\vec{V}_{l}}
%%%\newcommand{\v}{\vec{V}}
\newcommand{\g}{\vec{g}}
\newcommand{\dsdt}{\frac {\partial }{\partial t}}
\newcommand{\dst}{\partial_t}
\newcommand{\dsp}{\partial_p}
\newcommand{\dsz}{\partial_z}
\newcommand{\dsy}{\partial_y}
\newcommand{\dsx}{\partial_x}
\newcommand{\dsth}{\partial_{\theta}}
\newcommand{\dsq}{\partial_q}
\newcommand{\dsv}{\partial_V}
\newcommand{\dalpha}{\delta \alpha}
\newcommand{\dw}{\delta w}
\newcommand{\dqv}{\delta q_v}
\newcommand{\dtheta}{\delta \theta}
\newcommand{\domg}{\delta \omega}
\newcommand{\domgcv}{\delta \omega^{\text{cv}}}
%%\newcommand{\domgtop}{\delta \omega_{\text{top}}}
\newcommand{\domgtop}{\delta \omega_{w}}
%%\newcommand{\domgupper}{\delta \omega_{\text{upper}}}
\newcommand{\domgupper}{\delta \omega_{h}}
%%\newcommand{\ptop}{p_{\text{top}}}
\newcommand{\ptop}{p_{w}}
%%\newcommand{\pupper}{p_{\text{upper}}}
\newcommand{\pupper}{p_{h}}
%%\newcommand{\htop}{h_{\text{top}}}
\newcommand{\htop}{h_{w}}
%%\newcommand{\hupper}{h_{\text{upper}}}
\newcommand{\hupper}{h_{h}}
\newcommand{\intw}{\int_{\Sigma_w}}
\newcommand{\intu}{\int_{\Sigma_u}}
\newcommand{\intx}{\int_{\Sigma_x}}
\newcommand{\inte}{\int_{\Sigma_e}}
\newcommand{\intp}{\int_{\Sigma'}}
\newcommand{\intgw}{\int_{\Gamma_{w,in}}}
\newcommand{\intgwp}{\int_{\Gamma_{w,in}^+}}
\newcommand{\intgwm}{\int_{\Gamma_{w,in}^-}}
\newcommand{\intg}{\int_{\Gamma}}
\newcommand{\intgp}{\int_{\Gamma'}}
\newcommand{\intgpe}{\int_{\Gamma'_e}}
\newcommand{\intgpu}{\int_{\Gamma'_u}}
\newcommand{\intgpx}{\int_{\Gamma'_x}}
\newcommand{\intgpw}{\int_{\Gamma'_w}}
\newcommand{\kgw}{k_{\text{gw}}}
\newcommand{\taugw}{\tau_{\text{gw}}}
\newcommand{\klift}{k_{\text{lift}}}
%%\newcommand{\kliftw}{k_{\text{lift,wk}}}
\newcommand{\kliftw}{\mathcal{E}_{\text{lift}}^{\text{wk}}}
\newcommand{\Eliftw}{\mathcal{E}_{\text{lift}}^{\text{wk}}}
\newcommand{\Pb}{\ensuremath{P_{\text{buoy}}}}
\newcommand{\Pliftw}{\ensuremath{P_{\text{lift}}^{\text{wk}}}}
\newcommand{\Pliftwh}{\ensuremath{\widehat{P_{\text{lift}}^{\text{wk}}}}}
\newcommand{\Pliftwp}{\ensuremath{\widetilde{P_{\text{lift}}^{\text{wk}}}}}
\newcommand{\LG}{\ensuremath{L_{\Gamma}}}
\newcommand{\LGh}{\ensuremath{\widehat{L_{\Gamma}}}}
\newcommand{\LGp}{\ensuremath{\widetilde{L_{\Gamma}}}}
\newcommand{\LGt}{\ensuremath{\widetilde{L_{\Gamma}}}}
\newcommand{\Dw}{\ensuremath{D_{\text{wk}}}}
\newcommand{\QI}{\ensuremath{Q_1}}
\newcommand{\DQI}{\ensuremath{\delta Q_1}}
\newcommand{\QIsat}{\ensuremath{Q_{1}^{\text{sat}}}}
\newcommand{\QIcvsat}{\ensuremath{Q_{1}^{\text{cv,sat}}}}
\newcommand{\QIunsat}{\ensuremath{Q_{1}^{\text{unsat}}}}
\newcommand{\QIcvunsat}{\ensuremath{Q_{1}^{\text{cv,unsat}}}}
\newcommand{\QIIsat}{\ensuremath{Q_{2}^{\text{sat}}}}
\newcommand{\QIIcvsat}{\ensuremath{Q_{2}^{\text{cv,sat}}}}
\newcommand{\QIIunsat}{\ensuremath{Q_{2}^{\text{unsat}}}}
\newcommand{\QIIcvunsat}{\ensuremath{Q_{2}^{\text{cv,unsat}}}}
\newcommand{\QIw}{\ensuremath{Q_{1,w}}}
\newcommand{\QIcvw}{\ensuremath{Q_{1,w}^{\text{cv}}}}
\newcommand{\QIx}{\ensuremath{Q_{1,x}}}
\newcommand{\QIcvx}{\ensuremath{Q_{1,x}^{\text{cv}}}}
\newcommand{\QIcv}{\ensuremath{Q_{1}^{\text{cv}}}}
\newcommand{\DQIcv}{\ensuremath{\delta Q_{1}^{\text{cv}}}}
\newcommand{\QIbl}{\ensuremath{Q_{1}^{\text{bl}}}}
\newcommand{\QIcl}{\ensuremath{Q_{1}^{\text{cl}}}}
\newcommand{\QIwk}{\ensuremath{Q_{1}^{\text{wk}}}}
\newcommand{\DQIwk}{\ensuremath{\delta Q_{1}^{\text{wk}}}}
\newcommand{\QIp}{\ensuremath{Q_{1,p}}}
\newcommand{\QIS}{\ensuremath{\mathcal{Q}_{1}^S}}
\newcommand{\QID}{\ensuremath{\mathcal{Q}_{1}^D}}
\newcommand{\QIT}{\ensuremath{\mathcal{Q}_{1}^T}}
\newcommand{\QIE}{\ensuremath{\mathcal{Q}_{1}^E}}
\newcommand{\QII}{\ensuremath{Q_2}}
\newcommand{\DQII}{\ensuremath{\delta Q_2}}
\newcommand{\QIIw}{\ensuremath{Q_{2,w}}}
\newcommand{\QIIcvw}{\ensuremath{Q_{2,w}^{\text{cv}}}}
\newcommand{\QIIx}{\ensuremath{Q_{2,x}}}
\newcommand{\QIIcvx}{\ensuremath{Q_{2,x}^{\text{cv}}}}
\newcommand{\QIIcv}{\ensuremath{Q_{2}^{\text{cv}}}}
\newcommand{\DQIIcv}{\ensuremath{\delta Q_{2}^{\text{cv}}}}
\newcommand{\QIIbl}{\ensuremath{Q_{2}^{\text{bl}}}}
\newcommand{\QIIcl}{\ensuremath{Q_{2}^{\text{cl}}}}
\newcommand{\QIIwk}{\ensuremath{Q_{2}^{\text{wk}}}}
\newcommand{\DQIIwk}{\ensuremath{\delta Q_{2}^{\text{wk}}}}
\newcommand{\QIIp}{\ensuremath{Q_{2,p}}}
\newcommand{\QIIS}{\ensuremath{\mathcal{Q}_{2}^S}}
\newcommand{\QIID}{\ensuremath{\mathcal{Q}_{2}^D}}
\newcommand{\QIIT}{\ensuremath{\mathcal{Q}_{2}^T}}
\newcommand{\QIIE}{\ensuremath{\mathcal{Q}_{2}^E}}
\newcommand{\Sig}{\ensuremath{\Sigma\;}}
\newcommand{\Sigp}{\ensuremath{\Sigma'\;}}
\newcommand{\Sigw}{\ensuremath{\Sigma_w\;}}
\newcommand{\Sige}{\ensuremath{\Sigma_e}}
\newcommand{\Sigu}{\ensuremath{\Sigma_u}}
\newcommand{\Sigx}{\ensuremath{\Sigma_x}}
\newcommand{\Wb}{\ensuremath{W_{\text{buoy}}}}
\newcommand{\Wbb}{\bar{\bar{W}}}
\newcommand{\wa}{\langle w\rangle}
\newcommand{\wpa}{\langle w'\rangle}
\newcommand{\wwa}{\langle w^2\rangle}
\newcommand{\wwwa}{\langle w^3\rangle}
\newcommand{\wwpa}{\langle w^{\prime\prime\prime 2}\rangle}$$

Le but de cette note est d'établir les formules menant à la détermination des quantités
d'eau condensées sous forme d'eau liquide et de glace dans le sous-programme Fisrtilp. La
principale difficulté réside dans le fait que la condensation dégage de la chaleur latente
et que $$\qsat$$ dépend de la température. Dans tous ces calculs, on considère une couche
atmosphérique sans condensat de température $$T$$ et d'humidité spécifique $$q$$ et on se
donne la fraction $$f$$ de la maille où se produit la condensation. Le principe est d'écrire
que l'énergie statique humide et la quantité d'eau totale sont conservées dans le
processus de condensation.

Dans la première version créée par Laurent Li, on suppose que seule la partie
nuageuse de la maille voit sa température augmenter du fait de la condensation. Alors le
système que l'on décrit présente des températures différentes dans le nuage et hors du
nuage.

Laurent Li résoud :
\begin{equation}
C_p T' \;+\; L_v \qsat(T') \;=\; C_p T \;+\; L_v q
\label{eq.init}
\end{equation}
étant entendu que cela s'applique à la partie nuageuse de la maille.

Nous voulons, au contraire, écrire que les températures sont toujours uniformisées
horizontalement très rapidement (et même on va supposer instantanément). En outre nous
allons généraliser les calculs au cas où de la condensation solide est présente et nous
allons profiter de l'occasion pour écrire des formules complètes prenant en compte la
variation des chaleurs latentes avec la température (voir dans l'annexe les formules donnant
l'énergie statique humide $$h$$).

== Écriture des équations de conservation ==

Soit donc $$f$$ la fraction nuageuse. On écrit que sur $$f$$ il y a condensation et que sur
$$(1-f)$$ il n'y a pas de changement de phase. On suppose connu le contenu massique total
d'eau $$q_n$$ (i.e. vapeur + liquide + glace) dans le nuage. En désignant par $$T'$$ la
température (uniforme) après condensation et par $$q"$$ l'humidité spécifique hors des
nuages, les énergies statiques massiques initiale $$h$$, dans le nuage $$h_n$$ et hors du
nuage $$h"$$ s'écrivent :
\begin{equation}
\left \lbrace
\begin{array}{ll}
h_n &=\;[(1-q_n)C_{pd} + q_n C_l] T' + L_v \qsat - L_f q_i + gz \\
h" &=\;[(1-q")C_{pd} + q"C_l] T' + L_v q" + gz \\
h &=\;[(1-q)C_{pd} + q C_l] T + L_v q + gz \\
\end{array}
\right .
\end{equation}
où $$\qsat$$ désigne l'humidité saturante. Ce peut être l'humidité saturante à la
température $$T'$$ sur l'eau liquide ou sur la glace. A voir plus loin.

La conservation de l'énergie statique humide et de l'eau s'écrit :
\begin{equation}
\left \lbrace
\begin{array}{ll}
f h_n + (1-f)h" &=\;h \\
f q_n + (1-f)q" &=\;q \\
\end{array}
\right .
\end{equation}
Ce qui donne :
\begin{equation}
\left \lbrace
\begin{array}{l}
T'[C_{pd} + (f q_n + (1-f) q") (C_l-C_{pd})] + L_v [f \qsat + (1-f)q"] - L_f f q_i \\
\qquad \qquad = [C_{pd} + q (C_l-C_{pd})] T + L_v q \\
\\
q = f q_n + (1-f) q"\\
\end{array}
\right .
\end{equation}
Soit, en posant $$C_p = C_{pd} + (C_l-C_{pd}) q$$ :
\begin{equation}
\left \lbrace
\begin{array}{l}
C_p T' - f [L_v(q_n-\qsat) + L_f q_i] = C_p T \\
q = f q_n + (1-f) q"\\
\end{array}
\right .
\end{equation}

Maintenant nous supposons que la glace représente une fraction $$\gamma$$ du condensat :
$$q_i = \gamma (q_n - \qsat)$$. Finalement les équations de conservation s'écrivent :
\begin{equation}
\label{eq.conserv_fin}
\left \lbrace
\begin{array}{l}
C_p T' - f (L_v+\gamma L_f)(q_n-\qsat) = C_p T \\
q = f q_n + (1-f) q"\\
\end{array}
\right .
\end{equation}

== Mise en oeuvre ==
Nous supposons maintenant que $$\qsat$$ est égale à l'humidité saturante sur eau liquide
lorsqu'il n'y a pas de glace et à l'humidité saturante sur glace dès que la glace est
présente.

L'équation (\ref{eq.conserv_fin}) pour la température peut être écrite sous une forme
analogue à l'équation initiale (\ref{eq.init}) de Laurent Li :
\begin{equation}
\label{eq.conserv_fin_lli}
C_p T' \;+\; f (L_v+\gamma L_f) \qsat(T')] \;=\; C_p T \;+\; f (L_v+\gamma L_f) q_n
\end{equation}

On constate que l'on passe de l'équation (\ref{eq.init}) à l'équation
(\ref{eq.conserv_fin_lli}) en remplaçant $$L_v$$ par $$f (L_v+\gamma L_f)$$ et $$q$$ par $$q_n$$.
L'itération convergeant vers $$T'$$ s'obtient ainsi directement à partir des formules de
Laurent :
\begin{equation}
\displaystyle{
T'_{j+1}\;=\;T'_j\;+\;\frac{T-T'_j\;+\;\frac{f (L_v+\gamma L_f)}{C_p}[q_n\;-\;\qsat(T'_j)]}
{1\;+\;\frac{f (L_v+\gamma L_f)}{C_p}\partial_T\qsat(T'_j)} }
\label{eq.iter}
\end{equation}

L'expérience montre que cette itération converge assez vite, en 3 ou 4 passages. Elle va
fournir une valeur approchée $$T'_a$$ à partir de laquelle on peut calculer l'eau condensée :
$$
c\;=\;f[q_n\;-\;\qsat(T'_a)]
$$
et les estimations finales de la température $$T'$$ et de l'humidité $$q'$$ telles que
l'énergie et l'eau soient exactement conservées :
\begin{equation}
\left \lbrace
\begin{array}{l}
T'\;=\; T\;+\;\frac{L_v+\gamma L_f}{C_p}c \\
q'\;=\;q\;-\;c \\
\end{array}
\right .
\end{equation}

== Annexe: définitions des enthalpies et énergies statiques massiques ==
=== I - Enthalpies massiques: ===
Grandeurs conservées dans les transformations adiabatiques
isobares.
\begin{equation}
\begin{array}{ll}
\left \lbrace
\begin{array}{lll}
k &=& (1-q_t) k_d + q_v k_v + q_l k_l + q_i k_i \\
L_v(T) &=& k_v - k_l \\
L_s(T) &=& k_v - k_i\;, \quad L_f(T) = L_s(T)-L_v(T)\\
k_d &=& C_{pd} T \\
k_l &=& C_l T \\
k_i &=& C_i T \\
\end{array}
\right \rbrace & \left .
\begin{array}{l}
\rm{enthalpie} \\
k=[(1-q_t) C_{pd} +q_t C_l] T + L_v q_v - L_f q_i\\
\end{array} \right .
\end{array}
\end{equation}

Comme $$L_v-(C_{pv}-C_l)T = \rm{cste}$$,
l'enthalpie avec eau liquide $$k_w=k-q_t[L_v-(C_{pv}-C_l)T]$$
est aussi conservée dans les transformations adiabatiques isobares :
$$
k_w = [(1-q_t) C_{pd} + q_t C_{pv}] T - L_v q_l - L_s q_i
$$
On peut faire le même genre de chose pour obtenir une expression sans $$q_i$$,
ce qui donne : et comme $$L_f-(C_l-C_i)T = \rm{cste}$$, l'enthalpie (sans nom)
$$k_i=k-q_t[L_f-(C_l-C_i)T]$$ est aussi conservée dans les transformations
adiabatiques isobares :
$$
k_i = [(1-q_t) C_{pd} + q_t C_i] T+ L_s q_v + L_f q_l
$$

=== II - Energies statiques: ===
Grandeurs conservées dans les
transformations adiabatiques où le changement
de pression est purement hydrostatique.
\begin{equation}
\begin{array}{ll}
\left .
\begin{array}{lll}
\rm{Adiabatisme :} & dk &\displaystyle{=\frac{dp}{\rho} }\\
& & \\
\rm{Hydrostasie :} & \displaystyle{\frac{dp}{\rho} }&= -g dz
\end{array}
\right \rbrace & \left .
\begin{array}{l}
dk + g dz = 0 \\
\end{array} \right .
\end{array}
\end{equation}

{| border="1" style="border-collapse:collapse"
| Energie statique humide :
| $$h=[(1-q_t)C_{pd}+q_tC_l]T + L_v q_v - L_f q_i + gz$$
|-
| Energie statique avec eau liquide :
| $$h_w=[(1-q_t)C_{pd}+q_tC_{pv}]T - L_vq_l - L_s q_i + gz$$
|-
| Energie statique sans nom :
| $$h_i=[(1-q_t) C_{pd} + q_t C_i] T+ L_s q_v + L_f q_l + gz$$
|-
| Energie statique sèche :
| $$h_d=[(1-q)C_{pd}+qC_{pv}]T +gz$$ 
(= $$h_w$$ avec $$q_l\;=\;0$$ et $$q_i\;=\;0$$)
|}

[[Category:Physique]]
[[Category:Notes Techniques]]