渦度方程式
Theory and Physics
What is Vorticity?
Professor, the name "vorticity equation" sounds difficult, but first of all, what is "vorticity"?
Vorticity $\boldsymbol{\omega}$ is a vector quantity representing the local rotation of a fluid, defined as the curl of the velocity field.
In two dimensions, it becomes a scalar: $\omega = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$. It quantifies how much a fluid element is rotating.
I see, so the entire fluid doesn't need to be swirling; if there's local rotation, there is vorticity?
Exactly. For example, even in a linear shear flow $u = ky$, $v = 0$, we have $\omega = -k \neq 0$, so there is vorticity. In other words, vorticity is a concept separate from the existence of a "vortex".
Derivation of the Vorticity Equation
So, how do you derive the equation that describes the time evolution of vorticity?
Applying $\nabla \times$ to both sides of the Navier-Stokes equations yields the vorticity equation. For an incompressible fluid, it becomes:
The left side is the material derivative (Lagrangian derivative) of vorticity, representing the rate of change of vorticity seen from a moving fluid particle.
What is the meaning of each of the two terms on the right side?
The first term $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u}$ is the vortex stretching/tilting term. In 3D flow, when a vortex tube is stretched, its cross-section becomes thinner, and due to angular momentum conservation, the vorticity increases. The effect where a tornado or a bathtub vortex spins faster as it becomes thinner is due to this.
The second term $\nu \nabla^2 \boldsymbol{\omega}$ is the viscous diffusion term, representing the effect of vorticity diffusing due to molecular viscosity. The larger the kinematic viscosity coefficient $\nu$, the faster vorticity diffuses and dissipates.
In 2D, the vortex stretching term disappears, right?
Good observation. In 2D, $\boldsymbol{\omega}$ is only in the $z$ direction, so $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u} = 0$, and the vorticity equation becomes a simple advection-diffusion equation.
This means vorticity generation can only occur at boundary surfaces.
Vorticity Generation Mechanism
Doesn't vorticity just appear spontaneously?
In incompressible, barotropic fluids, vorticity is not generated within the fluid interior. The main sources of vorticity are as follows.
- Solid walls: The no-slip condition creates velocity gradients at the wall, generating vorticity. The wall vorticity flux shown by Lighthill (1963) is expressed as $\nu \frac{\partial \omega}{\partial n}\big|_{wall}$.
- Baroclinic effect: When density gradients and pressure gradients are not parallel, vorticity is generated from the term $\frac{1}{\rho^2}(\nabla \rho \times \nabla p)$. Ocean thermohaline circulation is a prime example.
- Non-uniform distribution of body forces: Such as the Coriolis force in a rotating coordinate system.
In CFD, how can we accurately capture vorticity at walls?
Mesh resolution near the wall is critically important. The choice is either to make the $y^+$ of the first cell near the wall less than 1, or to use wall functions to approximate the vorticity flux from the wall shear stress.
Relationship with Kelvin's Circulation Theorem
How are the vorticity equation and Kelvin's circulation theorem connected?
For an inviscid, barotropic fluid with only potential body forces acting, the circulation $\Gamma = \oint \mathbf{u} \cdot d\mathbf{l}$ along a closed curve is conserved. By Stokes' theorem, $\Gamma = \int_S \boldsymbol{\omega} \cdot d\mathbf{S}$, so this corresponds to the vorticity equation when the viscous and baroclinic terms are zero.
In actual CFD, viscosity exists, so circulation is not conserved. Is the vorticity equation like a generalization of the circulation theorem?
That understanding is correct. The vorticity equation is a more general description that includes all viscous and non-barotropic effects.
History of the Vorticity Equation—Helmholtz's Vortex Theorems (1858) and Vortex Tube Conservation
The concept of vorticity and the "Helmholtz vortex theorems" governing the behavior of vortex tubes were published by Hermann von Helmholtz in 1858. The essence of these theorems: ① A vortex tube moves as a line of fluid particles (vortex strength is conserved as it advects), ② A vortex tube cannot end in the fluid—it must form a closed loop or terminate at a fluid boundary. This also explains why a tornado cannot detach from the ground and why vortices are "sucked up" from the sea surface. Lord Kelvin developed this further to derive "Kelvin's circulation theorem (for a barotropic perfect fluid, the circulation Γ around a material curve is constant)". This theorem is the fundamental principle behind lift generation mechanisms (bound vortex) for propellers and wings, and forms the mathematical foundation for modern CFD codes like the vortex lattice method (VLM) and vortex panel methods.
Physical Meaning of Each Term
- Temporal term $\partial(\rho\phi)/\partial t$: Imagine the moment you turn on a faucet. At first, water comes out spluttering and unstable, but after a while, the flow becomes steady, right? This "period of change" is described by the temporal term. The pulsation of blood flow from a heartbeat, or the flow fluctuation each time an engine valve opens and closes—all are unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—meaning setting this term to zero. Since computational cost drops significantly, solving first in steady-state is a basic CFD strategy.
- Convection term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far corner of a room is also because the "carrier," air, transports heat via convection. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as the flow becomes faster, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar" → They are completely different! Convection is transport by flow, conduction is transmission by molecules. There's an order of magnitude difference in efficiency.
- Diffusion term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while it naturally mixes, right? That's molecular diffusion. Now a question—honey and water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. When viscosity is large, the diffusion term becomes strong, and the fluid moves in a "sluggish" manner. In low Reynolds number flows (slow, viscous), diffusion is dominant. Conversely, in high Re number flows, convection overwhelms and diffusion plays a minor role.
- Pressure term $-\nabla p$: When you push the plunger of a syringe, liquid shoots out forcefully from the needle tip, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are tightly packed? That's right, strong winds blow. "Where there is a pressure difference, flow is generated"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: In CFD, "pressure" is often gauge pressure, not absolute pressure. When results go wrong immediately after switching to compressible analysis, mixing up absolute/gauge pressure might be the cause.
- Source term $S_\phi$: Heated air rises—why? Because it becomes lighter (less dense) than its surroundings, so buoyancy pushes it upward. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat from a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... These are all actions that "inject energy or force into the fluid from the outside," expressed by the source term. What happens if you forget the source term? In a natural convection analysis, if you forget to include buoyancy, the fluid doesn't move at all—you get a physically impossible result like turning on a heater in a winter room but the warm air doesn't rise.
Assumptions and Applicability Limits
- Continuum assumption: Valid for Knudsen number Kn < 0.01 (mean free path ≪ characteristic length)
- Newtonian fluid assumption: Shear stress and strain rate have a linear relationship (non-Newtonian fluids require viscosity models)
- Incompressibility assumption (for Ma < 0.3): Treat density as constant. For Mach numbers above 0.3, compressibility effects must be considered.
- Boussinesq approximation (natural convection): Density variation is considered only in the buoyancy term; constant density is used in other terms.
- Non-applicable cases: Rarefied gases (Kn > 0.1), supersonic/hypersonic flows (shock capturing required), free surface flows (VOF/Level Set, etc. required)
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units. |
| Pressure $p$ | Pa | Distinguish between gauge and absolute pressure. Use absolute pressure for compressible analysis. |
| Density $\rho$ | kg/m³ | Air: approx. 1.225 kg/m³ @20°C, Water: approx. 998 kg/m³ @20°C |
| Viscosity coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity coefficient $\nu = \mu/\rho$ [m²/s] |
| Reynolds number $Re$ | Dimensionless | $Re = \rho u L / \mu$. Indicator for laminar/turbulent transition. |
| CFL number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time step stability. |
Numerical Methods and Implementation
Vorticity-Stream Function Method
When solving the vorticity equation in CFD, what is the most basic method?
For 2D incompressible flow, the vorticity-stream function method is a classical approach. It couples the vorticity equation with the Poisson equation for the stream function.
Velocity is recovered as $u = \frac{\partial \psi}{\partial y}$, $v = -\frac{\partial \psi}{\partial x}$. The benefit of this method is that pressure disappears from the unknowns.
What advantages does it have compared to solving the N-S equations directly?
In 2D, there are only two unknowns, $\omega$ and $\psi$, avoiding the pressure-velocity coupling problem (no need for iterations like the SIMPLE method). However, extension to 3D is complex because vorticity becomes a 3-component vector, so it's not very practical.
Spatial Discretization
How do you discretize in space?
Let's look at a typical discretization using the finite difference method (FDM). Using central differences with uniform grid spacing $h$, the Poisson equation becomes:
For the advection term, central differences can cause wiggles when the grid Reynolds number $Re_h = |u|h/\nu > 2$. To avoid this, upwind differencing or the QUICK method is often used.
How does it change with the finite volume method (FVM)?
In FVM, the vorticity flux at cell faces is evaluated. For advective flux, second-order upwind schemes (Second Order Upwind) or TVD schemes are used; for diffusive flux, central differences are used. OpenFOAM's scalarTransportFoam solver is a good reference for FVM implementation solving vorticity advection-diffusion.
Time Integration
How do you advance in time?
Let's summarize typical time integration schemes.
| Scheme | Accuracy | Stability | Features |
|---|---|---|---|
| Forward Euler | 1st order | Conditionally stable | Constraint: $\Delta t < h^2/(4\nu)$ |
| Crank-Nicolson | 2nd order | Unconditionally stable |