High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws

In [10], the authors have designed a new fifth order WENO finite-difference scheme by adding a higher order smoothness indicator which is obtained as a simple and inexpensive linear combination of the already existing low order smoothness indicators. Moreover, this new scheme, dubbed as WENO-Z, has a CPU cost which is equivalent to the one of the classical WENO-JS [2], and smaller than that of the mapped WENO-M, [5], since it involves no mapping of the nonlinear weights. In this article, we take a closer look at Taylor expansions of the Lagrangian polynomials of the WENO substencils and the related inherited symmetries of the classical lower order smoothness indicators to obtain a general formula for the higher order smoothness indicators that allows the extension of the WENO-Z scheme to all (odd) orders of accuracy. We further investigate the improved accuracy of the WENO-Z schemes at critical points of smooth solutions as well as their distinct numerical features as a result of the new sets of nonlinear weights and we show that regarding the numerical dissipation WENO-Z occupies an intermediary position between WENO-JS and WENO-M. Some standard numerical experiments such as the one dimensional Riemann initial values problems for the Euler equations and the Mach 3 shock density-wave interaction and the two dimensional double-Mach shock reflection problems are presented.     

An adaptive central-upwind weighted essentially non-oscillatory scheme

In this work, an adaptive central-upwind 6th-order weighted essentially non-oscillatory (WENO) scheme is developed. The scheme adapts between central and upwind schemes smoothly by a new weighting relation based on blending the smoothness indicators of the optimal higher order stencil and the lower order upwind stencils. The scheme achieves 6th-order accuracy in smooth regions of the solution by introducing a new reference smoothness indicator. A number of numerical examples suggest that the present scheme, while preserving the good shock-capturing properties of the classical WENO schemes, achieves very small numerical dissipation.     

Hybrid weighted essentially non-oscillatory schemes with different indicators

A key idea in finite difference weighted essentially non-oscillatory (WENO) schemes is a combination of lower order fluxes to obtain a higher order approximation. The choice of the weight to each candidate stencil, which is a nonlinear function of the grid values, is crucial to the success of WENO schemes. For the system case, WENO schemes are based on local characteristic decompositions and flux splitting to avoid spurious oscillation. But the cost of computation of nonlinear weights and local characteristic decompositions is very high. In this paper, we investigate hybrid schemes of WENO schemes with high order up-wind linear schemes using different discontinuity indicators and explore the possibility in avoiding the local characteristic decompositions and the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong shocks. The idea is to identify discontinuity by an discontinuity indicator, then reconstruct numerical flux by WENO approximation in discontinuous regions and up-wind linear approximation in smooth regions. These indicators are mainly based on the troubled-cell indicators for discontinuous Galerkin (DG) method which are listed in the paper by Qiu and Shu (J. Qiu, C.-W. Shu, A comparison of troubled-cell indicators for Runge–Kutta discontinuous Galerkin methods using weighted essentially non-oscillatory limiters, SIAM Journal of Scientific Computing 27 (2005) 995–1013). The emphasis of the paper is on comparison of the performance of hybrid scheme using different indicators, with an objective of obtaining efficient and reliable indicators to obtain better performance of hybrid scheme to save computational cost. Detail numerical studies in one- and two-dimensional cases are performed, addressing the issues of efficiency (less CPU time and more accurate numerical solution), non-oscillatory property.     

A Hermite upwind WENO scheme for solving hyperbolic conservation laws

This paper proposes a new WENO procedure to compute problems containing both discontinuities and a large disparity of characteristic scales.In a one-dimensional context, the WENO procedure is defined on a three-points stencil and designed to be sixth-order in regions of smoothness. We define a finite-volume discretization in which we consider the cell averages of the variable and its first derivative as discrete unknowns. The reconstruction of their point-values is then ensured by a unique sixth-order Hermite polynomial. This polynomial is considered as a symmetric and convex combination, by ideal weights, of three fourth-order polynomials: a central polynomial, defined on the three-points stencil, is combined with two polynomials based on the left and the right two-points stencils.The symmetric nature of such an interpolation has an important consequence: the choice of ideal weights has no influence on the properties of the discretization. This advantage enables to formulate the Hermite interpolation for non-uniform meshes. Following the methodology of the classic WENO procedure, nonlinear weights are then defined.To deal with the peculiarities of the Hermite interpolation near discontinuities, we define a new procedure in order for the nonlinear weights to smoothly evolve between the ideal weights, in regions of smoothness, and one-sided weights, otherwise.The resulting scheme is a sixth-order WENO method based on central Hermite interpolation and TVD Runge–Kutta time-integration. We call this scheme the HCWENO6 scheme.Numerical experiments in the scalar and the 1D Euler cases make it possible to check and to validate the options selected. In these experiments, we emphasize the resolution power of the method by computing test cases that model realistic aero-acoustic problems.     

