On the absolute instability of semi-implicit schemes for hydrostatic models

The dependence of the linear stability of two-time-level finite-difference semi-implicit schemes on the choice of reference temperature profile is studied. Particular vertical profiles of the temperature are considered to derive analytical conditions of stability. Analysis is made for general form of different model parameters such as the number of vertical levels and their distribution, the time step size, and the values of the viscosity coefficients. The derived conditions of stability are more restrictive than those for three-time-level schemes, but obtained necessary and sufficient condition for constant vertical lapse rates of the temperature has the form frequently applied to three-time-level schemes: the basic temperature profile should be warmer than the actual one. Performed numerical experiments show that the last restriction is neither necessary nor sufficient condition of stability for general temperature profiles.     

色散方程的一类本性并行的差分格式

对一维色散方程给出了本性并行的一般的交替差分格式，证明了该类格式的绝对稳定性已有的交替分组显格式(AGE)是该类格式的特例．作为特例，进一步得到交替分段显一隐格式(ASF-I)和交替分段Crank-Nicolson格式(ASC-N)．数值实验比较了这几个格式数值解的精确性．     

A modified Crank-Nicolson scheme with incremental unknowns for convection dominated diffusion equations

A modified Crank-Nicolson scheme based on one-sided difference approximations is proposed for solving time-dependent convection dominated diffusion equations in two-dimensional space. The modified scheme is consistent and unconditionally stable. A priori L² error estimate for the fully discrete modified scheme is derived. With the use of the incremental unknowns preconditioner at each time step, a comparison among several classical numerical schemes has been made and numerical results confirm stability and efficiency of the modified Crank-Nicolson scheme.     

Parallel difference schemes with interface extrapolation terms for quasi-linear parabolic systems

