利用通量限制的高分辨Godunov型格式的熵条件

1．引言考虑非线性双曲型守恒律方程的Cauchy问题式中f（w）∈C2（R）f＂,（w）≥0,初值。u0∈BV（R）．此问题通常只存在弱解,且需附加熵条件以保证解的唯一性．方程（1.1）的数值方法研究发展很快,但一阶精度格式（如Godunov格式）分辨率很低,而二阶精度格式在间断附近存在振荡;TVD格式则是一种成功的高分辨率无振荡格式．此外,双曲型守恒律数值方法的收敛性取决于差分格式的总变差稳定和离散熵条件．文献[2]中给出了利用通量限制构造TVD格式的方法,[1]则讨论了SOR-TVD格式的熵条件．本文第2节回顾了问的方法,具体导出了…     

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

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

一类交错网格的Gauss型格式

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

交错网格紧致差分格式和满足等价性的压力Poisson方程

1．引言随着电子计算机的发展，越来越多的实际问题数值模拟成为现实，但还有很多非线性问题数值计算时间太长，内存要求过大．数值方法的改进可使计算量和存储量大大减少，例如，对二维非定常问题，要使误差达到N－4量级，二阶格式计算点数为（N2）3，（包括时间方向），而四阶格式计算点数仅为N3，差N3倍！而计算量差的倍数更多．当N＝16时N3＝4096，当N＝256时，N31678万．紧致差分格式具有精度高，差分式基点少，＜线性）稳定性好，对高频波分辨率高，边界差分点少等优点【’，‘，’，’。’，’，‘’］，本文中的格式基点数为3，…     

广义非线性Schroedinger方程的一个新的守恒差分格式

对广义非线性Schroedinger方程提出了一种新的差分格式。揭示了该差分格式满足两个守恒律，并证明该格式的收敛性和稳定性．数值实验结果表明，新的差分格式优于Crank-Nicolson格式以及Zhang Fei等人提出的格式。     

二维非定常Euler方程的守恒迹线格式

双曲型守恒律的计算方法研究,得到了很大的发展,有许多优秀的差分格式.这些格式向多维的推广往往基于维数分步.     

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

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

非线性抛物组非均匀网格差分解的唯一性和稳定性

１．引言 １．对一维非线性抛物组,在文献山中已构造一般非均匀网格差分格式,其中差分逼近的组合系数对不同的网格点和不同的网格层可以不同,并且运用不动点原理证明了差分解的存在性和收敛性．在非均匀网格差分格式中差分逼近的组合系数为常数的情形,文献［２］证明了具有有界二阶差商的离散向量解的存在性、唯一性和稳定性．本文将对文献［１］中构造的一般非均匀网格差分格式,证明所得到的差分解的唯一性和稳定性． 考虑如下非线性抛物组其中是未知的ｍ－维向量函数是给定的矩阵函数,ｊ（ｘ,ｔ,ｕ,ｐ）。是给定的ｍ－维向量函数…     

一维单个守恒律的一类二阶、几乎熵耗散的格式

本文考虑一维单个守恒律方程，对其设计了一个基于熵耗散的非线性守恒型差分格式．本格式的数值流函数是Lax-Freidrichs格式和Lax-Wendroff格式数值流函数的凸组合，凸组合中的系数是由考虑耗散熵来决定的．这样在解的光滑区域内，格式几乎、甚至完全是Lax-Wendroff格式，而在解的间断处，格式几乎、甚至完全是Lax—Freidrichs格式．从而消除了间断附近的非物理振荡，实现了计算的非线性稳定性．理论分析表明本格式在解的非极值点处是二阶精度的，而在解的极值点处至少有一阶精度．数值试验表明格式是有效的．     

一个刚性守恒律方程组的全隐式差分方法

１．引言本文研究如下非线性刚性守恒律方程组的全隐式差分逼近． 方程（1．１）中的源项ｇ（ｕ,ｖ）定义为 ｇ（ｕ,ｖ）＝ｖ－（１－μ）ｆ（ｕ）,（１．２）其中ｆ是ｕ的一个给定函数,δ是一个小正参数,称为松弛时间,μ是参数．方程组（１．１）频繁出现于粘弹性力学中． 在零松弛时间限（δ→０）下,从（１．１）可得到如下方程组该方程组通常称为“平衡”模型,而方程组（１．１）称为“非平衡”模型． 文中将假设μ满足 ０＜ μ＜ １,（１．４）以便保证拟稳定性条件［１９,２０］和次特征条件［１１,２,３］： λ１≤λ＊…     

二维双曲守恒律的大时间步长Godunov方法(英文)

考虑了关于二维守恒律的大时间步长Godunov方法.该方法是关于一维问题的自然推广.证明了文中给出的数值流函数下,该方法是守恒的.进一步还给出了近似Riemann解应满足的条件,并且证明了利用满足这些条件的近似Riemann解的大时间步长Godunov方法守恒.最后,补充证明了满足这些条件的近似Riemann解是满足熵条件的.     

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.     

The third-order relaxation schemes for hyperbolic conservation laws

A class of semi-discrete third-order relaxation schemes are presented for relaxation systems which approximate systems of hyperbolic conservation laws. These schemes for the scalar conservation law are shown to satisfy the property of total variation diminishing (TVD) in the zero relaxation limit. A third-order TVD Runge–Kutta splitting method is developed for the temporal discretization of the semi-discrete schemes. Numerical results are given illustrating these schemes on one-dimensional nonlinear problems.     

ON THE EXPLICIT COMPACT SCHEMES II: EXTENSION OF THE STCE/CE METHOD ON NONSTAGGERED GRIDS ON THE EXPLICIT COMPACT SCHEMES II: EXTENSION OF THE STCE/CE METHOD ON NONSTAGGERED GRIDS

1. IntroductionTills paper is interested in the genuinely nonlinear conserVation lawswith initial data u(0, x) = "o(x), x = (x', ...l x').It is well known that the above problem may not always have a smooth global solutinn even if the initial data no is adequately smooth[6]. Thus, we consider its weaksolutinn so that the Problem (1.1) Ililght have a global solution allowing discontinuities(e.g. shock wave.etc.). Moreover, the elltrOPy conditinn should be deposed inorder to single out a phyS…     

$\delta-$Wave for 1-D and 2-D Hyperbolic Systems

Here a new kind of nonlinear wave, which is called $\delta-$wave, is described by some high resolution difference solutions for Riemann problems of one-dimensional (1-D) and two-dimensional (2-D) nonlinear hyperbolic systems in conservation laws. Some phenomena are numerically shown for the solutions of Riemann problems for 2-D gas dynamics systems.     

Explicit second-order accurate schemes for the nonlinear Schrödinger equations

We consider three-level explicit schemes for solving the nonlinear variable coefficient Schrödinger-type equation. Using spectral and energy methods we establish the stability and convergence of these schemes. The existence of discrete conservation laws is investigated. General results are applied for the DuFort-Frankel and leap-frog diffenrence schemes.     

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…     

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.     

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