Next Article in Journal
A Novel Fractional-Order Non-Convex TVα,p Model in Image Deblurring
Next Article in Special Issue
Enhancing Power Quality in Decentralized Hybrid Microgrids: Optimized DSTATCOM Performance Using Cascaded Fractional-Order Controllers and Hybrid Optimization Algorithms
Previous Article in Journal
Theoretical Analysis of Viscoelastic Friction System Characteristics of Robotic Arm Brake Based on Fractional Differential Theory
Previous Article in Special Issue
Leveraging the Performance of Integrated Power Systems with Wind Uncertainty Using Fractional Computing-Based Hybrid Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

State-Space Approach to the Time-Fractional Maxwell’s Equations under Caputo Fractional Derivative of an Electromagnetic Half-Space under Four Different Thermoelastic Theorems

by
Eman A. N. Al-Lehaibi
1 and
Hamdy M. Youssef
2,*
1
Mathematics Department, Jamoum University College, Umm Al-Qura University, Jamoum 25375, Saudi Arabia
2
Mechanical Engineering Department, College of Engineering and Architecture, Umm Al-Qura University, Makkah 21955, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(10), 566; https://doi.org/10.3390/fractalfract8100566
Submission received: 23 August 2024 / Revised: 19 September 2024 / Accepted: 24 September 2024 / Published: 28 September 2024

Abstract

:
This paper introduces a new mathematical modelling method of a thermoelastic and electromagnetic half-space in the context of four different thermoelastic theorems: Green–Naghdi type-I, and type-III; Lord–Shulman; and Moore–Gibson–Thompson. The bunding plane of the half-space surface is subjected to ramp-type heat and traction-free. We consider that Maxwell’s time-fractional equations have been under Caputo’s fractional derivative definition, which is the novelty of this work. Laplace transform techniques are utilized to obtain solutions using the state-space approach. Laplace transform’s inversions were calculated using Tzou’s iteration method. The temperature increment, strain, displacement, stress, induced electric field, and induced magnetic field distributions were obtained numerically and are illustrated in figures. The time-fraction parameter of Maxwell’s equations had a major impact on all the studied functions. The time-fractional parameter of Maxwell’s equations worked as resistant to the changing of temperature, particle movement, and induced magnetic field, while it acted as a catalyst to the induced electric field through the material. Moreover, all the studied functions have different values in the context of the four studied theorems.

1. Introduction

The link between strain and magnetic fields in thermoelastic substances is a subject of increasing interest due to its many important applications in geophysics, plasma physics, and related disciplines. Temperature gradients, very high temperatures, and internal magnetic fields of nuclear reactors all affect the field’s comprehension of how these devices operate. This is the domain of magneto-thermoelasticity. This hybrid theory combines the domains of thermoelasticity and electromagnetism [1]. As a result, the governing system of equations of the electromagnetic theory consists of hyperbolic partial differential equations, which ensure finite wave propagation rates. Biot’s theory of linked thermoelasticity is based on the parabolic partial differential equation of heat conduction and the hyperbolic partial differential equation of motion [2]. The features of the second equation suggest that heat waves may move at infinite speeds, despite evidence to the contrary. Lord and Shulman developed generalized thermoelasticity, with one relaxation time theory to overcome this restriction, and suggested a new and distinct rule of heat transfer to replace the traditional Fourier’s law of heat conduction (L-S) [3]. Both the heat flow vector and its time derivative are included in this law. Moreover, it has a new parameter, the relaxation time. There are limited propagation speeds for both heat and elastic waves since this theory’s heat equation is of the wave type. Like coupled and uncoupled theories, the other governing differential equations of this theory are the differential equations of motion and the constitutive relations of stress and strain. This theory was expanded upon by Dhaliwal and Sherief to include the general anisotropic case [4]. Green and Naghdi (G-N) presented three theories for investigating thermoelasticity. Instead of inequality, these models are based on entropy equality. Furthermore, the heat flow vector’s focus shifts [5,6,7]. The Green–Naghdi theorem is divided into three types: type I, type II, and type III, which represent three distinct assumptions. Type I uses the same traditional thermoelasticity mechanism and has a linear form. The waste of energy is prohibited under Type II. Type II is a subset of type III. Type III allows for the waste of energy [6,7,8,9,10,11,12]. In recent years, there has been significant interest in the Moore–Gibson–Thompson (MGT) heat conduction equation, and several articles have been published on it [13,14,15,16].
Whether one is discussing coupled or uncoupled versions, the generalized heat equation is often disregarded while discussing magneto-thermoelastic problems. This method is often justified since the outcomes of solving any of these equations vary very little quantitatively. However, when the short-term effects are considered, a significant degree of accuracy is lost unless the whole generalized system of differential equations is employed. This concept comes in quite handy for many problems involving heat gradients and for searching for fast-acting effects [1,17,18].
Leibniz presented the derivative of order half for the first time, initiating the extensive history of fractional calculus [19,20]. Leibniz, Liouville, Grunwald, Letnikov, and Riemann are credited with establishing the theory of fractional-order derivatives and integrals. Regarding fractional calculus and fractional differential equations with solutions, there are many fascinating sources and books available [19,21]. Recent physics research has used fractional-order derivatives and integrals, as well as fractional integro-differential equations, for a variety of purposes [22,23]. In scholarly circles, fractional-order electrodynamics is a relatively young subject. Among the many techniques, the addition of fractional-order derivatives to Maxwell’s equations in both space and time stands out [24,25,26,27,28]. These expansions of Maxwell’s equations, which make use of fractional-order derivatives, may be used to characterize the dynamics of electromagnetic systems with memory and energy dissipation [28]. The electric potential may be expressed in terms of fractional-order poles, which is an additional advantage of using fractional-order derivatives in electrodynamics [28,29]. Youssef developed the theory of fractional-order generalized thermoelasticity, which is based on the fractional-order non-Fourier heat conduction law [30] and the theory of generalized thermoelasticity with fractional-order strain based on fractional-order equations of motion, which is related to fractional-order stress–strain relations [31]. Neslihan examined the chaotic and regular dynamics of classical and fractional Gross–Pitaevskii equations (GPE) for interacting boson systems subjected to combined harmonic and optical lattice potentials through Poincaré sections of phase space, Lyapunov exponents, power spectrum analysis, and bifurcation techniques. They also studied chaotic and regular behaviours of classical and fractional Gross–Pitaevskii equations, including two-body, three-body, and higher-order interactions [32]. Al-Raeei simulated the spatial form of the fractional Schrödinger equation for the electrical screening potential using the Riemann–Liouville definition of the fractional derivatives and the numerical simulation methods [33].
State space approaches provide the fundamental basis of contemporary control theory, using differential equations rather than transfer functions to model processes. This technique is more encompassing than traditional Laplace and Fourier transform theory, allowing it to be applied to any system that can be examined using integral transforms in time. Additionally, it offers a distinct perspective on the temporal dynamics of linear systems. The state space approach is advantageous due to its ability to analyse linear systems with time-varying parameters; facilitate problem programming; analyse high-order linear systems; handle multiple input–multiple output systems; and serve as a basis for studying nonlinear systems, stochastic systems, and optimal control. The benefits stem from the use of state space, which enhances the generality and rigour of classical transform theory. In summary, state space methods are a crucial tool in modern control theory, providing a more general and rigorous approach to analysing the time behaviour of physical processes [34]. Recently, many authors have used the state-space approach in many applications [35,36,37,38,39,40,41].
The novelty of this study is that it addresses an application of electromagnetic generalized thermoelastic half-space in the context of the Green–Naghdi type I and type III, Lord–Shulman, and Moore–Gibson–Thompson theorems when the half-space surface is exposed to the time-fractional Maxwell electromagnetic effect. The time-fractional Maxwell’s equations are considered using the Caputo fractional derivative formulation. The main goal of this work is to investigate the effect of the fractional-order parameter of the time-fractional Maxwell’s equations, which has not been studied before, and to introduce new results where the time-fractional Maxwell’s equation has not been applied before. Moreover, this work introduces unified governing equations of four different theorems of thermoelasticity which are somehow close to each other, so one of the benefits of this work is that it compares these four theorems.

