# An Analytical and Approximate Solution for Nonlinear Volterra Partial Integro-Differential Equations with a Weakly Singular Kernel Using the Fractional Differential Transform Method.

1. IntroductionMany engineering and physical problems result in the analysis of the nonlinear weakly singular Volterra integral equations (WSVIEs). These equations are applied in many areas [1] such as reaction-diffusion problems in small cells [2], theory of elasticity, heat conductions, hydrodynamics, stereology [3], the radiation of heat from semi-infinite solids [4], and other applications. Such equations have been studied by several authors [5-14].

The aim of this paper is applying the fractional differential transform method (FDTM) for solving WSVIE. The fractional differential transform method has recently been developed for solving the differential and integral equations. For example, in [15], FDTM is applied for fractional differential equations and in [16] it is used for fractional integro-differential equations. This method is applied to nonlinear fractional partial differential equations in [17]. The use of the differential transform method (DTM) in electric circuit analysis was first proposed by Zhou [18].

The main challenge of partial integro-differential equations (PIDEs) with a weakly singular kernel is faced when we are looking for an analytical solution. By applying the differential transform method, the result mostly obtained is an analytical solution in the form of a polynomial. The differential transform method is different from the traditional high order Taylor series method, which requires symbolic competition of the necessary derivatives of the data functions. Making use of this method enables us to obtain highly accurate results or exact solutions for a partial integro-differential equation. The use of application of DTM and FDTM does not require linearization, discretization, or perturbation in contrast to the methods discussed in the literature [8, 9, 19].

The form of WSVIE that we will consider in this paper with FDTM is

[mathematical expression not reproducible], (1)

where 0 < p < 1, 0 [less than or equal to] [xi] [less than or equal to] x, 0 [less than or equal to] [eta] [less than or equal to] t, and (x, t) [member of] [0,1] x [0,1] with the initial condition

[phi](x, 0) = [[phi].sub.0](x), (2)

where [phi] is an unknown function in [LAMBDA](= [0, 1] x [0,1]) which should be determined and f(x, t), [[phi].sub.0](x) are known functions and N is a nonlinear operator. The given functions f(x, t) and [phi](x, t) are assumed to be sufficiently smooth in order to guarantee the existence and uniqueness of a solution [phi] [member of] C(A). It is assumed that the nonlinear term N[[phi]] satisfies the Lipschitz condition in [L.sup.2]([LAMBDA]).

The numerical treatment of (1) is not simple because the solutions of WSVIEs usually have a weak singularity at x = 0. Different numerical techniques have been developed for the solution of PIDEs [8, 10, 20-25]. In this article, FDTM is applied to solve (1) and the main theorem is proved on the two-dimensional FDTM, while the one-dimensional FDTM has been applied in [26].

The paper is organized as follows: In Section 2, Caputo and Riemann-Liouville fractional derivatives are introduced. In Section 3, the theorems of the fractional differential transform method, preliminaries, and notations are explained. In Section 4, we have proposed the main theorem, for which a WSVIE can be considered as a series of FDT. Further, some examples of the application of FDTM are demonstrated, which show the accuracy of the method, in Section 5. We conclude our discussion in Section 6.

2. Riemann-Liouville and Caputo Fractional Derivatives

There are different kinds of definitions for the fractional derivative of order q > 0; among various definitions of fractional derivatives of order q > 0, the Riemann-Liouville and Caputo formulas are the most common [27]. The Riemann-Liouville fractional integration of order q is defined as

[mathematical expression not reproducible]. (3)

The following equations define Riemann-Liouville and Caputo fractional derivatives of order q, respectively:

[mathematical expression not reproducible], (4)

[mathematical expression not reproducible]. (5)

where m - 1 [less than or equal to] q < m and m [member of] N. From (4) and (5), we have

[mathematical expression not reproducible]. (6)

3. The Fractional Differential Transform Method (FDTM)

There are some approaches to the generalization of the notion of differentiation to fractional orders. According to the Riemann-Liouville formula, the fractional differentiation is defined by (6). The analytical and continuous function f(x) is expended in terms of a fractional power series as follows:

f(x) = [[infinity].summation over (k=0)] F(k)[(x - [x.sub.0]).sup.k/[alpha]], (7)

