Unconditional stability of parallel alternating difference schemes for semilinear parabolic systems

The general alternating schemes with intrinsic parallelism for semilinear parabolic systems are studied. First we prove the a priori estimates in the discrete H¹ space of the difference solution for these schemes. Then the existence of the difference solution for these schemes follows from the fixed point principle. Finally the unconditional stability of the general alternating schemes is proved. The alternating group explicit scheme, the alternating segment explicit–implicit scheme and the alternating segment Crank–Nicolson scheme are the special cases of the general alternating schemes.     

Efficient discrete implementations for a dynamic problem of linear elasticity

For a dynamic three-dimensional linear elasticity problem in velocities-stresses, we construct efficient difference schemes on the basis of various additive decompositions of the original spatial operator. They include a difference scheme whose efficient implementation at the “predictor” stage has the property of complete conservativeness. Another class of efficient difference schemes is related to the representation of the operator as a product of triangular operators, that is, an operator analog of the LU-decomposition. The parallelism degree of these difference schemes is the same as of explicit schemes.     

三维两相渗流驱动问题迎风区域分裂显隐差分法

对三维两相渗流驱动问题提出了两种迎风区域分裂显隐差分格式．压力方程采用了七点差分格式,为了能达到实际并行计算的要求,对饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式．得到了离散l2模收敛性分析,最后给出数值试验,支撑了理论分析结果．     

对流扩散方程区域分裂特征差分方法

1引言对流扩散方程是许多物理问题的数学模型,研究其稳定的数值解法具有重要的应用价值．而标准的差分法和有限元法通常会失效,出现数值振荡．80年代,Douglas和Russel提出了特征线方法,在一定程度上克服了数值振荡,保证了数值的稳定,尤其对“对流占优”问题,更能突出特征法的优越性,并有了大量的理论成果[1,2,3]．区域分裂是一种解决大规模的科学与工程计算问题的有效方法,Dawson,Du和Dupont对热传导方程给出了非重叠区域分裂格式及分析,由于内边界的显格式,需要一定的稳定性条件Δt≤CH2;而Du等在[5]给出了抛物方程的几种区域分裂格式,对区域分裂法的     

Unconditional stability of alternating difference schemes with intrinsic parallelism for two‐dimensional parabolic systems

The difference method with intrinsic parallelism for two dimensional parabolic system is studied. The general alternating difference schemes, in particular those with variable time steplengthes, are constructed and proved to be unconditionally stable. The two dimensional alternating group explicit scheme, alternating block explicit‐implicit scheme, alternating block Crank‐Nicolson scheme and block ADI scheme are the special cases of the general schemes constructed here. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 625–636, 1999     

A parallel nearly implicit time-stepping scheme

Across-the-space parallelism still remains the most mature, convenient and natural way to parallelize large scale problems. One of the major problems here is that implicit time stepping is often difficult to parallelize due to the structure of the system. Approximate implicit schemes have been suggested to circumvent the problem (M.A. Botchev et al., Appl. Numer. Math. 31 (3) (1999) 239). These schemes have attractive stability properties and they are also very well parallelizable.The purpose of this article is to give an overall assessment of the parallelism of the method.     

高维抛物方程迎风区域分裂差分方法

结合迎风方法和区域分裂思想,采用一阶迎风、二阶修正迎风法逼近高维抛物方程的对流项.内边界处和子区域分别对应区域分裂显隐格式;并运用极值原理和嵌入定理给出了收敛性分析,最后给出数值试验,说明其实际意义.     

UNCONDITIONAL STABLE DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR SEMILINEAR PARABOLIC SYSTEMS OF DIVERGENCE TYPE

