不同类型的“追赶”法及其稳定性判别

“追赶”法是求解差分方程两点边值问题或条形矩阵(即只有几条对角线上的元素不为零的矩阵)线代数问题的有效解法。“追赶”法的主要问题是稳定性问题。在早期的工作[1,2]中,利用主对角线占优的性质证明了“追赶”法的稳定性。后来“追赶”法利用到求解比较一般的差分方程两点边值问题,并利用差分方程中系数矩阵的特征值性质证明了稳定性。在[3]中证明了:当差分方程两点边值问题是C-良态的,则正交“追赶”法是稳定的。直接利用问题的性态证明“追赶”法的稳定性是有意义的,因为有些差分方程     

广义非线性Schrdinger方程组的守恒差分格式

本文对两类广义非线性schrdinger方程组的初边值问题给出了守恒的差分格式,它保持了原来微分方程所具有的一个或两个守恒关系,文中综合地运用了差分算子的Coболев不等式,Gronwall不等式和嵌入定理对差分方程的解,作出了先验估计,在此基础上证明了差分格式的稳定性和收敛性。差分方程组是一个超越方程组,对它,我们给出追赶迭代法求解的公式,并证明了解的收敛性。     

带参数微分方程边值问题的追赶法在流体力学中的应用

本文首先介绍带参数微分方程边值问题的追赶法(包括共轭方法及拟线性化方法),然后用追赶法计算流体力学润滑方程定解问题的一个实例.其数值结果比较满意.     

二阶非线性q-差分微分方程两点边值问题两个迭代正解的存在性*

研究了二阶非线性q-差分微分方程两点边值问题,给出了系统两个正解存在的充分条件. 与其他文献中使用的不动点定理不同,文章不仅证明了该系统正解的存在性, 而且还利用单调迭代技巧给出了逼近正解的迭代格式.     

一类奇异方程两点边值问题的差分—样条校正解

本文考虑一类奇异方程两点边值问题的差分解和样条数值解法,证明了差分解,样条解分别从两侧逼近精确解,从而得到高精度的差分-样条校正解. 考虑如下形式的奇异边值问题:     

一类奇异方程两点边值问题的差分—样条校正解

本文考虑一类奇异方程两点边值问题的差分解和样条数值解法,证明了差分解,样条解分别从两侧逼近精确解,从而得到高精度的差分-样条校正解. 考虑如下形式的奇异边值问题:     

三维热传导型半导体问题的特征有限元方法和分析

本文研究三维热传导型半导体瞬态问题的特征有限元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出Ｇａｌｅｒｋｉｎ逼近;对电子,空穴浓度方程采用特征有限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近．应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

三维热传导型半导体问题的特征混合元方法和分析

本文研究三维热传导型半导体态问题的特征混合元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出混合元逼近,对电子,空穴浓度方程笔挺表限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近,应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

广义非线性Sine-Gordon方程的两个隐式差分格式

本文对一类非线性Sine-Gordon方程的初边值问题提出了两个隐式差分格式．两个隐式差分格式的精度均为O(τ~2 h~2)．我们用离散泛函分析的方法证明了格式的收敛性和稳定性,并证明了求解格式的追赶迭代法的收敛性,最后给出了数值结果．结果表明本文的格式是有效的和可靠的．     

具有积分边界条件的常微分方程边值问题的数值解

讨论了一类具有积分边界条件的二阶常微分方程非局部边值问题的数值解.对非局部积分边界条件采用了离散的多点边值问题进行逼近,通过常系数情况下解的局部性质,建立了这类边值问题的指数型差分格式,并且给出了格式的误差分析,证明了格式是一致收敛的.     

Parallel algorithms for solving tridiagonal and near-circulant systems

Many problems in mathematics and applied science lead to the solution of linear systems having circulant coefficient matrices. This paper presents a new stable method for the exact solution of non-symmetric tridiagonal circulant linear systems of equations. The method presented in this paper is quite competitive with Gaussian elimination both in terms of arithmetic operations and storage requirements. It is also competitive with the modified double sweep method. This method can be applied to solve the near-circulant tridiagonal system. In addition, the method is modified to allow for parallel processing.     

THE l~1-STABILITY OF A HAMILTONIAN-PRESERVING SCHEME FOR THE LIOUVILLE EQUATION WITH DISCONTINUOUS POTENTIALS

We study the l^1-stability of a Haxniltonian-preserving scheme, developed in [Jin and Wen, Comm. Math. Sci., 3 （2005）, 285-315], for the Liouville equation with a discontinuous potential in one space dimension. We prove that, for suitable initial data, the scheme is stable in the l^1-norm under a hyperbolic CFL condition which is in consistent with the l^1-convergence results established in [Wen and Jin, SIAM J. Numer. Anal., 46 （2008）, 2688-2714] for the same scheme. The stability constant is shown to be independent of the computational time. We also provide a counter example to show that for other initial data, in particular, the measure-valued initial data, the numerical solution may become l^1-unstable.     

非定常不可压缩Navier-Stokes方程的旋度形式压力校正格式分析