2. The Basic Steps of the Work

Now, we will describe the steps of the work to shaw the formulation of the problem and obtain the solutions in the Laplace transform domain, then figure out the results and discussions as in Scheme 1.

3. Formulation of the Problem

Assuming that a generalized thermoelastic and electromagnetic half-space occupies the space (see Figure 1):
Ω = x , 0 , 0 : 0 x < .
H is a magnetic vector field with constant intensity which acts tangent to the bounding plane of the half-space x = 0 . All the studied functions of the material are functions only of the distance x and the time variable t .
An induced magnetic vector field h = 0 , h x , t , 0 is the result of the effect of the primary magnetic vector field H 0 = 0 , H 0 , 0 . As a sequence, an induced electric vector field E = E x , E y , E z is generated as shown in Figure 1 [17,42,43,44,45]:
We consider h and E to have small magnitudes based on the linear theory of Green–Naghdi thermoelasticity. Therefore, the vector of displacement will possess the following components:
U = U x , U y , U z = u x , t , 0 , 0 .
The magnetic intensity vector will have the following components [17,42,43,44,45]:
H 0 = H x , H y , H z = 0 , H 0 + h x , t , 0 .
According to the left-hand rule, the electric intensity vector field must be perpendicular to both the magnetic intensity vector and the displacement vectors. Then, E has the following components:
E = 0 , 0 , E x , t .
The current density vector J must be parallel to the electric intensity vector E ; hence, we have:
J = 0 , 0 , J .
The time-fractional Maxwell’s equations in general vector forms are as follows [28]:
× h = J + ε 0 α D t α E ,
× E = μ 0 α D t α h ,
h = 0 ,
E = 0 ,
and
B = μ 0 H 0 + h ,
where any symbol which is written in bold is a vector. μ 0 and ε 0 are the magnetic and electric permeabilities, respectively [5,17,28,42,45,46,47,48].
The Equations (6)–(10) are supplemented by Ohm’s law, namely [5,17,28,42,45,46,47,48]:
J = σ 0 μ 0 U t × H 0 + E .
Ohm’s law gives [5,17,28,42,45,46,47,48]:
J = σ 0 0 ,   0 ,   H 0 μ 0 u t + E ,
where σ 0 is the electric conductivity.
The well-known Lorentz force F is given by the following law [5,17,28,42,45,46,47,48]:
F = J × B = σ 0 μ 0 H 0 + h H 0 μ 0 u t + E ,   0 ,   0 .
After linearization, we obtain:
F = σ 0 μ 0 H 0 μ 0 H 0 u t + E ,   0 ,   0 .
The strain components have the following forms [17]:
e x x = u x , e y y = e z z = e x z = e y z = e x y = 0 ,
and
e = e x x + e y y + e z z = u x
The stress components are given by the following constitutive relation [17]:
σ i j = 2 μ e i j + λ e δ i j γ T T 0 δ i j ,
which has the components [17]:
σ x x = σ = λ + 2 μ u x γ T T 0 ,
σ y y = σ z z = λ u x γ T T 0 ,
and
σ x y = σ y z = σ z x = 0 ,
where T is the absolute temperature, T 0 is a reference temperature such that T T 0 T 0 1 , λ and μ are Lamé’s moduli, γ is given by γ = 3 λ + 2 μ α T , and α T is the coefficient of the thermal linear expansion.
Equations of motion have the following form [1,2,3]:
σ i j , j + F i = ρ u ¨ i ,    i = 1 , 2 , 3 ,
where ρ is the density of the material.
Substituting from Equation (12) into Equation (9), we obtain the partial differential equation of the motion in the following form [5,17,42,43,45,46,47,48,49]:
λ + 2 μ 2 u x 2 γ T x σ 0 μ 0 H 0 E σ 0 H 0 2 μ 0 2 u t = ρ 2 u t 2 ,
By exciting the partial derivative with respect to the variable x and using Equation (16), Equation (22) takes the following new form:
λ + 2 μ 2 e x 2 γ 2 T x 2 σ 0 μ 0 H 0 E x σ 0 H 0 2 μ 0 2 e t = ρ 2 e t 2 .
The one-dimensional heat conduction equation of Green–Naghdi type I and type III, Lord–Shulman, and Moore–Gibson–Thompson takes the following unified form [50,51,52]:
3 T x 2 t + K * K 2 T x 2 = 1 + τ q t ρ C E K 2 T t 2 + T 0 γ K 2 e t 2 .
The unified Equation (24) can be applied to the MGT theory, the Lord and Shulman theory (LS), and the type I (GN-I), and type III (GN-III) Green–Naghdi equations [53]. Hence, we have:
(i)
K * = 0   and   τ q = 0 provides the Green–Naghdi (GN-I) heat conduction equation.
(ii)
K * = 0   and   τ q 0 provides the Lord–Shulman (LS) model heat conduction equation [54].
(iii)
K * 0   and   τ q = 0 provides the Green–Naghdi (GN-III) model heat conduction equation [53,54].
(iv)
K * 0   and   τ q 0 provides the Moore–Gibson–Thompson (MGT) model heat conduction equation [15].
We may take K * = C E λ + 2 μ 4 as the main character of the Green–Naghdi theory, where K gives the thermal conductivity and C E denotes the specific heat with constant deformation.
The above four theorems are used because they are well known in applications and they are close to each other in terms of explaining the propagation behaviour of thermal, mechanical, and magnetic waves and can be unified into an easily assembled thermal conductivity equation.
Equation (6) for the current model will take the following form:
h x = J + ε 0 α D t α E .
Substituting from Equation (12), we obtain [28]:
h x = σ 0 + ε 0 α D t α E + σ 0 H 0 μ 0 u t .
Moreover, Equation (7) will take the following form [28]:
E x = μ 0 α D t α h .
By executing the partial derivative with respect to the variable x for Equation (26) and using Equation (27), we obtain the following equation [28]:
2 h x 2 = σ 0 + ε 0 α D t α μ 0 α D t α h + σ 0 H 0 μ 0 e t .
Moreover, Equation (23) takes the following form [28]:
λ + 2 μ 2 e x 2 γ 2 T x 2 σ 0 H 0 μ 0 α + 1 D t α h σ 0 H 0 2 μ 0 2 e t = ρ 2 e t 2 .
The following dimensionless variables are used for simplifications, [17,42,43,45,46]:
x , u = x , u c 0 η ,    t = t c 0 2 η ,    σ i j = λ + 2 μ σ i j ,   T T 0 = λ + 2 μ γ θ ,   E = σ 0 H 0 c 0 μ 0 2 η E ,   h = σ 0 H 0 μ 0 η h .
Equations (24) and (28)–(29) are reduced to the following system of differential equations (dropping the primes for convenience):
2 e x 2 2 θ x 2 ν ε ε 1 D t α h ε 1 e t = 2 e t 2 ,
3 T x 2 t + ε 2 2 T x 2 = 1 + τ q t 2 T t 2 + ε 3 2 e t 2 ,
2 h x 2 = ν ε D t α h + V 2 ε ε 4 D t α D t α h + e t ,
and
E x = ε D t α h ,
where
η = ρ C E K ,   c 2 = 1 ε 0 μ 0 ,   ν = σ 0 μ 0 η , V = c 0 c ,    c 0 2 = λ + 2 μ ρ ,    ε = μ 0 c 0 2 η α 1 ,    ε 1 = H 0 2 μ 0 ν λ + 2 μ , ε 2 = K * K c 0 2 η ,    ε 3 = γ 2 T 0 η K λ + 2 μ ,   and   ε 4 = ε 0 c 0 2 η α 1 .
The constitutive equations are also reduced to the following forms:
σ x x = σ = e θ ,
and
σ y y = σ z z = β e θ
where β = λ λ + 2 μ .
Now, we define the operator D t α = d α d t α , which is a fractional-order derivative and is given by the normal derivative when α = 1 , and a Caputo fractional derivative when 1 < α < 2 , respectively, as in the following unified form [19,20,21,22,23,24,25,26,27,28,29,30,31]:
D t α f t = 1 Γ 2 α 0 t t ξ 1 α d 2 f ξ d t 2 d ξ ,     1 < α < 2 ,    t > 0 ,
and the Laplace transform of the above formula and its derivatives is given by [19,20,21,22,23,24,25,26,27,28,29,30,31]:
L D t α f t = s α f ¯ s s α 1 f 0 + s α f ( 1 ) 0 + ,     1 α < 2 .
The initial conditions of the current model are given as follows:
θ x , 0 + = e x , 0 + = h x , 0 + = θ x , 0 + t = e x , 0 + t = h x , 0 + t = 0 .
After applying the initial conditions, (37) is in the following simple form [19,20,21,22,23,24,25,26,27,28,29,30,31]:
L D t α f t = s α f ¯ s ,     1 α < 2 .
Then, we have:
d 2 d x 2 β 2 e ¯ d 2 θ ¯ d x 2 β 1 h ¯ = 0 ,
d 2 d x 2 β 3 θ ¯ β 4 e ¯ = 0 ,
d 2 d x 2 β 5 h ¯ β 6 e ¯ = 0 ,
and
d E ¯ d x = ε s α h ¯
where β 1 = ε ε 1 ν s α ,    β 2 = ε 1 s + s 2 ,   β 3 = s 2 1 + τ q s s + ε 2 ,   β 4 = s 2 1 + τ q s ε 3 s + ε 2 ,   β 5 = ε s α ν   + V 2 ε 4 s α , and β 6 = s .
The constitutive equations of stress–strain are given by:
σ ¯ x x = σ ¯ = e ¯ θ ¯ ,
and
σ ¯ y y = σ ¯ z z = β e ¯ θ ¯

