Estimation des termes source et puits de l'ozone dans un ... - Inra
Estimation des termes source et puits de l'ozone dans un ... - Inra
Estimation des termes source et puits de l'ozone dans un ... - Inra
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
où L (J kg -1 ) est la chaleur latente <strong>de</strong> vaporisation. K est le coefficient <strong>de</strong> transfert turbulent <strong>de</strong><br />
masse (en général égal à celui pour la chaleur sensible) ; q désigne l’humidité spécifique <strong>de</strong><br />
l’air (kg d’eau kg -1 d’air).<br />
Le flux conductif <strong>dans</strong> le sol (F s ) s’écrit sous la forme:<br />
F s = - λ ∂T/∂z<br />
λ étant la conductivité thermique du milieu (W m -1 K -1 ), z la profon<strong>de</strong>ur du sol, T la<br />
température du sol.<br />
L’évolution <strong><strong>de</strong>s</strong> profils <strong>de</strong> température <strong>dans</strong> le sol est pilotée par l’équation <strong>de</strong> la chaleur:<br />
∂F s /∂z = - Cp s ∂T/∂t<br />
ou Cp s est la chaleur spécifique du sol.<br />
On aboutit à la forme approximative:<br />
F S = λ’ (Ts – To)/∆z<br />
avec λ’ la conductivité thermique équivalente <strong>de</strong> la couche d’épaisseur ∆z, Ts la température<br />
<strong>de</strong> surface du sol <strong>et</strong> To la température <strong>de</strong> la couche d’épaisseur ∆z.<br />
2.2 Loi <strong>de</strong> conservation <strong>de</strong> la masse<br />
La loi <strong>de</strong> conservation <strong>de</strong> l’énergie présentée ci-<strong><strong>de</strong>s</strong>sus est similaire à la loi <strong>de</strong><br />
conservation <strong>de</strong> la masse. Celle-ci dit simplement que l’évolution <strong>de</strong> la concentration au cours<br />
du temps <strong>dans</strong> <strong>un</strong> volume élémentaire est égale à la somme <strong>de</strong> la divergence du flux (ici<br />
vertical, étant donné les hypothèses <strong>dans</strong> lesquelles nous nous plaçons), <strong>et</strong> <strong><strong>de</strong>s</strong> <strong>source</strong>s<br />
volumiques (S(z)). On a :<br />
dc(z) /dt = dF(z)/dz + S(z)<br />
La stationnarité <strong>de</strong> c nous donne :<br />
dF(z)/dz= S(z) (2)<br />
<strong>et</strong> pour estimer S(z) à partir <strong>de</strong> C(z) ou C(z) à partir <strong>de</strong> S(z), on utilise la relation <strong>de</strong> diffusion<br />
établie précé<strong>de</strong>mment :<br />
F= -K ∂C(z)/∂z (3)<br />
On a donc besoin <strong>de</strong> connaître K(z) <strong>dans</strong> tout le domaine que l’on souhaite étudier, c’est à<br />
dire le couvert végétal <strong>et</strong> la couche immédiatement au-<strong><strong>de</strong>s</strong>sus.<br />
8