Return to Article Details Solution to unsteady fractional heat conduction in the quarter-plane via the joint Laplace-Fourier sine transforms

Solution to unsteady fractional heat conductionin the quarter-plane via the joint Laplace-Fourier sine transforms

Arman Aghili $^{*}$

August 29, 2021; accepted: September 10, 2021; published online: November 8, 2021.

In this article, the author implemented the joint transform method, for solving the boundary value problems of time fractional heat equation. We also used methods of operational nature to solve a Fokker-Planck equation with non-constant coefficients. The results reveal that the integral transform method is reliable and efficient. Some illustrative non-trivial examples are also provided.

MSC. 44A10; 42B10; 33C10; 44A20.

Keywords. Laplace transform; Fourier transform; Hankel transform; Riemann-Liouville fractional derivative; Caputo fractional derivative; Modified Bessel functions; Stieltjes transform; Heat conduction; Fokker-Planck equation; Laguerre derivative.

$^{*}$ University of Guilan, Faculty of Mathematical Sciences, Department of Applied Mathematics, Iran - Rasht P.O. Box 1841, e-mail: arman.aghili@gmail.com.

1 Introduction and Preliminaries

In this study, the author provided mathematical results that are useful to the researchers in a variety of fields. This article is devoted to studying and application of the joint Laplace-Fourier sine transform for solving time fractional diffusion equation in the quarter-plane. We also consider a Fokker-Planck equation with variable coefficients. We provided methods and results for a partial fractional differential equations which arise in applications. So far, different methods of solution have been introduced to solve partial fractional differential equations, the Laplace transform method, [ 1 ] [ 2 ] , [ 3 ] , the Fourier transform method [ 10 ] , operational method [ 4 ] , [ 6 ] . We provided methods and results for a partial fractional differential equations which arise in applications. Different methods of solution have been introduced to solve partial fractional differential equations, the Laplace transform method, [ 1 ] , [ 2 ] , [ 3 ] , the Fourier transform method [ 10 ] , operational method [ 4 ] , [ 6 ] . The diffusion equation describes the flow of heat, or a concentration of particles. In [ 11 ] , the author considered the time fractional radial diffusion in a cylinder by using the joint Laplace-finite Hankel transforms https://www.overleaf.com/project/6148a908e65ae2a25915fa5d method. In [ 12 ] , the authors considered the time-fractional diffusion-wave equation. The corresponding Green’s function was obtained in closed form for arbitrary space dimension in terms of Fox H-functions.

Remark 1.1

In [ 9 ] , the authors summarize the essential definitions and notations for the Fox H-functions.They also provide for the general Green function a representation in terms of Mellin–Barnes integrals and, consequently, in terms of Fox H-functions. â–¡

1.1 Definitions and Notations

Definition 1.2

The left Riemann-Liouville fractional derivative of order $α$ $(0 < α < 1)$ of $ϕ (t)$ is defined as follows [ 10 ]

D_{a, t}^{R - L, α} ϕ (t) = \frac{1}{Γ (1 - α)} \frac{d}{d t} \int_{a}^{t} \frac{1}{(t - ξ)^{α}} ϕ (ξ) d ξ .

1.1

Definition 1.3

The left Caputo fractional derivative of order $α$ $(0 < α < 1)$ of $ϕ (t)$ is defined as follows [ 10 ]

