Solving the advection-diffusion equation with the Eulerian-Lagrangian localized adjoint method on unstructured meshes and non uniform time stepping

Authors:
A. Younes;P. Ackerer
Affiliations:
Laboratoire de Génie Industriel, Université de la Réunion, 15 Avenue René Cassin, BP 7151 - 97715, St. Denis Cedex 09, La Réunion, France and Institut de Mécanique de ...;Institut de Mécanique des Fluides et des Solides, Université Louis Pasteur, CNRS - UMR 7507, 2 Rue Boussingault, 67000 Strasbourg, France
Venue:
Journal of Computational Physics
Year:
2005

Citing 4
Cited 4

A characteristics-mixed finite element method for advection-dominated transport problems

SIAM Journal on Numerical Analysis
A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements

Journal of Computational Physics
A family of Eulerian-Lagrangian localized adjoint methods for multi-dimensional advection-reaction equations

Journal of Computational Physics
An ELLAM Scheme for Advection-Diffusion Equations in Two Dimensions

SIAM Journal on Scientific Computing

An efficient algorithm for characteristic tracking on two-dimensional triangular meshes

Computing
Random walk particle tracking simulations of non-Fickian transport in heterogeneous media

Journal of Computational Physics
A fourth-order numerical scheme for solving the modified Burgers equation

Computers & Mathematics with Applications
Impact of environmental factors on mosquito dispersal in the prospect of sterile insect technique control

Computers & Mathematics with Applications

Quantified Score

Hi-index	31.45

Visualization

Abstract

Eulerian-Lagrangian localized adjoint method (ELLAM) is used to solve the advection diffusion equation (ADE) which is a very common mathematical model in physics. In this work, ELLAM is extended to triangular meshes. Standard integration schemes, which perform well for rectangular grids, are improved to reduce oscillations with unstructured triangulations. Numerical experiments for grid Peclet numbers ranking from 1 to 100 show the efficiency of the developed scheme. A new algorithm is also developed in order to avoid excessive numerical diffusion when using many time steps with the ELLAM. The basic idea of this approach is to keep the same characteristics for all time steps and to interpolate only the concentration variations due to the dispersion process at the end of each time step. Although ELLAM requires a lot of integration points for unstructured meshes, it remains a competitive method when using a single or many time steps compared to explicit discontinuous Galerkin finite element method.