Mixed and hybrid finite element methods
Mixed and hybrid finite element methods
Convection enhanced diffusion for periodic flow
SIAM Journal on Applied Mathematics
A characteristics-mixed finite element method for advection-dominated transport problems
SIAM Journal on Numerical Analysis
SIAM Journal on Numerical Analysis
Journal of Computational Physics
SIAM Journal on Numerical Analysis
Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
SIAM Journal on Numerical Analysis
Computational solution of two-dimensional unsteady PDEs using moving mesh methods
Journal of Computational Physics
Locally Divergence-Free Discontinuous Galerkin Methods for MHD Equations
Journal of Scientific Computing
On Resistive MHD Models with Adaptive Moving Meshes
Journal of Scientific Computing
ELLAM for resolving the kinematics of two-dimensional resistive magnetohydrodynamic flows
Journal of Computational Physics
A Fully Mass and Volume Conserving Implementation of a Characteristic Method for Transport Problems
SIAM Journal on Scientific Computing
SIAM Journal on Scientific Computing
Hi-index | 31.46 |
Tracer transport is governed by a convection-diffusion problem modeling mass conservation of both tracer and ambient fluids. Numerical methods should be fully conservative, enforcing both conservation principles on the discrete level. Locally conservative characteristics methods conserve the mass of tracer, but may not conserve the mass of the ambient fluid. In a recent paper by the authors [T. Arbogast, C. Huang, A fully mass and volume conserving implementation of a characteristic method for transport problems, SIAM J. Sci. Comput. 28 (2006) 2001-2022], a fully conservative characteristic method, the Volume Corrected Characteristics Mixed Method (VCCMM), was introduced for potential flows. Here we extend and apply the method to problems with a solenoidal (i.e., divergence-free) flow field. The modification is a computationally inexpensive simplification of the original VCCMM, requiring a simple adjustment of trace-back regions in an element-by-element traversal of the domain. Our numerical results show that the method works well in practice, is less numerically diffuse than uncorrected characteristic methods, and can use up to at least about eight times the CFL limited time step.