D_{a, t}^{C, α} ϕ (t) = \frac{1}{Γ (1 - α)} \int_{a}^{t} \frac{1}{(t - ξ)^{α}} ϕ^{'} (ξ) d ξ .

It should be pointed out that in the literature, the Riemann-Liouville and the Caputo fractional derivatives generally mean the left Riemann-Liouville and the left Caputo derivatives, respectively.

Let us recall some definitions and properties that are related to the classical continuous Fourier transform.

Definition 1.4

The Fourier transform of the function $f (x)$ , and $- \infty < x < + \infty$ is defined as follows

\begin{aligned} F {f (x)} = \frac{1}{\sqrt{2 π}} \int_{- \infty}^{+ \infty} e^{i ω x} f (x) d x := F (ω) . \end{aligned}

If $F {f (x)} = F (ω)$ , then the inverse Fourier transform $F^{- 1} {F (ω)}$ is given by

\begin{array}{r} f (x) = \frac{1}{\sqrt{2 π}} \int_{- \infty}^{+ \infty} e^{- i x ω} F (ω) d ω . \end{array}

Lemma 1.5

The following identities hold true.

$F^{- 1} (\frac{1}{a \sqrt{2}} \exp (- \frac{ω^{2}}{4 a^{2}})) = \exp (- a^{2} x^{2})$ ,
$F^{- 1} (\sqrt{\frac{2}{π}} \frac{a}{a^{2} + ω^{2}}) = \exp (- a | x |)$ ,
$F^{- 1} (\sqrt{\frac{2}{π}} \frac{2 a i ω}{(a^{2} + ω^{2})^{2}}) = x \exp (- a | x |)$ ,
$F^{- 1} (i \sqrt{\frac{2}{π}} \frac{Γ (1 - δ)}{{| ω |}^{1 - δ}} \cos (\frac{π δ}{2})) = {| x |}^{- δ} s g n (x)$ .

Proof â–¼

See [ 8 ] .

In the sequel a new class of the inverse Fourier transforms of exponential functions involving square roots are determined. Inverse Fourier transforms involving square roots arise in many areas of applied mathematics and mathematical physics.

Lemma 1.6

The following identity holds true

\begin{aligned} F^{- 1} [\frac{e^{- t \sqrt{ω^{2} + 2 i λ ω + k^{2}}}}{\sqrt{ω^{2} + 2 i λ ω + k^{2}}}; ω \to x] = \sqrt{\frac{2}{π}} e^{- λ x} K_{0} (\sqrt{(x^{2} + t^{2}) (k^{2} + λ^{2})}) . \end{aligned}

Proof â–¼

Let us assume that

\begin{aligned} F (ω) = \frac{e^{- (t \sqrt{ω^{2} + 2 i λ ω + k^{2}}})}{\sqrt{ω^{2} + 2 i λ ω + k^{2}}} . \end{aligned}

Let us consider the following well-known elementary integral

\begin{aligned} \int_{0}^{+ \infty} e^{(- a^{2} x^{2} - \frac{b^{2}}{x^{2}})} d x = \frac{\sqrt{π}}{2 a} e^{- 2 a b} . \end{aligned}

In view of the above integral, $F (ω)$ can be written as follows

\begin{aligned} F (ω) = \frac{2}{\sqrt{π}} [\frac{\sqrt{π} e^{- 2 (\frac{t}{2} \sqrt{ω^{2} + 2 i λ ω + k^{2}}})}{2 \sqrt{ω^{2} + 2 i λ ω + k^{2}}}] = \frac{2}{\sqrt{π}} \int_{0}^{+ \infty} e^{- (ω^{2} + 2 i λ ω + k^{2}) ξ^{2} - \frac{t^{2}}{4 ξ^{2}}} d ξ . \end{aligned}

Upon using Fourier inversion formula, we arrive at

\begin{aligned} F^{- 1} (F (ω)) = f (x) = \frac{1}{\sqrt{2 π}} \int_{- \infty}^{+ \infty} e^{- i x ω} (\frac{2}{\sqrt{π}} \int_{0}^{+ \infty} e^{- (ω^{2} + 2 i λ ω + k^{2}) ξ^{2} - \frac{t^{2}}{4 ξ^{2}}} d ξ) d ω, \end{aligned}

changing the order of integration, we get the following

\begin{aligned} f (x) = \frac{2}{\sqrt{π}} \int_{0}^{+ \infty} e^{- k^{2} ξ^{2} - \frac{t^{2}}{4 ξ^{2}}} (\frac{1}{\sqrt{2 π}} \int_{- \infty}^{+ \infty} e^{- i (x + 2 λ ξ^{2}) ω - ξ^{2} ω^{2}} d ω) d ξ, \end{aligned}

after evaluating the inner integral by using first part of the lemma 1.5, we have

\begin{aligned} f (x) = \frac{2}{\sqrt{π}} \int_{0}^{+ \infty} e^{- k^{2} ξ^{2} - \frac{t^{2}}{4 ξ^{2}}} (\frac{1}{ξ \sqrt{2}}) e^{- \frac{(x + 2 λ ξ^{2})^{2}}{4 ξ^{2}}} d ξ, \end{aligned}

at this stage, if we make the simple change of variable $ξ^{2} = θ$ in the above integral, after simplifying we obtain

\begin{aligned} f (x) = \frac{e^{- λ x}}{π} \int_{0}^{+ \infty} e^{- (λ^{2} + k^{2}) θ - \frac{0.5 (t^{2} + x^{2})}{θ}} \frac{d θ}{2 θ}, \end{aligned}

by means of the following integral representation for the modified Bessel’s function of the second kind of order zero (Macdonald’s function)

\begin{aligned} K_{0} (2 \sqrt{p q}) = \int_{0}^{+ \infty} e^{- (p ξ + \frac{q}{ξ})} \frac{d ξ}{2 ξ}, \end{aligned}

we have finally

\begin{aligned} f (x) = \sqrt{\frac{2}{π}} e^{- λ x} K_{0} (\sqrt{(t^{2} + x^{2}) (λ^{2} + k^{2})}) . \end{aligned}

Corollary 1.7

The following identity holds true

\begin{aligned} F^{- 1} [\frac{e^{- t \sqrt{ω^{2} + k^{2}}}}{\sqrt{ω^{2} + k^{2}}}; ω \to x] = \sqrt{\frac{2}{π}} K_{0} (k \sqrt{(x^{2} + t^{2})}) . \end{aligned}

Proof â–¼

In the above lemma 1.6 let us choose

λ = 0

we get the desired identity.

Proof â–¼

Corollary 1.8

The following integral identity holds true

\begin{aligned} \int_{- \infty}^{+ \infty} \frac{e^{- (t \sqrt{ω^{2} + 2 i λ ω + k^{2}})}}{\sqrt{ω^{2} + 2 i λ ω + k^{2}}} d ω = 2 K_{0} (t (\sqrt{k^{2} + λ^{2}})) . \end{aligned}

Proof â–¼

In the above lemma 1.6, if we take

x = 0

we arrive at the result.

Proof â–¼

Corollary 1.9

The following identity holds true

\begin{aligned} F^{- 1} [\frac{e^{- t \sqrt{ω^{2} + 2 i λ ω + k^{2}}}}{\sqrt{ω^{2} + 2 i λ ω + k^{2}}}; ω \to x] = \sqrt{\frac{2}{π}} e^{- λ x} \int_{1}^{+ \infty} e^{- η \sqrt{(λ^{2} + k^{2}) (t^{2} + x^{2})}} \frac{d η}{\sqrt{η^{2} - 1}} . \end{aligned}

Proof â–¼

Let us recall an integral representation for the modified Bessel’s function of order zero

\begin{aligned} K_{0} (ξ) = \int_{0}^{+ \infty} e^{- ξ \cosh ϕ} d ϕ . \end{aligned}

At this point, let us choose $ξ = \sqrt{(k^{2} + λ^{2}) (t^{2} + x^{2})}$ and making a change of variable $\cosh ϕ = η$ , we get

\begin{aligned} K_{0} (ξ) = K_{0} (\sqrt{(t^{2} + x^{2}) (k^{2} + λ^{2})}) = \int_{1}^{+ \infty} e^{- η \sqrt{(t^{2} + x^{2}) (k^{2} + λ^{2})}} \frac{d η}{\sqrt{η^{2} - 1}} . \end{aligned}

By replacing the above integral on the right hand side of the corollary 1.7 we obtain

\begin{aligned} F^{- 1} [\frac{e^{- t \sqrt{ω^{2} + 2 i λ ω + k^{2}}}}{\sqrt{ω^{2} + 2 i λ ω + k^{2}}}; ω \to x] = \sqrt{\frac{2}{π}} e^{- λ x} \int_{1}^{+ \infty} e^{- η \sqrt{(λ^{2} + k^{2}) (t^{2} + x^{2})}} \frac{d η}{\sqrt{η^{2} - 1}} . \end{aligned}

Remark 1.10

In the above relation, let us choose $k = 0$ , then we get the following relation

\begin{aligned} F^{- 1} [\frac{e^{- t \sqrt{ω^{2} + 2 i λ ω}}}{\sqrt{ω^{2} + 2 i λ ω}}; ω \to x] = \sqrt{\frac{2}{π}} e^{- λ x} \int_{1}^{+ \infty} e^{- η λ \sqrt{(t^{2} + x^{2})}} \frac{d η}{\sqrt{η^{2} - 1}} . \end{aligned}

Note. To the best of the author’s knowledge, in the literature the same result is obtained by using complex integration around a complicated key-hole contour in the complex plane.

Definition 1.11

The Laplace transform of the function $f (t)$ , $0 < t < + \infty$
is defined as follows

\begin{aligned} L {f (t)} = \int_{0}^{\infty} e^{- s t} f (t) d t := F (s) . \end{aligned}

If $L {f (t)} = F (s)$ , then $L^{- 1} {F (s)}$ is given by

\begin{array}{r} f (t) = \frac{1}{2 π i} \int_{c - i \infty}^{c + i \infty} e^{s t} F (s) d s = \sum_{k = 1}^{n} Res [e^{t s} F (s); s = s_{k}] . \end{array}

where $F (s)$ is analytic in the region $Re (s) > c$ with finite number of isolated singularities in the complex-plane and $s$ is a complex number with positive real part. This complex integral can be evaluated by using an appropriate contour in the complex plane, with the parameter $c$ chosen to take into account the analytic structure of the integrand, such as the presence of any poles and branch cuts. The expression in 1.5 is the inverse Laplace transform for the function $F (s)$ , and is often called the Bromwich integral.

Lemma 1.12 Gross-Levi

Let us assume that $L [f (t); t \to s] = F (s)$ and $F (s) = F (r e^{i ϕ})$ , $| ϕ | < π$ and $\int_{0}^{+ \infty} {| F (r e^{i ϕ}) |}^{2} d r < + \infty$ , then we have the following inversion formula

\begin{aligned} f (t) = \frac{1}{π} \int_{0}^{+ \infty} e^{- t r} Im [lim_{ϕ \to - π} F (r e^{i ϕ})] d r . \end{aligned}

Proof â–¼

See [ 5 ] .

The most important use of the Caputo fractional derivative is treated in initial value problems where initial conditions are expressed in terms of integer order derivatives. In this respect, it is interesting to know the Laplace transform of this kind of derivative. In the following lemma, Laplace transform of the Caputo fractional derivatives of order non integer $α$ is given.

Lemma 1.13

We have the following relations

\begin{array}{r} L {D_{0, t}^{C, α} f (t)} = s^{α} F (s) - s^{α - 1} f (0 +), 0 < α < 1 . \end{array}

and generally

\begin{array}{r} L {D_{0, t}^{C, α} f (t)} = s^{α} F (s) - \sum_{k = 0}^{k = m - 1} s^{α - 1 - k} f^{k} (0 +), m - 1 < α < m . \end{array}

Proof â–¼

See [ 7 ] .

The Laplace transform provides a useful technique for the solution of such fractional singular integro-differential equations.

Lemma 1.14

Let $L {f (t)} = F (s)$ , then the following identities hold true.

$L^{- 1} [F (s^{α}); s \to t)] =$

