Conservation form

From HandWiki

Conservation form or Eulerian form refers to an arrangement of an equation or system of equations, usually representing a hyperbolic system, that emphasizes that a property represented is conserved, i.e. a type of continuity equation. The term is usually used in the context of continuum mechanics.

General form

Equations in conservation form take the form [math]\displaystyle{ \frac{d \xi}{d t} + \boldsymbol \nabla \cdot \mathbf f(\xi) = 0 }[/math] for any conserved quantity [math]\displaystyle{ \xi }[/math], with a suitable function [math]\displaystyle{ \mathbf f }[/math]. An equation of this form can be transformed into an integral equation [math]\displaystyle{ \frac d{d t} \int_V \xi ~ dV = -\oint_{\partial V} \mathbf f(\xi) \cdot \boldsymbol \nu ~ dS }[/math] using the divergence theorem. The integral equation states that the change rate of the integral of the quantity [math]\displaystyle{ \xi }[/math] over an arbitrary control volume [math]\displaystyle{ V }[/math] is given by the flux [math]\displaystyle{ \mathbf f(\xi) }[/math] through the boundary of the control volume, with [math]\displaystyle{ \boldsymbol \nu }[/math] being the outer surface normal through the boundary. [math]\displaystyle{ \xi }[/math] is neither produced nor consumed inside of [math]\displaystyle{ V }[/math] and is hence conserved. A typical choice for [math]\displaystyle{ \mathbf f }[/math] is [math]\displaystyle{ \mathbf f(\xi) = \xi \mathbf u }[/math], with velocity [math]\displaystyle{ \mathbf u }[/math], meaning that the quantity [math]\displaystyle{ \xi }[/math] flows with a given velocity field.

The integral form of such equations is usually the physically more natural formulation, and the differential equation arises from differentiation. Since the integral equation can also have non-differentiable solutions, the equality of both formulations can break down in some cases, leading to weak solutions and severe numerical difficulties in simulations of such equations.

Example

An example of a set of equations written in conservation form are the Euler equations of fluid flow: [math]\displaystyle{ \frac{\partial\rho}{\partial t} + \nabla\cdot(\rho\mathbf u) = 0 }[/math] [math]\displaystyle{ \frac{\partial\rho \mathbf u}{\partial t} + \nabla\cdot(\rho \mathbf u \otimes \mathbf u + p \mathbf I) = 0 }[/math] [math]\displaystyle{ \frac{\partial E}{\partial t} + \nabla\cdot(\mathbf u(E+pV)) = 0 }[/math]

Each of these represents the conservation of mass, momentum and energy, respectively.

See also

Further reading

  • Toro, E.F. (1999). Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer-Verlag. ISBN 3-540-65966-8. 
  • Randall J. LeVeque: Finite Volume Methods for Hyperbolic Problems. Cambridge University Press, Cambridge 2002, ISBN:0-521-00924-3 (Cambridge Texts in Applied Mathematics).