- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 679075, 13 pages
Nonlinear Stability and Convergence of Two-Step Runge-Kutta Methods for Volterra Delay Integro-Differential Equations
1Department of Mathematics, Heilongjiang Institute of Technology, Harbin 150050, China
2School of Management, Harbin University of Commerce, Harbin 150028, China
Received 29 January 2013; Revised 11 March 2013; Accepted 15 March 2013
Academic Editor: Changsen Yang
Copyright © 2013 Haiyan Yuan and Cheng Song. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This paper introduces the stability and convergence of two-step Runge-Kutta methods with compound quadrature formula for solving nonlinear Volterra delay integro-differential equations. First, the definitions of -algebraically stable and asymptotically stable are introduced; then the asymptotical stability of a -algebraically stable two-step Runge-Kutta method with is proved. For the convergence, the concepts of -convergence, diagonally stable, and generalized stage order are firstly introduced; then it is proved by some theorems that if a two-step Runge-Kutta method is algebraically stable and diagonally stable and its generalized stage order is , then the method with compound quadrature formula is -convergent of order at least , where depends on the compound quadrature formula.
Volterra delay integro-differential equations (VDIDEs) arise widely in the mathematical modeling of physical and biological phenomena. Significant advances in the theoretical analysis and in the numerical analysis for these problems have been made in the last few decades (see, e.g., [1, 2]). For the case of linear stability and convergence for these equations, numerical time-integration techniques of one-step collocation and Runge-Kutta type were investigated in [3–8]. Linear multistep-based methods were studied in [9–12]. De la Sen and Luo studied the uniform exponential stability of a wide class of linear time-delay systems in ; De la Sen considered the stability of impulsive time-varying systems in .
For the case of nonlinear stability and convergence, stability results were obtained in [15, 16], where the authors investigated the nonlinear stability of continuous Runge-Kutta methods, discrete Runge-Kutta methods, and backward differentiation (BDF) methods, respectively. However, most of these important results are based on the classical Lipschitz conditions, while the classical Lipschitz conditions are so strong that there are few equations satisfying them. Most of the Volterra delay integro-differential equations satisfy the one-sided Lipschitz condition, while the studies focusing on the stability and convergence of the numerical method for nonlinear VDIDEs based on a one-sided Lipschitz condition have not yet been seen in the literature. By means of a one-sided Lipschitz condition, we will discuss the stability and convergence of two-step Runge-Kutta (TSRK) methods for nonlinear VDIDEs in the present paper.
The paper is organized as follows. In Section 2, a fairly general class of VDIDEs is defined. We present a stability criterion for such problems, which generalizes the criteria in the above references. A class of two-step Runge-Kutta methods is also derived for solving VDIDEs. They are obtained by compound quadrature rules. In Sections 3 and 4, nonlinear stability and convergence of TSRK method for NDDEs are derived and proved. In Section 5 we present some numerical examples in order to illustrate the nonlinear stability and convergence of a two-step Runge-Kutta method. These numerical results show that the new methods are quite effective.
2. A Class of VDIDES and the Two-Step Runge-Kutta Methods
It is the purpose of this paper to investigate the nonlinear stability and convergence properties of the following initial-value problem VDIDEs: where , , and are smooth enough such that (1) has a unique solution and is a positive delay term. We assume that there exist some inner product and the induced norm in , such that and satisfy the following conditions: for , , and for all , and are nonnegative constants.
Throughout this paper, we assume that (1) has unique solution and denote the problem class that consist problem of type (1) with (2)–(4). In order to make the error analysis feasible, we always assume that (1) has a unique solution which is sufficiently differentiable and satisfies For function , we make the following assumptions, unless otherwise stated; all its partial derivatives used later exist and satisfy the following: Many numerical methods have been proposed for the numerical solution of (1).
In this paper, we are concerned with two-step Runge-Kutta (TSRK) methods of the formwhere , , is a stepsize, is an arbitrarily given positive integer, and . The above methods are studied in . Now we consider the adaptation of the two-step Runge-Kutta method to (1):In particular, , is the numerical approximation at to the analytic solution , the argument denotes an approximation to , and the argument denotes an approximation to which are obtained by a convergent compound quadrature (CQ) formula: using values with , .
The class of Runge-Kutta methods with CQ formula has been applied to delay-integro-differential equations by many authors (c.f. [18, 19]). For the CQ formula (9), we usually adopt the repeated trapezoidal rule, the repeated Simpson’s rule, or the repeated Newton-cotes rule, and so forth, denote . It should be pointed out that the adopted quadrature formula (9) is only a class of quadrature formula for , there also exist some other types of quadrature formula, such as Pouzet quadrature (PQ) formula and the quadrature formula based on Laguerre-Radau interpolations [20, 21]. It is the aim of our future research to investigate the adaptation of PQ formula and the quadrature formula based on Laguerre-Radau interpolations to VDIDEs.
3. The Nonlinear Stability Analysis
In this section, we will investigate the stability of the two-step Runge-Kutta methods for VDIDEs. In order to consider the stability property, we also need to consider the perturbed problem of (1): where is a given continuous function. The unique exact solution of the problem (10) is denoted as .
Applying the two-step Runge-Kutta method (7a)–(7c) to (10) leads towhere and denote approximations to and , respectively; the argument denotes an approximation to which are obtained by a convergent compound quadrature (CQ) formula: and denotes an approximation to , and for ,
3.1. Some Concepts
For the stability analysis, we need the compound quadrature formula (9) that satisfies the following condition: where is a positive constant.
Let be the internal stages and the external vectors and . Then we have a -stage partitioned general linear method:where , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and , , , .
We can identify the coefficient matrices
Definition 1 (see ). Let be real constants, a TSRK method is said to be algebraically stable if there exists a diagonal matrix and a diagonal nonnegative matrix such that is nonnegative, whereIn particular, a -algebraically stable method is called algebraically stable.
3.2. Numerical Stability of the Methods
Numerical stability is an important feature of an effective numerical method. An unstable numerical method may be consistent of high order, yet arbitrarily small perturbations will eventually cause large deviations from the true solution. In this section, we will focus on the asymptotic stability of the TSRK method.
Theorem 3. Assume that the TSRK method (7a)–(7c) is -algebraically stable with , suppose that the quadrature formula (9) satisfies the condition (14) and the conditions (2)–(4) hold, then, method (17a) and (17b) satisfies the following: when the following condition holds: where depends only on the method, the parameters , , , , , and .
Proof. It follows from a fairly straightforward (but tedious) computation and algebraically stability that (compare also [22, 23])
By means of algebraical stability of the method, we have the following:
It follows from (2)–(4) that
Substituting (25) into (24), using (14) we get the following:
By induction, we get the following:
Hence, , where .
The proof of Theorem 3 is completed.
Theorem 4. Assume that a TSRK method (7a)–(7c) is -algebraically stable with , then the TSRK method (7a)–(7c) with (9), (12) and (14) is asymptotically stable for the problem (1) with (2)–(4), when the following condition holds:
Then when , we have and .
The application of Theorem 3 yields the following: By induction, we get the following: The inequality together with the knowledge leads to the following: Because , we can get from (33) that The proof of the theorem is completed.
4. The Convergence of TSRK Method for NDDEs
4.1. Some Concepts
Definition 5. Method (7a)–(7c) with an approximation procedure (9) is said to be convergent of order if the global error satisfies a bound of the form where and are defined by where and depend on , , and .
Remark 7. The concepts of algebraic stability and diagonal stability of TSRK method are the generalizations of corresponding concepts of Runge-Kutta methods. Although it is difficult to examine these conditions, many results have been found, especially, there exist algebraically stable and diagonally stable multistep formulas of arbitrarily high order (cf. ).
Definition 8. TSRK Method (7a)–(7c) is said to have generalized stage order if is the largest integer which possesses the following properties.
For any given problem (1) and , there exists an abstract function , such that where the maximum stepsize and the constant depend only on the method and the bounds , and , they are defined by the following equations: The function is defined by the following: Particularly, when , generalized stage order is called stage order.
4.2. -Convergence and Proofs
In this section, we focus on the error analysis of TSRK method for (1). For the sake of simplicity, we always assume that all constants , , , and are dependent on the method, the bounds , , the parameters , and .
First, we give a preliminary result which will later be used several times. To simplify, we denote for , and for , , where . Define and by the following:
Proof. Since the method (7a)–(7c) is diagonally stable, there exists a positive definite diagonal matrix such that the matrix is positive definite. Therefore, the matrix is obviously nonsingular and there exists a depends only on the method such that the matrix
is also positive definite.
DefinethenUsing (2)–(4), (44), (46a), and (46b), we have for , where is the minimum eigenvalue of . Therefore, whereFrom (43a), (49a), and (49b), it follows that where , which completes the proof of Theorem 9.
Proof. Define , we get from (45a)–(45d) that With algebraic stability, the matrix is nonnegative definite. As in , we have Using (2)–(4), we further obtain the following: which gives (53). The proof is completed.
Proof. A combination of (36a) and (41a) leads to the following: It follows from Theorem 9 that A combination of (4) and (9) gives the following: where and depend only on the compound quadrature (CQ) formula (9), which on substitution into (63) gives the following: Therefore, there exist and such that (61) holds. The proof of Theorem 11 is completed.
Proof. A combination of (36b), (59b), and (53) leads to the following: Using Theorem 9, we have the following: which on substitution into (66) gives the following: where is the minimum characteristic value of . On the other hand, In view of (41a)–(41c) and (59a) and (59b), the application of Theorem 9 leads to the following: which gives where denotes the maximum eigenvalue of the matrix . A combination of (40) and (66)–(71) leads to the following: where