4. State-Space Formulation

We re-write the Equations (40)–(42) in the following forms:
d 2 e ¯ d x 2 = β 2 + β 4 e ¯ + β 3 θ ¯ + β 1 h ¯ ,
d 2 θ ¯ d x 2 = β 4 e ¯ + β 3 θ ¯ ,
d 2 h ¯ d x 2 = β 6 e ¯ + β 5 h ¯ .
Taking as the state variable the functions e ¯ , θ ¯ , h ¯ , d e ¯ d x , d θ ¯ d x , d h ¯ d x in the x-direction, Equations (46)–(48) can be written in matrix form using the Bahar–Hetnarski method [34,55,56]:
d ψ ¯ ( x , s ) d x = A ( s ) ψ ¯ ( x , s ) ,
where
ψ ¯ x , s = e ¯ x , s θ ¯ x , s h ¯ x , s e ¯ x , s θ ¯ x , s h ¯ x , s
and
A s = 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 β 2 + β 4 β 3 β 1 0 0 0 β 4 β 3 0 0 0 0 β 6 0 β 5 0 0 0
The formal and bounded solution in the domain (1) is given by:
ψ ¯ x , s = exp A s x ψ ¯ 0 , s ,
where
ψ ¯ 0 , s = e ¯ 0 , s θ ¯ 0 , s h ¯ 0 , s e ¯ 0 , s θ ¯ 0 , s h ¯ 0 , s
The characteristic equation of the matrix A s takes the following form:
k 6 L k 4 + M k 2 N = 0 ,
where L = β 2 + β 3 + β 4 + β 5 ,   M = β 3 β 5 + β 2 β 3 + β 2 β 5 + β 4 β 5 β 1 β 6 ,   N = β 3 β 2 β 5 β 1 β 6 .
The expansion of the Taylor series for the matrix exponential exp A s x is given by:
exp A s x = i = 0 A s x i i ! .
Applying the Cayley–Hamilton theorem [55,56,57], this infinite series can be truncated to the following form:
exp A s x = a 0 I 6 + a 1 A + a 2 A 2 + a 3 A 3 + a 4 A 4 + a 5 A 5 ,
where I 6 is the unit matrix of order 6 and a 0 a 5 are some parameters which depend on s and x to be determined.
Applying the Cayley–Hamilton theorem again, we obtain [55,56,57]:
exp k 1 x = a 0 I 6 + a 1 k 1 + a 2 k 1 2 + a 3 k 1 3 + a 4 k 1 4 + a 5 k 1 5 ,
exp k 1 x = a 0 I 6 a 1 k 1 + a 2 k 1 2 a 3 k 1 3 + a 4 k 1 4 a 5 k 1 5 ,
exp k 2 x = a 0 I 6 + a 1 k 2 + a 2 k 2 2 + a 3 k 2 3 + a 4 k 2 4 + a 5 k 2 5 ,
exp k 2 x = a 0 I 6 a 1 k 2 + a 2 k 2 2 a 3 k 2 3 + a 4 k 2 4 a 5 k 2 5 ,
exp k 3 x = a 0 I 6 + a 1 k 3 + a 2 k 3 2 + a 3 k 3 3 + a 4 k 3 4 + a 5 k 3 5 ,
exp k 3 x = a 0 I 6 a 1 k 3 + a 2 k 3 2 a 3 k 3 3 + a 4 k 3 4 a 5 k 3 5 ,
where ± k 1 , ± k 2 , ± k 3 are the roots of the characteristic Equation (54).
By solving the system of Equations (57)–(62), we obtain the parameters a 0 , a 1 , a 2 , a 3 , a 4 ,   a 5 as follows:
a 0 = ξ c 1 k 2 2 k 3 2 + c 2 k 1 2 k 3 2 + c 3 k 2 2 k 1 2 ,
a 1 = ξ s 1 k 2 2 k 3 2 + s 2 k 1 2 k 3 2 + s 3 k 2 2 k 1 2 ,
a 2 = ξ c 1 k 2 2 + k 3 2 + c 2 k 1 2 + k 3 2 + c 3 k 2 2 + k 1 2 ,
a 3 = ξ s 1 k 2 2 + k 3 2 + s 2 k 1 2 + k 3 2 + s 3 k 2 2 + k 1 2 ,
a 4 = ξ c 1 + c 2 + c 3 ,
a 5 = ξ s 1 + s 2 + s 3 ,
where c 1 = k 2 2 k 3 2 cosh k 1 x , c 2 = k 3 2 k 1 2 cosh k 2 x , c 3 = k 1 2 k 2 2 cosh k 3 x , s 1 = k 2 2 k 3 2 k 1 sinh k 1 x , s 2 = k 3 2 k 1 2 k 2 sinh k 2 x , s 3 = k 1 2 k 2 2 k 3 sinh k 3 x , ξ = k 1 2 k 2 2 1 k 2 2 k 3 2 1 k 3 2 k 1 2 1 .
Substituting Equations (63)–(68) into Equation (56), we obtain the matrix exponential in the form:
exp A s x = l x , s = l i j x , s ,    i , j = 1 , 2 , 3 , 4 , 5 , 6 ,
where the elements of the matrix l i j x , s ,    i , j = 1 , 2 , 3 , 4 , 5 , 6 are defined in the Appendix A.
Thus, the solution of the Equation (52) is given by:
ψ ¯ x , s = l i j x , s ψ ¯ 0 , s .
Now, to obtain the column matrix ψ ¯ 0 , s , we must apply the following boundary conditions:
(i)
The thermal boundary condition: the bounding surface of the half-space x = 0 is thermally loaded by a ramp-type heat, i.e.,
θ 0 , x = θ 0 t / t 0 0 t < t 0 1 t t 0 = g t ,
where t 0 is the parameter of the ramp-time heat. θ 0 is constant and gives the thermal loading intensity.
Using the Laplace transform, the above condition takes the form:
θ ¯ 0 , s = 1 e t 0 s t 0 s 2 = g ¯ s .
(ii)
The mechanical boundary condition: the bounding plane of the half-space x = 0 traction free, i.e.,
σ 0 , t = 0   which   gives   σ ¯ 0 , s = 0 .
Hence, from Equations (44) and (72), we obtain:
e ¯ 0 , s = g ¯ s .
(iii)
The electromagnetic boundary condition: the magnetic and electric functions h and E must satisfy the continuity conditions as follows [17,45]:
h x , t x = 0 = h 0 x , t x 0   and    E x , t x = 0 = E 0 x , t x 0 ,
where E 0 x , t   and   h 0 x , t are the electric and magnetic intensities in the free space, respectively.
In the free space, using the non-dimensional Maxwell’s equations in the Laplace transform domain, we can establish ν = 0   and   α = 1 , which give the following equations:
h ¯ 0 x = s V 2 E ¯ 0 ,
and
E ¯ 0 x = s h ¯ 0
Eliminating E ¯ 0 between the above two equations, we obtain:
d 2 h ¯ 0 d x 2 = δ 2 h ¯ 0    for   x 0 ,
where δ = V s .
The general and bounded solution of the Equation (78) is given by:
h ¯ 0 x , s = A e δ x   ,    x 0 .
From Equations (79) and (77), we obtain:
E ¯ 0 = s δ A e δ x   ,    x 0 .
From Equations (80) and (76), we obtain the condition when x = 0 as [17,45]:
d h ¯ x , s d x δ h ¯ x , s x = 0 = 0 .
Using the boundary conditions (73), (72), and (81), we obtain the following:
e ¯ x , s θ ¯ x , s h ¯ x , s e ¯ x , s θ ¯ x , s h ¯ x , s = l i j x , s 6 × 6 g ¯ s g ¯ s h ¯ 0 , s e ¯ 0 , s θ ¯ 0 , s h ¯ 0 , s .
To obtain e ¯ 0 , s , θ ¯ 0 , s , h ¯ 0 , s , we can use the Equations (81) and (82) when x = 0 in the following form:
e ¯ 0 , s θ ¯ 0 , s h ¯ 0 , s e ¯ 0 , s θ ¯ 0 , s δ h ¯ 0 , s = l i j 0 , s 6 × 6 g ¯ s g ¯ s h ¯ 0 , s e ¯ 0 , s θ ¯ 0 , s δ h ¯ 0 , s ,
which gives the following equations:
1 l 44 0 e ¯ 0 , s l 45 0 θ ¯ 0 , s l 43 0 + l 46 0 δ h ¯ 0 , s = l 41 0 + l 42 0 g ¯ s ,
l 54 0 e ¯ 0 , s 1 l 55 0 θ ¯ 0 , s + l 53 0 + l 56 0 δ h ¯ 0 , s = l 51 0 + l 52 0 g ¯ s ,
l 64 0 e ¯ 0 , s + l 65 0 θ ¯ 0 , s 1 l 63 0 δ l 66 0 h ¯ 0 , s = l 61 0 + l 62 0 g ¯ s ,
where l i j 0 = l i j 0 , s ,    i , j = 1 , 2 , 3 , 4 , 5 , 6 .
After obtaining the solution of the above system of equations and using Equation (81), we can obtain the known e ¯ 0 , s , θ ¯ 0 , s , h ¯ 0 , s , and h ¯ 0 , s .
Now, the complete solutions are obtained in the domain of Laplace transform.
The inversions of the Laplace transforms can be obtained by applying the following iteration form of Tzou [58]:
f ( t ) = L 1 f ¯ s e κ t t 1 2 f ¯ κ + Re r = 1 R 1 r f ¯ κ + i r π t ,
where i = 1 is the imaginary number unit, “ Re ” denotes the real part of a complex function, and R is an integer parameter that can be chosen such that:
g t R + 1 g t R < 10 6 .
To obtain a faster convergence for the above iteration, some experiments verify that the parameter “ κ ” may satisfy the following relation κ t 4.7 [57,58,59].

