# On Generalized Fractional Kinetic Equations Involving Generalized Bessel Function of the First Kind

Authors: Dinesh Kumar; S. D. Purohit; A. Secer; A. Atangana
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

We develop a new and further generalized form of the fractional kinetic equation involving generalized Bessel function of the first kind. The manifold generality of the generalized Bessel function of the first kind is discussed in terms of the solution of the fractional kinetic equation in the paper. The results obtained here are quite general in nature and capable of yielding a very large number of known and (presumably) new results.


---

## Body

## 1. Introduction and Preliminaries

Bessel functions are playing the important role in studying solutions of differential equations, and they are associated with a wide range of problems in important areas of mathematical physics, like problems of acoustics, radiophysics, hydrodynamics, and atomic and nuclear physics. These considerations have led various workers in the field of special functions to explore the possible extensions and applications of the Bessel functions. Among many properties of Bessel functions, they also have investigated some possible extensions of the Bessel functions.The generalized Bessel function of the first kindω
p
z is defined for z
∈
C
∖
{
0
} and b, c, p
∈
C
R
p
>
-
1 by the following series [1, page 10, (
1.15
)] (for recent work, see also [2–6]):
(1)
ω
p
,
b
,
c
z
=
ω
p
z
=
∑
k
=
0
∞
‍
-
1
k
c
k
k
!
Γ
p
+
b
+
1
/
2
+
k
z
2
2
k
+
p
,
where C denotes the set of complex numbers and Γ
a is the familiar Gamma function.The special cases of series (1) can be obtained as follows.(i)
If we putb
=
c
=
1 in (1), then we obtain the familiar Bessel function of the first kind [7] of order p for z
,
p
∈
C with R
p
>
-
1 defined and represented by the following expressions (see also [1, 8]):
(2)
J
p
z
=
∑
k
=
0
∞
-
1
k
k
!
Γ
p
+
k
+
1
z
2
2
k
+
p
,
z
∈
C
.
(ii)
Puttingb
=
1 and c
=
-
1 in series (1), we get the modified Bessel function of the first kind of order p defined by (see [1, 7])
(3)
I
p
z
=
∑
k
=
0
∞
1
k
!
Γ
p
+
k
+
1
z
2
2
k
+
p
,
z
∈
C
;
the series given by (3) is also a special case of Galué’s generalized modified Bessel function [9] depending on parameters a
=
0,1
,
2
,
… and p
>
-
1, given as follows:
(4)
I
a
p
z
=
∑
k
=
0
∞
z
/
2
2
k
+
p
k
!
Γ
p
+
a
k
+
1
,
z
∈
C
.
(iii)
Lettingb
=
2 and c
=
1 in series (1), we have the spherical Bessel function of the first kind of order p defined by (see [1])
(5)
j
p
z
=
π
2
∑
k
=
0
∞
-
1
k
k
!
Γ
p
+
k
+
3
/
2
z
2
2
k
+
p
,
z
∈
C
.Furthermore, Deniz et al. [10] considered the function φ
p
,
b
,
c
z, defined in terms of the generalized Bessel function ω
p
z, by the transformation
(6)
φ
p
,
b
,
c
z
=
2
p
Γ
p
+
b
+
1
2
z
1
-
p
/
2
ω
p
z
=
z
+
∑
k
=
1
∞
-
c
k
4
k
ν
k
z
k
+
1
k
!
,
where ν
=
p
+
(
b
+
1
)
/
2
∉
Z
0
-
≔
{
0
,
-
1
,
-
2
,
…
} and (
a
)
k is the Pochhammer symbol defined (for a
∈
C) by
(7)
a
k
≔
1
k
=
0
a
a
+
1
⋯
a
+
k
-
1
k
∈
N
≔
1,2
,
3
,
…
=
Γ
a
+
k
Γ
a
a
∈
C
∖
Z
0
-
.Fractional differential equations appear more and more frequently for modeling of relevant systems in several fields of applied sciences. These equations play important roles, not only in mathematics, but also in physics, dynamical systems, control systems, and engineering, to create the mathematical model of many physical phenomena. In particular, the kinetic equations describe the continuity of motion of substance and are the basic equations of mathematical physics and natural science. Therefore, in literature we found several papers that analyze extensions and generalizations of these equations involving various fractional calculus operators. One may, for instance, refer to such type of works by [11–23].Haubold and Mathai [13] have established a functional differential equation between rate of change of reaction, the destruction rate, and the production rate as follows:
(8)
d
N
d
t
=
-
d
N
t
+
p
N
t
,
where N
=
N
t is the rate of reaction, d
=
d
N is the rate of destruction, p
=
p
N is the rate of production, and N
t denotes the function defined by N
t
t
*
=
N
t
-
t
*, t
*
>
0.Haubold and Mathai studied a special case of (8), when spatial fluctuations or inhomogeneities in the quantity N
t are neglected, is given by the equation
(9)
d
N
i
d
t
=
-
c
i
N
i
t
,
together with the initial condition that N
i
t
=
0
=
N
0, is the number of density of species i at time t
=
0, c
i
>
0. If we decline the index i and integrate the standard kinetic equation (9), we have
(10)
N
t
-
N
0
=
-
c
D
0
t
-
1
N
t
,
where D
0
t
-
1 is the special case of the Riemann-Liouville integral operator D
0
t
-
ν defined as
(11)
D
0
t
-
ν
f
(
t
)
=
1
Γ
ν
∫
0
t
t
-
s
ν
-
1
f
s
d
s
,
t
>
0
,
R
(
ν
)
>
0
.
Haubold and Mathai [13] have given the fractional generalization of the standard kinetic equation (10) as
(12)
N
t
-
N
0
=
-
c
ν
D
0
t
-
ν
N
t
and have provided the solution of (12) as follows:
(13)
N
t
=
N
0
∑
k
=
0
∞
-
1
k
Γ
ν
k
+
1
c
t
ν
k
.
Further, Saxena and Kalla [17] considered the following fractional kinetic equation:
(14)
N
t
-
N
0
f
t
=
-
c
ν
D
0
t
-
ν
N
t
,
R
ν
>
0
,
where N
t denotes the number density of a given species at time t, N
0
=
N
0 is the number density of that species at time t
=
0, c is a constant, and f
∈
L
0
,
∞.By applying the Laplace transform to (14), we have
(15)
L
N
t
p
=
N
0
F
p
1
+
c
ν
p
-
ν
=
N
0
∑
n
=
0
∞
-
c
ν
n
p
-
n
ν
F
p
n
∈
N
0
,
c
p
<
1
,
where the Laplace transform [24] is defined by
(16)
F
p
=
L
f
t
=
∫
0
∞
e
-
p
t
f
t
d
t
,
R
p
>
0
.The aim of this paper is to develop a new and further generalized form of the fractional kinetic equation involving generalized Bessel function of the first kind. The manifold generality of the generalized Bessel function of the first kind is discussed in terms of the solution of the above fractional kinetic equation. Moreover, the results obtained here are quite capable of yielding a very large number of known and (presumably) new results.


