一类非线性离散系统的指数二分性及其在数值分析和计算中的应用

本文中我们考虑一类二阶非线性常微分方程的边值问题的迎风差分格式.我们运用奇异摄动方法构造了该迎风差分方程解的渐近近似,并利用指数二分性理论证明了有一个低阶方程其解是该迎风方程式的在边界外的一个良好近似.我们还构造了校正项,使校正项与低阶方程的解之和是一个渐近近似.最后一些数值例子用于显示本文方法的应用.     

Uniform methods for semilinear problems with attractive boundary turning points

An upwind finite‐difference scheme for the numerical solution of semilinear convection‐diffusion problems with attractive boundary turning points is considered. We show that the maximum nodal error is bounded by a special weighted ℓ₁‐type norm of the truncation error. This result is used to establish uniform convergence with respect to the perturbation parameter on Shishkin meshes.     

油藏数值模拟中动边值问题的迎风差分格式

考虑了三维油藏数值模拟中的动边值问题,对压力方程,给出中心差分格式;对饱和度方程给出隐式迎风差分格式及修正的迎风差分格式,并证明了格式的收敛性。数值算例与理论结果是一致的。     

Algorithm composition scheme for problems in composite domains based on the difference potential method

An algorithm composition scheme for the numerical solution of boundary value problems in composite domains is proposed and illustrated using an example. The scheme requires neither difference approximations of the boundary conditions nor matching conditions on the boundary between the subdomains. The scheme is suited for multiprocessor computers.     

An Unconditionally Stable Scheme for the Numerical Solution of Finite Difference Forms of Poisson's Equation

A procedure for an unconditionally stable non-iterative step-by-stepscheme for the solution of second order ordinary boundary valuedifference equations is extended to Poisson's equation withDirichlet boundary conditions. It is shown that owing to theparticular method of solution rounding error from a particulargrid point has practically no effect on other grid points, andit is absolutely stable for all finite difference forms of Poisson'sequation. The extension to other two dimensional partial differentialequations is discussed. It is also shown that the scheme isparticularly valuable in problems in which a high degree ofaccuracy is demanded and it is recommended for numerical weatherexperiments.     

三次多项式样条解一类四阶边值问题

应用三次多项式样条函数解一组四阶单侧、障碍、接触边值问题,取半结点为网格点,并增加了边界方程,应用方法解文献中的数值例子,说明了方法的高效性,数值结果也显示了方法的优越性.     

A parameter-uniform numerical method for a singularly perturbed two parameter elliptic problem

In this paper, a class of singularly perturbed elliptic partial differential equations posed on a rectangular domain is studied. The differential equation contains two singular perturbation parameters. The solutions of these singularly perturbed problems are decomposed into a sum of regular, boundary layer and corner layer components. Parameter-explicit bounds on the derivatives of each of these components are derived. A numerical algorithm based on an upwind finite difference operator and a tensor product of piecewise-uniform Shishkin meshes is analysed. Parameter-uniform asymptotic error bounds for the numerical approximations are established.     

Multigrid solution of high order discretisation for three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind

In this paper, we propose two compact finite difference approximations for three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind. In these methods there is no need to define special formulas near the boundaries and boundary conditions are incorporated with these techniques. The unknown solution and its second derivatives are carried as unknowns at grid points. We derive second-order and fourth-order approximations on a 27 point compact stencil. Classical iteration methods such as Gauss–Seidel and SOR for solving the linear system arising from the second-order and fourth-order discretisation suffer from slow convergence. In order to overcome this problem we use multigrid method which exhibit grid-independent convergence and solve the linear system of equations in small amount of computer time. The fourth-order finite difference approximations are used to solve several test problems and produce high accurate numerical solutions.     

一类反应扩散方程组的保持正性的差分格式

本文研究一类具有正解的反应扩散方程组的有限差分解法.构造了一个保持正性的差分格式.利用离散的最大值原理证明了差分格式解的非负性,有界性及差分格式的无条件稳定性.这些估计的证明不依赖于微分方程的解而仅仅与初边值条件有关.当微分方程的解适当光滑时,证明了差分格式的一致收敛性.最后给出了数值计算结果,并与以往方法进行了比较.计算结果说明了本文给出的方法的有效性.     

二维线性双曲型方程Neumann边值问题的紧交替方向隐格式

对二维Neumann边界条件的线性双曲型方程建立了紧交替方向的隐格式.利用方程和边界条件得到在空间上的三阶与五阶导数的边界值,进而在内点、边界内点和边界角点分别建立9点、6点和4点紧差分格式;通过引进新的范数和L₂范数估计L_∞范数;借助能量估计、Gronwall不等式和Schwarz不等式等技巧,详细分析了差分格式在无穷范数下关于时间和空间分别为二阶和四阶收敛性,并给出了稳定性结果;通过数值算例,验证了理论分析结果.     

Stochastic methods for Dirichlet problems

