# Modelling the Nonlinear Wave Motion within the Scope of the Fractional Calculus

Authors: Abdon Atangana; Seyma Tuluce Demiray; Hasan Bulut
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

The aim of this paper was to first extend the model describing the nonlinear wave movement to the concept of noninteger order derivative. The extended equation was investigated within the scope of an iterative method. The stability and convergence analysis of the iteration method for this extended equation was presented in detail. The uniqueness of the special solution was also investigated. A resume of the method for solving this equation was provided. The algorithm was used to derive the unique special solution for given initial conditions.


---

## Body

## 1. Introduction

The Kadomtsev-Petviashvili equation (or basically the KP equation) is a nonlinear partial differential equation in two spatial and one temporal coordinates that illustrate the development of nonlinear, extended waves of diminutive amplitude with sluggish reliance on the crosswise coordinate. There are two dissimilar descriptions of the KP equation. Nevertheless the version we are interesting in is given below as(1)
∂
x
(
∂
t
u
+
u
∂
x
u
+
ε
2
∂
x
x
x
u
)
+
λ
∂
y
y
u
=
0
.
The case λ
=
1 is known as the KPII equation and models, for instance, water waves with small surface tension. The case λ
=
i is known as the KPI equation and may be used to model waves in thin films with high surface tension [1–4]. The equation is frequently in black and white with dissimilar coefficients in front of the various terms, but the particular values are inessential, since they can be modified by appropriately rescaling the dependent and independent variables [5]. The KP equation is a worldwide integrable structure in two spatial dimensions in the similar line of attack that the KdV equation can be looked upon as a widespread integrable system in one spatial dimension, since many other integrable systems can be obtained as reductions [5].Per se, the KP equation has been comprehensively considered in the mathematical society in the last forty years. The KP equation is also one of the generally common models in nonlinear wave theory, which arises as a reduction of system with quadratic nonlinearity, which admits weakly dispersive waves, in a paraxial wave approximation [3]. The equation unsurprisingly comes into view as a well-known boundary in the asymptotic explanation of such systems in which merely the most important order terms are engaged and an asymptotic equilibrium flanked by weak scattering, quadratic nonlinearity and diffraction is supposed. The poles separately acted by the two spatial variables accounts for the asymmetric way in that they appear in the equation [4].In the last decade, many physical problems have been with great success described within differential equations with both integer order and fractional order. More importantly it was observed most of the time that, all physical problems described or modelled within the folder of the fractional order derivative were more predictable than the conventional order derivative, for instance, in the study of groundwater problems [6–10] and others [11–15]. One aspect of this study is to revert the conventional KP equation to the concept of fractional order derivative. However one of the big problems encountered in modelling real world problems is perhaps to find the solution of these nonlinear equations, especially when it came to fractional differential equation. Another aspect of this paper will therefore be devoted to the derivation of approximation of the generalized KP equation; we will also present the convergence and the uniqueness of the special solution. The generalised equation under study here is given as
(2)
∂
t
x
α
+
1
u
+
(
∂
x
u
)
2
+
u
∂
x
x
u
+
ε
2
∂
x
x
x
x
4
u
+
λ
∂
y
y
2
u
=
0
,
1111111111111111111111111111111
0
<
α
≤
1
.
We will present some basic formulas of the noninteger derivative in the following sections.


## 2. Basic Knowledge about Fractional Derivatives

