Hyperbolic two-temperature generalized thermoelasticity with fractional order strain of solid cylinder
Alaa A. El-Bary1
1Basic and Applied Science Institute, Arab Academy for Science, Technology and Maritime Transport, P.O. Box 1029, Alexandria, Egypt
1National Committee for Mathematics, Academy of Scientific Research and Technology, Cairo, Egypt
Journal of Engineering and Thermal Sciences, (in Press).
Received 26 March 2021; received in revised form 6 July 2021; accepted 6 August 2021; published 13 September 2021
A novel mathematical model of thermoelastic of a homogenous isotropic solid cylindrical infinite medium has been constructed in this paper. Thermally shocked is the bounding surface of the cylinder. In the sense of the hyperbolic two-temperature generalized thermoelasticity with fractional stress theory, the governing equations have been taken. Different values of the fractional order and two-temperature parameters have shown numerical results for the dynamical and conductive temperature increment, strain, displacement, and average stress, which are graphically applicable to all the functions studied. The fractional-order parameter has significant effects on stress and displacement distributions, while it has little effect on the dynamical and conductive temperatures increment and significant effects on all studied functions as well as on the two-temperature parameter. The two-temperature hyperbolic model is precious and effective.
Keywords: hyperbolic two-temperature generalized thermoelasticity, cylindrical body, fractional-order strain, two-temperature generalized thermoelasticity.
Finding the best mathematical model that simulates the behavior of these materials close to experience is one of the key problems in materials and solid sciences. The mathematical model produces results that are completely consistent with experimental findings, and researchers in this area are primarily interested in the physical behavior of these compounds. Various mathematical models have been proposed by scientists and researchers to explain the propagation of mechanical waves and heat in solid and elastic materials. However, not all these models are suitable, since one of the metrics of an ideal model is to include the speed of progression of mechanical and thermal waves with finite values as experimental results. Mathematical models of thermomechanical transformation through elastic materials involve a large amount of space that cannot be included in a single study. However, new models that seek to demonstrate new observations that are very close to laboratory studies and are generally consistent with the physical behavior of thermo-elastic materials should be discussed.
Chen and Gurtin  have introduced a thermoelastic theory based on two different temperatures, conductive and dynamic temperature. The difference between these two temperatures is proportional to the heat source. Warren and Chen  studied wave propagation in the two-temperature thermo-elasticity principle. However, there is no investigation of this theory until Youssef  has updated this theory and introduced a two-temperature generalized thermo-elasticity model. Youssef and Elsibai enhanced the two-temperature thermoelasticity principle without a Green-Naghdi model energy dissipation . Youssef and many other authors have used this model in many applications and research [5-8]. Youssef and El-Bary approved that the two-temperature generalized thermoelasticity model does not provide a finite speed of propagation of the thermal waves . Youssef and El-Bary have therefore updated the model and have introduced a new two-temperature model based on different heat piping rules, known as hyperbolic, generalized 2-temperature thermoelasticity .
Youssef suggested in this model that the difference between the conductive temperature acceleration and the acceleration of the thermal temperature during the material transition is proportionate to the heater supply. The speed of the spread of thermal waves is limited within this model.
Additional new models of thermoelasticity based on the fractional calculus principle were introduced. Magin and Royston  have introduced the first model that applies the fractional deformation derivative characterizing the material's behavior. The zero-order of the derivative is for a Hookean solid while the one-order is for a Newtonian fluid. The thermoelastic and thermoviscoelastic materials are sitting in the intermediate range with a fractional order parameter between zero and one . C Li et al. studied the generalized thermoelastic diffusion problems with fractional-order strain [11, 12]. Li, Chenlin, et al. studied the generalized piezoelectric thermoelasticity problems with strain rate and transient thermo‐electromechanical responses analysis .
Recently, many applications used the fractional-order definition with various construction and different models [14-22]. Depending on the fractional strain, Youssef  presented another new theory of generalized thermoelectricity. It is seen as another shift in the stress-strain relationship between Duhamel and Neumann. Youssef solved an application of thermoelasticity with fractional order strain for a one-dimensional half-space, in the context of models of Biot, Lord-Shulman, Green-Lindsay, and Green-Naghdi type II . Youssef solved many applications of thermoelasticity of infinite thermoelastic cylindrical medium [8, 24].
In the present investigation, we will study the induced conductive temperature, dynamic temperature, deformation, displacement, and stress fields in a thermoelastic one-dimensional infinite cylindrical body. The work under the purview of the theory of thermoelasticity with fraction order strain is based on hyperbolic two-temperature heat conduction laws.
The main target of this work is to discuss and study the effects of the fractional-order strain parameter and hyperbolic two-temperature parameter on thermomechanical waves through a thermoelastic body.
2. The governing equations
Consider a perfect thermal conducting, isotropic, elastic, and cylindrical body occupies the region . We apply a cylindrical coordinates system in which the -axis was lying along the axis of the cylinder, as in Fig. 1, and due to symmetry, the medium is one-dimensional. Moreover, all the state-functions will depend on the radial distance and the time . It is proposed that there are no body forces or heat sources in the medium and initially at rest.
Thus, the governing equations in cylindrical one-dimensional with fractional order strain can be written as [8, 24].
The equation of motion:
The heat equation:
where (m/s) is constant (two-temperature parameter).
Fig. 1. The isotropic homogeneous thermoelastic solid cylinder
The constitutive relations will take the forms:
and is the volume dilatation and satisfies the relation:
In the above equations, we apply the definition of the Riemann–Liouville fractional integral written in a convolution type form :
which gives Caputo fractional derivatives in the form:
For convenience, we will use the following non-dimensional variables [8, 24]:
where and .
Eqs. (1)-(6) take the following forms (the primes are suppressed for simplicity):
where , , , , , .
We apply the Laplace transform defined as:
The Laplace transform of the fractional derivative is defined as :
We assume the following initial conditions:
Thus, the Eqs. (12)-(17) have the forms:
We can re-write Eqs. (22)-(24) in the forms:
where , , , .
By inserting Eq. (31) into the Eq. (29), we obtain:
Substitute from Eq. (31) into the Eq. (30), we get:
where and .
Substitute from Eq. (33) into the Eq. (31), we get:
where and .
By eliminating between Eq. (32) and (33), we obtain:
By eliminating between Eqs. (32) and (33), we get:
where and .
The bounded solutions at of the Eqs. (35) and (36) take the forms:
where is the modified Bessel function of the first kind of order zero, and , are the roots of the characte, ristic equation:
For the thermal boundary condition, we consider the edge of the cylinder to be thermally loaded as follows:
where is the function of the thermal loading and is constant (the thermal loading strength). For the mechanical boundary conditions, we consider the edge of the cylinder is connected to a rigid foundation which can prevent any displacement which leads to zero volumetric deformation as follows:
Applying Laplace transform on the Eqs. (40) and (41), we obtain:
Applying the boundary conditions in the Eqs. (37) and (38), we get the following system:
By solving the system in Eq. (43) and Eq. (44), we get:
Hence, we have:
To get the displacement function, we use the Eqs. (28) and (46) as follows:
We use the well-known relation of Bessel integration as follows:
where is the modified Bessel function of the first kind of first order.
Hence, we obtain:
To obtain the dynamical temperature, we substitute from equations (45) and (46) in the Eqs. (34) to be in the form:
To get the stress distribution and to avoid the singularity when , we take the average of the three principal stresses components as follows:
We re-write the Eq. (52) to be in the form:
Substitute from Eqs. (46) and (51) in the Eq. (53), hence, we get:
The conducting and dynamical temperature increments, pressure, displacement, and average stress can both be measured using Riemann-sum approximation techniques over time. This method is used to invert some function’s Laplace transformation :
where “” is an imaginary number unit, and “” is the real part. For convergence with a faster process, many numerical experiments have shown that the value satisfies the relation .
3. Numerical results and discussion
Copper is known as the thermoelastic material for the numerical results and for the various physical constants we take the following values: 386 kg m k-1 s-3, 1.78×10-5 k-1, 8954 kg m-3, 293 k, 383.1 m2 k-1s-2, 3.86 1010 kg m-1 s-2, 7.76×1010 kg m-1 s-2.
Thus, we get the non-dimensional values of the problem as [8, 24]: 0.01047, 0.0419, 1.6086, 4, 1.0, 0.02, 0.01, 0.1 (assumed).
The numerical results of the conductive and dynamical temperature increments, the strain, the displacement, and the average stress distributions have been figured for a wide range of the dimensionless radial distance at the instant value of dimensionless time 1.0.
In Figs. 2-6, the calculations have been carried out for various values of the fractional-order parameter in the context of the hyperbolic two-temperature thermoelasticity model to stand on the effect of this parameter on all the studied functions.
It is noted from Figs. 2 and 3 that the fractional-order parameter has very limited effects on the conductive and dynamical temperature increment. It is expectable results where the fractional-order definition has been applied to the strain function.
Fig. 2. The conductive temperature distribution with the variance value of the fractional-order parameter
Fig. 3. The dynamical temperature distribution with variance value of fractional order parameter
Fig. 4 shows that the fractional-order parameter has a significant effect on the strain distributions. The absolute values of the peak points of the strain distribution take the following order:
Fig. 5 shows that the fractional-order parameter has a significant effect on the displacement distributions, and this effect is visible near the impact of the thermal shock. The absolute values of the displacement at the edge of the cylinder take the following order:
Fig. 6 shows that the fractional-order parameter has a significant effect on the average stress distributions. The absolute values of the average stress at the center of the cylinder take the following order:
Fig. 4. The strain distribution with variance value of fractional order parameter
Fig. 5. The displacement distribution with variance value of fractional order parameter
Fig. 6. The average stresses distribution with variance value of fractional order parameter
In Figs. 7-11, the calculations have been carried out for various models of thermoelasticity (one-temperature, two-temperature, and hyperbolic two-temperature) without fractional-order strain to stand on the two-temperature parameter effects on all the studied functions.
In Fig. 7, it is noticed that the behavior of the conductive temperature increment distributions in the context of the one-temperature and hyperbolic two-temperature models is closed but with different values. In contrast, the behavior of the classical two-temperature model is different. The curves of the one-temperature and hyperbolic two-temperature models fail down close to zero, more than the curve of the two-temperature model. The one-temperature and hyperbolic two-temperature models offer thermal waves with finite speed of propagation while the model of classical two-temperature parameter offers thermal waves with an infinite speed of propagation.
Fig. 8 shows that the dynamical temperature increment distributions in the context of the one-temperature and hyperbolic two-temperature models have the same behavior but with different values, while the behavior of the classical two-temperature model is different. The values of the dynamical temperature increment based on the one-temperature model in Fig. 8 coincide with the value of the conductive temperature increment in Fig. 7 based on the same model.
Fig. 7. The conductive temperature increment distribution with variance models and
Fig. 8. The dynamical temperature increment distribution with variance models
Fig. 9 shows that the strain distributions in the context of the one-temperature and hyperbolic two-temperature models have the same behavior but with different values, while its behavior based on the classical two-temperature model is different. The absolute values of the peak points of the strain take the following order:
Fig. 9. The strain increment distribution with variance models
Fig. 10 shows that the displacement distributions in the context of the one-temperature and hyperbolic two-temperature models have closed behavior, while its behavior based on the classical two-temperature model is different. The values of the peak points of the displacement take the following order:
Fig. 11 shows that the average stress distributions in the context of the one-temperature and hyperbolic two-temperature models have a closed behavior, while its behavior based on the classical two-temperature model is different. The absolute values of the stress in the edge of the cylinder take the following order:
We can see that for the special cases of zero fractional order of strain, one-temperature model, and classical two-temperature model, the results agree with the results in the publications [8, 27, 28].
Fig. 10. The displacement distribution with variance models
Fig. 11. The average stresses distribution with variance models
The issue of a thermoelastic homogeneous isotropic infinite cylindrical medium when the bounding surface is thermally shocked and no volume deformation is addressed in this article. One-temperature, classical two-temperature, and hyperbolic two-temperature models have also been used to solve the problem.
The fractional-order parameter has negligible impacts on conductive and dynamical temperature increments, but has a substantial influence on the strain, displacement, and average stresses, according to the numerical data.
The one-temperature model, like the hyperbolic two-temperature model, produces thermal waves with finite propagation speeds, while the classical two-temperature model generates thermal waves with an infinite propagation speed.
As a result, the hyperbolic two-temperature thermoelasticity model is a good fit for describing the thermodynamic behavior of thermoelastic materials.
- P. J. Chen and M. E. Gurtin, “On a theory of heat conduction involving two temperatures,” Zeitschrift für angewandte Mathematik und Physik ZAMP, Vol. 19, No. 4, pp. 614–627, Jul. 1968, https://doi.org/10.1007/bf01594969 [Publisher]
- W. E. Warren and P. J. Chen, “Wave propagation in the two temperature theory of thermoelasticity,” Acta Mechanica, Vol. 16, No. 1-2, pp. 21–33, Mar. 1973, https://doi.org/10.1007/bf01177123 [Publisher]
- H. M. Youssef, “Theory of two-temperature-generalized thermoelasticity,” IMA Journal of Applied Mathematics, Vol. 71, No. 3, pp. 383–390, Jun. 2006, https://doi.org/10.1093/imamat/hxh101 [Publisher]
- H. M. Youssef and K. A. Elsibai, “On the theory of two-temperature thermoelasticity without energy dissipation of Green-Naghdi model,” Applicable Analysis, Vol. 94, No. 10, pp. 1997–2010, Oct. 2015, https://doi.org/10.1080/00036811.2014.961920 [Publisher]
- I. A. Abbas and H. M. Youssef, “Two-temperature generalized thermoelasticity under ramp-type heating by finite element method,” Meccanica, Vol. 48, No. 2, pp. 331–339, Mar. 2013, https://doi.org/10.1007/s11012-012-9604-8 [Publisher]
- H. Youssef, “A two-temperature generalized thermoelastic medium subjected to a moving heat source and ramp-type heating: a state-space approach,” Journal of Mechanics of Materials and Structures, Vol. 4, No. 9, pp. 1637–1649, Jan. 2010, https://doi.org/10.2140/jomms.2009.4.1637 [Publisher]
- H. M. Youssef, “Two-dimensional problem of a two-temperature generalized thermoelastic half-space subjected to ramp-type heating,” Computational Mathematics and Modeling, Vol. 19, No. 2, pp. 201–216, Apr. 2008, https://doi.org/10.1007/s10598-008-0014-7 [Publisher]
- H. M. Youssef, “Two-temperature generalized thermoelastic infinite medium with cylindrical cavity subjected to moving heat source,” Archive of Applied Mechanics, Vol. 80, No. 11, pp. 1213–1224, Nov. 2010, https://doi.org/10.1007/s00419-009-0359-1 [Publisher]
- H. M. Youssef and A. A. El-Bary, “Theory of hyperbolic two-temperature generalized thermoelasticity,” Materials Physics and Mechanics, Vol. 40, No. 2, pp. 158–171, Jan. 2018, https://doi.org/10.18720/mpm.4022018_4 [Search CrossRef]
- R. L. Magin and T. J. Royston, “Fractional-order elastic models of cartilage: A multi-scale approach,” Communications in Nonlinear Science and Numerical Simulation, Vol. 15, No. 3, pp. 657–664, Mar. 2010, https://doi.org/10.1016/j.cnsns.2009.05.008 [Publisher]
- C. Li, H. Guo, X. Tian, and T. He, “Generalized thermoviscoelastic analysis with fractional order strain in a thick viscoelastic plate of infinite extent,” Journal of Thermal Stresses, Vol. 42, No. 8, pp. 1051–1070, Aug. 2019, https://doi.org/10.1080/01495739.2019.1587331 [Publisher]
- C. Li, H. Guo, X. Tian, and T. He, “Generalized thermoelastic diffusion problems with fractional order strain,” European Journal of Mechanics – A/Solids, Vol. 78, p. 103827, Nov. 2019, https://doi.org/10.1016/j.euromechsol.2019.103827 [Publisher]
- C. Li, H. Guo, X. Tian, and T. He, “Generalized piezoelectric thermoelasticity problems with strain rate and transient thermo-electromechanical responses analysis,” ZAMM – Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik, Vol. 100, No. 5, p. e201900067, May 2020, https://doi.org/10.1002/zamm.201900067 [Publisher]
- E. M. Hussein, “Fractional order thermoelastic problem for an infinitely long solid circular cylinder,” Journal of Thermal Stresses, Vol. 38, No. 2, pp. 133–145, Feb. 2015, https://doi.org/10.1080/01495739.2014.936253 [Publisher]
- S. Sheoran and P. Kundu, “Fractional order generalized thermoelasticity theories: A review,” International Journal of Advances in Applied Mathematics and Mechanics, Vol. 3, No. 4, pp. 76–81, 2016. [Search CrossRef]
- M. A. Ezzat, A. S. El-Karamany, and A. A. El-Bary, “Two-temperature theory in Green–Naghdi thermoelasticity with fractional phase-lag heat transfer,” Microsystem Technologies, Vol. 24, No. 2, pp. 951–961, Feb. 2018, https://doi.org/10.1007/s00542-017-3425-6 [Publisher]
- M. A. Ezzat and A. A. El-Bary, “Application of fractional order theory of magneto-thermoelasticity to an infinite perfect conducting body with a cylindrical cavity,” Microsystem Technologies, Vol. 23, No. 7, pp. 2447–2458, Jul. 2017, https://doi.org/10.1007/s00542-016-2976-2 [Publisher]
- N. Sarkar, D. Ghosh, and A. Lahiri, “A two-dimensional magneto-thermoelastic problem based on a new two-temperature generalized thermoelasticity model with memory-dependent derivative,” Mechanics of Advanced Materials and Structures, Vol. 26, No. 11, pp. 957–966, Jun. 2019, https://doi.org/10.1080/15376494.2018.1432784 [Publisher]
- S. Choudhary, S. Kumar, and J. S. Sikka, “Thermo-mechanical interactions in a fractional order generalized thermoelastic solid with diffusion,” Microsystem Technologies, Vol. 23, No. 12, pp. 5435–5446, Dec. 2017, https://doi.org/10.1007/s00542-017-3340-x [Publisher]
- M. A. Ezzat, A. S. El-Karamany, and A. A. El-Bary, “Application of fractional order theory of thermoelasticity to 3D time-dependent thermal shock problem for a half-space,” Mechanics of Advanced Materials and Structures, Vol. 24, No. 1, pp. 27–35, Jan. 2017, https://doi.org/10.1080/15376494.2015.1091532 [Publisher]
- M. H. Hendy, M. M. Amin, and M. A. Ezzat, “Application of fractional order theory to a functionally graded perfect conducting thermoelastic half space with variable Lamé’s Modulii,” Microsystem Technologies, Vol. 23, No. 10, pp. 4891–4902, Oct. 2017, https://doi.org/10.1007/s00542-017-3409-6 [Publisher]
- M. A. Ezzat and A. A. El-Bary, “Fractional order theory to an infinite thermo-viscoelastic body with a cylindrical cavity in the presence of an axial uniform magnetic field,” Journal of Electromagnetic Waves and Applications, Vol. 31, No. 5, pp. 495–513, Mar. 2017, https://doi.org/10.1080/09205071.2017.1285728 [Publisher]
- H. M. Youssef, “Theory of generalized thermoelasticity with fractional order strain,” Journal of Vibration and Control, Vol. 22, No. 18, pp. 3840–3857, Oct. 2016, https://doi.org/10.1177/1077546314566837 [Publisher]
- H. M. Youssef, “Generalized thermoelastic infinite medium with cylindrical cavity subjected to moving heat source,” Mechanics Research Communications, Vol. 36, No. 4, pp. 487–496, Jun. 2009, https://doi.org/10.1016/j.mechrescom.2008.12.004 [Publisher]
- Y. Povstenko, Fractional Thermoelasticity. Springer, 2015. [Publisher]
- D. Y. Tzou, “A unified field approach for heat conduction from macro – to micro-scales,” Journal of Heat Transfer, Vol. 117, No. 1, pp. 8–16, Feb. 1995, https://doi.org/10.1115/1.2822329 [Publisher]
- Hm Youssef, “Two-temperature generalized thermoelastic infinite medium with cylindrical cavity subjected to different types of thermal loading,” WSEAS Transactions on Heat and Mass Transfer, Vol. 1, No. 10, p. 769, Jan. 2006. [Search CrossRef]
- H. M. Youssef, “Two-temperature generalized thermoelastic infinite medium with cylindrical cavity subjected to non-gaussian laser beam,” Two-Temperature Generalized Thermoelastic Infinite Medium with Cylindrical Cavity Subjected to Non-Gaussian Laser Beam, Vol. 21, No. 2, pp. 222–229, 2016, https://doi.org/10.20855/ijav.2016.21.2415 [Publisher]