7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Editorial
Invited review
Letter to the Editor
Original Article
REVIEW
Review Article
SHORT COMMUNICATION
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Editorial
Invited review
Letter to the Editor
Original Article
REVIEW
Review Article
SHORT COMMUNICATION
View/Download PDF

Translate this page into:

Original article
31 (
1
); 110-113
doi:
10.1016/j.jksus.2017.07.006

The variational Adomian decomposition method for solving nonlinear two- dimensional Volterra-Fredholm integro-differential equation

Department of Mathematics Faculty of Science, King Abdul Aziz University, Jeddah, Saudi Arabia
Department of Mathematics Faculty of Science, King Khalid University, Abha, Saudi Arabia

⁎Corresponding author. malqrni@kku.edu.sa (M.M. Al-Qarni)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University.

Abstract

This paper outlines the coupling of variational iteration method (VIM) with Adomian decomposition method (ADM) for solving nonlinear mixed Volterra-Fredholm integro-diffierential equation (V-FIDE), this method is called variational Adomian decomposition method (VADM). Some numerical examples are introduced to verify that the method handles the difficulty of nonlinear term, red reduces the computational size and accelerates the convergence.

Keywords

Volterra-Fredholm integro-differential equation
Variational iteration method (VIM)
Variational Adomian decomposition method (VADM)
1

1 Introduction

In recent years, there has been a clear interest in the integro-differential equations which are a combination of differential and Volterra-Fredholm integral equations. Integro-differential equations play an important role in many branches of linear and non-linear functional analysis and their applications. The mentioned integro- differential equations are usually difficult to solve analytically, so approximation strategies are required to obtain the solution of the linear and nonlinear integro-differential equations (Huesin et al., 2008).

Several strategies are proposed to achieve this goal. The ADM is introduced by G. Adomian. It is based on the search for a solution in the form of a series and on decomposing the nonlinear operator into a series in which the terms are calculated recursively using the Adomian polynomials (Adomian, 1994). Many researchers studied and discussed the linear V-FIDE, Babolian et al. (2008) and Al-Jubory (2010) introduced some approximation method for solving Volterra-Fredholm integral and integro-differential equations. Dadkhah et al. (2010) utilized numerical solution of nonlinear V-FIDEs utilizing Legendre wavelets. Rabbani and Kiasodltani (2011) studied finding solution of nonlinear system of V-FIDE by utilizing discrete collocation method. Gherjalar and Mohammadikia, 2012 solved integral and integro-differential equations by utilizing B-splines function. The VIM (He, 2007, 1999) is a powerful device for solving various kinds of equations, linear and nonlinear. The technique has successfully been applied to many situations. For instance, He (2007) utilized the strategy to solve some integro-differential equations where he chose initial approximate solution in the form of exact solution with unknown constants.

In this work, the two-dimensional nonlinear V-FIDE is solved by the VADM. The two-dimensional nonlinear V-FIDE is given by

(1)
j = 0 k P j ( x 1 , ȷ 1 ) u j ( x 1 , ȷ 1 ) = f ̇ ( x 1 , ȷ 1 ) + a x 1 Ω F ( x 1 , ȷ 1 , y , τ ) γ ( u l ( y , τ ) ) dyd τ , ( x 1 , ȷ 1 ) J ́ = [ a , x 1 ] × Ω with initial conditions
(2)
u r ( a , ȷ 1 ) = g r , r = 0 , 1 , , k - 1 , Ω = [ a , b ]
where u j ( x 1 , ȷ 1 ) = d j u dx 1 j . The functions f ̇ ( x 1 , ȷ 1 ) , F ( x 1 , ȷ 1 , y , τ ) and γ ( u l ( y , τ ) ) , l 0 are analytic functions on J ́ , and functions P j ( x 1 , ȷ 1 ) , j = 0 , 1 , , k , P k ( x 1 , ȷ 1 ) 0 are given.

2

2 Adomian decomposition method

In what follows we display an outline for utilizing the ADM for solving the nonlinear V-FIDE. The Eq. (1) can be written as follows:

