Given a one-step numerical scheme,on which ordinary differential equations is it exact?

A necessary condition for a (non-autonomous) ordinary differential equation to be exactly solved by a one-step, finite difference method is that the principal term of its local truncation error be null. A procedure to determine some ordinary differential equations exactly solved by a given numerical scheme is developed. Examples of differential equations exactly solved by the explicit Euler, implicit Euler, trapezoidal rule, second-order Taylor, third-order Taylor, van Niekerk’s second-order rational, and van Niekerk’s third-order rational methods are presented.     

A second order backward difference method with variable steps for a parabolic problem

The numerical solution of a parabolic problem is studied. The equation is discretized in time by means of a second order two step backward difference method with variable time step. A stability result is proved by the energy method under certain restrictions on the ratios of successive time steps. Error estimates are derived and applications are given to homogenous equations with initial data of low regularity.     

On exponential stability of nonlinear Volterra difference equations in phase spaces

Using a novel approach, we present some new explicit criteria for global exponential stability of the zero solution of general nonlinear Volterra difference equations in phase spaces. In particular, this gives a solution to an open problem posed very recently by E. Braverman and I. M. Karabash in Journal of Difference Equations and Applications 18 , 909–939 (2012). As an application, we apply the obtained results to study asymptotic behavior of equilibriums of discrete time neural networks.     

On some explicit Adams multistep methods for fractional differential equations

In this paper we present a family of explicit formulas for the numerical solution of differential equations of fractional order. The proposed methods are obtained by modifying, in a suitable way, Fractional-Adams–Moulton methods and they represent a way for extending classical Adams–Bashforth multistep methods to the fractional case. The attention is hence focused on the investigation of stability properties. Intervals of stability for k$k$-step methods, k=1,…,5$k=1,\dots ,5$, are computed and plots of stability regions in the complex plane are presented.     

Global error estimation and control in linearly-implicit parallel two-step peer W-methods

The class of linearly-implicit parallel two-step peer W-methods has been designed recently for efficient numerical solutions of stiff ordinary differential equations. Those schemes allow for parallelism across the method, that is an important feature for implementation on modern computational devices. Most importantly, all stage values of those methods possess the same properties in terms of stability and accuracy of numerical integration. This property results in the fact that no order reduction occurs when they are applied to very stiff problems. In this paper, we develop parallel local and global error estimation schemes that allow the numerical solution to be computed for a user-supplied accuracy requirement in automatic mode. An algorithm of such global error control and other technical particulars are also discussed here. Numerical examples confirm efficiency of the presented error estimation and stepsize control algorithm on a number of test problems with known exact solutions, including nonstiff, stiff, very stiff and large-scale differential equations. A comparison with the well-known stiff solver RODAS is also shown.     

Doubly quasi-consistent parallel explicit peer methods with built-in global error estimation

Recently, Kulikov presented the idea of double quasi-consistency, which facilitates global error estimation and control, considerably. More precisely, a local error control implemented in such methods plays a part of global error control at the same time. However, Kulikov studied only Nordsieck formulas and proved that there exists no doubly quasi-consistent scheme among those methods.Here, we prove that the class of doubly quasi-consistent formulas is not empty and present the first example of such sort. This scheme belongs to the family of superconvergent explicit two-step peer methods constructed by Weiner, Schmitt, Podhaisky and Jebens. We present a sample of s-stage doubly quasi-consistent parallel explicit peer methods of order s−1 when s=3. The notion of embedded formulas is utilized to evaluate efficiently the local error of the constructed doubly quasi-consistent peer method and, hence, its global error at the same time. Numerical examples of this paper confirm clearly that the usual local error control implemented in doubly quasi-consistent numerical integration techniques is capable of producing numerical solutions for user-supplied accuracy conditions in automatic mode.     

Stability and error estimates for the $\theta$-method for strongly monotone and infinitely stiff evolution equations

Summary.   For evolution equations with a strongly monotone operator we derive unconditional stability and discretization error estimates valid for all . For the -method, with , we prove an error estimate , if , where is the maximal integration step for an arbitrary choice of sequence of steps and with no assumptions about the existence of the Jacobian as well as other derivatives of the operator , and an optimal estimate under some additional relation between neighboring steps. The first result is an improvement over the implicit midpoint method , for which an order reduction to sometimes may occur for infinitely stiff problems. Numerical tests illustrate the results. Received March 10, 1999 / Revised version received April 3, 2000 / Published online February 5, 2001     

Pseudospectra for matrix pencils and stability of equilibria

The concept of ε-pseudospectra for matrices, introduced by Trefethen and his coworkers, has been studied extensively since 1990. In this paper, ε-pseudospectra for matrix pencils, which are relevant in connection with generalized eigenvalue problems, are considered. Some properties as well as the practical computation of ε-pseudospectra for matrix pencils will be discussed. As an application, we demonstrate how this concept can be used for investigating the asymptotic stability of stationary solutions to time-dependent ordinary or partial differential equations; two cases, based on Burgers' equation, will be shown. This research has been supported by the Netherlands Organization for Scientific Research (N.W.O.)     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

Convergence and stability of the finite difference scheme for nonlinear parabolic systems with time delay

In this article, a class of nonlinear evolution equations – reaction–diffusion equations with time delay – is studied. By combining the domain decomposition technique and the finite difference method, the results for the existence, convergence and the stability of the numerical solution are obtained in the case of subdomain overlap and when the time-space is completely discretized. This revised version was published online in June 2006 with corrections to the Cover Date.     

On quasi-consistent integration by Nordsieck methods

In this paper we study advantages of numerical integration by quasi-consistent Nordsieck formulas. All quasi-consistent numerical methods possess at least one important property for practical use, which has not attracted attention yet, i.e. the global error of a quasi-consistent method has the same order as its local error. This means that the usual local error control will produce a numerical solution for the prescribed accuracy requirement if the principal term of the local error dominates strongly over remaining terms. In other words, the global error control can be as cheap as the local error control in the methods under discussion.     

Global dynamics of a competitive system of rational difference equations

In this paper, we study the qualitative behavior of a competitive system of second‐order rational difference equations. More precisely, we investigate the boundedness character, existence and uniqueness of positive equilibrium point, local asymptotic stability and global stability of the unique positive equilibrium point and rate of convergence of positive solutions of the system. Some numerical examples are given to verify our theoretical results. Copyright © 2014 John Wiley & Sons, Ltd.     

Explicit methods for fractional differential equations and their stability properties

The use of explicit methods in the numerical treatment of differential equations of fractional order is an area not yet widely investigated. In this paper stability properties of some multistep methods of explicit type are investigated and new methods with larger intervals of stability are proposed. Some numerical experiments are presented in order to validate theoretical results.     

Methods for parallel integration of stiff systems of odes

This paper presents a class of parallel numerical integration methods for stiff systems of ordinary differential equations which can be partitioned into loosely coupled sub-systems. The formulas are called decoupled backward differentiation formulas, and they are derived from the classical formulas by restricting the implicit part to the diagnonal sub-system. With one or several sub-systems allocated to each processor, information only has to be exchanged after completion of a step but not during the solution of the nonlinear algebraic equations.The main emphasis is on the formula of order 1, the decoupled implicit Euler formula. It is proved that this formula even for a wide range of multirate formulations has an asymptotic global error expansion permitting extrapolation. Besides, sufficient conditions for absolute stability are presented.     

Multirate Runge–Kutta schemes for advection equations

Explicit time integration methods can be employed to simulate a broad spectrum of physical phenomena. The wide range of scales encountered lead to the problem that the fastest cell of the simulation dictates the global time step. Multirate time integration methods can be employed to alter the time step locally so that slower components take longer and fewer time steps, resulting in a moderate to substantial reduction of the computational cost, depending on the scenario to simulate [S. Osher, R. Sanders, Numerical approximations to nonlinear conservation laws with locally varying time and space grids, Math. Comput. 41 (1983) 321–336; H. Tang, G. Warnecke, A class of high resolution schemes for hyperbolic conservation laws and convection-diffusion equations with varying time and pace grids, SIAM J. Sci. Comput. 26 (4) (2005) 1415–1431; E. Constantinescu, A. Sandu, Multirate timestepping methods for hyperbolic conservation laws, SIAM J. Sci. Comput. 33 (3) (2007) 239–278]. In air pollution modeling the advection part is usually integrated explicitly in time, where the time step is constrained by a locally varying Courant–Friedrichs–Lewy (CFL) number. Multirate schemes are a useful tool to decouple different physical regions so that this constraint becomes a local instead of a global restriction. Therefore it is of major interest to apply multirate schemes to the advection equation. We introduce a generic recursive multirate Runge–Kutta scheme that can be easily adapted to an arbitrary number of refinement levels. It preserves the linear invariants of the system and is of third order accuracy when applied to certain explicit Runge–Kutta methods as base method.     

Global asymptotic stability beyond 3/2 type stability for a logistic equation with piecewise constant arguments

In this paper, a logistic equation with multiple piecewise constant arguments is investigated in detail. We generalize the approach in two papers, [K. Uesugi, Y. Muroya, E. Ishiwata, On the global attractivity for a logistic equation with piecewise constant arguments, J. Math. Anal. Appl. 294 (2) (2004) 560-580] and [Y. Muroya, E. Ishiwata, N. Guglielmi, Global stability for nonlinear difference equations with variable coefficients, J. Math. Anal. Appl. 334 (1) (2007) 232-247], and establish a new condition for the global stability of the equation. Their results are given as one of the special cases. Moreover, we improve the 3/2 type stability condition under several dominance assumptions on the coefficients of the equation. Some examples and numerical simulations are also presented. All of these examples show that there are several conditions for the global stability of the equation, depending on the coefficients on the delay terms of the equation, beyond the 3/2 type stability condition.     

Sufficient conditions for the global stability of nonautonomous higher order difference equations

We present some explicit sufficient conditions for the global stability of the zero solution in nonautonomous higher order difference equations. The linear case is discussed in detail. We illustrate our main results with some examples. In particular, the stability properties of the equilibrium in a nonlinear model in macroeconomics is addressed.     

Stability of Runge-Kutta methods for the generalized pantograph equation

Summary. This paper deals with stability properties of Runge-Kutta (RK) methods applied to a non-autonomous delay differential equation (DDE) with a constant delay which is obtained from the so-called generalized pantograph equation, an autonomous DDE with a variable delay by a change of the independent variable. It is shown that in the case where the RK matrix is regular stability properties of the RK method for the DDE are derived from those for a difference equation, which are examined by similar techniques to those in the case of autonomous DDEs with a constant delay. As a result, it is shown that some RK methods based on classical quadrature have a superior stability property with respect to the generalized pantograph equation. Stability of algebraically stable natural RK methods is also considered. Received May 5, 1998 / Revised version received November 17, 1998 / Published online September 24, 1999     

Defect-based local error estimators for splitting methods,with application to Schrödinger equations,Part I: The linear case

We introduce a defect correction principle for exponential operator splitting methods applied to time-dependent linear Schrödinger equations and construct a posteriori local error estimators for the Lie–Trotter and Strang splitting methods. Under natural commutator bounds on the involved operators we prove asymptotical correctness of the local error estimators, and along the way recover the known a priori convergence bounds. Numerical examples illustrate the theoretical local and global error estimates.     

Beyond conventional Runge–Kutta methods in numerical integration of ODEs and DAEs by use of structures and local models

There are two parts in this paper. In the first part we consider an overdetermined system of differential-algebraic equations (DAEs). We are particularly concerned with Hamiltonian and Lagrangian systems with holonomic constraints. The main motivation is in finding methods based on Gauss coefficients, preserving not only the constraints, symmetry, symplecticness, and variational nature of trajectories of holonomically constrained Hamiltonian and Lagrangian systems, but also having optimal order of convergence. The new class of (s,s)$(s,s)$-Gauss–Lobatto specialized partitioned additive Runge–Kutta (SPARK) methods uses greatly the structure of the DAEs and possesses all desired properties. In the second part we propose a unified approach for the solution of ordinary differential equations (ODEs) mixing analytical solutions and numerical approximations. The basic idea is to consider local models which can be solved efficiently, for example analytically, and to incorporate their solution into a global procedure based on standard numerical integration methods for the correction. In order to preserve also symmetry we define the new class of symmetrized Runge–Kutta methods with local model (SRKLM).