$= \frac{1}{π} \int_{0}^{+ \infty} f (u) (\int_{0}^{+ \infty} e^{- t r - u r^{α} \cos α π} \sin (u r^{α} \sin α π) d r) d u .$
$L^{- 1} (e^{- k \sqrt{s}}) = \frac{k}{(2 \sqrt{π})} \int_{0}^{\infty} e^{- s ξ - \frac{k^{2}}{4 ξ}} d ξ$
$e^{- ω s^{α}} = \frac{1}{π} \int_{0}^{\infty} e^{- r^{α} (ω c o s α π)} \sin (ω r^{α} \sin α π) (\int_{0}^{\infty} e^{- s τ - r τ} d τ) d r$

Proof â–¼

See [ 1 ] , [ 2 ] .

Definition 1.15

The Stieltjes transform (i.e. the second iterate of the Laplace transform) of a function $ψ (t) : R_{+} \to C$ is defined by means of the following relation

\begin{aligned} S [ψ] (s) = L [L [ψ (t); p]; s] = \int_{0}^{+ \infty} \frac{ψ (t)}{t + s} d t = Ψ (s) . \end{aligned}

Provided that the integral exists.

We have the following inversion formula for the Stieltjes transform [ 5 ] ,

\begin{aligned} S^{- 1} [Ψ (s); s \to t] = \frac{1}{π} Im [lim_{s \to t e^{- i π}} Ψ (s)], \end{aligned}

