A FAMILY OF HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEMES WITH BRANCHING STABILITY FOR SOLVING 3-D PARABOLIC PARTIAL DIFFERENTIAL EQUATION

A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and the truncation error is 0(<Delta>t(2) + Deltax(4)).

High Order of Accuracy for Poisson Equation Obtained by Grouping of Repeated Richardson Extrapolation with Fourth Order Schemes

In this article,we improve the order of precision of the two-dimensional Poisson equation by combining extrapolation techniques with high order schemes.The high order solutions obtained traditionally generate non-sparse matrices and the calculation time is very high.We can obtain sparse matrices by applying compact schemes.In this article,we compare compact and exponential finite difference schemes of fourth order.The numerical solutions are calculated in quadruple precision(Real*16 or extended precision)in FORTRAN language,and iteratively obtained until reaching the round-off error magnitude around 1.0E−32.This procedure is performed to ensure that there is no iteration error.The Repeated Richardson Extrapolation(RRE)method combines numerical solutions in different grids,determining higher orders of accuracy.The main contribution of this work is based on a process that initializes with fourth order solutions combining with RRE in order to find solutions of sixth,eighth,and tenth order of precision.The multigrid Full Approximation Scheme(FAS)is also applied to accelerate the convergence and obtain the numerical solutions on the fine grids.

A－HIGH－ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THE EQUATION OF TWO－DIMENSIONAL PARABOLIC TYPE

In this paper. a three explicit difference shcemes with high order accuracy for solving the equations of two-dimensional parabolic type is proposed. The stability condition is r=△t/△x￣ 2=△t/△y￣2≤1/4 and the truncation error is O (△t￣2 + △x￣4 ).

A CLASS OF COMPACT UPWIND TVD DIFFERENCE SCHEMES

A new method was proposed for constructing total variation diminishing （TVD） upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD （CTVD） schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.

A NEW HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THREE-DIMENSIONAL PARABOLIC EQUATIONS

In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).

线性双曲型方程新的TVD格式

从双曲型方程的ＴＶＤ条件出发，分析了ＬａｘＷｅｎｄｒｏｆ格式和ＷａｒｍｉｎｇＢｅａｍ格式所存在的ＴＶＤ区间，构造了空间三点和四点新的二阶ＴＶＤ格式，使其在极值点也保持二阶精度。

双曲型线性方程三阶和四阶TVD格式的新构造

利用 Taylor级数理论和总变差减小 ( TVD)格式的充分条件构造了时间二阶、空间五点三阶和四阶新 TVD格式 .给出了新 TVD格式与传统 TVD格式及近期建立的二阶新 TVD格式用于线性双曲型方程的计算结果 ,表明本文新格式特别是四阶 TVD格式具有比二阶新 TVD格式和传统 TVD格式峰值衰减更慢、间断更陡 ,而计算工作量具有与传统二阶

双曲型守恒律的一种高精度TVD差分格式

构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例.

高精度两步TVD格式的构造及数值检验

基于通量分裂的思想,利用Taylor级数理论和第二步限制器函数的控制,构造了在光滑区域空间为三阶、时间为一阶精度的高分辨率全变差递减(TVD)格式。该格式在激波过渡点降为一阶迎风格式。并从单个线形方程推广到非线性方程及方程组情形。通过几个典型算例的计算,并与二阶TVD格式作了比较,表明该方法对流场中的激波有较高的分辨率,且是无波动的。

一类高精度TVD差分格式及其应用

构造了一维非线性双曲型守恒律的一个新的高精度、高分辨率的守恒型TVD差分格式。其构造思想是:首先,将计算区间划分为若干个互不相交的小区间,再根据精度要求等分小区间,通过各细小区间上的单元平均状态变量,重构各细小区间交界面上的状态变量,并加以校正;其次,利用近似R iemann解计算细小区间交界面上的数值通量,并结合高阶Runge-Kutta TVD方法进行时间离散,得到了高精度的全离散方法。证明了该格式的TVD特性。该格式适合于使用分量形式计算而无须进行局部特征分解。通过计算几个典型的问题,验证了格式具有高精度、高分辨率且计算简单的优点。

一种二阶TVD差分格式构造方法

采用一阶迎风差分格式,作Taylor展开,消去低阶项,给出了求解一维和二维等熵气动力学方程组的一种双参数二阶精度TVD差分格式的构造方法.

一类基于通量分裂的二阶精度TVD差分格式

基于通量分裂和逆风特性,选取单元交界面上的正、负数值通量,并结合高阶Runge-KuttaTVD时间离散方法,构造了一维非线性双曲型守恒律的一类二阶精度的差分格式,证明了差分格式的TVD特性。按分量形式推广到方程组。通过几个典型的数值算例验证了格式的有效性。

