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.
3. Formulation of the Problem
Assuming that a generalized thermoelastic and electromagnetic half-space occupies the space (see
Figure 1):
is a magnetic vector field with constant intensity which acts tangent to the bounding plane of the half-space . All the studied functions of the material are functions only of the distance and the time variable .
An induced magnetic vector field
is the result of the effect of the primary magnetic vector field
. As a sequence, an induced electric vector field
is generated as shown in
Figure 1 [
17,
42,
43,
44,
45]:
We consider
and
to have small magnitudes based on the linear theory of Green–Naghdi thermoelasticity. Therefore, the vector of displacement will possess the following components:
The magnetic intensity vector will have the following components [
17,
42,
43,
44,
45]:
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,
has the following components:
The current density vector
must be parallel to the electric intensity vector
; hence, we have:
The time-fractional Maxwell’s equations in general vector forms are as follows [
28]:
and
where any symbol which is written in bold is a vector.
and
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]:
Ohm’s law gives [
5,
17,
28,
42,
45,
46,
47,
48]:
where
is the electric conductivity.
The well-known Lorentz force
is given by the following law [
5,
17,
28,
42,
45,
46,
47,
48]:
After linearization, we obtain:
The strain components have the following forms [
17]:
and
The stress components are given by the following constitutive relation [
17]:
which has the components [
17]:
and
where
is the absolute temperature,
is a reference temperature such that
,
and
are Lamé’s moduli,
is given by
, and
is the coefficient of the thermal linear expansion.
Equations of motion have the following form [
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]:
By exciting the partial derivative with respect to the variable
and using Equation (16), Equation (22) takes the following new form:
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]:
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)
provides the Green–Naghdi (GN-I) heat conduction equation.
- (ii)
provides the Lord–Shulman (LS) model heat conduction equation [
54].
- (iii)
provides the Green–Naghdi (GN-III) model heat conduction equation [
53,
54].
- (iv)
provides the Moore–Gibson–Thompson (MGT) model heat conduction equation [
15].
We may take as the main character of the Green–Naghdi theory, where gives the thermal conductivity and 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:
Substituting from Equation (12), we obtain [
28]:
Moreover, Equation (7) will take the following form [
28]:
By executing the partial derivative with respect to the variable
for Equation (26) and using Equation (27), we obtain the following equation [
28]:
Moreover, Equation (23) takes the following form [
28]:
The following dimensionless variables are used for simplifications, [
17,
42,
43,
45,
46]:
Equations (24) and (28)–(29) are reduced to the following system of differential equations (dropping the primes for convenience):
and
where
The constitutive equations are also reduced to the following forms:
and
where
.
Now, we define the operator
, which is a fractional-order derivative and is given by the normal derivative when
, and a Caputo fractional derivative when
, respectively, as in the following unified form [
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31]:
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]:
The initial conditions of the current model are given as follows:
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]:
Then, we have:
and
where
and
The constitutive equations of stress–strain are given by:
and
4. State-Space Formulation
We re-write the Equations (40)–(42) in the following forms:
Taking as the state variable the functions
in the
x-direction, Equations (46)–(48) can be written in matrix form using the Bahar–Hetnarski method [
34,
55,
56]:
where
and
The formal and bounded solution in the domain (1) is given by:
where
The characteristic equation of the matrix
takes the following form:
where
.
The expansion of the Taylor series for the matrix exponential
is given by:
Applying the Cayley–Hamilton theorem [
55,
56,
57], this infinite series can be truncated to the following form:
where
is the unit matrix of order 6 and
are some parameters which depend on
s and
x to be determined.
Applying the Cayley–Hamilton theorem again, we obtain [
55,
56,
57]:
where
are the roots of the characteristic Equation (54).
By solving the system of Equations (57)–(62), we obtain the parameters
as follows:
where
,
,
,
,
,
,
.
Substituting Equations (63)–(68) into Equation (56), we obtain the matrix exponential in the form:
where the elements of the matrix
are defined in the
Appendix A.
Thus, the solution of the Equation (52) is given by:
Now, to obtain the column matrix , we must apply the following boundary conditions:
- (i)
The thermal boundary condition: the bounding surface of the half-space
is thermally loaded by a ramp-type heat, i.e.,
where
is the parameter of the ramp-time heat.
is constant and gives the thermal loading intensity.
Using the Laplace transform, the above condition takes the form:
- (ii)
The mechanical boundary condition: the bounding plane of the half-space
traction free, i.e.,
Hence, from Equations (44) and (72), we obtain:
- (iii)
The electromagnetic boundary condition: the magnetic and electric functions
and
must satisfy the continuity conditions as follows [
17,
45]:
where
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
, which give the following equations:
and
Eliminating
between the above two equations, we obtain:
where
.
The general and bounded solution of the Equation (78) is given by:
From Equations (79) and (77), we obtain:
From Equations (80) and (76), we obtain the condition when
as [
17,
45]:
Using the boundary conditions (73), (72), and (81), we obtain the following:
To obtain
, we can use the Equations (81) and (82) when
in the following form:
which gives the following equations:
where
.
After obtaining the solution of the above system of equations and using Equation (81), we can obtain the known , , , and .
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]:
where
is the imaginary number unit, “
” denotes the real part of a complex function, and
is an integer parameter that can be chosen such that:
To obtain a faster convergence for the above iteration, some experiments verify that the parameter “
” may satisfy the following relation
[
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]:
The dimensionless variables were taken as follows: and .
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:
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:
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:
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
take the following order:
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
take the following order:
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
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:
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:
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:
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:
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
:
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
:
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
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.