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.
In this paper we consider two commonly used classes of finite volume weighted essentially non-oscillatory (WENO) schemes in two dimensional Cartesian meshes. We compare them in terms of accuracy, performance for smooth and shocked solutions, and efficiency in CPU timing. For linear systems both schemes are high order accurate, however for nonlinear systems, analysis and numerical simulation results verify that one of them (Class A) is only second order accurate, while the other (Class B) is high order accurate. The WENO scheme in Class A is easier to implement and costs less than that in Class B. Numerical experiments indicate that the resolution for shocked problems is often comparable for schemes in both classes for the same building blocks and meshes, despite of the difference in their formal order of accuracy. The results in this paper may give some guidance in the application of high order finite volume schemes for simulating shocked flows.
f448fe82f3