One can find nowadays in the literature different definitions of fractional derivatives. However the furthermost well-liked ones are the Riemann-Liouville and the Caputo derivatives. For Caputo we have(3)
0
C
D
x
α
(
f
(
x
)
)
=
1
Γ
(
n
-
α
)
∫
0
x
(
x
-
t
)
n
-
α
-
1
d
n
f
(
t
)
d
t
n
d
t
.
For the case of Riemann-Liouville we have the following definition:
(4)
D
x
α
(
f
(
x
)
)
=
1
Γ
(
n
-
α
)
d
n
d
x
n
∫
0
x
(
x
-
t
)
n
-
α
-
1
f
(
t
)
d
t
.
Each one of these fractional derivatives presents some compensations and weakness [16–18].Definition 1.
LetΩ
=
[
a
,
b
]
(
-
∞
≤
a
<
b
≤
∞
) be a finite or infinite interval of real axis R
=
(
-
∞
,
∞
). We denote by L
p
(
a
,
b
)
(
1
≤
p
≤
∞
) the set of those Lebesgue complex-valued measurable functions f on Ω for which ∥
f
∥
p
<
∞, where
(5)
∥
f
∥
p
=
(
∫
a
b
|
f
(
t
)
|
p
d
t
)
1
/
p
(
1
≤
p
≤
∞
)
.
We will, in addition to the above definition, present the following useful theorem.Theorem 2.
Ifh
(
t
)
∈
L
1
(
R
) and h
1
(
t
)
∈
L
p
(
R
), then their convolution (
h
*
h
1
)
(
x
)
∈
L
p
(
R
)
(
1
≤
p
≤
∞
), and the following inequality holds [16]:
(6)
∥
f
(
h
*
h
1
)
(
x
)
∥
p
<
∥
h
∥
1
∥
h
1
∥
p
.
In particular, if h
(
t
)
∈
L
1
(
R
) and h
1
(
t
)
∈
L
2
(
R
), then their convolution (
h
*
h
1
)
(
x
)
∈
L
2
(
R
); then,
(7)
∥
f
(
h
*
h
1
)
(
x
)
∥
2
<
∥
h
∥
1
∥
h
1
∥
2
.Lemma 3 (see [16]).
The fractional integration operators withR
(
α
)
>
0 are bounded in L
p
(
a
,
b
)
(
1
≤
p
≤
∞
):
(8)
∥
I
a
α
f
∥
p
≤
K
∥
f
∥
p
,
K
=
(
b
-
a
)
R
(
α
)
R
(
α
)
|
Γ
(
α
)
|
.


## 3. Derivation of the Special Solution

The real problem with nonlinear equation is to find a suitable analytical method that can be used to derive their exact or special solutions. No wonder, many scholars have devoted their attention in developing methods to solve these equations analytically or numerically, but when the partial differential equation with mixed derivative is concerned no well-accurate approach has been proposed. We will in this paper propose an iteration method to derive the fractional Kadomtsev-Petviashvili equation.The method proposed here consists of applying first the inverse operator of∂
t
x
α
+
1 on both sides of (2) to obtain the following:
(9)
u
(
x
,
y
,
t
)
=
u
(
x
,
y
,
0
)
+
u
(
0
,
y
,
t
)
-
u
(
0
,
y
,
0
)
+
1
Γ
(
α
)
∫
0
x
∫
0
t
(
t
-
l
)
α
-
1
(
(
∂
v
u
)
2
+
u
∂
v
v
u
0000000000000000000000000
+
ε
2
∂
v
v
v
v
4
u
+
λ
∂
y
y
2
u
)
d
v
d
l
.
To solve the above equation, we assume that its solution is in series form as
(10)
u
(
x
,
y
,
t
)
=
∑
n
=
0
∞
u
n
(
x
,
y
,
t
)
.
Secondly we introduce an artificial parameter according to the concept of the homotopy, now after substituting (10) into (9) and putting together terms of same power of the artificial parameter p, we arrive at the following equations:
(11)
u
0
(
x
,
y
,
t
)
=
u
(
x
,
y
,
0
)
+
u
(
0
,
y
,
t
)
-
u
(
0
,
y
,
0
)
u
1
(
x
,
y
,
t
)
=
1
Γ
(
α
)
∫
0
x
∫
0
t
(
t
-
l
)
α
-
1
(
(
∂
v
u
0
)
2
+
u
∂
v
v
u
0
00000000000000000000000
+
ε
2
∂
v
v
v
v
4
u
0
+
λ
∂
y
y
2
u
0
)
d
v
d
l
,
u
n
(
x
,
y
,
t
)
=
1
Γ
(
α
)
∫
0
x
∫
0
t
(
t
-
l
)
α
-
1
(
H
n
-
1
1
+
H
n
-
1
2
+
ε
2
∂
v
v
v
v
4
u
n
-
1
00000000000000000000000
+
λ
∂
y
y
2
u
n
-
1
)
d
v
d
l
,
0
n
≥
1
,
H
n
-
1
1
=
∑
k
=
0
n
-
1
∂
v
u
k
∂
v
u
n
-
k
-
1
,
H
n
-
1
2
=
∑
k
=
0
n
-
1
∂
v
u
k
u
n
-
k
-
1
.
The iteration formula in (11) can be used to derive all terms providing that the initial conditions are given. We will put the above procedure in a resume usually called algorithm to help readers of the numerical method field implement method via computer.Algorithm 4.
Consider(i)
input:I
(
x
,
y
,
t
)
=
u
(
x
,
y
,
0
)
+
u
(
0
,
y
,
t
)
-
u
(
0
,
y
,
0
) as early boarder,
(ii)
k—number terms in the series calculation,
(iii)
output:u
part
(
x
,
y
,
t
), the special solution.
Step
1. Put u
0
(
x
,
y
,
t
)
=
I
(
x
,
y
,
t
) and u
part
(
x
,
y
,
t
)
=
u
0
(
x
,
y
,
t
).
Step
2. For k
=
0 to n
-
1 do Steps 3, 4, and 5.
Step
3. Compute
(12)
H
n
-
1
1
=
∑
k
=
0
n
-
1
∂
v
u
k
∂
v
u
n
-
k
-
1
,
H
n
-
1
2
=
∑
k
=
0
n
-
1
∂
v
u
k
u
n
-
k
-
1
,
b
n
=
1
Γ
(
α
)
×
∫
0
x
∫
0
t
(
t
-
l
)
α
-
1
(
H
n
-
1
1
+
H
n
-
1
2
+
ε
2
∂
v
v
v
v
4
u
n
-
1
00000000000000000
+
λ
∂
y
y
2
u
n
-
1
)
d
v
d
l
,
n
≥
1
.
Step
4. Compute
(13)
u
n
+
1
(
x
,
t
)
=
v
n
+
u
n
.
Step
5. Compute
(14)
u
part
(
x
,
y
,
t
)
=
u
part
(
x
,
y
,
t
)
+
u
n
+
1
(
x
,
y
,
t
)
.
Stop.


