Riemann solvers with evolved initial conditions

The scope of this paper is three fold. We first formulate upwind and symmetric schemes for hyperbolic equations with non‐conservative terms. Then we propose upwind numerical schemes for conservative and non‐conservative systems, based on a Riemann solver, the initial conditions of which are evolved non‐linearly in time, prior to a simple linearization that leads to closed‐form solutions. The Riemann solver is easily applied to complicated hyperbolic systems. Finally, as an example, we formulate conservative schemes for the three‐dimensional Euler equations for general compressible materials and give numerical results for a variety of test problems for ideal gases in one and two space dimensions. Copyright © 2006 John Wiley & Sons, Ltd.     

A Riemann solver and upwind methods for a two‐phase flow model in non‐conservative form

We present a theoretical solution for the Riemann problem for the five‐equation two‐phase non‐conservative model of Saurel and Abgrall. This solution is then utilized in the construction of upwind non‐conservative methods to solve the general initial‐boundary value problem for the two‐phase flow model in non‐conservative form. The basic upwind scheme constructed is the non‐conservative analogue of the Godunov first‐order upwind method. Second‐order methods in space and time are then constructed via the MUSCL and ADER approaches. The methods are systematically assessed via a series of test problems with theoretical solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Efficient construction of high‐resolution TVD conservative schemes for equations with source terms: application to shallow water flows

High‐resolution total variation diminishing (TVD) schemes are widely used for the numerical approximation of hyperbolic conservation laws. Their extension to equations with source terms involving spatial derivatives is not obvious. In this work, efficient ways of constructing conservative schemes from the conservative, non‐conservative or characteristic form of the equations are described in detail. An upwind, as opposed to a pointwise, treatment of the source terms is adopted here, and a new technique is proposed in which source terms are included in the flux limiter functions to get a complete second‐order compact scheme. A new correction to fix the entropy problem is also presented and a robust treatment of the boundary conditions according to the discretization used is stated. Copyright © 2001 John Wiley & Sons, Ltd.     

Local maximum principle satisfying high‐order non‐oscillatory schemes

The main contribution of this work is to classify the solution region including data extrema for which high‐order non‐oscillatory approximation can be achieved. It is performed in the framework of local maximum principle (LMP) and non‐conservative formulation. The representative uniformly second‐order accurate schemes are converted in to their non‐conservative form using the ratio of consecutive gradients. Using the local maximum principle, these non‐conservative schemes are analyzed for their non‐linear LMP/total variation diminishing stability bounds which classify the solution region where high‐order accuracy can be achieved. Based on the bounds, second‐order accurate hybrid numerical schemes are constructed using a shock detector. The presented numerical results show that these hybrid schemes preserve high accuracy at non‐sonic extrema without exhibiting any induced local oscillations or clipping error. Copyright © 2015 John Wiley & Sons, Ltd.     

Composite high resolution localized relaxation scheme based on upwinding for hyperbolic conservation laws

In this work we present an upwind‐based high resolution scheme using flux limiters. Based on the direction of flow we choose the smoothness parameter in such a way that it leads to a truly upwind scheme without losing total variation diminishing (TVD) property for hyperbolic linear systems where characteristic values can be of either sign. Here we present and justify the choice of smoothness parameters. The numerical flux function of a high resolution scheme is constructed using wave speed splitting so that it results into a scheme that truly respects the physical hyperbolicity property. Bounds are given for limiter functions to satisfy TVD property. The proposed scheme is extended for non‐linear problems by using the framework of relaxation system that converts a non‐linear conservation law into a system of linear convection equations with a non‐linear source term. The characteristic speed of relaxation system is chosen locally on three point stencil of grid. This obtained relaxation system is solved using composite scheme technique, i.e. using a combination of proposed scheme with the conservative non‐standard finite difference scheme. Presented numerical results show higher resolution near discontinuity without introducing spurious oscillations. Copyright © 2008 John Wiley & Sons, Ltd.     

High‐order ENO and WENO schemes for unstructured grids

This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd.     

MUSTA schemes for multi‐dimensional hyperbolic systems: analysis and improvements

We develop and analyse an improved version of the multi‐stage (MUSTA) approach to the construction of upwind Godunov‐type fluxes whereby the solution of the Riemann problem, approximate or exact, is not required. The new MUSTA schemes improve upon the original schemes in terms of monotonicity properties, accuracy and stability in multiple space dimensions. We incorporate the MUSTA technology into the framework of finite‐volume weighted essentially nonoscillatory schemes as applied to the Euler equations of compressible gas dynamics. The results demonstrate that our new schemes are good alternatives to current centred methods and to conventional upwind methods as applied to complicated hyperbolic systems for which the solution of the Riemann problem is costly or unknown. Copyright © 2005 John Wiley & Sons, Ltd.     

