关于“交错网格紧致差分格式和满足等价性的压力Poisson方程”一文的两点说明

“交错网格紧致差分格式和满足等价性的压力Poisson方程”（19：1（1997），83－90）作者注：1．四阶格式要求边条件至少三阶精度．而（2．13）只有二阶精度，要得到三阶，我们可以象（2．14）那样，在（2．13）左端加2．为了保证四阶Runge－Kutta方法（对非定常边条件）的精度，我们将在下一文中用如下中间层边条件，其中见（2．17）．##F56关于“交错网格紧致差分格式和满足等价性的压力Poisson方程”一文的两点说明@于欣$中国科学院力学研究所     

RLW-KdV方程的紧致有限差分格式

本文研究了RLW-KdV方程的一个三层线性紧致有限差分格式.该格式是质量守恒和能量守恒的,用离散能量法证明了差分格式的收敛性和稳定性.所建格式的收敛阶为O(τ~2+h~4).数值实验验证了该格式的有效性和可靠性.     

一维电磁波方程的四阶紧致差分格式

本文建立求解一维电磁波方程的四阶紧致差分格式,运用von Neumann法给出方法的稳定条件.运用能量法证明格式的收敛性.最后,数值例子验证了格式的有效性.     

一种求解Burgers方程的高精度紧致差分格式

针对Burgers方程,采用余项修正法和欧拉公式,推导了一种新的四层高精度紧致差分隐格式,其截断误差为O(τ~2+τh~2+h~4),即当τ=O(h~2)时,格式空间具有四阶精度;然后通过数值实验验证了格式的精确性和可靠性.     

线性Boussinesq方程的四阶紧致差分格式

基于紧致差分方法,推导出一个时间和空间方向均为四阶精度的三层隐式紧致格式,并采用Fourier分析法给出了格式的稳定性条件.最后,数值例子验证了所给出来的格式的精度和可靠性.     

非定常不可压Navier-Stokes方程的高效和稳健的差分格式Ⅱ

A second order accurate implicit finite difference scheme CNMT2 is proposed in this paper for the unsteady incompressible Navier-Stokes equations. It is proved that the scheme is unconditionally nonlinearly stable on smoothly nonuniform halfstaggered meshes; this stability also holds for this scheme with the pressure correction projection method. However, it is found that the pressure correction projection method may lead to deviation problems in practical simulation of high Re flow;the reason and the cure is given in this paper in terms of differential-algebraic equations.     

求解对流扩散方程的紧致pade'逼近差分格式

将指数变换u(x,t)=p(x,t)e~(k/(2ε)x),p(x,t)=v(x,t)e~(st)、pade'逼近与紧致差分方法相结合,对线性对流扩散问题提出了精度为o(τ~4+h~4)的差分格式,分析了稳定性.最后通过数值算例说明格式的有效性.     

RLW-KdV方程的高阶紧致有限差分格式（英文）

对RLW-KdV方程提出一种新的四阶精度紧致有限差分格式.用离散能量法证明差分格式的能量守恒性、可解性、收敛性和稳定性.在离散L~∞-范数下,所建格式在空间上四阶收敛且在时间上二阶收敛.通过两个数值算例验证了该格式的有效性和可靠性.     

发展型Navier-Stokes方程基于压力稳定化方法的谱逼近格式

本文采用压力稳定化方法近似模拟不可压缩条件,进而构造了发展型非周期NavierStokes方程的全离散Legendre谱逼近计算格式,严格分析了格式的广义稳定性与收敛性．本文建立的逼近结果也适用其它非周期问题．     

解四阶抛物型方程高精度紧致差分格式

针对四阶抛物型方程周期初值问题,提出了一个两层隐式差分格式和一个三层隐式差分格式.它们的局部截断误差分别为O((Δt)2+(Δx)4)和O((Δt)2+(Δt)(Δx)2+(Δx)4),其中Δt,Δx分别为时间步长和空间步长.误差分析和数值实验均表明,本文构造的差分格式比经典的Crank-Nicolson格式和Saul’ev构造的差分格式精度更高.从精度及稳定性方面考虑,本文构造的格式也比文[5]的显式格式要好.     

气动计算中的紧致格式与迎风紧致格式

§1.引言 巨型计算机的发展为解决流体力学问题提供了强有力的手段.由于我们所面临的问题日趋复杂,在今天如何提高数值模拟精度和提高求解效率仍是计算流体力学工作者面临的重要课题. 自Beam和Warming提出隐式求解法以来,NS方程的求解效率大为提高.此后人们在这方面作了很多工作,如在[2]中利用线Gauss-Seidel迭代法加快了迭代的收     

对流扩散方程的四阶紧凑迎风差分格式