Boundary closures for fourth-order energy stable weighted essentially non-oscillatory finite-difference schemes

A general strategy was presented in 2009 by Yamaleev and Carpenter for constructing energy stable weighted essentially non-oscillatory (ESWENO) finite-difference schemes on periodic domains. ESWENO schemes up to eighth order were developed that are stable in the energy norm for systems of linear hyperbolic equations. Herein, boundary closures are developed for the fourth-order ESWENO scheme that maintain, wherever possible, the WENO stencil biasing properties and satisfy the summation-by-parts (SBP) operator convention, thereby ensuring stability in an L₂ norm. Second-order and third-order boundary closures are developed that are stable in diagonal and block norms, respectively, and achieve third- and fourth-order global accuracy for hyperbolic systems. A novel set of nonuniform flux interpolation points is necessary near the boundaries to simultaneously achieve (1) accuracy, (2) the SBP convention, and (3) WENO stencil biasing mechanics.     

A central WENO scheme for solving hyperbolic conservation laws on non-uniform meshes

This paper proposes a new WENO procedure to compute multi-scale problems with embedded discontinuities, on non-uniform meshes.In a one-dimensional context, the WENO procedure is first defined on a five-points stencil and designed to be fifth-order accurate in regions of smoothness. To this end, we define a finite-volume discretization in which we consider the cell averages of the variable as the discrete unknowns. The reconstruction of their point-values is then ensured by a unique fifth-order polynomial. This optimum polynomial is considered as a symmetric and convex combination, by ideal weights, of four quadratic polynomials.The symmetric nature of the resulting interpolation has an important consequence: the choice of ideal weights has no influence on the accuracy of the discretization. This advantage enables to formulate the interpolation for non-uniform meshes. Following the methodology of the classic WENO procedure, non-oscillatory weights are then calculated from the ideal weights.We adapt this procedure for the non-linear weights to maintain the theoretical convergence properties of the optimum reconstruction, whatever the problem considered.The resulting scheme is a fifth-order WENO method based on central interpolation and TVD Runge–Kutta time-integration. We call this scheme the CWENO5 scheme.Numerical experiments in the scalar and the 1D Euler cases make it possible to check and to validate the options selected. In those experiments, we emphasize the resolution power of the method by computing test cases that model realistic aero-acoustic problems. Finally, the new algorithm is directly extended to bi-dimensional problems.     

On large time step TVD scheme for hyperbolic conservation laws and its efficiency evaluation

A large time step (LTS) TVD scheme originally proposed by Harten is modified and further developed in the present paper and applied to Euler equations in multidimensional problems. By firstly revealing the drawbacks of Harten’s original LTS TVD scheme, and reasoning the occurrence of the spurious oscillations, a modified formulation of its characteristic transformation is proposed and a high resolution, strongly robust LTS TVD scheme is formulated. The modified scheme is proven to be capable of taking larger number of time steps than the original one. Following the modified strategy, the LTS TVD schemes for Yee’s upwind TVD scheme and Yee–Roe–Davis’s symmetric TVD scheme are constructed. The family of the LTS schemes is then extended to multidimensional by time splitting procedure, and the associated boundary condition treatment suitable for the LTS scheme is also imposed. The numerical experiments on Sod’s shock tube problem, inviscid flows over NACA0012 airfoil and ONERA M6 wing are performed to validate the developed schemes. Computational efficiencies for the respective schemes under different CFL numbers are also evaluated and compared. The results reveal that the improvement is sizable as compared to the respective single time step schemes, especially for the CFL number ranging from 1.0 to 4.0.     

Long-time effect of relaxation for hyperbolic conservation laws

In processes such as invasion percolation and certain models of continuum percolation, in which a possibly random labelf(b) is attached to each bondb of a possibly random graph, percolation models for various values of a parameterr are naturally coupled: one can define a bondb to be occupied at levelr iff(b)r. If the labeled graph is stationary, then under the mild additional assumption of positive finite energy, a result of Gandolfi, Keane, and Newman ensures that, in lattice models, for each fixedr at which percolation occurs, the infinite cluster is unique a.s. Analogous results exist for certain continuum models. A unifying framework is given for such fixed-r results, and it is shown that if the site density is finite and the labeled graph has positive finite energy, then with probability one, uniqueness holds simultaneously for all values ofr. An example is given to show that when the site density is infinite, positive finite energy does not ensure uniqueness, even for fixedr. In addition, with finite site density but without positive finite energy, one can have fixed-r uniqueness a.s. for eachr, yet not have simultaneous uniqueness.Research supported by NSF grant DMS-9206139     

An efficient local time-stepping scheme for solution of nonlinear conservation laws

