任意精度的三点紧致显格式及其在CFD中的应用

通过在泰勒级数展开中运用逐阶迭代的方法,推导出了空间任意精度的三点紧致显格式的表达式,又由Fourier分析法得到了格式的数值弥散和耗散特性．与以往的高精度紧致差分格式不同,提出的格式不用隐式求解代数方程组并且可以达到任意精度．通过方波问题和顶盖方腔流的算例表明,格式在稀疏网格下可以得到很高的精度,不仅能节省计算量,而且易于编程,有很高的计算效率．     

一维非线性Schrödinger 方程的两个无条件收敛的守恒紧致差分格式

本文对一维非线性 Schrödinger 方程给出两个紧致差分格式, 运用能量方法和两个新的分析技 巧证明格式关于离散质量和离散能量守恒, 而且在最大模意义下无条件收敛. 对非线性紧格式构造了 一个新的迭代算法, 证明了算法的收敛性, 并在此基础上给出一个新的线性化紧格式. 数值算例验证 了理论分析的正确性, 并通过外推进一步提高了数值解的精度.     

Uniform convergent monotone iterates for semilinear singularly perturbed parabolic problems

This paper deals with discrete monotone iterative methods for solving semilinear singularly perturbed parabolic problems. Monotone sequences, based on the accelerated monotone iterative method, are constructed for a nonlinear difference scheme which approximates the semilinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. An analysis of uniform convergence of the monotone iterative method to the solutions of the nonlinear difference scheme and continuous problem is given. Numerical experiments are presented.     

Numerical performance of parallel group explicit solvers for the solution of fourth order elliptic equations

Many applications in applied mathematics and engineering involve numerical solutions of partial differential equations (PDEs). Various discretisation procedures such as the finite difference method result in a problem of solving large, sparse systems of linear equations. In this paper, a group iterative numerical scheme based on the rotated (skewed) five-point finite difference discretisation is proposed for the solution of a fourth order elliptic PDE which represents physical situations in fluid mechanics and elasticity. The rotated approximation formulas lead to schemes with lower computational complexities compared to the centred approximation formulas since the iterative procedure need only involve nodes on half of the total grid points in the solution domain. We describe the development of the parallel group iterative scheme on a cluster of distributed memory parallel computer using Message-Passing Interface (MPI) programming environment. A comparative study with another group iterative scheme derived from the centred difference formula is also presented. A detailed performance analysis of the parallel implementations of both group methods will be reported and discussed.     

解非定常不可压缩N-S方程的迭代压力Poisson方程法

１．引 言 数值求解不可压缩流体流动问题可以采用原始变量的方程作为控制方程,也可以用涡量一流函数方程作为控制方程．直接求解原始变量的不可压缩 Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程存在一个主要困难：速度向量在每一时刻都必须满足零散度约束条件,即不可压缩性连续方程．用涡量一流函数方程求解时,连续方程自动满足,所以不存在约束条件的问题,但涡量的边界条件比较难处理,且不易应用于三维问题和带有自由表面或其它流体交界面的问题． 解决上述速度向量必须满足零散度约束条件的困难的方法有：人工压缩法［３,１７］;压力Ｐｏｉｓ…     

Symbolic computation of high order compact difference schemes for three dimensional linear elliptic partial differential equations with variable coefficients

We present a symbolic computation procedure for deriving various high order compact difference approximation schemes for certain three dimensional linear elliptic partial differential equations with variable coefficients. Based on the Maple software package, we approximate the leading terms in the truncation error of the Taylor series expansion of the governing equation and obtain a 19 point fourth order compact difference scheme for a general linear elliptic partial differential equation. A test problem is solved numerically to validate the derived fourth order compact difference scheme. This symbolic derivation method is simple and can be easily used to derive high order difference approximation schemes for other similar linear elliptic partial differential equations.     

Uniform quadratic convergence of monotone iterates for nonlinear singularly perturbed parabolic problems

This paper deals with a monotone iterative method for solving nonlinear singularly perturbed parabolic problems. Monotone sequences, based on the method of upper and lower solutions, are constructed for a nonlinear difference scheme which approximates the nonlinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. The monotone sequences possess quadratic convergence rate. An analysis of uniform convergence of the monotone iterative method to the solutions of the nonlinear difference scheme and to the continuous problem is given. Numerical experiments are presented.     

Monotone iterates with quadratic convergence rate for solving weighted average approximations to semilinear parabolic problems

This paper deals with discrete monotone iterative methods for solving semilinear singularly perturbed parabolic problems. Monotone sequences, based on the accelerated monotone iterative method, are constructed for a nonlinear difference scheme which approximates the semilinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. An analysis of convergence of the monotone iterative method to the solutions of the nonlinear difference scheme is given. Numerical experiments are presented.     

交错网格紧致差分格式和满足等价性的压力Poisson方程

1．引言随着电子计算机的发展，越来越多的实际问题数值模拟成为现实，但还有很多非线性问题数值计算时间太长，内存要求过大．数值方法的改进可使计算量和存储量大大减少，例如，对二维非定常问题，要使误差达到N－4量级，二阶格式计算点数为（N2）3，（包括时间方向），而四阶格式计算点数仅为N3，差N3倍！而计算量差的倍数更多．当N＝16时N3＝4096，当N＝256时，N31678万．紧致差分格式具有精度高，差分式基点少，＜线性）稳定性好，对高频波分辨率高，边界差分点少等优点【’，‘，’，’。’，’，‘’］，本文中的格式基点数为3，…     

High order parameter-robust numerical method for singularly perturbed reaction-diffusion problems

