A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

Self-adjusting entropy-stable scheme for compressible Euler equations

In this work,a self-adjusting entropy-stable scheme is proposed for solving compressible Euler equations.The entropy-stable scheme is constructed by combining the entropy conservative flux with a suitable diffusion operator.The entropy has to be preserved in smooth solutions and be dissipated at shocks.To achieve this,a switch function,which is based on entropy variables,is employed to make the numerical diffusion term be automatically added around discontinuities.The resulting scheme is still entropy-stable.A number of numerical experiments illustrating the robustness and accuracy of the scheme are presented.From these numerical results,we observe a remarkable gain in accuracy.     

A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

An adaptive GRP scheme for compressible fluid flows

This paper presents a second-order accurate adaptive generalized Riemann problem (GRP) scheme for one and two dimensional compressible fluid flows. The current scheme consists of two independent parts: Mesh redistribution and PDE evolution. The first part is an iterative procedure. In each iteration, mesh points are first redistributed, and then a conservative interpolation formula is used to calculate the cell-averages and the slopes of conservative variables on the resulting new mesh. The second part is to evolve the compressible fluid flows on a fixed nonuniform mesh with the Eulerian GRP scheme, which is directly extended to two-dimensional arbitrary quadrilateral meshes. Several numerical examples show that the current adaptive GRP scheme does not only improve the resolution as well as accuracy of numerical solutions with a few mesh points, but also reduces possible errors or oscillations effectively.     

A fully discrete,kinetic energy consistent finite-volume scheme for compressible flows

A robust, implicit, low-dissipation method suitable for LES/DNS of compressible turbulent flows is discussed. The scheme is designed such that the discrete flux of kinetic energy and its rate of change are consistent with those predicted by the momentum and continuity equations. The resulting spatial fluxes are similar to those derived using the so-called skew-symmetric formulation of the convective terms. Enforcing consistency for the time derivative results in a novel density weighted Crank–Nicolson type scheme. The method is stable without the addition of any explicit dissipation terms at very high Reynolds numbers for flows without shocks. Shock capturing is achieved by switching on a dissipative flux term which tends to zero in smooth regions of the flow. Numerical examples include a one-dimensional shock tube problem, the Taylor–Green problem, simulations of isotropic turbulence, hypersonic flow over a double-cone geometry, and compressible turbulent channel flow.     

Positivity-preserving high order finite difference WENO schemes for compressible Euler equations

In Zhang and Shu (2010) [20], Zhang and Shu (2011) [21] and Zhang et al. (in press) [23], we constructed uniformly high order accurate discontinuous Galerkin (DG) and finite volume schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics. In this paper, we present an extension of this framework to construct positivity-preserving high order essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) finite difference schemes for compressible Euler equations. General equations of state and source terms are also discussed. Numerical tests of the fifth order finite difference WENO scheme are reported to demonstrate the good behavior of such schemes.     

A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry

The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes.     

A massively parallel multi-block hybrid compact–WENO scheme for compressible flows

A new multi-block hybrid compact–WENO finite-difference method for the massively parallel computation of compressible flows is presented. In contrast to earlier methods, our approach breaks the global dependence of compact methods by using explicit finite-difference methods at block interfaces and is fully conservative. The resulting method is fifth- and sixth-order accurate for the convective and diffusive fluxes, respectively. The impact of the explicit interface treatment on the stability and accuracy of the multi-block method is quantified for the advection and diffusion equations. Numerical errors increase slightly as the number of blocks is increased. It is also found that the maximum allowable time steps increase with the number of blocks. The method demonstrates excellent scalability on up to 1264 processors.     

A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes

We present a high-order cell-centered Lagrangian scheme for solving the two-dimensional gas dynamics equations on unstructured meshes. A node-based discretization of the numerical fluxes for the physical conservation laws allows to derive a scheme that is compatible with the geometric conservation law (GCL). Fluxes are computed using a nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The first-order scheme is conservative for momentum and total energy, and satisfies a local entropy inequality in its semi-discrete form. The two-dimensional high-order extension is constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess this new scheme. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new scheme.     

A second order numerical scheme for the improved Boussinesq equation

A finite-difference scheme arising from the use of rational approximants to the matrix-exponential term in a three-time level recurrence relation is used for the numerical solution of the improved Boussinesq equation (IBq). The resulting linear scheme, which is analyzed for local truncation error and stability, is tested numerically and conclusions with corresponding results known in the bibliography are derived.     

A multislope MUSCL method on unstructured meshes applied to compressible Euler equations for axisymmetric swirling flows

A finite volume method for the numerical solution of axisymmetric inviscid swirling flows is presented. The governing equations of the flow are the axisymmetric compressible Euler equations including swirl (or tangential) velocity. A first-order scheme is introduced where the convective fluxes at cell interfaces are evaluated by the Rusanov or the HLLC numerical flux while the geometric source terms are discretizated to provide a well-balanced scheme i.e. the steady-state solutions with null velocity are preserved. Extension to the second-order space approximation using a multislope MUSCL method is then derived. To test the numerical scheme, a stationary solution of the fluid flow following the radial direction has been established with a zero and nonzero tangential velocity. Numerical and exact solutions are compared for classical Riemann problems where we employ different limiters and effectiveness of the multislope MUSCL scheme is demonstrated for strongly shocked axially symmetric flows like in spherical bubble compression problem. Two other tests with axisymmetric geometries are performed: the supersonic flow in a tube with a cone and the axisymmetric blunt body with a free stream.     

Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms

In [16], [17], we constructed uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics with the ideal gas equation of state. The technique also applies to high order accurate finite volume schemes. For the Euler equations with various source terms (e.g., gravity and chemical reactions), it is more difficult to design high order schemes which do not produce negative density or pressure. In this paper, we first show that our framework to construct positivity-preserving high order schemes in [16], [17] can also be applied to Euler equations with a general equation of state. Then we discuss an extension to Euler equations with source terms. Numerical tests of the third order Runge–Kutta DG (RKDG) method for Euler equations with different types of source terms are reported.     

On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes

We construct uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for Euler equations of compressible gas dynamics. The same framework also applies to high order accurate finite volume (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO)) schemes. Motivated by Perthame and Shu (1996) [20] and Zhang and Shu (2010) [26], a general framework, for arbitrary order of accuracy, is established to construct a positivity preserving limiter for the finite volume and DG methods with first order Euler forward time discretization solving one-dimensional compressible Euler equations. The limiter can be proven to maintain high order accuracy and is easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the positivity property. Following the idea in Zhang and Shu (2010) [26], we extend this framework to higher dimensions on rectangular meshes in a straightforward way. Numerical tests for the third order DG method are reported to demonstrate the effectiveness of the methods.     

Implementation of the GRP scheme for computing radially symmetric compressible fluid flows

The study of radially symmetric compressible fluid flows is interesting both from the theoretical and numerical points of view. Spherical explosion and implosion in air, water and other media are well-known problems in application. Typical difficulties lie in the treatment of singularity in the geometrical source and the imposition of boundary conditions at the symmetric center, in addition to the resolution of classical discontinuities (shocks and contact discontinuities). In the present paper we present the implementation of direct generalized Riemann problem (GRP) scheme to resolve this issue. The scheme is obtained directly by the time integration of the fluid flows. Our new contribution is to show rigorously that the singularity is removable and derive the updating formulae for mass and energy at the center. Together with the vanishing of the momentum, we obtain new numerical boundary conditions at the center, which are then incorporated into the GRP scheme. The main ingredient is the passage from the Cartesian coordinates to the radially symmetric coordinates.     

A high order kinetic flux-vector splitting method for the reduced five-equation model of compressible two-fluid flows

We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) [21]. The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two space dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge–Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor [14]. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.     

Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows

The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant–Friedrichs–Lewy is kept near unity.     

A turbulence characteristic length scale for compressible flows

The current RANS models are generally established and calibrated under incompressible condition and these kinds of models could succeed in predicting many features of incompressible flows. However, these models extended to the high-speed, compressible flows are always less accurate. In the paper, a compressible von Kármán length scale is proposed for compressible flows considering the variable densities. It contains no empirical coefficients and is based on phenomenological theory. In the turbulent kinetic equation, the extra unclosed terms induced by non-constant densities are treated as dissipation terms and the equation is closed algebraically via the introduction of the von Kármán length scale. The original and the proposed von Kármán length scale lead to two different kinds of SAS (scale adaption simulation) models, KDO (turbulence kinetic energy dependent only) and CKDO (compressible KDO), respectively. Compressible mixing layer with significant compressibility is studied within standard k–?, k–ω, KDO turbulence models and their compressible versions. The compressibility effects such as the reduced mixing layer thickness, growth rate and turbulence intensity can be reproduced by CKDO. The new length scale can improve the performances of the model in predicting the mixing layer thickness, stream-wise velocity and Reynolds shear stresses when the convective Mach number is 0.8. Besides, the new length scale also leads to accurate computed growth rate when the convective Mach number ranges from 0.1 to 1.0.     

A smart nonstandard finite difference scheme for second order nonlinear boundary value problems

A new kind of finite difference scheme is presented for special second order nonlinear two point boundary value problems. An artificial parameter is introduced in the scheme. Symbolic computation is proposed for the construction of the scheme. Local truncation error of the method is discussed. Numerical examples are illustrated. Numerical results show that the method is very effective.     

The effect of shocks on second order sensitivities for the quasi-one-dimensional Euler equations

The effect of discontinuity in the state variables on optimization problems is investigated on the quasi-one-dimensional Euler equations in the discrete level. A pressure minimization problem and a pressure matching problem are considered. We find that the objective functional can be smooth in the continuous level and yet be non-smooth in the discrete level as a result of the shock crossing grid points. Higher resolution can exacerbate that effect making grid refinement counter productive for the purpose of computing the discrete sensitivities. First and second order sensitivities, as well as the adjoint solution, are computed exactly at the shock and its vicinity and are compared to the continuous solution. It is shown that in the discrete level the first order sensitivities contain a spike at the shock location that converges to a delta function with grid refinement, consistent with the continuous analysis. The numerical Hessian is computed and its consistency with the analytical Hessian is discussed for different flow conditions. It is demonstrated that consistency is not guaranteed for shocked flows. We also study the different terms composing the Hessian and propose some stable approximation to the continuous Hessian.