Next Article in Journal
Towards Ethical Engineering: Artificial Intelligence as an Ethical Governance Tool for Emerging Technologies
Previous Article in Journal
Research on the Uncertainty of Musical Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Fractional Order Thermoelastic Model with Voids in Three-Phase-Lag Thermoelasticity †

Department of Mathematics, Jadavpur University, Kolkata 700032, India
*
Author to whom correspondence should be addressed.
Presented at the 1st International Online Conference on Mathematics and Applications, 1–15 May 2023; Available online: https://iocma2023.sciforum.net/.
These authors contributed equally to this work.
Comput. Sci. Math. Forum 2023, 7(1), 57; https://doi.org/10.3390/IOCMA2023-14430
Published: 28 April 2023

Abstract

:
This paper deals with fractional order three-phase-lag (TPL) thermo-elasticity in a micropolar thermoelastic half space medium with voids. The subsequent non-dimensional coupled equations are solved by using the normal mode analysis and eigenvalue approach methods. By doing numerical computations of the physical fields for a substance that resembles a magnesium crystal in the presence of an electromagnetic field, the issue is proven to exist. The effect of the fractional order, the phase lags on the components of temperature, displacement, the stress, and changes in volume fraction field have been depicted graphically. Additionally, a graphic comparison of several types of models employing phase delays and the influence of the magnetic field is displayed.

1. Introduction

Lord and Shulman (1967) [1] derived a new model to overcome the infinite speed of heat transportation which was established by Biot (1956) [2]. The heat equation of this theory is of the wave type, it automatically ensures finite speeds of propagation for heat and elastic waves. Linear model of thermoelasticity is incapable to characterize the behaviour of many new synthetic solids of the elastomer and polymer type. Eringen used the micropolar elasticity term to describe such materials.
Fractional calculus was first used by Abel in the solution of an integral equation that arises in the tautochrone problem. This area has grown rapidly, and applications have been found in several fields, including solid mechanics, geophysics, physics, and mathematical biology. Many mathematical models in the fields of solid mechanics, bio-rheology, non-linear dynamical systems in ecology, and so on, have been successfully modified using fractional calculus. The fractional derivative exhibits non-local properties, and global dependency is among the main reason for its use. Kimmich considered anomalous diffusion and characterised it with the time-fractional diffusion wave equation using the Riemann–Liouville fractional integral. Povstenko [3,4] demonstrated the effect of fractional heat transportation in the presence of thermal stresses. Many authors [5,6,7,8,9] have also discussed fractional calculus in thermoelasticity.
Recently, Tzou [10,11] developed dual-phase-lag heat equation theory by incorporating of two-phase-lags associated with temperature gradient heat flux vector. Two phase lag parameters are related to the fastest effects of thermal inertia. This theory is known as the dual-phase-lag (DPL) model. The later three-phase-lag (TPL) model was developed by Roychoudhury [12].
In the present article, we consider two dimensional problem under TPL heat conduction theory in present of electric and magnetic field. The governing equation are transformed by normal modes. Finally the displacement component and temperature distribution are fixed out analytically from the vector matrix differential equation by eigenvalue method in transformed domain.

2. Basic Equation

The system of governing equations of a rotating micropolar thermoelastic solid is given by [12,13]
σ i j , j + F i = ρ [ u i ¨ + Ω × ( Ω × u ) i + ( 2 Ω × u ˙ ) i ]
m i j , i + ϵ i j r σ i r = J ρ [ ϕ ¨ + ( Ω × ϕ ) i ]
α ψ , i i ε ψ ω ψ ˙ β * u i , i + m θ = ρ ε 1 ψ ¨
The components of the force stress and couple stress tensors are
σ i j = λ u r , r δ i j + μ ( u i , j + u j , i ) + k ( u j , i ϵ i j r ϕ r ) β * ψ δ i j γ T ( θ + ν θ t ) δ i j
m i j = α ϕ r , r δ i j + β ϕ i , j + γ ϕ j , i
Fractional order Three-phase-lag Heat conduction without heat source as
K * 1 + τ ν α α ! α t α + K t 1 + τ θ α α ! α t α θ i i = 1 + τ α q α ! α t α + τ 2 α q 2 α ! 2 α t 2 α ( ρ C E 2 θ t 2 + β T 0 2 t 2 ( e k k ) + m T 0 2 ψ t 2 )
and the strain components are
e i j = 1 2 ( u i , j + u j i ) and e k k = ( u x + w z )

