Numerical Algorithms for Solving One Type of Singular Integro-Differential Equation Containing Derivatives of the Time Delay States

Shihchung Chiang^1*, Terry L. Herdman²

¹Department of Finance, Chung Hua University, Hsinchu.
²Department of Mathematics, Virginia Tech, Blacksburg, VA, USA.

DOI: 10.4236/am.2015.68123                

Abstract

This study presents numerical algorithms for solving a class of equations that partly consists of derivatives of the unknown state at previous certain times, as well as an integro-differential term containing a weakly singular kernel. These equations are types of integro-differential equation of the second kind and were originally obtained from an aeroelasticity problem. One of the main contributions of this study is to propose numerical algorithms that do not involve transforming the original equation into the corresponding Volterra equation, but still enable the numerical solution of the original equation to be determined. The feasibility of the proposed numerical algorithm is demonstrated by applying examples in measuring the maximum errors with exact solutions at every computed nodes and calculating the corresponding numerical rates of convergence thereafter.

Keywords

Integro-Differential Equation of the Second Kind, Weakly Singular Kernel, Numerical Algorithms, Rates of Convergence

1. Introduction

A dynamical system describing a two-dimensional physical thin airfoil moving inside an incompressible flow was introduced by Burns, Cliff, and Herdman [1] in 1983. The system contains a form of linear singular integro-differential equations with integration over a deterministic interval (i.e., equations not of the Volterra types). Other studies [2] [3] have presented the well-posedness of the problem regarding specific product spaces and the exact solutions of the original class of integro-differential equations of the first kind, and have reported numerical methods and corresponding numerical results [4] [5] . Associated optimal control problems are topics discussed in [6] . Another study [7] applied semigroup theory to this particular type of equation and constructed an associated abstract Cauchy problem. The current study presents a numerical algorithm for solving the type of equations containing not only the original aeroelastic integro-differential term as a part of the equation but also time-derivative states evaluated at different previous times. This new linear equation is in the category of “integro-differential equations of the second kind”. The main purpose of this study is to develop feasible numerical algorithms for solving this type of integro-differential equation. According to previous studies (for example, [8] ), all existing numerical methods can be used for solving only integro-differential equations of the second kind that can be transformed into Volterra integral equations of the second kind that linearly containing the state, and no numerical method (except the papers by current authors) has been proposed for solving the integro-differential equations of the second kind directly and the integro-differential equations of the second kind containing time delay states. The remainder of this paper is organized as follows: Section 2 presents the derivation of the associated Volterra integral equations of the second kind. Section 3 presents numerical algorithms used for directly solving singular integro-differential equations of the second kind. Section 4 presents the numerical results of test examples obtained by applying the numerical method described in Section 3. Finally, Section 5 presents a summary of this study.

2. Problem Description

Consider the class of an integro-differential equation of the second kind expressed as follows:

(1)

and the initial condition

(2)

where are constants and, are nonnegative constants. The term is the derivative of the delay state with respect to t, and the difference operator D is defined as

.

The second part of the integrand represents

,

and the first part is a weakly singular function

,

that is integrable, positive, nondecreasing, and weakly singular at. Assume the forcing term is locally integrable for Although a more general kernel is also suitable, this study focuses on the Abel- type kernel and considers and for A specific value of corresponds to the original aeroelastic problem. Assume that the initial condition, for is in space, a weighted space with weight.

If the differential part of the integro-differential term can be removed, that is, the term exists, then applying the integration to Equation (1) forms a new equation of the following form:

This equation can be developed into a Volterra integral equation of the second kind, provided that the function

is absolutely continuous with respect to, and the product of the kernel and initial functions, , belongs to. Therefore, the corresponding weakly singular Volterra integral equation of the second kind is

3. Numerical Algorithms

The proposed algorithms involve using the separating variables method to directly solve the numerical solution of Equations (1) and (2). Without loss of generality, assume that and, the equation is expressed as

(3)

with initial data

(4)

where, is a locally integrable function.

Let, then Equation (3) can be divided into two categories:, and.

3.1.

For this category, following study [6] , define a new functional such that

(5)

Reformulate Equation (3) as a first-order hyperbolic partial differential equation

(6)

with the condition

(7)

Next, assume that the solutions to Equations (6) and (7) have the form

, (8)

where the bases, , are

Specifically, and are piecewise linear functions. The mesh points, are defined by and, for One restriction for the

mesh points is, namely, the time lag terms coincide with some of the absolute

values of mesh points.

After substituting the form of previously defined in Equation (8) into Equations (6) and (7), the governing equations for and become the following equations:

(9)

and

(10)

By the property of the bases, rewrite Equation (10) as

(11)

where are the corresponding terms of with respect to,

Define

and Equations (9) and (11) thus become

(12)

and

(13)

This produces the following linear system of first-order ordinary differential equations:

, (14)

where

s represent certain values depending on the typical equation, and, in which T

is the transpose of the corresponding vector.

The procedure for obtaining the initial condition for the first-order ordinary differential system (14) is described as follows: For the initial condition, combine Equations (4), (5), and (8), and fix; the state thus becomes

The structure of and indicates that is equal to for

