一类线性双曲型方程Neumann边值问题的高阶差分格式

对一维Neumann边界条件的线性双曲方程,利用有限差分方法建立高阶差分格式.由方程和边界条件得到在空间边界点的三阶和五阶导数值,进而分别在内点和边界点建立三点和两点紧差分格式,其截断误差关于时间和空间分别为二阶和四阶;利用离散的能量估计方法,分析差分格式的收敛性和稳定性;通过数值算例,验证理论分析结果.     

一维线性Klein-Gordon方程Neumann边值问题的高阶差分格式

本文主要研究非线性Klein-Gordon方程Neumann边值问题的高阶差分格式.利用边界条件及非线性Klein-Gordon方程,得到其在空间上的三阶与五阶导数的边界值,进而分别在内点和边界点建立三点和两点紧差分格式.借助能量估计、Gronwall和Schwarz不等式、数学归纳法等技巧进行分析,得到截断误差是关于时间和空间上的二阶和四阶收敛.通过理论分析差分格式的收敛性和稳定性以及数值算例,验证了理论分析结果.     

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

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

一类非线性延迟抛物偏微分方程的紧致差分格式

对于一类带有Dirichlet边界条件的延迟非线性抛物型偏微分方程的初边值问题建立了一个紧差分格式,用能量分析法证明该差分格式在L_∞范数下是无条件收敛的,且收敛阶为O(τ~2+h~4).最后,通过数值算例验证了理论结果.     

基于不光滑边界的变系数抛物型方程的高精度紧格式

本文讨论基于不光滑边界的变系数抛物型方程求解的高精度紧格式.首先构造一般变系数抛物型方程的高精度紧格式,并在理论上证明格式具有空间方向四阶精度.然后针对非光滑边界条件,引入局部网格加密技巧在奇异点附近进行不均匀的网格加密.数值实验以期权定价中Black-Scholes偏微分方程的求解为例,验证高精度紧格式用于光滑边界条件的微分方程离散可以达到四阶精度.对于处理非光滑边界条件,网格局部加密技巧能有效的提高数值解精度,使得高精度紧格式用于定价欧式期权可以接近四阶精度.     

关于热传导方程半离散差分格式的一个注记

本文研究了热传导方程初边值问题的半离散化差分格式直接解算法.分别从Dirichlet和Neumann边界条件出发,直接由空间差分格式导出与时间相关的一阶常微分方程组,随后通过正/余弦变换获得了原方程的半解析解,并给出了相关收敛性分析.并对中心差分格式和紧差分格式的精度差异,通过矩阵特征值理论给出了相关原因分析.另外,对于二维热传导方程初边值问题,应用矩阵张量积运算,该直接解算法可直接演变成二重正(余)弦变换.该方法由于不涉及时间上的离散,从而具有较好的计算效率.     

拟线性抛物型方程奇异摄动问题的数值解法

本文讨论拟线性抛物型方程奇异摄动问题的差分解法,在非均匀网格上建立了线性三层差分格式,并证明了在离散的L₂范数意义下格式的一致收敛性,最后给出了一些数值例子.     

求解Riesz空间分数阶扩散方程的一种新的数值方法

本文利用Diethelm方法构造了一种逼近Riesz空间分数阶导数的O（△x^3-α）格式,其中1 < α < 2,△x是空间步长.进一步对一阶时间导数采用Crank-Nicolson方法离散,得到了求解Riesz空间分数阶扩散方程的一种新的有限差分格式,并用矩阵方法证明了稳定性和收敛性,其误差估计为O（△t²+△x^3-α）,其中△t为时间步长.最后,数值算例验证了差分格式的正确性和有效性.     

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

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

求解对流扩散方程的一致四阶紧致格式

目前,许多高精度差分格式,由于未成功地构造与其精度匹配的稳定的边界格式,不得不采用低精度的边界格式.本文针对对流扩散方程证明了存在一致四阶紧致格式,它的边界点的计算格式和内点的计算格式的截断误差主项保持一致,给出了具体内点和边界格式;并分析了此半离散格式的渐近稳定性.数值结果表明该格式是四阶精度;在对流占优情况下,本文边界格式的数值结果比四阶精度的显式差分格式的的数值结果的数值振荡小,取得了不错的效果,理论结果得到了数值验证;驱动方腔数值结果显示,本文对N-S方程的离散格式具有很好的可靠性,适合对复杂流体流动的数值模拟和研究.     

The Conservation and Convergence of Two Finite Difference Schemes for KdV Equations with Initial and Boundary Value Conditions

Korteweg-de Vries equation is a nonlinear evolutionary partial differential equation that is of third order in space. For the approximation to this equation with the initial and boundary value conditions using the finite difference method, the difficulty is how to construct matched finite difference schemes at all the inner grid points. In this paper, two finite difference schemes are constructed for the problem. The accuracy is second-order in time and first-order in space. The first scheme is a two-level nonlinear implicit finite difference scheme and the second one is a three-level linearized finite difference scheme. The Browder fixed point theorem is used to prove the existence of the nonlinear implicit finite difference scheme. The conservation, boundedness, stability, convergence of these schemes are discussed and analyzed by the energy method together with other techniques. The two-level nonlinear finite difference scheme is proved to be unconditionally convergent and the three-level linearized one is proved to be conditionally convergent. Some numerical examples illustrate the efficiency of the proposed finite difference schemes.     