We develop an efficient local time-stepping algorithm for the method of lines approach to numerical solution of transient partial differential equations. The need for local time-stepping arises when adaptive mesh refinement results in a mesh containing cells of greatly different sizes. The global CFL number and, hence, the global time step, are defined by the smallest cell size. This can be inefficient as a few small cells may impose a restrictive time step on the whole mesh. A local time-stepping scheme allows us to use the local CFL number which reduces the total number of function evaluations. The algorithm is based on a second order Runge–Kutta time integration. Its important features are a small stencil and the second order accuracy in the L² and L^∞ norms.     

Multiple-mode diffusion waves for viscous nonstrictly hyperbolic conservation laws

We study the large-time behaviors of solutions of viscous conservation laws whose inviscid part is a nonstrictly hyperbolic system. The initial data considered here is a perturbation of a constant state. It is shown that the solutions converge to single-mode diffusion waves in directions of strictly hyperbolic fields, and to multiple-mode diffusion waves in directions of nonstrictly hyperbolic fields. The multiple-mode diffusion waves, which are the new elements here, are the self-similar solutions of the viscous conservation laws projected to the nonstrictly hyperbolic fields, with the nonlinear fluxes replaced by their quadratic parts. The convergence rate to these diffusion waves isO(t ^–3/4+1/2p+) inL ^p , 1p, with >0 being arbitrarily small.This work was supported by the Applied Mathematical Sciences subprogram of the Office of Energy Research, U.S. Department of Energy, under Contract W-31-109-Eng-38     

Multi-dimensional limiting process for hyperbolic conservation laws on unstructured grids

The present paper deals with an efficient and accurate limiting strategy for the multi-dimensional hyperbolic conservation laws on unstructured grids. The multi-dimensional limiting process (MLP) which has been successfully proposed on structured grids is extended to unstructured grids. The basic idea of the proposed limiting strategy is to control the distribution of both cell-centered and cell-vertex physical properties to mimic multi-dimensional nature of flow physics, which can be formulated to satisfy so called the MLP condition. The MLP condition can guarantee high-order spatial accuracy and improved convergence without yielding spurious oscillations. Starting from the MUSCL-type reconstruction on unstructured grids followed by the efficient implementation of the MLP condition, MLP slope limiters on unstructured meshes are obtained.Thanks to its superior limiting strategy and maximum principle satisfying characteristics, the newly developed MLP on unstructured grids is quite effective in controlling numerical oscillations as well as accurate in capturing multi-dimensional flow features. Numerous test cases are presented to validate the basic features of the proposed approach.     

A Cartesian grid embedded boundary method for hyperbolic conservation laws

We present a second-order Godunov algorithm to solve time-dependent hyperbolic systems of conservation laws on irregular domains. Our approach is based on a formally consistent discretization of the conservation laws on a finite-volume grid obtained from intersecting the domain with a Cartesian grid. We address the small-cell stability problem associated with such methods by hybridizing our conservative discretization with a stable, nonconservative discretization at irregular control volumes, and redistributing the difference in the mass increments to nearby cells in a way that preserves stability and local conservation. The resulting method is second-order accurate in L¹ for smooth problems, and is robust in the presence of large-amplitude discontinuities intersecting the irregular boundary.     

Linear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence

The Osher–Chakrabarthy family of linear flux-modification schemes is considered. Improved lower bounds on the compression factors are provided, which suggest the viability of using the unlimited version. The LLF flux formula is combined with these schemes in order to obtain efficient finite-difference algorithms. The resulting schemes are applied to a battery of numerical tests, going from advection and Burgers equations to Euler and MHD equations, including the double Mach reflection and the Orszag–Tang 2D vortex problem. Total-variation-bounded (TVB) behavior is evident in all cases, even with time-independent upper bounds. The proposed schemes, however, do not deal properly with compound shocks, arising from non-convex fluxes, as shown by Buckley–Leverett test simulations.     

The multi-dimensional limiters for solving hyperbolic conservation laws on unstructured grids

Novel limiters based on the weighted average procedure are developed for finite volume methods solving multi-dimensional hyperbolic conservation laws on unstructured grids. The development of these limiters is inspired by the biased averaging procedure of Choi and Liu [10]. The remarkable features of the present limiters are the new biased functions and the weighted average procedure, which enable the present limiter to capture strong shock waves and achieve excellent convergence for steady state computations. The mechanism of the developed limiters for eliminating spurious oscillations in the vicinity of discontinuities is revealed by studying the asymptotic behavior of the limiters. Numerical experiments for a variety of test cases are presented to demonstrate the superior performance of the proposed limiters.     

Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functions