(3)
u x 1 , ȷ 1 = L - 1 f ̇ ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) + r = 0 k - 1 1 ( r ! ) ( x 1 - ȷ 1 ) r g r + L - 1 a x 1 Ω F ( x 1 , ȷ 1 , y , τ ) γ ( u l ( y , τ ) ) P k ( x 1 , ȷ 1 ) dyd τ - L - 1 j = 0 k - 1 P j ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) u j x 1 , ȷ 1 where L - 1 is the multiple integration operator, and Eq. (3) takes the form
(4)
u ( x 1 , ȷ 1 ) = L - 1 f ̇ ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) + r = 0 k - 1 1 ( r ! ) ( x 1 - ȷ 1 ) r g r + a x 1 Ω ( x 1 - ȷ 1 ) k ( k ) ! F ( x 1 , ȷ 1 , y , τ ) γ ( u l ( y , τ ) ) P k ( x 1 , ȷ 1 ) dyd τ - j = 0 k - 1 a x 1 ( x 1 - ȷ 1 ) k - 1 ( k - 1 ) ! P j ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) u j ( x 1 , ȷ 1 ) d ȷ 1
since j = 0 k - 1 L - 1 P j ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) u j ( x 1 , ȷ 1 ) = j = 0 k - 1 a x 1 ( x 1 - ȷ 1 ) k - 1 ( k - 1 ) ! P j ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) u j ( x 1 , ȷ 1 ) d ȷ 1

The ADM introduces the following expression

(5)
u ( x 1 , ȷ 1 ) = i = 0 u i ( x 1 , ȷ 1 ) for the solution u ( x 1 , ȷ 1 ) of Eq. (1) with initial conditions (2), where the components u i ( x 1 , ȷ 1 ) will be determined recurrently. In addition, the technique defines the nonlinear function γ ( u l ( y , τ ) ) by an infinite series of polynomials
(6)
γ ( u l ( y , τ ) ) = n = 0 Λ n
where Λ n , n 0 are defined by
(7)
Λ n = 1 n ! d n d α n γ i = 0 α i u i α = 0 , n = 0 , 1 , 2 ,

Substituting Eqs. (5) and (6) into Eq. (3) yields

(8)
i = 0 u i ( x 1 , ȷ 1 ) = L - 1 f ̇ ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) + r = 0 k - 1 1 ( r ! ) ( x 1 - ȷ 1 ) r g r + L - 1 a x 1 Ω F ( x 1 , ȷ 1 , y , τ ) P k ( x 1 , ȷ 1 ) i = 0 Λ i dyd τ - L - 1 j = 0 k - 1 P j ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) L i j

To determine the components u 0 , u 1 , u 2 , of the solution u ( x 1 , ȷ 1 ) , we set the recurrence relation

(9)
u 0 ( x 1 , ȷ 1 ) = L - 1 f ̇ ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) + r = 0 k - 1 1 ( r ! ) ( x 1 - ȷ 1 ) r g r , u i + 1 ( x 1 , ȷ 1 ) = L - 1 a x 1 Ω F ( x 1 , ȷ 1 , y , τ ) P k ( x 1 , ȷ 1 ) Λ i dyd τ - L - 1 j = 0 k - 1 P j ( x 1 , ȷ 1 ) P k ( x 1 , ȷ 1 ) L i j , i 0

The nonlinear terms given as follows: γ ( u l ( y , τ ) ) = i = 0 Λ i , D j ( u ( x 1 , ȷ 1 ) ) = i = 0 L i j . where D j = j u ( x 1 , ȷ 1 ) x 1 j is derivative operator .

3

3 Variational iteration method

According to the VIM (He, 2007; He, 1999; He, 2000; He and Wang, 2007; Nadjafi and Tamamgar, 2008; Sadigh Behzadi, 2012) we can write the correction functional for Eq. (1) as

(10)
u i + 1 ( x 1 , ȷ 1 ) = u i ( x 1 , ȷ 1 ) + 0 ȷ 1 λ ( ς ) [ u i k ( x 1 , ς ) - f ̇ ( x 1 , ς ) P k ( x 1 , ς ) - a x 1 Ω F ( x 1 , ς , y , τ ) P k ( x 1 , ς 1 ) γ ( u l ( y , τ ) ) dyd τ + j = 0 k - 1 P j ( x 1 , ς ) P k ( x 1 , ς ) u j ( x 1 , ς ) ] d ς where λ ( ς ) is a Lagrange multiplier which can be identified optimally via variational theory, u i is the ith approximate solution, and δ u i = 0 .

To find the optimal λ , we proceed as follows δ u i + 1 ( x 1 , ȷ 1 ) = δ u i ( x 1 , ȷ 1 ) + δ 0 ȷ 1 λ ( ς ) [ u i k ( x 1 , ς ) - f ̇ ( x 1 , ς ) P k ( x 1 , ς ) - a x 1 Ω F ( x 1 , ς , y , τ ) P k ( x 1 , ς 1 ) γ ( u l ( y , τ ) ) dyd τ + j = 0 k - 1 P j ( x 1 , ς ) P k ( x 1 , ς ) u j ( x 1 , ς ) ] d ς = 0

Then we apply the following stationary conditions 1 - λ = 0 , λ ( ς = ȷ 1 ) = 0 , λ = 0