where [alpha] is the order of fraction and F(k) is the fractional differential transform of f(x) [15, 28-30]. In Caputo sense [31], (6) is modified to handle integer-order initial conditions as follows:

[mathematical expression not reproducible]. (8)

Since the initial conditions are implemented to the integer-order derivatives, the transformation of the initial conditions also can be represented as follows:

[mathematical expression not reproducible], (9)

for k = 0, 1, 2, ..., (n[alpha] - 1), where n is the order of FDE that is considered.

Consider a function u(x, t) of two variables, and assume that it can be expressed as a product of two single-variable functions as u(x, t) = f(x)g(t). The expansion of the function u(x, t) in a Taylor series around a point ([x.sub.0], [t.sub.0]) is as follows:

[mathematical expression not reproducible]. (10)

If we take ([x.sub.0], [t.sub.0]) as (0, 0), then (10) can be illustrated as

[mathematical expression not reproducible], (11)

where [alpha], [beta] are the order of the fractions, [alpha], [beta] [member of] N, and U(k, h) = F(k)G(h) is called the spectrum of u(x, t) and defined by

[mathematical expression not reproducible]. (12)

If we choose [alpha] = 1 and [beta] = 1, the fractional two-dimensional differential transform reduces to the classical two-dimensional differential transform. Using (11) and (12), the theorems of FDTM are introduced as follows. The proofs of these theorems can be found in [15,17].

Theorem 1. Suppose that W(k, h), U(k, h), and V(k, h) are the differential transformations of the functions w(x, t), u(x, t), and v(x, t), respectively, with order of fraction a and [beta]; then

(1) if w(x, t) = u(x, t) [+ or -] v(x, t), then W(k, h) = U(k, h) [+ or -] V(k, h);

(2) if w(x, t) = [lambda]u(x, t), then W(k, h) = [lambda]U(k, h);

(3) if w(x, t) = [(x - [x.sub.0]).sup.p] [(t - [t.sub.0]).sup.q], then

[mathematical expression not reproducible]; (13)

(4) if w(x, t) = u(x, t)v(x, t), then

W(k, h) = U(k, h)V(k, h)

= [k.summation over (r=0)] [h.summation over (s=0)] (r, h - s)v(k - r, s); (14)

(5) if w(x, t) = [x.sup.p] sin(at + b), then

W(k, h) = [a.sup.h]/h! [delta](k - p) sin (h[pi]/2 + b); (15)

(6) if w(x, t) = [x.sup.p] cos(at + b), then

W(k, h) = [a.sup.p] [delta](k - p) cos (h[pi]/2 + b); (16)

(7) if w(x, t) = [x.su.p] exp([lambda]t), then

W(k, h) = [delta](k - p)([[lambda].sup.h])/h!, (17)

where p and q are positive and [lambda], a, b are scalars.

Theorem 2. If u(x, t) = [D.sup.q.sub.x] v(x, t) and [beta] [member of] N is order of fractional, then

U(k, h) = [GAMMA](q + 1 + k/[beta])/[GAMMA](1 + k/[beta]) V(k + [beta]q, h). (18)

Definition 3. The Beta function B(a,b) of two variables is defined by

B(a, b) = [[integral].sup.1.sub.0] [(1 - x).sup.a-1] [x.sup.(b-1)]dx. (19)

The following is proved easily:

B(a, b) = [GAMMA] (a) [GAMMA](b)/[GAMMA] (a + b). (20)

Definition 4. The Kronecker delta function is given by

[mathematical expression not reproducible]. (21)

4. Main Theorem

Now, we represent the main theorem of this study, through which a weakly singular Volterra integral equation can be expressed as a series of fractional differential transform for

[mathematical expression not reproducible].

Theorem 5. Suppose that [PHI](k, h) and W(k, h) are the fractional differential transforms of the functions [phi](x, t) and w(x, t), respectively, such that

[mathematical expression not reproducible]. (22)

Then, by choosing a suitable [beta] [member of] [z.sup.+] such that [beta]p [member of] [z.sup.+], we have

[mathematical expression not reproducible], (23)

where B(x, x) is the Beta function and [delta] is the Kronecker delta function.

Proof. By putting [mathematical expression not reproducible] in (22), it will change into

[mathematical expression not reproducible]. (24)