High-Order Accurate Entropy Stable Finite Difference Schemes for One- and Two-Dimensional Special Relativistic Hydrodynamics

This paper develops the high-order accurate entropy stable finite difference schemes for one-and two-dimensional special relativistic hydrodynamic equations.The schemes are built on the entropy conservative flux and the weighted essentially non-oscillatory(WENO)technique as well as explicit Runge-Kutta time discretization.The key is to technically construct the affordable entropy conservative flux of the semi-discrete second-order accurate entropy conservative schemes satisfying the semi-discrete entropy equality for the found convex entropy pair.As soon as the entropy conservative flux is derived,the dissipation term can be added to give the semidiscrete entropy stable schemes satisfying the semi-discrete entropy inequality with the given convex entropy function.The WENO reconstruction for the scaled entropy variables and the high-order explicit Runge-Kutta time discretization are implemented to obtain the fully-discrete high-order entropy stable schemes.Several numerical tests are conducted to validate the accuracy and the ability to capture discontinuities of our entropy stable schemes.

双曲方程的一种二阶TVD差分格式构造方法

本文采用一阶迎风差分格式作Ｔａｙｌｏｒ展开，消去低阶项，给出求解守恒型双曲方程初值问题的一种二阶ＴＶＤ 差分格式的构造方法。

一类椭圆型Dirichlet边值问题的高精度Richardson外推法

针对椭圆型偏微分方程,先建立四阶和六阶精度的紧致差分格式,在此基础上用Richardson外推法,得到其六阶和八阶精度的外推差分格式。并通过两个Poisson方程算例,验算已建立的差分格式。数值算例结果表明,基于紧致差分格式的Richardson外推法能够得到有效的、健壮的高精度数值解。

High Order Finite Difference WENO Methods with Unequal-Sized Sub-Stencils for the Degasperis-Procesi Type Equations

In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,which contain nonlinear high order derivatives,and possibly peakon solutions or shock waves.By introducing auxiliary variable(s),we rewrite the DP equation as a hyperbolic-elliptic system,and theµDP equation as afirst order system.Then we choose a linearfinite difference scheme with suitable order of accuracy for the auxiliary variable(s),and twofinite difference WENO schemes with unequal-sized sub-stencils for the primal variable.One WENO scheme uses one large stencil and several smaller stencils,and the other WENO scheme is based on the multi-resolution framework which uses a se-ries of unequal-sized hierarchical central stencils.Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil,both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights.Another advantage is that thefinal reconstructed polynomial on the target cell is a polynomial of the same de-gree as the polynomial over the big stencil,while the classicalfinite difference WENO reconstruction can only be obtained for specific points inside the target interval.Nu-merical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构，通常采用高阶精度的差分格式进行数值模拟。为能抑制计算过程中和在激波附近产生的虚假波动，本文从物理构思出发，提出了建立高阶格式的两个基本原则，作者称之为抑制波动的原则和稳定性原则。通过对典型问题的计算表明，这些原则是有效的，由此建立的二阶格式即ＮＮＤ格式无虚假波动，由此建立的三、四阶以上的格式仅在激波附近有微小波动。算例还表明，不符合文中原则的高阶格式的结果则要差得多。为了进一步克服激波附近的微小波动，文中还建议了两种方法，一种是基于熵增原则的方法，另一种是采用激波附近为二阶ＮＮＤ格式、其它区域为高阶格式的混合方法。

解抛物型方程的一个高精度显式差分格式

引入耗散项的方法,构造一个条件稳定的显格式,其稳定性条件为r≤1/2,截断误差可达到O(τ2+h4+hτ22).当τ=O(h)时,此格式可逼近精度,特别当τ=O(h2)时,格式达到二阶精度.数值例子表明,所建立的差分格式是有效的.

双曲型守恒律的一类局部化的高效差分格式

构造了一维非线性双曲型守恒律的一类局部化的高效全离散差分格式,并将该格式推广到一维守恒方程组及二维守恒方程(组).最后,给出了几个标准算例.数值计算结果表明此格式具有高精度高分辨激波、稀疏波和接触间断,且边界条件易于处理等优点.

一个基于通量分裂的高精度MmB差分格式

本文研究双曲型守恒律的高精度差分方法。在新的算法中 ,首先将计算区间划分为互不相交的小区间 ,再根据精度要求等分小区间 ;其次 ,根据流动方向进行通量分裂 ,重构小区间交界面上的正、负数值通量 ,并进行校正 ;然后 ,采用高阶Runge KuttaTVD方法进行时间离散 ,构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的守恒型差分格式。推广到二维双曲型守恒律方程 ,证明了格式的MmB特性。进而推广到二维守恒型方程组情形。最后对二维Burgers方程及Euler方程进行了数值试验 。