Schanuel's lemma: Difference between revisions
en>ZéroBot m r2.7.1) (Robot: Adding de:Lemma von Schanuel |
|||
Line 1: | Line 1: | ||
The '''omega equation''' is of great importance in [[meteorology]] and [[atmospheric physics]]. It is a [[partial differential equation]] for the vertical velocity, <math>\omega</math>, which is defined as a [[Lagrangian]] rate of change of pressure with time, that is, <math>\omega = \frac{dp}{dt}</math>. | |||
<br /> | |||
The equation reads: | |||
:{{NumBlk|:|<math> \sigma\nabla^2_H\omega + f^2\frac{\partial^2\omega}{\partial p^2} = f \frac{\partial}{\partial p} \left[ \mathbf{V}_g \cdot \nabla_H (\zeta_g + f) \right] - \nabla^2_H \left( \mathbf{V}_g\cdot\nabla_H \frac{\partial \phi}{\partial p}\right) </math>|{{EquationRef|1}}}} | |||
where <math> f </math> is the [[Coriolis parameter]], <math> \sigma </math> is the static stability, <math> \mathbf{V}_g </math> is the geostrophic velocity vector, <math> \zeta_g </math> is the geostrophic relative vorticity, <math> \phi </math> is the [[geopotential]], <math> \nabla^2_H </math> is the horizontal Laplacian operator and <math> \nabla_H </math> is the horizontal [[del]] operator.<ref>Holton, J.R., 1992, ''An Introduction to Dynamic Meteorology'' Academic Press, 166-175</ref> | |||
==Derivation== | |||
The derivation of the <math>\omega</math> equation is based on the [[vorticity equation]] and the thermodynamic equation. The [[vorticity equation]] for a frictionless atmosphere may be written as: | |||
:{{NumBlk|:|<math> \frac{\partial \xi}{\partial t} + V \cdot \nabla\eta - f \frac{\partial \omega}{\partial p} = \left( \xi \frac{\partial \omega}{\partial p} - \omega \frac{\partial \xi}{\partial p} \right) + k \cdot \nabla\omega \times \frac{\partial V}{\partial p} </math>|{{EquationRef|2}}}} | |||
Here <math>\xi</math> is the relative vorticity, <math>V</math> the horizontal wind velocity vector, whose components in the <math>x</math> and <math>y</math> directions are <math>u</math> and <math>v</math> respectively, <math>\eta</math> the absolute vorticity, <math>f</math> the [[Coriolis frequency|Coriolis parameter]], <math>\omega = \frac{dp}{dt}</math> the individual rate of change of pressure <math>p</math>. <math>k</math> is the unit vertical vector, <math>\nabla</math> is the isobaric Del (grad) operator, <math>\left( \xi \frac{\partial \omega}{\partial p} - \omega \frac{\partial \xi}{\partial p} \right)</math> is the vertical | |||
advection of vorticity and <math>k \cdot \nabla\omega \times \frac{\partial V}{\partial p} </math> represents the transformation of horizontal vorticity into vertical vorticity.<ref>Singh & Rathor, 1974, Reduction of the Complete Omega Equation to the Simplest Form, Pure and Applied Geophysics, 112, 219-223</ref> | |||
The thermodynamic equation may be written as: | |||
:{{NumBlk|:|<math> \frac{\partial}{\partial t} \left( - \frac{\partial Z}{\partial p} \right) + V \cdot \nabla \left( - \frac{\partial Z}{\partial p} \right) - k\omega = \frac{R}{C_p \cdot g} \cdot \frac{q}{p} </math> |{{EquationRef|3}}}} | |||
<br /> | |||
where <math> k \equiv \left( \frac{\partial Z}{\partial p}\right) \frac{\partial}{\partial p} \ln\theta</math>, in which <math>q</math> is the supply of heat per unit-time and mass, <math>C_p</math>the specific heat of dry air, <math>R</math> the gas constant for dry air, <math>\theta</math> is the potential temperature and <math>\phi</math> is geopotential <math>(gZ)</math>. | |||
The <math>\omega</math> equation ({{EquationNote|1}}) is then obtained from equation ({{EquationNote|2}}) and ({{EquationNote|3}}) by substituting values: | |||
:<math>\xi = \frac{g}{f}\nabla^2 Z </math> | |||
and | |||
<br /> | |||
:<math>\hat k \cdot \nabla\omega \times \frac{\partial V}{\partial p} = \frac{\partial \omega}{\partial y}\frac{\partial u}{\partial p} - \frac{\partial \omega}{\partial x}\frac{\partial v}{\partial p}</math> | |||
<br /> | |||
into ({{EquationNote|2}}), which gives: | |||
<br /> | |||
:{{NumBlk|:|<math>\frac{\partial}{\partial t}\left(\frac{g}{f}\nabla^2 Z \right) + V \cdot \nabla\eta - f \frac{\partial \omega}{\partial p} = \left(\xi \frac{\partial \omega}{\partial p } - \omega \frac{\partial \xi}{\partial p} \right) + \left(\frac{\partial \omega}{\partial x}\frac{\partial v}{\partial p}\right)</math>|{{EquationRef|4}}}} | |||
<br /> | |||
<br /> | |||
Differentiating ({{EquationNote|4}}) with respect to <math>p</math> gives: | |||
<br /> | |||
:{{NumBlk|:|<math>\frac{g}{f}\frac{\partial}{\partial t} \nabla^2 \left(\frac{\partial Z}{\partial p} \right) + \frac{\partial}{\partial p} (V \cdot \nabla\eta) - f \frac{\partial^2 \omega}{\partial p^2} - \frac{\partial f}{\partial p}\frac{\partial \omega}{\partial p} = \frac{\partial}{\partial p}\left(\xi \frac{\partial \omega}{\partial p } - \omega \frac{\partial \xi}{\partial p} \right) + \frac{\partial}{\partial p} \left(\frac{\partial \omega}{\partial y} \cdot \frac{\partial u}{\partial p} - \frac{\partial \omega}{\partial x}\cdot \frac{\partial v}{\partial p}\right)</math>|{{EquationRef|5}}}} | |||
<br /> | |||
<br /> | |||
Taking the Laplacian (<math> \nabla^2 </math>) of ({{EquationNote|3}}) gives: | |||
<br /> | |||
:{{NumBlk|:|<math>\nabla^2 \left(-\frac{\partial Z}{\partial p} \right) + \nabla^2 V \cdot \nabla \left(-\frac{\partial Z}{\partial p} \right) - \nabla^2 k \omega = \frac{R}{C_p \cdot g} \cdot \frac{\nabla^2 q}{p}</math>|{{EquationRef|6}}}} | |||
<br /> | |||
Adding ({{EquationNote|5}}) and ({{EquationNote|6}}), simplifying and substituting <math>gk = \sigma</math>, gives: | |||
:{{NumBlk|:|<math>\nabla^2\omega + \frac{f^2}{\sigma} \frac{\partial^2\omega}{\partial p^2} = \frac{1}{\sigma} \left[ \frac{\partial}{\partial p} J(\phi,\eta) + \frac{1}{f}\nabla^2 J \left(\phi, -\frac{\partial \phi}{\partial p} \right) \right] - \frac{f}{\sigma} \frac{\partial}{\partial p} \left( \frac{\partial \omega}{\partial y} \cdot \frac{\partial u}{\partial p} - \frac{\partial \omega}{\partial x} \cdot \frac{\partial v}{\partial p} \right) - \frac{f}{\sigma} \frac{\partial}{\partial p} \left( \xi \frac{\partial \omega}{\partial p} - \omega \frac{\partial \xi}{\partial p} \right) \frac{R \cdot \nabla^2 q}{C_p \cdot S \cdot p}</math> |{{EquationRef|7}}}} | |||
<br /> | |||
<br /> | |||
Equation ({{EquationNote|7}}) is now a linear differential equation in <math>\omega</math>, such that it can be split into two part, namely <math>\omega_1</math> and <math>\omega_2</math>, such that: | |||
:{{NumBlk|:|<math>\nabla^2\omega_1 + \frac{f^2}{\sigma} \frac{\partial^2\omega_1}{\partial p^2} =\frac{1}{\sigma} \left[ \frac{\partial}{\partial p} J(\phi,\eta) + \frac{1}{f}\nabla^2 J \left(\phi, -\frac{\partial \phi}{\partial p} \right) \right] - \frac{f}{\sigma} \frac{\partial}{\partial p} \left( \frac{\partial \omega}{\partial y} \cdot \frac{\partial u}{\partial p} - \frac{\partial \omega}{\partial x} \cdot \frac{\partial v}{\partial p} \right) - \frac{f}{\sigma} \frac{\partial}{\partial p} \left( \xi \frac{\partial \omega}{\partial p} - \omega \frac{\partial \xi}{\partial p} \right)</math>|{{EquationRef|8}}}} | |||
<br /> | |||
and | |||
<br /> | |||
:{{NumBlk|:|<math>\nabla^2\omega_2 + \frac{f^2}{\sigma} \frac{\partial^2\omega_2}{\partial p^2} =\frac{R \cdot \nabla^2 q}{C_p \cdot \sigma \cdot p}</math>|{{EquationRef|9}}}} | |||
<br /> | |||
<br /> | |||
where <math>\omega_1</math> is the vertical velocity due to the mean baroclinicity in the atmosphere and <math>\omega_2</math> is the vertical velocity due to the non-adiabatic heating, which includes the latent heat of condensation, sensible heat radiation, etc. (Singh & Rathor, 1974). | |||
==Interpretation== | |||
Physically, the omega equation combines the effects of vertical differential of geostrophic absolute vorticity advection (first term on the right-hand side) and three-dimensional Laplacian of thickness thermal advection (second term on the right-hand side) and determines the resulting vertical motion (as expressed by the dependent variable <math>\omega</math>.) | |||
The above equation is used by meteorologists and operational weather forecasters to assess development from synoptic charts. In rather simple terms, positive vorticity advection (or PVA for short) and no thermal advection results in a negative <math>\omega</math>, that is, ascending motion. Similarly, warm advection (or WA for short) also results in a negative <math>\omega</math> corresponding to ascending motion. Negative vorticity advection (NVA) or cold advection (CA) both result in a positive <math>\omega</math> corresponding to descending motion. | |||
==References== | |||
<references/> | |||
==External links== | |||
* [http://amsglossary.allenpress.com/glossary/search?id=omega-equation1 American Meteorological Society definition] | |||
[[Category:Atmospheric dynamics]] | |||
[[Category:Partial differential equations]] |
Revision as of 08:50, 22 January 2014
The omega equation is of great importance in meteorology and atmospheric physics. It is a partial differential equation for the vertical velocity, , which is defined as a Lagrangian rate of change of pressure with time, that is, .
The equation reads:
where is the Coriolis parameter, is the static stability, is the geostrophic velocity vector, is the geostrophic relative vorticity, is the geopotential, is the horizontal Laplacian operator and is the horizontal del operator.[1]
Derivation
The derivation of the equation is based on the vorticity equation and the thermodynamic equation. The vorticity equation for a frictionless atmosphere may be written as:
Here is the relative vorticity, the horizontal wind velocity vector, whose components in the and directions are and respectively, the absolute vorticity, the Coriolis parameter, the individual rate of change of pressure . is the unit vertical vector, is the isobaric Del (grad) operator, is the vertical advection of vorticity and represents the transformation of horizontal vorticity into vertical vorticity.[2]
The thermodynamic equation may be written as:
where , in which is the supply of heat per unit-time and mass, the specific heat of dry air, the gas constant for dry air, is the potential temperature and is geopotential .
The equation (Template:EquationNote) is then obtained from equation (Template:EquationNote) and (Template:EquationNote) by substituting values:
and
into (Template:EquationNote), which gives:
Differentiating (Template:EquationNote) with respect to gives:
Taking the Laplacian () of (Template:EquationNote) gives:
Adding (Template:EquationNote) and (Template:EquationNote), simplifying and substituting , gives:
Equation (Template:EquationNote) is now a linear differential equation in , such that it can be split into two part, namely and , such that:
and
where is the vertical velocity due to the mean baroclinicity in the atmosphere and is the vertical velocity due to the non-adiabatic heating, which includes the latent heat of condensation, sensible heat radiation, etc. (Singh & Rathor, 1974).
Interpretation
Physically, the omega equation combines the effects of vertical differential of geostrophic absolute vorticity advection (first term on the right-hand side) and three-dimensional Laplacian of thickness thermal advection (second term on the right-hand side) and determines the resulting vertical motion (as expressed by the dependent variable .)
The above equation is used by meteorologists and operational weather forecasters to assess development from synoptic charts. In rather simple terms, positive vorticity advection (or PVA for short) and no thermal advection results in a negative , that is, ascending motion. Similarly, warm advection (or WA for short) also results in a negative corresponding to ascending motion. Negative vorticity advection (NVA) or cold advection (CA) both result in a positive corresponding to descending motion.