To calculate [[integral].sup.x.sub.0] [(x - [xi]).sup.p-1] [[xi].sup.k]/[beta]] d[xi], we change variable x/[xi], = v and according to Definition 3 we get

[mathematical expression not reproducible]. (25)

The following equation is obtained by using Theorem 1 and replacing (25) into (24):

[mathematical expression not reproducible]. (26)

Hence

[mathematical expression not reproducible]. (27)

According to w(x, t) = [[summation].sup.[infinity].sub.k=0] [[summation].sup.[infinity].sub.h=0] W(k, h) [x.sup.k/[beta]] [t.sup.h/[alpha]], one can conclude that

[mathematical expression not reproducible]. (28)

The proof is completed.

5. Description of Method

In this section, we try to describe the FDTM for (1) and initial condition (2). Based on Theorems 1, 2, and 5, FDTM for (1) and (2) would result as follows:

[mathematical expression not reproducible]. (29)

in which [PHI](k, h) and F(k, h) are the FDTM of [phi](x, t) and f(x, t), respectively. Then according to the recurrence relation (11) the unknown function would result.

6. Applications

In this section, we take some examples to clarify the advantages and the accuracy of the fractional differential transform method (FDTM) for solving a kind of nonlinear partial integro-differential equation with a weakly singular kernel. For each of these examples, we obtain a recurrence relation. In all of the examples, we choose [alpha] = 1 and [beta] is chosen in a way where [beta]p [member of] [z.sup.+].

Example 1. Consider the following nonlinear partial integro-differential equation with a weakly singular kernel with p = 1/2 [23]:

[mathematical expression not reproducible], (30)

with the initial condition [phi](x, 0) = [x.sup.2] and f(x, t) = 2t - [x.sup.2] - t2 - (2/315)[x.sup.1/2] t(128[x.sup.4] + 112[x.sup.2][t.sup.2] + 63[t.sup.4]).

For solving (30), we employ the FDTM, to get

[mathematical expression not reproducible], (31)

where h, k [greater than or equal to] 1 in the upper bound of the sigmas and the differential transform of initial condition is as follows:

[mathematical expression not reproducible]. (32)

Also we have

[mathematical expression not reproducible]. (33)

By using the recurrence relation (31) and the transform initial condition (32), we get the following:

[mathematical expression not reproducible]. (34)

and by applying the same calculations, the following can be concluded:

[PHI](k, 1) = 0 k [greater than or equal to] 5. (35)

Also we put

[mathematical expression not reproducible]. (36)

By Definition 3, we have

B(5, 1/2) = 2 x 128/315. (37)

So [PHI](9, 2) = 0.

And by applying the same calculations, we can conclude that

[PHI](k, 2) = 0 k > 9. (38)

By continuing this process, we can also conclude the following:

[PHI](k, h) = 0 k [greater than or equal to] 0, h [greater than or equal to] 3. (39)

Therefore, by substituting the above values into (33), the exact solution is obtained in the following form:

[PHI](x, t) = [x.sup.2] + [t.sup.2] (40)

which is the particular solution obtained in [23].

Example 2. Consider the following nonlinear partial integro-differential equation with a weakly singular kernel:

[mathematical expression not reproducible] (41)

with the initial condition [phi](x, 0) = x and f(x, t) = x sin t - x cos t - (64/231)[x.sup.11/4]t - (32/231)[x.sup.11/4] sin 2t.

Taking into consideration the two-dimensional transform for (41) and the related theorems, we have

[mathematical expression not reproducible]. (42)

The differential transform of the initial condition is as follows:

[mathematical expression not reproducible]. (43)

Also we have

[mathematical expression not reproducible]. (44)

By using the recurrence relation (42) and the differential transform of initial condition (43), we get

[mathematical expression not reproducible] (45)

and by applying the same calculations, it can be concluded that

[mathematical expression not reproducible], (46)

and by applying the same calculations, we conclude the following:

[mathematical expression not reproducible]. (47)

Therefore, by substituting the above values into (44), the exact solution is obtained in the following form:

[mathematical expression not reproducible]. (48)

Example 3. Consider the following nonlinear partial integro-differential equation with a weakly singular kernel:

[mathematical expression not reproducible] (49)

with the initial condition [phi](x, 0) = x and f(x, t) = (125/168)[x.sup.12/5][e.sup.2t].

