American Journal of Computational and Applied Mathematics

p-ISSN: 2165-8935    e-ISSN: 2165-8943

2012;  2(3): 101-104

doi: 10.5923/j.ajcam.20120203.06

Solving Partial Integro-Differential Equations Using Laplace Transform Method

Jyoti Thorwe, Sachin Bhalekar

Department of Mathematics, Shivaji University, Kolhapur, 416004, India

Correspondence to: Sachin  Bhalekar, Department of Mathematics, Shivaji University, Kolhapur, 416004, India.

Email:

Copyright © 2012 Scientific & Academic Publishing. All Rights Reserved.

Abstract

Partialintegro-differential equations (PIDE) occur naturally in various fields of science, engineering and social sciences. In this article, we propose a most general form of a linear PIDE with a convolution kernel. We convert the proposed PIDE to an ordinary differential equation (ODE) using a Laplace transform (LT). Solving this ODE and applying inverse LT an exact solution of the problem is obtained. It is observed that the LT is a simple and reliable technique for solving such equations. A variety of numerical examples are presented to show the performance and accuracy of the proposed method.

Keywords: Partialintegro-differential equations, Laplace transform

1. Introduction

Real life phenomena are often modelled by ordinary/partial differential equations. Due to the local nature of ordinary differential operator(ODO), the models containing merely ODOs do not help in modelling memory and hereditary properties. One of the best remedies to overcome this drawback is the introduction of integral term in the model. The ordinary/partial differential equation along with the weighted integral of unknown function gives rise to an integro-differential equation (IDE) or a partial integro-differential equation (PIDE) respectively. Analysis of such equations can be found in[1-4].
Applications of PIDEs can be found in various fields. Dehghan and shakeri[5] have used variational iteration method (VIM) to solve PIDEs arising in heat conduction of materials with memory. Various numerical schemes are proposed by Dehghan[6] to solve PIDEs arising in viscoelasticity. Nonlinear PIDEs arising in nuclear reactor dynamics are solved by Pao[7] and Pachapatte[8]. PIDEs have been used in jump-diffusion models for pricing of derivatives in finance[9]. Abergel[10] used a nonlinear PIDE in financial modelling. Hepperger[11] proposed a PIDE in the model of electricity swaptions. A PIDE governing biofluid flow in fractured biomaterials is proposed by Zadeh in[12].
The most promising tool for solving linear equations is the Laplace transform (LT) method[13,14]. LT is used in[16] for calculations of water flow and heat transfer in fractured rocks. Alquran et al.[17] used LT to solve non -homogeneous partial differential equations. Merdan et al. [18] proposed a new method for nonlinear oscillatory systems using LT.
Stiff systems of ODEs are solved by Aminikhah[19] using a combined LT and HPM. Kexue and Jiger[20] have utilized LT to solve problems arising in fractional differential equations.
In this article we propose a most general form of a linear PIDE in two independent variables with a convolution kernel. In Section 2 we provide some preliminaries regarding LT. Section 3 is devoted to the proposed method and Section 4 provides an ample number of examples of various types.

2. Preliminaries

2.1. Laplace Transform method:

Definition: TheLaplace transform of a function f(x), is defined by
,(whenever integral on RHS exists)where, x ≥0, p is real and ℒis the Laplace transform operator.Convolution Theorem:
If and then
where, .

3. Solving PIDEs using Laplace Transform Method

Consider PIDE,
(*)
(with prescribed conditions)
where are known functions. and c are constants or the functions of x.
Taking Laplace transform on both sides of PIDE (*) with respect to t we get,
(**)
where,
Equation (**) is an ordinary differential equation in . Solving this ordinary differential equation and taking inverse Laplace transform of , we get a solution of (*).

4. Illustrative examples

Figure 1. Solution of (1.1)
Example 1. Consider the PIDE
(1)
with initial conditions
(2)
and boundary condition
(3)
Taking Laplace transform with respect to t on both sides of (1),
(4)
Solution of (1) is
(5)
Where, C is a constant to be determined.
From (3),
(6)
Taking inverse Laplace transform on both the sides of (6), we get
The solution of (1) is plotted in the Fig.1.
Example 2. Consider the PIDE
(7)
with initial condition
(8)
and boundary condition
(9)
Taking Laplace transform w.r.t. t on both sides of (7),
(10)
Therefore the solution of (10) is
(11)
From the boundary condition (9)
(12)
Using (11) and (12), we get C = 0.
∴ Equation (11) becomes,
(13)
Taking inverse Laplace transform of (13)
(14)
The solution (14) is plotted in Figure 2.
Figure 2. Solution of
Example 3. Consider
(15)
(16)
(17)
Taking Laplace transform of (15) w.r.t. t we get
(18)
(19)
Solving (19) we get,
(20)
Now,
(22)
Using (21) and (22) in (20) we get,
(23)
(24)
Solving (23) and (24) we get,
Equation (20) becomes,
Taking inverse Laplace transform we get,
This is an exact solution of (15).
Figure 3 represents the graph of
Figure 3. Solution of (15)
Example 4. Consider
(25)
(26)
(27)
Taking Laplace transform of (25) with respect to t and using equation (26) we get,
(28)
Solving (28) we get,
(29)
Using (27) we get,
Therefore equation (29) becomes,
Taking inverse Laplace transform we get,
(30)
Figure 4 shows the graph of (30).
Figure 4. Solution of (25)

