# Extension of Matched Asymptotic Method to Fractional Boundary Layers Problems

Authors: Abdon Atangana; Emile Franc Doungmo Goufo
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

We were concerned with the description of the boundary layers problems within the scope of fractional calculus. However, we will note that one of the main methods used to solve these problems is the matched asymptotic method. We should mention that this was not achievable via the existing fractional derivative definitions, because they do not obey the chain rule.  In order to accommodate the matched asymptotic method to the scope of fractional derivative, we proposed a relatively new derivative called thebeta-derivative. We presented some useful information for this operator. With the reward of this operator, we presented the idea of matched asymptotic method in finding solutions of the fractional boundary layers problems. The method was illustrated with an example.


---

## Body

## 1. Introduction

In a large class of singular perturbed problems, the domain may be divided into two or more subdomains. In one of these, often the largest, the solution is accurately approximated by an asymptotic series found by treating as a regular perturbation. The other subdomains consist of one or more small areas in which that approximation is inaccurate, generally because the perturbation terms in the problem are not negligible there. These areas are referred to as transition layers and as boundary or interior layers depending on whether they occur at the domain boundary or inside the domain [1–6]. An approximation in the form of an asymptotic series is obtained in transition layers by treating that part of the domain as a separate perturbation problem. This approximation is called the “inner solution,” and the other is the “outer solution” named for their relationship to transition layers. The outer and inner solutions are then combined through a process called “matching” in such a way that an approximate solution for the whole domain is obtained [1–6]. Therefore, the matched asymptotic method is a common approach to finding an accurate approximation to solution to an equation in particular when solving perturbed differential equation with conventional order derivative. This class of differential equations is used to describe real world problems; for instance, in physic and fluid mechanics, a boundary layer is the layer of fluid in the immediate neighbourhood of a bouncing surface where the consequences of viscidness are noteworthy. In the earth’s atmosphere, the environmental borderline layer is the air layer neighbouring the earth affected by diurnal heat, moisture, or momentum transfer to or from the surface. On an aircraft wing, the boundary layer is the part of the flow close to the wing, where viscous forces distort the surrounding nonviscous flow.In the recent decades, attention has been paid by several scholars to modeling the real world problems with the concept of fractional order derivatives. It was revealed by many proofs and results that the modelling of these real world problems with the concept of fractional order derivatives gives better prediction rather than using conventional derivative, which is regarded as integer order derivatives.It is perhaps important to mention that the matched asymptotic method has never been used to solve any kind of fractional differential equations because of the nature and properties of the fractional derivative. In particular, the most commonly used fractional derivative (Caputo fractional derivative) for modelling the real world problem does not obey the chain rule which is one of the key elements of the matched asymptotic method. Recently, the so-called conformable fractional derivative was proposed. This fractional derivative is theoretically very easier to handle and also obeys some conventional properties that cannot be satisfied by the existing fractional derivatives, for instance, the chain rule. However this fractional derivative has a very big weakness, which is the fractional derivative of any differentiable function at the point zero and this does not satisfy any physical problem or cannot for the moment have any physical interpretation. A modified version was proposed in order to extend the limitation of the conformable derivative; however this derivative depends on the interval on which the function is being differentiated which is also a true problem for some physical problem.The aim of this work is to further extend the so-called matched asymptotic method in the scope of the fractional differential equation on the one hand. On the other hand, we will propose a suitable derivative using fractional order, that is easier and respects the chain rule. This version of fractional derivative will be used to further describe the boundary layers problems within the folder of fractional calculus. We will start with some existing definitions of fractional derivatives in the coming section.


## 2. Some Definitions of Fractional Derivatives