Lemma 1.16

The following integral identity holds true.

\begin{aligned} \int_{0}^{+ \infty} \frac{e^{t} K_{ν} (t)}{(t + ξ) \sqrt{t}} d t = \frac{e^{- ξ}}{\sqrt{ξ}} K_{ν} (ξ) . \end{aligned}

Proof â–¼

The left hand side of the above identity can be written in terms of the Stieltjes transform (i.e. the second itrate of the Laplace transform) of a function as below

\begin{aligned} L [L [\frac{e^{- t} K_{ν} (t)}{\sqrt{t}}; t \to η]; η \to ξ] = S [\frac{e^{- t} K_{ν} (t)}{\sqrt{t}}; t \to ξ] = \frac{π e^{- ξ}}{\cos (π ν) \sqrt{ξ}} K_{ν} (ξ) . \end{aligned}

Equivalently, we need to show that

\begin{aligned} \frac{e^{- t} K_{ν} (t)}{\sqrt{t}} = S^{- 1} [\frac{π e^{- ξ}}{\cos (π ν) \sqrt{ξ}} K_{ν} (ξ); ξ \to t] . \end{aligned}

At this point, let us evaluate the right hand side by means of the inversion formula for the Stieltjes transform [ 5 ] as follows

\begin{aligned} S^{- 1} [\frac{e^{- ξ}}{\sqrt{ξ}} K_{ν} (ξ); ξ \to t] = \frac{1}{π} Im [lim_{ξ \to t e^{- i π}} \frac{π e^{- ξ} K_{ν} (ξ)}{\cos (π ν) \sqrt{ξ}}], \end{aligned}

after simplifying we get

\begin{aligned} S^{- 1} [\frac{π e^{- ξ}}{\cos (π ν) \sqrt{ξ}} K_{ν} (ξ); ξ \to t] = \frac{1}{π} Im [\frac{π e^{- t} K_{ν} (- t)}{- i \cos (π ν) \sqrt{t}}] = Im [\frac{i e^{- t} K_{ν} (- t)}{\cos (π ν) \sqrt{t}}] . \end{aligned}

Let us use the well-known identity for the Macdonal’s function as below

\begin{aligned} K_{ν} (- t) = e^{i π ν} K_{ν} (t) = (\cos π ν + i \sin π ν) K_{ν} (t), \end{aligned}

hence, we have

\begin{aligned} S^{- 1} [\frac{π e^{- ξ}}{\cos (π ν) \sqrt{ξ}} K_{ν} (ξ); ξ \to t] = Im [\frac{e^{- t} (i \cos π ν - \sin π ν) K_{ν} (t)}{\cos (π ν) \sqrt{t}}] = \frac{e^{- t} K_{ν} (t)}{\sqrt{t}} . \end{aligned}

Proof â–¼

The operational methods provide a fast and universal mathematical tool for obtaining the solution of partial differential equations. The most commonly exponential operators which act on the function $Ψ (t)$ are as follows.

Lemma 1.17

The following exponential identities hold true.

$\exp (\pm λ \frac{d}{d t}) Ψ (t) = Ψ (t \pm λ)$
$\exp (\pm λ t \frac{d}{d t}) Ψ (t) = Ψ (t e^{\pm λ})$
$\exp (λ q (t) \frac{d}{d t}) Ψ (t) = Ψ (Q (F (t) + λ))$

where $F (t)$ is the primitive function of $\frac{1}{q (t)}$ and $Q (t)$ is the inverse function of $F (t)$ .

Proof â–¼

See [ 6 ] .

Example 1.18

Let us consider the following Fokker-Planck equation with variable coefficients as

\begin{aligned} \frac{\partial u}{\partial t} - α t^{α - 1} (\frac{\partial}{\partial x} x \frac{\partial u}{\partial x}) = - λ u . λ, γ > 0. \end{aligned}

\begin{aligned} u (x, 0) = e^{- γ x} . \end{aligned}

Note: The operator $_{L} D_{x} (.) = \frac{\partial}{\partial x} x \frac{\partial}{\partial x}$ is known as Laguerre derivative [ 7 ] .

Solution. The above partial differential equation can be written as below,

\begin{aligned} \frac{\partial u}{\partial t} = (- λ + α t^{α - 1} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}) u, \end{aligned}

in order to solve the above differential equation, we separate the variables and rewrite the above equation as follows

\begin{aligned} \frac{d u}{u} = (- λ + α t^{α - 1} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}) d t, \end{aligned}

by integrating the above equation, we get

\begin{aligned} \ln u = (- λ t + t^{α} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}) + c (x), \end{aligned}

\begin{aligned} u (x, t) = \exp [(- λ t + t^{α} \frac{\partial}{\partial x} x \frac{\partial}{\partial x})] c_{1} (x), \end{aligned}

by using the fact that $u (x, 0) = e^{- γ x} = c_{1} (x) = e^{c (x)}$ , we arrive at

\begin{aligned} u (x, t) = e^{(- λ t + t^{α} (\frac{\partial}{\partial x} x \frac{\partial u}{\partial x}))} e^{- γ x} . \end{aligned}

After simplifying, we obtain

\begin{aligned} u (x, t) = e^{- λ t} [e^{t^{α} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}}] e^{- γ x} . \end{aligned}

Let us define the operators A and B as follows

A = t^{α} \frac{\partial}{\partial x}, B = t^{α} x \frac{\partial^{2}}{\partial x^{2}} .

