Convergence of rational multistep methods of Adams-Padé type

Rational generalizations of multistep schemes, where the linear stiff part of a given problem is treated by an A-stable rational approximation, have been proposed by several authors, but a reasonable convergence analysis for stiff problems has not been provided so far. In this paper we directly relate this approach to exponential multistep methods, a subclass of the increasingly popular class of exponential integrators. This natural, but new interpretation of rational multistep methods enables us to prove a convergence result of the same quality as for the exponential version. In particular, we consider schemes of rational Adams type based on A-acceptable Padé approximations to the matrix exponential. A numerical example is also provided.     

非自治刚性随机微分方程正则EM分裂方法的收敛性和稳定性

本文研究了数值求解非自治随机微分方程的正则Euler-Maruyama分裂（CEMS）方法,该方程的漂移项系数带有刚性且允许超线性增长,扩散项系数满足全局Lipschitz条件.首先,证明了CEMS方法的强收敛性及收敛速度.其次,证明了在适当条件下CEMS方法是均方稳定的.进一步,利用离散半鞅收敛定理,研究了CEMS方法的几乎必然指数稳定性.结果表明,CEMS方法在漂移系数的刚性部分满足单边Lipschitz条件下可保持几乎必然指数稳定性.最后通过数值实验,检验了CEMS方法的有效性并证实了我们的理论结果.     

Implicit Euler and Lie splitting discretizations of nonlinear parabolic equations with delay

A convergence analysis is presented for the implicit Euler and Lie splitting schemes when applied to nonlinear parabolic equations with delay. More precisely, we consider a vector field which is the sum of an unbounded dissipative operator and a delay term, where both point delays and distributed delays fit into the framework. Such equations are frequently encountered, e.g. in population dynamics. The main theoretical result is that both schemes converge with an order (of at least) \(q=1/2\) , without any artificial regularity assumptions. We discuss implementation details for the methods, and the convergence results are verified by numerical experiments demonstrating both the correct order, as well as the efficiency gain of Lie splitting as compared to the implicit Euler scheme.     

High order splitting methods for analytic semigroups exist

In this paper, we are concerned with the construction and analysis of high order exponential splitting methods for the time integration of abstract evolution equations which are evolved by analytic semigroups. We derive a new class of splitting methods of orders three to fourteen based on complex coefficients. An optimal convergence analysis is presented for the methods when applied to equations on Banach spaces with unbounded vector fields. These results resolve the open question whether there exist splitting schemes with convergence rates greater then two in the context of semigroups. As a concrete application we consider parabolic equations and their dimension splittings. The sharpness of our theoretical error bounds is further illustrated by numerical experiments.     

刚性多滞量积分微分方程的Runge—Kutta方法

本文研究了求解刚性多滞量积分微分方程的Runge-Kutta方法的非线性稳定性和计算有效性.经典Runge—Kutta方法连同复合求积公式和Pouzet求积公式被改造用于求解一类刚性多滞量Volterra型积分微分方程.其分析导出了:在适当条件下,扩展的Runge-Kutta方法是渐近稳定和整体稳定的.此外,数值试验表明所给出的方法是高度有效的.     

Convergence analysis of high-order time-splitting pseudo-spectral methods for rotational Gross–Pitaevskii equations

A convergence analysis of time-splitting pseudo-spectral methods adapted for time-dependent Gross–Pitaevskii equations with additional rotation term is given. For the time integration high-order exponential operator splitting methods are studied, and the space discretization relies on the generalized-Laguerre–Fourier spectral method with respect to the $(x,y)$ -variables as well as the Hermite spectral method in the $z$ -direction. Essential ingredients in the stability and error analysis are a general functional analytic framework of abstract nonlinear evolution equations, fractional power spaces defined by the principal linear part, a Sobolev-type inequality in a curved rectangle, and results on the asymptotical distribution of the nodes and weights associated with Gauß–Laguerre quadrature. The obtained global error estimate ensures that the nonstiff convergence order of the time integrator and the spectral accuracy of the spatial discretization are retained, provided that the problem data satisfy suitable regularity requirements. A numerical example confirms the theoretical convergence estimate.     

PARALLEL COMPOUND METHODS FOR SOLVING PARTITIONED STIFF SYSTEMS

1. IntroductionMany stiff systems occurring in practice have a special structu-re, and they can be split illtotwo coupied subsystemsyb(x) = fS(x, YS(x), YN(X)), YS(X0) = YS0, YS E R"s(1)yk(X) = fN(X, YS(X), YN(X)), YN(X0) = YN0, YS 6 R"Nwhere ys denotes the vector of sto f comPonents and yN denotes the vector of nonstiff componeots. FOr such partitioned systems a partitioned discretization method is used, i.e. thestAf subsystem is sothed by a "stiff" method and the nonstffi by a…     