There exist a quite few definitions of fractional derivatives in the literature; we will present the definition of a few [7–12].(1)
The Riemann-Liouville fractional derivative of a function saysf is defined as
(1)
D
x
α
f
x
=
1
Γ
n
-
α
d
d
x
n
∫
0
x
x
-
t
n
-
α
-
1
f
t
d
t
,
kkkkkkkkkkkkkkkkkkkk
kk
n
-
1
<
α
≤
n
.
(2)
The Caputo fractional derivative of a differentiable function saysf is defined as
(2)
D
x
α
f
(
x
)
=
1
Γ
n
-
α
∫
0
x
x
-
t
n
-
α
-
1
d
d
t
n
f
t
d
t
,
kkkkkkkkkkkkkkkkkkkkk
k
n
-
1
<
α
≤
n
.
(3)
The modified Liouville fractional derivative of a function saysf is defined as
(3)
D
x
α
f
(
x
)
=
1
Γ
n
-
α
d
d
x
n
×
∫
0
x
x
-
t
n
-
α
-
1
f
t
-
f
(
0
)
d
t
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
n
-
1
<
α
≤
n
.
(4)
The local fractional derivative of a function saysf is defined as
(4)
L
x
α
f
(
x
)
=
lim
⁡
x
→
x
0
⁡
Γ
1
+
α
f
x
-
f
x
0
x
-
x
0
α
.
(5)
The conformable fractional derivative of a function saysf is given as
(5)
T
α
f
(
x
)
=
lim
⁡
ε
→
0
⁡
f
x
+
ε
x
1
-
α
-
f
(
x
)
ε
.
(6)
The modified conformable fractional derivative of a given functionf defined in the interval (
a
,
b
) is given as
(6)
T
α
f
(
x
)
=
lim
⁡
ε
→
0
⁡
f
x
+
ε
x
-
a
1
-
α
-
f
(
x
)
ε
.The last two seem to satisfy some common properties of the standard concept of derivative, but as we said before they have some limitations that will not allow them to be used in modelling real world problems. We will therefore propose a suitable fractional derivative that will allow us to escape the lack of these fractional derivatives.Definition 1.
Letf be a function, such that f
:
[
a
,
∞
)
→
R. Then, the beta derivative of a function f is defined as
(7)
D
0
A
x
β
f
x
=
lim
⁡
ε
→
0
⁡
f
x
+
ε
x
+
1
/
Γ
β
1
-
β
-
f
(
x
)
ε
,
for all x
≥
a, β
∈
(
0,1
]. Then if the limit of the above exists, f is said to be β-differentiable.
One can remark that the above definition does not depend on the interval on which the function is defined as in the case of Definition6. If the function is differentiable, our definition at a point zero is different to zero, which was not the case for the 5 definitions.Theorem 2.
Assuming that a given function saysf
:
[
a
,
∞
)
→
R is β-differentiable at the points says x
0
≥
a, β
∈
(
0
,
1
], then f is also continuous at x
0.Proof.
Assuming thatf is β-differentiable, then
(8)
D
0
A
x
β
f
x
0
=
lim
⁡
ε
→
0
⁡
f
x
0
+
ε
x
0
+
1
/
Γ
β
1
-
β
-
f
(
x
0
)
ε
exists. Therefore,
(9)
lim
⁡
ε
→
0
⁡
f
x
0
+
ε
x
0
+
1
Γ
β
1
-
β
-
f
x
0
=
lim
⁡
ε
→
0
⁡
f
x
0
+
ε
x
0
+
1
/
Γ
β
1
-
β
-
f
(
x
0
)
ε
ε
=
D
A
x
β
f
x
0
·
0
=
0
,
thus giving
(10)
lim
⁡
ε
→
0
⁡
f
x
0
+
ε
x
0
+
1
Γ
β
1
-
β
=
lim
⁡
ε
→
0
⁡
f
x
0
+
ε
x
0
+
1
Γ
β
1
-
β
-
f
(
x
0
)
+
f
x
0
=
0
+
f
x
0
.
However, taking x
=
x
0
+
ε
x
0
+
1
/
Γ
β
1
-
β so that ε
=
x
-
x
0
/
x
0
+
1
/
Γ
β
1
-
β, thus
(11)
lim
⁡
(
x
-
x
0
)
/
(
x
0
+
1
/
Γ
(
β
)
)
1
-
β
→
0
⁡
f
x
=
f
x
0
;
since x
0
+
1
/
Γ
β
1
-
β
≠
0, then the above can be rewritten as
(12)
lim
⁡
x
-
x
0
→
0
⁡
f
x
=
f
x
0
or
lim
⁡
x
→
x
0
⁡
f
x
=
f
x
0
,
which completes the proof.Theorem 3.
Assuming thatf is β-differentiable on an open interval (
a
,
b
), then(1)
ifD
0
A
x
β
f
(
x
)
<
0 for all x
∈
a
,
b, then f is decreasing there;
(2)
ifD
0
A
x
β
f
(
x
)
>
0 for all x
∈
a
,
b, then f is increasing there;
(3)
ifD
0
A
x
β
f
(
x
)
=
0 for all x
∈
a
,
b, then f is constant there.Proof.
Assume thatf is β-differentiable on an open interval (
a
,
b
) and D
0
A
x
β
f
(
x
)
>
0; then, consider x
1
>
x
2
(13)
f
x
1
+
ε
x
2
+
1
/
Γ
β
1
-
β
-
f
x
2
ε
>
0
.
Then,
(14)
f
x
1
+
ε
x
2
+
1
Γ
β
1
-
β
-
f
x
2
>
0
.
Now, taking the limit on both sides, we have that
(15)
lim
⁡
ε
→
0
⁡
f
x
1
+
ε
x
2
+
1
Γ
β
1
-
β
-
f
x
2
=
f
x
1
-
f
x
2
>
0
.
This completes the proof. Note that, by using a similar approach, we can also prove 2 and 3.Theorem 4.
Assuming thatg
≠
0 and f are two functions β-differentiable with β
∈
(
0,1
], then the following relations can be satisfied:(1)
D
0
A
x
α
a
f
x
+
b
g
(
x
)
=
a
D
0
A
x
α
f
x
+
b
D
0
A
x
α
f
x for all a and b real number;
(2)
D
0
A
x
α
c
=
0 for c any given constant;
(3)
D
0
A
x
α
f
x
g
(
x
)
=
g
x
D
0
A
x
α
f
x
+
f
x
D
0
A
x
α
g
x;
(4)
D
0
A
x
α
f
x
/
g
(
x
)
=
g
x
D
0
A
x
α
f
x
-
f
x
D
0
A
x
α
g
x
/
g
2
(
x
).The proofs of the above relations are the same as the one in [12].Theorem 5.
Letf
:
[
a
,
∞
)
→
R be a function such that f is differentiable and also α-differentiable. Let g be a function defined in the range of f and also differentiable; then, one has the following rule:
(16)
D
0
A
x
β
g
o
f
(
x
)
=
x
+
1
Γ
β
1
-
β
f
′
x
g
′
f
x
.Proof.
We have from the definition that(17)
lim
⁡
ε
→
0
⁡
g
o
f
x
+
ε
x
+
1
/
Γ
β
1
-
β
-
g
o
f
x
ε
.
Let h
=
ε
x
+
1
/
Γ
β
1
-
β such that
(18)
lim
⁡
ε
→
0
⁡
g
o
f
x
+
ε
x
+
1
/
Γ
β
1
-
β
-
g
o
f
x
ε
=
x
+
1
Γ
β
1
-
β
lim
⁡
h
→
0
⁡
g
o
f
x
+
h
-
g
o
f
x
h
.
According to chain rule, we have that
(19)
lim
⁡
h
→
0
⁡
g
o
f
x
+
h
-
g
o
f
x
h
=
f
′
x
g
′
f
x
.
Therefore, we have
(20)
D
A
x
β
g
o
f
(
x
)
=
x
+
1
Γ
β
1
-
β
f
′
x
g
′
f
(
x
)
.
The above rule is referred to as Atangana beta-rule. We will note that if β
=
1, we recover the Chain rule. We will present new derivative for some special functions.(1)
The beta-derivative of the Mittag-Leffler is given as(21)
D
0
A
x
β
E
α
(
x
)
=
x
+
1
Γ
β
1
-
β
∑
n
=
1
∞
n
x
n
-
1
Γ
α
n
+
1
.
(2)
The beta-derivative of Bessel function of first kind is given as(22)
D
0
A
x
β
J
v
(
x
)
=
1
2
x
+
1
Γ
β
1
-
β
×
∑
n
=
1
∞
-
1
n
n
-
1
!
2
n
+
v
x
/
2
2
n
+
v
-
1
n
!
Γ
v
+
n
+
1
.
(3)
The beta-derivative of Error function is given as(23)
D
0
A
x
β
E
r
f
(
x
)
=
x
+
1
Γ
β
1
-
β
2
π
e
-
x
2
.We can carry on and find all the fractional derivatives of the special functions; however, we will notice that if β
=
1, we obtain without any problem their derivatives with the standard derivative. We will now look at the inverse operator of the proposed fractional derivative.Definition 6.
Letf
:
[
a
,
∞
)
→
R be given function; then one proposes that Atangana’s beta-integral of f is
(24)
I
a
A
x
β
f
(
x
)
=
∫
a
x
t
+
1
Γ
β
β
-
1
f
t
d
t
.
The above operator is the inverse operator of the proposed fractional derivative. One will present to underpin this statement by the following theorem.Theorem 7.
D
0
A
x
α
I
0
A
x
α
f
(
x
)
=
f
x for all x
≥
a with f being a given continuous and differentiable function.Proof.
Letf be a continuous function; then, by definition, if we let I
0
A
x
α
f
x
=
F
(
x
), we have
(25)
D
0
A
x
α
I
0
A
x
α
f
(
x
)
=
lim
⁡
ε
→
0
⁡
F
x
+
ε
x
+
1
/
Γ
β
1
-
β
-
F
(
x
)
ε
=
x
+
1
Γ
β
1
-
β
d
F
(
x
)
d
x
=
x
+
1
Γ
β
1
-
β
d
d
x
×
∫
a
x
t
+
1
Γ
β
β
-
1
f
(
t
)
d
t
=
x
+
1
Γ
β
1
-
β
x
+
1
Γ
β
β
-
1
f
x
=
f
x
.
This completes the proof.Theorem 8.
I
a
A
x
α
D
0
A
x
α
f
(
x
)
=
f
x
-
f
(
a
) for all x
≥
a with f being a given continuous and differentiable function.Proof.
Sincef is continuous and differentiable, by definition, we have that
(26)
I
a
A
x
α
D
0
A
x
α
f
(
x
)
=
∫
a
x
t
+
1
Γ
β
β
-
1
D
A
t
α
f
(
t
)
d
t
=
∫
a
x
t
+
1
Γ
β
β
-
1
t
+
1
Γ
β
1
-
β
kkkk
×
d
f
x
d
x
d
t
=
∫
a
x
d
f
(
x
)
d
x
d
t
=
f
x
-
f
(
a
)
.
This completes the proof.We will present beta-integral of some simple functions.(1)
Forf
x
=
c as given constant, we have
(27)
I
0
A
x
α
c
=
∫
0
x
t
+
1
Γ
β
β
-
1
c
d
t
=
c
Γ
α
-
α
-
1
+
1
+
x
Γ
α
α
α
.
(2)
Forf
x
=
c
x
n for any given natural number n and real number a, we have
(28)
∫
0
x
t
+
1
Γ
β
β
-
1
c
t
n
d
t
=
-
c
-
x
-
n
x
n
Beta
-
x
Γ
α
,
1
+
n
,
α
Γ
α
-
n
-
α
.
(3)
Forf
x
=
cos
⁡
⁡
(
x
), we have the following:
(29)
∫
0
x
t
+
1
Γ
β
β
-
1
cos
⁡
⁡
t
d
t
=
-
1
2
ⅇ
-
1
/
2
ⅈ
π
α
×
-
ⅇ
ⅈ
π
α
Γ
α
,
-
ⅈ
Γ
α
k
k
k
k
k
k
k
×
cos
⁡
1
Γ
α
-
ⅈ
sin
1
Γ
α
k
k
k
k
k
k
k
-
Γ
α
,
ⅈ
Γ
α
cos
⁡
1
Γ
α
+
ⅈ
sin
1
Γ
α
-
1
2
1
+
x
Γ
α
α
1
+
x
Γ
α
2
-
α
×
ⅇ
α
Conjugate
Log
ⅈ
1
+
x
Γ
α
-
Log
Γ
α
Γ
α
ⅈ
1
+
x
Γ
α
Γ
α
α
Γ
α
α
Γ
α
,
-
ⅈ
1
+
x
Γ
α
Γ
α
kkkkkkk
×
ⅇ
α
Conjugate
Log
ⅈ
1
+
x
Γ
α
-
Log
Γ
α
Γ
α
ⅈ
1
+
x
Γ
α
Γ
α
α
kkkkkkk
×
cos
⁡
1
Γ
α
-
ⅈ
sin
1
Γ
α
kkkkkkk
+
-
ⅈ
1
+
x
Γ
α
α
Γ
kkkkkkk
×
α
,
ⅈ
1
+
x
Γ
α
Γ
α
kkkkkkk
×
cos
⁡
1
Γ
α
+
ⅈ
sin
1
Γ
α
ⅇ
α
Conjugate
Log
ⅈ
1
+
x
Γ
α
-
Log
Γ
α
Γ
α
ⅈ
1
+
x
Γ
α
Γ
α
α
.
(4)
Forf
x
=
exp
⁡
⁡
[
t
], we have the following:
(30)
∫
0
x
t
+
1
Γ
β
β
-
1
cos
⁡
⁡
[
t
]
d
t
=
ⅇ
-
1
/
Γ
[
α
]
1
Γ
α
-
1
-
x
Γ
α
-
α
1
+
x
Γ
α
α
k
k
k
k
k
k
k
k
k
×
Γ
α
-
Γ
α
,
-
x
-
1
Γ
α
k
k
k
k
k
k
k
k
k
+
-
1
-
α
-
Γ
α
+
Γ
α
,
-
1
Γ
α
.Other properties of the proposed operator will be investigated in the next work. However, we will now present the application of this operator in the scope of fractional boundary layers problems.