Then we have

t^{α} \frac{\partial}{\partial x} x \frac{\partial}{\partial x} = t^{α} \frac{\partial}{\partial x} + t^{α} x \frac{\partial^{2}}{\partial x^{2}} = A + B,

with

[A, B] = A B - B A = t^{2 α} \frac{\partial^{2}}{\partial x^{2}} = A^{2} .

In this case we have the following decomposition [ 6 ]

e^{A + B} = (1 + A) e^{B} .

From the above relation, we arrive at

\begin{aligned} u (x, t) & = e^{- λ t} [e^{t^{α} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}}] e^{- γ x} = \\ = e^{- λ t} (1 + A) [e^{B} e^{- γ x}] = e^{- λ t} (1 + t^{α} \frac{\partial}{\partial x}) [e^{t^{α} x \frac{\partial^{2}}{\partial x^{2}}} e^{- γ x}], \end{aligned}

at this point let us recall the following integral identity

e^{ξ} = \frac{1}{\sqrt{π}} \int_{- \infty}^{+ \infty} e^{- (s^{2} + 2 s \sqrt{ξ})} d s .

In the above integral let us set $ξ = t^{α} x \frac{\partial^{2}}{\partial x^{2}}$ we obtain the following operational identity

e^{t^{α} x \frac{\partial^{2}}{\partial x^{2}}} = \frac{1}{\sqrt{π}} \int_{- \infty}^{+ \infty} e^{- (s^{2} + 2 s \sqrt{x t^{α}} \frac{\partial}{\partial x})} d s = \frac{1}{\sqrt{π}} \int_{- \infty}^{+ \infty} e^{- s^{2}} [e^{- 2 \sqrt{t^{α}} s \sqrt{x} \frac{\partial}{\partial x}}] d s,

thus, we have

\begin{aligned} u (x, t) & = e^{- λ t} [e^{t^{α} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}}] e^{- γ x} \\ = e^{- λ t} (1 + t^{α} \frac{\partial}{\partial x}) [\frac{1}{\sqrt{π}} \int_{- \infty}^{+ \infty} e^{- s^{2}} [e^{(- 2 s \sqrt{t^{α}}) \sqrt{x} \frac{\partial}{\partial x}} e^{- γ x}] d s] . \end{aligned}

Then apply the identity

e^{ϕ \sqrt{x} \frac{\partial}{\partial x}} f (x) = f ([\sqrt{x} + \frac{ϕ}{2}]^{2}),

leads to

\begin{aligned} u (x, t) = e^{- λ t} (1 + t^{α} \frac{\partial}{\partial x}) [\frac{1}{\sqrt{π}} \int_{- \infty}^{+ \infty} e^{- s^{2}} e^{- γ (\sqrt{x} - \sqrt{t^{α}} s)^{2}} d s] . \end{aligned}

After evaluation of the inner integral we have

\begin{aligned} u (x, t) = \frac{e^{- λ t}}{\sqrt{π}} (1 + t^{α} \frac{\partial}{\partial x}) [\frac{\sqrt{π}}{\sqrt{1 + γ t^{α}}} e^{- (γ - \frac{γ^{2} t^{α}}{1 + γ t^{α}}) x}], \end{aligned}

after simplification

\begin{aligned} u (x, t) = \frac{e^{- λ t}}{(1 + γ t^{α})^{\frac{3}{2}}} e^{- (γ - \frac{γ^{2} t^{α}}{1 + γ t^{α}}) x} . \end{aligned}

Note: It is easy to verify that $u (x, 0) = e^{- γ x}$ .

In the next lemma, let us illustrate the use of Bromwich integral and residues theorem.

Lemma 1.19

Using Bromwich complex inversion formula to show that

\begin{aligned} L^{- 1} [\frac{K_{0} (a \sqrt{s})}{\sqrt{s} - b}] = 2 b K_{0} (a b) e^{b^{2} t} - \int_{0}^{+ \infty} e^{- t ξ^{2}} \frac{b J_{0} (a ξ) + ξ Y_{0} (a ξ)}{b^{2} + ξ^{2}} ξ d ξ . \end{aligned}

Note. The techniques introduce for evaluating the inverse Laplace transforms are adequate for a wide variety of routine applications involving the Laplace transform.

Proof â–¼

Direct application of the complex inversion formula and 1.5 leads to the following

\begin{array}{r} L^{- 1} [\frac{K_{0} (a \sqrt{s})}{\sqrt{s} - b}] = g (t) = \frac{1}{2 π i} \int_{c - i \infty}^{c + i \infty} e^{s t} \frac{K_{0} (a \sqrt{s})}{\sqrt{s} - b} d s = \end{array}

\begin{aligned} = \sum_{k = 1}^{n} Res [e^{t s} \frac{K_{0} (a \sqrt{s})}{\sqrt{s} - b}; s = s_{k}] . \end{aligned}

The transform $G (s) = [\frac{K_{0} (a \sqrt{s})}{\sqrt{s} - b}]$ has a simple pole at $s = b^{2}$ and a branch point at $s = 0$ . Then, the inverse Laplace transform will be obtained by the Gross-Levi method as follows [ 5 ]

\begin{array}{r} g (t) = lim_{s \to b^{2}} [(s - b^{2}) \frac{K_{0} (a \sqrt{s}) e^{t s}}{\sqrt{s} - b}] + \frac{1}{π} \int_{0}^{+ \infty} e^{- t r} Im [lim_{ϕ \to - π} G (r e^{i ϕ})] . \end{array}

or,

