Eau Condensee Fisrtilp : Différence entre versions

Version actuelle en date du 14 septembre 2021 à 12:20

$$\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}