Conservative and Easily Implemented Finite Volume Semi-Lagrangian WENO Methods for 1D and 2D Hyperbolic Conservation Laws

The paper is devised to propose finite volume semi-Lagrange scheme for approximating linear and nonlinear hyperbolic conservation laws. Based on the idea of semi-Lagrangian scheme, we transform the integration of flux in time into the integration in space. Compared with the traditional semi-Lagrange scheme, the scheme devised here tries to directly evaluate the average fluxes along cell edges. It is this difference that makes the scheme in this paper simple to implement and easily extend to nonlinear cases. The procedure of evaluation of the average fluxes only depends on the high-order spatial interpolation. Hence the scheme can be implemented as long as the spatial interpolation is available, and no additional temporal discretization is needed. In this paper, the high-order spatial discretization is chosen to be the classical 5thorder weighted essentially non-oscillatory spatial interpolation. In the end, 1D and 2D numerical results show that this method is rather robust. In addition, to exhibit the numerical resolution and efficiency of the proposed scheme, the numerical solutions of the classical 5th-order WENO scheme combined with the 3rd-order Runge-Kutta temporal discretization (WENOJS) are chosen as the reference. We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS, but with less CPU time.


Introduction
The characteristic of hyperbolic conservation laws is that they can generate discontinuous solutions even if the initial condition is smooth.It is this property that leads to the challenge in developing the high-order numerical methods.The main difficulty of designing high-order methods is how to maintain the highorder accuracy around the smooth regions and meanwhile suppress the spurious oscillations around the regions with large gradients.In the past few decades, many high-order numerical methods were developed to evolve the hyperbolic problem, e.g. the essentially non-oscillatory schemes (ENO) [1] [2] [3] [4], the weighted essentially non-oscillatory schemes (WENO) [5] [6] and the discontinuous Galerkin finite element method [7] [8], and so on.
The discretization procedure of the classical WENO schemes [5] can be divided into two stages.First, through spatial discretization implemented by WENO reconstruction, leaving the equations continuous in time; this leads to a system of ordinary differential equations in time.Then, any numerical methods for ordinary differential equations, e.g.Runge-Kutta methods, can be used to evolve this system.In order to make the solution stable in total variation norm, Shu and Osher [3] and Shu [9] devised a set of TVD Runge-Kutta methods as ODE solvers.Taking three-stage TVD Runge-Kutta methods for example, three WENO reconstruction processes are needed to evolve one time step.And they also suffer from a CFL time step stability restriction.
To fix the problem above, the semi-Lagrangian methods [10] [11] [12] [13] avoid the time discretization by characteristic tracing technique and have more alleviative restriction for CFL number.Recently, in the papers [14] [15] [16] [17], the authors combined the idea of semi-Lagrangian with WENO reconstruction for solving advection problems.The semi-Lagrangian methods above depend on the characteristic tracing.However, for generally nonlinear cases, it is impossible to find the trace points exactly (even finding the trace points with high accuracy is very hard).It is this restriction that makes the semi-Lagrangian methods very hard to apply to nonlinear problems.The authors Huang and Arbogast in [18] devised a semi-Lagrangian method for nonlinear conservation laws.However, a costly flux correction is needed for keeping high accuracy and numerical stability.
In this paper, we try to directly evaluate the average flux ( ) ( ) where n t is the time step and x + is the mesh boundary.The main idea of evaluation of the average fluxes is the transformation of the integration of flux function in time into the integration in space.A similar transformation had appeared in [16].The computation for integration in space only depends on the spatial interpolation polynomial.Hence, the scheme here can be implemented easily as long as the procedure of spatial interpolation is available.Several high-order polynomial interpolations can be chosen, e.g. the PPM [19], ENO [1] [2] or WENO [5] [6], etc.In this paper, we choose the classical 5th-order WENO interpolation to reconstruct the polynomial in space and we call the scheme here as WENOEL.It is noted that the scheme is obviously conservative since we directly evaluate the average flux at each cell edge.As the semi-Lagrangian methods presented in literature, we also need the tracking back points along cha-F.X. Hu racteristic lines to evaluate the average flux.Since the characteristic points can be tracked back exactly (or with high-order accuracy) only for some simple linear equations, we do not try to find the characteristic points with high-order accuracy.For linear problems, the scheme can achieve optimal accuracy because the average flux can be evaluated with high order of accuracy.For nonlinear problems, in testing the order of accuracy for smooth solution, no optimal accuracy can be guaranteed.For testing the resolution in nonlinear problems of the scheme, we choose the classical WENOJS (the 5th-order WENO reconstruction with 3rd-order Runge-Kutta time discretization) as benchmark and compare them for several 1D and 2D examples.We find that the scheme proposed in the paper generates comparable solutions with that of WENOJS but needs less than half computing time.The free CFL restriction is another advantage of semi-Lagrangian methods.For advection problems, the CFL number can be chosen generally larger than 1.However, for nonlinear hyperbolic conservation laws, the CFL 1 > will lead to numerical instability generally if there exist shocks.In this paper, the CFL 0.6 = is chosen for the WENOEL scheme in all the numerical examples.
In the paper that follows, we will present the evaluation procedure of the average flux in Section 2 in detail.In Section 3, to show the high resolution and efficiency of our scheme, we compare the schemes WENOEL proposed here and WENOJS for several 1D and 2D examples.And from the comparisons of resolution and efficiency for these 1D and 2D examples, we can find that the proposed method in this paper can present comparable results with the classical WENOJS scheme, but with less CPU time.

The Finite Volume WENOEL Scheme
Consider the conservation laws subject to initial condition Rearranging this equation and dividing by x ∆ yields Denote the i th cell average by ( ) and average flux at cell edge then (3) can be written as ( ) The value n i U will approximate the average value j i u .Given the approx- imated cell average n i U , the formula (5) tells us that be updated if the time integrals on the right of (4) can be evaluated efficiently.Denoting the approximated flux function then we obtain ( )

The Linear Equation
Here, we mainly consider the 5th-order finite volume version because of its popularity and other cases can be obtained similarly.Consider the linear advection equation In this case, the flux function . For this equation, the initial condition ( ) For evaluating the average flux in (4), we can firstly apply a 5th-order reconstruction based on piecewise constant average fluxes { } 1 , .
Here, we assume the advection velocity 0 c > , then ( ) simply spreads to the right with velocity c , which gives , .
In this case, the average flux Omitting the high-order term ( ) , we obtain the numerical flux ( ) The last equality in (11) is obtained by integration of substitution ( ) − .From the Equation (11), the integration in time [ ] Due to the integrand ( ) i P x is reconstruction polynomial, the last integration in ( 12) can be com- puted exactly.
Remark 1.This transformation is the main idea of this paper.And it is similar to the equation (3.9) in [16] but with a slightly different way.The main difference is that the operation in [16] transforms the flux integration in time at cell center into the integration of ( ) where i x * is the backward characteristic point of cell center i x .And here we transform the flux integration in time at cell edge into the integration of interpolation polynomial of flux function in space.This difference makes the extension of the scheme to nonlinear cases much easier.
In the following of this subsection, we will present the 5th-order WENO reconstruction process to approximate the integral in (12) ( ) The 5th-order WENO reconstruction procedure is represented as the convex combination of three 3rd-order reconstructions.First, we intend to reconstruct the 3rd-order conservative polynomials on cell It is noted that there are 3 three-point stencils contained the cell i I can be used to reconstruct the polynomials.We denote these stencils respectively by , , , , , , , , .
The corresponding polynomial of each stencil is ( ) where the subscript i denotes the polynomial on cell i I and superscript j denotes the reconstruction based on stencil j i S .The coefficients in ( 13) are given by (for simplicity, we omit superscript in where 0,1, 2 j = and So far, we have obtained the conservative interpolation polynomials ( ) In the end, the integral in ( 12) can be expressed as ( ) ( ) where the linear optimal weights j d are .For alleviating the effect of the non-smooth stencils, the nonlinear weights can be constructed as follows The β is the indicator of smoothness based on the stencil j i S ,

(
) ( ) Finally, the numerical flux should be expressed as ( ) Substituting the Formulas ( 13) into (16) gives In contrast, when the advection velocity 0 c < , ( ) , n u x t spreads to the left with velocity c , which gives , , and the flow rate , .
Similarly, the average flux The nonlinear weight 1 j i w + can be constructed similarly as (14) and (15).The linear optimal weights are Substituting the formula ( 13) into (20) gives The scheme proposed in the paper is conservative for linear and nonlinear cases.The property of conservation is obvious since the scheme ( 6) is represented in conservation form.Remark 2. In [16], the authors Qiu and Shu presented a semi-Lagrangian finite difference WENO method for advection in incompressible flow.Actually, there are two reconstruction procedures needed to obtain the numerical flux at cell edge x + .It is noted that, for linear equation with constant coefficient, the formulas of average flux ( 17) and ( 21) are totally same as ones in [16].However, these two methods are devised by totally different procedure.They mainly have the following differences: the method present here is finite volume; our method can be easily extended to nonlinear cases by freezing locally and detailed description is given in Subsection 2.3; there are always exist linear optimal weights to construct the WENO schemes for our scheme which is not satisfied the schemes in [16].

Algorithm
Here, we conclude the algorithm for computing numerical flux 1) Distinguish the propagation direction of solution at 2) The propagation velocity can be chosen to be ( ) 3) Insert the propagation velocity into (17) or (21) to get numerical flux

The Nonlinear Equations
In this subsection, we will extent the method above to nonlinear case.Firstly, it is noted that for nonlinear case the method here cannot approach the optimal accuracy since the solution no longer simply translates uniformly.Instead, it deforms as it evolves and the characteristic curves are no longer parallel straight lines, which lead to the evaluation of average fluxes is hard to obtain.And generally the tracking back points cannot be found exactly (even cannot find the points with high accuracy).Hence for nonlinear case, rather than trying to find the tracking back points, we freeze the equation to linear formation locally and apply the procedure in Subsection 2.1 to it.For solving the nonlinear case, the propagation direction is firstly distinguished by Rankine-Hugoniot jump conditions and propagation velocity is chosen to be Remark 3. The semi-Lagrangian finite volume method proposed here is largely dependent on the reconstruction from the cell flux average ( ) ( ) For finite volume method, only the cell average ( ) , the cell flux average can be apparently chosen to be However, for nonlinear problem, the cell flux average must be computed by some numerical quadrature.In this paper, we choose Gauss quadrature to evaluate the cell flux average where g k A and g k x are the Gauss weights and Gauss points, respectively.For simple, the value of ( ) x can be obtained by the 5thorder Lagrangian interpolation.
In addition, for solving Euler equations, the Roe speed is used to identify the upwind direction which may lead to some numerical instabilities, especially for 2D cases.In Subsection 3.3, in computing the numerical tests: double Mach reflection and Mach 3 wind tunnel with a step, there indeed exist numerical instabilities.In addition, it is noted that this numerical instability not only appears in the WENOEL scheme, but also appears in the WENOJS scheme.It is said that the numerical instability is mainly due to the choice of Roe speed.The detailed introductions for the numerical instabilities of upwind schemes are referred to [20] [21] [22].In order to eliminate such kind of instabilities, the H-correction procedure in [21] is taken to solve these flaws.For these two 2D problems, the numerical flux at the cell edge ( ) min , , then the numerical flux F + can be chosen to be (17) or (21) and the average fluxes k f in these two equations are corres- pondingly substituted by ( ) , 2, , 3 λ is the eigenvalue and the value , and , i j a is the speed of sound.
2) Otherwise, a more dissipative Lax-Friedrich flux splitting method is used to split the flux , k j f into positive and negative fluxes ( ) For the positive and negative fluxes, the propagation velocities are chosen to be  17) and ( 21), respectively, we can get the positive and negative numerical fluxes

Numerical Results
In this section, we will apply this method to 1D and 2D hyperbolic problems.At present, the Strang dimensional splitting method [23] [24] is applied to 2D problems, which is used to split the 2D equation into two 1D equation.In each coordinate direction, the equation is evolved by the proposed 1D WENOEL method.The detailed comparisons of resolution and efficiency are shown between our method and WENOJS method (5th-order spatial and 3rd-order TVD RK temporal discretization).In the end, it is noted that, in the WENOJS method, the upwinding (not flux vector splitting) is used to construct the numerical flux.
So the WENOJS method in this paper can also be called by WENOJS-Roe.The detailed description of WENOJS-Roe is referred to the procedure 2.4 on page 33 in [25].
For the choice of CFL number, It is found that when 0 CFL 1 < < the WENOEL scheme is numerically stable in all the examples we have tested.Especially, we had chose CFL = 0.9 to simulate all the examples appeared in this paper and no unstable phenomenon had been found.For comparing the resolution and computing time between the WENOEL and WENOJS schemes, we choose CFL = 0.6 for these two schemes for all the tests contained the discontinuities.
This choice of CFL number for the WENOJS scheme is also adopted in many classical literatures, e.g.[25].For the problem (22) (23), which is used to test the order of accuracy of schemes, we choose for WENOJS in order that spatial error is dominant but still let CFL = 0.6 for WENOEL.

1D Scalar Examples
In this subsection, we will consider linear advection, Burger's and Buckley-Leverett equations.For the linear advection equation, two initial conditions are used to test the schemes.We use a smooth initial condition to test the order of accuracy of the scheme and a condition contained discontinuities and high-frequency waves to test the performance of the scheme in simulating discontinuous and large-gradient solution.For the Burger equation, we also consider two initial conditions which commonly used to validate the proposed scheme.For the nonconvex Buckley-Leverett equation, we study a simple model for two-phase fluid flow in a porous medium.Except for the comparison of resolution between the WENOJS and WENOEL schemes, the computing time is also shown in this subsection.

Linear Advection Equation
Consider the linear advection equation subject to two initial conditions and periodic boundary conditions.One initial condition is ( ) and the other is where The problem ( 22) (23) here is used to test the accuracy of our method in linear case.The Table 1 and Table 2 present the comparison of errors and accuracy between the WENOEL and WENOJS methods.The output time is chosen to be 2 t = and the time steps are chosen to be 0.6 t x ∆ = ∆ and ( ) for the WENOEL and WENOJS, respectively.As the classical WENOJS method, the WENOEL also arrives 5th-order of accuracy.For the 1 l and l ∞ errors, the WENOEL method shows much better results than the WENOJS method.This  23) with WENOJS method.The time step is chosen to be ( ) and final time is result is expected since the tracking point can be found exactly and the numerical flux is obtained much accurately.The solution of (22) (24) contains a discontinuous square pulse and several continuous but high-gradient profiles.This problem is conventionally employed to test high-resolution schemes.The output time is chosen to be 2 t = and the time step is chosen to be 0.6 t x ∆ = ∆ for the WENOEL and WENOJS schemes.
Figure 1 gives the solutions of the WENOEL and WENOJS methods with cells 200 N = . From Figure 1, we can find that the WENOEL and WENOJS methods can both capture the discontinuous and steep solutions.And for the first and third high-frequency waves, the WENOEL performs slightly better than WENOJS method.Figure 2 is used to test the long-time behavior of these two methods  after 10 cycles, i.e. 20 t = .Obviously, following long-time evolution, the WENOEL gives higher-resolution solution around four waves than that of WENOJS method.In addition, we present the comparison of CPU time (in seconds) and 1 l error for this problem between the WENOEL and WENOJS methods.

Burger Equation
In this subsection, we consider the Burger's equation subject to two initial conditions.one initial condition is ( ) and the other is For the nonlinear scalar equation, we first use Rankine-Hugoniot jump condition to identify the upwinding direction, then compute the average flux by (17) or (21).In addition, an entropy fix is used to modify the average flux when rarefaction wave appears.For the problems (25) ( 26) and ( 25) (27), we both set for the WENOEL and WENOJS schemes, where For the problem (25) (26), the numerical solutions are computed with cells 100 N = until time 1 t = .From Figure 3, the solution of this problem contains a rarefaction and a shock, and these two methods arrive almost same resolution.In Figure 4, the problem (25) (27) has the same set but the output time 0.37 t = , and has the same conclusion with last problem.In a word, for nonlinear Burger's  .The blue "  " and red " ◊ " de- note the solutions of WENOEL and WENOJS methods, respectively.equation, the WENOEL scheme can also obtain comparable numerical results and only need half CPU time.The presentation of CPU time is omitted for saving space.

Buckley-Leverett Equation
As an example where nonconvex flux functions arise, we study a simple model for two-phase fluid flow in a porous medium.Consider the scalar conservation law ( ) with the initial condition ( ) Figure 5 shows the solutions computed by the WENOEL and WENOJS methods with cells 100 N = until time 0.4 t = . Also, we choose the time steps for WENOEL and WENOJS schemes, where From Figure 5 the WENOEL method generates the solution which is comparable with that of the WENOJS method.For the nonconvex flux problem, the WENOEL scheme can also have a robust performance as WENOJS scheme.

The 1D Euler Equation
In this subsection, we consider 1D Euler equations since one of the main application areas of high-resolution scheme is compressible gas dynamics,  ) where ρ , u , p , E are density, velocity, pressure and total energy, respec- tively.The system of equations is closed by the equation of state for an ideal polytropic gas: where the ratio of specific heats 1.4 γ = .
For solving Euler equations, Roe linearization (i.e.Roe average) is used to locally freeze a nonlinear system to a linear system.Then this system is decoupled into three advection equations and each equation is solved by procedure in Section 0. It is noted that the propagation direction of solution of each advection equation is distinguished by corresponding eigenvalue of Roe matrix.
For the problem (30) (31), called Lax problem [26], we solve it with 0.6 t xλ ∆ = ∆ for the WENOEL and WENOJS methods.Figure 6 is computed with cells 200 N = and output time 1.3 t = .Apparently, the WENOEL me- thod performs almost the same as WENOJS method and presents a robust solution for this shock-tube problem.The exact solution is obtained by the solver presented in [24].
The problem (30) (32), called shock entropy wave interaction [4], is usually used to test the high-order methods for capturing the high-frequency waves.The CFL number is also set to be 0.6 t xλ ∆ = ∆ for the WENOEL and WENOJS methods.We evolve the equations up to .It is plainly seen from Figure 7 that, as WENOJS, the WENOEL method performs with the high resolution.Furthermore, from Figure 8, we can find that in the part of high-frequency waves the WENOEL method performs less dissipative than the WENOJS method.
The problem (30) (33) was originally proposed as a benchmark for testing several numerical methods by Woodward and Colella [19].The reflective boundary conditions are applied to both boundaries.The CFL number is also set to be 0.6 t xλ ∆ = ∆ for the WENOEL and WENOJS methods.We evolve the equations up to 0.38 t = with cells 400 N = . The reference solution in Figure 9 is computed by WENOJS method with cells 2000 N = . As the last two problems, the WENOEL method gives almost the same resolution as WENOJS method.

The 2D Euler System
Finally, we consider a numerical experiment for 2D Euler equations for gas dynamics, ( ) ( ) where the equation of state is ( ) ( ) .The blue "  " and red " ◊ " denote the solutions of WENOEL and WENOJS methods, respectively.

Double Mach Reflection
Here we firstly apply the WENOEL and WENOJS schemes to the 2D double-Mach shock reflection problem where a strong vertical shock moves horizontally into a wedge which inclined with some angle with the ratio of specific heats 1.4 γ = .Initially, this problem was proposed by Woodward and Colella [19] and had been taken extensively as a test example for high-order schemes.The computational domain is chosen to be [ ] [ ] 0, 4 0,1 × and the reflective wall lies on the bottom of the computational domain for 1 4 6 x ≤ ≤ .In the beginning, a Mach 10 shock, moving right, is located at 1 6 x = , 0 y = and makes an angle 60  with the x -axis.For the boundary conditions, the exact postshock condition is imposed for bottom boundary from 0 x = to

Mach 3 Tunnel with a Step
This wind problem is also originally presented in [19].This problem is set up as follows.The wind tunnel is 1 length unit wide and 3 length units long.The step is 0.2 length units high and is located 0.6 length units from the left-hand end of the tunnel.The problem is initialized by a uniform right-going Mach 3 flow.Reflective boundary conditions are applied along the wall of the tunnel.The inflow and outflow boundary conditions are applied at the entrance and exit, respectively.The corner of is a singular point and we use the same technique to correct it [19], which is based on the assumption of a nearly steady flow in the region near the corner.In addition, to eliminate the numerical instability originated from the upwind scheme, the same H-correction procedure as last problem is used to eliminate this flaw.We evolve the initial data until time 4 with a grid of 720 240 × cells by the WENOEL and WENOJS methods.The CFL number is chosen to be 0.6 for these methods.The contour lines of the density are displayed in Figure 11.We observe the good resolution and strong reflective waves in this test for these methods.Moreover, a better resolution can be observed for the WENOJS method.

2D Riemann Problem
This is a simple 2D Riemann problem on square [ ] [ ] 0,1 0,1 × is used to test our method which is originally proposed in [27] 12 gives the density profiles computed by the WENOEL and WENOJS methods, respectively, and these two high-resolution scheme both perform well.

Rayleigh-Taylor Instability Problem
Finally, we consider Rayleigh-Taylor instability problem [28].This problem describes the flow motions on the interface between fluids with different density.The heavy fluid moves into the region of the light fluid with a fingering nature, which lead to the bubbles of light fluid rising into the heavy fluid and the spikes of heavy fluid falling into the light fluid.The computational domain is  (34), respectively.For this problem, from Figure 13, we can find that the WENOEL and WENOJS methods both generate high-resolution numerical solutions.
From the performance of these four 2D problems, as in the 1D problems, we can find that the similar resolution with the WENOJS scheme can be observed.

Conclusion
In this paper, we propose a simple and easily implemented method to solve hyperbolic conversation laws.The main idea of this method is the transformation of integration of flux function in time into integration in space.It is this procedure that leads to the average flux at cell edge which can be directly evaluated.In addition, the evaluation of the average flux is easily implemented and can be combined with any non-oscillatory spatial reconstruction.Through the performance on lots of classical examples, we can find this method is rather robust.We compare our scheme with the classical WENOJS scheme and almost the same performance on resolution is observed.And from the comparisons of resolution and efficiency for these 1D and 2D examples, we can find that the proposed method in this paper can present comparable results with the classical WENOJS scheme, but with less CPU time.

Figure 3 .
Figure 3.The numerical solutions of (25) (26) are computed by the WENOEL and WENOJS methods with 100 N = until time 1 t = .The blue "  " and red " ◊ " denote the solutions of WENOEL and WENOJS methods, respectively.

Figure 4 .
Figure 4.The numerical solutions of (25) (27) are computed by the WENOEL and WENOJS methods with 100 N = until time 0.37 t =

Figure 5 .
Figure 5.The numerical solutions of (28) (29) are computed by the WENOEL and WENOJS methods with 100 N = until time 0.4 t = .The blue "  " and red " ◊ " denote the solutions of WENOEL and WENOJS methods, respectively.

(
problem, a Mach 3 shock wave moved right interacts with sine waves in a density which lead to a field obtained smooth and discontinuous structures.The reference solution in Figure 7 is computed by WENOJS method with cells 2000 N =

Figure 7 .
Figure 7.The shock entropy wave problem (30) (32) is computed by the WENOEL and WENOJS methods with cells 400 N = until time 1.8 t =.The blue "  " and red " ◊ " denote the solutions of WENOEL and WENOJS methods, respectively.

Figure 8 .
Figure 8.The zoomed figure of Figure 7 around the high-frequency waves.

Figure 9 .
Figure 9.The blast wave problem (30) (33) is computed by the WENOEL and WENOJS methods with cells 400 N = until time 0.38 t =.The blue "  " and red " ◊ " denote the solutions of WENOEL and WENOJS methods, respectively.

1 6 x
= , and the reflective boundary condition is imposed for the rest; the flows are imposed on the top boundary such that there is no interaction with the Mach 10 shock; inflow and outflow boundary conditions are set for the left and right boundaries respectively.The unshocked fluid has a density of 1.4, a pressure of 1 and this problem is run until 0.2 t = with cells 960 240 × .In addition, there exists numerical instability in computing this problem by unwind methods.A H-correction procedure presented in Subsection 2.3 is used to eliminate this instability.The H-correction procedure injects much dissipation in unstable region by the way of shifting the upwind scheme into flux splitting scheme.The contour lines of the density are shown Figure 10.These two methods both obtain high-resolution solutions.In region of complicated structure, the WENOJS scheme performs better than our scheme.

Figure 11 .
Figure 11.Mach 3 tunnel with a step problem.The density ρ is shown with meshes 20 240× .30 equally spaced contour lines are plotted from 0.12 to 6.41.The top: WENOEL scheme.The bottom: WENOJS scheme.

Figure 12 .
Figure 12.The solutions of 2D Riemann problem computed by the WENOEL and WENOJS methods.The left: WENOEL scheme.The right: WENOJS scheme.
For the boundary conditions, the reflective boundary conditions are imposed on the left and right boundaries; the flows are set as ( boundaries, respectively.In addition, the source terms ρ and v ρ are added to the right of the third and fourth equations of Euler equations

Figure 13 .
Figure 13.The solutions of Rayleigh-Taylor instability problem computed by the WENOEL and WENOJS methods with cells 100 400 × .The left: WENOEL scheme.The right: WENOJS scheme.

Table 1 .
The test of order of accuracy for the 1 l and l ∞ errors for initial value problem

Table 2 .
(22)test of order of accuracy for the 1 l and l ∞ errors for initial value problem(22)( . The square is divided into four quadrants by straight lines 7 0.8