Grid approximation of singularly perturbed parabolic reaction-diffusion equations on large domains with respect to the space and time variables

In an unbounded (with respect to x and t) domain (and in domains that can be arbitrarily large), an initial-boundary value problem for singularly perturbed parabolic reaction-diffusion equations with the perturbation parameter ε² multiplying the higher order derivative is considered. The parameter ε takes arbitrary values in the half-open interval (0, 1]. To solve this problem, difference schemes on grids with an infinite number of nodes (formal difference schemes) are constructed that converge ε-uniformly in the entire unbounded domain. To construct these schemes, the classical grid approximations of the problem on the grids that are refined in the boundary layer are used. Schemes on grids with a finite number of nodes (constructive difference schemes) are also constructed for the problem under examination. These schemes converge for fixed values of ε in the prescribed bounded subdomains that can expand as the number of grid points increases. As ε → 0, the accuracy of the solution provided by such schemes generally deteriorates and the size of the subdomains decreases. Using the condensing grid method, constructive difference schemes that converge ε-uniformly are constructed. In these schemes, the approximation accuracy and the size of the prescribed subdomains (where the schemes are convergent) are independent of ε and the subdomains may expand as the number of nodes in the underlying grids increases.     

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.     

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.     

THE FINITE DIFFERENCE STREAMLINE DIFFUSION METHODS FOR TIME-DEPENDENT CONVECTION-DIFFUSION EQUATIONS

In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L~2-norm.     

高阶Schrdinger方程的一族恒稳差分格式

1 引言 高阶Schrdinger方程在量子力学、非线性光学及流体力学中有着广泛的应用,其最简单的模型方程为,即 (其中,i=(-1)~(1/2),m为正整数) (1) 文[2-3]研究了方程(1)的辛算法及差分解法。文[2]先将方程(1)写成Hamilton形     

二维抛物型方程的高精度多重网格解法

提出了数值求解二维抛物型方程的一种新的高精度加权平均紧隐格式，利用Fourier分析方法证明了该格式是无条件稳定的，为了克服传统迭代法在求解隐格式是收敛速度慢的缺陷，利用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的精确性和可靠性。     

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.     

Monotonicity criteria for difference schemes designed for hyperbolic equations

Previously formulated monotonicity criteria for explicit two-level difference schemes designed for hyperbolic equations (S.K. Godunov’s, A. Harten’s (TVD schemes), characteristic criteria) are extended to multileveled, including implicit, stencils. The characteristic monotonicity criterion is used to develop a universal algorithm for constructing high-order accurate nonlinear monotone schemes (for an arbitrary form of the desired solution) based on their analysis in the space of grid functions. Several new fourth-to-third-order accurate monotone difference schemes on a compact three-level stencil and nonexpanding (three-point) stencils are proposed for an extended system, which ensures their monotonicity for both the desired function and its derivatives. The difference schemes are tested using the characteristic monotonicity criterion and are extended to systems of hyperbolic equations.     

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.     

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     

Unconditional and optimal <Emphasis Type="Bold">H</Emphasis><Superscript>2</Superscript>-error estimates of two linear and conservative finite difference schemes for the Klein-Gordon-Schrödinger equation in high dimensions

The focus of this paper is on the optimal error bounds of two finite difference schemes for solving the d-dimensional (d = 2, 3) nonlinear Klein-Gordon-Schrödinger (KGS) equations. The proposed finite difference schemes not only conserve the mass and energy in the discrete level but also are efficient in practical computation because only two linear systems need to be solved at each time step. Besides the standard energy method, an induction argument as well as a ‘lifting’ technique are introduced to establish rigorously the optimal H ²-error estimates without any restrictions on the grid ratios, while the previous works either are not rigorous enough or often require certain restriction on the grid ratios. The convergence rates of the proposed schemes are proved to be at O(h ² + τ ²) with mesh-size h and time step τ in the discrete H ²-norm. The analysis method can be directly extended to other linear finite difference schemes for solving the KGS equations in high dimensions. Numerical results are reported to confirm the theoretical analysis for the proposed finite difference schemes.     

New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

