6.2. Математическая постановка модели атмосферного пограничного слоя

Математическая модель включает следующие уравнения:

Уравнение неразрывности


\[\frac{\partial (\rho u)}{\partial x}+\frac{\partial (\rho v)}{\partial y}+\frac{\partial (\rho w)}{\partial z}=0;\] (6.1)

Уравнения движения

\[\rho \left( \frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}+w\frac{\partial u}{\partial z} \right)=-\frac{\partial p}{\partial x}+\rho fv+{{K}_{H}}\left( \frac{{{\partial }^{2}}u}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}u}{\partial {{y}^{2}}} \right)+\frac{\partial }{\partial z}\left( K_{Z}^{m}\frac{\partial u}{\partial z} \right);\] (6.2)
\[\rho \left( \frac{\partial v}{\partial t}+u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}+w\frac{\partial v}{\partial z} \right)=-\frac{\partial p}{\partial y}-\rho fu+{{K}_{H}}\left( \frac{{{\partial }^{2}}v}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}v}{\partial {{y}^{2}}} \right)+\frac{\partial }{\partial z}\left( K_{Z}^{m}\frac{\partial v}{\partial z} \right);\] (6.3)
\[\rho \left( \frac{\partial v}{\partial t}+u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}+w\frac{\partial v}{\partial z} \right)=-\frac{\partial p}{\partial y}-\rho fu+{{K}_{H}}\left( \frac{{{\partial }^{2}}v}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}v}{\partial {{y}^{2}}} \right)+\frac{\partial }{\partial z}\left( K_{Z}^{m}\frac{\partial v}{\partial z} \right);\] (6.4)

Здесь \(t\) - время, \(u\),\(v\),\(w\) – продольная, поперечная и вертикальная компоненты вектора осредненной скорости ветра в направлении декартовых координат (\x\),(\y\),(\z\) соответственно, \(\rho \) – плотность, \(f\) – параметр Кориолиса, \({{K}_{H}}\) – коэффициент горизонтальной диффузии, \(K_{z}^{m}\) – коэффициент вертикальной диффузии количества движения.

Уравнение баланса энергии

\[\rho \left( \frac{\partial \theta }{\partial t}+u\frac{\partial \theta }{\partial x}+v\frac{\partial \theta }{\partial y}+w\frac{\partial \theta }{\partial z} \right)={{K}_{H}}\left( \frac{{{\partial }^{2}}\theta }{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}\theta }{\partial {{y}^{2}}} \right)+\frac{\partial }{\partial z}\left( K_{Z}^{h}\frac{\partial \theta }{\partial z} \right)+\frac{\theta }{{{c}_{p}}T}\left( {{Q}_{rad}}-\rho {{L}_{w}}{{\Phi }_{v}} \right);\]  (6.5)

Здесь \(T\) – абсолютная температура, \(\theta \) – потенциальная температура, \({{c}_{p}}\) – теплоемкость воздуха при постоянном давлении, \({{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{Q}}_{rad}}\) – нагрев (охлаждение) атмосферы за счет лучистых потоков тепла, распространяющихся в атмосфере,  \(\rho {{L}_{w}}{{\Phi }_{v}}\) – изменение температуры за счет фазовых переходов воды в атмосфере,  \(K_{z}^{h}\) – коэффициент вертикальной диффузии тепла и влаги, \({{L}_{w}}=2,501\cdot {{10}^{6}}Дж/кг\) – теплота парообразования.