二维半线性反应扩散方程的交替方向隐格式

本文研究一类二维半线性反应扩散方程的差分方法.构造了一个二层线性化交替方向隐格式.利用离散能量估计方法证明了差分格式解的存在唯一性、差分格式在离散H~1模下的二阶收敛性和稳定性.最后给出两个数值例子验证了理论分析结果.     

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

提出了求解三维抛物型方程的一个高精度显式差分格式.首先,推导了一个特殊节点处一阶偏导数(■u)/(■/t)的一个差分近似表达式,利用待定系数法构造了一个显式差分格式,通过选取适当的参数使格式的截断误差在空间层上达到了四阶精度和在时间层上达到了三阶精度.然后,利用Fourier分析法证明了当r1/6时,差分格式是稳定的.最后,通过数值试验比较了差分格式的解与精确解的区别,结果说明了差分格式的有效性.     

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.     

New conservative difference schemes for a coupled nonlinear Schrödinger system

In this paper, two conservative difference schemes for solving a coupled nonlinear Schrödinger (CNLS) system are numerically analyzed. Firstly, a nonlinear implicit two-level finite difference scheme for CNLS system is studied, then a linear three-level difference scheme for CNLS system is presented. An induction argument and the discrete energy method are used to prove the second-order convergence and unconditional stability of the linear scheme. Numerical examples show the efficiency of the new scheme.     

Long-Time Behavior of Finite Difference Solutions of a Nonlinear Schrödinger Equation with Weakly Damped

A weakly damped Schrödinger equation possessing a global attractor are considered. The dynamical properties of a class of finite difference scheme are analysed. The existence of global attractor is proved for the discrete system. The stability of the difference scheme and the error estimate of the difference solution are obtained in the autonomous system case. Finally, long-time stability and convergence of the class of finite difference scheme also are analysed in the nonautonomous system case.     

二维热传导方程的三层显式差分格式

对二维热传导方程构造了一个稳定的三层显式差分格式求其数值解,其背景源于高维热力学反问题迭代算法中对正问题小计算量算法的需求。首先建立一个含参数的一般差分格式去逼近微分方程,并得到了最优截断误差。然后导出了参数应满足的条件以保证差分格式的稳定性。最后给出了数值的例子并和其它算法进行比较,说明了格式在精度上的有效性和计算量上的优越性。     

LONG-TIME BEHAVIOR OF FINITE DIFFERENCE SOLUTIONS OF A NONLINEAR SCHRODINGER EQUATION WITH WEAKLY DAMPED

1. IntroductionThe nonlinear schr~r equation with weakly dampedwhere t = N, o > 0, together with appropriate boUndary and hatal condition, is ared inmany physical fields. The echtence of an attractor is one of the most boortant ~eristiCSfor a dissipative system. The long-tabs dynamics is completely determined by the attractorof the system. J.M. Ghidaglia[1] studied the lOng-the behavior of the nonlineaz Sequation (1.1) and proved the eAstence of a compact global attractor A in H'(n) which…     

Numerical bifurcation of a delayed diffusive food‐limited model with Dirichlet boundary condition

This paper studies the Neimark–Sacker bifurcation of a diffusive food‐limited model with a finite delay and Dirichlet boundary condition by the backward Euler difference scheme, Crank‐Nicolson difference scheme, and nonstandard finite‐difference scheme. The existence of Neimark‐Sacker bifurcation at the equilibrium is obtained. Our results show that Crank‐Nicolson and nonstandard finite‐difference schemes are superior to the backward Euler difference scheme under the means of describing approximately the dynamics of the original system. Finally, numerical examples are provided to illustrate the analytical results. Copyright © 2015 John Wiley & Sons, Ltd.     

Analysis of new conservative difference scheme for two-dimensional Rosenau-RLW equation

In this article, numerical solution for the Rosenau-RLW equation in 2D is considered and a conservative Crank–Nicolson finite difference scheme is proposed. Existence of the numerical solutions for the difference scheme has been shown by Browder fixed point theorem. A priori bound and uniqueness as well as conservation of discrete mass and discrete energy for the finite difference solutions are discussed. Unconditional stability and a second-order accuracy on both space and time of the difference scheme are proved. Numerical experiments are given to support our theoretical results.     

Numerical solutions of diffusion mathematical models with delay

In this paper an explicit numerical difference scheme for mixed problems for the delay diffusion equation is proposed, as a generalization of the classic difference scheme for the diffusion problem. A sufficient condition for the asymptotic stability of the new scheme is proved. Consistence, convergence and some properties of stability for this scheme are studied. Illustrative examples of numerical results are also included.     

Estimation of 1-dimensional nonlinear stochastic differential equations based on higher-order partial differential equation numerical scheme and its application

A method based on higher-order partial differential equation (PDE) numerical scheme are proposed to obtain the transition cumulative distribution function (CDF) of the diffusion process (numerical differentiation of the transition CDF follows the transition probability density function (PDF)), where a transformation is applied to the Kolmogorov PDEs first, then a new type of PDEs with step function initial conditions and 0, 1 boundary conditions can be obtained. The new PDEs are solved by a fourth-order compact difference scheme and a compact difference scheme with extrapolation algorithm. After extrapolation, the compact difference scheme is extended to a scheme with sixth-order accuracy in space, where the convergence is proved. The results of the numerical tests show that the CDF approach based on the compact difference scheme to be more accurate than the other estimation methods considered; however, the CDF approach is not time-consuming. Moreover, the CDF approach is used to fit monthly data of the Federal funds rate between 1983 and 2000 by CKLS model.     