\begin{aligned} g (t) = lim_{s \to b^{2}} [(\sqrt{s} - b) (\sqrt{s} + b) \frac{K_{0} (a \sqrt{s}) e^{t s}}{\sqrt{s} - b}] + \frac{1}{π} \int_{0}^{+ \infty} e^{- t r} Im [\frac{K_{0} (a \sqrt{r e^{- i π}})}{\sqrt{r e^{- i π}} - b}] . \end{aligned}

after simplifying we get

\begin{aligned} g (t) = 2 b K_{0} (a b) e^{b^{2} t} - \frac{1}{π} \int_{0}^{+ \infty} e^{- t r} Im [\frac{K_{0} (- i a \sqrt{r})}{i \sqrt{r} + b}] d r . \end{aligned}

At this stage, we use the following well-known identity for the Bessel’s functions of the first and second kind [ 5 ]

\begin{aligned} K_{0} (- i a \sqrt{r}) = \frac{π i}{2} [J_{0} (a \sqrt{r}) + i Y_{0} (a \sqrt{r})] . \end{aligned}

therefore,

\begin{aligned} g (t) = 2 b K_{0} (a b) e^{b^{2} t} - \frac{1}{π} \int_{0}^{+ \infty} e^{- t r} Im [\frac{\frac{i π}{2} (b - i \sqrt{r}) (J_{0} (a \sqrt{r} + i Y_{0} (a \sqrt{r})}{r + b^{2}}] d r, \end{aligned}

after taking imaginary part of the fraction under integral sign, followed by making a change of variable $r = ξ^{2}$ and simplifying, we arrive at

\begin{aligned} g (t) = 2 b K_{0} (a b) e^{b^{2} t} - \int_{0}^{+ \infty} e^{- t ξ^{2}} \frac{b J_{0} (a ξ) + ξ Y_{0} (a ξ)}{b^{2} + ξ^{2}} ξ d ξ . \end{aligned}

Definition 1.20

The Hankel transform of order $ν$ of a function $f (t)$ is given by

\begin{aligned} H_{ν} [f (t); ρ] = \int_{0}^{+ \infty} f (t) t J_{ν} (ρ t) d t = F (ρ) . \end{aligned}

It is well to note that

ν

is not specified at this point and can be chosen to fit best the particular problem under consideration. In order for a transformation to be useful in solving boundary value problems, it must have an inverse. The inverse Hankel transform of a function

F (ρ)

is given by [ 4 , 5 ]

\begin{aligned} H_{ν}^{- 1} [F (ρ); t] = \int_{0}^{+ \infty} F (ρ) ρ J_{ν} (t ρ) d ρ = f (t) . \end{aligned}

Lemma 1.21

Let us define the function

h_{ν} (r) = \frac{1}{r^{ν + 1}} \frac{d}{d r} [\frac{1}{r^{ν}} y (r)] .

and apply the Hankel transform of order $ν$ , then we obtain

1. H_{ν} [h_{ν} (r); r \to ρ] = - ρ^{2} H_{ν} [y (r); r \to ρ] .

2. H_{0} [h_{0} (r); r \to ρ] = - ρ^{2} H_{0} [y (r); r \to ρ] .

As an example of the Hankel transform we will consider the problem of heat conduction formulated as follows.

Solution to time fractional heat equation via the joint Laplace-Hankel transform.

Problem 1.22

Let us solve the time fractional heat conduction equation with boundary conditions as follows

\begin{aligned} D_{t}^{C, α} u = k^{2} \frac{1}{r} \frac{\partial}{\partial r} (r \frac{\partial u}{\partial r}), \end{aligned}

\begin{aligned} u (r, 0) = f (r), 0 < α \leq 1, 0 < r < + \infty . \end{aligned}

\begin{aligned} lim_{r \to + \infty} u (r, t) = 0, \end{aligned}

\begin{array}{r} lim_{r \to 0} | u (r, t) | < + \infty . \end{array}

Note. Notice that the fractional derivative is in the Caputo sense. The constant $k^{2}$ is the coefficient of diffusion for the substance under consideration. In the most general case the coefficient of diffusion will depend on the concentration and the coordinates of the point in question.

Solution: In order to obtain a solution for equations 1.11–1.14, let us define the joint Laplace-Hankel transform as follows

\begin{aligned} H_{ν} (L [u (r, t); t \to s], r \to ρ]) = U (ρ, s) = \int_{0}^{+ \infty} r J_{0} (ρ r) (\int_{0}^{+ \infty} e^{- s t} u (r, t) d t) d r, \end{aligned}

taking the joint Laplace-Hankel transform to each term of equation 1.11 followed by second part of the lemma 1.19 and using boundary conditions 1.12, 1.13, 1.14 leads to

\begin{aligned} s^{α} U (ρ, s) + (k^{2} ρ^{2}) U (ρ, s) = s^{α - 1} F (ρ), \end{aligned}

from which we obtain

\begin{aligned} U (ρ, s) = \frac{s^{α - 1} F (ρ)}{s^{α} + k^{2} ρ^{2}} . \end{aligned}

At this point, taking the inverse joint Laplace-Hankel transform to obtain

\begin{aligned} u (r, t) = \int_{0}^{+ \infty} ρ J_{0} (r ρ) F (ρ) [L^{- 1} (\frac{s^{α - 1}}{s^{α} + k^{2} ρ^{2}})] d ρ . \end{aligned}

But, using the fact that $L^{- 1} (\frac{s^{α - 1}}{s^{α} + k^{2} ρ^{2}}) = E_{α, 1} (- k^{2} ρ^{2} t^{α})$ , we get the formal solution as follows

\begin{aligned} u (r, t) = \int_{0}^{+ \infty} ρ J_{0} (r ρ) F (ρ) E_{α, 1} (- k^{2} ρ^{2} t^{α}) d ρ . \end{aligned}

Note. $E_{α, 1} (z)$ stands for the Mittag-Leffler function in one-parameter and we have

E_{α, 1} (z) = E_{α} (z) = \sum_{k = 0}^{+ \infty} \frac{z^{k}}{Γ (α k + 1)} .

Let us consider the special case $α = 0.5$ , then we have

\begin{aligned} L^{- 1} [\frac{s^{α - 1}}{s^{α} + k^{2} ρ^{2}}] = L^{- 1} [\frac{1}{\sqrt{s} \sqrt{s} + k^{2} ρ^{2}}] = e^{k^{4} ρ^{4} t} e r f c (k^{2} ρ^{2} \sqrt{t}) . \end{aligned}

Therefore, our formal solution becomes

\begin{aligned} u (r, t) = \int_{0}^{+ \infty} ρ J_{0} (r ρ) E r f c (k^{2} ρ^{2} \sqrt{t}) F (ρ) d ρ . \end{aligned}

Let us recall that $H_{0} [f (r); r \to ρ] = \int_{0}^{+ \infty} f (ξ) ξ J_{0} (ρ ξ) ξ = F (ρ)$ then we get finally

\begin{aligned} u (r, t) = \int_{0}^{+ \infty} ξ f (ξ) [\int_{0}^{+ \infty} ρ J_{0} (r ρ) J_{0} (ξ ρ) E r f c (\sqrt{t} k^{2} ρ^{2}) d ρ] d ξ . \end{aligned}

The last step is to verify that $u (r, 0) = f (r)$ , we have

u (r, 0) = \int_{0}^{+ \infty} ξ f (ξ) (\int_{0}^{+ \infty} ρ J_{0} (ξ ρ) J_{0} (r ρ) d ρ) d ξ .

But the value of the inner integral is

\int_{0}^{+ \infty} ρ J_{0} (ξ ρ) J_{0} (r ρ) d ρ = \frac{1}{2} δ (\frac{r^{2} - ξ^{2}}{4}) .

From which we deduce that

u (r, 0) = \int_{0}^{+ \infty} ξ f (ξ) [\frac{1}{2} δ (\frac{r^{2} - ξ^{2}}{4})] d ξ .

In order to evaluate the above integral, we introduce a new change of variable $\frac{r^{2} - ξ^{2}}{4} = η$ , we get

u (r, 0) = \int_{- \infty}^{\frac{r^{2}}{4}} \sqrt{r^{2} - 4 η} f (\sqrt{r^{2} - 4 η}) \frac{δ (η) d η}{\sqrt{r^{2} - 4 η}} =

= \int_{- \infty}^{\frac{r^{2}}{4}} f (\sqrt{r^{2} - 4 η}) δ (η) d η = f (r) .

Lemma 1.23

We have the following integral identity for the Bessel’s functions

J_{0} (2 | \sinh \frac{t}{2} | ξ) = \frac{4}{π^{2}} \int_{0}^{+ \infty} \sin (t τ) \sinh (π τ) K_{i τ} (ξ) d τ .

Proof â–¼

Let us start with the well-known identity for the product of the modified Bessel’s functions

K_{ν} (x) K_{ν} (y) = \frac{π}{2 \sin π ν} \int_{0}^{+ \infty} J_{0} (\sqrt{2 x y \cosh t - (x^{2} + y^{2})}) \sinh (ν t) d t .

Let us take $x = y = ξ, ν = i τ$ then after simplifying we get

K_{i ν}^{2} (ξ) = \frac{π}{2 i \sinh π τ} \int_{0}^{+ \infty} J_{0} (\sqrt{2 ξ^{2} \cosh t - 2 ξ^{2}}) i \sin (τ t) d t,

or,

\begin{aligned} \frac{2}{π} \sinh π τ K_{i τ}^{2} (ξ) & = \int_{0}^{+ \infty} J_{0} (ξ \sqrt{2 (\cosh t - 1)}) \sin τ t d t = \\ = \int_{0}^{+ \infty} J_{0} (2 | \sinh (\frac{t}{2}) | ξ) \sin τ t d t . \end{aligned}

At this stage, by taking the inverse Fourier-sine transform, we have

\frac{4}{π^{2}} \int_{0}^{+ \infty} K_{i τ} (ξ) \sinh (π τ) \sin (t τ) d τ = J_{0} (2 | \sinh \frac{t}{2} | ξ) .

Proof â–¼

2 Solution to unsteady fractional heat conduction in the quarter-plane via the joint Laplace-Fourier sine transform.

During the last three decade, many mathematical methods were widely applied in fractal analysis. Fractional calculus used to investigate fractal functions is an important tool in this fields [ 13 ] . The main physical purpose for investigating fractional diffusion equations is to describe phenomena of anomalous diffusion.

Problem 2.1

Let us solve the time fractional heat conduction equation in two dimensions with boundary conditions as follows

\begin{aligned} D_{t}^{c, α} u = Δ u = \frac{\partial^{2} u}{\partial x^{2}} + \frac{\partial^{2} u}{\partial y^{2}}, \end{aligned}

\begin{aligned} u (x, 0, t) = 0, \end{aligned}

\begin{array}{r} u (x, y, 0) = 0, \end{array}

\begin{aligned} u (0, y, t) - u_{x} (0, y, t) = \frac{t^{α - 1}}{Γ (α)}, 0 < α \leq 1, 0 < x, y < + \infty . \end{aligned}

\begin{aligned} lim_{y \to + \infty} u (x, y, t) = 0, \end{aligned}

\begin{array}{r} lim_{x \to + \infty} | u (x, y, t) | = 0. \end{array}

Solution: In order to obtain a solution for equations 2.16–2.21, let us define the joint Laplace-Fourier sine transform as follows

\begin{aligned} F_{s} (L [u (x, t); t \to s], y \to ω]) & = U (x, w, s) \\ = \sqrt{\frac{2}{π}} \int_{0}^{+ \infty} \sin (w y) (\int_{0}^{+ \infty} e^{- s t} u (x, y, t) d t) d y, \end{aligned}

