# Analytical Solutions of the One-Dimensional Heat Equations Arising in Fractal Transient Conduction with Local Fractional Derivative

Authors: Ai-Ming Yang; Carlo Cattani; Hossein Jafari; Xiao-Jun Yang
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

The one-dimensional heat equations with the heat generation arising in fractal transient conduction associated with local fractional derivative operators are investigated. Analytical solutions are obtained by using the local fractional Adomian decomposition method via local fractional calculus theory. The method in general is easy to implement and yields good results. Illustrative examples are included to demonstrate the validity and applicability of the new technique.


---

## Body

## 1. Introduction

The Adomian decomposition method [1–3] was applied to process linear and nonlinear problems in the fields of science and engineering. Tatari and Dehghan [4] applied Adomian decomposition method to process the multipoint boundary value problem. Wazwaz [5] used Adomian decomposition method to deal with the Bratu-type equations. Daftardar-Gejji and Jafari [6] considered Adomian decomposition method to analyze the Bagley Torvik equation. Larsson [7] presented the solution for Helmholtz equation by using the Adomain decomposition method. Tatari and coworkers [8] investigated solution for the Fokker-Planck equation by Adomain decomposition method.Fractional calculus [9–12] was applied to model the physical and engineering problems for expressions of stress-strain constitutive relations of different viscoelastic fractional order properties of materials, diffusion processes with fractional order properties, fractional order flows, analytical mechanics of fractional order discrete system vibrations [13–15], and so on. Recently, the application of Adomian decomposition method for solving the linear and nonlinear fractional partial differential equations in the fields of the physics and engineering had been established in [16, 17]. Adomian decomposition method was applied to handle the time-fractional Navier-Stokes equation [18], fractional space diffusion equation [19], fractional KdV-Burgers equation [20], linear and nonlinear fractional diffusion and wave equations [21], KdV-Burgers-Kuramoto equation [22], fractional Burgers’ equation [23], and so on. For more details on some methods for solving fractional differential equations, see [24–28].Recently, local fractional calculus theory was applied to model some nondifferentiable problems for mathematical physics (see [29–36] and the references therein). The Adomian decomposition method, as one of efficient tools for solving the linear and nonlinear differential equations, was extended to find the solutions for local fractional differential equations [37–40] and nondifferentiable solutions were obtained.The partial differential equationfs describing thermal process of fractal heat conduction were suggested in [30, 38] in the following form:
(1)∂αu(x,t)∂tα-∂2αu(x,t)∂x2α=0.
The initial and boundary conditions are
(2)u(0,t)=f(t),∂αu(0,t)∂xα=g(t),
where the operator is the local fractional differential operator [29, 30, 34, 37, 38], which is applied to model the heat conduction problems in fractal media, fractal materials, fractal fracture mechanics, fractal wave behavior, Navier-Stokes equations on Cantor sets, Schrödinger equation with local fractional derivative, and diffusion equations on cantor space-time.The one-dimensional heat equations with the heat generation arising in fractal transient conduction were considered in [30] as
(3)∂αu(x,t)∂tα-∂2αu(x,t)∂x2α=ϕ(x,t),
where ϕ(x,t) is the heat generation term.We use initial and boundary conditions as follows:(4)u(0,t)=f(t),∂αu(0,t)∂xα=g(t).
The aim of this paper is to investigate the one-dimensional heat equations with the heat generation arising in fractal transient conduction by using the local fractional Adomian decomposition method.This paper is structured as follows. In Section2, we give the basic notations and definitions of local fractional operators. In Section 3, local fractional Adomian decomposition method for heat generation arising in fractal transient conduction is presented. Three examples are shown in Section 4. Finally, Section 5 presents conclusions.


## 2. Preliminaries

In this section we present some basic definitions and notations of the local fractional operators which are used further through the paper.Let us denote local fractional continuity off(x) as
(5)f(x)∈Cα(a,b).Definition 1.
Local fractional derivative operator off(x) at the point x0 is given by [29, 30, 34–38]:
(6)f(α)(x0)=dαf(x)dxα|x=x0=limx→x0Δα(f(x)-f(x0))(x-x0)α,
where Δα(f(x)-f(x0))≅Γ(1+α)Δ(f(x)-f(x0)) and f(x)∈Cα(a,b).
Local fractional derivative of high order and local fractional partial derivative of high order are written in the form [29, 30, 38]
(7)f(kα)(x)=Dx(α)⋯Dx(α)︷ktimesf(x),(8)∂kα∂xkαf(x,y)=∂α∂xα⋯∂α∂xα︷ktimesf(x,y),
respectively.
As inverse of local fractional differential operator, the local fractional integral operator off(x) in the interval [a,b] is defined as [29, 30, 36–38]
(9)Iab(α)f(x)=1Γ(1+α)∫abf(t)(dt)α=1Γ(1+α)limΔt→0∑j=0j=N-1f(tj)(Δtj)α,
where a partition of the interval [a,b] is denoted as Δtj=tj+1-tj, Δt=max{Δt0,Δt1,Δtj,…} and j=0,…,N-1, t0=a, and tN=b.
The properties are only presented as follows [29, 30, 37]:
(10)Dx(α)[f(x)g(x)]=(Dx(α)f(x))g(x)+f(x)(Dx(α)g(x)),Iax(α)f(x)g(α)(x)=[f(x)g(x)]|ax-Iax(α)f(α)(x)g(x),Dx(α)xkαΓ(1+kα)=x(k-1)αΓ[1+(k-1)α],I0b(α)xkαΓ(1+kα)=x(k+1)αΓ[1+(k+1)α].


