Earth:Semi-Lagrangian scheme

From HandWiki

The Semi-Lagrangian scheme (SLS) is a numerical method that is widely used in numerical weather prediction models for the integration of the equations governing atmospheric motion. A Lagrangian description of a system (such as the atmosphere) focuses on following individual air parcels along their trajectories as opposed to the Eulerian description, which considers the rate of change of system variables fixed at a particular point in space. A semi-Lagrangian scheme uses Eulerian framework but the discrete equations come from the Lagrangian perspective.

Some background

The Lagrangian rate of change of a quantity [math]\displaystyle{ F }[/math] is given by

[math]\displaystyle{ \frac{DF}{Dt} = \frac{\partial F}{\partial t} + (\mathbf{v}\cdot\vec\nabla)F, }[/math]

where [math]\displaystyle{ F }[/math] can be a scalar or vector field and [math]\displaystyle{ \mathbf{v} }[/math] is the velocity field. The first term on the right-hand side of the above equation is the local or Eulerian rate of change of [math]\displaystyle{ F }[/math] and the second term is often called the advection term. Note that the Lagrangian rate of change is also known as the material derivative.

It can be shown that the equations governing atmospheric motion can be written in the Lagrangian form

[math]\displaystyle{ \frac{D\mathbf{V}}{Dt} = \mathbf{S}(\mathbf{V}), }[/math]

where the components of the vector [math]\displaystyle{ \mathbf{V} }[/math] are the (dependent) variables describing a parcel of air (such as velocity, pressure, temperature etc.) and the function [math]\displaystyle{ \mathbf{S}(\mathbf{V}) }[/math] represents source and/or sink terms.

In a Lagrangian scheme, individual air parcels are traced but there are clearly certain drawbacks: the number of parcels can be very large indeed and it may often happen for a large number of parcels to cluster together, leaving relatively large regions of space completely empty. Such voids can cause computational problems, e.g. when calculating spatial derivatives of various quantities. There are ways round this, such as the technique known as Smoothed Particle Hydrodynamics, where a dependent variable is expressed in non-local form, i.e. as an integral of itself times a kernel function.

Semi-Lagrangian schemes avoid the problem of having regions of space essentially free of parcels.

The Semi-Lagrangian scheme

Semi-Lagrangian schemes use a regular (Eulerian) grid, just like finite difference methods. The idea is this: at every time step the point where a parcel originated from is calculated. An interpolation scheme is then utilized to estimate the value of the dependent variable at the grid points surrounding the point where the particle originated from. The references listed contain more details on how the Semi-Lagrangian scheme is applied.

See also

External links

  • ctraj: C++ trajectory library, including semi-Lagrangian tracer codes.

References