taking the joint Laplace-Fourier sine transform of 2.16 and using boundary conditions 2.17, 2.18, 2.20 yields

\begin{aligned} U_{x x} (x, ω, s) - (s^{α} + ω^{2}) U (x, ω, s) = 0, \end{aligned}

with boundary conditions as follows

\begin{aligned} U (0, ω, s) - U_{x} (0, ω, s) = \frac{1}{ω s^{α}}, \end{aligned}

\begin{array}{r} lim_{x \to + \infty} | U (x, ω, s) | = 0. \end{array}

The solution to equation 2.22 that satisfies the boundary conditions 2.23, 2.22 is as follows

\begin{aligned} U (x, ω, s) = \frac{e^{- x \sqrt{s^{α} + ω^{2}}}}{ω s^{α} (\sqrt{s^{α} + ω^{2}} + 1)} . \end{aligned}

Upon inverting the joint Laplace-Fourier sine transform, we obtain

\begin{aligned} u (x, y, t) = \sqrt{\frac{2}{π}} \int_{0}^{+ \infty} \frac{\sin y ω}{ω} (\frac{1}{2 π i} \int_{c - i \infty}^{c + i \infty} \frac{e^{t s - x \sqrt{s^{α} + ω^{2}}}}{s^{α} (\sqrt{s^{α} + ω^{2}} + 1)} d s) d ω . \end{aligned}