The general finite difference schemes with intrinsic parallelism for the boundary value problem of the semilinear parabolic system of divergence type with bounded measurable coefficients is studied. By the approach of the discrete functional analysis, the existence and uniqueness of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. Moreover the unconditional stability of the general difference schemes with intrinsic parallelism justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete initial data of the original problems in the discrete W_2~(2,1) (Q△) norms. Finally the convergence of the discrete vector solutions of the certain difference schemes with intrinsic parallelism to the unique generalized solution of the original semilinear parabolic problem is proved.     

Discrete models in the dynamic problem of linear elasticity and conservation laws

The dynamic problem of linear elasticity (a continuous model) is based on the momentum and moment of momentum conservation laws and hence has an additional conservation law for the total energy (the kinetic plus potential energy). Continuous models with this property are referred to as entropy models (S.K. Godunov). For difference schemes, this property is described with the use of the notion of complete conservativeness) (A.A. Samarskii). For the considered problem, we make a complete analysis of a two-parameter family of “conjugatecoordinated” two-layer difference schemes including a conservative scheme. For the latter, we construct efficient triangular-factorized implementations with the same parallelism degree as in ordinary explicit schemes. For difference schemes without the property of complete conservativeness, we discuss the role of the law of passage to the limit for the minimization of imbalance in the total energy.     

The unconditional stable and convergent difference methods with intrinsic parallelism for quasilinear parabolic systems

A kind of the general finite difference schemes with intrinsic parallelism forthe boundary value problem of the quasilinear parabolic system is studied without assum-ing heuristically that the original boundary value problem has the unique smooth vectorsolution. By the method of a priori estimation of the discrete solutions of the nonlineardifference systems, and the interpolation formulas of the various norms of the discretefunctions and the fixed-point technique in finite dimensional Euclidean space, the exis-tence and uniqueness of the discrete vector solutions of the nonlinear difference systemwith intrinsic parallelism are proved. Moreover the unconditional stability of the generalfinite difference schemes with intrinsic parallelism is justified in the sense of the continu-ous dependence of the discrete vector solution of the difference schemes on the discretedata of the original problems in the discrete w_2~(2,1) norms. Finally the convergence of thediscrete vector solutions of the certain differe     

Numerical methods for Hamilton Jacobi functional differential equations

Initial and initial boundary value problems for first order partial functional differential equations are considered. Explicit difference schemes of the Euler type and implicit difference methods are investigated. The following theoretical aspects of the methods are presented. Sufficient conditions for the convergence of approximate solutions are given and comparisons of the methods are presented. It is proved that assumptions on the regularity of given functions are the same for both the methods. It is shown that conditions on the mesh for explicit difference schemes are more restrictive than suitable assumptions for implicit methods. There are implicit difference schemes which are convergent and corresponding explicit difference methods are not convergent. Error estimates for both the methods are construted.     

THE UNCONDITIONAL STABLE DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR TWO DIMENSIONAL SEMILINEAR PARABOLIC SYSTEMS

In this paper we are going to discuss the difference schemes with intrinsic parallelism for the boundary value problem of the two dimesional semilinear parabolic systems.The unconditional stability of the general finite difference schemes with intrinsic parallelism is justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete data of the original problems in the discrete W2^(2,1) norms.Then the uniqueness of the discrete vector solution of this difference scheme follows as the consequence of the stability.     

拟线性抛物型方程组的主对角隐格式

拟线性抛物型方程组的主对角隐格式韩臻，沈隆钧，符鸿源（北京应用物理与计算数学研究所）ＡＤＩＡＧＯＮＡＬＩＭＰＬＩＣＩＴＳＣＨＥＭＥＦＯＲＱＵＡＳＩ－ＬＩＮＥＡＲＰＡＲＡＢＯＬＩＣＳＹＳＴＥＭ￥ＨａｎＺｈｅｎ；ＳｌｉｅｎＬｏｎｇ－ｊｕｎ；ＦｕＨｏｎｇ－...     

扩散方程一种无条件稳定的保正并行有限差分方法