将求解不可压缩流动的旋度形式压力校正格式从Stokes方程延拓到非定常不可压缩Navier-Stokes方程.在第一步需要求解一个线性的对流-扩散方程,在第二步求解一个Stokes问题.首先给出格式的稳定性分析,然后采用P_2—P_1元分别使用标准形式的压力校正格式和旋度形式的压力校正格式进行了数值模拟,模拟结果表明,对于速度的L~2,H~1误差,两种格式几乎一样,对于压力的L~2误差,旋度形式的压力校正格式略有改进.     

海水入侵数值模拟的特征块中心差分法

本文研究二维海水入侵数值模拟的有限差分法．对关于压力的流动方程本文采用块中心差分法,对关于含盐浓度方程的对流扩散方程采用基于完全二次矩形插值的特征差分法,运用先验估计的理论和技巧得到了最佳阶L~2误差估计的结果．     

Reconstruction of a time-dependent coefficient in nonlinear Klein–Gordon equation using Bernstein spectral method

In this work, we present a spectral method for recovering an unknown time-dependent lower-order coefficient and unknown wave displacement in a nonlinear Klein–Gordon equation with overdetermination at a boundary condition. We apply the initial and boundary conditions to construct the satisfier function and use this function in a transformation to convert the main problem to a nonclassical hyperbolic equation with homogeneous initial and boundary conditions. Then, we utilize the orthonormal Bernstein basis functions to approximate the solution of the reformulated problem and use a direct technique based on the operational matrices of integration and differentiation of these basis functions together with the collocation technique to reduce the problem to a system of nonlinear algebraic equations. Regarding the perturbed measurements, the method takes advantage of the mollification method in order to derive stable numerical derivatives. Numerical simulations for solving several test examples are presented to show the applicability of the proposed method for obtaining accurate and stable results.     

Mixed covolume method for parabolic problems on triangular grids

In this paper, we present a mixed covolume method for parabolic equations on triangular grids. This method use the lowest order Raviart–Thomas (R–T) mixed finite element space as the trial space. We prove the optimal order of convergence for the approximate pressure and velocity in L²-norm. Furthermore, we obtain the quasi-optimal error estimates for the approximate pressure in L^∞-norm.     

双变量矩阵方程组对称最小二乘解的迭代算法

本文研究了求双变量线性矩阵方程组的对称最小二乘解的问题.利用求解线性代数方程组的共轭梯度法的基本思想,通过对有关矩阵和系数的变形与近似处理,建立了一种迭代算法.拓宽了共轭梯度法的适用范围.算例表明,迭代算法是有效的.     

Optimal error estimates and modified energy conservation identities of the ADI-FDTD scheme on staggered grids for 3D Maxwell’s equations

This paper is concerned with the optimal error estimates and energy conservation properties of the alternating direction implicit finite-difference time-domain(ADI-FDTD) method which is a popular scheme for solving the 3D Maxwell’s equations.Precisely,for the case with a perfectly electric conducting(PEC) boundary condition we establish the optimal second-order error estimates in both space and time in the discrete H 1-norm for the ADI-FDTD scheme,and prove the approximate divergence preserving property that if the divergence of the initial electric and magnetic fields are zero,then the discrete L 2-norm of the discrete divergence of the ADI-FDTD solution is approximately zero with the second-order accuracy in both space and time.The key ingredient is two new discrete modified energy norms which are second-order in time perturbations of two new energy conservation laws for the Maxwell’s equations introduced in this paper.Furthermore,we prove that,in addition to two known discrete modified energy identities which are second-order in time perturbations of two known energy conservation laws,the ADI-FDTD scheme also satisfies two new discrete modified energy identities which are second-order in time perturbations of the two new energy conservation laws.This means that the ADI-FDTD scheme is unconditionally stable under the four discrete modified energy norms.Experimental results which confirm the theoretical results are presented.     

Algorithms for solving inverse geophysical problems on parallel computing systems

For solving inverse gravimetry problems, efficient stable parallel algorithms based on iterative gradient methods are proposed. For solving systems of linear algebraic equations with block-tridiagonal matrices arising in geoelectrics problems, a parallel matrix sweep algorithm, a square root method, and a conjugate gradient method with preconditioner are proposed. The algorithms are implemented numerically on a parallel computing system of the Institute of Mathematics and Mechanics (PCS-IMM), NVIDIA graphics processors, and an Intel multi-core CPU with some new computing technologies. The parallel algorithms are incorporated into a system of remote computations entitled “Specialized Web-Portal for Solving Geophysical Problems on Multiprocessor Computers.” Some problems with “quasi-model” and real data are solved.     

Bernstein polynomial basis for numerical solution of boundary value problems

The purpose of this paper is to propose a computational method for the approximate solution of linear and nonlinear two-point boundary value problems. In order to approximate the solution, the expansions in terms of the Bernstein polynomial basis have been used. The properties of the Bernstein polynomial basis and the corresponding operational matrices of integration and product are utilized to reduce the given boundary value problem to a system of algebraic equations for the unknown expansion coefficients of the solution. On this approach, the problem can be solved as a system of algebraic equations. By considering a special case of the problem, an error analysis is given for the approximate solution obtained by the present method. At last, five examples are examined in order to illustrate the efficiency of the proposed method.