This in turn gives λ = ς - ȷ 1

The solution is given by

(11)
u ( x 1 , ȷ 1 ) = lim i u i ( x 1 , ȷ 1 ) .

4

4 Variational Adomian decomposition method (VADM) for solving Eqs. (1) and (2)

This modified version of VADM which is obtained by the elegant coupling of VIM and ADM.

In VIM, from the correction functional for Eq. (1) we can write

(12)
u i + 1 ( x 1 , ȷ 1 ) = u i ( x 1 , ȷ 1 ) + 0 ȷ 1 λ ( ς ) u i k ( x 1 , ς ) - f ̇ ( x 1 , ς ) P k ( x 1 , ς ) - a x 1 Ω F ( x 1 , ς , y , τ ) P k ( x 1 , ς 1 ) γ ( u l ( y , τ ) ) dyd τ + j = 0 k - 1 P j ( x 1 , ς ) P k ( x 1 , ς ) u j ( x 1 , ς ) d ς

Substituting from Eq. (6) in Eq. (12)we have

(13)
u i + 1 ( x 1 , ȷ 1 ) = u i ( x 1 , ȷ 1 ) + 0 ȷ 1 λ ( ς ) u i k ( x 1 , ς ) - f ̇ ( x 1 , ς ) P k ( x 1 , ς ) - a x 1 Ω F ( x 1 , ς , y , τ ) P k ( x 1 , ς 1 ) Λ i dyd τ + j = 0 k - 1 P j ( x 1 , ς ) P k ( x 1 , ς ) u i j ( x 1 , ς ) d ς

Then the recurrence relation is given by

(14)
u 0 ( x 1 , ȷ 1 ) = r = 0 k - 1 1 ( r ! ) ( x 1 - ȷ 1 ) r g r , u i + 1 ( x 1 , ȷ 1 ) = u i ( x 1 , ȷ 1 ) + 0 ȷ 1 λ ( ς ) u i k ( x 1 , ς ) - f ̇ ( x 1 , ς ) P k ( x 1 , ς ) - a x 1 Ω F ( x 1 , ς , y , τ ) P k ( x 1 , ς 1 ) Λ i dyd τ + j = 0 k - 1 P j ( x 1 , ς ) P k ( x 1 , ς ) u i j ( x 1 , ς ) d ς , i 0

The approximate solution is given by

(15)
u ( x 1 , ȷ 1 ) = lim i u i ( x 1 , ȷ 1 ) .

To demonstrate the efficiency of the VADM we have considered the following examples.

5

5 Numerical examples

Example (5.1.) (Aghazadeh and Khajehnasiri, 2013; Poorfattah and Shaerlar, 2015; Darania et al., 2011): Consider the nonlinear integro - differential equation

(16)
2 u ( x 1 , ȷ 1 ) x 1 2 + sin ( x 1 ȷ 1 ) u ( x 1 , ȷ 1 ) - 0 x 1 0 1 x 1 ȷ 1 u ( y , τ ) x 1 dyd τ = f ( x 1 , ȷ 1 ) where
(17)
f ( x 1 , ȷ 1 ) = ( x 1 ȷ 1 ) sin ( x 1 ȷ 1 ) - 1 3 ȷ 1 3 .
with the initial conditions
(18)
I . Cs : u ( 0 , ȷ 1 ) = 0 , u ( 0 , ȷ 1 ) x 1 = ȷ 1 ,
which has exact solution u ( x 1 , ȷ 1 ) = x 1 ȷ 1 .

This example is solved by the VADM formule (14), the results are shown in (Table 1).

Table 1 Exact solution, approximate solution and Error by using VADM.
x 1 ȷ 1 Exact App.-VADM Err.-VADM
1.00E−02 1.00E−02 1.00000E−04 1.00041E−04 4.10000E−08
4.00E−02 4.00E−02 1.60000E−03 1.60245E−03 2.45000E−06
6.00E−02 6.00E−02 3.60000E−03 3.60790E−03 7.90000E−06
8.00E−02 8.00E−02 6.40000E−03 6.41774E−03 1.77400E−05
1.00E−01 1.00E−01 1.00000E−02 1.00326E−02 3.26000E−05

Example (5.2) (Aghazadeh and Khajehnasiri, 2013; Poorfattah and Shaerlar, 2015; Darania et al., 2011): Consider the nonlinear integro-differential equation

