Abstract
We develop a family of Eulerian-Lagrangian localized adjoint methods for the solution of the initial-boundary value problems for first-order advection-reaction equations on general multi-dimensional domains. Different tracking algorithms, including the Euler and Runge-Kutta algorithms, are used. The derived schemes, which are fully mass conservative, naturally incorporate inflow boundary conditions into their formulations and do not need any artificial outflow boundary conditions. Moreover, they have regularly structured, well-conditioned, symmetric, and positive-definite coefficient matrices, which can be efficiently solved by the conjugate gradient method in an optimal order number of iterations without any preconditioning needed. Numerical results are presented to compare the performance of the ELLAM schemes with many well studied and widely used methods, including the upwind finite difference method, the Galerkin and the Petrov-Galerkin finite element methods with backward-Euler or Crank-Nicolson temporal discretization, the streamline diffusion finite element methods, the monotonic upstream-centered scheme for conservation laws (MUSCL), and the Minmod scheme.
Original language | English |
---|---|
Pages (from-to) | 120-163 |
Number of pages | 44 |
Journal | Journal of Computational Physics |
Volume | 152 |
Issue number | 1 |
DOIs | |
Publication status | Published - Jun 10 1999 |
Keywords
- Advection-reaction equations
- Characteristic methods
- Comparison of numerical methods
- Eulerian-Lagrangian methods
- Linear hyperbolic problems
- Numerical solutions of advection-reaction equations
ASJC Scopus subject areas
- Numerical Analysis
- Modelling and Simulation
- Physics and Astronomy (miscellaneous)
- Physics and Astronomy(all)
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics