Order conditions for partitioned Runge-Kutta methods

We illustrate the use of the recent approach by P. Albrecht to the derivation of order conditions for partitioned Runge-Kutta methods for ordinary differential equations.     

Quadratic invariants and multi-symplecticity of partitioned Runge-Kutta methods for Hamiltonian PDEs

In this paper, we study the preservation of quadratic conservation laws of Runge-Kutta methods and partitioned Runge-Kutta methods for Hamiltonian PDEs and establish the relation between multi-symplecticity of Runge-Kutta method and its quadratic conservation laws. For Schrödinger equations and Dirac equations, it reveals that multi-symplectic Runge-Kutta methods applied to equations with appropriate boundary conditions can preserve the global norm conservation and the global charge conservation, respectively.     

Hamilton系统MSRK算法的数值色散关系

Numerical dispersion relation of the multi-symplectic Runge-Kutta (MSRK) method for linear Hamiltonian PDEs is derived in the present paper, which is shown to be a discrete counterpart to that possessed by the differential equation. This provides further understanding of MSRK methods. However, much still remains to be investigated further.     

Convergence of Runge-Kutta methods for neutral Volterra delay-integro-differential equations

In this paper, we focus on the error behavior of Runge-Kutta methods for nonlinear neutral Volterra delay-integro-differential equations (NVDIDEs) with constant delay. The convergence properties of the Runge-Kutta methods with two classes of quadrature technique, compound quadrature rule and Pouzet type quadrature technique, are investigated.      

The multi-symplectic Fourier pseudospectral method for solving two-dimensional Hamiltonian PDEs

In this paper, the multi-symplectic Fourier pseudospectral (MSFP) method is generalized to solve two-dimensional Hamiltonian PDEs with periodic boundary conditions. Using the Fourier pseudospectral method in the space of the two-dimensional Hamiltonian PDE (2D-HPDE), the semi-discrete system obtained is proved to have semi-discrete multi-symplectic conservation laws and a global symplecticity conservation law. Then, the implicit midpoint rule is employed for time integration to obtain the MSFP method for the 2D-HPDE. The fully discrete multi-symplectic conservation laws are also obtained. In addition, the proposed method is applied to solve the Zakharov-Kuznetsov (ZK) equation and the Kadomtsev-Petviashvili (KP) equation. Numerical experiments on soliton solutions of the ZK equation and the KP equation show the high accuracy and effectiveness of the proposed method.     

Convergence Aspects of Step-Parallel Iteration of Runge-Kutta Methods for Delay Differential Equations

Implicit Runge-Kutta methods are known as highly accurate and stable methods for solving differential equations. However, the iteration technique used to solve implicit Runge-Kutta methods requires a lot of computational efforts. To lessen the computational effort, one can iterate simultaneously at a number of points along the t-axis. In this paper, we extend the PDIRK (Parallel Diagonal Iterated Runge-Kutta) methods to delay differential equations (DDEs). We give the region of convergence and analyze the speed of convergence in three parts for the P-stability region of the Runge-Kutta corrector. It is proved that PDIRK methods to DDEs are efficient, and the diagonal matrix D of the PDIRK methods for DDES can be selected in the same way as for ordinary differential equations (ODEs).     

The stability of modified Runge-Kutta methods for the pantograph equation

In the present paper, the modified Runge-Kutta method is constructed, and it is proved that the modified Runge-Kutta method preserves the order of accuracy of the original one. The necessary and sufficient conditions under which the modified Runge-Kutta methods with the variable mesh are asymptotically stable are given. As a result, the -methods with , the odd stage Gauss-Legendre methods and the even stage Lobatto IIIA and IIIB methods are asymptotically stable. Some experiments are given.

     

On higher-order semi-explicit symplectic partitioned Runge-Kutta methods for constrained Hamiltonian systems

Summary. In this paper we generalize the class of explicit partitioned Runge-Kutta (PRK) methods for separable Hamiltonian systems to systems with holonomic constraints. For a convenient analysis of such schemes, we first generalize the backward error analysis for systems in to systems on manifolds embedded in . By applying this analysis to constrained PRK methods, we prove that such methods will, in general, suffer from order reduction as well-known for higher-index differential-algebraic equations. However, this order reduction can be avoided by a proper modification of the standard PRK methods. This modification increases the number of projection steps onto the constraint manifold but leaves the number of force evaluations constant. We also give a numerical comparison of several second, fourth, and sixth order methods. Received May 5, 1995 / Revised version received February 7, 1996     

Three-stage Stiffly Accurate Runge-Kutta Methods for Stiff Stochastic Differential Equations

In this paper we discuss diagonally implicit and semi-implicit methods based on the three-stage stiffly accurate Runge-Kutta methods for solving Stratonovich stochastic differential equations（SDEs）.Two methods,a three-stage stiffly accurate semi-implicit（SASI3） method and a three-stage stiffly accurate diagonally implicit （SADI3） method,are constructed in this paper.In particular,the truncated random variable is used in the implicit method.The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of stiff SDEs.     

Maximal order for second derivative general linear methods with Runge-Kutta stability

An extension of general linear methods (GLMs), so-called SGLMs (GLMs with second derivative), was introduced to the case in which second derivatives, as well as first derivatives, can be calculated. SGLMs are divided into four types, depending on the nature of the differential system to be solved and the computer architecture that is used to implement these methods. In this paper, we obtain maximal order for two types of SGLMs with Runge-Kutta stability (RKS) property. Also, we construct methods of these types which possess RKS property and A-stability. Efficiency of the constructed methods is shown by numerical experiments.     

An improvement on the positivity results for 2-stage explicit Runge-Kutta methods

In this paper, we investigate the positivity property for a class of 2-stage explicit Runge-Kutta (RK2) methods of order two when applied to the numerical solution of special nonlinear initial value problems (IVPs) for ordinary differential equations (ODEs). We also pay particular attention to monotonicity property. We obtain new results for positivity which are important in practical applications. We provide some numerical examples to illustrate our results.     

Dissipativity of Runge-Kutta methods for Volterra functional differential equations

This paper is concerned with the numerical dissipativity of nonlinear Volterra functional differential equations (VFDEs). We give some dissipativity results of Runge-Kutta methods when they are applied to VFDEs. These results provide unified theoretical foundation for the numerical dissipativity analysis of systems in ordinary differential equations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs), Volterra delay integro-differential equations (VDIDEs) and VFDEs of other type which appear in practice. Numerical examples are given to confirm our theoretical results.     

On Runge-Kutta-type methods for two-dimensional stochastic differential equations

In the multidimensional case, second-order weak Runge-Kutta methods for stochastic differential equation (SDE) need simulation of correlated random variables, unless the diffusion matrix of SDE satisfies the commutativity condition. In this paper, we show that this can be avoided for some types of diffusion matrices and test functions important for applications. Published in Lietuvos Matematikos Rinkinys, Vol. 46, No. 3, pp. 403–412, July–September, 2006.     

Parallel continuous Runge-Kutta methods and vanishing lag delay differential equations

We present an explicit Runge-Kutta scheme devised for the numerical solution ofdelay differential equations (DDEs) where a delayed argument lies in the current Runge-Kutta interval. This can occur when the lag is small relative to the stepsize, and the more obvious extensions of the explicit Runge-Kutta method produce implicit equations. It transpires that the scheme is suitable forparallel implementation for solving both ODEs and more general DDEs. We associate our method with a Runge-Kutta tableau, from which the order of the method can be determined. Stability will affect the usefulness of the scheme and we derive the stability equations of the scheme when applied to the constant-coefficient test DDEu(t)=u(t) +u(t –), where the lag and the Runge-Kutta stepsizeH _n H are both constant. (The case=0 is treated separately.) In the case that 0, we consider the two distinct possibilities: (i) H and (ii)<H.In memory of Professor Leslie Fox, Balliol College, OxfordWork performed in part at The University of Auckland, New Zealand.This paper is presented as an outcome of the LMS Durham Symposium convened by Professor C.T.H. Baker on 4th–14th July 1992 with support from the SERC under Grant reference number GR/H03964.     

基于龙格—库塔方法的一阶微分方程组的初值问题

本文在R unge-K u tta方法的基础上,讨论了一阶微分方程组当其初始状态具有模糊不确定性时,运用模糊仿真的近似推理规则,求其数值解的方法.     

B-Theory of Runge-Kutta methods for stiff Volterra functional differential equations

B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations(VFDEs)are established which provide unified theoretical foundation for the studyof Runge-Kutta methods when applied to nonlinear stiff initial value problems(IVPs)in ordinary differentialequations(ODEs),delay differential equations(DDEs),integro-differential equatioons(IDEs)and VFDEs of     

Rational Runge-Kutta methods are not suitable for stiff systems of ODEs

The present paper shows that rational RK-methods are not very appropriate to solve stiff differential equations. The CA₀-stability (i.e. componentwise contractivity) is defined and the non-existence of CA₀-stable rational RK-methods is demonstrated. Furthermore it is shown that the stepsizes which can be expected when solving a stiff differential system with a rational or with an explicit linear RK-method are of the same order of magnitude.     

An $h$-Adaptive Runge-Kutta Discontinuous Galerkin Method for Hamilton-Jacobi Equations

In [35, 36], we presented an $h$-adaptive Runge-Kutta discontinuous Galerkin method using troubled-cell indicators for solving hyperbolic conservation laws. A tree data structure (binary tree in one dimension and quadtree in two dimensions) is used to aid storage and neighbor finding. Mesh adaptation is achieved by refining the troubled cells and coarsening the untroubled "children". Extensive numerical tests indicate that the proposed $h$-adaptive method is capable of saving the computational cost and enhancing the resolution near the discontinuities. In this paper, we apply this $h$-adaptive method to solve Hamilton-Jacobi equations, with an objective of enhancing the resolution near the discontinuities of the solution derivatives. One- and two-dimensional numerical examples are shown to illustrate the capability of the method.     

Preservation of Oscillations in the Runge-Kutta Method for a Type of Advanced Differential Equation

This article is concerned with the preservation of oscillations for differential equations with piecewise constant arguments of advanced type. By using the Runge-Kutta method, new oscillation conditions for numerical solution are established. We prove that oscillations of the analytic solution are preserved by the numerical solution in the Runge-Kutta method under some conditions. Some experiments are given.     

A stability property of A-stable collocation-based Runge-Kutta methods for neutral delay differential equations

We consider a linear homogeneous system of neutral delay differential equations with a constant delay whose zero solution is asymptotically stable independent of the value of the delay, and discuss the stability of collocation-based Runge-Kutta methods for the system. We show that anA-stable method preserves the asymptotic stability of the analytical solutions of the system whenever a constant step-size of a special form is used.