## 4. Convergence and Uniqueness Analysis of the Proposed Method

In recent years, many papers have been published in which iterations methods are used to give approximate solution to the nonlinear equations. However in these papers, no study of stability, the convergence, and the uniqueness analysis has been done. This has become nowadays a routine whereas the study of convergence is very important and very difficult because the proof of convergence shows the strength of the method. We therefore devote this section to the study of the convergence of the scheme used here and applied to the fractional nonlinear wave motion. To achieve this we will think about the following fractional sub-Hilbert spaceH
α of the Hilbert space H
=
L
2
(
(
a
,
b
)
×
(
0
,
T
)
) [19] that can be defined as the set of those functions:
(15)
v
:
(
a
,
b
)
×
[
0
,
T
]
⟶
R
,
1
Γ
(
α
)
∫
∫
(
t
-
l
)
α
-
1
v
2
d
l
d
s
<
∞
.
We correspondingly undertake that the differential operatives are circumscribed under the L
2 norms. We consider the fractional Kadomtsev-Petviashvili equation in the light of the above assumptions and then the operator form of the equation is
(16)
B
(
u
)
=
-
∂
t
x
α
+
1
(
u
)
=
(
∂
x
u
)
2
+
u
∂
x
x
u
+
ε
2
∂
x
x
x
x
4
u
+
λ
∂
y
y
2
u
.
The homotopy decomposition method used here is convergent if the shadowing two suggestions are fulfilled (see [20] and references therein).(H1)
(
B
(
u
)
-
B
(
v
)
,
u
-
v
)
≤
q
∥
u
-
v
∥,  q
>
0  for all  (
u
,
v
)
∈
H
α.
(H2)
For any positive constantM
>
0 there exists a constant P
>
0 such that, for u
,
v
∈
H with ∥
u
∥
≤
M,  ∥
v
∥
≤
M.We have
(17)
(
B
(
u
)
-
B
(
v
)
,
w
)
≤
P
∥
u
-
v
∥
∥
w
∥
∀
w
∈
H
α
.
We will present some properties of the inner product.


### 4.1. Properties of the Inner Product [21, 22]