## 3. Analysis of the Method

Let us rewrite the heat equations with the heat generation arising in fractal transient conduction in the form(11)Lt(α)u-Lxx(2α)u=ϕ,
subject to the initial and boundary conditions
(12)u(0,t)=f(t),∂αu(0,t)∂xα=g(t),
where ∂α/∂tα and ∂2α/∂x2α symbolize Lt(α) and Lxx(2α), respectively.By defining the twofold local fractional integral operator asLxx(-2α), we have
(13)Lxx(-2α)[Lt(α)u-ϕ]=Lxx(-2α)Lxx(2α)u,
so that
(14)u=Lxx(-2α)Lt(α)u-Lxx(-2α)ϕ+xαΓ(1+α)g(t)+f(t).
Hence, we get
(15)u(x,t)=u0(x,t)+Lxx(-2α)[Lt(α)u(x,t)],
where
(16)u0(x,t)=-Lxx(-2α)ϕ+xαΓ(1+α)g(t)+f(t).
So, from (15) we have iterative formula as follows:
(17)un+1(x,t)=Lxx(-2α)[Lt(α)un(x,t)],n≥0,
where u0(x,t)=-Lxx(-2α)ϕ+(xα/Γ(1+α))g(t)+f(t).Finally, the exact solution can be constructed as follows:(18)u(x,t)=limn→∞∑i=1nui(x,t).


## 4. Illustrative Examples

Example 1.
In view of (3), we consider ϕ(x,t)=1, f(t)=tα/Γ(1+α), and g(t)=tα/Γ(1+α).
We have(19)∂αu(x,t)∂tα-∂2αu(x,t)∂x2α=1,
subject to the initial value condition
(20)u0(x,t)=-x2αΓ(1+2α)+xαΓ(1+α)tαΓ(1+α)+tαΓ(1+α).
From (19) we have the following recursive relations:
(21)un+1(x,t)=Lxx(-2α)[Lt(α)un(x,t)].
In view of (21), the first few terms of the decomposition series read
(22)u0(x,t)=-x2αΓ(1+2α)+xαΓ(1+α)tαΓ(1+α)+tαΓ(1+α),u1(x,t)=x3αΓ(1+3α)+x2αΓ(1+2α).
From (25) we get
(23)u2(x,t)=u3(x,t)=⋯=un(x,t)=0.
Therefore, the exact solution of (19) can be written as
(24)u(x,t)=x3αΓ(1+3α)+xαΓ(1+α)tαΓ(1+α)+tαΓ(1+α).
The value of the fractal-dimension order α=ln2/ln3 of the behavior of the solution is shown in Figure 1.Figure 1
Solution for the one-dimensional heat equations with a fixed valueα=ln2/ln3.Example 2.
Whenϕ(x,t)=1, f(t)=tα/Γ(1+α), and g(t)=0, we get
(25)∂αu(x,t)∂tα-∂2αu(x,t)∂x2α=1.
We give the initial value condition as follows:
(26)u0(x,t)=-x2αΓ(1+2α)+tαΓ(1+α).
From (19) we have the following recursive relations:
(27)un+1(x,t)=Lxx(-2α)[Lt(α)un(x,t)].
From (27), we have the first few terms of the decomposition series as follows:
(28)u0(x,t)=-x2αΓ(1+2α)+tαΓ(1+α),u1(x,t)=x2αΓ(1+2α).
Hence, we get
(29)u2(x,t)=u3(x,t)=⋯=un(x,t)=0.
So, the exact solution of (19) reads
(30)u(x,t)=tαΓ(1+α).
The solution with fractal-dimension order α=ln2/ln3 is shown in Figure 2.Figure 2
Solution for the one-dimensional heat equations with a fixed valueα=ln2/ln3.Example 3.
Whenϕ(x,t)=1, f(t)=0, and g(t)=tα/Γ(1+α), we get
(31)∂αu(x,t)∂tα-∂2αu(x,t)∂x2α=1.
The initial value condition is presented as follows:
(32)u0(x,t)=-x2αΓ(1+2α)+xαΓ(1+α)tαΓ(1+α).
From (19) the recursive relations follow
(33)un+1(x,t)=Lxx(-2α)[Lt(α)un(x,t)].
In view of (27), we get the few terms of the series; namely,
(34)u0(x,t)=-x2αΓ(1+2α)+xαΓ(1+α)tαΓ(1+α),u1(x,t)=x3αΓ(1+3α).
Hence, we get
(35)u2(x,t)=u3(x,t)=⋯=un(x,t)=0.
So, the exact solution of (19) reads
(36)u(x,t)=x3αΓ(1+3α)-x2αΓ(1+2α)+xαΓ(1+α)tαΓ(1+α).
Figure 3 shows the exact solution when α=ln2/ln3.Figure 3
The surface shows the exact solutionu(x,t) with a fixed value α=ln2/ln3.


## 5. Conclusions

In this work, analytical solutions for the one-dimensional heat equations with the heat generation arising in fractal transient conduction associated with local fractional derivative operators were discussed. The obtained solutions are nondifferentiable functions, which are Cantor functions and they discontinuously depend on the local fractional derivative. It is shown that the local fractional Adomian decomposition method is an efficient and simple tool for solving local fractional differential equations.



---

*Source: 462535-2013-12-17.xml*