5. The Numerical Results

The copper material was used to obtain the numerical calculations. The parameters and the material properties’ constants were taken as follows [4,5,17,42,43,44,45,46,47,49]:
λ = 7.76 × 10 10 ,   μ = 3.86 × 10 10 ,   T 0 = 293 ,   ρ = 9854 , α T = 1.78 × 10 5 ,   K = 386 ,    C E = 381.0 ,   σ 0 = 5.7 × 10 7 ,   ε 0 = 10 9 36 π ,    μ 0 = 4 π × 10 7 .
The dimensionless variables were taken as follows: t = t 0 = 1.0 and 0.0 x 2.0 .
The calculations by using MAPLE software will be go on as in the Scheme 2 as follows:
Figure 2 contains four figures showing the temperature increment distributions under the theorems GN-I in Figure 2a, LS in Figure 2b, GN-III in Figure 2c, and MGT in Figure 2d, with variance values of the time-fractional parameter of Maxwell’s equations to establish this parameter of the temperature increment. We can see that the temperature increment distributions almost are the same under the four studied models. Moreover, the effect of the time-fractional parameter of Maxwell’s equations is limited.
Figure 3 contains four figures showing the volumetric dilatation distributions under the theorems GN-I in Figure 3a, LS in Figure 3b, GN-III in Figure 3c, and MGT in Figure 3d, along with variance values of the time-fractional parameter of Maxwell’s equations to establish this parameter of the volumetric dilatation. We can see that the volumetric dilatation distributions have the same behaviour with different values under the four studied models. Moreover, the effect of the time-fractional parameter of Maxwell’s equations is significant; increasing the values of the time-fractional parameter of Maxwell’s equations leads to a decrease in the value of the volumetric dilatation. In other words, the time-fractional parameter of Maxwell’s equations works as a resistor to the volumetric dilatation. In addition, the maximum values of the volumetric dilatation distributions take the following order:
e Max LS > e Max GN - I > e Max GN - III = e Max MGT
Also, we can say that the behaviour of the volumetric dilatation distributions under the two models of Green–Naghdi type III and Moore–Gibson–Thompson are too close to each other.
Figure 4 presents four figures showing the displacement distributions under the theorems GN-I in Figure 4a, LS in Figure 4b, GN-III in Figure 4c, and MGT in Figure 4d, along with variance values of the time-fractional parameter of Maxwell’s equations to establish this parameter of the displacement. We can see that the displacement distributions have the same behaviour with different values under the four studied models. Moreover, the effects of the time-fractional parameter of Maxwell’s equations are significant; increasing the values of the time-fractional parameter of Maxwell’s equations leads to a decrease in the absolute value of the displacement. This means that the time-fractional parameter of Maxwell’s equations works as a resistor to the displacement. In addition, the maximum values of the absolute value of the displacement distributions take the following order:
u Max LS > u Max GN - I > u Max GN - III = u Max MGT
Thus, the behaviour of the displacement distributions under the two models of Green–Naghdi type III and Moore–Gibson–Thompson are too close to each other.
Figure 5 presents four figures showing the stress distributions under the theorems GN-I in Figure 5a, LS in Figure 5b, GN-III in Figure 5c, and MGT in Figure 5d, along with variance values of the time-fractional parameter of Maxwell’s equations to establish this parameter of the stress. We can see that the stress distributions have the same behaviour with different values under the four studied models. Also, the effects of the time-fractional parameter of Maxwell’s equations are significant; increasing the values of the time-fractional parameter of Maxwell’s equations leads to a decrease in the value of the stress. This means that the time-fractional parameter of Maxwell’s equations works as a resistor to the stress. In addition, the maximum values of the stress distributions take the following order:
σ Max LS > σ Max GN - I > σ Max GN - III = σ Max MGT
Thus, the behaviour of the stress distributions under the two models of Green–Naghdi type III and Moore–Gibson–Thompson are too close to each other.
Figure 6 presents four figures showing the induced magnetic field distributions under the theorems GN-I in Figure 6a, LS in Figure 6b, GN-III in Figure 6c, and MGT in Figure 6d, along with variance values of the time-fractional parameter of Maxwell’s equations to establish this parameter on the induced magnetic field. We can see that the induced magnetic field distributions have the same behaviour with different values under the four studied models. Also, the effects of the time-fractional parameter of Maxwell’s equations are significant; increasing the values of the time-fractional parameter of Maxwell’s equations leads to an increase in the absolute value of the induced magnetic field. In other words, the time-fractional parameter of Maxwell’s equations works as a resistor to the induced magnetic field. In addition, the maximum values of the absolute value of the induced magnetic field distributions in the boundary x = 0 take the following order:
h Max LS = h Max GN - I > h Max GN - III > h Max MGT
Thus, the behaviour of the displacement distributions under the two models of Lord–Shulman and Green–Naghdi type I are too close to each other.
Figure 7 presents four figures showing the induced electric field distributions under the theorems GN-I in Figure 7a, LS in Figure 7b, GN-III in Figure 7c, and MGT in Figure 7d, along with variance values of the time-fractional parameter of Maxwell’s equations to establish this parameter on the induced electric field. We can see that the induced electric field distributions show the same behaviour with different values under the four studied models. Also, the effects of the time-fractional parameter of Maxwell’s equations are significant; increasing the values of the time-fractional parameter of Maxwell’s equations leads to an increase in the value of the induced electric field. In other words, the time-fractional parameter of Maxwell’s equations works as a catalyst to the induced electric field. In addition, the maximum values of the induced electric field distributions in the boundary x = 0 take the following order:
E Max LS = E Max GN - I > E Max GN - III = E Max MGT
Thus, the behaviour of the stress distributions under the two models of Lord–Shulman and Green–Naghdi type I are too close to each other.
Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 show the temperature increment, volumetric dilatation, displacement, stress, induced magnetic field, and induced electric field distributions, respectively, when α = 1.2 under the theorems of GN-I, LS, GN-III, and MGT.
Figure 8 shows that the temperature increment distributions under the four studied theorems have the same behaviour, but with different values. The temperature increment values have the following order:
θ GN - III > θ MGT > θ GN - I > θ LS
This means that the speed of propagation of the thermal wave under the Lord–Shulman theorem has a limited value, more than the other theorems, and goes to zero faster.
Figure 9 shows that the volumetric dilatation distributions under the four studied theorems have the same behaviour, but with different values. The volumetric dilatation values have the following order:
e LS > e GN - I > e MGT > e GN - III
Figure 10 shows that the displacement distributions under the four studied theorems have the same behaviour, but with different values. The absolute values of displacement have the following order:
u LS > u GN - I > u MGT > u GN - III
Figure 11 shows that the stress distributions under the four studied theorems have the same behaviour, but with different values. The values of stress have the following order:
ρ LS > ρ GN - I > σ MGT > σ GN - III
Figure 9, Figure 10 and Figure 11 show that the speed of the mechanical wave under the theorem of GN-III has a limited value, which is higher than its values under the other theorems.
Figure 12 shows that the induced magnetic field distributions under the four studied theorems have the same behaviour, but with different values. The values of the induced magnetic field have the following order on the bounding plane of the half-space x = 0 :
h LS > h GN - I > h MGT > h GN - III
Figure 13 shows that the induced electric field distributions under the four studied theorems have the same behaviour, but with different values. The values of the induced electric field have the following order on the bounding plane of the half-space x = 0 :
E LS > E GN - I > E MGT > E GN - III
Figure 12 and Figure 13 show that the speed of the electromagnetic waves under the theorem of GN-III has a limited value, which is higher than its values under the other theorems.
For the validation of the current results, the current results of the case α = 1.0 show the same behaviour regarding the temperature increment, volumetric dilatation, displacement, stress, induced magnetic field, and induced electric field distributions seen in previously published papers [5,10,17,24,43,44,45,46,47,48].