广义Rosenau-Burgers方程的差分算法

对广义Rosenau-Burgers方程的初边值问题进行了数值研究,提出了新的两层隐式差分格式,得到了差分解的存在唯一性,并利用能量方法分析了该格式的二阶收敛性与无条件稳定性,并且给出数值算例进行验证.     

美式看跌期权定价的两种有限差分格式

基于Black-Scholes模型,采用指数拟合有限差分法与外推的指数拟合有限差分法对美式看跌期权价值进行了数值计算,对这两种数值方法及其与已往的显式、隐式、C-N等有限差分的优缺点进行了比较,并给出数值算例,通过对此算例做的一系列数值试验,验证了算法的有效性,并得到了一些在期权交易的实际操作中有用的结果.     

带有多项式非线性项的高维反应扩散方程有限差分格式的长时间行为

本文考虑了一类带有多项式非线性项的高维反应扩散方程．建立了一个全离散的有限差分格式,并证明了差分解的存在唯一性．分析了由差分格式生成的离散系统的动力性质,在对差分解先验估计的基础上得到了离散动力系统的整体吸引子的存在性．最后证明了差分格式的长时间稳定性和收敛性．     

A high-order compact difference scheme for 2D Laplace and Poisson equations in non-uniform grid systems

In this study, a high-order compact scheme for 2D Laplace and Poisson equations under a non-uniform grid setting is developed. Based on the optimal difference method, a nine-point compact difference scheme is generated. Difference coefficients at each grid point and source term are derived. This is accomplished through the consideration of compatibility between the partial differential equation and its difference discretization. Theoretically, the proposed scheme has third- to fourth-order accuracy; its fourth-order accuracy is achieved under uniform grid settings. Two examples are provided to examine performance of the proposed scheme. Compared with the traditional five-point difference scheme, the proposed scheme can produce more accurate results with faster convergence. Another reference scheme with the same nine-point grid stencil is derived based on the five-point scheme. The two nine-point schemes have the same coefficients for each grid points; however, their coefficients for the source term are different. The overall accuracy level of the solution resulting from the proposed scheme is higher than that of the nine-point reference scheme. It is also indicated that the smoothness of grids has significant effects on accuracy and convergence of the solutions; efforts in optimizing the grid configuration and allocation can improve solution accuracy and efficiency. Consequently, with the proposed method, solution under the non-uniform grid setting with appropriate grid allocation would be more accurate than that under the uniform-grid manipulation, with the same number of grid points.     

A three-time-level explicit difference scheme of the second order of accuracy for parabolic equations

The difference schemes of Richardson [1] and of Crank-Nicolson [2] are schemes providing second-order approximation. Richardson's three-time-level difference scheme is explicit but unstable and the Crank-Nicolson two-time-level difference scheme is stable but implicit. Explicit numerical methods are preferable for parallel computations. In this paper, an explicit three-time-level difference scheme of the second order of accuracy is constructed for parabolic equations by combining Richardson's scheme with that of Crank-Nicolson. Restrictions on the time step required for the stability of the proposed difference scheme are similar to those that are necessary for the stability of the two-time-level explicit difference scheme, but the former are slightly less onerous.Translated fromMatematicheskie Zametki, Vol. 60, No. 5, pp. 751–759, November, 1996.This research was supported by the Russian Foundation for Basic Research under grant No. 95-01-00489 and by the International Science Foundation under grants No. N8Q300 and No. JBR100.     

非线性Kundu方程的弱守恒差分格式

本文考察了一类非线性Kundu方程的周期初值问题,提出了一种弱守恒的差分格式,对其差分解作了先验估计,证明了格式的收敛性与稳定性,最后,通过数值计算检验了格式的可信性.     

On finite difference approximation of a matrix-vector product in the Jacobian-free Newton–Krylov method

The Jacobian-free Newton–Krylov (JFNK) method is a special kind of Newton–Krylov algorithm, in which the matrix-vector product is approximated by a finite difference scheme. Consequently, it is not necessary to form and store the Jacobian matrix. This can greatly improve the efficiency and enlarge the application area of the Newton–Krylov method. The finite difference scheme has a strong influence on the accuracy and robustness of the JFNK method. In this paper, several methods for approximating the Jacobian-vector product, including the finite difference scheme and the finite difference step size, are analyzed and compared. Numerical results are given to verify the effectiveness of different finite difference methods.     

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.     

A higher-order hybrid spline difference method on adaptive mesh for solving singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions

We propose a hybrid numerical scheme to discretize a class of singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions on an equidistributed grid. The hybrid difference scheme is developed by using a modified backward difference scheme in time, a combination of the cubic spline and exponential spline difference scheme in space. The proposed scheme uses a cubic spline difference scheme for the discretization of robin-boundary conditions. For the time discretization of the problem, we use the standard uniform mesh while a layer adapted equidistributed grid is generated for the spatial discretization. By equidistributing a curvature-based monitor function, the spatial adaptive grid is able to capture the presence of parabolic boundary layers without using any prior information about the solution. Parameter uniform error estimates are derived to illustrate an optimal convergence of first-order in time and second-order in space for the proposed discretization. The accuracy of the proposed scheme is confirmed by the numerical experiments that underpin the theoretical analysis.