Using an equivalent expression for solutions of second order Dirichlet problems in terms of Ito type stochastic differential equations, we develop a numerical solution method for Dirichlet boundary value problems. It is possible with this idea to solve for solution values of a partial differential equation at isolated points without having to construct any kind of mesh and without knowing approximations for the solution at any other points. Our method is similar to a recently published approach, but differs primarily in the handling of the boundary. Some numerical examples are presented, applying these techniques to model Laplace and Poisson equations on the unit disk. Visiting Professor, Universidad de Salamanca.     

自适应迎风格式求解奇异摄动两点边值问题的高精度算法(二)

0引言 考虑与文[1]相同的奇异摄动两点边值问题的数值解法: Tu(x):=-εu″(x)-p(x)u′(x)=f(x),x∈(0,1); (1) u(0)=0,u(1)=1. (2) 其中ε是一个常数,0<ε≤1,f∈C2[0,1].假定P∈C3[0,1]且存在常数β和-β使得0<β≤p(x)≤-β,|p′(x)|≤-β,(V)x∈[0,1] (3) 成立.     

A parameter-uniform numerical method for time-dependent singularly perturbed differential-difference equations

A numerical study is made for solving a class of time-dependent singularly perturbed convection–diffusion problems with retarded terms which often arise in computational neuroscience. To approximate the retarded terms, a Taylor’s series expansion has been used and the resulting time-dependent singularly perturbed differential equation is approximated using parameter-uniform numerical methods comprised of a standard implicit finite difference scheme to discretize in the temporal direction on a uniform mesh by means of Rothe’s method and a B-spline collocation method in the spatial direction on a piecewise-uniform mesh of Shishkin type. The method is shown to be accurate of order O(M⁻¹ + N⁻² ln³N), where M and N are the number of mesh points used in the temporal direction and in the spatial direction respectively. An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. Numerical results are given to illustrate the parameter-uniform convergence of the numerical approximations. Comparisons of the numerical solutions are performed with an upwind and midpoint upwind finite difference scheme on a piecewise-uniform mesh to demonstrate the efficiency of the method.     

Efficient computational methods for singular free boundary problems using smoothing variable substitutions

We study a class of singular free boundary problems for the degenerate m$m$-Laplacian. Taking into account the behavior of the solution in the neighborhood of the singular points, a variable substitution is introduced, which makes the solution smooth in all the domain. Then, a standard finite difference scheme is used to discretize the problem. The numerical results suggest that in this way the second order convergence of the finite difference scheme is recovered, in spite of the singularities.     

Finite‐difference schemes for transport‐dominated equations using special collocation nodes

We introduce finite‐difference schemes based on a special upwind‐type collocation grid, in order to obtain approximations of the solution of linear transport‐dominated advection‐diffusion problems. The method is well suited when the diffusion parameter is very small compared to the discretization parameter. A theory is developed and many numerical experiments are shown. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Uniform Convergence Analysis of Finite Difference Scheme for Singularly Perturbed Delay Differential Equation on an Adaptively Generated Grid

<正>Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.     

Parameter-uniform finite difference schemes for singularly perturbed parabolic diffusion-convection-reaction problems

In this paper, parameter-uniform numerical methods for a class of singularly perturbed parabolic partial differential equations with two small parameters on a rectangular domain are studied. Parameter-explicit theoretical bounds on the derivatives of the solutions are derived. The solution is decomposed into a sum of regular and singular components. A numerical algorithm based on an upwind finite difference operator and an appropriate piecewise uniform mesh is constructed. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are given to illustrate the parameter-uniform convergence of the numerical approximations.

     

Uniformly convergent difference schemes for a singularly perturbed third order boundary value problem

In this paper we consider a numerical approximation of a third order singularly perturbed boundary value problem by an upwind finite difference scheme on a Shishkin mesh. The behavior of the solution, and the stability of the continuous problem are discussed. The proof of the uniform convergence of the proposed numerical method is based on the strongly uniform stability and a weak consistency property of the discrete problem. Numerical experiments verify our theoretical results.     

A new pseudospectral method with upwind features

A new pseudospectral method is presented for numerical solutionsof singular perturbation problems without turning points. Boththeoretical and numerical analyses show that this new methodis an upwind scheme. It is shown that when the perturbationparameter is fixed the computed solution converges spectrallyto the exact solution as the number of collocation points tendsto infinity. ^* Supported by a Royal Fellowship Award from the Royal Societyof London.     

l 1-error estimates on the immersed interface upwind scheme for linear convection equations with piecewise constant coefficients: A simple proof

Abstract A linear convection equation with discontinuous coefcients arises in wave propagation through interfaces.An interface condition is needed at the interface to select a unique solution.An upwind scheme that builds this interface condition into its numerical flux is called the immersed interface upwind scheme.An l1-error estimate of such a scheme was frst established by Wen et al.(2008).In this paper,we provide a simple analysis on the l1-error estimate.The main idea is to formulate the solution to the underline initial-value problem into the sum of solutions to two convection equations with constant coefcients,which can then be estimated using classical methods for the initial or boundary value problems.