This paper presents a new approach, so-called boundary variation diminishing (BVD), for reconstructions that minimize the discontinuities (jumps) at cell interfaces in Godunov type schemes. It is motivated by the observation that diminishing the jump at the cell boundary might effectively reduce the dissipation in numerical flux. Different from the existing practices which seek high-order polynomials within mesh cells while assuming discontinuities being always at the cell interfaces, we proposed a new strategy that combines a high-order polynomial-based interpolation and a jump-like reconstruction that allows a discontinuity being partly represented within the mesh cell rather than at the interface. It is shown that new schemes of high fidelity for both continuous and discontinuous solutions can be devised by the BVD guideline with properly-chosen candidate reconstruction schemes. Excellent numerical results have been obtained for both scalar and Euler conservation laws with substantially improved solution quality in comparison with the existing methods. This work provides a simple and accurate alternative with great practical significance to the current Godunov paradigm which overly pursues the smoothness within mesh cell under the questionable premiss that discontinuities only appear at cell interfaces.

No abstract

Two-dimensional, oblique detonations induced by a wedge are simulated using the reactive Euler equations with a detailed chemical reaction model. The focus of this study is on the oblique shock-to-detonation transition in a stoichiometric hydrogenair mixture. A combustible, gas mixture at low pressure and high temperature, corresponding to the realistic, inflow conditions applied in oblique detonation wave engines, is presented in this study. At practical flight conditions, the present numerical results illustrate that oblique detonation initiation is achieved through a smooth transition from a curved shock, which differs from the abrupt transition depicted in the previous studies. The formation mechanism of this smooth transition is discussed and a quantitative analysis is carried out by defining a characteristic length for the initiation process. The dependence of the initiation length on different parameters including the wedge angle, flight Mach number, and inflow Mach number is discussed. Despite the hypothetical nature of the simulation configuration, the present numerical study uses parameters we deem relevant to practical conditions and provides important observations for which future investigations can benefit from in reaching toward a rigorous theory of the formation and self-sustenance of oblique detonation waves. C 2015 AIP Publishing LLC. [http://dx

SUMMARYNumerical oscillation has been an open problem for high order numerical methods with increased local degrees of freedom (DOFs). Current strategies mainly follow the limiting projections derived originally for conventional finite volume methods, and thus are not able to make full use of the sub-cell information available in the local high order reconstructions. This paper presents a novel algorithm which introduces a nodal-value based weighted essentially non-oscillatory limiter for constrained interpolation profile/multimoment finite volume method (CIP/MM FVM) (Ii and Xiao, J. Comput. Phys., 222(2007), 849-871 ) as an effort to pursue a better suited formulation to implement the limiting projection in schemes with local DOFs. The new scheme, CIP-CSL-WENO4 scheme, extends the CIP/MM FVM method by limiting the slope constraint in the interpolation function using the WENO reconstruction which makes use of the subcell information available from the local DOFs and is built from the point values (PVs) at the solution points within three neighboring cells, thus resulting a more compact WENO stencil. The proposed WENO limiter matches well the original CIP/MM FVM, which leads to a new scheme of high accuracy, algorithmic simplicity and computational efficiency. We present the numerical results of benchmark tests for both scalar and Euler conservation laws to manifest the 4th order accuracy and oscillation-suppressing property of the proposed scheme.

In this work we propose a new formulation for high-order multi-moment constrained finite volume (MCV) method. In the one-dimensional building-block scheme, three local degrees of freedom (DOFs) are equidistantly defined within a grid cell. Two candidate polynomials for spatial reconstruction of third-order are built by adopting one additional constraint condition from the adjacent cells, i.e. the DOF at middle point of left or right neighbour. A boundary gradient switching (BGS) algorithm based on the variation-minimization principle is devised to determine the spatial reconstruction from the two candidates, so as to remove the spurious oscillations around the discontinuities. The resulted non-oscillatory MCV3-BGS scheme is of fourth-order accuracy and completely free of case-dependent ad hoc parameters. The widely used benchmark tests of one-and two-dimensional scalar and Euler hyperbolic conservation laws are solved to verify the performance of the proposed scheme in this paper. The MCV3-BGS scheme is very promising for the practical applications due to its accuracy, non-oscillatory feature and algorithmic simplicity.