The commentary on Hilbert space entertains numerous illustrations of the inner product spaces wherein the metric produced by the inner product profits a complete metric space. Inner product spaces have an instinctively outlined norm grounded upon the inner product of the space itself that does not the parallelogram equality:(18)
∥
x
∥
=
(
x
,
x
)
.
It is well defined by the no negativity axiom of the definition of the inner product space. The following properties can be observed:
(19)
|
(
x
,
y
)
|
≤
∥
x
∥
∥
y
∥
.
The above is the well-known Cauchy-Schwarz inequality. Also the following can be obtained:
(20)
∥
a
·
x
∥
=
|
a
|
·
∥
x
∥
.
The above is called homogeneity. The last interesting one for this paper will be given as
(21)
∥
x
+
y
∥
≤
∥
x
∥
+
∥
y
∥
.
The above is called triangular inequality.Proof.
Using the definition of our operatorB, we have the following:
(22)
B
(
u
)
-
B
(
v
)
=
(
∂
x
u
)
2
-
(
∂
x
v
)
2
+
u
∂
x
x
u
-
v
∂
x
x
v
+
ε
2
∂
x
x
x
x
4
(
u
-
v
)
+
λ
∂
y
y
2
(
u
-
v
)
=
∂
x
(
u
∂
x
u
-
v
∂
x
v
)
+
ε
2
∂
x
x
x
x
4
(
u
-
v
)
+
λ
∂
y
y
2
(
u
-
v
)
=
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
+
ε
2
∂
x
x
x
x
4
(
u
-
v
)
+
λ
∂
y
y
2
(
u
-
v
)
.
With the above reduction in hand, we can now proceed to compute
(23)
(
B
(
u
)
-
B
(
v
)
,
u
-
v
)
=
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
+
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
+
(
λ
∂
y
y
2
(
u
-
v
)
,
u
-
v
)
.
We will consider case by case the components of the above equation:
(24)
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
.
Regarding the bodily problem under examination, we are required to point out that the wave can only be propagated in finite space; therefore, u, v are bordered, as a result, we can find a positive constant N such that (
u
,
u
)
,
(
v
,
v
)
<
N
2. It follows by making use of Schwartz inequality that
(25)
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
≤
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
∥
u
-
v
∥
.
Using the properties of the norm and the derivative, we can possibly find a positive real number, say, for instance, O
1, satisfying
(26)
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
≤
O
1
∥
∂
x
u
2
-
∂
x
v
2
2
∥
.
Using the same properties, we can get hold of a second positive real constant, say, for instance, O
2, that will allow us to have
(27)
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
≤
O
1
O
2
2
∥
u
2
-
v
2
∥
.
We can now further extend this to the following:
(28)
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
≤
O
1
O
2
2
∥
u
-
v
∥
∥
u
+
v
∥
.
Now we can use the fact that the nonlinear wave is propagating in the finite space to obtain ∥
u
-
v
∥
∥
u
+
v
∥
≤
4
N
4 such that
(29)
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
≤
4
N
4
O
1
O
2
2
∥
u
-
v
∥
.
We will continue our investigation with the following case:
(30)
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
.
Without doubt, we can find four position real numbers, say, for instance, O
3, O
4, O
5 and O
4, allowing us to have the following relation:
(31)
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
≤
ε
2
O
3
O
4
O
5
O
4
∥
u
-
v
∥
∥
u
-
v
∥
.
We can further use the fact that (
u
,
u
)
,
(
v
,
v
)
<
N
2 to obtain
(32)
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
≤
2
N
2
ε
2
O
3
O
4
O
5
O
4
∥
u
-
v
∥
.
And also
(33)
(
λ
∂
y
y
2
(
u
-
v
)
,
u
-
v
)
≤
2
N
2
λ
O
6
O
7
∥
u
-
v
∥
.
However, by inserting (33), (32), and (29) into (23), we arrive at the following result:
(34)
(
B
(
u
)
-
B
(
v
)
,
u
-
v
)
≤

(
2
N
2
ε
2
O
3
O
4
O
5
O
4

+

2
N
2
λ
O
6
O
7

+

4
N
4
O
1
O
2
2
)
∥
u
-
v
∥
.
It is obvious that we can take
(35)
q
=
2
N
2
ε
2
O
3
O
4
O
5
O
4
+
2
N
2
λ
O
6
O
7
+
4
N
4
O
1
O
2
2
.
And hypothesis (H1) is verified.
It is also easier to prove hypothesis two just by realizing that(36)
(
B
(
u
)
-
B
(
v
)
,
w
)
≤
P
∥
u
-
v
∥
∥
w
∥
,
P
=
2
M
2
O
8
O
9
2
+
ε
2
O
10
O
11
O
12
O
13
+
λ
O
14
O
15
.
And then hypothesis (H2) is also verified. Therefore following with hypotheses (H1) and (H2) verified, we can without reservation establish the following theorem.Theorem 5.
Let us consider(37)
B
(
u
)
=
(
∂
x
u
)
2
+
u
∂
x
x
u
+
ε
2
∂
x
x
x
x
4
u
+
λ
∂
y
y
2
u
and consider the initial and boundary condition for (2); then the proposed method leads to a special solution of (2).The proof of the above theorem follows directly from hypotheses (H1) and (H2) and also Lemma3.Theorem 6.
Taking into account the initial conditions for (2), then the special solution of (2) u
e
s
p to which u converges is unique.