A variable relaxation scheme for multiphase,multicomponent flow

We propose a variable relaxation scheme for the simulation of 1D, two-phase, multicomponent flow in porous media. For these strongly nonlinear systems, traditional high order upwind schemes are impractical: Riemann solutions are not directly available when the phase behavior is complex, and the systems are weakly hyperbolic at isolated points. Relaxation schemes avoid the dependency on the eigenstructure and nonlinear Riemann solvers by approximating the original system with a strongly hyperbolic linear system. We exploit the known information about the eigenvalues to construct first order and second order variable relaxation schemes with much reduced numerical diffusion as compared to the standard relaxation formulations. The proposed second order variable relaxation scheme is competitive in accuracy and efficiency with a third order component-wise ENO reconstruction, and performs at least as well as second order component-wise TVD schemes.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

Implicit weighted essentially non‐oscillatory schemes for the incompressible Navier–Stokes equations

A class of lower–upper/approximate factorization (LUAF) implicit weighted essentially non‐oscillatory (ENO; WENO) schemes for solving the two‐dimensional incompressible Navier–Stokes equations in a generalized co‐ordinate system is presented. The algorithm is based on the artificial compressibility formulation, and symmetric Gauss–Seidel relaxation is used for computing steady state solutions while symmetric successive overrelaxation is used for treating time‐dependent flows. WENO spatial operators are employed for inviscid fluxes and central differencing for viscous fluxes. Internal and external viscous flow test problems are presented to verify the numerical schemes. The use of a WENO spatial operator not only enhances the accuracy of solutions but also improves the convergence rate for the steady state computation as compared with using the ENO counterpart. It is found that the present solutions compare well with exact solutions, experimental data and other numerical results. Copyright © 1999 John Wiley & Sons, Ltd.     

Developing shock-capturing difference methods

A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the characteristic decompositions are carefully controlled by limiters to prevent non-physical oscillations. Lastly, the schemes are converted into conservative forms, and the oscillation-free shock-capturing schemes are acquired. Two explicit upwind schemes (2nd-order and 3rd-order) and three compact upwind schemes (3rd-order, 5th-order and 7th-order) are modified by the method for hyperbolic systems and the modified schemes are checked on several one-dimensional and two-dimensional test cases. Some numerical solutions of the schemes are compared with those of a WENO scheme and a MP scheme as well as a compact-WENO scheme. The results show that the method with high order accuracy and high resolutions can capture shock waves smoothly.     

Adaptive mesh refinement techniques for high‐order shock capturing schemes for multi‐dimensional hydrodynamic simulations

The numerical simulation of physical phenomena represented by non‐linear hyperbolic systems of conservation laws presents specific difficulties mainly due to the presence of discontinuities in the solution. State of the art methods for the solution of such equations involve high resolution shock capturing schemes, which are able to produce sharp profiles at the discontinuities and high accuracy in smooth regions, together with some kind of grid adaption, which reduces the computational cost by using finer grids near the discontinuities and coarser grids in smooth regions. The combination of both techniques presents intrinsic numerical and programming difficulties. In this work we present a method obtained by the combination of a high‐order shock capturing scheme, built from Shu–Osher's conservative formulation (J. Comput. Phys. 1988; 77 :439–471; 1989; 83 :32–78), a fifth‐order weighted essentially non‐oscillatory (WENO) interpolatory technique (J. Comput. Phys. 1996; 126 :202–228) and Donat–Marquina's flux‐splitting method (J. Comput. Phys. 1996; 125 :42–58), with the adaptive mesh refinement (AMR) technique of Berger and collaborators (Adaptive mesh refinement for hyperbolic partial differential equations. Ph.D. Thesis, Computer Science Department, Stanford University, 1982; J. Comput. Phys. 1989; 82 :64–84; 1984; 53 :484–512). Copyright © 2006 John Wiley & Sons, Ltd.     

Augmented upwind numerical schemes for the groundwater transport advection‐dispersion equation with local operators

When solute transport is advection‐dominated, the advection‐dispersion equation approximates to a hyperbolic‐type partial differential equation, and finite difference and finite element numerical approximation methods become prone to artificial oscillations. The upwind scheme serves to correct these responses to produce a more realistic solution. The upwind scheme is reviewed and then applied to the advection‐dispersion equation with local operators for the first‐order upwinding numerical approximation scheme. The traditional explicit and implicit schemes, as well as the Crank‐Nicolson scheme, are developed and analyzed for numerical stability to form a comparison base. Two new numerical approximation schemes are then proposed, namely, upwind–Crank‐Nicolson scheme, where only for the advection term is applied, and weighted upwind‐downwind scheme. These newly developed schemes are analyzed for numerical stability and compared to the traditional schemes. It was found that an upwind–Crank‐Nicolson scheme is appropriate if the Crank‐Nicolson scheme is only applied to the advection term of the advection‐dispersion equation. Furthermore, the proposed explicit weighted upwind‐downwind finite difference numerical scheme is an improvement on the traditional explicit first‐order upwind scheme, whereas the implicit weighted first‐order upwind‐downwind finite difference numerical scheme is stable under all assumptions when the appropriate weighting factor (θ) is assigned.     