3. Formulation of the Problem

( 1 + ϵ 5 ) 2 u x 2 + ( 1 a 2 + ϵ 5 ) 2 w x z + a 2 2 u z 2 ϵ 3 ϕ 2 z 1 ( ϵ 1 + ϵ 2 ) ( 1 + v t ) θ x + β 1 ψ x = ( ϵ 4 + 1 ( ϵ 1 + ϵ 2 ) ) 2 u t 2
( 1 + ϵ 5 ) 2 w z 2 + ( 1 a 2 + ϵ 5 ) 2 u x z + a 2 2 w x 2 ϵ 3 ϕ 2 x 1 ( ϵ 1 + ϵ 2 ) ( 1 + v t ) θ z + β 1 ψ z = ( ϵ 4 + 1 ( ϵ 1 + ϵ 2 ) ) 2 w t 2
2 ϕ 2 x 2 + 2 ϕ 2 z 2 2 ϵ 7 ϕ 2 + ϵ 7 ( u z w x ) = ϵ 8 2 ϕ 2 t 2
K 1 ( 1 + τ v α α ! α t α ) + t ( 1 + τ θ α α ! α t α ) ( 2 θ x 2 + 2 θ z 2 ) = p l 2 θ t 2 + ϵ 9 2 t 2 ( u x + w z ) + β 1 2 ψ t 2
The components of stress tensor and force tensor are written in the below form:
σ x x = ( ϵ 1 + ϵ 2 ) u x + ( ϵ 1 1 ) w x ( 1 + ν t ) θ
σ z z = ( ϵ 1 + ϵ 2 ) w z + ( ϵ 1 1 ) u x ( 1 + ν t ) θ
σ x z = u z + ϵ 2 w x ( ϵ 2 1 ) ϕ 2
σ z x = w x + ϵ 2 u z ( ϵ 2 1 ) ϕ 2
m z y = ( ϵ 2 1 ) ϵ 4 ϕ 2 z
m x y = ( ϵ 2 1 ) ϵ 4 ϕ 2 x
where ϵ i , i = 1 ( 1 ) 9 are given in Appendix A.
To simplify the Equations (1)–(7), we use the non-dimensional variables as below
( x * , z * ) = η 0 c 0 ( x , z ) , ( t * , t 0 * , τ 0 * , τ q * , τ T * , υ * ) = η 0 ( t , t 0 , τ 0 , τ q , τ T , υ ) , Ω * = Ω η 0 ( u * , w * ) = ρ c 0 η 0 γ T T 0 ( u , w ) , θ * = θ T 0 , σ i j * = σ i j γ T T 0 , m i j * = η 0 m i j c 0 γ T T 0 , ( ψ * , ϕ 2 * ) = ρ c 0 2 γ T T 0 ( ψ , ϕ 2 ) .
where
c 0 2 = μ ρ , β = λ + 2 μ μ , η 0 = ρ c E c 0 2 k , p l = 1 + τ α q α ! α t α + τ 2 α q 2 α ! 2 α t 2 α .
Equations (8)–(17) takes the above form, dropping the * for convenience.

4. Normal Mode Analysis

The decomposition of the solution of the physical variables under consideration has the following form in the normal mode
( u , w , ψ , θ , ϕ 2 , σ j k ) ( x , z , t ) = ( u ¯ , w ¯ , ψ ¯ , θ ¯ , ϕ 2 ¯ , σ j k ¯ ) ( x ) e ( s t + i q z ) ;
where u ¯ , w ¯ , θ ¯ , ϕ 2 ¯ , σ j k ¯ , Q ¯ are the amplitudes of the functions, s is a complex constant, i = 1 and q is the wave numbers in the z-direction.
Using above equations we obtain the vector matrix differential equation as
d z ¯ d x = A z ¯
where
z ¯ = ( u ¯ , w ¯ , ϕ 2 ¯ , θ ¯ , ψ ¯ , d u ¯ d x , d w ¯ d x , d ϕ 2 ¯ d x , d θ ¯ d x , d ψ ¯ d x ) T
A = 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 a 51 a 52 a 53 0 0 0 a 56 0 a 58 a 59 a 61 a 62 0 a 64 0 a 65 0 a 67 0 0 a 71 0 a 73 0 0 0 a 76 a 77 0 0 0 a 82 0 a 84 0 a 85 a 86 a 87 a 88 0 0 a 92 0 a 94 a 95 a 96 0 0 0 0
and p l s = 1 + τ q α α ! s α + τ q 2 α 2 α ! s 2 α K 1 ( 1 + τ v α α ! s α ) + s ( 1 + τ θ α α ! s α ) and all the components of matrix A in Equation (22) are given in Appendix A.