Fourth‐order compact difference schemes for 1D nonlinear Kuramoto–Tsuzuki equation

In this article, first, we establish some compact finite difference schemes of fourth‐order for 1D nonlinear Kuramoto–Tsuzuki equation with Neumann boundary conditions in two boundary points. Then, we provide numerical analysis for one nonlinear compact scheme by transforming the nonlinear compact scheme into matrix form. And using some novel techniques on the specific matrix emerged in this kind of boundary conditions, we obtain the priori estimates and prove the convergence in norm. Next, we analyze the convergence and stability for one of the linearized compact schemes. To obtain the maximum estimate of the numerical solutions of the linearized compact scheme, we use the mathematical induction method. The treatment is that the convergence in norm is obtained as well as the maximum estimate, further the convergence in norm. Finally, numerical experiments demonstrate the theoretical results and show that one of the linearized compact schemes is more accurate, efficient and robust than the others and the previous. It is worthwhile that the compact difference methods presented here can be extended to 2D case. As an example, we present one nonlinear compact scheme for 2D Ginzburg–Landau equation and numerical tests show that the method is accurate and effective. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2080–2109, 2015     

A FINITE DIFFERENCE METHOD FOR THE HEAT EQUATION WITH A NONLINEAR BOUNDARY CONDITION

This article is concerned with the numerical solution to a parabolic equation with a kind of nonlinear boundary conditions. A difference scheme is constructed by the method of reduction of order on uniform mesh to solve the problem. It is proved that the difference scheme is uniquely solvable and uncon-ditionaUy convergent with the convergence order 2 in both space and time in an energy norm. An effective iterative algorithm is given and a numerical example is presented to demonstrate the theoretical results.     

热方程非线性边界条件问题的有限差分方法

This article is concerned with the numerical solution to a parabolic equation with a kind of nonlinear boundary conditions. A difference scheme is constructed by the method of reduction of order on uniform mesh to solve the problem. It is proved that the difference scheme is uniquely solvable and uncon-ditionaUy convergent with the convergence order 2 in both space and time in an energy norm. An effective iterative algorithm is given and a numerical example is presented to demonstrate the theoretical results.     

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

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

Compact difference schemes for heat equation with Neumann boundary conditions

In this article, two recent proposed compact schemes for the heat conduction problem with Neumann boundary conditions are analyzed. The first difference scheme was proposed by Zhao, Dai, and Niu (Numer Methods Partial Differential Eq 23, (2007), 949–959). The unconditional stability and convergence are proved by the energy methods. The convergence order is O(τ² + h^2.5) in a discrete maximum norm. Numerical examples demonstrate that the convergence order of the scheme can not exceeds O(τ² + h³). An improved compact scheme is presented, by which the approximate values at the boundary points can be obtained directly. The second scheme was given by Liao, Zhu, and Khaliq (Methods Partial Differential Eq 22, (2006), 600–616). The unconditional stability and convergence are also shown. By the way, it is reported how to avoid computing the values at the fictitious points. Some numerical examples are presented to show the theoretical results. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A parameter robust numerical method for a two dimensional reaction-diffusion problem

In this paper a singularly perturbed reaction-diffusion partial differential equation in two space dimensions is examined. By means of an appropriate decomposition, we describe the asymptotic behaviour of the solution of problems of this kind. A central finite difference scheme is constructed for this problem which involves an appropriate Shishkin mesh. We prove that the numerical approximations are almost second order uniformly convergent (in the maximum norm) with respect to the singular perturbation parameter. Some numerical experiments are given that illustrate in practice the theoretical order of convergence established for the numerical method.

     

Numerical method for a singularly perturbed convection-diffusion problem with delay

This paper deals with the singularly perturbed boundary value problem for a linear second-order delay differential equation. For the numerical solution of this problem, we use an exponentially fitted difference scheme on a uniform mesh which is accomplished by the method of integral identities with the use of exponential basis functions and interpolating quadrature rules with weight and remainder term in integral form. It is shown that one gets first order convergence in the discrete maximum norm, independently of the perturbation parameter. Numerical results are presented which illustrate the theoretical results.     

A linearized high‐order difference scheme for the fractional Ginzburg–Landau equation

The numerical solution for the one‐dimensional complex fractional Ginzburg–Landau equation is considered and a linearized high‐order accurate difference scheme is derived. The fractional centered difference formula, combining the compact technique, is applied to discretize fractional Laplacian, while Crank–Nicolson/leap‐frog scheme is used to deal with the temporal discretization. A rigorous analysis of the difference scheme is carried out by the discrete energy method. It is proved that the difference scheme is uniquely solvable and unconditionally convergent, in discrete maximum norm, with the convergence order of two in time and four in space, respectively. Numerical simulations are given to show the efficiency and accuracy of the scheme. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 105–124, 2017