(19)
u ( x 1 , ȷ 1 ) 2 u ( x 1 , ȷ 1 ) ȷ 1 2 - 4 u ( x 1 , ȷ 1 ) 2 u ( x 1 , ȷ 1 ) x 1 2 + 4 0 x 1 0 1 u 2 ( y , τ ) dyd τ = f ̇ ( x 1 , ȷ 1 ) , x 1 [ 0 , 1 ] , where
(20)
f ̇ ( x 1 , ȷ 1 ) = x 1 - 1 2 π sin ( 2 π x 1 ) ȷ 1 - 1 4 π sin ( 4 π ȷ 1 ) .
with the boundary conditions
(21)
B . Cs : u ( 0 , ȷ 1 ) = u ( 1 , ȷ 1 ) = 0 ,
and the initial conditions
(22)
I . Cs : u ( x 1 , 0 ) = sin ( π x 1 ) , 0 x 1 1 u ( x 1 , 0 ) ȷ 1 = 0 , 0 x 1 1 .
which has exact solution u ( x 1 , ȷ 1 ) = sin ( π x 1 ) cos ( 2 π ȷ 1 ) .

This example is solved by the VADM formule (14), the results are written in (Table 2).

Table 2 Exact solution, approximate solution and Error by using VADM.
x 1 ȷ 1 Exact App.-VADM Err.-VADM
1.00E−03 1.00E−03 3.14152E−03 3.14300E−03 1.48000E−06
4.00E−03 4.00E−03 1.25621E−02 1.25333E−02 2.88000E−05
6.00E−03 6.00E−03 1.88350E−02 1.87290E−02 1.06000E−04
8.00E−03 8.00E−03 2.50984E−02 2.48510E−02 2.47400E−04
1.00E−02 1.00E−02 3.13487E−02 3.08400E−02 5.08700E−04

6

6 Conclusion

In this work, we applied the VADM for solving nonlinear mixed V-FIDE. The approximate solutions of V-FIDE are obtained by two powerful methods VIM and ADM in VADM. The given numerical examples showed the efficiency and accuracy of the introduced method, it reduces the size of computation without the restrictive assumption to handle nonlinear terms.

References

  1. , . Solving Frontier Problems of Physics: The Decomposition Method. Boston, MA: Kluwer Academic Publishers; .
  2. , , . Solving nonlinear two-dimensional volterra integro-differential equations by block-pulse functions. Math. Sci. 2013:1-6.
    [Google Scholar]
  3. , . Some Approximation Methods for Solving Volterra-Fredholm Integral and Integro-differential Equations. University of Technology; . (Ph.D. Thesis)
  4. Babolian, E., Masouri, Z., Hatamzadeh-Varmazyar, S., 2008. New Direct Method to Solve Nonlinear Volterra-Fredholm Integral and Integro-Differential Equation Using Operational Matrix with Block-pulse Functions. 8, 59–79.
  5. , , , . Numerical solution of nonlinear fredholm-volterra integro- differential equations using legendre wavelets. In: Proceedings of the 6th IMT-GT Conference on Mathematics, Statistics and its Applications (ICMSA). . p. :738-744.
    [Google Scholar]
  6. , , , . New computational method for solving some 2-dimensional nonlinear Volterra integro-differential equations. Numer. Algorithms. 2011;57:125-147.
    [Google Scholar]
  7. , , . Numerical solution of functional integral and integro-differential equations by using B-splines. Appl. Math.. 2012;3:1940-1944.
    [Google Scholar]
  8. , . Variational iteration method: a kind of nonlinear analytical technique, some examples. Int. J. Nonilnear Mech.. 1999;34(4):699-708.
    [Google Scholar]
  9. , . Variational iteration method for autonomous ordinary differential systems. Appl. Math. Comput.. 2000;114(2–3):115-123.
    [Google Scholar]
  10. , . Variational iteration method for eight order initial boundary value problem. Phys. Sci.. 2007;76:680-682.
    [Google Scholar]
  11. , , . Variational iteration method for solving integro-differential equations. Phys. Lett. A. 2007;367:188-191.
    [Google Scholar]
  12. , , , . Numerical solution of linear integro – differential equations. J. Math. Stat.. 2008;4(4):250-254.
    [Google Scholar]
  13. , , . The variational iteration method: a highly promising method for solving the system of integro-differential equations. Comput. Math. Appl.. 2008;56:346-351.
    [Google Scholar]
  14. , , . Direct method for solving nonlinear two-dimensional volterra-fredholm integro-differential equations by block-pulse functions. Int. J. Inf. Secur. Syst. Manage.. 2015;4:418-423.
    [Google Scholar]
  15. , , . Solving of nonlinear system of Volterra-Fredholm integro-differential equations by using discrete collocation method. J. Math. Comput. Sci.. 2011;5(4):282-289.
    [Google Scholar]
  16. , . The use of iterative methods to solve two-dimensional nonlinear volterra-fredholm integro-differential equations. ISPACS 2012:1-20.
    [Google Scholar]
Show Sections