We present a high order parameter-robust finite difference method for singularly perturbed reaction-diffusion problems. The problem is discretized using a suitable combination of fourth order compact difference scheme and central difference scheme on generalized Shishkin mesh. The convergence analysis is given and the method is proved to be almost fourth order uniformly convergent in maximum norm with respect to singular perturbation parameter ε. Numerical experiments are conducted to demonstrate the theoretical results.     

From the fitting techniques to accurate schemes for the Liouville‐Bratu‐Gelfand problem

We are interested in numerical methods for the Liouville‐Bratu‐Gelfand problem. The ideas and techniques developed here to construct the schemes are inspired from the fitted method and the so‐called compact exponentially fitted method. Some of those schemes can be viewed as extensions of both the Buckmire scheme and the standard scheme which results from the use of the standard finite‐difference procedures. We study and compare computationally the accuracy of methods introduced here. It is also mentioned that the Buckmire's techniques and the standard scheme are a particular case of the fitted method. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Strong Convergence of an Iterative Scheme by a New Type of Projection Method for a Family of Quasinonexpansive Mappings

We deal with a common fixed point problem for a family of quasinonexpansive mappings defined on a Hilbert space with a certain closedness assumption and obtain strongly convergent iterative sequences to a solution to this problem. We propose a new type of iterative scheme for this problem. A feature of this scheme is that we do not use any projections, which in general creates some difficulties in practical calculation of the iterative sequence. We also prove a strong convergence theorem by the shrinking projection method for a family of such mappings. These results can be applied to common zero point problems for families of monotone operators.     

A NEW APPROACH TO SOLVE SYSTEMS OF LINEAR EQUATIONS

1. IntroductionThe new aPProaCh is based on the analysis of the motion of a damped harmonic oscillatorin the gravitational field 11]. The associated equation of motion ismXtt + oXt + aX = b (1)where X = X(t), is the one dimensions di8PlaCement of Of a mass m under a dissipation(o > 0), a ~nic potential (a > 0) and a constant acceleration (b, gravitational field). Thetotal energy variation is given by the equationwhereThe solution of the motion equation (1) is given by the sum of two contr…     

Compact difference methods applied to initial-boundary value problems for mixed systems

Summary. An initial--boundary value problem to a system of nonlinear partial differential equations, which consists of a hyperbolic and a parabolic part, is taken into consideration. The problem is discretised by a compact finite difference method. An approximation of the numerical solution is constructed, at which the difference scheme is linearised. Nonlinear convergence is proved using the stability of the linearised scheme. Finally, a computational experiment for a noncompact scheme is presented. Received May 20, 1995     

High-order compact finite difference scheme for option pricing in stochastic volatility models

We derive a new high-order compact finite difference scheme for option pricing in stochastic volatility models. The scheme is fourth order accurate in space and second order accurate in time. Under some restrictions, theoretical results like unconditional stability in the sense of von Neumann are presented. Where the analysis becomes too involved we validate our findings by a numerical study. Numerical experiments for the European option pricing problem are presented. We observe fourth order convergence for non-smooth payoff.     

求解一维定常对流扩散反应方程的混合型紧致差分格式

借助显式紧致格式和隐式紧致格式的思想,基于截断误差余项修正,并结合原方程本身,构造出了一种求解一维定常对流扩散反应方程的高精度混合型紧致差分格式.格式仅用到三个点上的未知函数值及一阶导数值,而一阶导数值利用四阶Pade格式进行计算,格式整体具有四阶精度.数值实验结果验证了格式的精确性和可靠性.     

Weak Convergence of an Iterative Method for Pseudomonotone Variational Inequalities and Fixed-Point Problems

We consider an iterative scheme for finding a common element of the set of solutions of a pseudomonotone, Lipschitz-continuous variational inequality problem and the set of common fixed points of N nonexpansive mappings. The proposed iterative method combines two well-known schemes: extragradient and approximate proximal methods. We derive a necessary and sufficient condition for weak convergence of the sequences generated by the proposed scheme.     

Strong Convergence Theorems for Variational Inequalities and Fixed Point Problems of Asymptotically Strict Pseudocontractive Mappings in the Intermediate Sense

The purpose of this paper is to investigate the problem of finding a common element of the set of fixed points of an asymptotically strict pseudocontractive mapping in the intermediate sense and the set of solutions of the variational inequality problem for a monotone, Lipschitz continuous mapping. We introduce a modified hybrid Mann iterative scheme with perturbed mapping which is based on well-known CQ method, Mann iteration method and hybrid (or outer approximation) method. We establish a strong convergence theorem for three sequences generated by this modified hybrid Mann iterative scheme with perturbed mapping. Utilizing this theorem, we also design an iterative process for finding a common fixed point of two mappings, one of which is an asymptotically strict pseudocontractive mapping in the intermediate sense and the other taken from the more general class of Lipschitz pseudocontractive mappings.     

An iterative approach to solve a nonlinear moving boundary problem describing the solvent diffusion within glassy polymers

This paper consists in studying a mathematical model of solvent diffusion through the glassy polymers as a one-dimensional moving boundary problem with kinetic undercooling. We establish an iterative variable time-step method based on a nonstandard finite difference (NSFD) scheme to solve the considered moving boundary problem. The monotonicity and positivity of the numerical solution are proved. The numerical approach is investigated for three test problems composed of constant and inconstant diffusion coefficients for different values of parameters to demonstrate the validity and ability of the method.     

New methods for linear inequalities

The iterative method of Cimmino for solving linear equations is generalized to linear inequalities. We also present a Richardson-type iterative method for solving the inequality problem, which includes the generalized Cimmino scheme. Convergence proofs are provided.