next up previous contents
Next: Boundary conditions Up: Non-oscillatory central schemes Previous: The third-order central scheme

Central schemes in two space dimensions

Following the one dimensional setup, one can derive a non-oscillatory, two-dimensional central scheme. Here we sketch the construction of the second-order two-dimensional scheme following [16] (see also [2],[1]). For the two-dimensional third-order accurate scheme, we refer to [24].

We consider the two-dimensional hyperbolic system of conservation laws
To approximate a solution to (1.4.34), we start with a two-dimensional linear reconstruction
Here, the discrete slopes in the x and in the y direction approximate the corresponding derivatives, tex2html_wrap_inline5589, and tex2html_wrap_inline5591 is the characteristic function of the cell tex2html_wrap_inline5593. Of course, it is essential to reconstruct the discrete slopes, tex2html_wrap_inline5595 and tex2html_wrap_inline5597, with built in limiters, which guarantee the non-oscillatory character of the reconstruction; the family of min-mod limiters is a prototype example tex2html_wrap_inline5351tex2html_wrap_inline5351


An exact evolution of this reconstruction, which is based on integration of the conservation law over the staggered volume yields

The exact averages at tex2html_wrap_inline5509 - consult the floor plan in Figure 1.4.6 yields


So far everything is exact. We now turn to approximate the four fluxes on the right of (1.4.37), starting with the one along the East face, consult figure 1.4.7, tex2html_wrap_inline5609. We use the midpoint quadrature rule for second-order approximation of the temporal integral, tex2html_wrap_inline5611; and, for reasons to be clarified below, we use the second-order rectangular quadrature rule for the spatial integration across the y-axis, yielding
In a similar manner we approximate the remaining fluxes.


These approximate fluxes make use of the midpoint values,
tex2html_wrap_inline5615, and it is here that we take advantage of utilizing these midvalues for the spatial integration by the rectangular rule. Namely, since these midvalues are secured at the smooth center of their cells, tex2html_wrap_inline5617, bounded away from the jump discontinuities along the edges, we may use Taylor expansion, tex2html_wrap_inline5619. Finally, we use the conservation law (1.4.34) to express the time derivative, tex2html_wrap_inline5621, in terms of the spatial derivatives, tex2html_wrap_inline5623 and tex2html_wrap_inline5625,
Here, tex2html_wrap_inline5627 and tex2html_wrap_inline5629, are one-dimensional discrete slopes in the x- and y-directions, of the type reconstructed in (1.4.36)-(1.4.36); for example, multiplication by the corresponding Jacobians A and B yields
Equipped with the midvalues (1.4.40), we can now evaluate the approximate fluxes, e,g., (1.4.39). Inserting these values, together with the staggered average computed in (1.4.39), into (1.4.37), we conclude with new staggered averages at tex2html_wrap_inline5639, given by

In summary, we end up with a simple two-step predictor-corrector scheme which could be conveniently expressed in terms on the one-dimensional staggered averaging notations

Our scheme consists of a predictor step
followed by the corrector step

In figures 1.4.8 taken from [16], we present the two-dimensional computation of a double-Mach reflection problem; in figure 1.4.9 we quote from [45] the two-dimensional computation of MHD solution of Kelvin-Helmholtz instability due to shear flow. The computations are based on our second-order central scheme. It is remarkable that such a simple 'two-lines' algorithm, with no characteristic decompositions and no dimensional splitting, approximates the rather complicated double Mach reflection problem with such high resolution. Couple of remarks are in order.

Figure 1.4.8: Double Mach reflection problem computed with the central scheme using tex2html_wrap_inline5641 limiter with CFL=0.475 at t=0.2 (a) density computed with tex2html_wrap_inline5649 cells (b) density computed with tex2html_wrap_inline5651 cells (c) x-velocity computed with tex2html_wrap_inline5651 cells


We conclude this section with brief remarks on further results related to central schemes. Remarks.

  1. Simplicity.

    Again, we would like to highlight the simplicity of the central schemes, which is particularly evident in the multidimensional setup: no characteristic information is required - in fact, even the exact Jacobians of the fluxes are not required; also, since no (approximate) Riemann solvers are involved, the central schemes require no dimensional splitting; as an example we refer to the approximation of the incompressible equations by central schemes, §1.5; the results in [7] provide another example of a weakly hyperbolic multidimensional system which could be efficiently solved in term of central schemes, by avoiding dimensional splitting.

  2. Non-staggering. We refer to [15] for a non-staggered version of the central schemes.
  3. Stability.

    The following maximum principle holds for the nonoscillatory scalar central schemes:


  4. Third-order accuracy. Extensions to third-order accuracy in two space dimensions can be found in [24].

next up previous contents
Next: Boundary conditions Up: Non-oscillatory central schemes Previous: The third-order central scheme

Eitan Tadmor
Mon Dec 8 17:34:34 PST 1997