This paper addresses the convergence properties of implicit numerical solution algorithms for nonlinear hyperbolic transport problems. It is shown that the Newton–Raphson (NR) method converges for any time step size, if the flux function is convex, concave, or linear, which is, in general, the case for CFD problems. In some problems, e.g., multiphase flow in porous media, the nonlinear flux function is S-shaped (not uniformly convex or concave); as a result, a standard NR iteration can diverge for large time steps, even if an implicit discretization scheme is used to solve the nonlinear system of equations. In practice, when such convergence difficulties are encountered, the current time step is cut, previous iterations are discarded, a smaller time step size is tried, and the NR process is repeated. The criteria for time step cutting and selection are usually based on heuristics that limit the allowable change in the solution over a time step and/or NR iteration. Here, we propose a simple modification to the NR iteration scheme for conservation laws with S-shaped flux functions that converges for any time step size. The new scheme allows one to choose the time step size based on accuracy consideration only without worrying about the convergence behavior of the nonlinear solver. The proposed method can be implemented in an existing simulator, e.g., for CO₂ sequestration or reservoir flow modeling, quite easily. The numerical analysis is confirmed with simulation studies using various test cases of nonlinear multiphase transport in porous media. The analysis and numerical experiments demonstrate that the modified scheme allows for the use of arbitrarily large time steps for this class of problems.     

Lattice Boltzmann method for n-dimensional nonlinear hyperbolic conservation laws with the source term

It is important for nonlinear hyperbolic conservation laws (NHCL) to own a simulation scheme with high order accuracy, simple computation, and non-oscillatory character. In this paper, a unified and novel lattice Boltzmann model is presented for solving n-dimensional NHCL with the source term. By introducing the high order source term of explicit lattice Boltzmann method (LBM) and the optimum dimensionless relaxation time varied with the specific issues, the effects of space and time resolutions on the accuracy and stability of the model are investigated for the different problems in one to three dimensions. Both the theoretical analysis and numerical simulation validate that the results by the proposed LBM have second-order accuracy in both space and time, which agree well with the analytical solutions.     

A stable high-order Spectral Difference method for hyperbolic conservation laws on triangular elements

Numerical schemes using piecewise polynomial approximation are very popular for high order discretization of conservation laws. While the most widely used numerical scheme under this paradigm appears to be the Discontinuous Galerkin method, the Spectral Difference scheme has often been found attractive as well, because of its simplicity of formulation and implementation. However, recently it has been shown that the scheme is not linearly stable on triangles. In this paper we present an alternate formulation of the scheme, featuring a new flux interpolation technique using Raviart–Thomas spaces, which proves stable under a similar linear analysis in which the standard scheme failed. We demonstrate viability of the concept by showing linear stability both in the semi-discrete sense and for time stepping schemes of the SSP Runge–Kutta type. Furthermore, we present convergence studies, as well as case studies in compressible flow simulation using the Euler equations.     

The finite volume local evolution Galerkin method for solving the hyperbolic conservation laws

This paper presents a finite volume local evolution Galerkin (FVLEG) scheme for solving the hyperbolic conservation laws. The FVLEG scheme is the simplification of the finite volume evolution Galerkin method (FVEG). In FVEG, a necessary step is to compute the dependent variables at cell interfaces at t_n + τ (0 < τ ? Δt). The FVLEG scheme is constructed by taking τ → 0 in the evolution operators of FVEG. The FVLEG scheme greatly simplifies the evaluation of the numerical fluxes. It is also well suited with the semi-discrete finite volume method, making the flux evaluation being decoupled with the reconstruction procedure while maintaining the genuine multi-dimensional nature of the FVEG methods. The derivation of the FVLEG scheme is presented in detail. The performance of the proposed scheme is studied by solving several test cases. It is shown that FVLEG scheme can obtain very satisfactory numerical results in terms of accuracy and resolution.     

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Luká?ová-Medvid’ová, J. Saibertov’a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533– 562; M. Luká?ová-Medvid’ová, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1–30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor–corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.     

A class of large time step Godunov schemes for hyperbolic conservation laws and applications

A large time step (LTS) Godunov scheme firstly proposed by LeVeque is further developed in the present work and applied to Euler equations. Based on the analysis of the computational performances of LeVeque’s linear approximation on wave interactions, a multi-wave approximation on rarefaction fan is proposed to avoid the occurrences of rarefaction shocks in computations. The developed LTS scheme is validated using 1-D test cases, manifesting high resolution for discontinuities and the capability of maintaining computational stability when large CFL numbers are imposed. The scheme is then extended to multidimensional problems using dimensional splitting technique; the treatment of boundary condition for this multidimensional LTS scheme is also proposed. As for demonstration problems, inviscid flows over NACA0012 airfoil and ONERA M6 wing with given swept angle are simulated using the developed LTS scheme. The numerical results reveal the high resolution nature of the scheme, where the shock can be captured within 1–2 grid points. The resolution of the scheme would improve gradually along with the increasing of CFL number under an upper bound where the solution becomes severely oscillating across the shock. Computational efficiency comparisons show that the developed scheme is capable of reducing the computational time effectively with increasing the time step (CFL number).