Adomian’s decomposition method to modeling power functionally graded thermoelastic materials in heat transfer and thermal stress analysis

Hamdy M. Youssef1 , Eman A. N. Al-Lehaibi2

1Mathematics Department, Faculty of Education, Alexandria University, Alexandria, Egypt

1Mechanics Department, Faculty of Engineering, Umm Al-Qura University, Makka, Saudi Arabia

2Mathematics Department, Al-Lith University College, Umm Al-Qura University, Makka, Saudi Arabia

1Corresponding author

Vibroengineering PROCEDIA, Vol. 22, 2019, p. 188-193.
Received 21 January 2019; accepted 28 January 2019; published 15 March 2019

36th International Conference on Vibroengineering in Dubai, United Arab Emirates, March 15-17th, 2019

Copyright © 2019 Hamdy M. Youssef, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF References
Cite this article
Views 89
Reads 24
Downloads 1035

This work deals with an iteration method for numerical solving the problem of one-dimensional coupled thermoelasticity under given boundary conditions. This iteration based on the Adomian’s decomposition method. All the material properties have been considered variable on position with a power law. The numerical results have been calculated for different cases of the gradient parameter and the gradient index. The numerical results have been shown in figures. The gradient parameter and the gradient index have significant effects on the temperature increment, the strain, the stress, and the displacement.

Keywords: Adomian’s decomposition method, thermoelasticity, iteration method, power law, functionally graded.

1. Introduction

Recently, much attention has been devoted to the numerical methods in which do not require discretization of time-space variables or to the linearization of the non-linear equations [1]. Adomian introduced the decomposition method for solving linear and non-linear ordinary and partial differential equations [2-4]. This method offers accurate and computable with approximately convergent solutions to linear and non-linear partial and ordinary differential equations [5-16]. Adomian got the solutions of many bio-mathematical models of viruses, bacteria, antigens, and tumor tissues [4]. Adomian’s decomposition method (ADM) is to divide the given equation into linear and nonlinear parts of the equation [1, 13].

2. Mathematical idealizations of a FGM by using power law

This particular idealization for FGM modeling is well-known in the fracture mechanics science. The thickness h, the typical material properties Px at any point at a distance x from the reference surface has been modeled to this equation [17]:

Px=PTop1-Rxhn,   R=1-PBotPTop,   PTop=Pxx=0,   PBot=Pxx=L.

R is the gradient parameter when it vanishes all the material properties are in the standard case with constant values. n is the material gradient index which depends upon the design requirements [17].

3. Formulation the problem in general form

Consider an isotropic and thermo-elastic body in one-dimensional fill the region which is defined by Ψ=x: 0x<h  where h is the thickness of the body, and it is initially at rest and has been loaded by the harmonic thermal wave, and the surface is traction free.

The displacement components for one-dimension medium have the form [18]:

ux,t=uxx,t,      uy=uz=0.

The equation of motion:

ρx    2e t2=λx+2 μx  2e x2-γx 2θ x2.

The generalized equation of heat conduction has the form:

xKxθx=ρx CExθ t+γx T0 e t,
Kx2θx2+Kxxθx=ρx CExθ t+γx T0 e t.

The constitutive relation takes the form:

σ=λx+2μxe- γx θ,   e= u x.

In the above equations, θ=T-T0 is the temperature increment, ρx is the density, λx and μx are Lame’s parameters, Kx is the thermal conductivity, γx is a material constant given by γ(x)=3λ(x)+2μ(x)αT(x), αT(x) being the coefficient of linear thermal expansion, and CEx is the specific heat at constant strain.

4. Formulation of the problem by using the exponential law

Substitute from Eq. (7) into Eqs. (4)-(6), we get:

Kx,λx, μx, CEx, ρx, αTx=K0, λ0, μ0,C0E, ρ0,α0T1-Rxhn,
ρ0    2e t2=λ0+2 μ0-γ01-Rxhn 2θ x2,
1-Rxhn2θx2-nRh1-Rxhn-1θx=ρ0 C0EK01-Rxh2nθ t
     +γ0T0K01-Rxh2ne t,

For simplicity, we use the non-dimensional variables (we will drop the primes) [18]:

x',u',h'=c0ηx,u,h,   t'=c02ηt,   θ'=θT0,   σ=σ'λ0+2μ0,


η=ρ0C0EK0,   γ0=3λ0+2μ0α0T,   c0=λ0+2μ0ρ0,   ε1=γ0T0λ0+2μ0,   ε2=γ0ρ0C0E,

thus, we obtain:

  2e x2=  2e t2+ε11-Rxhn 2θ x2,  
