摘要
We consider constraint preserving multidimensional evolution equations.A prototypical example is provided by the magnetic induction equation of plasma physics.The constraint of interest is the divergence of the magnetic field.We design finite volume schemes which approximate these equations in a stable manner and preserve a discrete version of the constraint.The schemes are based on reformulating standard edge centered finite volume fluxes in terms of vertex centered potentials.The potential-based approach provides a general framework for faithful discretizations of constraint transport and we apply it to both divergence preserving as well as curl preserving equations.We present benchmark numerical tests which confirm that our potential-based schemes achieve high resolution,while being constraint preserving.
基金
E.Tadmor Research was supported in part by NSF grant 07-07949 and ONR grant N00014-091-0385.