Low Regularity Exponential-Type Integrators for Semilinear Schrödinger Equations

We introduce low regularity exponential-type integrators for nonlinear Schrödinger equations for which first-order convergence only requires the boundedness of one additional derivative of the solution. More precisely, we will prove first-order convergence in \(H^r\) for solutions in \(H^{r+1}\) (with \(r > d/2\)) of the derived schemes. This allows us lower regularity assumptions on the data than for instance required for classical splitting or exponential integration schemes. For one-dimensional quadratic Schrödinger equations, we can even prove first-order convergence without any loss of regularity. Numerical experiments underline the favorable error behavior of the newly introduced exponential-type integrators for low regularity solutions compared to classical splitting and exponential integration schemes.     

求解带刚性源项标量双曲型守恒律方程的保有界WCNS格式

带刚性源项的双曲守恒律方程是很多物理问题,特别是化学反应流的数学模型.本文考虑带刚性源项的标量双曲型守恒律方程,通过时空分离的方式,发展了一类保有界的WCNS格式.对于空间离散,我们将参数化的通量限制器推广到WCNS框架,使得方程对流项离散后满足极值原理.对于时间离散,我们将半离散的WCNS改写成指数形式,采用三阶修正指数型Runge-Kutta格式来控制方程的刚性,保持数值解的界.可以证明,本文格式对带刚性源项的一维标量守恒律方程具有保有界性和弱渐近保持性.数值试验验证了方法的有效性.     

On interpolation variants of Newton’s method for functions of several variables

A generalization of the variants of Newton’s method based on interpolation rules of quadrature is obtained, in order to solve systems of nonlinear equations. Under certain conditions, convergence order is proved to be 2d+1, where d is the order of the partial derivatives needed to be zero in the solution. Moreover, different numerical tests confirm the theoretical results and allow us to compare these variants with Newton’s classical method, whose convergence order is d+1 under the same conditions.     

A Class of Explicit Exponential General Linear Methods

In this paper, we consider a class of explicit exponential integrators that includes as special cases the explicit exponential Runge–Kutta and exponential Adams–Bashforth methods. The additional freedom in the choice of the numerical schemes allows, in an easy manner, the construction of methods of arbitrarily high order with good stability properties. We provide a convergence analysis for abstract evolution equations in Banach spaces including semilinear parabolic initial-boundary value problems and spatial discretizations thereof. From this analysis, we deduce order conditions which in turn form the basis for the construction of new schemes. Our convergence results are illustrated by numerical examples. AMS subject classification (2000) 65L05, 65L06, 65M12, 65J10     

Smoothing schemes for reaction-diffusion systems with nonsmooth data

Cox and Matthews [S.M. Cox, P.C. Matthews, Exponential time differencing for stiff systems, J. Comput. Phys. 176 (2002) 430–455] developed a class of Exponential Time Differencing Runge–Kutta schemes (ETDRK) for nonlinear parabolic equations; Kassam and Trefethen [A.K. Kassam, Ll. N. Trefethen, Fourth-order time stepping for stiff pdes, SIAM J. Sci. Comput. 26 (2005) 1214–1233] have shown that these schemes can suffer from numerical instability and they proposed a modified form of the fourth-order (ETDRK4) scheme. They use complex contour integration to implement these schemes in a way that avoids inaccuracies when inverting matrix polynomials, but this approach creates new difficulties in choosing and evaluating the contour for larger problems. Neither treatment addresses problems with nonsmooth data, where spurious oscillations can swamp the numerical approximations if one does not treat the problem carefully. Such problems with irregular initial data or mismatched initial and boundary conditions are important in various applications, including computational chemistry and financial engineering. We introduce a new version of the fourth-order Cox–Matthews, Kassam–Trefethen ETDRK4 scheme designed to eliminate the remaining computational difficulties. This new scheme utilizes an exponential time differencing Runge–Kutta ETDRK scheme using a diagonal Padé approximation of matrix exponential functions, while to deal with the problem of nonsmooth data we use several steps of an ETDRK scheme using a sub-diagonal Padé formula. The new algorithm improves computational efficiency with respect to evaluation of the high degree polynomial functions of matrices, having an advantage of splitting the matrix polynomial inversion problem into a sum of linear problems that can be solved in parallel. In this approach it is only required that several backward Euler linear problems be solved, in serial or parallel. Numerical experiments are described to support the new scheme.     

Stable strong order 1.0 schemes for solving stochastic ordinary differential equations