5. Solution

The following equation represents the matrix A’s characteristic equation:
det ( A Λ I ) = 0
The eigenvalues of the characteristic Equation (20) are taken in the form Λ = ± Λ i ( i = 1 , 2 , 3 , 4 , 5 )
The eigen vector X corresponding to the eigen value Λ calculated as
X = Γ 1 Γ 2 Γ 3 Γ 4 Γ 5 Λ Γ 1 Λ Γ 2 Λ Γ 3 Λ Γ 4 Λ Γ 5 T
where Γ i ( i , j = 1 , 2 , 3 , 4 , 5 ) are given in Appendix A. We construct the inverse of the matrix V = X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , X 7 , X 8 , X 9 , X 10 = ( X i j ) 10 × 10 , i , j = 1 ( 1 ) 10 as V 1 = ω i j , i , j = 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 10 .
Then the solution of the differential Equation (23) is [14,15]
z ¯ = j = 1 5 X j y j
y r = C r e λ r x + e λ r x Q r e λ r x a n d   Q r = j = 1 10 ω r j f j
where C r is an arbitrary constant which is to be evaluated using initial and boundary conditions.
u = e ( s t + i q y ) j = 1 5 X 1 j ( C j e Λ j x Q j Λ j )
w = e ( s t + i q y ) j = 1 5 X 2 j ( C j e Λ j x Q j Λ j )
θ = e ( ω t + i b y ) j = 1 5 X 3 j ( C j e Λ j x Q j Λ j )
Using Equations (27)–(29) simplifying the above equations we obtain the stress components as follows
σ z z = e ( ω t + i b z ) j = 1 5 C j S 1 j ( x ) e Λ j x M 1 , σ z x = e ( ω t + i b z ) j = 1 5 C j S 2 j ( x ) e Λ j x M 2
σ x x = e ( ω t + i b z ) j = 1 5 C j S 4 j ( x ) e Λ j x M 4 , σ x z = e ( ω t + i b z ) j = 1 5 C j S 3 j ( x ) e Λ j x M 3
m x y = e ( ω t + i b z ) j = 1 5 C j S 6 j e Λ j x M 6 , m z y = e ( ω t + i b z ) j = 1 5 C j S 5 j e Λ j x M 5
where C j , S i j (x) and M k (i = 1, 2, 3, 4, 5, 6, 7; j = 1, 2, 3, 4; k = 1, 2, 3, 4, 5, 6, 7) are given in the Appendix A.

6. Discussion

We discuss TPL model in Figure 1, Figure 2, Figure 3 and Figure 4. In Figure 1 the fractional order parameter α has a significant effect on the temperature distribution, where increasing on α causes increasing on and the rate of change of with respect to x also increases when α increases at fixed t = 0.1 .
In Figure 2, when fractional order α = 1 , time t = 0.1 are fixed and other constants are same, then u displacement distribution starts with negative values for TPL, LS, CT theory. The u-distribution up to x = 1.2 for the LS model is greater than for the TPL model and CT model. After x = 1.2 the result is reversed.
In Figure 3 and Figure 4, the fractional parameter has significant effects on the stress σ x x and σ x z distributions. Both starts from zero initially which shows that they satisfy boundary condition. Increasing of fractional parameter α causes decreasing the absolute values of the stresses, and the rate of change of them with respect to x also increase when α increases. For the fixed value of x the stress σ x z distribution has the lowest value for the fractional parameter α = 0.2 in the range 0 x 1.0 . For x 1.0 the distribution gradually increases towards zero.