2θx2=1-Rxhnθ t+ε21-Rxhne t+nRh1-Rxh-1θx,

Eq. (13) has been reduced to the form:

2θx2=1-Rxhnθ t+ε21-Rxhne t+nRh1+Rxhθx.

5. Adomian’s decomposition method (ADM)

The differential operator L is defined as following [8, 10, 15]:

Lxxθx,t=1-RxhnLtθx,t+ε2  Ltex,t+nRh1+RxhLxθx,t.

The appeared operators in the above equations are defined as:

Lt=t,   Ltt=2t2,    Lx=x,    Lxx=2x2.

Assuming that the inverse of the operator Lx-1 , Lxx-1 exists in the forms [8, 10, 15]:

Lx-1fx=0xfx1dx1,    Lxx-1fx=0x0x2fx1dx1dx2.

Thus, applying the inverse operator on both the sides of Eqs. (17) and (18), we obtain:

ex,t=e0,t+ex,txx=0+Lxx-1Lttex,t+ε11-RxhnLxxθx,t ,

We decompose the functions θx,t and ex,t as following [8, 10, 15]:

θx,t=k=0θkx,t=θ0+k=1θkx,t,    ex,t=k=0ekx,t=e0+k=1ekx,t,
e0=e0,t+ex,txx=0,    θ0=θ0,t+θx,txx=0,
     +Lxx-1Lttk=0ekx,t+ε11-Rxhn Lxxk=0θkx,t ,
     +Lxx-11-Rxhn Ltk=0θkx,t+ε21-RxhnLtk=0ekx,t+nRh1+RxhLxk=0θkx,t.

We obtain these components by ekx,t and θkx,t the recursive formulas [8, 10, 15]:

ek+1x,t=Lxx-1Lttekx,t+ε11-RxhnLxxθkx,t ,    k1,
θk+1x,t=Lxx-11-RxhnLtθkx,t+ε2  Ltekx,t
     +nRhLxx-11+RxhLxθkx,t,    k1.

The bounding plane x= 0 is thermally loaded by harmonic heat and traction free as follows:

θx,tx=0=θ0sinωt,   θx,txx=0=0,    σx,tx=0=0,    ex,txx=0=0.

Moreover, the series solutions of Eq. (22) are convergent very rapidly in real physical problems as in [12, 13]. The convergence of the series has investigated by several authors in [1, 5, 8, 9, 12-15]. In an algorithmic form, the suitable value for the tolerance Tol=  10-6 [10]:

6. The numerical results

The constants of the material properties were taken as follows [18]: K0= 386 W/(mK), α0T= 1.78×10-5K-1, C0E= 383.1 J/(kgK), η= 8886.73 s/m2, T0= 293 K, μ0=  3.86×1010 N/m2, λ0=  7.76×1010 N/m2, ρ0=  8954 kg/m3, ε1= 0.0104443, ε2= 1.60862, ω=π, θ0= 1.0, h= 1.0.

Fig. 1. θx,2.0 distribution for R= 0.0, 0.1

θx,2.0 distribution for R= 0.0, 0.1

Fig. 2.ex,2.0 distribution when R= 0.0, 0.1

ex,2.0 distribution when R= 0.0, 0.1

Fig. 3.σx,2.0 distribution when R= 0.0, 0.1

σx,2.0 distribution when R= 0.0, 0.1

Fig. 4.ux,2.0 distribution when R= 0.0, 0.1

ux,2.0 distribution when R= 0.0, 0.1

Fig. 5.θx,2.0 distribution when R= 0.0, 0.5

θx,2.0 distribution when R= 0.0, 0.5

Fig. 6.ex,2.0 distribution when R= 0.0, 0.5

ex,2.0 distribution when R= 0.0, 0.5

Fig. 7.σx,2.0 distribution when R= 0.0, 0.5

σx,2.0 distribution when R= 0.0, 0.5

Fig. 8.ux,2.0 distribution when R= 0.0, 0.5

ux,2.0 distribution when R= 0.0, 0.5

7. Conclusions

