1. Introduction
The integrodifferential equations occur in the mathematical models of many real and scientific phenomena. The exact solutions of equations of this kind are usually difficult to find. Therefore, the various approximate methods are applied for finding their solutions. In this paper, we propose the application of the homotopy analysis method for solving these kind of equations.
The homotopy analysis method was developed in the 1990s [
1,
2]. This method gives the possibility to determine the solution of the wide class of problems described by means of various operator equations. In particular, it can be the ordinary and partial differential equations [
3,
4], the fractional differential or integrodifferential equations [
5,
6], and also the integral equations [
7,
8]. Theoretical results concerning the convergence of the homotopy analysis method in the case of differential equations can be found, among others, in References [
2,
9,
10], whereas, in the case of integral equations, can be found in References [
8,
11].
The homotopy analysis method was already applied for solving the integrodifferential equations in some papers (see, for example, References [
12,
13]). In most of these papers, presented are only the examples of applications without any theoretical results concerning the convergence of the method. The theoretical result, the most frequently cited in the other papers, is the classical theorem stating that, if the series created in the homotopy analysis method is convergent, then its sum is the solution of the input equation. There are only two papers [
14,
15] that contain also some other theoretical results. In Reference [
14], the theorem defining the sufficient conditions for convergence and divergence of the created series is formulated. Next, in Reference [
15], the authors used the homotopy analysis method combined with the parametrization methods to solve the optimal control problems governed by the integrodifferential equations. This paper includes also the sufficient condition for convergence of the obtained series (in addition to the classical theorem mentioned before).
In the current paper, except the classical theorem, we also formulate and prove the sufficient condition for convergence of the created series. Moreover, we give the estimation of error of the approximate solution obtained by reducing the sum of the series to the partial sum. Additionally, we present the example of using the examined method for determining the vibrations of the freely supported reinforced concrete beam as well as for solving the equation of movement of the electromagnet jumper mechanical system.
2. Homotopy Analysis Method
The homotopy analysis method serves to solve the following operator equations:
where
N denotes the operator (in particular, it may be the nonlinear operator), while
u is the unknown function. In the first step, we specify the homotopy operator
as given below:
where
is the embedding parameter,
is the convergence control parameter [
2,
9,
16,
17],
denotes an initial approximation of the solution of Equation (
1), and
L is the auxiliary linear operator with property
.
Considering the equation
, we obtain the so-called zero-order deformation equation:
Substituting
, we get
, so we have
. However, if we assume
, we obtain
. Therefore,
, where
u is the desired solution of Equation (
1). Thus, the change of parameter
p from zero to one involves a change from the trivial problem to the original problem (and, thus, the solution from
to
u). If the operator
N has a unique kernel, then Equation (
1) has a single solution, which could be determined by using the described method. If the kernel is not unique, the method allows to determine one of the existing solutions of Equation (
1).
Taking the Maclaurin series of function
(with respect to parameter
p), we get
where
and
If the above series is convergent for
, then we obtain the required solution:
In order to determine the function
, we differentiate
m times with respect to parameter
p the left- and the right-hand sides of dependence of Equation (
3). Then, we divide the result by
and we substitute
, obtaining the so-called
mth-order deformation equation (
):
where
, and
and
If we are not able to determine the sum of series in Equation (
6), then we can accept the partial sum of this series as the approximate solution of considered equation:
The proper selection of the convergence control parameter
h affects the area of convergence of the created series as well as the rate of this convergence [
2,
10,
18]. In order to determine the value of this parameter, we can use the so-called “optimization method” [
2], in which we define the following squared residual of governing equation:
Next, the optimal value of the convergence control parameter is obtained by minimizing the above squared residual.
3. Integrodifferential Equations
We consider the integrodifferential equations of the form
with condition
, where
,
, and where, for
,
is the bounded linear operators
,
,
and
, whereas the function
u is sought. In this elaboration, we will use the norm in space
, where
is the compact set, defined as follows:
and the norm in
is defined as follows:
We assume
or
. The norm in space
can be also defined as given below:
where
w is a fixed element of set
.
Integrating both sides of Equation (
12), we get the equivalent form of this equation:
where we used the condition
and we assumed that
The above equation gives the ground for defining operators
L and
N:
By setting
and by applying the homotopy analysis method, we get the following formula for functions
:
where
and operator
are determined by Equations (
8) and (
9).
In the considered case, the operator
takes the following form (under assumption that the series is convergent, which will be discussed later):
Hence, we obtain the following formulae for function
:
and for
:
Now, we intend to present and to prove the following three theorems describing the application of the homotopy analysis method for solving linear integrodifferential equations.
Theorem 1. Let the functions and , be defined by Equations (17) and (18). If the series in Equation (6) is uniformly convergent, then the sum of this series is the solution of Equation (12). Proof. Let us assume that the series in Equation (
6) is convergent. For any
, according to the necessary condition for convergence of the series, we have
Using the definition of operator
L, we can write
From Equation (
7), we get
And because
, we obtain
By using the convergence of series in Equation (
6) and the continuity of respective operators, we get
This means that the function
satisfies Equation (
13), and therefore, it satisfies Equation (
12). □
In the next theorem, we prove the sufficient condition for convergence of the series created in the discussed method.
Theorem 2. If the inequalityis satisfied, then with a suitable choice of the control parameter h, the series in Equation (6) is uniformly convergent in the interval . Proof. Let
be the function of class
. Let us find the boundaries of functions
in interval
:
where
It can be shown by induction that, for
and
, we have
Thus, for the considered series in Equation (
6), we get
In the above estimation, the last series is the geometric series with common ratio . Thus, if only (of course we have ), then based on the comparison test, we state that the investigated series is uniformly convergent in the interval .
Let us investigate whether one can choose such value of parameter
h, that
, that is
We consider the above inequality in three cases: for
, for
and for
. In the first case, we get the conditions
In the second case, we have
, whereas the third case leads to contradiction. It means that, if the inequality of Equation (
19) is fulfilled, then we can select such a value of the parameter controlling convergence (it would be enough to take any
) that
. □
The estimations introduced in the previous proof will be also used in the proof of the following theorem.
Theorem 3. If the inequality (19) is fulfilled and , then for the approximate solution , we have the following error estimation:for all . Thus, we havewhere and . Proof. For all
, we have
Now, we proceed to the proof of the second inequality. For this purpose, we apply the homotopy analysis method for Equation (
12) by defining the operators
L and
N in the following way:
Then, we get the following relations for the derivatives of functions
:
and for
:
Now, we determine the boundaries of functions
:
Hence, we get in general for
By using the above estimations, we have for any
which proves the inequality of Equation (21).
The last inequality follows from the definition of the taken norm:
and from the inequalities already proven. □
4. Examples
Example 1. We start with the solution of the following Volterra–Fredholm equation:with condition . In this equation, we have , , , , and , , where , and . The solution of Equation (14) is then the function . Let us mention that all the calculations presented here were carried out with the aid of Mathematica software. In the considered equation, we have Hence, we getwhich means that the sum of the created series is the solution of the discussed equation. Assuming the initial approximation , in the first step, we get The optimal value of the convergence control parameter equals . Figure 1 displays the graph of logarithm of the squared residual for . Table 1 compiles the errors of reconstruction of the exact solution for the successive approximate solutions , . As revealed by the above results, together with an increase of the number of components in Equation (10), the errors quickly decrease. The errors decrease most quickly for the optimal value of parameter h. For this value, the approximate solution provides the approximation of the sought function with the error not higher than , whereas the approximate solution gives the error not higher than . Also, for example, in the case of , we obtain, respectively, and . Table 1 contains also the estimation of error resulting from the inequalityof Equation (22). In the considered example, we have , for the optimal value of convergence control parameter . Equation (22) gives the estimation of error of the approximate solution (the worst possible case). In fact, the errors of approximate solutions are, in general, much smaller than the value determined in the right-hand side of this inequality. In Figure 2, the absolute error for and is plotted. The obtained results show that the investigated method is very quickly convergent and that the calculation of only a few first terms of the series ensures very good approximation of the exact solution. Example 2. Vibrations of the freely supported reinforced concrete beam of span m and cross section () m/m can be described by means of the linear differential equation. Evenly distributed mass of the beam per length unit, by assuming that the density of beam material (reinforced concrete) is of value kg/m, equals The discussed beam is considered as the set of segments of length , and the kinetic energy of such an elementary segment, that is, the element of mass (see Figure 3a) [19,20], is expressed as follows:where denotes the velocity of point of coordinate x (in the beam). In case of the freely supported beam wrapped by the mass , the deflection curve is symmetrical with respect to the centre of gravity and, for , is expressed by the following relation:where denotes the largest deflection of the beam (see Figure 3a). The velocity of each point of coordinate x is described by the following formula:where and the kinetic energy of the whole beam is equal to By substituting Equation (26) into Equation (27) and by integrating, we get The kinetic energy described by Equation (28) corresponds to the kinetic energy of the system with one degree of freedom, and the substitute mass equalsattached in the vibration arrow, that is, at the point where the displacement of beam is the highest (see Figure 3b). Stiffness of the substitute system is expressed by the following formula:where GPa is the modulus of elasticity of the beam material and I denotes the moment of inertia of the cross section . Thus, we have MN/m. The equation describing the first free vibration of the reinforced concrete beam of span m and cross section () m/m has then the following form:where kg is the beam substitute mass and MN/m describes the beam stiffness in the center of span (see Figure 4). Equation (31), after dividing by m, takes the formwhere denotes the circular frequency of undamped vibrations. The above equation must be completed by the initial conditions Now, we will demonstrate how to apply the homotopy analysis method to solve Equation (32) with the conditions of Equation (33). For this purpose, let us integrate Equation (32) in limits from 0 to t, obtainingwith condition . For the initial conditions and , the exact solution takes the form Equation (34) can be solved by using the procedure described in this paper. If we take the function as the initial approximation, then in the first steps of the method, we get Figure 5 displays the logarithm of squared residual for . The optimal value of the convergence control parameter, determined numerically, is equal to . In Table 2, there are collected the errors (in norm ) of reconstructing the exact solution by the successive approximate solutions , . At the beginning, that is, until , the errors increase, but after that, the errors start quickly decrease. So, for , the errors are at the level of ; for , they are at the level of ; and for , the errors are at the level of . The main part of these errors is generated by the approximation of derivative. For , the error in norm is equal to , whereas in norm , it is only equal to . A similar situation can be observed in the other cases. As revealed by the above results, starting from some moment, we can observe that an increase of the number of components in Equation (10) causes the errors to quickly decrease (see also Figure 6). The investigated equation shows that the homotopy analysis method can be successfully applied also for solving the equations not satisfying the condition in Equation (19). In the above considered case, we have because the condition of Equation (19) is only the sufficient condition and it is not the necessary condition for convergence of the method. Example 3. The equation of movement of the electromagnet jumper mechanical system can be presented in the following way:where m is the mass, denotes the damping, y describes the displacement of the electromagnet jumper, and f denotes the force occurring in consequence to the electrical system action. In this example, we take . The third term in the left-hand side of the equation is the result of taking into account in the mechanical system the nonlinear spring, which generates the force component depending on the displacement. The above differential equation is completed by the following initial conditions: Now: we present the application of the homotopy analysis method for solving Equation (32) with the conditions of Equation (33). Equation (35), after integration in limits from 0 to t and some transformations, can be written in the following form:where is the function depending on the initial conditions. In the considered example, we take , , , , , , . In order to test the usage of investigated method, we select the function so that the exact solution of above equation is known ().
By applying the homotopy analysis method for solving Equation (38) and by taking , we obtain in the first step The optimal value of the convergence control parameter, determined numerically, is equal to . In Table 3, there are compiled the errors (in norm ) of reconstructing the exact solution by the successive approximate solutions , . As revealed by the above results, starting from some moment, the errors quickly decrease together with an increase of the number of components in Equation (10). Next, in Figure 7, the absolute error is plotted for and . Again, we may conclude, on the basis of the obtained results, that the method is very quickly convergent to the exact solution. 5. Conclusions
This paper presents the application of the homotopy analysis method for solving linear integrodifferential equations. In this method, we create a series, the successive elements of which are calculated by integrating properly the previous elements. In this paper, it is proven that, if the created series is convergent, then its sum is the solution of the input equation. The sufficient condition of this convergence is also formulated and proven. Moreover, the estimation of error of approximate solution, obtained by reducing the considered series to its partial sum, is also given.
The series created with the discussed method is, in general, very quickly convergent. After selecting the optimal value of the convergence control parameter, calculation of only a few or several first terms of the series ensures very good approximation of the sought solution. The executed calculations show that the homotopy analysis method is effective in solving the problems under consideration.
The results obtained in this paper can be easily generalized to the case of equations having the following form:
where all the functions are defined analogously to Equation (
12). In this case, the condition of Equation (
19) should be replaced by the condition in the following form:
The homotopy analysis method can be also applied for solving the equation of the following form:
in a similar way as described above, with the appropriate initial conditions. One has to only integrate
l times the given equation and to repeat the procedure described in this paper.
It is also possible to find another formulation of the homotopy operator [
2], namely
where
H is an auxiliary function. Also in this case, all the results obtained and presented in this work are true. The only change concerns the formula for constant
.
Author Contributions
Conceptualization, D.S. and R.W.; formal analysis, D.S.; investigation, D.S., E.H., R.W., K.G. and T.T.; methodology, D.S., E.H. and R.W.; writing—original draft, D.S. and E.H.; writing—review and editing, R.W., K.G. and T.T.
Funding
This research received no external funding.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Liao, S. Homotopy analysis method: A new analytic method for nonlinear problems. Appl. Math. Mech. 1998, 19, 957–962. [Google Scholar]
- Liao, S. Homotopy Analysis Method in Nonlinear Differential Equations; Springer: Berlin, Germany; Higher Education Press: Beijing, China, 2012. [Google Scholar]
- Van Gorder, R. Control of error in the homotopy analysis of semi-linear elliptic boundary value problems. Numer. Algor. 2012, 61, 613–629. [Google Scholar] [CrossRef]
- Hetmaniok, E.; Słota, D.; Wituła, R.; Zielonka, A. An analytical method for solving the two-phase inverse Stefan problem. Bull. Pol. Acad. Sci. Tech. Sci. 2015, 63, 583–590. [Google Scholar] [CrossRef] [Green Version]
- Zurigat, M.; Momani, S.; Odibat, Z.; Alawneh, A. The homotopy analysis method for handling systems of fractional differential equations. Appl. Math. Model. 2010, 34, 24–35. [Google Scholar] [CrossRef]
- Zhang, X.; Tang, B.; He, Y. Homotopy analysis method for higher-order fractional integro-differential equations. Comput. Math. Appl. 2011, 62, 3194–3203. [Google Scholar] [CrossRef] [Green Version]
- Abbasbandy, S.; Shivanian, E. A new analytical technique to solve Fredholm’s integral equations. Numer. Algor. 2011, 56, 27–43. [Google Scholar] [CrossRef]
- Hetmaniok, E.; Nowak, I.; Słota, D.; Wituła, R. Convergence and error estimation of homotopy analysis method for some type of nonlinear and linear integral equations. J. Numer. Math. 2015, 23, 331–344. [Google Scholar] [CrossRef]
- Van Gorder, R.; Vajravelu, K. On the selection of auxiliary functions, operators, and convergence control parameters in the application of the homotopy analysis method to nonlinear differential equations: A general approach. Commun. Nonlinear Sci. Numer. Simulat. 2009, 14, 4078–4089. [Google Scholar] [CrossRef]
- Odibat, Z. A study on the convergence of homotopy analysis method. Appl. Math. Comput. 2010, 217, 782–789. [Google Scholar] [CrossRef]
- Hetmaniok, E.; Słota, D.; Trawiński, T.; Wituła, R. Usage of the homotopy analysis method for solving the nonlinear and linear integral equations of the second kind. Numer. Algor. 2014, 67, 163–185. [Google Scholar] [CrossRef]
- Behzadi, S.; Abbasbandy, S.; Allahviranloo, T.; Yildirim, A. Application of homotopy analysis method for solving a class of nonlinear Volterra-Fredholm integro-differential equations. J. Appl. Anal. Comput. 2012, 2, 127–136. [Google Scholar]
- Mohamed, M.; Gepreel, K.; Alharthi, M.; Alotabi, R. Homotopy analysis transform method for integro-differential equations. Gen. Math. Notes 2016, 32, 32–48. [Google Scholar]
- Ghanbari, B. The convergence study of the homotopy analysis method for solving nonlinear Volterra-Fredholm integrodifferential equations. Sci. World J. 2014, 2014. [Google Scholar] [CrossRef] [PubMed]
- Alipouru, M.; Vali, M.; Borzabadi, A. A direct approach for approximate optimal control of integro-differential equations based on homotopy analysis and parametrization method. IMA J. Math. Control Inf. 2017, 34, 611–630. [Google Scholar] [CrossRef]
- Liao, S. Notes on the homotopy analysis method: Some definitions and theorems. Commun. Nonlinear Sci. Numer. Simulat. 2009, 14, 983–997. [Google Scholar] [CrossRef]
- Russo, M.; Van Gorder, R. Control of error in the homotopy analysis of nonlinear Klein-Gordon initial value problems. Appl. Math. Comput. 2013, 219, 6494–6509. [Google Scholar] [CrossRef]
- Turkyilmazoglu, M. Convergence of the homotopy analysis method. arXiv 2010, arXiv:1006.4460v1. [Google Scholar]
- Gromysz, K. Investigations of Reinforced-Concrete Composite Slabs Loaded Temporarily, Cyclically And Kinematically; Publishing House of the Silesian University of Technology: Gliwice, Poland, 2013. [Google Scholar]
- Gromysz, K. Nonlinear analytical model of composite concrete slab free and forced vibrations. Procedia Eng. 2017, 193, 281–288. [Google Scholar] [CrossRef]
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).