6. Conclusions

This work introduced a new mathematical model of a thermoelastic and electromagnetic half-space based on Green–Naghdi theory type I (GN-I) and type III (GN-III), the Lord–Shulman (LS) theorem, and the Moore–Gibson–Thompson (MGT) theorem of thermoelasticity based on time-fractional Maxwell’s equations. The Caputo fractional derivative is the fractional derivative which was applied in this model.
The solutions were obtained directly using the state-space approach, with the Laplace transform and the general solutions for any set of boundary conditions obtained in the Laplace transform domain. Tzou’s iteration method was used to compute the inverse Laplace transforms. The distributions of temperature increment, volumetric dilatation, displacement, stress, induced magnetic field, and induced electric field were discovered and discussed.
The time-fractional Maxwell’s equation parameter based on the Caputo fractional derivative has significant influence on all the functions under investigation.
The time-fractional parameter of Maxwell’s equations functions as a hindrance to the particle’s displacement, deformation, and the induced magnetic field, but it catalyses the generated electric field inside the material.
The speed of propagation of the thermal wave under the Lord–Shulman theorem has a limited value, which is higher than the other theorems, and goes to zero faster, which makes the Lord–Shulman model the most successful model for the thermal transfer.
The speed of the mechanical and electromagnetic waves under the theorem of GN-III has a limited value, which is higher than its values under the other theorems. This means that the theorem of Green–Naghdi type III is the most successful model in explaining the mechanical and electromagnetic behaviour.
The time-fractional parameter of Maxwell’s equations works as a resistor to the volumetric dilatation, displacement, stress, and induced magnetic field, while it works as catalyst to the induced electric field.
Therefore, the time-fractional parameter of Maxwell’s equations could be used to control the thermomechanical and electromagnetic waves through thermoelastic materials.

Author Contributions

Conceptualization, H.M.Y.; Validation, H.M.Y. and E.A.N.A.-L.; Formal analysis, H.M.Y. and E.A.N.A.-L.; Investigation, H.M.Y.; Resources, E.A.N.A.-L.; Writing—original draft, H.M.Y.; Writing—review & editing, H.M.Y.; Visualization, E.A.N.A.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