7. Conclusions

In this work, the effect of the heat source, fractional order and phase lag parameters on the temperature distribution, displacement components, the stress components have been studying for a two-dimensional problem in a half space micro-polar thermoelastic material is considered with the context of the fractional order thermoelasticity theory caring a heat source. We found that, the fractional order parameter has significant effects on all the studied fields and the results supporting the definition of the classification of the thermal conductivity of the materials. We also found significant effects of phase lag parameters.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data took places in this paper.

Acknowledgments

Sourov Roy would like to acknowledge the support of Council of Scientific and Industrial Research, New Delhi, India (Fileno.-09/096(0949)/2018-EMR-I).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
TPLThree phase lag
LSLord–Shulman
CTCoupled Thermoelasticity theory

Appendix A

a 51 = ϵ 4 s 2 + a 2 q 2 + s 2 ϵ 1 + ϵ 2 1 + ϵ 5 , a 52 = i q ϵ 3 1 + ϵ 5 , a 53 = i q ( a 2 1 ) i q ϵ 5 1 + ϵ 5 , a 56 = ( 1 + ν s ) ( ϵ 1 + ϵ 2 ) ( 1 + ϵ 5 ) , a 58 = b 3
a 61 = q 2 + ϵ 4 s 2 + ϵ 5 q 2 + s 2 ϵ 1 + ϵ 2 a 2 , a 62 = i q ( 1 + ν s ) a 2 ( ϵ 1 + ϵ 2 ) , a 64 = i q b 3 a 2 , a 65 = i q ( a 2 1 ) ϵ 5 a 2 , a 67 = ϵ 3 a 2
a 71 = i q ϵ 7 , a 73 = q 2 + ϵ 8 s 2 + 2 ϵ 7 , a 77 = ϵ 7
a 82 = i q p l s ϵ 9 , a 84 = ( p l s + q 2 ) , a 85 = p l s ϵ 10 , a 86 = p l s ϵ 9 s ,
a 92 = i q ϵ 11 , a 93 = 0 , a 94 = ϵ 13 , a 95 = q 2 + ϵ 14 s 2 + s ϵ 6 + ϵ 12 , a 96 = ϵ 11
ϵ 1 = λ + μ ρ c 0 2 , ϵ 2 = μ + k ρ c 0 2 , ϵ 3 = k / ρ c 0 2 ( ϵ 1 + ϵ 2 ) , ϵ 4 = ϵ 0 μ 0 2 H 0 2 ρ ( ϵ 1 + ϵ 2 ) , ϵ 5 = μ 0 2 H 0 2 ρ c 0 2 ( ϵ 1 + ϵ 2 ) , ϵ 7 = k c 0 2 γ η 0 2 ,
ϵ 8 = ρ J η 0 2 γ , ϵ 9 = β T 0 γ T ρ K η 0 , ϵ 10 = m T 0 γ T ρ K η 0
ϵ 6 = ω c 0 2 β η 0 , ϵ 11 = b c 0 2 β η 0 2 ϵ 12 = ϵ c 0 2 β η 0 ϵ 13 = m ρ c 0 4 β γ T η 0 2 ϵ 14 = ρ ξ c 0 2 β

