求解一维定常对流扩散反应方程的混合型紧致差分格式

借助显式紧致格式和隐式紧致格式的思想,基于截断误差余项修正,并结合原方程本身,构造出了一种求解一维定常对流扩散反应方程的高精度混合型紧致差分格式.格式仅用到三个点上的未知函数值及一阶导数值,而一阶导数值利用四阶Pade格式进行计算,格式整体具有四阶精度.数值实验结果验证了格式的精确性和可靠性.     

二维系数不连续Helmholtz方程的高阶紧致差分格式

针对二维系数不连续Helmholtz方程,提出和研究了高阶紧致差分格式,在波数跳跃位置引入局部网格加密技巧进行网格加密.数值实验验证,该高阶紧致差分格式用于求解二维系数不连续Helmholtz方程可以达到四阶精度,局部网格加密技巧能够有效地提高数值解的精度.     

对论文"任意精度的三点紧致显格式及其在CFD中的应用"的评论

证明了林建国等(林建国,谢志华,周俊陶,任意精度的三点紧致显格式及其在CFD中的应用.应用数学和力学,2007,28(7):843-852)提出的紧致显格式与传统的差分格式实质相同,是传统差分格式的另一表达形式,并不具有紧致格式的优点.尽管如此,但这种表达形式更紧凑,推导获得高精度的差分表达式相对于传统的Taylor展开求待定系数的方法也更加简单.     

具有三阶精度的数值微分紧致格式及其应用

本文从微分代数精度概念出发 ,引入了数值微分的紧致性概念 ,并且构造了在一般意义下的具有三点三阶精度的数值微分格式 ,以及在等距节点这种特殊情况下的计算格式 .最后 ,通过数值实验证实了格式的精度 .     

高阶精度耗散加权紧致非线性格式

基于构造耗散紧致线性格式的方法,发展了新的高阶精度耗散加权紧致非线性格式(DWCNS). 通过Fourier分析方法讨论了DWCNS的耗散与色散特性. 从修正波数来看,DWCNS 在光滑区的精度与显式迎风偏斜的5阶精度格式接近. 发展了边界格式和靠近边界格式,分析了均匀网格和拉伸网格上的渐近稳定性,并讨论了在多维Euler方程和Navier-Stokes方程中的应用. 几个典型的无黏/有黏算例表明DWCNS对间断有很好的捕捉能力,对边界层的分辨精度也很高,并具有很好的收敛性（含有强激波流场计算,其均方根残差可以收敛到机器零）,尤其是在网格设计合理的情况下能抑制TVD和ENO格式均无法克服的涡量场的非物理波动.     

任意精度的三点紧致显格式及其在CFD中的应用

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间任意精度的三点紧致显格式的表达式,又由Fourier分析法得到了格式的数值弥散和耗散特性．与以往的高精度紧致差分格式不同,提出的格式不用隐式求解代数方程组并且可以达到任意精度．通过方波问题和顶盖方腔流的算例表明,格式在稀疏网格下可以得到很高的精度,不仅能节省计算量,而且易于编程,有很高的计算效率．     

二维分数阶发展型方程的正式的二阶BDF交替方向隐式紧致差分格式

该文将研究二维分数阶发展型方程的正式的二阶向后微分公式(BDF)的交替方向隐式(ADI)紧致差分格式.在时间方向上用二阶向后微分公式离散一阶时间导数,积分项用二阶卷积求积公式近似,在空间方向上用四阶精度的紧致差分离散二阶空间导数得到全离散紧致差分格式.基于与卷积求积相对应的实二次型的非负性,利用能量方法研究了差分格式的稳定性和收敛性,理论结果表明紧致差分格式的收敛阶为O(k~(a+1)+h_1~4+h_2~4),其中k为时间步长,h_1和h_2分别是空间x和y方向的步长.最后,数值算例验证了理论分析的正确性.     

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

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

A High-Order Finite Difference Scheme for 3D Unsteady Convection Diffusion Reaction Equations北大核心CSCD

针对三维非稳态对流扩散反应方程,构造了一种高精度紧致有限差分格式,对空间的离散采用四阶紧致差分方法,对时间的离散采用Taylor级数展开和余项修正技术,所提格式在时间上的精度为二阶、在空间上的精度为四阶。利用Fourier稳定性分析法证明了该格式是无条件稳定的。最后给出数值算例验证了理论结果。     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.     

A fourth-order method of the convection-diffusion equations with Neumann boundary conditions

