A New Technique of Initial Boundary Value Problems Using Adomian Decomposition Method
A New Technique of Initial Boundary Value Problems Using Adomian Decomposition Method
A New Technique of Initial Boundary Value Problems Using Adomian Decomposition Method
n=0
(2.4)
Furthermore, suppose that the nonlinear term Nu can be written as infinite series in terms
of the Adomian polynomials A
n
of the form
Nu = A
n
n=0
, (2.S)
where the Adomian polynomials A
n
of Nu are evaluated using the formula
A
n
=
1
n!
J
n
Jz
n
N((z
))
=0
_
x=0
, n = u,1,2, (2.6)
Then substituting (2.4) and (2.5) in (2.3) gives
u
n
= +I
-1
(g) - I
-1
(R u
n
)
n=0
-I
-1
(A
n
)
n=0
n=0
(2.7)
Each term of series (2.4) is given by the recurrent relation
u
0
= + +I
-1
(g), n = u
u
n+1
= I
-1
(Ru
n
) - I
-1
(A
n
). n u
_ (2.8)
The convergence of the decomposition series have investigated by several authors
[Cherrualt Y. (1989), Cherrualt Y. and Adomian G. (1993)].
802 Elaf Jaafar Ali
3. Adomian decomposition method for solving initial boundary value
problems
To convey the basic idea for treatment of initial and boundary conditions by Adomian
decomposition method for solving initial boundary value problems, we consider the
following one dimensional differen-tial equation
F(u(x, t)) = g(x, t), u < x < 1, t > u, (S.1)
the initial conditions associated with (3.1) are of the form
u(x, u) =
0
(x),
ou(x, u)
ot
=
1
(x), u x 1, (S.2)
and the boundary conditions are given by
u(u, t) = q
0
(t), u(1, t) = q
1
(t), t > u, (S.S)
where
0
(x),
1
(x), q
0
(t) onJ q
1
(t) are given functions. The initial solution can be
written as
u
0
(x, t) =
0
(x) + t
1
(x). (S.4)
We construct a new successive initial solutions u
n
-
at every iteration for Eq. (3.1) by
applying a new technique
u
n
-
(x, t) = u
n
(x, t) + (1 -x)|q
0
(t) - u
n
(u, t)] + x |q
1
(t) -u
n
(1, t)], (S.S)
where n = u,1,2, , clearly that the new successive initial solutions u
n
-
of Eq. (3.1)
satisfying the initial and boundary conditions together as follows
i t = u tbcn u
n
-
(x, u) = u
n
(x, u),
i x = u tbcn u
n
-
(u, t) = q
0
(t),
i x = 1 tbcn u
n
-
(1, t) = q
1
(t). (S.6)
generally, sometimes a new technique of formula (3.5) can be applying for higher
dimensional initial boundary value problems by mixed initial and boundary conditions
together to have a new successive initial solutions u
n
-
as follows
Adomian decomposition method 803
Consider the k dimensional equation
F(u(x
1
, x
2
, , x
k
, t)) = g(x
1
, x
2
, , x
k
, t), u < x
1
, x
2
, , x
k
< 1, t > u, (S.7)
Then we have
u
n
-
(x
1
, x
2
, , x
k
, t)
= u
n
(x
1
, x
2
, , x
k
, t)
+(1 -x
1
)|q
01
(x
2
x
3
, , x
k
, t) - u
n
(u, x
2
, , x
k
, t)]
+x
1
|q
11
(x
2
x
3
, , x
k
, t) -u
n
(1, x
2
, , x
k
, t)] +
+(1 -x
k
)|q
0k
(x
1
x
2
, , x
k-1
, t) -u
n
(x
1
, x
2
, , x
k-1
, u, t)]
+x
k
|q
1k
(x
1
, x
2
, , x
k-1
, t) - u
n
(x
1
, x
2
, , x
k-1
, 1, t)], (S.8)
where n = u,1,2, , and the boundary conditions associated with Eq. (3.7) are of the
form
u(u, x
2
, , x
k
, t) = q
01
(x
2
, x
3
, , x
k
, t), u(1, x
2
, , x
k
, t) = q
11
(x
2
, x
3
, , x
k
, t),
u(x
1
, u, , x
k
, t) = q
02
(x
1
, x
3
, , x
k
, t), u(x
1
, 1, , x
k
, t) = q
12
(x
1
, x
3
, , x
k
, t),
.
u(x
1
, , x
k-1
, u, t) = q
0k
(x
1
, , x
k-1
, t), u(x
1
, , x
k-1
, 1, t) = q
1k
(x
1
, , x
k-1
, t),
(S.9)
and the initial conditions are given by
u(x
1
, x
2
, , x
k
, u) =
0
(x),
ou(x
1
, x
2
, , x
k
, u)
ot
=
1
(x). (S.1u)
Such as treatment of formula (3.8) when k = 2 is a very effective as shown in this paper.
4. Applications and results
Example 1: Consider the following one-dimensional heat-like problem [Momani
S.(2005)]
ou
ot
-
1
2
x
2
o
2
u
ox
2
= u, u < x < 1, t > u, (4.1)
804 Elaf Jaafar Ali
subject to the initial conditions
u(x, u) = x
2
, u < x < 1, (4.2)
and the boundary conditions
u(u, t) = u, u(1, t) =e
t
, t > u. (4.S)
By applying a new approximations u
n
-
in Eq. (3.5) we obtain
u
n
-
(x, t) = u
n
(x, t) + (1 -x)|u-u
n
(u, t)] +x|e
t
-u
n
(1, t)], (4.4)
where n = u,1,2, . The initial approximation is
u
0
(x, t) = u(x, u) = x
2
.
Now, we begin with a new initial approximationu
0
-
(wbcn n = u)
u
0
-
(x, t) = x
2
+ x (e
t
- 1). (4.S)
According to the Adomian decomposition method [Adomian G., 1994; Wazwaz A.M.,
2001], we have an operater form for Eq.(4.1) as
Iu =
1
2
x
2
o
2
u
ox
2
, u < x < 1, t > u, (4.6)
where the differential operator is I =
t
, so that I
-1
is integral operator
I
-1
(. ) = _(. )Jt. (4.7)
t
0
By operating with I
-1
on both sides of Eq.(4.6) and using a new technique of initial
solutions u
n
-
we have
u
n+1
(x, t) = __
1
2
x
2
o
2
u
n
-
(x, t)
ox
2
_Jt. (4.8)
t
0
By Eq. (4.5) we have
u
1
(x, t) = x
2
t ,
And by Eq. (4.4) we have consequently the following approximants
u
2
(x, t) =
x
2
t
2
2!
,
Adomian decomposition method 805
u
3
(x, t) =
x
2
t
3
S!
, (4.9)
and so on, the series solution is given by
u(x, t) = x
2
_1 + t +
t
2
2!
+
t
3
S!
+ _, (4.1u)
and in a closed form by
u(x, t) = x
2
c
t
. (4.11)
Which is an exact solution.
Example 2: We next consider the one-dimensional wave-like equation[Momani
S.(2005)]
o
2
u
ot
2
-
1
2
x
2
o
2
u
ox
2
= u, u < x < 1, t > u, (4.12)
subject to the initial conditions
u(x, u) = x,
ou(x, u)
ot
= x
2
, u < x < 1, (4.1S)
and the boundary conditions
u(u, t) = u , u(1, t) =1 +sinh t , t > u. (4.14)
By applying a new approximations u
n
-
in Eq. (3.5) we obtain
u
n
-
(x, t) = u
n
(x, t) + (1 -x)|u -u
n
(u, t)] +x|1 +sinht - u
n
(1, t)], (4.1S)
where n = u,1,2, . The initial approximation is
u
0
(x, t) = u(x, u) +
u(x,0)
t
t = x + x
2
t. (4.16)
Now, we begin with a new initial approximationu
0
-
(wbcn n = u)
u
0
-
(x, t) = x + x
2
t + x (sinht -t). (4.17)
By Adomian decomposition method , we have an operator form for Eq.(4.12) as
806 Elaf Jaafar Ali
Iu =
1
2
x
2
o
2
u
ox
2
, u < x < 1, t > u, (4.18)
where the differential operator is I =
2
t
2
, so that I
-1
is a two-fold integral operator
I
-1
(. ) = __(. )JtJt. (4.19)
t
0
t
0
By operating with I
-1
on both sides of (4.18) and using a new technique of initial
solutions u
n
-
we have
u
n+1
(x, t) = ___
1
2
x
2
o
2
u
n
-
(x, t)
ox
2
_JtJt ,
t
0
t
0
n u. (4.2u)
By Eq. (4.17) we have
u
1
(x, t) =
x
2
t
3
S!
,
And by Eq. (4.15) we have consequently the following approximants
u
2
(x, t) =
x
2
t
5
S!
,
u
3
(x, t) =
x
2
t
7
7!
, (4.21)
and so on, the series solution is given by
u(x, t) = x + x
2
_t +
t
3
S!
+
t
5
S!
+
t
7
7!
+_, (4.22)
and in a closed form by
u(x, t) = x + x
2
sinht. (4.2S)
Which is an exact solution.
Example 3: Consider the two-dimensional initial boundary value problem [Noor M.A.
and Mohyud-Din S.T.(2008)]:
o
2
u
ot
2
=
1
2
y
2
o
2
u
ox
2
+
1
2
x
2
o
2
u
oy
2
, u < x, y < 1, t > u, (4.24)
and the initial conditions
Adomian decomposition method 807
u(x, y, u) = x
2
+ y
2
,
ou(x, y, u)
ot
= -(x
2
+y
2
), (4.2S)
with boundary conditions
u(u, y, t) = y
2
c
-t
, u(1, y, t) = (1 + y
2
)c
-t
,
u(x, u, t) = x
2
c
-t
, u(x, 1, t) = (1 +x
2
)c
-t
. (4.26)
By applying a new approximations u
n
-
in Eq. (3.8) we obtain
u
n
-
(x, y, t) = u
n
(x, y, t) +(1 -x)|y
2
c
-t
-u
n
(u, y, t)] + x|(1 +y
2
)c
-t
-u
n
(1, y, t)] +
(1 - y)|x
2
c
-t
-u
n
(x, u, t)] +y|(1 +x
2
)c
-t
-u
n
(x, 1, t)], (4.27)
where n = u,1,2, . The initial approximation is
u
0
(x, y, t) = x
2
+y
2
- (x
2
+ y
2
)t. (4.28)
Now, we begin with a new initial approximationu
0
-
(wbcn n = u)
u
0
-
(x, y, t) = x
2
+y
2
- (x
2
+ y
2
)t + (1 -x)|y
2
c
-t
-y
2
+y
2
t]
+x|(1 +y
2
)c
-t
-(1 +y
2
) + (1 + y
2
)t] +(1 -y)|x
2
c
-t
-x
2
+x
2
t]
+y|(1 +x
2
)c
-t
-(x
2
+ 1) + (x
2
+ 1)t],
= -x -y +xc
-t
+ yc
-t
+x
2
c
-t
+y
2
c
-t
+ (x + y)t. (4.29)
According to the Adomian decomposition method, we have an operator form for
Eq.(4.24) as
Iu =
1
2
y
2
o
2
u
ox
2
+
1
2
x
2
o
2
u
oy
2
, u < x < 1, t > u, (4.Su)
where the differential operator is I =
2
t
2
, so that I
-1
is a two-fold integral operator
I
-1
(. ) = __(. )JtJt. (4.S1)
t
0
t
0
By operating with I
-1
on both sides of (4.30) and using a new technique of initial
solutions u
n
-
we have the following iteration formula
u
n+1
-
(x, y, t) = ___
1
2
y
2
o
2
u
n
-
(x, y, t)
ox
2
+
1
2
x
2
o
2
u
n
-
(x, y, t)
oy
2
_JtJt. (4.S2)
t
0
t
0
808 Elaf Jaafar Ali
By Eq. (4.29) we obtain a first iteration
u
1
(x, y, t) = -x
2
- y
2
+ (x
2
+ y
2
)t + (x
2
+ y
2
)c
-t
. (4.SS)
We can readily check
u(x, y, t) = u
0
(x, y, t) + u
1
(x, y, t) = (x
2
+y
2
)c
-t
. (4.S4)
Which yields an exact solution of Eq. (4.24).
Example 4: Consider the two-dimensional nonlinear inhomogeneous initial boundary
value problem [Noor M.A. and Mohyud-Din S.T. (2008)]:
o
2
u
ot
2
=
1S
2
(x u
xx
2
+ y u
2
) +2x
2
+ 2y
2
, u < x, y < 1, t > u, (4.SS)
and the initial conditions
u(x, y, u) = u,
ou(x, y, u)
ot
= u, (4.S6)
with boundary conditions
u(u, y, t) = y
2
t
2
+ yt
6
, u(1, y, t) = (1 +y
2
)t
2
+ (1 +y)t
6
,
u(x, u, t) = x
2
t
2
+xt
6
, u(x, 1, t) = (1 +x
2
)t
2
+(1 +x)t
6
. (4.S7)
By applying a new approximations u
n
-
in Eq. (3.8) we obtain
u
n
-
(x, y, t) = u
n
(x, y, t) +(1 -x)|y
2
t
2
+ yt
6
-u
n
(u, y, t)]
+x|(1 +y
2
)t
2
+ (1 + y)t
6
- u
n
(1, y, t)]
+(1 -y)|x
2
t
2
+ xt
6
-u
n
(x, u, t)]
+y|(1 +x
2
)t
2
+ (1 + x)t
6
-u
n
(x, 1, t)], n = u,1,2, (4.S8)
The initial approximation is
u
0
(x, y, t) = u(x, y, u) +
u(x,,0)
t
t = u.
Now, we begin with a new initial approximationu
0
-
(wbcn n = u)
u
0
-
(x, y, t) = (x
2
+y
2
)t
2
+ (1 -x)yt
6
+ x(1 +y)t
6
+ (1 -y)xt
6
+y(1 +x)t
6
, (4.S9)
Adomian decomposition method 809
By Adomian decomposition method, we have an operator form for Eq.(4.35) as
Iu =
1S
2
_x
o
2
u
ox
2
+ y
o
2
u
oy
2
_ + 2x
2
+2y
2
, u < x < 1, t > u, (4.4u)
where the differential operator is I =
2
t
2
, so that I
-1
is a two-fold integral operator
I
-1
(. ) = __(. )JtJt. (4.41)
t
0
t
0
By operating with I
-1
on both sides of (4.40) the components u
n
(x, y, t) can be
determind as
u
0
(x, y, t) = u(x, y, u) +
ou(x, y, u)
ot
t + I
-1
(2x
2
+2y
2
)
= __(2x
2
+ 2y
2
) Jt Jt = (x
2
+y
2
)t
2
,
t
0
t
0
(4.42)
u
n+1
(x, y, t) = ___
1S
2
_x
o
2
u
n
(x, y, t)
ox
2
+ y
o
2
u
n
(x, y, t)
oy
2
__JtJt ,
t
0
t
0
n u. (4.4S)
By a new technique of initial solutionsu
n
-
, we have the following iteration formula
u
n+1
-
(x, y, t) = ___
1S
2
_x
o
2
u
n
-
(x, y, t)
ox
2
+ y
o
2
u
n
-
(x, y, t)
oy
2
__JtJt. (4.44)
t
0
t
0
By Eq. (4.39) we obtain a first iteration
u
1
(x, y, t) = (x + y)t
6
. (4.4S)
We can readily check
u(x, y, t) = u
0
(x, y, t) + u
1
(x, y, t) = (x
2
+y
2
)t
2
+ (x + y)t
6
. (4.46)
Which yields an exact solution of Eq. (4.35).
Example 5: Consider the two-dimensional initial boundary value problem [Wazwaz A.M.
(2009)]:
810 Elaf Jaafar Ali
o
2
u
ot
2
=
1
2
x
2
o
2
u
ox
2
+
1
2
y
2
o
2
u
oy
2
, u < x, y < 1, t > u, (4.47)
and the initial conditions
u(x, y, u) = x
2
+ y
2
,
ou(x, y, u)
ot
= y
2
- x
2
, (4.48)
with boundary conditions
u(u, y, t) = y
2
c
t
, u(1, y, t) = c
-t
+y
2
c
t
,
u(x, u, t) = x
2
c
-t
, u(x, 1, t) = c
t
+x
2
c
-t
. (4.49)
By applying a new approximations u
n
-
in Eq. (3.8) we obtain
u
n
-
(x, y, t) = u
n
(x, y, t) +(1 -x)|y
2
c
t
-u
n
(u, y, t)] + x|c
-t
+y
2
c
t
-u
n
(1, y, t)] +
(1 - y)|x
2
c
-t
-u
n
(x, u, t)] +y|c
t
+ x
2
c
-t
-u
n
(x, 1, t)], (4.Su)
where n = u,1,2, . The initial approximation is
u
0
(x, y, t) = x
2
+y
2
+ (y
2
-x
2
)t. (4.S1)
Now, we begin with a new initial approximationu
0
-
(wbcn n = u)
u
0
-
(x, y, t) = x
2
+y
2
+ (y
2
-x
2
)t + (1 -x)|y
2
c
t
-y
2
- y
2
t] + x|c
-t
+ y
2
c
t
-1 -
y
2
-(y
2
- 1)t] +(1 -y)|x
2
c
-t
-x
2
+x
2
t] +y|c
t
+x
2
c
-t
-x
2
-1 - (1 -
x
2
)
t]. (4.S2)
According to the Adomian decomposition method, we have an operator form for
Eq.(4.47) as
Iu =
1
2
x
2
o
2
u
ox
2
+
1
2
y
2
o
2
u
oy
2
, u < x < 1, t > u, (4.SS)
where the differential operator is I =
2
t
2
, so that I
-1
is a two-fold integral operator
I
-1
(. ) = __(. )JtJt. (4.S4)
t
0
t
0
By operating with I
-1
on both sides of (4.53) and using a new technique of initial
solutions u
n
-
we have the following iteration formula
Adomian decomposition method 811
u
n+1
-
(x, y, t) = ___
1
2
x
2
o
2
u
n
-
(x, y, t)
ox
2
+
1
2
y
2
o
2
u
n
-
(x, y, t)
oy
2
_JtJt. (4.SS)
t
0
t
0
By Eq. (4.52) we obtain a first iteration
u
1
(x, y, t) = -x
2
- y
2
- (y
2
-x
2
)t + x
2
c
-t
+y
2
c
t
. (4.S6)
We can readily check
u(x, y, t) = u
0
(x, y, t) + u
1
(x, y, t) = x
2
c
-t
+ y
2
c
t
. (4.S7)
Which yields an exact solution of Eq. (4.47).
Example 6: Consider the two-dimensional initial boundary value problem [Wazwaz A.M.
(2009)]:
o
2
u
ot
2
=
1
2
x
2
o
2
u
ox
2
+
1
2
y
2
o
2
u
oy
2
, u < x, y < 1, t > u, (4.S8)
and the initial conditions
u(x, y, u) = y
2
,
ou(x, y, u)
ot
= x
2
, (4.S9)
with boundary conditions
u(u, y, t) = y
2
cosht , u(1, y, t) = sinht +y
2
cosht,
u(x, u, t) = x
2
sinht , u(x, 1, t) = x
2
sinh t +cosht. (4.6u)
By applying a new approximations u
n
-
in Eq. (3.8) we obtain
u
n
-
(x, y, t) =
u
n
(x, y, t) + (1 - x)|y
2
cosht -u
n
(u, y, t)] +x|sinh t +y
2
cosht -u
n
(1, y, t)] +
(1 - y)|x
2
sinht -u
n
(x, u, t)] + y|x
2
sinht +cosht -u
n
(x, 1, t)], (4.61)
where n = u,1,2, . The initial approximation is
u
0
(x, y, t) = y
2
+x
2
t. (4.62)
Now, we begin with a new initial approximationu
0
-
(wbcn n = u)
812 Elaf Jaafar Ali
u
0
-
(x, y, t) = y
2
+x
2
t + (1 - x)|y
2
cosht -y
2
] + x|sinht +y
2
cosht -y
2
- t]
+(1 -y)|x
2
sinht -x
2
t] +y|x
2
sinht +cosht -1 - x
2
t]. (4.6S)
According to the Adomian decomposition method, we have an operator form for
Eq.(4.58) as
Iu =
1
2
x
2
o
2
u
ox
2
+
1
2
y
2
o
2
u
oy
2
, u < x < 1, t > u, (4.64)
where the differential operator is I =
2
t
2
, so that I
-1
is a two-fold integral operator
I
-1
(. ) = __(. )JtJt. (4.6S)
t
0
t
0
By operating with I
-1
on both sides of (4.64) and using a new technique of initial
solutions u
n
-
we have the following iteration formula
u
n+1
-
(x, y, t) = ___
1
2
x
2
o
2
u
n
-
(x, y, t)
ox
2
+
1
2
y
2
o
2
u
n
-
(x, y, t)
oy
2
_JtJt. (4.66)
t
0
t
0
By Eq. (4.63) we obtain a first iteration
u
1
(x, y, t) = -y
2
-x
2
t + x
2
sinht +y
2
cosht. (4.67)
We can readily check
u(x, y, t) = u
0
(x, y, t) + u
1
(x, y, t) = x
2
sinht +y
2
cosht. (4.68)
Which yields an exact solution of Eq. (4.58).
5. Conclusions
In this paper, a very effective to construct a new initial successive solutions u
n
-
by
mixed initial and boundary conditions together which explained in formula (3.5) and
(3.8) when k = 2 which is used to find successive approximations u
n
of the solution u by
Adomian decomposition method to solve initial boundary value problems. These
technique to construct of a new successive initial solutions can give a more accurate
solution. Some examples (linear and nonlinear, homogeneous and nonhomogeneous
problems) are given in this paper to illustrate the effectiveness and convenience of a new
Adomian decomposition method 813
technique. It is important and obvious that the exact solutions have found directly from a
first iteration for the last four examples but if we use initial conditions only we will have
exact solution by calculating infinite successive solutions u
n
which closed form by Eq.
(2.4).
References
[1] A.M. Wazwaz, On the solution of the fourth order parabolic equation by the
decomposition method, Intern. J . ComputerMath., 57 (1995)213-217.
[2] A.M. Wazwaz, The decomposition method for solving higher dimensional initial
boundary value problems of variable coefficients, International J ournal of
ComputerMathematics 76 (2) (2000), 159-172.
[3] A.M. Wazwaz, Analytic treatment for variable coefficient fourth-order parabolic
partial differential equations, Applied Mathematics and Computation, 123 (2) (2001)
219-227.
[4] A.M. Wazwaz, Partial differential equations and solitary waves theory, Higher
education press, Beijing and Springer-Verlag Berlin Heidelberg (2009).
[5] A. Shidfar and A. Molabahrami, A weighted algorithm based on the homotopy
analysis method: Application to inverse heat conduction problems, Commun Nonlinear
Sci Numer Simulat 15 (2010), 2908-2915.
[6] C. Chun and R. Sakthivel, Homotopy perturbation technique for solving two-point
boundary value problems comparison with other methods, Computer Physics
Communications 181 (2010), 1021-1024.
[7] G. Adomian, Solving Frontier Problems on Physics: The Decomposition Method,
Kluwer Academic Publisher, Boston, 1994.
[8] G. Wu and E. W Lee, Fractional variational iteration method and its application,
Physics Letters A 374 (2010) 2506-2509.
[9] J . Biazar and H. Aminikhah, Exact and numerical solutions for non-linear Burger's
equation by VIM, Mathematical and Computer Modelling 49 (2009), 1394-1400.
[10] M.A. Noor and S.T. Mohyud-Din, Solving Higher Dimensional Initial Boundary
Value Problems by Variational Iteration Decomposition Method, Applications and
Applied Mathematics An International J ournal 3(2) (2008), 254-266.
[11] R. Yulita Molliq, M.S.M. Noorani, I. Hashim, Variational iteration method for
fractional heat- and wave-like equations, Nonlinear Analysis: Real World Applications
10 (2009) 1854-1869.
814 Elaf Jaafar Ali
[12] S. Momani, Analytical approximate solution for fractional heat-like and wave-like
equations with variable coefficients using the decomposition method. Applied
Mathematics and Computation 165 (2005) 459472.
[13] Y. Cherrualt, Convergence of Adomians method, Kybernetes 18 (1989) 31-38.
[14] Y. Cherrualt and G. Adomian, Decomposition methods: a new proof of
convergence, Math. Comput. Model. 18 (1993) 103-106.
[15] Z. Niu and C. Wang, A one-step optimal homotopy analysis method for nonlinear
differential equations, Commun Nonlinear Sci Numer Simulat 15 (2010), 2026-2036.
Received: September, 2011