The vorticity equation of fluid dynamics describes evolution of the vorticity ω→ of a particle of a fluid as it moves with its flow, that is, the local rotation of the fluid (in terms of vector calculus this is the curl of the velocity field). The equation is:
where d/dt the total time derivative operator, also denoted by in capital D notation as D/Dt, v→ is the velocity field, ρ is the local fluid density, p is the local pressure, τ is the viscous stress tensor and B→ represents the sum of the external body forces. The first source term on the right hand side represents vortex stretching. The equation is valid in the absence of any concentrated torques and line forces, for a compressible Newtonian fluid.
- The term dω→/dt on the left-hand side is the material derivative of the vorticity vector ω→. It describes the rate of change of vorticity (that is, the angular acceleration) of the fluid particle. This change can be attributed to unsteadiness in the flow (∂ω→/∂t, the unsteady term) or due to the motion of the fluid particle as it moves from one point to another (v→ ∙ (∇→ω→), the convection term).
- The term (ω→∙ ∇→) v→ on the right-hand side describes the stretching or tilting of vorticity due to the velocity gradients. Note that ∇→v→ is an order-2 tensor with nine components.
- The term ω→(∇→ ∙ v→) describes stretching of vorticity due to flow compressibility. It follows from the Navier-Stokes equation for continuity,namely
- where V = 1/ρ is the specific volume of the fluid element. One can think of ∇→ ∙ V→ as a measure of flow compressibility. Sometimes the negative sign is included in the term.
- The term (1/ρ2)∇→ρ × ∇→p is the baroclinic term. It accounts for the changes in the vorticity due to the intersection of density and pressure surfaces.
- The term ∇→ × (∇→ ∙ τ/ρ), accounts for the diffusion of vorticity due to the viscous effects.
- The term ∇→ × B→ provides for changes due to external body forces. These are forces that are spread over a three-dimensional region of the fluid, such as gravity or electromagnetic forces. (As opposed to forces that act only over a surface (like drag on a wall) or a line (like surface tension around a meniscus).
- In case of conservative body forces, ∇→ × B→ = 0.
- For a barotropic fluid, ∇→ρ × ∇→p = 0. This is also true for a constant density fluid (including incompressible fluid) where ∇→ρ = 0. Note that this is not the same as an incompressible flow, for which the barotropic term cannot be neglected.
- For inviscid fluids, the viscosity tensor τ is zero.
Thus for an inviscid, barotropic fluid with conservative body forces, the vorticity equation simplifies to
Alternately, in case of incompressible, inviscid fluid with conservative body forces,
For a brief review of additional cases and simplifications, see also.
Now, vorticity is defined as the curl of the velocity vector. Taking curl of momentum equation yields the desired equation.
The following identities are useful in derivation of the equation,
- , where ϕ is any scalar field.
In specific sciences
In the atmospheric sciences, the vorticity equation can be stated in terms of the absolute vorticity of air with respect to an inertial frame, or of the vorticity with respect to the rotation of the Earth. The absolute version is
- K. P. Burr, Marine Hydrodynamics, Lecture 9, MIT Lectures. http://web.mit.edu/2.20/www/lectures/Lecture-2014/lecture9-2014.pdf
- Utpal Manna and S. S. Sritharan (2007), "Lyapunov Functionals and Local Dissipativity for the Vorticity Equation in Lp and Besov spaces,", Differential and Integral Equations 20 (5): 581–598.
- V. Barbu and S. S. Sritharan, “M-accretive Quantization of the Vorticity Equation", in Semi-groups of Operators: Theory and Applications, edited by A. V. Balakrishnan, Birkhauser, Boston, 2000, pp. 296–303. http://www.nps.edu/Academics/Schools/GSEAS/SRI/BookCH12.pdf
- A. M. Krigel, "Vortex evolution", Geophysical, Astrophysical Fluid Dynamics, 1983, 24, pp. 213–223.
|This article needs additional citations for verification. (May 2009)|