Constraint-Preserving Upwind Methods for Multidimensional Advection Equations

Authors:
M. Torrilhon;M. Fey
Affiliations:
-;-
Venue:
SIAM Journal on Numerical Analysis
Year:
2004

Citing 0
Cited 3

Increasing the accuracy in locally divergence-preserving finite volume schemes for MHD

Journal of Computational Physics
Splitting based finite volume schemes for ideal MHD equations

Journal of Computational Physics
Constraint Preserving Schemes Using Potential-Based Fluxes. II. Genuinely Multidimensional Systems of Conservation Laws

SIAM Journal on Numerical Analysis

Quantified Score

Hi-index	0.01

Visualization

Abstract

A general framework for constructing constraint-preserving numerical methods is presented and applied to a multidimensional divergence-constrained advection equation. This equation is part of a set of hyperbolic equations that evolve a vector field while locally preserving either its divergence or its curl. We discuss the properties of these equations and their relation to ordinary advection. Due to the constraint, such equations form model equations for general evolution equations with intrinsic constraints which appear frequently in physics.The general framework allows the construction of numerical methods that preserve \emph{exactly} the discretized constraint by special flux distribution. Assuming a rectangular, two-dimensional grid as a first approach, application of this framework leads to a locally constraint-preserving multidimensional upwind method. We prove consistency and stability of the new method and present several numerical experiments. Finally, extensions of the method to the three-dimensional case are described.