A family of Eulerian-Lagrangian localized adjoint methods for multi-dimensional advection-reaction equations
- Univ. of South Carolina, Columbia, SC (United States). Dept. of Mathematics
- Texas A and M Univ., College Station, TX (United States). Inst for Scientific Computation
- Mobil Technology Co., Dallas, TX (United States). Upstream Strategic Research Center
- Sultan Qaboos Univ., Muscat (Oman). Dept. of Mathematics and Statistics
Many difficult problems arise in the numerical simulation of fluid flow processes within porous media in petroleum reservoir simulation and in subsurface contaminant transport and remediation. The authors 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 full 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.
- OSTI ID:
- 361771
- Journal Information:
- Journal of Computational Physics, Vol. 152, Issue 1; Other Information: PBD: 10 Jun 1999
- Country of Publication:
- United States
- Language:
- English
Similar Records
Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. II. Implicit Methods
Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems