Search

This paper serves to simulate compressible multi-fluid flows using an energy-splitting method based on the generalized Riemann problem (GRP) solver. The emphasis is put on two-dimensional problems to address the fundamental role of transversal effect in the construction of numerical fluxes, which exhibits the genuine multi-dimensionality of the underlying numerical schemes. In contrast with the popularly-used methods of line which use exact or approximate Riemann solvers, the GRP solver, as the discontinuous version of the Lax–Wendroff method, takes the temporal-spatial coupling process to include the transversal effect. Combined with the energy-splitting technique developed earlier in Lei and Li [24], the resolution of flow fields around material interfaces is improved significantly.

This paper serves to treat boundary conditions numerically with high order accuracy in order to match the two-stage fourth-order finite volume schemes for hyperbolic problems developed in [{\em J. Li and Z. Du, A two-stage fourth order time-accurate discretization {L}ax--{W}endroff type flow solvers, {I}. {H}yperbolic conservation laws, SIAM, J. Sci. Comput., 38 (2016), pp.~A3046--A3069}]. As such, it is significant when capturing small scale structures near physical boundaries. Different from previous contributions in literature, the current approach constructs a fourth order accurate approximation to boundary conditions by only using the Jacobian of the flux function (characteristic information) instead of its successive differentiation leading to tensors of high ranks in the inverse Lax-Wendroff method. Technically, data in several ghost cells are constructed with interpolation so that the interior scheme can be implemented over boundary cells, and theoretical boundary condition has to be modified properly at intermediate stages so as to make the two-stage scheme over boundary cells fully consistent with that over interior cells. This highlights the fact that {\em continuous boundary conditions only match continuous partial differential equations (PDEs), and they must be approximated in a consistent way (even though it could be exactly valued) when the PDEs are discretized.} Several numerical examples are provided to illustrate the performance of the current approach when dealing with general boundary conditions.

This paper develops a new fifth order accurate Hermite WENO (HWENO) reconstruction method for hyperbolic conservation schemes in the framework of the two-stage fourth order accurate temporal discretization in [{\em J. Li and Z. Du, A two-stage fourth order time-accurate discretization {L}ax--{W}endroff type flow solvers, {I}. {H}yperbolic conservation laws, SIAM, J. Sci. Comput., 38 (2016), pp.~A3046--A3069}]. Instead of computing the first moment of the solution additionally in the conventional HWENO or DG approach, we can directly take the {\em interface values}, which are already available in the numerical flux construction using the generalized Riemann problem (GRP) solver, to approximate the first moment. The resulting scheme is fourth order temporal accurate by only invoking the HWENO reconstruction twice so that it becomes more compact. Numerical experiments show that such compactness makes significant impact on the resolution of nonlinear waves.

In this paper we develop a novel two-stage fourth order time-accurate discretization for time-dependent flow problems, particularly for hyperbolic conservation laws. Different from the classical Runge-Kutta (R-K) temporal discretization for first order Riemann solvers as building blocks, the current approach is solely associated with Lax-Wendroff (L-W) type schemes as the building blocks. As a result, a two-stage procedure can be constructed to achieve a fourth order temporal accuracy, rather than using well-developed four stages for R-K methods. The generalized Riemann problem (GRP) solver is taken as a representative of L-W type schemes for the construction of a two-stage fourth order scheme., 24 pages