§1.引言 流动和传热传质的基本方程均是对流扩散型的.对流扩散方程的高阶紧凑差分格式,作为提高计算可靠性和节省计算量的一条有效途径,已引起相当的重视.作为该领域的一大进展,新近由Dennis推出的对流扩散方程四阶紧凑格式,在二维情形下呈九点式且勿须引入中间变量,只涉及对流扩散量本身,能在较粗网格下获取较为准确的数值结果.从本质上说,该格式系指数型四阶紧凑格式的多项式型翻版.它与指数型紧凑格     

一个新的差分格式的稳定性和收敛性

Douglas给出了三种差分格式并讨论其收敛性和稳定性.[2]中给出了四种既稳定又收敛的差分格式,本文则构造了一个高精度的差分格式,并讨论该格式的收敛性和稳定性.     

OPTIMAL POINT-WISE ERROR ESTIMATE OF A COMPACT FINITE DIFFERENCE SCHEME FOR THE COUPLED NONLINEAR SCHRODINGER EQUATIONS

In this paper, we analyze a compact finite difference scheme for computing a coupled nonlinear SchrSdinger equation. The proposed scheme not only conserves the totM mass and energy in the discrete level but also is decoupled and linearized in practical computa- tion. Due to the difficulty caused by compact difference on the nonlinear term, it is very hard to obtain the optimal error estimate without any restriction on the grid ratio. In order to overcome the difficulty, we transform the compact difference scheme into a special and equivalent vector form, then use the energy method and some important lemmas to obtain the optimal convergent rate, without any restriction on the grid ratio, at the order of O（h4 ＋r2） in the discrete L∞ -norm with time step - and mesh size h. Finally, numerical results are reported to test our theoretical results of the proposed scheme.     

A COMPACT UPWIND SECOND ORDER SCHEME FOR THE EIKONAL EQUATION

We present a compact upwind second order scheme for computing the viscosity solution of the Eikonal equation. This new scheme is based on： 1. the numerical observation that classical first order monotone upwind schemes for the Eikonal equation yield numerical upwind gradient which is also first order accurate up to singularities; 2. a remark that partial information on the second derivatives of the solution is known and given in the structure of the Eikonal equation and can be used to reduce the size of the stencil. We implement the second order scheme as a correction to the well known sweeping method but it should be applicable to any first order monotone upwind scheme. Care is needed to choose the appropriate stencils to avoid instabilities.     

Camassa-Holm方程的守恒有限差分格式

1引言Camassa-Holm(C-H)方程是一类十分重要而又特别的新型浅水波方程.1981年,C-H方程由Fuchssteniner和Fokas作为具有双Hamilton结构的例子给出,随后在1993年,Comassa和Holm将其作为浅水波方程重新提出[1],发现了其具有的一些特殊性质—尖峰孤波解和blow-up解等,由此引发了人们对C-H方程的极大兴趣.关于其解的各种性质已有许多工作:1998年,Constantin研究了C-H方程周期整体解的存在性,谱与逆谱问题,     

A FIFTH-ORDER ACCURATE WEIGHTED ENN DIFFERENCE SCHEME AND ITS APPLICATIONS

1. IntroductionIn the paPer [1], where Zhang Hed et al. preseated the nonoocUlatory 3rdorder ENNWence scheme. The idea of ENN scheme is to compare the 1st-order dmerence and 2ndorderWence to attain 3rdorder accurate scheme and to avod spurious oscillatious neax shocks.However the ENN stwe has certain drawbaCks. One Problem is ouly 3rdords accuracyeven in the very smooth regions. AIlOther is to use a 1ot of logical statemeats which dst theconvergence rate and the efficiency Of parallel…     

A MIXED FINITE ELEMENT METHOD ON A STAGGERED MESH FOR NAVIER-STOKES EQUATIONS

In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are defined on three different meshes. This method is a conforming quadrilateral Q1 × Q1 - P0 element approximation for the Navier-Stokes equations. First-order error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to illustrate the effectiveness of the proposed method.     

反应扩散方程的紧交替方向差分格式

本文研究二维常系数反应扩散方程的紧交替方向隐式差分格式．首先综合应用降阶法和降维法导出了紧差分格式,并给出了差分格式截断误差的表达式．其次引进过渡层变量,给出了紧交替方向隐式差分格式算法．接着用能量分析方法给出了紧交替方向隐式差分格式的解在离散H^1范数下的先验估计式,证明了差分格式的可解性、稳定性和收敛性,在离散H^1范数下收敛阶为O(r^2 H^4)．然后将Rechardson外推法应用于紧交替方向隐式差分格式,外推一次得到具有O(r^4 H^6)阶精度的近似解．最后给出了数值例子,数值结果和理论结果是吻合的．