本文针对扩散方程提出了一种保正的并行差分格式,并且这个格式为无条件稳定的.我们在每个时间层将计算区域分成许多个子区域以便于实施并行计算.格式构造中首先我们使用前两个时间层的计算结果在分区界面处通过一种非线性的保正外插来预估子区域界面值.然后在每个子区域内部使用经典的全隐格式进行计算.最后在界面处使用全隐格式进行校正（本质上这一步计算是显式计算）.我们给出了一维与二维情形下的保正并行差分格式,并相应的给出了无条件稳定性证明.数值实验显示此并行格式具有二阶数值精度,而且无条件稳定性与保正性也均在数值实验中得到验证.     

Burgers方程的一类交替分组方法

对于Burgers方程给出了一组新的Saul'yev型非对称差分格式，并用这些差分格式构造了求解非线性Burgers方程的交替分组四点方法．该算法把剖分节点分成若干组，在每组上构造能够独立求解的差分方程．因此算法具有并行本性，能直接在并行计算机上使用．章还证明了所给算法线性绝对稳定．数值试验表明，该方法使用简便，稳定性好，有很好的精度。     

THE UNCONDITIONAL CONVERGENT DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR QUASILINEAR PARABOLIC SYSTEMS WITH TWO DIMENSIONS

In the present work we are going to solve the boundary value problem for the quasilinear parabolic systems of partial differential equations with two space dimensions by the finite difference method with intrinsic parallelism.Some fundamental behaviors of general finite difference schemes with intrinsic parallelism for the mentioned problems are studied.By the method of a priori estimation of the discrete solutions of the nonlinear difference systems,and the interpolation formulas of the various norms of the discrete functions and the fixed-point technique in finite dimensional Euclidean space,the existennce of the discrete vector solutions of the nonliear difference system with intrinsic parallelism are proved .Moreover the convergence of the discrete vector solutions of these difference schemes to the unique generalizd solution of the original quasilinear parabolic problem is proved.     

ON THE CONVERGENCE OF IMPLICIT DIFFERENCE SCHEMES FOR HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are are interested in the fOllowing Cauchy problem for scalar conservation lawswhere the initial data uo E BV(R) and the flux function f 6 C'(n).It is well known that this problem may not always have a smooth global solution even ifthe i…     

Implicit Solution of a Two-Dimensional Parabolic Inverse Problem with Temperature Overspecification

Three different implicit finite difference schemes for solving the two-dimensional parabolic inverse problem with temperature overspecification are considered. These schemes are developed for indentifying the control parameter which produces, at any given time, a desired temperature distribution at a given point in the spatial domain. The numerical methods discussed, are based on the second-order (5,1) Backward Time Centered Space (BTCS) implicit formula, and the second-order (5,5) Crank-Nicolson implicit finite difference formula and the fourth-order (9,9) implicit scheme. These finite difference schemes are unconditionally stable. The (9,9) implicit formula takes a huge amount of CPU time, but its fourth-order accuracy is significant. The results of a numerical experiment are presented, and the accuracy and central processor (CPU) times needed for each of the methods are discussed and compared. The implicit finite difference schemes use more central processor times than the explicit finite difference techniques, but they are stable for every diffusion number.     

General difference schemes with intrinsic parallelism for nonlinear parabolic systems

The boundary value problem for nonlinear parabolic system is solved by the finite difference method with intrinsic parallelism. The existence of the discrete vector solution for the general finite difference schemes with intrinsic parallelism is proved by the fixed-point technique in finite-dimensional Euclidean space. The convergence and stability theorems of the discrete vector solutions of the nonlinear difference system with intrinsic parallelism are proved. The limitation vector function is just the unique generalized solution of the original problem for the parabolic system.     

带扩散项四阶抛物方程的一类本性并行差分格式

给出逼近带扩散项四阶抛物方程一组非对称差分格式,对此组非对称格式重新组合,得到了一类新的具有并行本性的算法.随后,利用矩阵法证明了算法的绝对稳定性.最后给出数值实验.