MUSTA: A multi-stage numerical flux

Numerical methods for conservation laws constructed in the framework of finite volume and discontinuous Galerkin finite elements require, as the building block, a monotone numerical flux. In this paper we present some preliminary results on the MUSTA approach [E.F. Toro, Multi-stage predictor–corrector fluxes for hyperbolic equations, Technical Report NI03037-NPA, Isaac Newton Institute for Mathematical Sciences, University of Cambridge, UK, 17th June, 2003] for constructing upwind numerical fluxes. The scheme may be interpreted as an un-conventional approximate Riemann solver that has simplicity and generality as its main features. When used in its first-order mode we observe that the scheme achieves the accuracy of the Godunov method used in conjunction with the exact Riemann solver, which is the reference first-order method for hyperbolic systems. At least for the scalar model hyperbolic equation, the Godunov scheme is the best of all first-order monote schemes, it has the smallest truncation error. Extensions of the scheme of this paper are realized in the framework of existing approaches. Here we present a second-order TVD (TVD for the scalar case) extension and show numerical results for the two-dimensional Euler equations on non-Cartesian geometries. The schemes find their best justification when solving very complex systems for which the solution of the Riemann problem, in the classical sense, is too complex, too costly or is simply unavailable.     

守恒型双曲方程的中心差分TVD格式研究

研究了如何利用迎风格式的耗散性构造中心差分ＴＶＤ格式的方法，给 相应的定理，构造出新的耗散表达式。新格式既保留了二阶中心差分格式灵活方便的优点，又吸收了迎风格式耗散项比较精细的特点，同时具有ＴＶＤ性质，使得新格式具有较同的激波分辨率。     

一类交错网格的Gauss型格式

本文在交错网格的情况下 ,利用 Gauss型求积公式构造了一类不需解 Riemann问题的求解一维单个双曲守恒律的二阶显式 Gauss型差分格式 ,证明了该格式在CFL条件限制下为 TVD格式 ,并证明了这类格式的收敛性 ,然后将格式推广到方程组的情形 .由于在交错网格的情况下构造的这类差分格式 ,不需要求解 Riemann问题 ,因此这类格式与诸如 Harten等的 TVD格式相比具有如下优点 :由于不需要完整的特征向量系 ,因此可用于求解弱双曲方程组 ,计算更快、编程更加简便等 .     

Second order scheme for scalar conservation laws with discontinuous flux

Burger, Karlsen, Torres and Towers in [9] proposed a flux TVD (FTVD) second order scheme with Engquist–Osher flux, by using a new nonlocal limiter algorithm for scalar conservation laws with discontinuous flux modeling clarifier thickener units. In this work we show that their idea can be used to construct FTVD second order scheme for general fluxes like Godunov, Engquist–Osher, Lax–Friedrich, … satisfying (A, B)-interface entropy condition for a scalar conservation law with discontinuous flux with proper modification at the interface. Also corresponding convergence analysis is shown. We show further from numerical experiments that solutions obtained from these schemes are comparable with the second order schemes obtained from the minimod limiter.     

A composite semi-conservative scheme for hyperbolic conservation laws

In this work a first order accurate semi-conservative composite scheme is presented for hyperbolic conservation laws. The idea is to consider the non-conservative form of conservation law and utilize the explicit wave propagation direction to construct semi-conservative upwind scheme. This method captures the shock waves exactly with less numerical dissipation but generates unphysical rarefaction shocks in case of expansion waves with sonic points. It shows less dissipative nature of constructed scheme. In order to overcome it, we use the strategy of composite schemes. A very simple criteria based on wave speed direction is given to decide the iterations. The proposed method is applied to a variety of test problems and numerical results show accurate shock capturing and higher resolution for rarefaction fan.     

ON THE CENTRAL RELAXING SCHEMES I:SINGLE CONSERVATION LAWS

1. IntroductionIn [6], Jin and adn constructed a class of uPWind relaxing schemes for nonlinearconservation lawswith initial data u(0, x) ~ "o(x), x ~ (xl, ...t -cd), by using the idea of the local relaxation approximation [2,3,6,10].The relaxing scheme is obtained in the following way: A linear hyperbolic systemwith a stiff source term is first constructed to approximate the original equation (1.1)with a small dissipative correction. Then this linear hyperbolic system is solved easilyby und…     