l 11 x , s = β 2 2 + 2 β 2 + β 3 β 4 + β 1 β 6 + β 2 2 a 4 + a 2 β 2 + β 4 + a 0 ,
l 12 x , s = β 2 + β 3 + β 4 a 4 + a 2 β 3 ,
l 13 x , s = β 2 + β 4 + β 5 a 4 + a 2 β 1 ,
l 14 x , s = β 2 2 + 2 β 2 + β 3 β 4 + β 1 β 6 + β 2 2 a 5 + a 3 β 2 + β 4 + a 1 ,
l 15 x , s = β 2 + β 4 + β 5 a 5 + a 3 β 3 ,
l 16 x , s = β 2 + β 4 + β 5 a 5 + a 3 β 1 ,
l 21 x , s = β 2 + β 3 + β 4 a 4 + a 2 β 4 ,
l 22 x , s = a 4 β 3 2 + β 3 a 4 β 4 + a 2 β 3 + a 0 ,
l 23 x , s = a 4 β 1 β 4 ,
l 24 x , s = β 2 + β 3 + β 4 a 5 + a 3 β 4 ,
l 25 x , s = a 5 β 3 2 + β 3 a 5 β 4 + a 3 β 3 + a 1 ,
l 26 x , s = a 5 β 1 β 4 ,
l 31 x , s = β 2 + β 4 + β 5 a 4 + a 2 β 6 ,
l 32 x , s = a 4 β 3 β 6 ,
l 33 = a 4 β 1 β 6 + β 5 2 + a 2 β 5 + a 0 ,
l 34 x , s = β 2 + β 4 + β 5 a 5 + a 3 β 6 ,
l 35 x , s = a 5 β 3 β 6 ,
l 36 x , s = a 5 β 1 β 6 + β 5 2 + a 3 β 5 + a 1 ,
l 41 x , s = a 5 β 4 3 + 3 β 2 + 2 β 3 a 5 + a 3 β 4 2 + 2 β 1 β 6 + 3 β 2 2 + 2 β 3 β 2 + β 3 2 a 5 + 2 a 3 β 2 + a 3 β 3 + a 1 β 4        + 2 β 1 β 2 β 6 + β 1 β 5 β 6 + β 2 3 a 5 + a 3 β 1 β 6 + a 3 β 2 2 + a 1 β 2
l 42 x , s = β 3 β 3 2 + β 2 + 2 β 4 β 3 + β 1 β 6 + β 2 2 + 2 β 4 β 2 + β 4 2 a 5 + a 3 β 2 + β 3 + β 4 + a 1 ,
l 43 x , s = β 1 β 4 2 + 2 β 2 + β 3 + β 5 β 4 + β 1 β 6 + β 2 2 + β 5 β 2 + β 5 2 a 5 + a 3 β 2 + β 3 + β 4 + a 1 ,
l 44 x , s = β 4 2 + 2 β 2 + β 3 β 4 + β 1 β 6 + β 2 2 a 4 + a 2 β 2 + a 2 β 4 + a 0 ,
l 45 x , s = β 2 + β 3 + β 4 a 4 + a 2 β 3 ,
l 46 x , s = β 2 + β 4 + β 5 a 4 + a 2 β 1 ,
l 51 x , s = β 4 β 2 2 + 2 β 2 + 2 β 3 β 4 + β 1 β 6 + β 2 2 + β 3 β 2 + β 3 2 a 5 + a 3 β 2 + β 3 + β 4 + a 1 ,
l 52 x , s = β 3 a 5 β 3 2 + 2 a 5 β 4 + a 3 β 3 + a 5 β 4 2 + a 5 β 2 + a 3 β 4 + a 1 ,
l 53 x , s = β 2 + β 3 + β 4 + β 5 a 5 + a 3 β 1 β 4 ,
l 54 x , s = β 2 + β 3 + β 4 a 4 + a 2 β 4 ,
l 55 x , s = a 4 β 3 2 + a 4 β 4 + a 2 β 3 + a 0 ,
l 56 x , s = a 4 β 1 β 4 ,
l 61 x , s = β 6 β 4 2 + 2 β 2 + β 3 + β 5 β 4 + β 1 β 6 + β 2 2 + β 5 β 2 + β 5 2 a 5 + a 3 β 2 + β 4 + β 5 + a 1 ,
l 62 x , s = β 2 + β 3 + β 4 + β 5 a 5 + a 3 β 3 β 6 ,
l 63 x , s = β 1 β 6 2 a 5 β 5 + β 2 + β 4 a 5 + a 3 + β 5 a 5 β 5 2 + a 3 β 5 + a 1 ,
l 64 x , s = β 2 + β 4 + β 5 a 4 + a 2 β 6 ,
l 65 x , s = a 4 β 3 β 6 ,
l 66 x , s = a 4 β 1 β 6 + β 5 2 + a 2 β 5 + a 0