## 4.1. Properties of the Inner Product [21, 22]

The commentary on Hilbert space entertains numerous illustrations of the inner product spaces wherein the metric produced by the inner product profits a complete metric space. Inner product spaces have an instinctively outlined norm grounded upon the inner product of the space itself that does not the parallelogram equality:(18)
∥
x
∥
=
(
x
,
x
)
.
It is well defined by the no negativity axiom of the definition of the inner product space. The following properties can be observed:
(19)
|
(
x
,
y
)
|
≤
∥
x
∥
∥
y
∥
.
The above is the well-known Cauchy-Schwarz inequality. Also the following can be obtained:
(20)
∥
a
·
x
∥
=
|
a
|
·
∥
x
∥
.
The above is called homogeneity. The last interesting one for this paper will be given as
(21)
∥
x
+
y
∥
≤
∥
x
∥
+
∥
y
∥
.
The above is called triangular inequality.Proof.
Using the definition of our operatorB, we have the following:
(22)
B
(
u
)
-
B
(
v
)
=
(
∂
x
u
)
2
-
(
∂
x
v
)
2
+
u
∂
x
x
u
-
v
∂
x
x
v
+
ε
2
∂
x
x
x
x
4
(
u
-
v
)
+
λ
∂
y
y
2
(
u
-
v
)
=
∂
x
(
u
∂
x
u
-
v
∂
x
v
)
+
ε
2
∂
x
x
x
x
4
(
u
-
v
)
+
λ
∂
y
y
2
(
u
-
v
)
=
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
+
ε
2
∂
x
x
x
x
4
(
u
-
v
)
+
λ
∂
y
y
2
(
u
-
v
)
.
With the above reduction in hand, we can now proceed to compute
(23)
(
B
(
u
)
-
B
(
v
)
,
u
-
v
)
=
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
+
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
+
(
λ
∂
y
y
2
(
u
-
v
)
,
u
-
v
)
.
We will consider case by case the components of the above equation:
(24)
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
.
Regarding the bodily problem under examination, we are required to point out that the wave can only be propagated in finite space; therefore, u, v are bordered, as a result, we can find a positive constant N such that (
u
,
u
)
,
(
v
,
v
)
<
N
2. It follows by making use of Schwartz inequality that
(25)
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
≤
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
∥
u
-
v
∥
.
Using the properties of the norm and the derivative, we can possibly find a positive real number, say, for instance, O
1, satisfying
(26)
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
≤
O
1
∥
∂
x
u
2
-
∂
x
v
2
2
∥
.
Using the same properties, we can get hold of a second positive real constant, say, for instance, O
2, that will allow us to have
(27)
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
≤
O
1
O
2
2
∥
u
2
-
v
2
∥
.
We can now further extend this to the following:
(28)
∥
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
∥
≤
O
1
O
2
2
∥
u
-
v
∥
∥
u
+
v
∥
.
Now we can use the fact that the nonlinear wave is propagating in the finite space to obtain ∥
u
-
v
∥
∥
u
+
v
∥
≤
4
N
4 such that
(29)
(
∂
x
(
∂
x
u
2
-
∂
x
v
2
2
)
,
u
-
v
)
≤
4
N
4
O
1
O
2
2
∥
u
-
v
∥
.
We will continue our investigation with the following case:
(30)
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
.
Without doubt, we can find four position real numbers, say, for instance, O
3, O
4, O
5 and O
4, allowing us to have the following relation:
(31)
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
≤
ε
2
O
3
O
4
O
5
O
4
∥
u
-
v
∥
∥
u
-
v
∥
.
We can further use the fact that (
u
,
u
)
,
(
v
,
v
)
<
N
2 to obtain
(32)
(
ε
2
∂
x
x
x
x
4
(
u
-
v
)
,
u
-
v
)
≤
2
N
2
ε
2
O
3
O
4
O
5
O
4
∥
u
-
v
∥
.
And also
(33)
(
λ
∂
y
y
2
(
u
-
v
)
,
u
-
v
)
≤
2
N
2
λ
O
6
O
7
∥
u
-
v
∥
.
However, by inserting (33), (32), and (29) into (23), we arrive at the following result:
(34)
(
B
(
u
)
-
B
(
v
)
,
u
-
v
)
≤