In this paper some new parallel difference schemes with interface extrapolation terms for a quasi-linear parabolic system of equations are constructed. Two types of time extrapolations are proposed to give the interface values on the interface of sub-domains or the values adjacent to the interface points, so that the unconditional stable parallel schemes with the second accuracy are formed. Without assuming heuristically that the original boundary value problem has the unique smooth vector solution, the existence and uniqueness of the discrete vector solutions of the parallel difference schemes constructed are proved. Moreover the unconditional stability of the parallel difference schemes is justified in the sense of the continuous dependence of the discrete vector solution of the schemes on the discrete known data of the original problems in the discrete W2(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the parallel difference schemes with interface extrapolation terms to the unique generalized solution of the original quasi-linear parabolic problem is proved. Numerical results are presented to show the good performance of the parallel schemes, including the unconditional stability, the second accuracy and the high parallelism.     

色散方程的一类新的并行交替分段隐格式

本文给出了一组逼近色散方程的非对称差分格式，并用这组格式和对称的Crank-Nicolson型格式构造了求解色散方程的并行交替分段差分隐格式．这个格式是无条件稳定的，能直接在并行计算机上使用．数值试验表明，这个格式有很好的精度．     

KdV方程的一类本性并行差分格式

对三阶KdV方程给出了—组非对称的差分公式,并用这些差分公式和对称的Crank-Nicolson型公式构造了一类具有本性并行的交替差分格式．证明了格式的线性绝对稳定性．对—个孤立波解、二个孤立波解和三个孤立波解的情况分别进行了数值试验,并对—个孤立波解的数值解的收敛阶和精确性进行了试验和比较．     

Monotone difference schemes stabilized by discrete mollification for strongly degenerate parabolic equations

The discrete mollification method is a convolution‐based filtering procedure suitable for the regularization of ill‐posed problems and for the stabilization of explicit schemes for the solution of PDEs. This method is applied to the discretization of the diffusive terms of a known first‐order monotone finite difference scheme [Evje and Karlsen, SIAM J Numer Anal 37 (2000) 1838–1860] for initial value problems of strongly degenerate parabolic equations in one space dimension. It is proved that the mollified scheme is monotone and converges to the unique entropy solution of the initial value problem, under a CFL stability condition which permits to use time steps that are larger than with the unmollified (basic) scheme. Several numerical experiments illustrate the performance and gains in CPU time for the mollified scheme. Applications to initial‐boundary value problems are included. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 38–62, 2012     

The stability and convergence of two linearized finite difference schemes for the nonlinear epitaxial growth model

The numerical simulation of the dynamics of the molecular beam epitaxy (MBE) growth is considered in this article. The governing equation is a nonlinear evolutionary equation that is of linear fourth order derivative term and nonlinear second order derivative term in space. The main purpose of this work is to construct and analyze two linearized finite difference schemes for solving the MBE model. The linearized backward Euler difference scheme and the linearized Crank‐Nicolson difference scheme are derived. The unique solvability, unconditional stability and convergence are proved. The linearized Euler scheme is convergent with the convergence order of O(τ + h²) and linearized Crank‐Nicolson scheme is convergent with the convergence order of O(τ² + h²) in discrete L²‐norm, respectively. Numerical stability with respect to the initial conditions is also obtained for both schemes. Numerical experiments are carried out to demonstrate the theoretical analysis. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

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

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

非线性波动方程的显式差分法

利用有界延拓法，研究了非线性波动方程周期初边值问题的显式差分解的收敛性与稳定性，避免了较难的先验估计，并放宽了非线性项的条件。     

关于色散方程的一类二阶恒稳显格式

1　引　　言对于具有高阶空间导数的发展方程 ,其显格式因结构简单 ,易于计算 ,具有明显的计算优越性 ,但已有的绝大多数显格式的稳定性条件都十分苛刻 (见 [6 ] -[1 5] ) ,远不如一般隐格式 ,使其应用受到限制 .1 994年《计算物理》中关于“色散方程的一类具任意稳定性的显格式”一文 (见 [1 4 ] ) ,把色散方程显格式的稳定性条件提高到了可以任意选择的程度 ,但截断误差仅为 O(τ+h) .本文构造了新一类双参数显式差分格式 ,它是绝对稳定的 ,且其截断误差是 O(τ+h2 ) ,它结构简单 ,易于实现计算 ,利于实际应用 .我们用数值例子验证了理论…     

PARALLEL COMPUTING METHOD OF PURE ALTERNATIVE SEGMENT EXPLICIT-IMPLICIT DIFFERENCE SCHEME FOR NONLINEAR LELAND EQUATION

The research on the numerical solution of the nonlinear Leland equation has important theoretical significance and practical value. To solve nonlinear Leland equation, this paper offers a class of difference schemes with parallel nature which are pure alternative segment explicit-implicit(PASE-I) and implicit-explicit(PASI-E) schemes. It also gives the existence and uniqueness,the stability and the error estimate of numerical solutions for the parallel difference schemes. Theoretical analysis demonstrates that PASE-I and PASI-E schemes have obvious parallelism, unconditionally stability and second-order convergence in both space and time. The numerical experiments verify that the calculation accuracy of PASE-I and PASI-E schemes are better than that of the existing alternating segment Crank-Nicolson scheme, alternating segment explicit-implicit and implicit-explicit schemes. The speedup of PASE-I scheme is 9.89, compared to classical Crank-Nicolson scheme. Thus the schemes given by this paper are high efficient and practical for solving the nonlinear Leland equation.     

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.     

三阶非线性KdV方程的交替分段显-隐差分格式

对三阶非线性KdV方程给出了一组非对称的差分公式,用这些差分公式与显、隐差分公式组合,构造了一类具有本性并行的交替分段显-隐格式A·D2证明了格式的线性绝对稳定性．对1个孤立波解、2个孤立波解的情况分别进行了数值试验．数值结果显示,交替分段显-隐格式稳定,有较高的精确度．     

Two difference schemes for the numerical solution of Maxwell’s equations as applied to extremely and super low frequency signal propagation in the Earth-ionosphere waveguide

Two explicit two-time-level difference schemes for the numerical solution of Maxwell’s equations are proposed to simulate propagation of small-amplitude extremely and super low frequency electromagnetic signals (200 Hz and lower) in the Earth-ionosphere waveguide with allowance for the tensor conductivity of the ionosphere. Both schemes rely on a new approach to time approximation, specifically, on Maxwell’s equations represented in integral form with respect to time. The spatial derivatives in both schemes are approximated to fourth-order accuracy. The first scheme uses field equations and is second-order accurate in time. The second scheme uses potential equations and is fourth-order accurate in time. Comparative test computations show that the schemes have a number of important advantages over those based on finite-difference approximations of time derivatives. Additionally, the potential scheme is shown to possess better properties than the field scheme.     

On a difference scheme of fourth order of accuracy for the Bitsadze–Samarskii type nonlocal boundary value problem

The Bitsadze–Samarskii type nonlocal boundary value problem for the differential equation in a Hilbert space H with the self‐adjoint positive definite operator A with a closed domain D(A) ? H is considered. Here, f(t) be a given abstract continuous function defined on [0,1] with values in H, φ and ψ be the elements of D(A), and λ_j are the numbers from the set [0,1]. The well‐posedness of the problem in Hölder spaces with a weight is established. The coercivity inequalities for the solution of the nonlocal boundary value problem for elliptic equations are obtained. The fourth order of accuracy difference scheme for approximate solution of the problem is presented. The well‐posedness of this difference scheme in difference analogue of Hölder spaces is established. For applications, the stability, the almost coercivity, and the coercivity estimates for the solutions of difference schemes for elliptic equations are obtained. Mathematical Methods in the Applied Sciences. Copyright © 2012 John Wiley & Sons, Ltd.     

提高反应—扩散方程有限差分格式的稳定性问题

This paper deals with the special nonlinear reaction-diffusion equation.The finite difference scheme with incremental unknowns approximating to the differential equation (2.1) is set up by means of introducing incremental unknowns methods.Through the stability analyzing for the scheme,it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.     

耗散对称正则长波方程的有限差分逼近

本文对耗散对称正则长波方程的初边值问题进行了数值研究,提出了一个两层隐式Crank-Nicolson差分格式,讨论了差分解的存在唯一性,并利用能量方法分析了该格式的二阶收敛性与稳定性,数值算例表明本文的格式是可靠的.     

A second order difference scheme with nonuniform rectangular meshes for nonlinear parabolic system

In this paper, a difference scheme with nonuniform meshes is proposed for the initial-boundary problem of the nonlinear parabolic system. It is proved that the difference scheme is second order convergent in both space and time.     

五维热传导方程的一族两层显格式

构造了五维热传导方程的一族两层显格式,证明了当截断误差阶为O（τ＋h2）时,其稳定性条件为网比r=hτ2≤21,优于同类的其它显格式,当截断误差阶为O（τ2＋h2）时,可以得到一个简洁而实用的二阶精度的两层显格式.