Numerical simulation of unsteady multidimensional free surface motions by level set method

This paper presents a numerical method that couples the incompressible Navier–Stokes equations with the level set method in a curvilinear co‐ordinate system for study of free surface flows. The finite volume method is used to discretize the governing equations on a non‐staggered grid with a four‐step fractional step method. The free surface flow problem is converted into a two‐phase flow system on a fixed grid in which the free surface is implicitly captured by the zero level set. We compare different numerical schemes for advection of the level set function in a generalized curvilinear format, including the third order quadratic upwind interpolation for convective kinematics (QUICK) scheme, and the second and third order essentially non‐oscillatory (ENO) schemes. The level set equations of evolution and reinitialization are validated with benchmark cases, e.g. a stationary circle, a rotating slotted disk and stretching of a circular fluid element. The coupled system is then applied to a travelling solitary wave, and two‐ and three‐dimensional dam breaking problems. Some interesting free surface phenomena are revealed by the computational results, such as, the large free surface vortices, air entrapment and splashing of the water surge front. The computational results are in excellent agreement with theoretical predictions and experimental data, where they are available. Copyright © 2003 John Wiley & Sons, Ltd.     

Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

Third-order modified coefficient scheme based on essentially non-oscillatory scheme

A third-order numerical scheme is presented to give approximate solutions to multi-dimensional hyperbolic conservation laws only using modified coefficients of an essentially non-oscillatory （MCENO） scheme without increasing the base points during construction of the scheme. The construction process shows that the modified coefficient approach preserves favourable properties inherent in the original essentially nonoscillatory （ENO） scheme for its essential non-oscillation, total variation bounded （TVB）, etc. The new scheme improves accuracy by one order compared to the original one. The proposed MCENO scheme is applied to simulate two-dimensional Rayleigh-Taylor （RT） instability with densities 1：3 and 1：100, and solve the Lax shock-wave tube numerically. The ratio of CPU time used to implement MCENO, the .third-order ENO and fifth-order weighed ENO （WENO） schemes is 0.62：1：2.19. This indicates that MCENO improves accuracy in smooth regions and has higher accuracy and better efficiency compared to the original ENO scheme.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

On the discontinuous Galerkin computation of N‐waves focusing

Sonic boom focusing phenomenon can be predicted using the solution to the nonlinear Tricomi equation which is a hybrid (hyperbolic‐elliptic) second‐order partial differential equation. In this paper, the hyperbolic conservation law form is derived, which is valid in the entire domain. In this manner, the presence of two regions where the equation behaves differently (hyperbolic in the upper and elliptic in the lower half‐plane) is avoided. On the upper boundary, a new mixed boundary condition for the acoustic pressure is employed. The discretization is carried out using a discontinuous Galerkin (DG) method combined with a Runge–Kutta total‐variation diminishing scheme. The results show the accuracy of DG methods to solve problems involving sharp gradients and discontinuities. Comparisons with analytical results for the linear case, and other numerical results using classical explicit and compact finite difference schemes and weighted essentially non‐oscillatory schemes are included. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluation of some high‐order shock capturing schemes for direct numerical simulation of unsteady two‐dimensional free flows

The present study addresses the capability of a large set of shock‐capturing schemes to recover the basic interactions between acoustic, vorticity and entropy in a direct numerical simulation (DNS) framework. The basic dispersive and dissipative errors are first evaluated by considering the advection of a Taylor vortex in a uniform flow. Two transonic cases are also considered. The first one consists of the interaction between a temperature spot and a weak shock. This test emphasizes the capability of the schemes to recover the production of vorticity through the baroclinic process. The second one consists of the interaction of a Taylor vortex with a weak shock, corresponding to the framework of the linear theory of Ribner. The main process in play here is the production of an acoustic wave. The results obtained by using essentially non‐oscillatory (ENO), total variation diminishing (TVD), compact‐TVD and MUSCL schemes are compared with those obtained by means of a sixth‐order accurate Hermitian scheme, considered as reference. The results are as follows; the ENO schemes agree pretty well with the reference scheme. The second‐order accurate Upwind‐TVD scheme exhibits a strong numerical diffusion, while the MUSCL scheme behavior is very sensitive to the value on the parameter β in the limiter function minmod. The compact‐TVD schemes do not yield improvement over the standard TVD schemes. Copyright © 2000 John Wiley & Sons, Ltd.