Let us evaluate first the complex inner integral, we have

\begin{aligned} h (x, ω, t) = L^{- 1} [\frac{e^{- x \sqrt{s^{α} + ω^{2}}}}{s^{α} (\sqrt{s^{α} + ω^{2}} + 1)}] = L^{- 1} [G (s^{α})], \end{aligned}

in order to evaluate $h (x, t)$ , let us assume that

\begin{aligned} G (s) = \frac{e^{- x \sqrt{s + ω^{2}}}}{s (\sqrt{s + ω^{2}} + 1)} . \end{aligned}

We first evaluate $L^{- 1} [G (s); s \to t] = g (x, ω, t)$ by means of the Gross-Levi method as follows

\begin{aligned} L^{- 1} [G (s); s \to t] = g (x, ω, t) = \frac{e^{- ω x}}{ω + 1} + \frac{e^{- ω^{2} t}}{π} \int_{0}^{+ \infty} e^{- t ξ} (\frac{\sqrt{ξ} \cos x \sqrt{ξ} + \sin x \sqrt{ξ}}{(1 + ξ) (ξ + ω^{2})}) d ξ . \end{aligned}

Then we evaluate $L^{- 1} [G (s^{α}); s \to t] = h (x, ω, t)$ , by means of the lemma 1.13, we get

\begin{aligned} h (x, ω, t) = \frac{1}{π} \int_{0}^{+ \infty} g (x, ω, η) (\int_{0}^{+ \infty} e^{- t r - η r^{α} \cos α π} \sin (η r^{α} \sin α π) d r) d η . \end{aligned}

Finally, we get the exact solution to diffusion equation as follows

\begin{aligned} u (x, y, t) = \sqrt{\frac{2}{π}} \int_{0}^{+ \infty} \frac{\sin y ω}{ω} h (x, ω, t) d ω . \end{aligned}

Let us consider the special case $α = 0.5$ , then we have

\begin{aligned} h (x, ω, t) = \frac{1}{π} \int_{0}^{+ \infty} g (x, ω, η) (\int_{0}^{+ \infty} e^{- t r} \sin (η \sqrt{r}) d r) d η \end{aligned}

\begin{aligned} = \frac{1}{2 t \sqrt{π t}} \int_{0}^{+ \infty} g (x, ω, η) [η e^{- \frac{η^{2}}{4 t}}] d η . \end{aligned}

Finally, we obtain

\begin{aligned} u (x, y, t) = \frac{1}{π t \sqrt{2 t}} \int_{0}^{+ \infty} \frac{\sin (y ω)}{ω} (\int_{0}^{+ \infty} η g (x, ω, η) e^{- \frac{η^{2}}{4 t}} d η) d ω, \end{aligned}

with

\begin{aligned} g (x, ω, t) = \frac{e^{- ω x}}{ω + 1} + \frac{e^{- ω^{2} t}}{π} \int_{0}^{+ \infty} e^{- t ξ} (\frac{\sqrt{ξ} \cos x \sqrt{ξ} + \sin x \sqrt{ξ}}{(1 + ξ) (ξ + ω^{2})}) d ξ . \end{aligned}

3 conclusion

The Fourier, Laplace and Hankel transforms provide a powerful method for solving certain linear differential and integral equations, and can be used for evaluating certain definite integrals. The paper is devoted to studying and application of the joint Laplace-Fourier sine transform for solving time fractional diffusion equation in the first quadrant. As a mathematical tool, the proposed method is extremely simple, attractive and concise.

Acknowledgements

The author would like to express his sincere thanks to the anonymous referee for useful comments and suggestions that lead to a significant improvement in the paper.