## 2. Solution of Generalized Fractional Kinetic Equations

In this section, we will investigate the solution of the generalized fractional kinetic equations. The results are as follows.Theorem 1.
Ifd
>
0, ν
>
0, c
,
b
,
l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(17)
N
t
-
N
0
ω
l
,
b
,
c
t
=
-
d
ν
D
0
t
-
ν
N
t
,
there holds the formula:
(18)
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
t
2
2
k
+
l
×
E
ν
,
2
k
+
l
+
1
-
d
ν
t
ν
,
where E
ν
,
2
k
+
l
+
1
· is the generalized Mittag-Leffler function [25].Proof.
The Laplace transform of the Riemann-Liouville fractional integral operator is given by [26, 27]
(19)
L
D
0
t
-
ν
f
t
;
p
=
p
-
ν
F
p
,
where F
p is defined in (16). Now, applying the Laplace transform to both sides of (17), we get
(20)
L
N
t
;
p
=
N
0
L
ω
l
,
b
,
c
t
;
p
-
d
ν
L
D
0
t
-
ν
N
t
;
p
,
N
p
=
N
0
∫
0
∞
e
-
p
t
∑
k
=
0
∞
‍
-
c
k
k
!
Γ
l
+
k
+
b
+
1
/
2
t
2
2
k
+
l
-
d
ν
p
-
ν
N
p
,
N
p
1
+
d
ν
p
-
ν
=
N
0
∑
k
=
0
∞
-
c
k
2
-
2
k
+
l
k
!
Γ
l
+
k
+
b
+
1
/
2
∫
0
∞
e
-
p
t
t
2
k
+
l
d
t
=
N
0
∑
k
=
0
∞
-
c
k
2
-
2
k
+
l
k
!
Γ
l
+
k
+
b
+
1
/
2
Γ
2
k
+
l
+
1
p
2
k
+
l
+
1
,
N
p
=
N
0
∑
k
=
0
∞
-
c
k
2
-
2
k
+
l
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
p
-
2
k
+
l
+
1
∑
r
=
0
∞
‍
1
r
-
p
/
d
-
ν
r
r
!
.
Taking Laplace inverse of (20) and using L
-
1
p
-
ν
=
t
ν
-
1
/
Γ
ν, R
ν
>
0, we have
(21)
L
-
1
N
p
=
N
0
∑
k
=
0
∞
-
c
k
2
-
2
k
+
l
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
L
-
1
∑
r
=
0
∞
-
1
r
d
ν
r
p
-
2
k
+
l
+
ν
r
+
1
,
N
t
=
N
0
∑
k
=
0
∞
-
c
k
2
-
2
k
+
l
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
∑
r
=
0
∞
-
1
r
d
ν
r
t
2
k
+
l
+
ν
r
Γ
ν
r
+
2
k
+
l
+
1
=
N
0
∑
k
=
0
∞
-
c
k
2
-
2
k
+
l
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
t
2
k
+
l
×
∑
r
=
0
∞
-
1
r
d
ν
r
t
ν
r
Γ
ν
r
+
2
k
+
l
+
1
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
t
2
2
k
+
l
×
∑
r
=
0
∞
-
1
r
d
ν
r
t
ν
r
Γ
ν
r
+
2
k
+
l
+
1
,
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
t
2
2
k
+
l
×
E
ν
,
2
k
+
l
+
1
-
d
ν
t
ν
.
This completes the proof of Theorem 1.If we setb
=
c
=
1 in (17), then generalized Bessel function ω
l
,
b
,
c
z reduces to Bessel function of the first kind J
l
z given by (2), and we arrive at the following result.Corollary 2.
Ifd
>
0, ν
>
0, l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(22)
N
t
-
N
0
J
l
t
=
-
d
ν
D
0
t
-
ν
N
t
,
there holds the formula:
(23)
N
t
=
N
0
∑
k
=
0
∞
-
1
k
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
1
t
2
2
k
+
l
×
E
ν
,
2
k
+
l
+
1
-
d
ν
t
ν
.Further, takingb
=
1 and c
=
-
1 in (17), then we obtain result of generalized fractional kinetic equation having modified Bessel function of the first kind.Corollary 3.
Ifd
>
0, ν
>
0, l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(24)
N
t
-
N
0
I
l
t
=
-
d
ν
D
0
t
-
ν
N
t
,
there holds the formula:
(25)
N
t
=
N
0
∑
k
=
0
∞
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
1
t
2
2
k
+
l
×
E
ν
,
2
k
+
l
+
1
-
d
ν
t
ν
.Lettingb
=
2 and c
=
1 in (17), then generalized Bessel function ω
l
,
b
,
c
z reduces to the spherical Bessel function of the first kind j
l
z given by (5), and we obtain the following interesting result.Corollary 4.
Ifd
>
0, ν
>
0, l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(26)
N
t
-
N
0
j
l
t
=
-
d
ν
D
0
t
-
ν
N
t
,
there holds the solution of (22)(27)
N
t
=
N
0
π
2
∑
k
=
0
∞
-
1
k
Γ
2
k
+
l
+
1
k
!
Γ
l
+
k
+
3
/
2
t
2
2
k
+
l
×
E
ν
,
2
k
+
l
+
1
-
d
ν
t
ν
.Theorem 5.
Ifd
>
0, ν
>
0, c
,
b
,
l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(28)
N
t
-
N
0
ω
l
,
b
,
c
d
ν
t
ν
=
-
d
ν
D
0
t
-
ν
N
t
,
there holds the formula:
(29)
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
d
ν
t
ν
,
where E
ν
,
2
k
ν
+
ν
l
+
1
· is the generalized Mittag-Leffler function.Proof.
The Laplace transform of the Riemann-Liouville fractional integral operator is given by [26]
(30)
L
D
0
t
-
ν
f
t
;
p
=
p
-
ν
F
p
,
where F
p is defined in (16). Now, applying the Laplace transform to both sides of (28), we get
(31)
L
N
t
;
p
=
N
0
L
ω
l
,
b
,
c
d
ν
t
ν
;
p
-
d
ν
L
D
0
t
-
ν
N
t
;
p
,
(32)
N
p
=
N
0
∫
0
∞
e
-
p
t
∑
k
=
0
∞
-
c
k
k
!
Γ
l
+
k
+
b
+
1
/
2
×
d
ν
t
ν
2
2
k
+
l
-
d
ν
p
-
ν
N
p
,
N
p
1
+
d
ν
p
-
ν
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
k
!
Γ
l
+
k
+
b
+
1
/
2
∫
0
∞
e
-
p
t
t
2
k
ν
+
ν
l
+
1
-
1
d
t
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
k
!
Γ
l
+
k
+
b
+
1
/
2
Γ
2
k
ν
+
ν
l
+
1
p
2
k
ν
+
ν
l
+
1
,
N
p
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
p
-
2
k
ν
+
ν
l
+
1
∑
r
=
0
∞
1
r
-
p
/
d
-
ν
r
r
!
.
Taking Laplace inverse of (32) and using L
-
1
p
-
ν
=
t
ν
-
1
/
Γ
ν, R
ν
>
0, we have
(33)
L
-
1
N
p
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
L
-
1
∑
r
=
0
∞
-
1
r
d
ν
r
p
-
2
k
ν
+
ν
l
+
ν
r
+
1
,
N
t
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
∑
r
=
0
∞
-
1
r
d
ν
r
t
ν
2
k
+
l
+
r
Γ
ν
r
+
2
k
ν
+
ν
l
+
1
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
t
ν
2
k
+
l
×
∑
r
=
0
∞
-
1
r
d
ν
r
t
ν
r
Γ
ν
r
+
2
k
ν
+
ν
l
+
1
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
d
ν
t
ν
2
2
k
+
l
×
∑
r
=
0
∞
-
1
r
d
ν
r
t
ν
r
Γ
ν
r
+
2
k
+
l
+
1
,
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
d
ν
t
ν
.
This completes the proof of Theorem 5.If we setb
=
c
=
1 in Theorem 5, then generalized Bessel function ω
l
,
b
,
c
z reduces to Bessel function of the first kind J
l
z, and we arrive at the special case of (28).Corollary 6.
Ifd
>
0, ν
>
0, l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(34)
N
t
-
N
0
J
l
d
ν
t
ν
=
-
d
ν
D
0
t
-
ν
N
t
,
the following result holds:
(35)
N
t
=
N
0
∑
k
=
0
∞
-
1
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
1
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
d
ν
t
ν
.On takingb
=
1 and c
=
-
1 in (28), then generalized Bessel function ω
l
,
b
,
c
z reduces to Bessel function of the first kind J
l
z, and we get the following result.Corollary 7.
Ifd
>
0, ν
>
0, l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(36)
N
t
-
N
0
I
l
d
ν
t
ν
=
-
d
ν
D
0
t
-
ν
N
t
,
the following result holds:
(37)
N
t
=
N
0
∑
k
=
0
∞
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
1
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
d
ν
t
ν
.Further, if we putb
=
2 and c
=
1 in (28), then we arrive at the following interesting result.Corollary 8.
Ifd
>
0, ν
>
0, l
,
t
∈
C, and R
l
>
-
1, then for the solution of the equation
(38)
N
t
-
N
0
j
l
d
ν
t
ν
=
-
d
ν
D
0
t
-
ν
N
t
,
the following result holds:
(39)
N
t
=
N
0
π
2
∑
k
=
0
∞
-
1
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
3
/
2
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
d
ν
t
ν
,
where j
l
z is the spherical Bessel function of the first kind.Theorem 9.
Ifa
>
0, d
>
0, ν
>
0, c
,
b
,
l
,
t
∈
C, a
≠
d, and R
l
>
-
1, then for the solution of the equation
(40)
N
t
-
N
0
ω
l
,
b
,
c
d
ν
t
ν
=
-
a
ν
D
0
t
-
ν
N
t
,
there holds the formula:
(41)
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
a
ν
t
ν
.Proof.
Applying the Laplace transform to both sides of (40), we get
(42)
L
N
t
;
p
=
N
0
L
ω
l
,
b
,
c
d
ν
t
ν
;
p
-
a
ν
L
D
0
t
-
ν
N
t
;
p
,
N
p
1
+
a
ν
p
-
ν
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
k
!
Γ
l
+
k
+
b
+
1
/
2
Γ
2
k
ν
+
ν
l
+
1
p
2
k
ν
+
ν
l
+
1
,
N
p
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
p
-
2
k
ν
+
ν
l
+
1
∑
r
=
0
∞
1
r
-
1
r
p
-
ν
r
a
ν
r
r
!
.
Taking Laplace inverse of (42), we arrive at
(43)
L
-
1
N
p
=
N
0
∑
k
=
0
∞
-
c
k
d
ν
/
2
2
k
+
l
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
×
L
-
1
∑
r
=
0
∞
-
1
r
a
ν
r
p
-
2
k
ν
+
ν
l
+
ν
r
+
1
,
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
d
ν
t
ν
2
2
k
+
l
×
∑
r
=
0
∞
-
1
r
a
ν
r
t
ν
r
Γ
ν
r
+
2
k
+
l
+
1
,
N
t
=
N
0
∑
k
=
0
∞
-
c
k
Γ
2
k
ν
+
ν
l
+
1
k
!
Γ
l
+
k
+
b
+
1
/
2
d
ν
t
ν
2
2
k
+
l
×
E
ν
,
2
k
+
l
ν
+
1
-
a
ν
t
ν
.
This completes the proof of Theorem 9.Remark 10.
The special cases for Theorem9 can be developed on similar lines to that of Corollaries 6–8, but we do not state here due to lack of space.


## 3. Conclusion

In this paper we have studied a new fractional generalization of the standard kinetic equation and derived solutions for it. It is not difficult to obtain several further analogous fractional kinetic equations and their solutions as those exhibited here by main results. Moreover, by the use of close relationships of the generalized Bessel function of the first kindω
p
z with many special functions, we can easily construct various known and new fractional kinetic equations.



---

*Source: 289387-2015-02-18.xml*