The Balanced method was introduced as a class of quasi-implicit methods, based upon the Euler-Maruyama scheme, for solving stiff stochastic differential equations. We extend the Balanced method to introduce a class of stable strong order 1.0 numerical schemes for solving stochastic ordinary differential equations. We derive convergence results for this class of numerical schemes. We illustrate the asymptotic stability of this class of schemes is illustrated and is compared with contemporary schemes of strong order 1.0. We present some evidence on parametric selection with respect to minimising the error convergence terms. Furthermore we provide a convergence result for general Balanced style schemes of higher orders.     

Radial basis collocation method and quasi‐Newton iteration for nonlinear elliptic problems

This work presents a radial basis collocation method combined with the quasi‐Newton iteration method for solving semilinear elliptic partial differential equations. The main result in this study is that there exists an exponential convergence rate in the radial basis collocation discretization and a superlinear convergence rate in the quasi‐Newton iteration of the nonlinear partial differential equations. In this work, the numerical error associated with the employed quadrature rule is considered. It is shown that the errors in Sobolev norms for linear elliptic partial differential equations using radial basis collocation method are bounded by the truncation error of the RBF. The combined errors due to radial basis approximation, quadrature rules, and quasi‐Newton and Newton iterations are also presented. This result can be extended to finite element or finite difference method combined with any iteration methods discussed in this work. The numerical example demonstrates a good agreement between numerical results and analytical predictions. The numerical results also show that although the convergence rate of order 1.62 of the quasi‐Newton iteration scheme is slightly slower than rate of order 2 in the Newton iteration scheme, the former is more stable and less sensitive to the initial guess. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Development and Application of an Exponential Method for Integrating Stiff Systems Based on the Classical Runge–Kutta Method

We study numerical methods for solving stiff systems of ordinary differential equations. We propose an exponential computational algorithm which is constructed by using an exponential change of variables based on the classical Runge–Kutta method of the fourth order. Nonlinear problems are used to prove and demonstrate the fourth order of convergence of the new method.     

Convergence in variation of solutions of nonlinear Fokker–Planck–Kolmogorov equations to stationary measures

We study convergence in variation of probability solutions of nonlinear Fokker–Planck–Kolmogorov equations to stationary solutions. We obtain sufficient conditions for the exponential convergence of solutions to the stationary solution in case of coefficients that can have an arbitrary growth at infinity and depend on the solutions through convolutions with unbounded discontinuous kernels. In addition, we study a more difficult case where the nonlinear equation has several stationary solutions and convergence to a stationary solution depends on initial data. Finally, we obtain sufficient conditions for solvability of nonlinear Fokker–Planck–Kolmogorov equations.     

On convergence of numerical schemes for hyperbolic conservation laws with stiff source terms

We deal in this study with the convergence of a class of numerical schemes for scalar conservation laws including stiff source terms. We suppose that the source term is dissipative but it is not necessarily a Lipschitzian function. The convergence of the approximate solution towards the entropy solution is established for first and second order accurate MUSCL and for splitting semi-implicit methods.

     

Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules

The Gauss product quadrature rules and collocation method are applied to reduce the second-kind nonlinear two-dimensional Fredholm integral equations (FIE) to a nonlinear system of equations. The convergence of the proposed numerical method is proved under certain conditions on the kernel of the integral equation. An iterative method for approximating the solution of the obtained nonlinear system is provided and its convergence is proved. Also, some numerical examples are presented to show the efficiency and accuracy of the proposed method.     

Exponential Time Differencing-Padé Finite Element Method for Nonlinear Convection-Diffusion-Reaction Equations with Time Constant Delay

In this paper, ETD3-Padé and ETD4-Padé Galerkin finite element methods are proposed and analyzed for nonlinear delayed convection-diffusion-reaction equations with Dirichlet boundary conditions. An ETD-based RK is used for time integration of the corresponding equation. To overcome a well-known difficulty of numerical instability associated with the computation of the exponential operator, the Padé approach is used for such an exponential operator approximation, which in turn leads to the corresponding ETD-Padé schemes. An unconditional $L^2$ numerical stability is proved for the proposed numerical schemes, under a global Lipshitz continuity assumption. In addition, optimal rate error estimates are provided, which gives the convergence order of $O(k^{3}+h^{r})$ (ETD3-Padé) or $O(k^{4}+h^{r})$ (ETD4-Padé) in the $L^2$ norm, respectively. Numerical experiments are presented to demonstrate the robustness of the proposed numerical schemes.     

非线性抛物型方程有限元法数值积分的有效性

Abstract. The effect of numerical integration in finite element methods applied to a class of nonlinear parabolic equations is considered and some sufficient conditions on the quadrature scheme to ensure that the order of convergence is unaltered in the presence of numerical integration are given. Optimal Lz and H1 estimates for the error and its time derivative are established.