一种新的求解双曲型守恒律的三阶中心差分格式

本文提出了一种求解双曲型守恒律新的三阶中心差分格式，主要是引入了一种推广的三阶重构，并证明了这种重构在网格边界无振荡．所提的格式保持了中心差分格式简单的优点，不需用Riemann解算器，避免了进行特征解耦．数值试验结果表明本文格式是高精度、高分辨率的。     

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

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

守恒型双曲型方程组的差分方法(Ⅱ)

本文是[1]的继续,将介绍守恒型双曲型方程组的各种其他差分方法,例如基于Riemann间断分解的 格式,Glimm格式和Chorin的随机选取法,人工粘性法,人工压缩法,特征型格式和质点法等。本文所采用的记号同[1]。 本文继续介绍下列守恒型双曲型方程组的差分方法     

具有初始跳跃的双曲型方程奇异摄动问题的数值解法

本文考虑具有初始跳跃的二阶双曲型方程初边值问题.首先给出解的导数估计.然后在一非均匀网格上构造了一个差分格式,最后在能量范数意义下证明了差分格式解的一致收敛性.     

非定常自由面流激波解的二阶守恒算法

将计算双曲型守恒律弱解的Lax-Wendroff型TVD格式推广到断面形状沿程任意变化的一般浅水方程组，构造了二阶精度的差分格式．新格式适用于模拟天然河道中溃坝洪水波的传播．提供了表明方法性能的算例，实际天然梯级水库溃坝问题的数值实验表明格式稳定，适应性强．     

一阶双曲型方程组初边值问题的差分格式及其稳定性

本文比较系统地讨论了有关数值求解两个自变量的一阶双曲型方程组初边值问题的某些问题,给出了几种能用于任何类型的初边值问题的差分格式,并在很宽的条件下证明了其中的某些变系数的初边值问题的差分格式对初值和边值是稳定的、差分格式所立出的方程组是良态的.其中的某些格式已用于解决某些复杂的实际问题(应用部分见[16]).     

流体力学方程组的总熵增量小的守恒型差分格式(续)

近年来,国外许多学者对求解双曲守恒律组的高分辨率、高精度差分格式进行了深入的研究.例如MUSCL方法、TVD格式、PPM方法、各种限流的方法以及ENO格式等等.将这些方法应用于流体力学方程组,其数值实践的结果表明,在消除波后振荡、提高激波间断分辨率、提高计算精度等方面有明显的效果.在设计这些计算格式时,通常都是研究单个标量方程的计算格式,再推广到方程组的情形.同时,或者对数值解的总变差提出某种要求(不增或基本不增),或者采用修正数值流措施,或者采用插值或重构的方法,在网格内部用线性分布和更高阶的分布取代Godunov方法中的常数分布,以及处理相应的小范围的解的算法.     

进化型方程的差分格式(Ⅰ)