Schémas compacts d'ordre élevé: application aux problèmes bidimensionnels de diffusion-convection instationnaire I

For unsteady 2D diffusion-convection problems, we present two classes of compact difference schemes of order 2 in time and 4 in space. These FDS are essentially derived from 1D schemes, extensively analyzed in our previous paper [6]. We propose two approaches: construction of 2D schemes as product of 1D schemes and global formulation of 2D schemes. Part II (by Michel Fournié) focuses on the development and analysis of global schemes with the assistance of symbolic computation software (AXIOM).     

求解二阶波动方程的三次样条差分方法

有限差分法在求解二阶波动方程初边值问题过程中通常受到精度和稳定性的限制.本文对二阶波动方程的时间、空间项分别采用三次样条公式进行离散,推导出精度分别为O(τ2+h2),0(τ2+h4),O(τ4+h2)和O(τ4+h4)的四种三层隐式差分格式,以及与之相匹配的第一个时间步的同阶离散格式,并采用Fourier方法分析了格...     

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.     

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

提出了一族三维热传导方程的两层显式差分格式,当截断误差阶为O(Δt+(Δx)2)时,稳定性条件为网格比r=Δt/(Δx)2=Δt/(Δy)2=Δt/(Δz)2≤1/2,优于其他显式差分格式。而当截断误差阶为O((Δt)2+(Δx)4)时,稳定性条件为r≤1/6,包含了已有的结果。     

一类时空二阶精度高分辨率MmB差分格式的构造及数值试验

1．引言考虑如下二维双曲型守恒律初值问题的数值解．H．M．Wu和S．L．Yang在文山中给出了MmB差分格式的定义如下：给定（．1）M差分格式定义．若则称格式（1．2）为MmB差分格式．这里BmB表示局部MaximumandminimumBounds．由定义可知，若差分格式（1．2）可写为形式且。＼P’三0，＞。：r’一1．则格式（1．4）为MmB差分格式．j＝l文山构造了二维双曲型守恒律的二类二阶精度的MmB差分格式，使构造二维高分辨格式有了新的突破，但他们是从标量线性双曲型守恒律出发，然后把结果推广到非线性情形．本文直接从二维非线性双曲型守恒律…     

空间四阶的时间亚扩散方程的有限差分方法

提出了两个求解空间四阶的时间亚扩散方程的数值方法,其误差阶分别为O（τ＋h2）和O（τ2＋h2）.通过Fourier方法,发现两个差分格式均为无条件稳定的.最后,通过数值例子,验证了两个算法的有效性.     

Convergence of Iterative Difference Method with Nonuniform Meshes for Quasilinear Parabolic Systems

In this paper, we study the general difference schemes with nonuniform meshes for the following problem: u_t = A(x,t,u,u_x)u_{xx}, + f(x,t,u,u_x), 0 < x < l, 0 < t ≤ T \qquad (1) u(0,t) = u(l ,t) = 0, 0 < t ≤ T \qquad\qquad (2) u(x,0) = φ(x), 0 ≤ x ≤ l \qquad\qquad (3) where u, φ, and f are m-dimensional vector valued functions, u_t = \frac{∂u}{∂t}, u_x = \frac{∂u}{∂x}, u_{xx} = \frac{∂²u}{∂_x²}. In the practical computation, we usually use the method of iteration to calculate the approximate solutions for the nonlinear difference schemes. Here the estimates of the iterative sequence constructed from the iterative difference schemes for the problem (1)-(3) is proved. Moreover, when the coefficient matrix A = A(x, t, u) is independent of u_x, t he convergence of the approximate difference solution for the iterative difference schemes to the unique solution of the problem (1)-(3) is proved without imposing the assumption of heuristic character concerning the existence of the unique smooth solution for the original problem (1)-(3).     

Algorithm for solving a class of one-dimensional variational inequalities

We consider one-dimensional variational inequalities with end constraints. An exact difference scheme and truncated difference schemes of any order of accuracy are constructed for this problem. The accuracy of the rank-m truncated scheme in the grid norm of C is O(h^2m+2). An algorithm for the implementation of the difference schemes is proposed. The algorithm reduces to two sweeps.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 64, pp. 24–30, 1988.