Figs. 1-8 represent the temperature increment, the strain, the stress, and the displacement distribution with various values of the gradient parameter. R= 0.0 gives the normal case of non-functionally graded material, while R0.0 performs a functionally graded material of power law with different valuesR= (0.1, 0.2). From the consideration in Eq. (1), R= 0.1 means that the ratio PBot= 90 % PTop, and R= 0.2 gives that the ratio PBot= 80 % PTop. According to the results and the figures, the parameter R has significant effects on all the stat-functions θ=θx,t, e=ex,t, σ=σxxx,t, and u=ux,t. Moreover, n has been assumed with different values through the calculations (n= 1, 5, 10) which give various materials’ designs. Figs.  1, 3, 5 and 7 show that when the value of the parameters n and R increase, the value of the temperature increment and the absolute value of the stress increase, while Figs. 2, 4, 6, and 8 show that when the value of the parameters n and R increases the value of the strain and the displacement decrease.


  1. Sweilam N. Harmonic wave generation in non linear thermoelasticity by variational iteration method and Adomian’s method. Journal of Computational and Applied Mathematics, Vol. 207, Issue 1, 2007, p. 64-72. [Publisher]
  2. Adomian G. Solving Frontier Problems of Physics: The Decomposition Method. Klumer, Boston, 1994. [Publisher]
  3. Adomian G. Nonlinear Stochastic Systems Theory and Applications to Physics. Springer Science and Business Media, 1988. [Search CrossRef]
  4. Adomian G., Cherruault Y., Abbaoui K. A nonperturbative analytical solution of immune response with time-delays and possible generalization. Mathematical and Computer Modelling, Vol. 24, Issue 10, 1996, p. 89-96. [Publisher]
  5. Ciarlet P., Jamelot E., Kpadonou F. Domain decomposition methods for the diffusion equation with low-regularity solution. Computers and Mathematics with Applications, Vol. 74, Issues 10, 2017, p. 2369-2384. [Publisher]
  6. Duz M. Solutions of complex equations with Adomian decomposition method. TWMS Journal of Applied and Engineering Mathematics, Vol. 7, Issue 1, 2017, p. 66-73. [Search CrossRef]
  7. El-Sayed S.-M., Kaya D. On the numerical solution of the system of two-dimensional Burgers’ equations by the decomposition method. Applied Mathematics and Computation, Vol. 158, Issue 1, 2004, p. 101-109. [Publisher]
  8. Górecki H., Zaczyk M. Decomposition method and its application to the extremal problems. Archives of Control Sciences, Vol. 26, Issue 1, 2016, p. 49-67. [Publisher]
  9. Kaya D., El-Sayed S.-M. On the solution of the coupled Schrödinger-KdV equation by the decomposition method. Physics Letters A, Vol. 313, Issue 1, 2003, p. 82-88. [Publisher]
  10. Kaya D., Inan I. E. A convergence analysis of the ADM and an application. Applied Mathematics and Computation, Vol. 161, Issue 3, 2005, p. 1015-1025. [Publisher]
  11. Kaya D., Yokus A. A decomposition method for finding solitary and periodic solutions for a coupled higher-dimensional Burgers equations. Applied Mathematics and Computation, Vol. 164, Issue 3, 2005, p. 857-864. [Publisher]
  12. Lesnic D. Convergence of Adomian’s decomposition method: periodic temperatures. Computers and Mathematics with Applications, Vol. 44, Issues 1-2, 2002, p. 13-24. [Publisher]
  13. Lesnic D. Decomposition methods for non-linear, non-characteristic Cauchy heat problems. Communications in Nonlinear Science and Numerical Simulation, Vol. 10, Issue 6, 2005, p. 581-596. [Publisher]
  14. Li L., Jiao L., Stolkin R., Liu F. Mixed second order partial derivatives decomposition method for large scale optimization. Applied Soft Computing, Vol. 61, 2017, p. 1013-1021. [Publisher]
  15. Mustafa I. Decomposition method for solving parabolic equations in finite domains. Journal of Zhejiang University-Science A, Vol. 6, Issue 10, 2005, p. 1058-1064. [Publisher]
  16. Vadasz P., Olek S. Convergence and accuracy of Adomian’s decomposition method for the solution of Lorenz equations. International Journal of Heat and Mass Transfer, Vol. 43, Issue 10, 2000, p. 1715-1734. [Publisher]
  17. Toudehdehghan A., Lim J., Foo K. E., Ma’arof M., Mathews J. A brief review of functionally graded materials. MATEC Web Conferences: EDP Sciences, 2017. [Publisher]
  18. Youssef H. M., El-Bary A.-A. Thermoelastic material response due to laser pulse heating in context of four theorems of thermoelasticity. Journal of Thermal Stresses, Vol. 37, Issue 12, 2014, p. 1379-1389. [Publisher]

Cited By

Waves in Random and Complex Media
Eman A. N. Al-Lehaibi
JVE Journals is rebranding to Extrica

Inspired by innovations from the previous century and the rapid growth during the last years, we are improving for excellence in your publishing experience

Read to know more
JVE Journals