5. Conclusions

PIDEs are used in modelling various phenomena in science, engineering and social sciences. The LT technique is successfully used to solve a general linear PIDE involving a convolution kernel. We get exact solutions of such PIDEs after a few steps of calculations. We hope some other types of PIDEs and these equations can be used in modelling real life phenomena.

ACKNOWLEDGEMENTS

Authors are thankful to the anonymous referees for their insightful comments leading to the improvement of the paper.

References

[1]  Appell, J.M., Kalitvin,A.S.andZabrejko, P.P., Partial integral operators and integro-differential equations,M.Dekkar, New York(2000).
[2]  Bahuguna, D. and Dabas, J., Existence and uniqueness of a solution to a PIDE by the method of lines, Electronic Journal of Qualatative Theory of Differential Equations, 4(2008)1-12.
[3]  Pachapatte, B.G., On some new integral and integro-differntial inequalities in two independent variables and their applications, Journal of Differential Equations, 33(1979)249-272.
[4]  Yanik, E.G. and Fairweather, G., Finite element methods for parabolic and hyperbolic partial integro-differential equations, Nonlinear Analysis: Theory, Method and Applications, 12(1988)785-809.
[5]  Dehghan, M. and Shakeri, F., Solution of parabolic integro-differential equations arising in heat conduction in materials with memory via He’s variational iteration technique, International Journal For Numerical Methods In Biomedical Engineering, 26(2010)705-715.
[6]  Dehghan, M., Solution of a partial integro-differential equation arising from viscoelasticity, International Journal of Computational Mathematics,83(2006)123-129.
[7]  Pao, C.V., Solution of a nonlinear integro-differential system arising in Nuclear reactor dynamics, Journal Of Mathematical Analysis And Application, 48(1974)470-492.
[8]  Pachapatte, B.G., On a nonlinear diffusion system arising in reactor dynamics, Journal Of Mathematical Analysis And Applications, 94(1983)501-508.
[9]  Sachs, E.W. and Strauss, A.K., Efficient solution of a partial integro-differential equation in finance, Applied Numerical Mathematics, 58(2008)1687-1703.
[10]  Abeergel, F. and Tachet, R., A nonlinear partial integro-differential equation from mathematical finance, AIMS Journal, 10(2010)10-20.
[11]  Hepperger, P., Hedging electricity swaptions using partial integro-differential equations, Stochastic Processes And Their Applications, 122(2012)600-622.
[12]  Zadeh, K.S., An integro-partial differential equation for modeling biofluids flow in fractured biomaterials, Journal Of Theoretical Biology, 273(2011)72-79.
[13]  Debnath, L. and Bhatt, D., Integral transforms and their applications, CRC Press (2007).
[14]  Rehman, M., Integral equations and their applications,WIT Press(2007).
[15]  Schiff, J.L., The Laplace transform theory and applications, Springer, New York (1999).
[16]  Xiang, Tan-yong, Guo,Jia-qi, A Laplace transform and Green function method for calculation of water flow and heat transfer in fractured rocks, Rock And Soil Mechanics, 32(2)(2011)333-340.
[17]  Alquran, M.T., AL-khaled, K., Ali, M. and Taany, A., The combined Laplace transform- differential transform method for solving linear non-homogeneous partial differential equations, Journal Of Mathematics Computer science, 2(2012)690-701.
[18]  Merdan, M., Yildirim, A., Gokdogan, A. and Mohyud-dins, T., Coupling of homotopy perturbation, Laplace transform and pade Approximants for nonlinear oscillatory systems, World Applied Sciences Journal, 16(3)(2012)320-328.
[19]  Aminikhah, H., The combined Laplace transform and new homotopy perturbation methods for stiff systems of ordinary differential equations, Applied Mathematical Modeling, 36(2012)3638-3644.
[20]  Kexue, Li, Jigen, P., Laplace transform and fractional differential equations, Applied Mathematics Letters, 24(2011)2019-2023.