Background

The advection equation is the partial differential equation that governs the motion of a conserved scalar as it is advected by a known velocity field. It is derived using the scalar's conservation law, together with Gauss's theorem, and taking the infinitesimal limit. Perhaps the best image to have in mind is the transport of dissolved salt in water. The advection equation expressed as:

$$\frac{\partial\psi}{\partial t}+\nabla \ (\psi {\b v})=0$$ (2)

where $\nabla$ is the divergence. Frequently, it is assumed that the velocity field is solenoidal, that is, that . If this is so, the above equation reduces to

$$\frac{\partial\psi}{\partial t}+\b{v}\cdot \nabla\psi =0$$ (3)

In particular, if the flow is steady, $\b{v}\cdot \nabla\psi =0$ which shows that $\psi$ is constant along a streamline.

The advection equation is not simple to solve numerically: the system is a hyperbolic partial differential equation, and interest typically centers on discontinuous "shock" solutions (which are notoriously difficult for numerical schemes to handle). Even in one space dimension and constant velocity, the system remains difficult to simulate (it is a standard test for advection schemes known as the pigpen problem ). The equation becomes

$$\frac{\partial q}{\partial t}+{\b v}\frac{\partial q}{\partial t} =0$$ (4)

where $q=q(x,t)$.