Next, to determine apply an ordinary differential equation solver (Matlab “ode45”)

to the system (14). Two methods can be used to solve, , depending on the setting of variables: fix or in Equation (8). According to the property at the two choices become two cases for the solution,:

Case 1:

and Case 2:

In Case 1, solve for based on Equation (14) and set Thus, yields the

corresponding solutions In Case 2, solve for by using Equation (14).

Subsequently, set to obtain for Therefore, is the solution for

A similar procedure can be extended to solve, for.

3.2.

For this category, Equation (3) can be rewritten as

then it becomes

a similar form of Equation (3) except for the integral interval of the second term on the left hand side, but this new equation can be treated by reconsidering the discretization interval to be; namely, by resetting the mesh points as, and then follow the procedures introduced in Section 3.1.

4. Numerical Examples

Consider examples involving, initial conditions, and forcing terms

, for

Example 1:;, , , , ,;

, or 1000;

Exact solution:

Example 2:;, , , ,;

, or 1000;

,

Exact solution:

Example 3:;, , , ,;

, or 1000;

Exact solution:

Example 4:;, , , ,;

, or 1000;

where is a gamma function.

Exact solution:

Example 5:;, , , ,;

, or 1000;

Exact solution:

Example 6:;, , , ,;

, or 1000;

Exact solution:

The feasibility of the proposed methods are determined by the maximum errors at every computed nodes after applying different number of mesh points, the formula is

for is the number of mesh points. is the computed solution and is the exact solution. The rate of convergence is defined as

,

for is a positive integer.

Table 1 and Table 2 contain the maximum errors at every computed nodes and mean rates of convergence evaluated at for the examples.

Although the mean rates of convergence for the linear cases (solutions are linear: and initial

Table 1. The maximum errors at every computed nodes for the examples.

Table 2. Mean rates of convergence evaluated at for the examples.

conditions:) such as Example 2 and Example 3 in Table 2 have some vibration phenomena, the maximum errors in Table 1 provide sufficient evidence for the correctness of the numerical solutions.

Remark

This study presents a numerical method for directly solving the integro-differential equations of the second kind. The method involves discretizing the space s, and retains the variable t. The unknown states are repre- sented by To solve system (14), which is a semi-discretized scheme, the authors suggest using an ordinary differential equation solver. The (mean) rates of convergence can be determined, although it depends on the separating variable form of the state as well as on the accuracy of the ordinary differential equation solver applied (shown in the coming papers). Another approach to determining the rate of convergence in this observed study is to discretize both variables s and t, and this process results in a full-discretized scheme, as described in [3] .

5. Summary

This study presents a numerical method for solving a class of singular integro-differential equations of the second kind that contain derivatives of the states at previous certain times of the finite history interval, as well as an integro-differential term containing a weakly singular kernel. The proposed equations can be transformed into Volterra integral equations of the second kind if the integro-differential term is integrable. This study presents direct numerical methods to the proposed equation. The tables of corresponding maximum errors and the mean rates of convergence show the feasibility of using the proposed numerical method for the equations.

NOTES

^*Corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Burns, J.A., Cliff, E.M. and Herdman, T.L. (1983) A State-Space Model for an Aeroelastic System. Proceedings: 22nd IEEE Conference on Decision and Control, San Antonio, December 1983, 1074-1077. http://dx.doi.org/10.1109/CDC.1983.269685
[2]	Burns, J.A. and Ito, K. (1995) On Well-Posedness of Solutions to Integro-Differential Equations of Neutral-Type in a Weighted-L2 Spaces. Differential and Integral Equations, 8, 627-646.
[3]	Herdman, T.L. and Turi, J. (1991) An Application of Finite Hilbert Transforms in the Derivation of a State Space Model for an Aeroelastic System. Journal of Integral Equations and Applications, 3, 271-287. http://dx.doi.org/10.1216/jiea/1181075618
[4]	Herdman, T.L. and Turi, J. (1991) On the Solutions of a Class of Integral Equations Arising in Unsteady Aerodynamics. In: Elaydi, S., Ed., Differential Equations: Stability and Control, Department of Mathematics, Trinity University, San Antonio, Texas, 241-248.
[5]	Kappel, F. and Zhang, K.P. (1986) Equivalence of Functional Equations of Neutral Type and Abstract Cauchy Problems. Monatshefte für Mathematik, 101, 115-133. http://dx.doi.org/10.1007/BF01298925
[6]	Chiang, S. and Herdman, T.L. (2013) Revised Numerical Methods on the Optimal Control Problem for a Class of Singular Integro-Differential Equations. Mathematics in Engineering, Science and Aerospace, 4, 176-189.
[7]	Ito, K. and Turi, J. (1991) Numerical Methods for a Class of Singular Integro-Differential Equations Based on Semigroup Approximation. SIAM Journal on Numerical Analysis, 28, 1698-1722. http://dx.doi.org/10.1137/0728085
[8]	Lubich, Ch. (1985) Fractional Linear Multistep Methods for Abel-Volterra Integral Equations of the Second Kind. Mathematics of Computation, 45, 463-469. http://dx.doi.org/10.1090/S0025-5718-1985-0804935-7