## 3. A Matched Asymptotic Method for Fractional Boundary Layer Problems

In this section, we will present the discussion that underpins the extension of the so-called matched asymptotic method in the scope of fractional boundary layer problem. The fractional derivative used here will be the beta-derivative because it provides us with useful rule called the new chain rule that cannot be obtained by the existing fractional derivative operator. We will then present the general methodology. Many singular perturbed differential equations have solutions that change rapidly in a narrow region. This may occur in an initial layer where there is a rapid adjustment of initial differential conditions to a quasisteady state, in a boundary layer where the solution away from the boundary adjusts to a boundary condition, or in an interior layer such as a propagation wave front. This class of problems will be further investigated in the scope of fractional calculus by using the fractional matched asymptotic method.Consider the following fractional singular perturbation problem:(31)
P
α
u
(
x
,
ε
)
=
0
,
u
a
=
A
,
u
b
=
B
,
where α
≥
1 is the fractional order derivative. To solve the above equation, we assume that
(32)
u
x
,
ε
=
∑
n
=
0
∞
ε
n
u
n
x
.
And then we can find the outer solution using the same routine as in the case of ordinary differential equation with integer order derivatives.To find the inner solution, we also follow the same routine as in the case of differential equation with integer order derivative, with the difference that after changing the variable we make use of Atangana’s beta-rule provided in Theorem5. We will illustrate this with an example.Example 9.
Consider the following fractional singular perturbation problem:(33)
ε
y
′′
+
1
+
ε
D
0
A
x
α
y
x
+
y
=
0
,
0
<
ε
≪
1
,
subject to the following initial and boundary conditions:
(34)
y
0
=
0
,
y
1
=
1
.
Outer solution, valid for t
=
O
(
1
), implying equation (33) is reduced to
(35)
D
0
A
x
α
y
x
+
y
=
0
.
But the solution is two times differentiable; this will imply that the fraction version will be rewritten as
(36)
x
+
1
Γ
β
1
-
α
y
′
+
y
=
0
,
by which the solution is given as
(37)
y
x
,
β
=
c
Exp
-
x
+
1
/
Γ
α
2
-
α
2
-
α
.
Applying the condition y
0
=
0, we have c
=
0; applying the boundary condition y
1
=
1, we have c
=
Exp
1
/
Γ
(
β
)
2
-
α
/
2
-
α. This is not possible; therefore, according to the matched asymptotic method, we consider the outer solution c
=
Exp
1
/
Γ
(
β
)
2
-
α
/
2
-
α. Therefore, the following solution is the leading-order solution in the outer region.
The inner solution is valid fort
=
O
(
ε
). In the inner region, t and ε are both tiny, but of comparable size, so define the new O
(
1
) time variable ξ
=
x
/
ε. Rescaling the original fractional boundary value problem by replacing t with ξ
ε, the problem becomes
(38)
1
ε
y
′′
(
ξ
)
+
1
+
ε
D
A
x
α
y
ξ
ε
+
y
(
ξ
)
=
0
.
The fractional part can be converted using first the fact that y is differentiable and secondly using theorem to obtained
(39)
1
ε
y
′′
(
ξ
)
+
ξ
+
1
Γ
α
1
-
α
1
ε
1
+
ε
y
′
ξ
+
y
(
ξ
)
=
0
.
Of course, by multiplying by the small parameter and taking it to be zero, we will obtain
(40)
y
′′
ξ
+
ξ
+
1
Γ
α
1
-
α
y
′
ξ
=
0
.
The solution is given as
(41)
y
ξ
,
β
,
ε
=
C
∫
Exp
-
ξ
+
1
/
Γ
α
2
-
α
2
-
α
d
ξ
=
C
1
Γ
α
+
ξ
kkkkk
×
ExpIntegral
E
1
+
1
-
2
+
α
,
ξ
+
1
/
Γ
α
2
-
α
2
-
α
kkkkk
×
-
2
+
α
-
1
1
Γ
α
+
B
,
where ExpIntegral
E is the exponential integral function defined as
(42)
ExpIntegral
E
n
,
z
=
∫
x
∞
exp
⁡
⁡
[
-
x
t
]
t
n
d
t
.
Therefore, the general solution for inner part is given as
(43)
y
ξ
,
β
,
ε
=
B
+
C
×
1
Γ
α
+
ξ
kkkkkkkk
×
ExpIntegral
E
1
+
1
-
2
+
α
,
ξ
+
1
/
Γ
α
2
-
α
2
-
α
kkkkkkkk
×
-
2
+
α
-
1
1
Γ
α
.
To find the constants B and C, we apply the boundary conditions. First, for t
=
0, which corresponds to the inner region, we have
(44)
B
+
C
1
Γ
α
k
k
k
k
k
k
×
ExpIntegral
E
1
+
1
-
2
+
α
,
1
/
Γ
α
2
-
α
2
-
α
k
k
k
k
k
k
×
-
2
+
α
-
1
1
Γ
α
=
0
,
such that the outer solution becomes
(45)
y
ξ
,
β
,
ε
=
-
C
1
Γ
α
k
k
k
k
k
k
k
×
ExpIntegral
E
1
+
1
-
2
+
α
,
1
/
Γ
α
2
-
α
2
-
α
k
k
k
k
k
k
k
×
-
2
+
α
-
1
1
Γ
α
+
C
1
Γ
α
+
ξ
kkkkkkkk
×
ExpIntegral
E
1
+
1
-
2
+
α
,
ξ
+
1
/
Γ
α
2
-
α
2
-
α
kkkkkkkk
×
-
2
+
α
-
1
1
Γ
α
.
Currently working to find the constant C, we employ the harmonising method also called matching method. The principal idea here is that the inner and the outer solutions should harmonise for values of x in an intermediary region, where ε
≪
t
≪
1. We demand the outer limit of the inner solution to match the limit of the outer solution, which in mathematical formula implies
(46)
lim
⁡
ξ
→
∞
⁡
y
I
=
lim
⁡
x
→
0
⁡
y
O
which corresponds to
(47)
C
=
-
1
×
1
Γ
α
kkkkk
×
ExpIntegral
E
1
+
1
-
2
+
α
,
1
/
Γ
α
2
-
α
2
-
α
kkkkk
×
-
2
+
α
-
1
1
Γ
α
-
1
.Composite Solution. To have the final matched or composite solution, which is valid, on the whole domain, one of the commonly used methods is the so-called uniform method. This method consists of an addition of the inner and the outer approximations and subtracting their overlapping value y
overlap, which would otherwise be counted twice. The overlapping value is the limit of the inner boundary layer solution and the inner limit of the outer solution. In mathematical formula, we have the following:
(48)
y
x
=
y
I
x
+
y
O
x
-
y
overlap
.


## 4. Conclusion

Boundary layers problems are a class of problems that can be found in many areas of real world situations. In particular, numerous extraordinary disconcerted differential equations have explanations that transform precipitously in a contracted section. This may happen in a preliminary layer everywhere there is a speedy modification of original discrepancy surroundings to a quasisteady state, in a boundary layer where the solution away from the borderline adjusts to a boundary condition, or in an interior layer such as a propagation wave front. The method in charge of solving these equations is the so-called matched asymptotic method. Due to the nature of this method, it was not possible to use it to solve the fractional boundary layer problems, because the existing fractional order derivative does not obey the so-called chain rule. In order to solve this problem, we have presented a relatively novel derivative, which is more suitable than the existing one. The properties of this derivative allowed us to extend the matched asymptotic method to the scope of fractional calculus. We have presented in detail an example to exemplify the method.



---

*Source: 107535-2014-11-20.xml*