1. Introduction
In the generalized theories, the governing equations involve thermal relaxation times and they are of the hyperbolic type. Green and Naghdi [
1,
2,
3] considered a new extend theory by including the thermal displacement gradient between the constitutive variables. As we know, the classically coupled thermoelasticity includes the temperature gradient as one of the constitutive variables.
An important feature of this theory is that it does not accommodate the dissipation of thermal energy. In paper by Sharma and Chauhan [
4], we find an approach regarding the elastic interactions without considering the energy dissipation due to heat sources and body forces.
An important step in evolution of the classical theory of elasticity was made through the appearance of the theory of poroelasticity, which consider the volume of void, in an elastic body with pores, as a kinematics variable.
This gave the opportunity to investigate some concrete types of biological and geological solids and their useful applications. See, for instance, the applications in the fuel-cell industry [
5,
6,
7,
8,
9,
10].
We have to point out that the theory of linear elastic bodies with pores allows for the approach of such properties of biological and geological medium that could not be studied in the context of classical theory. It is very important to note that, when the volume of the pores tends to zero, we can see that the poroelastic theory reduces to the theory of classical elasticity.
In the paper, Nunziato and Cowin [
11] first established a theory of elastic bodies with pores in the non-linear theory case.
This theory of porous media has gained a great extension over the last period of time, and many authors consider different mathematical models for the mechanical behavior of solids with pores, by combining the poroelasticity theory with other different theories, in other words, combining different effects, [
12,
13,
14,
15,
16].
The consideration of the dynamic reaction of a thermoelastic body with additional parameters is very helpful in solving many concrete applications. For instance, the initial stresses are considered in a thermoelastic body with pores due to different reasons, such as the gravity variations, the difference of temperature, the process of quenching, etc.
Clearly, the earth is constantly under the influence of high initial stresses. As such, the researchers have allocated great importance to the study the effect of initial stresses regarding the thermal and mechanical state of a solid. For instance, Montanaro in [
17] investigated a thermoelastic isotropic body with hydrostatic initial stress.
Of course, the laser pulse has an effect on thermal loading in an elastic body with voids. Othman and Abd-Elaziz studied this effect in the paper [
18]. Marin investigated Cesaro means in the thermoelasticity of dipolar bodies [
19]. Marin and Oechsner studied the effect of a dipolar structure on the Holder stability in Green-Naghdi thermoelasticity [
20].
Other effects, such as the effect of the Earth’s electromagnetic field on seismic propagations, the designing of different elements of machine, emissions of electromagnetic radiations from nuclear devices, plasma physics, etc., can be found in [
21,
22,
23,
24,
25,
26,
27].
In our present study, we approach of the plane strain problem of a half-space body consisting of an electro-magneto-thermoelastic material that possesses voids and is subjected to some initial stress and to the Thomson effect. Our mathematical model is regarding the Green–Naghdi theory of type II and III of thermoelasticity. We assume that the Thomson effect is a constant coefficient and the density of charges that are induced by electric current is a function that depends on time variable.
In order to obtain the expressions for the considered parameters, it used the known normal mode technique. We also have obtained some graphic representations for the repartition of the considered variables.
2. Formulation of the Problem
An isotropic and homogeneous elastic body with pores (voids) is considered, with the temperature , in the reference state, and the half space . The motion referred to a rectangular Cartesian system of coordinates with origins in the surface . Additionally, the X-axis is pointing vertically into the body. In the of a two-dimensional problem, we suppose that the evolution of the body will be characterized by the displacement vector u, with components . The functions that are considered in this context are dependent on the time variable t and of the spatial variables x and y.
We consider a magnetic field with components , having a constant intensity, which acts parallel to the direction of the Z-axis.
It is known that a magnetic field of the form
produces an induced electric field of components
, and an induced magnetic field, as denoted by
h, and these satisfy the electromagnetism equations, in the linearized form. We will use the Maxwell’s equations [
24] in order to characterize the evolution of the electric field and for variation of the magnetic field, as follows:
The modified Ohm’s law for a medium with finite conductivity supplements the above system of coupled equations, namely
where
is the magnetic permeability,
is the magnetic displacement vector,
is the electric permeability,
is the current density vector,
is the charge density,
is the electric displacement vector, and
is the induced electric field vector.
For an isotropic and homogeneous thermoelastic body having pores, the constitutive equations receive the following form:
The strain-displacement relation is
The tensor of rotation has the components:
In Green-Naghdi (G-N) theories we take into account the Thomson effect, so that the Fourier’s law becomes
which gives
If we take into account Equations (1) and (3), then from Equation (13), we deduce
where
is the temperature above the reference temperature
is chosen so that
,
are the counterparts of Lame’ constants,
is the time,
are the components of the stress tensor,
is the equilibrated stress vector,
is the equilibrated inertia,
is the intrinsic equilibrated body force,
are constants of material that are due to the presence of the pores,
, such that
is the coefficient of thermal expansion,
is the Kronecker delta,
is the mass density,
is the specific heat at the constant strain,
is the thermal conductivity,
is entropy per unit mass,
is a constant, and
are the components of the first heat flux moment vector, we write the equation of continuity for the charges in the body in the form
where the velocity of the charges has the components
.
Let us now consider that the charge density is a function that does not depend on spatial variables, but only on time variable. Thus, Equation (15) will reduce to
We will assume that the charges have the speed of components
, which are proportional to the components of the velocity for particles
, so that we can write
which gives
where
is a positive constant (non-dimensional).
If we take into account Equation (18), from Equation (16), we are led to
which gives
Hence, we obtain
where
is the charge density when the strain vanishes.
While taking into account the Equation (22), from Equation (14) we deduce that the Fourier’s law, in its generalized form, receives the form:
In the case of null heat supply, the balance energy becomes
Taking into account Equations (9) and (23), from Equation (24), we deduce that the equation of heat conduction can be written in the form
This equation can be substitute by an approximate form
As a consequence, we can obtain the stress components in a simplified form. Accordingly, from Equations (6), (10), and (11), we are led to
The equations of motion, taking into account the Lorentz force
The Lorentz force is given by
The current density vector is parallel to electric intensity vector , thus .
The Ohm’s law (5) after linearization gives
Equations (1), (4), and (34) give
From Equations (2) and (5), we get
From Equations (33) and (34), we obtain
From Equations (27)–(32) and (38), we get
in which we used the notation
For the equation of the equilibrated forces, we obtain
Also, while taking into account Equations (7), (8), and (41), we are led to
Let us define the non-dimensional sizes
where
For dimensionless sizes that are defined in Equation (43), we can write the above basic equations in the following from
by dropping the dashed, for convenience. Here,
is the Peltier coefficient at
and
,
From Equations (44) and (45), we obtain
From Equations (48) and (49), we obtain
5. Numerical Results and Discussion
For numerical computations, following Dhaliwal and Singh [
28], magnesium material was chosen for the purposes of numerical evaluations. All of the units of parameters that were used in the calculation are given in SI units. The constants were taken as
N/m
2,
N/m
2,
W/m·deg,
Kg/m
3,
N/m
2·deg,
J/Kg·deg,
/s,
N/m
2 and
T0 = 298 K.
The voids parameters are m2, N, N/m2, N/m2, N/m2·deg and N/m2s.
The Magnetic field parameters are , , and Col2/Cal·cm·sec.
The comparisons were carried out for , , , , , , and .
Since, we have , and for small values of time we can take , which is a real constant.
The above comparisons have been made in the context of two (G-N) theories of type II and III, in three situations:
- (i)
whether we have an initial stress or not [L* =0 and 105 at M0=0.5 and H0=105];
- (ii)
whether we have a Thomson effect or not [M0 =0 and 0.5 at H0=105 and L*=105];
- (iii)
whether we have some void parameters or not
Case i: In the
Figure 1,
Figure 2,
Figure 3 and
Figure 4, we made the calculations for
, at
. The values of the deformation
e, the values of the temperature
the electromagnetic field
, and the values of the voids function
are graphically represented, for different values of
y in some graphs of two-dimensional space. In these figures, we use the solid lines for the results in the case without initial stress for the (G-N) theory of type II. For the results in the case of the (G-N) theory of type II with initial stress, we have used the large dashes line. In the case without initial stress for the (G-N) theory of type III, we have used the small dashes. Finally, for the results in the case with initial stress for the (G-N) theory of type III, we have used the small dashes line with dot.
Figure 1 depicts the variation of the strain
e versus
y. The magnitude of the strain is found to be large for the G-N theory of type III. It can be seen that the initial stress shows an increasing effect on the magnitude of strain. In
Figure 2, the parameter for the initial stress is decreasing as an effect. Additionally, when the initial stress parameter
is increasing, the value of the temperature
is decreasing.
Figure 3 shows the variation of the induced magnetic field
versus
. The value of strain is found to be large for the theory of G-N of type III. It can be seen that the presence of the initial stress shows a decreasing effect on the magnitude of the induced magnetic field.
Figure 4 expresses the distribution of the change in the volume fraction field
versus
. It was observed that the initial stress has a great effect on the distribution of
.
Case ii: In the
Figure 5,
Figure 6,
Figure 7 and
Figure 8, calculations were made for
at
, and
. The strain
e, the temperature
, the electromagnetic field
h and the voids function
are graphically represented in some of the graphs for different values of y. Here, the solid lines is for results in the G-N theory of type II at
, which gives the classical Fourier’s law of heat conduction, the large dashes line is for results in the type II G-N theory at
, which gives the generalized Fourier’s law of heat conduction. We also use a small dashes line for the results for the type III G-N theory in the case
, while the line with small dashes and the dot is for results for the type III G-N theory for
.
Figure 5 is for the effect of
in the case that it exists and we can see that the value of the deformation
e is increasing when
increases with the corresponding difference.
Figure 6 is for the value of the temperature
, which is decreasing for the parameter
, which is increasing. In
Figure 7, the effect of parameter
exists and the value of the induced magnetic field
h increases when the parameter
increases. In
Figure 8, the value of the voids function
is increasing for the case that the parameter
is increasing.
Case iii:
Figure 9,
Figure 10 and
Figure 11 present the evolution of the physical sizes with regards to the distance
in 2D during
,
,
, and with and without void parameter, in the context of G-N theory of type II and type III G-N theories. In this case of comparison, the solid lines is for results of pore effect
in the type II G-N theory, the small dashes line is also the voids effect, but in the case of the type III G-N theory. We use a large dashes line for results in the G-N theory of type II by neglecting the effect of pores, while the small dashes line with dot is used for results in the type III G-N theory, neglecting the effect of pores. In
Figure 9, we find the repartition of the strain
, and we have a comparison between the values of the strain in the case of the presence of the pores to those in the case of neglecting the voids, in the range
; while, the values are the same for two cases at
.
Figure 10 illustrates the repartition of the temperature
and a comparison between the temperature in the case of presence of pores to those in the case of neglecting the voids, for
y in the range 0 <
y < 1.9. In the case
y > 1.9, the values are the same for two cases.
Figure 11 depicts the repartition of the magnetic field
h, an the values of the magnetic field
h in the case of the presence of pores are compared to those in the case of neglecting the voids, for
y in the range
; in the case
. The values are the same for the two cases.
Figure 12,
Figure 13,
Figure 14 and
Figure 15 are giving 3D surface curves for the physical quantities i.e., the strain
e, the temperature
, the magnetic field
, and the voids fuction
for the thermoelastic theory of electromagnetic bodies with pores, by taking into account the Thomson effect and the effect of the initial stress. The importance of these figures is that they give the dependence of the above physical sizes regarding the vertical component of distance.