In this paper, we have developed a fourth-order compact finite difference scheme for solving the convection-diffusion equation with Neumann boundary conditions. Firstly, we apply the compact finite difference scheme of fourth-order to discrete spatial derivatives at the interior points. Then, we present a new compact finite difference scheme for the boundary points, which is also fourth-order accurate. Finally, we use a Padé approximation method for the resulting linear system of ordinary differential equations. The presented scheme has fifth-order accuracy in the time direction and fourth-order accuracy in the space direction. It is shown through analysis that the scheme is unconditionally stable. Numerical results show that the compact finite difference scheme gives an efficient method for solving the convection-diffusion equations with Neumann boundary conditions.     

Iterative methods and high-order difference schemes for 2D elliptic problems with mixed derivative

We derive a fourth-order compact finite difference scheme for a two-dimensional elliptic problem with a mixed derivative and constant coefficients. We conduct experimental study on numerical solution of the problem discretized by the present compact scheme and the traditional second-order central difference scheme. We study the computed accuracy achieved by each scheme and the performance of the Gauss-Seidel iterative method, the preconditioned GMRES iterative method, and the multigrid method, for solving linear systems arising from the difference schemes.     

Compact difference scheme for two-dimensional fourth-order nonlinear hyperbolic equation

High-order compact finite difference method for solving the two-dimensional fourth-order nonlinear hyperbolic equation is considered in this article. In order to design an implicit compact finite difference scheme, the fourth-order equation is written as a system of two second-order equations by introducing the second-order spatial derivative as a new variable. The second-order spatial derivatives are approximated by the compact finite difference operators to obtain a fourth-order convergence. As well as, the second-order time derivative is approximated by the central difference method. Then, existence and uniqueness of numerical solution is given. The stability and convergence of the compact finite difference scheme are proved by the energy method. Numerical results are provided to verify the accuracy and efficiency of this scheme.     

一维Allen-Cahn方程紧差分格式的离散最大化原则和能量稳定性研究

本文主要研究相场模拟中的Allen-Cahn模型,考虑一维Allen-Cahn方程紧差分方法的数值逼近.建立具有O(∫τ2+h4)精度的全离散紧差分格式,证明在合理的步长比和时间步长的约束下,其数值解满足离散最大化原则,在此基础上,研究了全离散格式的能量稳定性.最后给出数值算例.     

A high‐order ADI finite difference scheme for a 3D reaction‐diffusion equation with neumann boundary condition

In this article, we extend the fourth‐order compact boundary scheme in Liao et al. (Numer Methods Partial Differential Equations 18 (2002), 340–354) to a 3D problem and then combine it with the fourth‐order compact alternating direction implicit (ADI) method in Gu et al. (J Comput Appl Math 155 (2003), 1–17) to solve the 3D reaction‐diffusion equation with Neumann boundary condition. First, the reaction‐diffusion equation is solved with a compact fourth‐order finite difference method based on the Padé approximation, which is then combined with the ADI method and a fourth‐order compact scheme to approximate the Neumann boundary condition, to obtain fourth order accuracy in space. The accuracy in the temporal dimension is improved to fourth order by applying the Richardson extrapolation technique, although the unconditional stability of the numerical method is proved, and several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed new algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Integrated fast and high‐accuracy computation of convection diffusion equations using multiscale multigrid method

We present an explicit sixth‐order compact finite difference scheme for fast high‐accuracy numerical solutions of the two‐dimensional convection diffusion equation with variable coefficients. The sixth‐order scheme is based on the well‐known fourth‐order compact (FOC) scheme, the Richardson extrapolation technique, and an operator interpolation scheme. For a particular implementation, we use multiscale multigrid method to compute the fourth‐order solutions on both the coarse grid and the fine grid. Then, an operator interpolation scheme combined with the Richardson extrapolation technique is used to compute a sixth‐order accurate fine grid solution. We compare the computed accuracy and the implementation cost of the new scheme with the standard nine‐point FOC scheme and Sun–Zhang's sixth‐order method. Two convection diffusion problems are solved numerically to validate our proposed sixth‐order scheme. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

一维非线性Schrödinger 方程的两个无条件收敛的守恒紧致差分格式

本文对一维非线性 Schrödinger 方程给出两个紧致差分格式, 运用能量方法和两个新的分析技 巧证明格式关于离散质量和离散能量守恒, 而且在最大模意义下无条件收敛. 对非线性紧格式构造了 一个新的迭代算法, 证明了算法的收敛性, 并在此基础上给出一个新的线性化紧格式. 数值算例验证 了理论分析的正确性, 并通过外推进一步提高了数值解的精度.