A class of semi-Lagrangian approximations for fluids

Abstract

Values of a fluid variable ψ at any two points of a space-time continuum are related through the integral of the space-time gradient of ψ along an arbitrary contour connecting these two points (Stokes' theorem). Noting that spatial and temporal components of the gradient are related through the fluid equations, and selecting the contour composed of a parcel trajectory and an appropriate residual, leads to the integral form of the fluid equations, which is particularly convenient for finite-difference approximations. In these equations, the inertial and forcing terms are separated such that forces are integrated along a parcel trajectory (the Lagrangian aspect), whereas advection of the variable is evaluated along the residual contour (the Eulerian aspect). The virtue of this method is an extreme simplicity of the resulting solver; the entire model for a fluid may be essentially built upon a single one-dimensional Eulerian advection scheme while retaining the formal accuracy of its constant-coefficient limit. The Lagrangian aspect of the approach allows for large-Courant-number (>1) computations in a broad spectrum of dynamic applications. -from Authors