An ELLAM Scheme for Advection-Diffusion Equations in Two Dimensions

Abstract

We develop an Eulerian-Lagrangian localized adjoint method (ELLAM) to solve two-dimensional advection-diffusion equations with all combinations of inflow and outflow Dirichlet, Neumann, and flux boundary conditions. The ELLAM formalism provides a systematic framework for implementation of general boundary conditions, leading to mass-conservative numerical schemes. The computational advantages of the ELLAM approximation have been demonstrated for a number of one-dimensional transport systems; practical implementations of ELLAM schemes in multiple spatial dimensions that require careful algorithm development are discussed in detail in this paper. Extensive numerical results are presented to compare the ELLAM scheme with many widely used numerical methods and to demonstrate the strength of the ELLAM scheme.

An Efficient Characteristic Method for the Magnetic Induction Equation with Various Resistivity Scales

“…Standard finite difference or finite element methods produce either excessive nonphysical oscillations or extra numerical diffusion, which smears these physical features. The Eulerian-Lagrangian method developed in [5] (see also the references therein) is a characteristic FEM for these problems. This method naturally incorporates all types of boundary conditions in its formulation, is not subject to the severe restrictions imposed by the Courant-Friedrichs-Lewy (CFL) condition, and generates accurate numerical solutions even if large time steps are used.…”

A Mass-Conservative Characteristic Finite Element Scheme for Convection-Diffusion Problems

“…where I n 1 and I n 2 are defined in (21) and (22). Set e n h = φ n h − Π h φ n and η n = φ n − Π h φ n .…”

“…Among them the procedure of the characteristic method is natural from the physical point of view since it approximates particle movements, and it is attractive from the mathematical point of view since it symmetrizes the problem. Many authors have contributed to develop, analyse and apply characteristic finite element schemes; see [1], [4], [6], [7], [8], [13], [14], [15], [16], [21] and references therein.…”

“…In the framework of characteristic methods it it not trivial to maintain this property. Some schemes have been proposed and studied from this point [1], [6], [8], [21].…”

A singularity-avoiding moving least squares scheme for two-dimensional unstructured meshes

“…This scheme is intended to allow the calculation of a function f and its spatial derivatives based upon the values and/or directional derivatives of that function (or a linear combination of these) at various scattered points in the local region. This amount of flexibility is necessary for handling the kinds of boundary conditions that are typically required when solving partial differential equations such as the convection diffusion equation [7]. In order to allow for such general boundary conditions, the following linear relationship is specified for each data point j…”