To solve (49), by applying FDTM, we have

[mathematical expression not reproducible], (50)

where k [greater than or equal to] 2, h [greater than or equal to] 1 in the upper bound of the sigmas and differential transform of initial condition is as follows:

[mathematical expression not reproducible]. (51)

Also we have

[mathematical expression not reproducible]. (52)

By using the recurrence relation (50), the differential transform of initial condition (51), and the same calculations of the above-mentioned examples, it is concluded that

[mathematical expression not reproducible]. (53)

Of course this solution is an analytical solution.

7. Conclusion

In this paper, we have described the definition and operation of two-dimensional fractional differential transform; fractional derivatives have been considered in the Caputo and Riemann-Liouville sense and the main theorem on fractional differential transform method. Using the fractional differential transform method, a kind of nonlinear partial integro-differential equation with a singular kernel was solved approximately and analytically. We have used FDTM in this paper to solve (30) which was solved by operational matrices in [23]. The advantages of this method are that one obtains satisfactory results in less time, there is no need to calculate anyrepeated integral, and there is no discretization.

https://doi.org/10.1155/2018/7237680

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1] G. Capobianco, D. Conte, and I. Del Prete, "High performance parallel numerical methods for Volterra equations with weakly singular kernels," Journal of Computational and Applied Mathematics, vol. 228, no. 2, pp. 571-579, 2009.

[2] J. A. Dixon, "A nonlinear weakly singular Volterra integro-differential equation arising from a reaction-diffusion study in a small cell," Journal of Computational and Applied Mathematics, vol. 18, no. 3, pp. 289-305, 1987.

[3] P. Linz, Analytical and Numerical Methods for Volterra Equations, vol. 7 of SIAM Studies in Applied Mathematics, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, Pa, USA, 1985.

[4] J. B. Keller and W. E. Olmstead, "Temperature of a nonlinearly radiating semi-infinite solid," Quarterly of Applied Mathematics, vol. 29, pp. 559-566, 1971/72.

[5] A. Palamara Orsi, "Product integration for Volterra integral equations of the second kind with weakly singular kernels," Mathematics of Computation, vol. 65, no. 215, pp. 1201-1212, 1996.

[6] Y. Chen and T. Tang, "Spectral methods for weakly singular Volterra integral equations with smooth solutions," Journal of Computational and Applied Mathematics, vol. 242, pp. 53-69, 2013.

[7] H. Brunner, A. Pedas, and G. Vainikko, "The piecewise polynomial collocation method for nonlinear weakly singular Volterra equations," Mathematics of Computation, vol. 68, no. 227, pp. 1079-1095, 1999.

[8] J. C. Lopez Marcos, "A difference scheme for a nonlinear partial Integrodifferential equation," SIAM Journal on Numerical Analysis, vol. 27, no. 1, pp. 20-31, 1990.

[9] L. Bougoffa, R. C. Rach, and A. Mennouni, "An approximate method for solving a class of weakly-singular Volterra integro-differential equations," Applied Mathematics and Computation, vol. 217, no. 22, pp. 8907-8913, 2011.

[10] M. Dehghan, "Solution of a partial integro-differential equation arising from viscoelasticity," International Journal of Computer Mathematics, vol. 83, no. 1, pp. 123-129, 2006.

[11] P. Baratella and A. P. Orsi, "A new approach to the numerical solution of weakly singular Volterra integral equations," Journal of Computational and Applied Mathematics, vol. 163, no. 2, pp. 401-418, 2004.

[12] T. Diogo and P. Lima, "Superconvergence of collocation methods for a class of weakly singular Volterra integral equations," Journal of Computational and Applied Mathematics, vol. 218, no. 2, pp. 307-316, 2008.

[13] T. Tang, "A finite difference scheme for partial integro-differential equations with a weakly singular kernel," Applied Numerical Mathematics, vol. 11, no. 4, pp. 309-319, 1993.

[14] Y. Chen and T. Tang, "Spectral methods for weakly singular Volterra integral equations with smooth solutions," Journal of Computational and Applied Mathematics, vol. 233, no. 4, pp. 938-950, 2009.