(
2
N
2
ε
2
O
3
O
4
O
5
O
4

+

2
N
2
λ
O
6
O
7

+

4
N
4
O
1
O
2
2
)
∥
u
-
v
∥
.
It is obvious that we can take
(35)
q
=
2
N
2
ε
2
O
3
O
4
O
5
O
4
+
2
N
2
λ
O
6
O
7
+
4
N
4
O
1
O
2
2
.
And hypothesis (H1) is verified.
It is also easier to prove hypothesis two just by realizing that(36)
(
B
(
u
)
-
B
(
v
)
,
w
)
≤
P
∥
u
-
v
∥
∥
w
∥
,
P
=
2
M
2
O
8
O
9
2
+
ε
2
O
10
O
11
O
12
O
13
+
λ
O
14
O
15
.
And then hypothesis (H2) is also verified. Therefore following with hypotheses (H1) and (H2) verified, we can without reservation establish the following theorem.Theorem 5.
Let us consider(37)
B
(
u
)
=
(
∂
x
u
)
2
+
u
∂
x
x
u
+
ε
2
∂
x
x
x
x
4
u
+
λ
∂
y
y
2
u
and consider the initial and boundary condition for (2); then the proposed method leads to a special solution of (2).The proof of the above theorem follows directly from hypotheses (H1) and (H2) and also Lemma3.Theorem 6.
Taking into account the initial conditions for (2), then the special solution of (2) u
e
s
p to which u converges is unique.


## 5. Derivation of the Special Solution

