The Fourier law of one-dimensional heat conduction equation in fractal media is investigated in this paper. An approximate solution to one-dimensional local fractional Volterra integral equation of the second kind, which is derived from the transformation of Fourier flux equation in discontinuous media, is considered. The Picard successive approximation method is applied to solve the temperature field based on the given Mittag-Leffler-type Fourier flux distribution in fractal media. The nondifferential approximate solutions are given to show the efficiency of the present method.
1. Introduction
Engineering problems can be mathematically described by differential equations. Many initial and boundary value problems associated with differential equations can be transformed into problems of solving some approximate integral equations. Heat transfer is described by theory of integral equations. Integral equation arising in heat transfer with smooth condition is valid for continuous media [1–4]. The common methods for solving the equations of heat transfer are purely mathematical are among them; the finite difference techniques (FDT) [5], the regression analysis (RA) [6], the Adomian decomposition method (ADM) [7], the combined Laplace-Adomian method (CLAM) [8], the homotopy analysis method (HAM) [9, 10], the differential transformation method (DTM) [11], the spline-wavelets techniques (SWT) [12], the boundary element method (BEM) [13], the heat-balance integral method (HBIM) [14, 15], the variational iteration method (VIM) [16], the local fractional variational iteration method (LFVIM) [17], and the Picard successive approximation method (PSAM) [18].
On the other hand, the nanoscale heat problem can be characterized as fractal behaviors. As usual, the materials are called the Cantor materials. Heat transfer in fractal media with nonsmooth conditions is a hot topic. For example, the heat transfer equations in a medium with fractal geometry [19] and fractal domains [20] were considered. The local fractional transient heat conduction equations based upon the Fourier law within local fractional derivative arising in heat transfer from discontinuous media were presented in [21–24].
Fractional calculus was successfully used to deal with the real world problems [25–30]. There is its limit that the operators do not deal with the local fractional continuous functions (nondifferential functions). Hence, the local fractional Fourier flux [21] is not handled by using some approaches from the classical and fractional operators. This paper focuses on analytical solution to local fractional Fourier flux in fractal media by using Picard’s successive approximation method [18]. This paper is organized as follows. In Section 2, we give notations to local fractional derivative and integrals and investigate the heat transfer in fractal media. Section 3 is devoted to Picard’s successive approximation method based upon local fractional integrals. Analysis solution is shown in Section 4. Conclusions are in Section 5.
2. Heat Transfer in Fractal Media with Local Fractional Derivative
In order to study the non-differential solution for the heat problem in fractal media with local fractional derivative, we here begin with the Fourier flux equation in discontinuous media.
The temperature field reads as [21]
(1)T(x,y,z,τ)=f(x,y,z,τ)atτ>τ0andinΩ,
where f(x,y,z,τ) is local fractional continuous at fractal domain Ω.
For a given temperature field T, a local fractional temperature gradient [21] can be written as follows:
(2)∇αT=∂αT∂u1αe⃑1α+∂αT∂u2αe⃑2α+∂αT∂u3αe⃑3α,
where the local fractional partial derivative is defined by [21–24]
(3)fx(α)(x0,y)=∂αf(x,y)∂xα|x=x0=limx→x0Δα(f(x,y)-f(x0,y))(x-x0)α,
where Δα(f(x,y)-f(x0,y))≅Γ(1+α)Δ(f(x,y)-f(x0,y)).
Here, the local fractional derivative is defined on the fractal set like a Cantor set. For example, when we consider the Cantor set, we can find the local fractional derivative of discontinuous function T (however, T is a local fractional continuous function).
We consider the heat flux per unit fractal area q⃑ is proportional to the temperature gradient in fractal medium. Fourier law of heat conduction in fractal medium with local fractional derivative is expressed by [21]
(4)q⃑(x,y,z,t)=-K2α∇αT(x,y,z,t),
where K2α denotes the thermal conductivity of the fractal material, and it is related to fractal dimensions of materials. It is shown that the fractal dimensions of materials are an important characteristic value. Here, we consider the fractal Fourier flow, which is discontinuous; however, it is found that it is local fractional continuous. Like classical Fourier flow, its thermal conductivity is an approximate value for fractal one when α=1.
Fourier law of one-dimensional heat conduction equation in fractal media reads as [21]
(5)q(x,t)=-K2αdαT(x,t)dxα,atτ>τ0andinA,
where K2α denotes the thermal conductivity of the fractal materials.
When τ=τ0, from (5) we have
(6)q(x)=-K2αdαT(x)dxα,
at τ>τ0 and in A, where K2α is the thermal conductivity of the fractal materials. Namely, T is a bi-Lipschitz mapping, and shows the fractal characteristic behavior [21].
Local fractional heat conduction equation with heat generation in fractal media can be written as [21]
(7)K2α∇2αT+g-ραcα∂αT∂tα=0atτ>τ0andinΩ.
Local fractional heat conduction equation with no heat generation in fractal media is suggested as [21, 22]
(8)K2α∇2αT-ραcα∂αT∂tα=0atτ>τ0andinΩ,
where ∇2α is a local fractional Laplace operator [21].
3. The Method
In this section, we discuss the Picard successive approximation method. Meanwhile, we transfer the Fourier law of one-dimensional heat conduction equation in fractal media into the local fractional Volterra integral equation of the second kind.
3.1. Picard’s Successive Approximation Method
This method is first proposed in [18]. Here we will give a short introduction to Picard’s successive approximation method within the local fractional calculus.
In this method, we set
(9)u0(x)=f(x).
We give the first approximation u1(x) by
(10)u1(x)=f(x)+λαΓ(1+α)∫0xK(x,t)u0(x)(dt)α,
where the local fractional integral of f(x) of order α in the interval [a,b] is defined as follows [21–24]:
(11)Ib(α)af(x)=1Γ(1+α)∫abf(t)(dt)α=1Γ(1+α)limΔt→0∑j=0j=N-1f(tj)(Δtj)α
with Δtj=tj+1-tj, Δt=max{Δt0,Δt1,Δt2,…}, and [tj,tj+1], j=0,…,N-1,t0=a, andtN=b.
Here, we find that the equality u1(x) is a local fractional continuous function if f(x), K(x,t) and u0(x) are local fractional continuous functions.
Continuing in this manner, we have the infinite sequences of functions
(12)u0(x),u1(x),u2(x),…,un(x),…
such that the recurrence equations are given by
(13)un(x)=f(x)-1Γ(1+α)∫0xK(x,t)un-1(x)(dt)α,222222222222222222222222222(n=1,2,3,…),
where u0(x) is equivalent to any selected function, which is the local fractional continuous function.
Hence, we have successive approximation as follows:
(14)u1(x)=f(x)+λαΓ(1+α)∫0xK(x,t)f(t)(dt)α,⋮un(x)=f(x)+λαΓ(1+α)∫0xK(x,t)un-1(x)(dt)α.
Thus, at the limit, the solution u(x) is written as
(15)u(x)=limn→∞un(x).
3.2. An Alternative Method from Local Fractional Derivative to Local Fractional Volterra Integral Equations
We directly observe that the local fractional differential equation of α order
(16)K2αdαTdxα=q(x)(0≤x≤b)
can be written immediately as the local fractional Volterra integral equation of thee second kind in the form
(17)T(x)=T(a)+λαΓ(1+α)∫0xq(t)(dt)α,
where λ=1/K2.
The Mittag-Leffler type Fourier flux distribution in fractal media can be written as follows:
(18)q(x)=qt(t,x)=Eα(x-t)αT(t)(0≤t≤b).
Making use of (18), we can get the local fractional Volterra integral equation of the second kind in the form:
(19)T(x)=T(a)+λαΓ(1+α)∫0xEα(x-t)αT(t)(dt)α.
4. Approximate Solutions for Local Fractional Volterra Integral Equation of the Second Kind
Let us assume that the zeroth approximation is
(20)u0(x)=0.
Then the first approximation can be written as follows:
(21)u1(x)=f(x)=T(a).
Here we obtain the second approximation, which reads as
(22)u2(x)=f(x)+λαΓ(1+α)∫0xEα(x-t)αu1(t)(dt)α.
Proceeding in this manner, we have the third approximation in the following form:
(23)u3(x)=f(x)+λαΓ(1+α)∫0xEα(x-t)αu2(t)(dt)α.
Hence, continuing in this manner, we obtain
(24)un(x)=f(x)+λαΓ(1+α)∫0xEα(x-t)αun-1(t)(dt)α=f(x)+λαΓ(1+α)×∫0x(1+⋯+λ(n-1)α(x-t)(n-1)αΓ(1+(n-1)α))×Eα(x-t)αf(t)(dt)α.
Taking the limit, we have
(25)u(x)=limn→∞un(x)=T(a)+λαΓ(1+α)×∫0xEα[(x-t)α(1+λα)]T(a)(dt)α=T(a){2λα1+λαEα[xα(1+λα)]}=T(a){2K2α+1Eα[xα(1+1K2α)]},
where the term Eα[xα(1+λα)] is a Mittag-Leffler type Fourier flux distribution in fractal media, which is related to the fractal coarse-grained mass function [21, 24]. When K=1, we get u(x)=T(a)Eα[2xα]. The nondifferentiable solution of (25) for parameters K=1, α=ln2/ln3, and T(a)=1 is shown in Figure 1; the non-differentiable solution of (25) for parameters K=1, α=ln2/ln3, and T(a)=2 is shown in Figure 2; the non-differentiable solution of (25) for parameters K=1, α=ln2/ln3, and T(a)=3 is shown in Figure 3; the non-differentiable solution of (25) for parameters K=1, α=ln2/ln3, and T(a)=4 is shown in Figure 4.
The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1, α=ln2/ln3, and T(a)=1.
The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1, α=ln2/ln3, and T(a)=2.
The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1, α=ln2/ln3, and T(a)=3.
The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1, α=ln2/ln3, and T(a)=4.
5. Conclusions
This work studied the Fourier law of one-dimensional heat conduction equation in fractal media. Mittag-Leffler type Fourier flux distribution in fractal media with temperature field effect was considered. An approximation solution for the local fractional Volterra integral equation of the second kind derived from Fourier law of one-dimensional heat conduction equation for heat conduction in discontinuous media was studied by using Picard’s successive approximation method. The non-differential approximate solutions were given to show the efficiency of the present method.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by National Scientific and Technological Support Projects (no. 2012BAE09B00), the National Natural Science Foundation of China (no. 51274270), and the National Natural Science Foundation of Hebei Province (no. E2013209215).
CrankJ.NicolsonP.A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type1947435067MR0019410ZBL0029.05901ShihT. M.A literature survey on numerical heat transfer1982543694202-s2.0-002019516410.1080/10407788208913456ChoiB. C.ChurchillS. W.A technique for obtaining approximate solutions for a class of integral equations arising in radiative heat transfer19856142482-s2.0-000725804510.1016/0142-727X(85)90030-XPovstenkoY. Z.Fractional heat conduction equation and associated thermal stress20042818310210.1080/014957390523741MR2119353DehghanM.The one-dimensional heat equation subject to a boundary integral specification200732266167510.1016/j.chaos.2005.11.010MR2280112ZBL1139.35352IoannouY.FyrillasM. M.DoumanidisC.Approximate solution to Fredholm integral equations using linear regression and applications to heat and mass transfer20123681278128310.1016/j.enganabound.2012.02.006MR2913114NadeemS.AkbarN. S.Effects of heat transfer on the peristaltic transport of MHD Newtonian fluid with variable viscosity: application of adomian decomposition method200914113844385510.1016/j.cnsns.2008.09.010WazwazA.-M.MehannaM. S.The combined Laplace-Adomian method for handling singular integral equation of heat transfer2010102248252MR2745250ZBL1215.65206AbbasbandyS.The application of homotopy analysis method to nonlinear equations arising in heat transfer2006360110911310.1016/j.physleta.2006.07.065MR2288118ZBL1236.80010RaftariB.VajraveluK.Homotopy analysis method for MHD viscoelastic fluid flow and heat transfer in a channel with a stretching wall201217114149416210.1016/j.cnsns.2012.01.032MR2930314ZBL06081939JoneidiA. A.GanjiD. D.BabaelahiM.Differential transformation method to determine fin efficiency of convective straight fins with temperature dependent thermal conductivity20093677577622-s2.0-6764909520810.1016/j.icheatmasstransfer.2009.03.020CattaniC.LaserraE.Spline-wavelets techniques for heat propagation200324348549610.1080/02522667.2003.10699579MR1984452ZBL1049.65104SimõesN.TadeuA.AntónioJ.MansurW.Transient heat conduction under nonzero initial conditions: a solution using the boundary element method in the frequency domain201236456256710.1016/j.enganabound.2011.10.006MR2880089HristovJ.Approximate solutions to fractional sub-diffusion equations: the heat-balance integral method201119312292432-s2.0-7995367959210.1140/epjst/e2011-01394-2HristovJ.Heat-balance integral to fractional (half-time) heat diffusion sub-model20101422913162-s2.0-7795589511410.2298/TSCI1002291HGanjiD. D.SajjadiH.New analytical solution for natural convection of Darcian fluid in porous media prescribed surface heat flux201115222122710.2298/TSCI100424001VYangX. J.BaleanuD.Fractal heat conduction problem solved by local fractional variation iteration method201317262562810.2298/TSCI121124216YYangX.-J.Picard’s approximation method for solving a class of local fractional Volterra integral equations2012136770NigmatullinR. R.The realization of the generalized transfer equation in a medium with fractal geometry198613314254302-s2.0-002249294310.1002/pssb.2221330150DaveyK.ProsserR.Analytical solutions for heat transfer on fractal and pre-fractal domains2013371-255456910.1016/j.apm.2012.02.047MR2994201YangX. J.2012New York, NY, USAWorld ScienceZhangY. Z.YangA. M.YangX.-J.1-D heat conduction in a fractal medium: a solution by the local fractional Fourier series method201317395395610.2298/TSCI130303041ZHuM.-S.BaleanuD.YangX.-J.One-phase problems for discontinuous heat transfer in fractal media201320133358473MR303037410.1155/2013/358473YangX. J.2011Hong KongAsian AcademicLiM.ZhaoW.On bandlimitedness and lag-limitedness of fractional Gaussian noise201339291955196110.1016/j.physa.2012.12.035MR3038268LiM.Approximating ideal filters by systems of fractional order201220126365054MR2874494ZBL1233.9204810.1155/2012/365054LiM.ZhaoW.On 1/f noise2012201223673648MR3007794ZBL1264.9406010.1155/2012/673648HeydariM. H.HooshmandaslM. R.Maalek GhainiF. M.LiM.Chebyshev wavelets method for solution of nonlinear fractional integro-differential equations in a large interval201320131248208310.1155/2013/482083RayS. S.BeraR. K.Analytical solution of a fractional diffusion equation by Adomian decomposition method2006174132933610.1016/j.amc.2005.04.082MR2204183ZBL1089.65108WangQ.Numerical solutions for fractional KdV-Burgers equation by Adomian decomposition method200618221048105510.1016/j.amc.2006.05.004MR2282548ZBL1107.65124