We consider a Cauchy problem of unidimensional
nonlinear diffusion equation on finite interval. This problem is ill-posed
and its approximate solution is unstable. We apply the He's homotopy
perturbation method (HPM) and obtain the third-order asymptotic expansion. We show that if the conductivity term in diffusion equation has
a specified condition, the above solution can be estimated. Finally, a
numerical experiment is provided to illustrate the method.
1. Introduction
The diffusion equation, one of the classical partial differential equations (PDEs), describes the process of diffusivity propagation. It has a great deal of application in different branches of sciences which have found a considerable amount of interest in recent years. This kind of equation arises naturally in a variety of models from theoretical physics, chemistry, and biology [1–8]. For instance, diffusion equations are used to investigate heat conduction, steady states and hysteresis, spatial patterns, blood oxygenation, moving fronts, pulses, and oscillations phenomena. Without any excessive simplification, these problems are all nonlinear. Therefore one needs to use a variety of different methods from different areas of mathematics such as numerical analysis, bifurcation and stability theory, similarity solutions, perturbations, topological methods, and many others, in order to study them [9–16].
Recently HPM is widely applied to linear and nonlinear problems. The method was proposed first by He in 1997 and systematical description in 2000 which is, in fact a coupling of traditional perturbation method and homotopy in topology. The application of the HPM to nonlinear problems has been developed, because this method continuously deforms the difficult problem under study into a simple one which is easy to solve. The method yields a very rapid convergence of the solution series in the most cases. Because of this rapid convergency, HPM has become a powerful mathematical tool, when it is successfully coupled with the perturbation theory. Also, HPM was used to solve variational problems by different investigators before [17–29]. One can find the recent developments of the HPM in [30–33].
This work is concerned to the nonlinear Cauchy diffusion problem and the HPM is applied to solve it. The organization of this paper is as follows. Section 2 is devoted to introduce the statement of Cauchy problem. In Section 3, we give the concepts of HPM. In Section 4, we derive the solution of Cauchy equation of nonlinear diffusion problem by HPM. In Section 5, we present an experiment wherein its numerical results illustrate the accuracy and efficiency of the proposed method, and finally in Section 6, some conclusions are considered.
2. Statement of the Cauchy Problem
Let ϕ(x,t) be a smooth function in Ω≡[0,l]×[0,T] where l and T are constant values, and f(t), g(t), a(t), and b(t) are known functions in [0,T]. Now, we assume that u(x,t) satisfies the nonlinear Cauchy diffusion equation:Au=ϕ(x,t)inΩ0≡(0,l)×(0,T),
subject to the initial conditions:u(0,t)=f(t),0≤t≤T,ux(0,t)=g(t),0≤t≤T,
where A is defined asA(u(x,t))=∂tu(x,t)-∂x{(a(t)u(x,t)+b(t))∂xu(x,t)},
such that a(t)u(x,t)+b(t) is positive [3–6], and u(x,t) is an unknown.
According to [34], we express HPM for the nonlinear problems in general case. Then, we apply this method to approximate the solution of the problem (2.1)–(2.3).
3. Description of the HPM
Suppose that A, a, b, Φ, f, and g satisfy to the above conditions. The operator A can be generally divided into two parts L and N, where L is a linear operator, and N is a nonlinear one. Therefore (2.1) can be rewritten as follows:L(u)+N(u)-ϕ(x,t)=0.
He [35] constructed a homotopy H:Ω×[0,1]→ℝ𝕖 which satisfiesH(v,p)=(1-p)[L(v)-L(v0)]+p[A(v)-ϕ(x,t)]=0,
orH(v,p)=L(v)-L(v0)+pL(v0)+p[N(v)-ϕ(x,t)]=0.
where p∈[0,1], that is called a homotopy parameter, and v0 is an initial approximation of (2.1) which satisfies initial conditions.
Hence, it is obvious thatH(v,0)=L(v)-L(v0)=0,H(v,1)=A(v)-ϕ(x,t)=0.
Now, the changing process of p from 0 to 1 is just that of H(v,p) from L(v)-L(v0) to A(v)-ϕ(x,t).
Applying the perturbation technique due to the fact that 0≤p≤1 can be considered as a small parameter, we can assume that the solution of (3.2) or (3.3) can be expressed as a series in p, as follows:v=v0+pv1+p2v2+p3v3⋯,
when p→1; (3.2) or (3.3) corresponds to (3.1) and becomes the approximate solution of (3.1). That is,u=limp→1v=v0+v1+v2+v3+⋯.
The series solution (3.6) is convergent for different terms of v, and the rate of convergence depends on A(v) [36–38].
4. Solution of Cauchy Equation of Nonlinear Diffusion Problem by HPM
Consider the nonlinear differential equation (2.1), with the indicated initial conditions (2.2). From (2.1) we have∂2u∂x2-{1b(t)∂u∂t-a(t)b(t)(∂∂x(u∂u∂x))}=-1b(t)Φ(x,t).
Then we can write (4.1) as follows:Lxu-Nu=Ψ(x,t),
where Ψ(x,t)=(-1/b(t))Φ(x,t),Lxu=∂2u/∂x2 and Nu=(1/b(t))(∂u/∂t)-(a(t)/b(t))((∂/∂x)(u∂u/∂x)) are the linear and nonlinear parts of Au, respectively.
By twice integration of (4.1) with respect to x, and applying the initial conditions (2.3), we obtain:u(x,t)-xg(t)-f(t)-∫0x∫0xNudxdx=∫0x∫0xΨ(x,t)dxdx=-1b(t)∫0x∫0xΦ(x,t)dxdx.
Consequently, we obtainu(x,t)=xg(t)+f(t)+∫0x∫0xNudxdx-1b(t)∫0x∫0xΦ(x,t)dxdx.
By HPM, let F(u)=u(x,t)-h(x,t)=0, where h(x,t)=xg(t)+f(t)+∬0xΨ(x,t)dxdx. That is, h(x,t)=xg(t)+f(t)-(1/b(t))∬0xΦ(x,t)dxdx.
Hence, we may choose a convex homotopy such that [23]H(v,p)=v(x,t)-h(x,t)-p∫0x∫0xNvdxdx=0,
whereF(u)=u(x,t)-h(x,t)=0,h(x,t)=xg(t)+f(t)+∫0x∫0xΨ(x,t)dxdx.
By using (4.5), we findv(x,t)=h(x,t)+p∫0x∫0xNvdxdx.
By combining (4.1) and (4.7), we obtainv(x,t)=xg(t)+f(t)-1b(t)∫0x∫0xΦ(x,t)dxdx+p∫0x∫0x(1b(t)∂∂tv(x,t)-a(t)b(t)∂∂x(v(x,t)∂∂xv(x,t)))dxdx,
orv0(x,t)=h(x,t)=xg(t)+f(t)-1b(t)∫0x∫0xΦ(x,t)dxdx,v1(x,t)=∫0x∫0x{1b(t)∂∂tv0-a(t)b(t)∂∂x(v0∂∂xv0)}dxdx,v2(x,t)=∫0x∫0x{1b(t)∂∂tv1-a(t)b(t)∂∂x(v0∂∂xv1+v1∂∂xv0)}dxdx,v3(x,t)=∫0x∫0x{1b(t)∂∂tv2-a(t)b(t)∂∂x(v0∂∂xv2+v1∂∂xv1+v2∂∂xv0)}dxdx,
where the above relations are obtained by equating the terms with identical powers of p in (4.8).
Therefore, the approximation solution isu(x,t)≃v0+v1+v2+v3.
In Section 5, we explain a numerical experiment. By using the HPM, an approximate solution for nonlinear diffusion equation is obtained.
5. Numerical Experiment
Let us consider the following nonlinear differential equationut-∂∂x{(16e-tu+(t+5)e-t)∂u∂x}=-73t-9,(x,t)∈[0,1]×[0,1],
with initial conditions:u(0,t)=t,0≤t≤1,ux(0,t)=0,0≤t≤1.
If we want to use our last notation, we haveΦ(x,t)=-73t-9,a(t)=16e-t,b(t)=(t+5)e-t.
Obviously, the above assumptions satisfy to consideration of aforesaid conditions. In addition, the exact solution of the problem is: u(x,t)=x2et+t.
In this experiment, we have obtained the solution of Cauchy problem at the points x = 0.1, 0.2, 0.3,…,1, where t = 0.25, 0.50, 0.75 and 1.
We construct a homotopy in the same form as we have described in Section 3: H(v,p)=u(x,t)-h(x,t)-p∫0x∫0x{et(t+5)∂u∂t-16(t+5)∂∂x(u∂u∂x)}dxdx=0.
By substituting (3.5) into the above equation, and equating the terms with identical powers of p, we have v0(x,t)=16(t+5)(6t2+30t+7etx2t+27etx2),v1(x,t)=-etx2432(t+5)3(-129etx2+7etx2t2+6etx2t+528t2-540t-5400+84t3),v2(x,t)=etx277760(t+5)5(-12177e2tx4+705e2tx4t2-2865e2tx4t+161e2tx4t3+6930etx2t2-255150etx2t-526500etx2+16830etx2t3+1890etx2t4+2520t5+28440t4+63000t3-243000t2-810000tx4),v3(x,t)=-etx278382080(t+5)7(187092e3tx6t3-131550e3tx6t2-3508884e3tx6t+2457e3tx6t4-6173667e3tx6-49829472e2tx4t2-219217320e2tx4t-217954800e2tx4+357268e2tx4t4+713160e2tx4t3+31164e2tx4t5+84672etx2t6+1155168etx2t5+35592480etx2t4-131997600etx2t3-805140000etx2t2-578340000etx2t+1360800000etx2+423360t7+6894720t6+3447360t5+1209600t4-34020000t3-680400000t2).
The exact solution, approximate solution, absolute error, relative error, L2-norm error, maximum absolute error, and maximum relative error at some time levels are presented in Tables 1 and 2.
Exact solution, approximate solution, absolute error, and relative error of u(x,t) at the time t = 0.25
x
Exact solution
Approximate solution
Absolute error
Relative error
0.1
0.2628402542
0.2628402542
0.00
0.00
0.2
0.3013610167
0.3013610083
8.20×10-9
2.75×10-8
0.3
0.3655622875
0.3655622312
5.62×10-8
1.54×10-7
0.4
0.4554440667
0.4554438471
2.20×10-7
4.81×10-7
0.5
0.5710063542
0.5710057031
6.51×10-7
1.14×10-6
0.6
0.7122491501
0.7122475258
1.62×10-6
2.28×10-6
0.7
0.8791724543
0.8791688716
3.58×10-6
4.07×10-6
0.8
1.071776267
1.071769074
7.19×10-6
6.71×10-6
0.9
1.290060588
1.290047189
1.34×10-6
1.03×10-5
1
1.534025417
1.534001959
2.34×10-5
1.52×10-5
Exact solution, approximate solution, absolute error, and relative error of u(x,t) at the time t = 0.50
x
Exact solution
Approximate solution
Absolute error
Relative error
0.1
0.5164872127
0.5164872161
3.40×10-9
6.19×10-9
0.2
0.5659488508
0.5659488481
2.70×10-9
4.77×10-9
0.3
0.6483849144
0.6483848411
7.33×10-8
1.13×10-7
0.4
0.7637954034
0.7637950723
3.31×10-7
4.33×10-7
0.5
0.9121803178
0.9121793115
1.01×10-6
1.10×10-6
0.6
1.093539658
1.093537168
2.49×10-6
2.27×10-6
0.7
1.307873423
1.307868043
5.38×10-6
4.11×10-6
0.8
1.555181613
1.555171074
1.05×10-5
6.77×10-6
0.9
1.835464230
1.835445111
1.91×10-5
1.04×10-5
1
2.148721271
2.148688717
3.26×10-5
1.51×10-5
Exact solution, approximate solution, absolute error, and relative error of u(x,t) at the time t = 0.75
x
Exact solution
Approximate solution
Absolute error
Relative error
0.1
0.7711700002
0.7711700127
1.19×10-8
1.54×10-8
0.2
0.8346800007
0.8346800302
2.88×10-8
3.45×10-8
0.3
0.9405300015
0.9405299770
2.45×10-8
2.60×10-8
0.4
1.088720003
1.088719693
3.11×10-7
2.85×10-7
0.5
1.279250004
1.279248891
1.11×10-6
8.70×10-7
0.6
1.512120006
1.512117104
2.90×10-6
1.91×10-6
0.7
1.787330008
1.787323655
6.35×10-6
3.55×10-6
0.8
2.104880011
2.104867651
1.24×10-5
5.87×10-6
0.9
2.464770014
2.464748039
2.20×10-5
8.91×10-6
1
2.867000017
2.866963750
3.62×10-5
1.26×10-5
Exact solution, approximate solution, absolute error, and relative error of u(x,t) at the time t = 1
x
Exact solution
Approximate solution
Absolute error
Relative error
0.1
1.027182818
1.027182849
3.07×10-8
2.98×10-8
0.2
1.108731273
1.108731374
1.00×10-7
9.04×10-8
0.3
1.244645364
1.244645495
1.31×10-7
1.04×10-7
0.4
1.434925092
1.434925050
4.23×10-8
2.94×10-8
0.5
1.679570457
1.679569755
7.01×10-7
4.17×10-7
0.6
1.978581458
1.978579197
2.26×10-6
1.14×10-6
0.7
2.331958096
2.331952871
5.22×10-6
2.24×10-6
0.8
2.739700370
2.739690321
1.00×10-5
3.66×10-6
0.9
3.201808281
3.201791426
1.69×10-5
5.26×10-6
1
3.718281828
3.718256914
2.49×10-5
6.70×10-6
L2-norm error ∥ue(x,t)-ua(x,t)∥2, maximum absolute error, and maximum relative error at the times t= 0.25, 0.50, 0.75, and 1.
t
∥ue(x,t)-ua(x,t)∥2
Maximum absolute error
Maximum relative error
0.25
0.001934504062
0.00002345900000
0.00001529179356
0.5
0.002320886866
0.00003255400000
0.00001515040617
0.75
0.002486049062
0.00003626700000
0.00001264980809
1
0.002172409395
0.00002491464924
0.00000670058118
In the tables fortunately, we do not have any diametrical sharp changes in our error bounds and it has no common difficulties that may appear in numerical approaches like Runge's phenomenon [39]. This means that our method works steady at all time levels. Also, the computed relative errors magnitude is acceptable and it makes our approach admissible.
Figure 1 represents the exact solution of the nonlinear diffusion problem on the interval [0,1]. As it is illustrated in Figure 2, the approximate solution gives the solution in function form. We would like to emphasize that we have presented the results in tables at some points, in order to compare our computed values with exact solution easily.
Exact solution of the nonlinear diffusion problem on the interval [0,1].
Approximate solution of the nonlinear diffusion problem on the interval [0,1].
In addition, it is possible to draw the absolute error graph because it is yield in function form too. We drew absolute error function in Figure 3, to show how little its magnitude is.
Absolute error of the nonlinear diffusion problem on the interval [0,1].
6. Conclusions
In this study, we consider the Cauchy problem of unidimensional nonlinear diffusion equation. This problem is inherently ill-posed and unsteady. If the analytical solution exists, it needs some rigid and sophisticated computation in practice. We investigate this problem with a very modern acclaimed powerful method called HPM. Our simple rapid exact approach yields good results as we have reported in Section 5. We have computed an approximate solution with acceptable error bounds which are at least of order 10-5. That makes our technique remarkable and convenient. We have used Maple 11 Packages on common home PC for all of our computations.
MehrerH.2007Berlin, GermanySpringerVázquezJ. L.2007Oxford, UKThe Clarendon Pressxxii+624Oxford Mathematical MonographsMR2286292ZBL1177.70029BurgersJ. M.1973Reidel Publishing CompanyKolokoltsovV.2000Berlin, GermanySpringerZBL1021.60011ChenZ.HuanG.MaY.2006Philadelphia, Pa, USASociety for Industrial and Applied Mathematicsxxx+531Computational Science & EngineeringMR2217767ZBL1184.11033CunninghamR. E. WilliamsR. J. J.1980New York, NY, USAPlenum PressZBL0442.90095JägerW.RannacherR.WarnatzJ.2007Berlin, GermanySpringerx+67610.1007/978-3-540-28396-6MR2265385KärgerJ.HeitjansP.2005Berlin, GermanySpringerAminataeiA.SharanM.SinghM. P. A numerical solution for the nonlinear convective facilitated-diffusion reaction problem for the process of blood oxygenation in the lungs19853182187ZBL0616.65086AminataeiA.SharanM.SinghM. P.A numerical model for the process of gas exchange in the pulmonary capillaries19871810401060ZBL0626.92007AminataeiA.SharanM.SinghM. P.Two-layer model for the process of blood oxygenation in the pulmonary capillaries-parabolic profiles in the core as well as in the plasma layer1988126601609ZBL0659.76130AminataeiA.A numerical two layer model for blood oxygenation in lungs200112456385AminataeiA.Comparison of explicit and implicit approaches to numerical solution of uni dimensional equation of diffusion2002152120 & 57AminataeiA.Blood oxygenation in the pulmonary circulation. a review20051025571AminataeiA.A numerical simulation of the unsteady convective-diffusion equation2008127387ZBL1179.65149SaiedE. A.HusseinM. M.New classes of similarity solutions of the inhomogeneous nonlinear diffusion equations1994271448674874MR129518610.1088/0305-4470/27/14/015ZBL0842.35002HeJ.-H.The homotopy perturbation method nonlinear oscillators with discontinuities20041511287292MR203796710.1016/S0096-3003(03)00341-2ZBL1039.65052HeJ.-H.Homotopy perturbation method for solving boundary value problems20063501-28788MR219932210.1016/j.physleta.2005.10.005SiddiquiA. M.MahmoodR.GhoriQ. K.Herschel-Bulkley fluid model for blood flow through a stenosed artery20063524540441610.1016/j.physleta.2005.12.033ZBL1175.76178ArielP. D.Optimal homotopy perturbation method for strongly nonlinear differential equations2010124352MarincaV.HerisanuN.Homotopy perturbation method and the natural convection flow of a third grade fluid through a circular tube201013273280El-ShahedM.Application of he’s homotopy perturbation method to volterra's integro-differential equation200562163GhasemiM.Tavassoli KhanjaniM.DavariA.Numerical solution of two-dimensional nonlinear differential equation by homotopy perturbation method200718934110.1016/j.amc.2006.11.164ZBL1170.91384DehghanM.ShakeriF.Solution of a partial differential equation subject to temperature overspecification by He' s homotopy perturbation method200775677878710.1088/0031-8949/75/6/007MR2342225ZBL1117.35326JafariM. A.AminataeiA.Numerical solution of problems in calculus of variations by homotopy perturbation method2008104828228510.1063/1.2990913ZBL1183.68133AminataeiA.JannatiQ.Using homotopy perturbation method in the numerical solution of partial differential equation2008224556ZakeriA.azakeri@kntu.ac.irJannatiQ.janati@sina.kntu.ac.irAn inverse problem for parabolic partial differential equations with nonlinear conductivity term200920094685702-s2.0-003196503310.3814/2009/468570JafariM. A.AminataeiA.Homotopy perturbation method for computing eigenelements of Sturm-Liouville two point boundary value problem200933115191524MR2544645JafariM. A.m-jafari@dena.kntu.ac.irAminataeiA.ataei@kntu.ac.irApplication of homotopy perturbation method in the solution of Fokker-Planck equation20098050550012-s2.0-003267277810.1088/0031-8949/80/05/055001ZBL1183.82049HeJ.-H.An elementary introduction to recently developed asymptotic methods and nanomechanics in textile engineering20082221348735782-s2.0-3474887716510.1142/S0217979208048668HeJ.-H.Recent development of the homotopy perturbation method2008312205209MR2432078ZBL1159.34333NoorM. A.Mohyud-DinS. T.Homotopy perturbation method for solving nonlinear higher-order boundary value problems200894395408ZBL1177.65160YildirimA.Exact solutions of nonlinear differential-difference equations by he's homotopy perturbation method200892111114FriedmanA.1964Englewood Cliffs, NJ, USAPrentice-Hallxiv+347MR0181836ZBL0178.44803HeJ.-H.Homotopy perturbation technique19991783-4257262MR171104110.1016/S0045-7825(99)00018-3ZBL0956.70017HeJ.-H.Some asymptotic methods for strongly nonlinear equations200620101141119910.1142/S0217979206033796MR2220565ZBL1102.34039LiaoS. J.2003Boca Raton, Fla, USACRC PressNayfehA. H.1981New York, NY, USAWiley-Intersciencexiv+519MR597894BerrutJ.-P.TrefethenL. N.Barycentric Lagrange interpolation200446350151710.1137/S0036144502417715MR2115059ZBL1061.65006