References

  1. Hetnarski, R.B.; Eslami, M.R.; Hetnarski, R.B.; Eslami, M.R. Basic laws of thermoelasticity. In Thermal Stresses—Advanced Theory and Applications; Springer: Cham, Switzerland, 2019; pp. 1–43. [Google Scholar]
  2. Biot, M.A. Thermoelasticity and irreversible thermodynamics. J. Appl. Phys. 1956, 27, 240–253. [Google Scholar] [CrossRef]
  3. Lord, H.W.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  4. Sherief, H.H.; Darwish, A.A. A short time solution for a problem in thermoelasticity of an infinite medium with a spherical cavity. J. Therm. Stress. 1998, 21, 811–828. [Google Scholar] [CrossRef]
  5. Sadeghi, M.; Kiani, Y. Generalized magneto-thermoelasticity of a layer based on the Lord–Shulman and Green–Lindsay theories. J. Therm. Stress. 2022, 45, 319–340. [Google Scholar] [CrossRef]
  6. Youssef, H.M.; El-Bary, A.A. The influence of the mechanical damage on a viscothermoelastic nanobeam due to ramp-type heating under Green-Naghdi theory type-II. J. Strain Anal. Eng. Des. 2024, 59, 44–55. [Google Scholar] [CrossRef]
  7. Hendy, M.M.; Ezzat, M.A. A modified Green-Naghdi fractional order model for analyzing thermoelectric MHD. Int. J. Numer. Methods Heat Fluid Flow 2024, 34, 2376–2398. [Google Scholar] [CrossRef]
  8. Sarkar, N.; Lahiri, A. A three-dimensional thermoelastic problem for a half-space without energy dissipation. Int. J. Eng. Sci. 2012, 51, 310–325. [Google Scholar] [CrossRef]
  9. Rezazadeh, M.; Tahani, M.; Hosseini, S.M. Thermoelastic damping in a nonlocal nano-beam resonator as NEMS based on the type III of Green–Naghdi theory (with energy dissipation). Int. J. Mech. Sci. 2015, 92, 304–311. [Google Scholar] [CrossRef]
  10. Atwa, S.Y. Generalized magneto-thermoelasticity with two temperature and initial stress under Green–Naghdi theory. Appl. Math. Model. 2014, 38, 5217–5230. [Google Scholar] [CrossRef]
  11. Othman, M.I.; Atwa, S.Y.; Farouk, R. The effect of diffusion on two-dimensional problem of generalized thermoelasticity with Green–Naghdi theory. Int. Commun. Heat Mass Transf. 2009, 36, 857–864. [Google Scholar] [CrossRef]
  12. Sharma, S.; Sharma, K.; Bhargava, R.R. Effect of viscosity on wave propagation in anisotropic thermoelastic with Green-Naghdi theory type-II and type-III. Mater. Phys. Mech. 2013, 16, 144–158. [Google Scholar]
  13. Lizama, C.; Warma, M.; Zamorano, S. Exterior controllability properties for a fractional Moore–Gibson–Thompson equation. Fract. Calc. Appl. Anal. 2022, 25, 887–923. [Google Scholar] [CrossRef]
  14. Liu, W.; Chen, Z. General decay rate for a Moore–Gibson–Thompson equation with infinite history. Z. Für Angew. Math. Und Phys. 2020, 71, 43. [Google Scholar] [CrossRef]
  15. Quintanilla, R. Moore–Gibson–Thompson thermoelasticity. Math. Mech. Solids 2019, 24, 4020–4031. [Google Scholar] [CrossRef]
  16. Iqbal, K.; Kulvinder, S. The two-temperature effect on a semiconducting thermoelastic solid cylinder based on the modified Moore–Gibson–Thompson heat transfer. St. Petersburg State Polytehnical Univ. J. Phys. Math. 2023, 16, 65–81. [Google Scholar]
  17. Sherief, H.H.; Yosef, H.M. Short time solution for a problem in magnetothermoelasticity with thermal relaxation. J. Therm. Stress. 2004, 27, 537–559. [Google Scholar] [CrossRef]
  18. Youssef, H.M.; Alghamdi, A.A. Influence of the fractional-order strain on an infinite material with a spherical cavity under Green-Naghdi hyperbolic two-temperature thermoelasticity theory. J. Eng. Therm. Sci. 2023, 3, 11–24. [Google Scholar] [CrossRef]
  19. Daftardar-Gejji, V. Fractional Calculus; Alpha Science International Limited: New Delhi, India, 2013. [Google Scholar]
  20. Li, C.; Cai, M. Theory and Numerical Approximations of Fractional Integrals and Derivatives; SIAM: Philadelphia, PA, USA, 2019. [Google Scholar]
  21. Almeida, R.; Tavares, D.; Torres, D.F. The Variable-Order Fractional Calculus of Variations; Springer: Cham, Switzerland, 2019. [Google Scholar]
  22. Petrás, I. Fractional Derivatives, Fractional Integrals, and Fractional Differential Equations in Matlab; IntechOpen: London, UK, 2011. [Google Scholar]
  23. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  24. Baleanu, D.; Golmankhaneh, A.K.; Golmankhaneh, A.K.; Baleanu, M.C. Fractional electromagnetic equations using fractional forms. Int. J. Theor. Phys. 2009, 48, 3114–3123. [Google Scholar] [CrossRef]
  25. Lazo, M.J. Gauge invariant fractional electromagnetic fields. Phys. Lett. A 2011, 375, 3541–3546. [Google Scholar] [CrossRef]
  26. Jaradat, E.; Hijjawi, R.; Khalifeh, J. Maxwell’s equations and electromagnetic Lagrangian density in fractional form. J. Math. Phys. 2012, 53, 033505. [Google Scholar] [CrossRef]
  27. Ortigueira, M.D.; Rivero, M.; Trujillo, J.J. From a generalised Helmholtz decomposition theorem to fractional Maxwell equations. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 1036–1049. [Google Scholar] [CrossRef]
  28. Stefański, T.P.; Gulgowski, J. Formulation of time-fractional electrodynamics based on Riemann-Silberstein vector. Entropy 2021, 23, 987. [Google Scholar] [CrossRef] [PubMed]
  29. Machado, J.T.; Jesus, I.S.; Galhano, A.; Cunha, J.B. Fractional order electromagnetics. Signal Process. 2006, 86, 2637–2644. [Google Scholar] [CrossRef]
  30. Youssef, H.M. Theory of fractional order generalized thermoelasticity. J. Heat Transf. 2010, 132, 061301. [Google Scholar] [CrossRef]
  31. Youssef, H.M. Theory of generalized thermoelasticity with fractional order strain. J. Vib. Control. 2016, 22, 3840–3857. [Google Scholar] [CrossRef]
  32. Üzar, N. Chaotic and regular behaviours of classical and fractional Gross–Pitaevskii equations including two-body, three-body and higher-order interactions. Pramana 2023, 97, 36. [Google Scholar] [CrossRef]
  33. Al-Raeei, M. Applying fractional quantum mechanics to systems with electrical screening effects. Chaos Solitons Fractals 2021, 150, 111209. [Google Scholar] [CrossRef]
  34. Ezzat, M.A. State space approach to solids and fluids. Can. J. Phys. 2008, 86, 1241–1250. [Google Scholar] [CrossRef]
  35. Mahato, C.S.; Biswas, S. State space approach to study thermal shock problem in nonlocal thermoelastic medium with double porosity. J. Therm. Stress. 2023, 46, 415–443. [Google Scholar] [CrossRef]
  36. Othman, M.I.; Abd-Elaziz, E.M.; Hilal, M.I. State-space approach to a 2-D generalized thermoelastic medium under the effect of inclined load and gravity using a dual-phase-lag model. Mech. Based Des. Struct. Mach. 2022, 50, 395–411. [Google Scholar] [CrossRef]
  37. Ezzat, M.A.; Ezzat, S.M.; Alkharraz, M.Y. State-space approach to nonlocal thermo-viscoelastic piezoelectric materials with fractional dual-phase lag heat transfer. Int. J. Numer. Methods Heat Fluid Flow 2022, 32, 3726–3750. [Google Scholar] [CrossRef]
  38. Mahato, C.S.; Biswas, S. State Space Approach to Characterize Rayleigh Waves in Nonlocal Thermoelastic Medium with Double Porosity under Three-Phase-Lag Model. Comput. Math. Math. Phys. 2024, 64, 555–584. [Google Scholar] [CrossRef]
  39. Youssef, H.M.; Al-Lehaibi, E.A. State-space approach to three-dimensional generalized thermoelasticity with fractional order strain. Mech. Adv. Mater. Struct. 2019, 26, 878–885. [Google Scholar] [CrossRef]
  40. Ismail, M.; Khamis, A.; El-Bary, A.; Youssef, H. Effect of the rotation of generalized thermoelastic layer subjected to harmonic heat: State-space approach. Microsyst. Technol. 2017, 23, 3381–3388. [Google Scholar] [CrossRef]
  41. Al-Lehaibi, E.A.; Youssef, H.M. State-space approach to nano-beam with variable material properties. Adv. Sci. Eng. Med. 2016, 8, 412–420. [Google Scholar] [CrossRef]
  42. Ezzat, M.A.; Youssef, H.M. Generalized magneto-thermoelasticity in a perfectly conducting medium. Int. J. Solids Struct. 2005, 42, 6319–6334. [Google Scholar] [CrossRef]
  43. Othman, M.I. Generalized Electromagneto-Thermoelastic Plane Waves by Thermal Shock Problem in a Finite Conductivity Half-Space with One Relaxation Time. Multidiscip. Model. Mater. Struct. 2005, 1, 231–250. [Google Scholar] [CrossRef]
  44. Abbas, I.A. Generalized magneto-thermoelasticity in a nonhomogeneous isotropic hollow cylinder using the finite element method. Arch. Appl. Mech. 2009, 79, 41–50. [Google Scholar] [CrossRef]
  45. Tiwari, R. Magneto-thermoelastic interactions in generalized thermoelastic half-space for varying thermal and electrical conductivity. Waves Random Complex Media 2024, 34, 1795–1811. [Google Scholar] [CrossRef]
  46. Othman, M.I. Relaxation effects on thermal shock problems in an elastic half-space of generalized magneto-thermoelastic waves. Mech. Mech. Eng. 2004, 7, 165–178. [Google Scholar]
  47. Ezzat, M.A.; Youssef, H.M. Generalized magneto-thermoelasticity for an infinite perfect conducting body with a cylindrical cavity. Mater. Phys. Mech. 2013, 18, 156–170. [Google Scholar]
  48. Sarkar, N. Generalized magneto-thermoelasticity with modified Ohm’s Law under three theories. Comput. Math. Model. 2014, 25, 544–564. [Google Scholar] [CrossRef]
  49. Ezzat, M.A.; El-Karamany, A.S.; El-Bary, A.A. Electro-magnetic waves in generalized thermo-viscoelasticity for different theories. Int. J. Appl. Electromagn. Mech. 2015, 47, 95–111. [Google Scholar] [CrossRef]
  50. Hobiny, A.; Abbas, I.A. A GL photo-thermal theory upon new hyperbolic two-temperatures in a semiconductor material. Waves Random Complex Media 2021, 34, 4799–4812. [Google Scholar] [CrossRef]
  51. Youssef, H.M.; El-Bary, A. Characterization of the photothermal interaction of a semiconducting solid sphere due to the mechanical damage, ramp-Type heating, and rotation under LS theory. Waves Random Complex Media 2024, 34, 547–572. [Google Scholar] [CrossRef]
  52. Hobiny, A.D.; Abbas, I.A. The influences of thermal relaxation time and varying thermal conductivity in thermoelastic media. Case Stud. Therm. Eng. 2024, 56, 104263. [Google Scholar] [CrossRef]
  53. Green, A.; Naghdi, P. Thermoelasticity without energy dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  54. Kumar, H.; Mukhopadhyay, S. Thermoelastic damping analysis in microbeam resonators based on Moore–Gibson–Thompson generalized thermoelasticity theory. Acta Mech. 2020, 231, 3003–3015. [Google Scholar] [CrossRef]
  55. Bahar, L.Y.; Hetnarski, R.B. State space afproach to thermoelasticity. J. Therm. Stress. 1978, 1, 135–145. [Google Scholar] [CrossRef]
  56. Modi, V. Proceedings of the Sixth Canadian Congress of Applied Mechanics (CANCAM 77), Vancouver, BC, Canada, 30 May–3 June 1977; Volume 2. Available online: https://ui.adsabs.harvard.edu/abs/1977STIA...7741951M/abstract (accessed on 23 August 2024).
  57. Özişik, M.; Tzou, D. On the wave theory in heat conduction. J. Heat Transf. 1994, 116, 526–535. [Google Scholar] [CrossRef]
  58. Tzou, D.; Guo, Z.-Y. Nonlocal behavior in thermal lagging. Int. J. Therm. Sci. 2010, 49, 1133–1137. [Google Scholar] [CrossRef]
  59. Tzou, D.Y. The generalized lagging response in small-scale and high-rate heating. Int. J. Heat Mass Transf. 1995, 38, 3231–3240. [Google Scholar] [CrossRef]
Scheme 1. The basic steps of the work.
Scheme 1. The basic steps of the work.
Fractalfract 08 00566 sch001
Figure 1. The electromagnetic half-space subjected to time-fractional Maxwell’s equations.
Figure 1. The electromagnetic half-space subjected to time-fractional Maxwell’s equations.
Fractalfract 08 00566 g001
Scheme 2. The calculation of numerical results.
Scheme 2. The calculation of numerical results.
Fractalfract 08 00566 sch002
Figure 2. The temperature increment distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The temperature increment distributions under the GN-I theorem. (b) The temperature increment distributions under the LS theorem. (c) The temperature increment distributions under the GN-III theorem. (d) The temperature increment distributions under the MGT theorem.
Figure 2. The temperature increment distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The temperature increment distributions under the GN-I theorem. (b) The temperature increment distributions under the LS theorem. (c) The temperature increment distributions under the GN-III theorem. (d) The temperature increment distributions under the MGT theorem.
Fractalfract 08 00566 g002
Figure 3. The volumetric dilatation distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The volumetric dilatation distributions under GN-I. (b) The volumetric dilatation distributions under LS. (c) The volumetric dilatation distributions under GN-III. (d) The volumetric dilatation distributions under MGT.
Figure 3. The volumetric dilatation distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The volumetric dilatation distributions under GN-I. (b) The volumetric dilatation distributions under LS. (c) The volumetric dilatation distributions under GN-III. (d) The volumetric dilatation distributions under MGT.
Fractalfract 08 00566 g003
Figure 4. The displacement distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The displacement distributions under GN-I. (b) The displacement distributions under LS. (c) The displacement distributions under GN-III. (d) The displacement distributions under MGT.
Figure 4. The displacement distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The displacement distributions under GN-I. (b) The displacement distributions under LS. (c) The displacement distributions under GN-III. (d) The displacement distributions under MGT.
Fractalfract 08 00566 g004
Figure 5. The stress distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The stress distributions under GN-I. (b) The stress distributions under LS. (c) The stress distributions under GN-III. (d) The stress distributions under MGT.
Figure 5. The stress distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The stress distributions under GN-I. (b) The stress distributions under LS. (c) The stress distributions under GN-III. (d) The stress distributions under MGT.
Fractalfract 08 00566 g005
Figure 6. The induced magnetic field distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The induced magnetic field distributions under GN-I. (b) The induced magnetic field distributions under LS. (c) The induced magnetic field distributions under GN-III. (d) The induced magnetic field distributions under MGT.
Figure 6. The induced magnetic field distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The induced magnetic field distributions under GN-I. (b) The induced magnetic field distributions under LS. (c) The induced magnetic field distributions under GN-III. (d) The induced magnetic field distributions under MGT.
Fractalfract 08 00566 g006
Figure 7. The induced electric field distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The induced electric field distributions under GN-I. (b) The induced electric field distributions under LS. (c) The induced electric field distributions under GN-III. (d) The induced electric field distributions under MGT.
Figure 7. The induced electric field distributions under four different thermoelastic theorems with variance values of the time-fractional parameter of Maxwell’s equations when t = t 0 = 1.0 . (a) The induced electric field distributions under GN-I. (b) The induced electric field distributions under LS. (c) The induced electric field distributions under GN-III. (d) The induced electric field distributions under MGT.
Fractalfract 08 00566 g007
Figure 8. The temperature increment distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Figure 8. The temperature increment distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Fractalfract 08 00566 g008
Figure 9. The volumetric dilatation distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Figure 9. The volumetric dilatation distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Fractalfract 08 00566 g009
Figure 10. The displacement distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Figure 10. The displacement distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Fractalfract 08 00566 g010
Figure 11. The stress distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and t = t 0 = 1.0 .
Figure 11. The stress distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and t = t 0 = 1.0 .
Fractalfract 08 00566 g011
Figure 12. The induced magnetic field distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Figure 12. The induced magnetic field distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Fractalfract 08 00566 g012
Figure 13. The induced electric field distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Figure 13. The induced electric field distribution under the four theorems of thermoelasticity when the value of the time-fractional parameter of Maxwell’s equations is α = 1.2 and when t = t 0 = 1.0 .
Fractalfract 08 00566 g013
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

Al-Lehaibi, E.A.N.; Youssef, H.M. State-Space Approach to the Time-Fractional Maxwell’s Equations under Caputo Fractional Derivative of an Electromagnetic Half-Space under Four Different Thermoelastic Theorems. Fractal Fract. 2024, 8, 566. https://doi.org/10.3390/fractalfract8100566

AMA Style

Al-Lehaibi EAN, Youssef HM. State-Space Approach to the Time-Fractional Maxwell’s Equations under Caputo Fractional Derivative of an Electromagnetic Half-Space under Four Different Thermoelastic Theorems. Fractal and Fractional. 2024; 8(10):566. https://doi.org/10.3390/fractalfract8100566

Chicago/Turabian Style

Al-Lehaibi, Eman A. N., and Hamdy M. Youssef. 2024. "State-Space Approach to the Time-Fractional Maxwell’s Equations under Caputo Fractional Derivative of an Electromagnetic Half-Space under Four Different Thermoelastic Theorems" Fractal and Fractional 8, no. 10: 566. https://doi.org/10.3390/fractalfract8100566

APA Style

Al-Lehaibi, E. A. N., & Youssef, H. M. (2024). State-Space Approach to the Time-Fractional Maxwell’s Equations under Caputo Fractional Derivative of an Electromagnetic Half-Space under Four Different Thermoelastic Theorems. Fractal and Fractional, 8(10), 566. https://doi.org/10.3390/fractalfract8100566

Article Metrics

Back to TopTop