We investigate the local fractional linear transport equations arising in fractal porous media by using the local fractional variational iteration method. Their approximate solutions within the nondifferentiable functions are obtained and their graphs are also shown.
1. Introduction
Transport equations [1–3] had successful applications in aeronomy [4], semiconductors [5], superconductor [6], turbulence [7], QCD [8], plasma [9], gas mixture [10], and biology [11]. The linear transport equation was written as follows [12]:
(1)∂u(x,t)∂t+a(x,t)∂u(x,t)∂x=0
subject to the initial condition
(2)u(x,0)=f(x),
where f(x) is a continuous function (differentiable function).
With the development of the transport theory in porous media [13], the disordered materials, which show fractal characteristics [14], had been investigated [15–17]. Recently, during anomalous transport [18] and evidence of its existence in point vortex flow [19], fractional theory of transport was developed by many researchers. Zaslavsky [20] and Tarasov [21] investigated the anomalous transport of fractional dynamics. Lutz presented the transport equations for Lévy stable processes with fractional derivative [22]. Uchaikin and Sibatov considered the applications in disordered semiconductors [23]. Metzler and Klafter reported developments in the description of anomalous transport based upon fractional derivative [24]. Kadem et al. studied the solutions for the fractional transport equation by using the spectral method [25].
The above results via fractional calculus are set upon the differentiable functions. However, there are many nondifferentiable functions, which do not deal with the classical and fractional calculus. More recently, the local fractional calculus developed in [26–36] was the best candidate for scientists to deal with the nondifferentiable functions. In this paper, we consider the local fractional linear transport equations arising in fractal porous media in one-dimensional case [28]:
(3)∂αu(x,t)∂tα+a(x,t)∂αu(x,t)∂xα=0,
subject to the initial condition
(4)u(x,0)=f(x),
where the velocity term a(x,t) and the quantity f(x) may be nondifferentiable functions. The purpose of the current paper is to find the nondifferentiable solutions for the local fractional linear transport equations arising in fractal porous media in one-dimensional case by using the local fractional variational iteration method [31–36]. The plan of the paper is as follows. In Section 2, the conceptions of local fractional derivatives and local fractional integrals are given. In Section 3, the idea of local fractional variational iteration method is presented. In Section 4, the nondifferentiable behaviors for solutions of local fractional linear transport equations are studied. In Section 5, the conclusions are given.
2. On the Local Fractional Calculus
In this section, we introduce the definitions of local fractional derivatives and integrals which are used in the paper.
We set the function [26, 27]
(5)f(x)∈Cα(a,b),
where
(6)|f(x)-f(x0)|<εα,
with |x-x0|<δ, for ε>0, 0<α<1 and ε∈R.
Definition 1.
Let f(x)∈Cα(a,b). We define the local fractional derivative of f(x) of order α by [26–36]
(7)dαf(x0)dxα=Δα(f(x)-f(x0))(x-x0)α,
where
(8)Δα(f(x)-f(x0))≅Γ(1+α)[f(x)-f(x0)].
Definition 2.
Let f(x)∈Cα[a,b]. We define the local fractional integral of f(x) of order α in the interval [a,b] by [26, 27, 29–36]
(9)Iab(α)f(x)=1Γ(1+α)∫abf(t)(dt)α=1Γ(1+α)limΔt→0∑j=0j=N-1f(tj)(Δtj)α,
where the partitions of the interval [a,b] are denoted as (tj,tj+1), j=0,…,N-1, t0=a and tN=b with Δtj=tj+1-tj and Δt=max{Δt0,Δt1,Δtj,…}.
From (8) and (9), there are some properties in the following form [26, 27, 31–35]:
(10)dαdxαxnαΓ(1+nα)=x(n-1)αΓ(1+(n-1)α),dαdxαEα(xα)=Eα(xα),dαdxαsinα(xα)=cosα(xα),dαdxαsinhα(xα)=coshα(xα),I0x(α)xnαΓ(1+nα)=x(n+1)αΓ(1+(n+1)α),
where
(11)Eα(xα)=∑k=0∞xαkΓ(1+kα),sinα(xα)=∑k=0∞(-1)kx(2k+1)αΓ[1+(2k+1)α],cosα(xα)=∑k=0∞(-1)kx2αkΓ(1+2αk),sinhα(xα)=∑k=0∞x(2k+1)αΓ[1+(2k+1)α],coshα(xα)=∑k=0∞x2αkΓ(1+2αk).
3. The Local Fractional Variational Iteration Method
In this section, the local fractional variational iteration method first proposed in [31] is applied to deal with the local fractional linear differential equations of order α.
Let us consider the following local fractional operator equation:
(12)Lαu(s)+Nαu(s)=0,
where the linear local fractional differential operator is defined as Lα=dα/dsα, a nonlinear local fractional operator Nα.
From (12), a correction local fractional functional can be structured as
(13)un+1(t)=un(t)+I0t(α){ξα[Lαun(s)+Nαun(s)]},
where ξα denotes a fractal Lagrange multiplier; that is, δαu~n=0 [26].
Making use of (18), the new iteration formula reads as
(14)δαun+1(t)=δαun(t)+δαI0t(α)×{ξα[Lαun(s)+Nαu~n(s)]},
which leads to
(15)(ξα(τ))(α)=0,1+ξα(τ)|τ=t=0,
where δαu~n denotes a restricted local fractional variation; that is, δαu~n=0 [26].
Therefore, from (15) the fractal Lagrange multiplier can be identified as
(16)ξα(τ)=-1.
Submitting (17) into (13), we have that
(17)un+1(t)=un(t)+I0t(α){-[Lαun(s)+Nαun(s)]}.
Consequently, the nondifferentiable solution can be written as
(18)u=limn→∞un.
For more results, see [31–36].
4. Approximate Solutions Example 3.
Consider sample local fractional linear transport equations arising in fractal porous media in the form
(19)∂αu(x,t)∂tα+∂αu(x,t)∂xα=0,
with the initial condition
(20)u(x,0)=xαΓ(1+α).
From (17) we derive the following iterative formula:
(21)un+1(x,t)=un(x,t)+I0t(α){-[∂αun(x,s)∂sα+∂αun(x,s)∂xα]},
where the initial value condition is
(22)u0(x,t)=u(x,0)=xαΓ(1+α).
From (21) the first approximate term reads as
(23)u1(x,t)=u0(x,t)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=xαΓ(1+α)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=xαΓ(1+α)-tαΓ(1+α).
In a like manner, the second approximate term is
(24)u2(x,t)=u1(x,t)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=xαΓ(1+α)-tαΓ(1+α)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=xαΓ(1+α)-tαΓ(1+α).
The third approximate term can be written as
(25)u3(x,t)=u2(x,t)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=xαΓ(1+α)-tαΓ(1+α)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=xαΓ(1+α)-tαΓ(1+α).
Continuing to calculate them in this manner, for n>1, we have that
(26)un+1(x,t)=un(x,t).
Hence, we obtain the nondifferentiable solution given by
(27)u(x,t)=xαΓ(1+α)-tαΓ(1+α),
and its graph is shown in Figure 1.
The plot of solution of (19) when α=ln2/ln3.
Example 4.
The sample local fractional linear transport equations arising in fractal porous media take the form
(28)∂αu(x,t)∂tα+∂αu(x,t)∂xα=0,
subject to the initial condition
(29)u(x,0)=Eα(xα).
In view of (17) we structure the following iterative formula:
(30)un+1(x,t)=un(x,t)+I0t(α){-[∂αun(x,s)∂sα+∂αun(x,s)∂xα]},
where the initial value condition is expressed by
(31)u0(x,t)=u(x,0)=Eα(xα).
In view of (31), we give the first approximation given by
(32)u1(x,t)=u0(x,t)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=Eα(xα)[1-tαΓ(1+α)].
Similarly, we obtain that
(33)u2(x,t)=u1(x,t)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=Eα(xα)[1-tαΓ(1+α)]+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=Eα(xα)[1-tαΓ(1+α)+t2αΓ(1+2α)],u3(x,t)=u2(x,t)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=Eα(xα)[1-tαΓ(1+α)+t2αΓ(1+2α)]+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=Eα(xα)[1-tαΓ(1+α)+t2αΓ(1+2α)-t3αΓ(1+3α)],u4(x,t)=u3(x,t)+I0t(α){-[∂αu3(x,s)∂sα+∂αu3(x,s)∂xα]}=Eα(xα)-Eα(xα)tαΓ(1+α)+I0t(α){-[∂αu3(x,s)∂sα+∂αu3(x,s)∂xα]}=Eα(xα)[1-tαΓ(1+α)+t2αΓ(1+2α)iiiiiiiiiiiiiiii-t3αΓ(1+3α)+t4αΓ(1+4α)],⋮
and so on.
Hence, the nondifferentiable solution is given by
(34)u(x,t)=limn→∞un(x,t)=limn→∞Eα(xα)[∑k=0∞t2αkΓ(1+2αk)-∑k=0∞t(2k+1)αΓ[1+(2k+1)α]]=Eα(xα)[coshα(tα)-sinhα(tα)].
together with the graph shown in Figure 2.
The plot of solution of (28) when α=ln2/ln3.
Example 5.
We now focus on the following local fractional linear transport equations arising in fractal porous media:
(35)∂αu(x,t)∂tα+∂αu(x,t)∂xα=0,
subject to the initial condition
(36)u(x,0)=cosα(xα).
In view of (17) the following iterative formula is given by
(37)un+1(x,t)=un(x,t)+I0t(α){-[∂αun(x,s)∂sα+∂αun(x,s)∂xα]},
where the initial value condition is expressed by
(38)u0(x,t)=u(x,0)=cosα(xα).
Using (37) gives the approximate terms
(39)u1(x,t)=u0(x,t)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=cosα(xα)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=cosα(xα)+sinα(xα)tαΓ(1+α),u2(x,t)=u1(x,t)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=cosα(xα)+sinα(xα)tαΓ(1+α)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=cosα(xα)+sinα(xα)tαΓ(1+α)-cosα(xα)t2αΓ(1+2α),u3(x,t)=u2(x,t)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=cosα(xα)+sinα(xα)tαΓ(1+α)-cosα(xα)t2αΓ(1+2α)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=sinα(xα)∑k=01(-1)kt(2k+1)αΓ(1+(2k+1)α)+cosα(xα)∑k=01(-1)kt2kαΓ(1+2kα),u4(x,t)=u3(x,t)+I0t(α){-[∂αu3(x,s)∂sα+∂αu3(x,s)∂xα]}=sinα(xα)∑k=01(-1)kt(2k+1)αΓ(1+(2k+1)α)+cosα(xα)∑k=01(-1)kt2kαΓ(1+2kα)+I0t(α){-[∂αu3(x,s)∂sα+∂αu3(x,s)∂xα]}=sinα(xα)∑k=02(-1)kt(2k+1)αΓ(1+(2k+1)α)+cosα(xα)∑k=02(-1)kt2kαΓ(1+2kα),⋮
and so on.
Therefore, the nondifferentiable solution of (35) can be written as
(40)u(x,t)=limn→∞un(x,t)=limn→∞[sinα(xα)∑k=0∞(-1)kt(2k+1)αΓ(1+(2k+1)α)iiiiiiiiiiiii+cosα(xα)∑k=0∞(-1)kt2kαΓ(1+2kα)]=limn→∞[sinα(xα)sinα(tα)+cosα(xα)cosα(tα)]
and its graph is given in Figure 3.
The plot of solution of (35) when α=ln2/ln3.
Example 6.
We discuss the following local fractional linear transport equations arising in fractal porous media:
(41)∂αu(x,t)∂tα+∂αu(x,t)∂xα=0,
subject to the initial condition
(42)u(x,0)=sinα(xα).
From (17) the following iterative formula can be written as
(43)un+1(x,t)=un(x,t)+I0t(α){-[∂αun(x,s)∂sα+∂αun(x,s)∂xα]},
where the initial value condition is given by
(44)u0(x,t)=u(x,0)=sinα(xα).
Making use of (43), we obtain the approximate terms
(45)u1(x,t)=u0(x,t)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=sinα(xα)+I0t(α){-[∂αu0(x,s)∂sα+∂αu0(x,s)∂xα]}=sinα(xα)-cosα(xα)tαΓ(1+α),u2(x,t)=u1(x,t)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=sinα(xα)-cosα(xα)tαΓ(1+α)+I0t(α){-[∂αu1(x,s)∂sα+∂αu1(x,s)∂xα]}=sinα(xα)-cosα(xα)tαΓ(1+α)-sinα(xα)t2αΓ(1+2α),u3(x,t)=u2(x,t)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=sinα(xα)-cosα(xα)tαΓ(1+α)-sinα(xα)t2αΓ(1+2α)+I0t(α){-[∂αu2(x,s)∂sα+∂αu2(x,s)∂xα]}=sinα(xα)-cosα(xα)tαΓ(1+α)-sinα(xα)t2αΓ(1+2α)+cosα(xα)t3αΓ(1+3α),u4(x,t)=u3(x,t)+I0t(α){-[∂αu3(x,s)∂sα+∂αu3(x,s)∂xα]}=sinα(xα)-cosα(xα)tαΓ(1+α)-sinα(xα)t2αΓ(1+2α)+cosα(xα)t3αΓ(1+3α)+I0t(α){-[∂αu3(x,s)∂sα+∂αu3(x,s)∂xα]}=sinα(xα)∑k=02(-1)kt2kαΓ(1+2kα)-cosα(xα)∑k=01(-1)kt(2k+1)αΓ(1+(2k+1)α),⋮
and so on.
Therefore, the nondifferentiable solution of (41) reads as
(46)u(x,t)=limn→∞un(x,t)=limn→∞[sinα(xα)∑k=0∞(-1)kt2kαΓ(1+2kα)iiiiiiiiiiiii-cosα(xα)∑k=0∞(-1)kt(2k+1)αΓ(1+(2k+1)α)]=sinα(xα)cosα(tα)-cosα(xα)sinα(tα)
and its graph is given in Figure 4.
The plot of solution of (41) when α=ln2/ln3.
5. Conclusions
Local fractional calculus theory is a tool for modeling the nondifferentiable problems for science and engineering. In this work we studied the local fractional linear transport equations arising in fractal porous media by using the local fractional variational iteration method. The solutions with nondifferentiable functions were also obtained and some examples were also discussed. These results show the reliabilities and efficiencies of the proposed local fractional variational iteration method.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work was supported by the National Natural Science Foundation of China under Grant no. 71173060. It was also supported by the China Postdoctoral Science Foundation under Grant no. 2013M541351.
AgoshkovV.1998Boston, Mass, USABirkhäuserxviii+278Modeling and Simulation in Science, Engineering and TechnologyMR1638817DumasL.GolseF.Homogenization of transport equations20006041447147010.1137/S0036139997332087MR1760042ZBL0964.350162-s2.0-0033721064WeinanE.Homogenization of linear and nonlinear transport equations199245330132610.1002/cpa.3160450304MR1151269ZBL0794.35014SchunkR. W.Transport equations for aeronomy19752334374852-s2.0-0001761745BlotekjaerK.Transport equations for electrons in two-valley semiconductors19701713847Betbeder-MatibetO.NozieresP.Transport equations in clean superconductors19695133924172-s2.0-2342648426DalyB. J.HarlowF. H.Transport equations in turbulence19701311263426492-s2.0-36849100388ElzeH.-T.GyulassyM.VasakD.Transport equations for the QCD quark Wigner operator19862763-47067282-s2.0-26144449149MikhailovskiiA. B.TsypinV. S.Transport equations of plasma in a curvilinear magnetic field1984244335354TanenbaumB. S.Transport equations for a gas mixture1965846836862-s2.0-22644437733PerthameB.2006Berlin, GermanySpringerPetrovaG.PopovB.Linear transport equations with discontinuous coefficients1999249-101849187310.1080/03605309908821484MR1708110ZBL0992.351042-s2.0-0033474746InghamD. B.PopI.20053Oxford, UKElsevierStanleyH. E.Application of fractal concepts to polymer statistics and to anomalous transport in randomly porous media1984365-68438602-s2.0-000164502310.1007/BF01012944SahimiM.Fractal and superdiffusive transport and hydrodynamic dispersion in heterogeneous porous media19931313402-s2.0-002768115710.1007/BF00613269RigbyS. P.GladdenL. F.NMR and fractal modelling studies of transport in porous media19965110226322722-s2.0-0030131264ColotinM.PompilianG. O.NicaP.GurluiS.PaunV.AgopM.Fractal transport phenomena through the scale relativity model200911621571642-s2.0-69849103385KlagesR.RadonsG.SokolovI. M.2008New York, NY, USAJohn Wiley & SonsLeonciniX.KuznetsovL.ZaslavskyG. M.Evidence of fractional transport in point vortex flow20041922592732-s2.0-004259299210.1016/S0960-0779(03)00040-7MR2007456ZBL1068.76012ZaslavskyG. M.Chaos, fractional kinetics, and anomalous transport200237164615802-s2.0-003688793610.1016/S0370-1573(02)00331-9MR1937584ZBL0999.82053TarasovV. E.Transport equations from Liouville equations for fractional systems20062033413532-s2.0-3264444004210.1142/S0217979206033267MR2206188ZBL1130.82335LutzE.Fractional transport equations for Lévy stable processes20018611220822112-s2.0-003584815310.1103/PhysRevLett.86.2208UchaikinV. V.SibatovR. T.Fractional theory for transport in disordered semiconductors2008134715727MR2381497ZBL1221.821412-s2.0-3554897345910.1016/j.cnsns.2006.07.008MetzlerR.KlafterJ.The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics20043731R161R208MR2090004ZBL1075.820182-s2.0-404315147710.1088/0305-4470/37/31/R01KademA.LuchkoY.BaleanuD.Spectral method for solution of the fractional transport equation2010661103115MR2732183ZBL1237.820412-s2.0-7995372050110.1016/S0034-4877(10)80026-6YangX.-J.2012New York, NY, USAWorld Science PublisherYangX. J.2011Hong KongAsian Academic PublisherYangX. J.BaleanuD.HeJ. H.Transport equations in fractal porous media within fractional complex transform method2013144287292ZhaoC.-G.YangA.-M.JafariH.HaghbinA.The Yang-Laplace transform for solving the IVPs with local fractional derivative20142014538645910.1155/2014/386459MR3166607YangA.-M.ZhangY.-Z.CattaniC.XieG.-N.RashidiM. M.ZhouY.-J.YangX.-J.Application of local fractional series expansion method to solve Klein-Gordon equations on Cantor sets20142014637274110.1155/2014/372741MR3182278YangX. J.BaleanuD.Fractal heat conduction problem solved by local fractional variation iteration method2013172625628YangX. J.BaleanuD.Local fractional variational iteration method for Fokker-Planck equation on a Cantor set201323238YangA. M.ChenZ. S.SrivastavaH. M.YangX. J.Application of the local fractional series expansion method and the variational iteration method to the Helmholtz equation involving local fractional derivative operators20132013610.1155/2013/259125259125BaleanuD.Tenreiro MachadoJ. A.CattaniC.BaleanuM. C.YangX.-J.Local fractional variational iteration and decomposition methods for wave equation on Cantor sets within local fractional operators20142014653504810.1155/2014/535048MR3166628HeJ. H.Local fractional variational iteration method for fractal heat transfer in silk cocoon hierarchy2013411520LiuC. F.KongS. S.YuanS. J.Reconstructive schemes for variational iteration method within Yang-Laplace transform with application to fractal heat conduction problem2013173715721