Thayer Watkins Silicon Valley & Tornado Alley USA |
---|
Vorticity Equation in Meteorology |
The purpose of this material is to show the relationships among the solenoids terms which arise in meteorological analysis. A solenoid term is of the form ∇f×∇g, where f and g are scalar field variables. The remarkable fact is that when f and g are field variables such as temperature, pressure, density, specific volume and entropy the solenoid terms are proportional.
The magnitude of such a term is equal to |∇f||∇g|sin(φ), where φ is the angle between the two gradients. The solenoid vector lies in the line intersection of the two surfaces of constant levels of f and g.
The blue vector perpendicular to the blue-outlined plane is the gradient vector for the f field. The red vector perpendicular to the red-outlined plane is the gradient vector for the g field. The green vector in the line of intersection of the two planes is their vector cross product. If the two constant surfaces coincide then the angle between the gradients is zero and the solenoid term vanishes. This condition prevails when the atmosphere is barotropic; i.e., when the atmospheric density is a function only of the pressure.
For an inertial frame of reference the equations of motion for a parcel of air are, in vector form,:
where v is the velocity vector, ρ the density, p pressure, g the acceleration due to gravity, k the unit vertical vector and f the vector of friction forces.
The pressure gradient term
is especially important.
This term can be put into an interesting form by noting that from the definition of potential temperature θ:
and when the gradient operator ∇ is applied to this equation the result is
From the defintion of entropy s it follows that
Multiplying through by T and noting that cp∇T is the same as ∇h, where h stands for enthalpy, results in
Thus if the pressure gradient term in the equations of motion is replaced with T∇s - ∇h the result is
Since k is the same as ∇z the above equation is equivalent to
The motion-following derivative dv/dt is composed of an instaneous rate of change at a point and an advection term; i.e.,
The advection term v·∇v can be expressed1 as
Thus the equations of motion for the atmosphere can be expressed in vector form as
The curl operator ∇× can be applied to this equation. The curl of any gradient of a scalar field vanishes; i.e., ∇×∇γ=0 for any scalar field γ because of the equality of cross derivatives. Therefore under the curl operation ∇(v2/2 + h + gz) vanishes.
Also, because the curl of a curl vanishes,
The result of applying the curl operator to the left-hand side of the above equation of motion (1) and taking into account the interchangeability of the time and space derivatives is
Equating this to the result of applying the curl operation to the right-hand side of the equation (1) gives
This form of the vorticity equation points out the role of the intersection or non-intersection of the isothermal surface and the isoentropic surface through the term ∇T×∇s.
Note that since ∇s = cp∇T/T - R∇p/p
Since by the ideal gas law
Thus the ∇T×∇s term in the vorticity equation can be replaced by a term involving ∇p×∇ρ. Generally all of these cross product terms, called solenoid terms, are proportional and they all vanish when the atmosphere is barotropic; i.e.; when ∇ρ always has the same direction as ∇p.
Combining the last two derivations gives
The −∇×(v×q) term in the vorticity equation reduces2 to
These two terms are known as the divergence term and the twisting term, respectively.
Not much can be done analytically with the friction term except note that it is directed opposite to the velocity vector. It is of a low order of magnitude than the other terms and is usually neglected in the analysis.
The term vorticity in meterology usually refers to the vertical component of the vorticity vector. The vertical component of a vector may be obtained analytically by take the dot product of the vector with the unit vector in the vertical direction k; i.e., ζ = k·q. Thus the time rate of change of vorticity ζ is given by
The terms on the right-hand side of the equation are known, respectively, as the divergence term, the twisting or tilting term and the solenoid term.
1This follows from the vector identity
with A=B=v; i.e.,
2 The vector identity
when A=v and B=q=∇×v reduces to
because the divergence of the curl of a vector vanishes.
HOME PAGE OF Thayer Watkins |