References

  1. Lord, H.W.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  2. Biot, M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range. J. Acoust. Soc. Am. 1956, 28, 168–178. [Google Scholar]
  3. Povstenko, Y.Z. Fractional heat conduction and associated thermal stresses. J. Therm. Stress. 2004, 28, 83–102. [Google Scholar]
  4. Povstenko, Y.Z. Fractional Catteneo-type equations and generalized thermoelasticity. J. Therm. Stress. 2011, 34, 94–114. [Google Scholar]
  5. Youssef, H.M.; Al-Lehaibi, E.A. Fractional order generalized thermoelastic infinite medium with cylindrical cavity subjected to harmonically varying heat. Sci. Res. 2010, 3, 32–37. [Google Scholar]
  6. Sherief, H.H.; El-Sayed, A.M.A.; El-Latief, A.M.A. Fractional order theory of thermoelasticity. Int. J. Solids Struct. 2010, 47, 269–275. [Google Scholar] [CrossRef]
  7. Roy, S.; Lahiri, A. A Study on Fractional Order Thermoelastic Half Space. Int. J. Appl. Mech. Eng. 2020, 25, 191–202. [Google Scholar] [CrossRef]
  8. Abbas, I.; Youssef, H. Two-Dimensional Fractional Order Generalized Thermoelastic Porous Material. Lat. Am. J. Solids Struct. 2015, 12, 1415–1431. [Google Scholar] [CrossRef]
  9. Peng, W.; Ma, Y.; Li, C.; He, T. Dynamic analysis to the fractional order thermoelastic diffusion problem of an infinite body with a spherical cavity and variable material properties. J. Therm. Stress. 2020, 43, 38–54. [Google Scholar] [CrossRef]
  10. Tzou, D.Y. A unified field approach for heat conduction from micro to macroscale. J. Heat Transfer. 1995, 117, 8–16. [Google Scholar]
  11. Tzou, D.Y. Macro to Microscale Heat Transfer: The Lagging Behavior; Taylor & Francis: Washington, DC, USA, 1996. [Google Scholar]
  12. Roychoudhuri, S.K. On a thermoelastic three-phase-lag model. J. Therm. Stress. 2007, 30, 231–238. [Google Scholar]
  13. Abouelregal, A.E. Modified Fractional Photo-Thermoelastic Model for a Rotating Semiconductor Half-Space Subjected to a Magnetic Field. Silicon 2020, 12, 2837–2850. [Google Scholar] [CrossRef]
  14. Das, N.C.; Lahiri, A.; Sarkar, S. Eigenvalue value approach three dimensional coupled thermoelasticity in a rotating transversely isotropic medium. Tamsui Oxf. J. Mathimatical Sci. 2009, 25, 237–257. [Google Scholar]
  15. Othman, M.I.A.; Abbas, I.A. Eigenvalue approach for generalized thermoelastic porous medium under the effect of thermal loading due to a laser pulse in DPL model. Indian J. Phys. 2019, 93, 1567–1578. [Google Scholar] [CrossRef]
Figure 1. Impact of fractional order parameter on temperature distribution at t = 0.1 and z = 1.
Figure 1. Impact of fractional order parameter on temperature distribution at t = 0.1 and z = 1.
Csmf 07 00057 g001
Figure 2. Impact of different thermoelasticity theory on u distribution at t = 0.1 and z = 1.
Figure 2. Impact of different thermoelasticity theory on u distribution at t = 0.1 and z = 1.
Csmf 07 00057 g002
Figure 3. Impact of fractional order parameter on σ x x distribution at t = 0.1 and z = 1.
Figure 3. Impact of fractional order parameter on σ x x distribution at t = 0.1 and z = 1.
Csmf 07 00057 g003
Figure 4. Impact of fractional order parameter on σ x z distribution at t = 0.1 and z = 1.
Figure 4. Impact of fractional order parameter on σ x z distribution at t = 0.1 and z = 1.
Csmf 07 00057 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Roy, S.; Lahiri, A. Fractional Order Thermoelastic Model with Voids in Three-Phase-Lag Thermoelasticity. Comput. Sci. Math. Forum 2023, 7, 57. https://doi.org/10.3390/IOCMA2023-14430

AMA Style

Roy S, Lahiri A. Fractional Order Thermoelastic Model with Voids in Three-Phase-Lag Thermoelasticity. Computer Sciences & Mathematics Forum. 2023; 7(1):57. https://doi.org/10.3390/IOCMA2023-14430

Chicago/Turabian Style

Roy, Sourov, and Abhijit Lahiri. 2023. "Fractional Order Thermoelastic Model with Voids in Three-Phase-Lag Thermoelasticity" Computer Sciences & Mathematics Forum 7, no. 1: 57. https://doi.org/10.3390/IOCMA2023-14430

APA Style

Roy, S., & Lahiri, A. (2023). Fractional Order Thermoelastic Model with Voids in Three-Phase-Lag Thermoelasticity. Computer Sciences & Mathematics Forum, 7(1), 57. https://doi.org/10.3390/IOCMA2023-14430

Article Metrics

Back to TopTop