MULTIDIMENSIONAL RELAXATION APPROXIMATIONS FOR HYPERBOLIC SYSTEMS OF CONSERVATION LAWS

We construct and implement a non-oscillatory relaxation scheme for multidimensional hyperbolic systems of conservation laws. The method transforms the nonlinear hyperbolic system to a semilinear model with a relaxation source term and linear characteristics which can be solved numerically without using either Riemann solver or linear iterations. To discretize the relaxation system we consider a high-resolution reconstruction in space and a TVD Runge-Kutta time integration. Detailed formulation of the scheme is given for problems in three space dimensions and numerical experiments are implemented in both scalar and system cases to show the effectiveness of the method.     

Convergence of relaxation schemes for hyperbolic conservation laws with stiff source terms

We focus in this study on the convergence of a class of relaxation numerical schemes for hyperbolic scalar conservation laws including stiff source terms. Following Jin and Xin, we use as approximation of the scalar conservation law, a semi-linear hyperbolic system with a second stiff source term. This allows us to avoid the use of a Riemann solver in the construction of the numerical schemes. The convergence of the approximate solution toward a weak solution is established in the cases of first and second order accurate MUSCL relaxed methods.

     

Relaxation WENO schemes for multidimensional hyperbolic systems of conservation laws

We present a class of high‐order weighted essentially nonoscillatory (WENO) reconstructions based on relaxation approximation of hyperbolic systems of conservation laws. The main advantage of combining the WENO schemes with relaxation approximation is the fact that the presented schemes avoid solution of the Riemann problems due to the relaxation approach and high‐resolution is obtained by applying the WENO approach. The emphasis is on a fifth‐order scheme and its performance for solving a wide class of systems of conservation laws. To show the effectiveness of these methods, we present numerical results for different test problems on multidimensional hyperbolic systems of conservation laws. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

High-order accurate monotone compact running scheme for multidimensional hyperbolic equations

Monotone absolutely stable conservative difference schemes intended for solving quasilinear multidimensional hyperbolic equations are described. For sufficiently smooth solutions, the schemes are fourth-order accurate in each spatial direction and can be used in a wide range of local Courant numbers. The order of accuracy in time varies from the third for the smooth parts of the solution to the first near discontinuities. This is achieved by choosing special weighting coefficients that depend locally on the solution. The presented schemes are numerically efficient thanks to the simple two-diagonal (or block two-diagonal) structure of the matrix to be inverted. First the schemes are applied to system of nonlinear multidimensional conservation laws. The choice of optimal weighting coefficients for the schemes of variable order of accuracy in time and flux splitting is discussed in detail. The capabilities of the schemes are demonstrated by computing well-known two-dimensional Riemann problems for gasdynamic equations with a complex shock wave structure.     

A high order central-upwind scheme for hyperbolic conservation laws

A high order central-upwind scheme for approximating hyperbolic conservation laws is proposed. This construction is based on the evaluation of the local propagation speeds of the discontinuities and Peer's fourth order non-oscillatory reconstruction. The presented scheme shares the simplicity of central schemes, namely no Riemann solvers are involved. Furthermore, it avoids alternating between two staggered grids, which is particularly a challenge for problems which involve complex geometries and boundary conditions. Numerical experiments demonstrate the high resolution and non-oscillatory properties of our scheme.     

求解双曲型方程的一致二阶守恒型差分格式的构造及数值实验

由于TVD格式具有激波分辨率高与非物理振荡小的特点,在气体动力学问题的求解中得到了广泛的应用.但现有的TVD格式受其构造方式所限,在解的局部极值点附近只能达到一阶精度. 考虑以下单个双曲型方程:     

The relaxation schemes for systems of conservation laws in arbitrary space dimensions

We present a class of numerical schemes (called the relaxation schemes) for systems of conservation laws in several space dimensions. The idea is to use a local relaxation approximation. We construct a linear hyperbolic system with a stiff lower order term that approximates the original system with a small dissipative correction. The new system can be solved by underresolved stable numerical discretizations without using either Riemann solvers spatially or a nonlinear system of algebraic equations solvers temporally. Numerical results for 1-D and 2-D problems are presented. The second-order schemes are shown to be total variation diminishing (TVD) in the zero relaxation limit for scalar equations. ©1995 John Wiley & Sons, Inc.     

Development and Comparison of Numerical Fluxes for LWDG Methods