[15] A. Arikoglu and I. Ozkol, "Solution of fractional differential equations by using differential transform method," Chaos, Solitons & Fractals, vol. 34, no. 5, pp. 1473-1481, 2007.

[16] A. Arikoglu and I. Ozkol, "Solution of fractional integro-differential equations by using fractional differential transform method," Chaos, Solitons & Fractals, vol. 40, no. 2, pp. 521-529, 2009.

[17] S. Momani and Z. Odibat, "A novel method for nonlinear fractional partial differential equations: combination of DTM and generalized Taylor's formula," Journal of Computational and Applied Mathematics, vol. 220, no. 1-2, pp. 85-95, 2008.

[18] J. K. Zhou, differential transformation and its applications for electrical circuits, Huazhong University Press, Wuha, China, 1986.

[19] R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific, Singapore, 2000.

[20] E. Tohidi and F. Toutounian, "Convergence analysis of Bernoulli matrix approach for one-dimensional matrix hyperbolic equations of the first order," Computers & Mathematics with Applications, vol. 68, no. 1-2, pp. 1-12, 2014.

[21] F. Toutounian and E. Tohidi, "A new Bernoulli matrix method for solving second order linear partial differential equations with the convergence analysis," Applied Mathematics and Computation, vol. 223, pp. 298-310, 2013.

[22] G. Fairweather, "Spline collocation methods for a class of hyperbolic partial integro-differential equations," SIAM Journal on Numerical Analysis, vol. 31, no. 2, pp. 444-460, 1994.

[23] S. Singh, V K. Patel, V. K. Singh, and E. Tohidi, "Numerical solution of nonlinear weakly singular partial integro-differential equation via operational matrices," Applied Mathematics and Computation, vol. 298, pp. 310-321, 2017.

[24] S. Singh, V. K. Patel, and V. K. Singh, "Operational matrix approach for the solution of partial integro-differential equation," Applied Mathematics and Computation, vol. 283, pp. 195207, 2016.

[25] W. McLean, I. H. Sloan, and V. Thomee, "Time discretization via Laplace transformation of an integro-differential equation of parabolic type," Numerische Mathematik, vol. 102, no. 3, pp. 497-522, 2006.

[26] E. Rahimi, H. Taghvafard, and G. H. Erjaee, "Fractional differential transform method for solving a class of weakly singular Volterra integral equations," Iranian Journal of Science & Technology, vol. 38, no. 1, pp. 69-73, 2014.

[27] D. Nazari and S. Shahmorad, "Application of the fractional differential transform method to fractional-order integro-differential equations with nonlocal boundary conditions," Journal of Computational and Applied Mathematics, vol. 234, no. 3, pp. 883-891, 2010.

[28] A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, New York, NY, USA, Elsevier, 2006.

[29] M. Sen, Introduction to fractional-order operators and their engineering applications, Department of Aerospace and Mechanical Engineering University of Notre Dame, 2014.

[30] Z. M. Odibat and N. T. Shawagfeh, "Generalized Taylor's formula," Applied Mathematics and Computation, vol. 186, no. 1, pp. 286-293, 2007.

[31] M. Caputo, "Linear models of dissipation whose Q is almost frequency independent-II," The Geophysical Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529-539, 1967.

Rezvan Ghoochani-Shirvan, (1) Jafar Saberi-Nadjafi (iD), (2) and Morteza Gachpazan (2)

(1) Department of Applied Mathematics, Ferdowsi University ofMashhad, International Campus, Mashhad, Iran

(2) Department of Applied Mathematics, School of Mathematical Sciences, Ferdowsi University ofMashhad, Mashhad, Iran

Correspondence should be addressed to Jafar Saberi-Nadjafi; najafi141@gmail.com

Received 17 August 2017; Revised 24 November 2017; Accepted 18 January 2018; Published 1 April 2018

Academic Editor: Jaume Gine

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Ghoochani-Shirvan, Rezvan; Saberi-Nadjafi, Jafar; Gachpazan, Morteza |

Publication: | International Journal of Differential Equations |

Article Type: | Report |

Date: | Jan 1, 2018 |

Words: | 3334 |

Previous Article: | Global Stability of an Economic Model with a Continuous Delay of Kaldor Type Modified. |

Next Article: | On the Control of Coefficient Function in a Hyperbolic Problem with Dirichlet Conditions. |

Topics: |