本文对一些常系数进化型方程证明了除常微分方程和一阶双曲型方程特殊的差分格式之外不存在绝对相容、绝对稳定的显式格式.作为推论,我们指出Hadjidimos关于他的格式是绝对相容的结论是错误的,以及对于一般常微分方程组的显式格式,能用的步长一定要满足条件△t=o(R~(-0.5),其中R为方程组右函数的.Jacobi阵的谱半径.对于热传导方程,本文利用差分算子的分裂技巧构造了两类可以显式求解的格式.它们是绝对     

流体力学方程组的总熵增量小的守恒型差分格式

1.引言 近年来,国外许多学者对求解双曲守恒律组的高分辨率、高精度差分格式进行了深入的研究。例如MUSCL方法、TVD格式、PPM方法、各种限流的方法以及ENO格式等等。将这些方法应用于流体力学方程组,其数值实践的结果表明,在消除波后振荡、提高激波间断分辨率、提高计算精度等方面有明显的效果。在设计这些     

一阶拟线性方程差分解法的误差估计

利用网格单元精确解结合守恒积分而导出差分格式这一途径,对于一阶拟线性方程和一阶拟线性双曲型方程组初始值问题有着理论意义和现实意义。早在五十年代,著名的Lax格式,格式,格式等实际上都可以通过网格单元精确解结合守恒积分而导出。本文企图通过这一离散化途径推导出一阶拟线性方程初值问题的差分格式,并讨论此差分格式的误差估计。     

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.     

Central ADER schemes for hyperbolic conservation laws

In this paper we first briefly review the very high order ADER methods for solving hyperbolic conservation laws. ADER methods use high order polynomial reconstruction of the solution and upwind fluxes as the building block. They use a first order upwind Godunov and the upwind second order weighted average (WAF) fluxes. As well known the upwind methods are more accurate than central schemes. However, the superior accuracy of the ADER upwind schemes comes at a cost, one must solve exactly or approximately the Riemann problems (RP). Conventional Riemann solvers are usually complex and are not available for many hyperbolic problems of practical interest. In this paper we propose to use two central fluxes, instead of upwind fluxes, as the building block in ADER scheme. These are the monotone first order Lax-Friedrich (LXF) and the third order TVD flux. The resulting schemes are called central ADER schemes. Accuracy of the new schemes is established. Numerical implementations of the new schemes are carried out on the scalar conservation laws with a linear flux, nonlinear convex flux and non-convex flux. The results demonstrate that the proposed scheme, with LXF flux, is comparable to those using first and second order upwind fluxes while the scheme, with third order TVD flux, is superior to those using upwind fluxes. When compared with the state of art ADER schemes, our central ADER schemes are faster, more accurate, Riemann solver free, very simple to implement and need less computer memory. A way to extend these schemes to general systems of nonlinear hyperbolic conservation laws in one and two dimensions is presented.     

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.

     

以中心- 迎风数值格式研究具有外力项的p 系统

本文以半离散中心- 迎风数值格式研究具有外力项的p 系统. 中心型数值格式用来处理双曲型守恒律或系统的优势是快速且简单, 因为不需要使用近似Riemann 解, 也不需要做特征分解. 我们的数值模拟验证了理论研究结果: 具有外力项的p 系统的解的收敛及爆破行为, 同时也指出一些尚待理论研究的问题.     

Sonic and kinetic phase transitions with applications to Chapman–Jouguet deflagrations

We consider an n × n system of hyperbolic conservation laws and focus on the case of strongly underdetermined sonic phase boundaries. We propose a Riemann solver that singles out solutions uniquely. This Riemann solver has two features: it selects phase boundaries by means of an exterior function and it allows compound waves. Then we prove the global existence of weak solutions to the Cauchy problem. Applications to Chapman–Jouguet deflagrations are given. Copyright © 2004 John Wiley & Sons, Ltd.     

ON THE CENTRAL RELAXING SCHEME Ⅱ: SYSTEMS OF HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are interested in construction of the central reltalng sChemes for system of noIilinearhyperbolic conservation lawswith initial data U(0, x) = Uo(x), x = (x1 ? ...! xd), based on the local relaJxation approkimationof Eq.(1.1) [2, 3, 6, 8, 9, 12].To i11ustrate the basic idea of the relaalng schemes, for the sake of simplicity in the presentation, we restrict our attention to onedimensional scalar conservaioll lawsFirst, introduce a linear hyperbollc system with a stiff sourc…     

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.     

Hybrid WENO schemes with different indicators on curvilinear grids

In {J. Comput. Phys. 229 (2010) 8105-8129}, we studied hybrid weighted essentially non-oscillatory (WENO) schemes with different indicators for hyperbolic conservation laws on uniform grids for Cartesian domains. In this paper, we extend the schemes to solve two-dimensional systems of hyperbolic conservation laws on curvilinear grids for non-Cartesian domains. Our goal is to obtain similar advantageous properties as those of the hybrid WENO schemes on uniform grids for Cartesian domains. Extensive numerical results strongly support that the hybrid WENO schemes with discontinuity indicators on curvilinear grids can also save considerably on computational cost in contrast to the pure WENO schemes. They also maintain the essentially non-oscillatory property for general solutions with discontinuities and keep the sharp shock transition.     

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     

Kinetic approximation of a boundary value problem for conservation laws

Summary. We design numerical schemes for systems of conservation laws with boundary conditions. These schemes are based on relaxation approximations taking the form of discrete BGK models with kinetic boundary conditions. The resulting schemes are Riemann solver free and easily extendable to higher order in time or in space. For scalar equations convergence is proved. We show numerical examples, including solutions of Euler equations.Mathematics Subject Classification (2000): 65M06, 65M12, 76M20Correspondence to: D. Aregba-Driollet