The discontinuous Galerkin （DO） or local discontinuous Galerkin （LDG） method is a spatial discretization procedure for convection-diffusion equations, which employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes and limiters. The Lax- Wendroff time discretization procedure is an altemative method for time discretization to the popular total variation diminishing （TVD） Runge-Kutta time discretizations. In this paper, we develop fluxes for the method of DG with Lax-Wendroff time discretization procedure （LWDG） based on different numerical fluxes for finite volume or finite difference schemes, including the first-order monotone fluxes such as the Lax-Friedfichs flux, Godunov flux, the Engquist-Osher flux etc. and the second-order TVD fluxes. We systematically investigate the performance of the LWDG methods based on these different numerical fluxes for convection terms with the objective of obtaining better performance by choosing suitable numerical fluxes. The detailed numerical study is mainly performed for the one-dimensional system case, addressing the issues of CPU cost, accuracy, non-oscillatory property, and resolution of discontinuities. Numerical tests are also performed for two dimensional systems.     

An efficient TVD‐WENO method for conservation laws

In this article, we present a high‐resolution hybrid scheme for solving hyperbolic conservation laws in one and two dimensions. In this scheme, we use a cheap fourth order total variation diminishing (TVD) scheme for smooth region and expensive seventh order weighted nonoscillatory (WENO) scheme near discontinuities. To distinguish between the smooth parts and discontinuities, we use an efficient adaptive multiresolution technique. For time integration, we use the third order TVD Runge‐Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Finite-difference methods for one-dimensional hyperbolic conservation laws

This article contains a survey of some important finite-difference methods for one-dimensional hyperbolic conservation laws. Weak solutions of hyperbolic conservation laws are introduced and the concept of entropy stability is discussed. Furthermore, the Riemann problem for hyperbolic conservation laws is solved. An introduction to finite-difference methods is given for which important concepts such as, e.g., conservativity, stability, and consistency are introduced. Godunov-type methods are elaborated for general systems of hyperbolic conservation laws. Finally, flux limiter methods are developed for the scalar nonlinear conservation law. © 1994 John Wiley & Sons, Inc.     

一个求解双曲型守恒律方程的高分辨率GVC格式

1　引　　言本文研究双曲型方程:ut+fx=0,　t>0,-∞相似文献   

An optimal error estimate for upwind Finite Volume methods for nonlinear hyperbolic conservation laws

The purpose of this paper is to show that the cell-centered upwind Finite Volume scheme applied to general hyperbolic systems of m conservation laws approximates smooth solutions to the continuous problem at order one in space and time. As it is now well understood, there is a lack of consistency for order one upwind Finite Volume schemes: the truncation error does not tend to zero as the time step and the grid size tend to zero. Here, following our previous papers on scalar equations, we construct a corrector that allows us to prove the expected error estimate for nonlinear systems of equations in one dimension.     

Error estimate for the upwind finite volume method for the nonlinear scalar conservation law

In this paper we estimate the error of upwind first order finite volume schemes applied to scalar conservation laws. As a first step, we consider standard upwind and flux finite volume scheme discretization of a linear equation with space variable coefficients in conservation form. We prove that, in spite of their lack of consistency, both schemes lead to a first order error estimate. As a final step, we prove a similar estimate for the nonlinear case. Our proofs rely on the notion of geometric corrector, introduced in our previous paper by Bouche et al. (2005) [24] in the context of constant coefficient linear advection equations.     

同位非结构和结构网格摄动有限体积法（PFV）求解二维Navier-Stokes方程组

根据NS方程组的一阶迎风和二阶中心有限体积（UFV和CFV）格式,导出NS方程组迎风和中心摄动有限体积(UPFV和CPFV)格式．该格式通过把控制体界面质量通量摄动展开成网格间距的幂级数,并由守恒方程本身求得幂级数系数而获得．迎风和中心摄动有限体积格式使用了与一阶迎风和二阶中心格式相同的基点数和相同的表达形式,宜于计算机编程．顶盖驱动方腔流和驻点流标量输运的数值实验证明,迎风PFV格式比一阶UFV、二阶CFV格式有更高的精度,更高的分辨率．尤其是良好的鲁棒特性．对顶盖驱动方腔流,在Re数从102到104范围内,亚松弛系数可在0.3～0.8任取,收敛性能良好．