We make use of Algorithm4 to derive the special solution for (2). Here we assume that the initial condition is given as
(38)
u
(
x
,
y
,
0
)
=
x
y
.
Then using the algorithm we obtain the following:
(39)
u
1
(
x
,
y
,
t
)
=
t
α
(
-
Sin
⁡
[
2
y
]
+
Sin
⁡
[
2
(
x
+
y
)
]
00000000000000
+
2
(
x
+
Sin
⁡
[
x
]
Sin
⁡
[
x
+
2
y
]
)
)
00000000
000
×
(
4
α
Gamma
[
α
]
)
-
1
,
u
2
(
x
,
y
,
t
)
=
1
24
Gamma
[
1
+
α
]
2
t
2
α
×
(
-
3
cos
⁡
⁡
[
x
-
y
]
+
12
cos
⁡
⁡
[
y
]
+
4
cos
⁡
⁡
[
3
y
]
0000000
-
9
cos
⁡
⁡
[
x
+
y
]
-
7
cos
⁡
⁡
[
3
(
x
+
y
)
]
0000000
+
3
cos
⁡
⁡
[
x
+
3
y
]
+
3
sin
⁡
[
x
-
y
]
0000000
-
36
sin
⁡
[
y
]
-
4
sin
⁡
[
3
y
]
+
39
sin
⁡
[
x
+
y
]
0000000
+

12
x
(

-

2
cos
⁡
⁡
[
2
y
]

+

2
sin
⁡
[
2
y
]

+

sin
⁡
[
x

+

y
]
)
0000000
+
sin
⁡
[
3
(
x
+
y
)
]
0000000
+
72
sin
⁡
[
x
]
(
-
cos
⁡
⁡
[
x
+
2
y
]
+
sin
⁡
[
x
+
2
y
]
)
0000000
+
3
sin
⁡
[
x
+
3
y
]
)
,
u
3
(
x
,
y
,
t
)
=
1
96
Gamma
[
1
+
α
]
3
t
3
α
×
6
(
15
+
38
x
+
6
x
2
-
12
cos
⁡
⁡
[
x
]
-
3
cos
⁡
⁡
[
2
x
]
+
4
sin
⁡
[
x
]
x
2
)
+
48
cos
⁡
⁡
[
y
]
(
(
-
4
+
x
)
(
-
1

+

cos
⁡
⁡
[
x
]
)

-

(
9
+
x
)
sin
⁡
[
x
]
)
+
4
cos
⁡
⁡
[
2
y
]
(
-
189
-
135
x
+
48
x
2
000000000000
0000
+
16
cos
⁡
⁡
[
x
]
+
(
173
-
9
x
)
cos
⁡
⁡
[
2
x
]
000000000000
0000
+
8
sin
⁡
[
x
]
+
(
-
190
+
9
x
)
sin
⁡
[
2
x
]
x
2
)
-
24
cos
⁡
⁡
[
3
y
]
(
-
8
+
6
x
+
(
3
+
2
x
)
cos
⁡
⁡
[
x
]

+

5
cos
⁡
⁡
[
3
x
]
0000000000000
0000
+
(
5
+
2
x
)
sin
⁡
[
x
]
+
35
sin
⁡
[
3
x
]
)
+
cos
⁡
⁡
[
4
y
]
(
21
+
8
cos
⁡
⁡
[
x
]
-
29
cos
⁡
⁡
[
4
x
]
+
8
sin
⁡
[
x
]
00000000000
000
+
18
sin
⁡
[
2
x
]
-
14
sin
⁡
[
4
x
]
)
+
48
(
8

+

3
x

+

(
-
8

+

x
)
cos
⁡
⁡
[
x
]

+

(
11

+

x
)
sin
⁡
[
x
]
)
sin
⁡
[
y
]
+
4
(
186
-
153
x
-
48
x
2
+
4
cos
⁡
⁡
[
x
]
000000
00
+
(
-
190
+
9
x
)
cos
⁡
⁡
[
2
x
]
-
20
sin
⁡
[
x
]
0000
0000
+
(
-
173
+
9
x
)
sin
⁡
[
2
x
]
x
2
)
sin
⁡
[
2
y
]
-
24
(
-
40
-
6
x
+
(
5
+
2
x
)
cos
⁡
⁡
[
x
]
+
35
cos
⁡
⁡
[
3
x
]
0000000
00
-
(
3
+
2
x
)
sin
⁡
[
x
]
-
5
sin
⁡
[
3
x
]
)
sin
⁡
[
3
y
]
+
(
8
cos
⁡
⁡
[
x
]

+

18
cos
⁡
⁡
[
2
x
]

-

2
(
6

+

7
cos
⁡
⁡
[
4
x
]

+

4
sin
⁡
[
x
]
)
00000
00
+
29
sin
⁡
[
4
x
]
)
sin
⁡
[
4
y
]
,
u
4
(
x
,
y
,
t
)
=
1
2880
Gamma
[
1
+
α
]
4
×
t
4
α
(
-
120
cos
⁡
[
2
y
]
×
(
1173
+
2
x
(
-
351
+
x
(
-
123
+
32
x
)
)
-
96
cos
⁡
[
x
]
+
3
(
-
359
+
24
x
)
cos
⁡
[
2
x
]
+
24
(
2
+
x
+
6
x
cos
⁡
[
x
]
)
sin
[
x
]
-
489
sin
[
2
x
]
)
+
360
(
-
39
-
36
x
+
(
26
+
6
x
)
cos
⁡
[
x
]
0000000000000000
+
(
13
+
6
x
)
cos
⁡
[
2
x
]
0000000000000000
-
2
(
1
+
x
)
sin
[
x
]
-
5
sin
[
2
x
]
)
-
30
cos
⁡
[
y
]
(
x
2
892
-
24
(
-
5
+
x
)
x
000000000000000000000
+
12
(
-
83
+
66
x
+
8
x
2
)
cos
⁡
[
x
]
00000000000000000000
0
+
84
cos
⁡
[
2
x
]
+
20
cos
⁡
[
3
x
]
00000000000000000000
0
-
3
(
-
405

+

4
x
(
62

+

13
x
)
)
sin
[
x
]
00000000000000000000
0
+
36
sin
[
2
x
]
+
15
sin
[
3
x
]
x
2
)
+
5
cos
⁡
[
3
y
]
(
8
(
-
15883
+
18
x
(
-
71
+
27
x
)
)
000000000000000
00000
+
36
(
45
+
x
(
53
+
16
x
)
)
cos
⁡
[
x
]
00000000000000
000000
+
216
cos
⁡
[
2
x
]
0000000000000
0000000
+
4
(
31307
-
105
x
)
cos
⁡
[
3
x
]
000000000000
00000000
+
9
(
-
37
+
4
x
(
85
+
16
x
)
)
sin
[
x
]
00000000000
000000000
+
432
sin
[
2
x
]
0000000000
0000000000
+
5
(
-
3665
+
12
x
)
sin
[
3
x
]
)
-
180
cos
⁡
[
4
y
]
(
111

+
56
x

+

12
(
7
+
x
)
cos
⁡
[
x
]
00000000000000000000000
+
6
cos
⁡
[
2
x
]
-
201
cos
⁡
[
4
x
]
00000000000000000000000
+
4
(
7
+
3
x
)
sin
[
x
]
00000000000000000000000
+
2
(
13
+
6
x
)
sin
[
2
x
]
0000000000000000000000
0
+
417
sin
[
4
x
]
)
+
3
cos
⁡
[
5
y
]
(
288
+
90
cos
⁡
[
x
]
-
150
cos
⁡
[
3
x
]
000000000000000000000
-
228
cos
⁡
[
5
x
]
+
105
sin
[
x
]
000000000000000000000
+
120
sin
[
2
x
]
+
200
sin
[
3
x
]
000000000000000000000
-
369
sin
[
5
x
]
)
-
30
(
X
2
(
609
+
12
x
(
-
41
+
3
x
)
)
cos
⁡
[
x
]
+
36
cos
⁡
[
2
x
]
+
15
cos
⁡
[
3
x
]
+
12
(
-
55
+
50
x
+
6
x
2
000000000000000000
+
(
105
+
x
(
51
+
8
x
)
)
sin
[
x
]
000000000000000000
-
5
sin
[
2
x
]
X
2
)
-
20
sin
[
3
x
]
X
2
)
sin
[
y
]
+
120
(
-
453
-
630
x
+
330
x
2
+
64
x
3
00000000000000
+
12
(
-
3
+
x
)
cos
⁡
[
x
]
00000000000000
+
(
489
-
72
x
)
cos
⁡
[
2
x
]
00000000000000
+
12
(
19
+
3
x
)
sin
[
x
]
00000000000000
+
3
(
-
359
+
24
x
)
sin
[
2
x
]
x
3
)
sin
[
2
y
]
+
5
(
9
(
-
51
+
4
x
(
85
+
16
x
)
)
cos
⁡
[
x
]
+
432
cos
⁡
[
2
x
]
+
5
(
-
3665
+
12
x
)
cos
⁡
[
3
x
]
+
4
(
4588
-
36
x
(
361
+
27
x
)
-
9
(
48
+
x
(
53
+
16
x
)
)
sin
[
x
]
-
54
sin
[
2
x
]
+
(
-
31307
+
105
x
)
sin
[
3
x
]
)
)
sin
[
3
y
]
-
180
(
4
(
7
+
3
x
)
cos
⁡
[
x
]
+
2
(
13
+
6
x
)
cos
⁡
[
2
x
]
00000000000000
-
3
(
157
+
12
x
-
139
cos
⁡
[
4
x
]
000000000000000000
+
4
(
7
+
x
)
sin
[
x
]
+
2
sin
[
2
x
]
000000000000000000
-
67
sin
[
4
x
]
)
)
sin
[
4
y
]
+
3
(
-
56
+
105
cos
⁡
[
x
]
+
120
cos
⁡
[
2
x
]
+
200
cos
⁡
[
3
x
]
-
369
cos
⁡
[
5
x
]
-
90
sin
[
x
]
0
+
150
sin
[
3
x
]
+
228
sin
[
5
x
]
)
sin
[
5
y
]
)
.
By using the proposed algorithm, one can obtain the remaining terms but for simplicity we stopped at 5 terms. The numerical solution is depicted in Figures 1 and 2.Figure 1
Solution for alpha equals 0.5 andy
=
0.Figure 2
Solution for alpha equals 0.05 andy
=
Pi.


## 6. Conclusion

In the recent years several physical phenomena were explained with great success in the light of the concept of noninteger order derivatives. More importantly to truthfully photocopy the nonlocal, frequency- and history-dependent possessions of power law phenomena, selected dissimilar modelling apparatuses constructed upon fractional operators have to be initiated. Especially, the remunerations of fractional calculus and fractional order models and their applications in the field of nonlinear wave motion have previously been intensively reconsidered during the last few eras with outstanding outcome. With this in mind, we have investigated and considered the motion of nonlinear wave within the folder of fractional derivative. A careful investigation of the stability, convergence, and uniqueness analysis has been done.



---

*Source: 481657-2014-05-22.xml*