Next Article in Journal
On the Structure of Finite Groupoids and Their Representations
Next Article in Special Issue
Multiparametric Analytical Solution for the Eigenvalue Problem of FGM Porous Circular Plates
Previous Article in Journal
Stability Analysis of Darcy-Forchheimer Flow of Casson Type Nanofluid Over an Exponential Sheet: Investigation of Critical Points
Previous Article in Special Issue
New Exact Solutions of the Generalized Benjamin–Bona–Mahony Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Effect of Thomson and Initial Stress in a Thermo-Porous Elastic Solid under G-N Electromagnetic Theory

by
Elsayed M. Abd-Elaziz
1,
Marin Marin
2,* and
Mohamed I. A. Othman
3
1
Zagazig Higher Institute of Engineering & Technology, Ministry of Higher Education, P.O. Box 44519, Zagazig, Egypt
2
Department of Mathematics and Computer Science, Transilvania University of Brasov, Brasov 500036, Romania
3
Department of Mathematics, Faculty of Science, Zagazig University, P.O. Box 44519, Zagazig, Egypt
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(3), 413; https://doi.org/10.3390/sym11030413
Submission received: 22 January 2019 / Revised: 6 March 2019 / Accepted: 14 March 2019 / Published: 20 March 2019
(This article belongs to the Special Issue Symmetry in Applied Continuous Mechanics)

Abstract

:
The present work investigated the effect of Thomson and initial stress in a thermo-porous elastic solid under G-N electromagnetic theory. The Thomson coefficient affects the heat condition equation. A constant Thomson coefficient, instead of traditionally a constant Seebeck coefficient, is assumed. The charge density of the induced electric current is taken as a function of time. A normal mode method is proposed to analyze the problem and to obtain numerical solutions. The results that were obtained for all physical sizes are graphically illustrated and we offer a comparison between the type II G-N theory and the G-N theory of type III, both in the present case and in the absence of specific parameters, as initial stress, pores and the Thomson effect. Some particular cases are also discussed in the context of the problem. The results indicate that the effect of initial stress, Thomson coefficient effect, and magnetic field are very pronounced.

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 T 0 , in the reference state, and the half space ( y 0 ) . The motion referred to a rectangular Cartesian system of coordinates ( x , y , z ) with origins in the surface ( z = 0 ) . 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 u = ( u , v , 0 ) . 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 H = ( 0 , 0 , H 3 ) , having a constant intensity, which acts parallel to the direction of the Z-axis.
It is known that a magnetic field of the form H = ( 0 , 0 , H 0 + h ( x , y , t ) ) produces an induced electric field of components E = ( E 1 , E 2 , 0 ) , 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:
×    h = J + D ˙ ,
×    E = B ˙ ,
.    B = 0 ,   .    D = ρ e ,
B = μ 0 ( H + h ) ,   D = ε 0 E .
The modified Ohm’s law for a medium with finite conductivity supplements the above system of coupled equations, namely
J = σ 0 ( E + μ 0 u ˙ × H ) ,
where μ 0 is the magnetic permeability, B is the magnetic displacement vector, ε 0 is the electric permeability, J is the current density vector, ρ e is the charge density, D is the electric displacement vector, and E is the induced electric field vector.
For an isotropic and homogeneous thermoelastic body having pores, the constitutive equations receive the following form:
σ i j = 2 μ e i j + ( λ e r r + λ 0    ϕ β T ) δ i j L ( δ i j + m i j ) ,
h i = α ϕ , i ,
g = λ 0 e r r ξ 1 φ + m T ,
ρ η * = β e r r + a 0 T + m ϕ .
The strain-displacement relation is
e i j = 1 2 ( u i , j + u j , i ) .
The tensor of rotation has the components:
m i j = 1 2 ( u j , i u i , j ) ,   i , j = 1 , 2 , 3 .
In Green-Naghdi (G-N) theories we take into account the Thomson effect, so that the Fourier’s law becomes
q i = [ k T , i + k T ˙ , i ] + M J i ,
which gives
q i , i = [ k T , i i + k T ˙ , i i ] + M J i , i .
If we take into account Equations (1) and (3), then from Equation (13), we deduce
q i , i = [ k T , i i + k T ˙ , i i ] + M ρ ˙ e .
where T is the temperature above the reference temperature T 0 is chosen so that | ( T T 0 ) / T 0 | < 1 , λ , μ are the counterparts of Lame’ constants, t is the time, σ i j are the components of the stress tensor, h , i is the equilibrated stress vector, ψ is the equilibrated inertia, g is the intrinsic equilibrated body force, α , λ 0 , ξ 1 , ω 0 , m are constants of material that are due to the presence of the pores, β = ( 3 λ + 2 μ ) α t , such that α t is the coefficient of thermal expansion, δ i j is the Kronecker delta, ρ is the mass density, C E is the specific heat at the constant strain, k is the thermal conductivity, η is entropy per unit mass, k * is a constant, and q i j 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
ρ ˙ e + ( ρ e v i ) = 0 ,
where the velocity of the charges has the components v i .
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
ρ ˙ e + ρ e ( v i ) = 0 .
We will assume that the charges have the speed of components v i , which are proportional to the components of the velocity for particles u ˙ i , so that we can write
v i = p 0 u ˙ i ,
which gives
v i = p 0 u ˙ i = p 0 e ˙ ,
where p 0 is a positive constant (non-dimensional).
If we take into account Equation (18), from Equation (16), we are led to
ρ ˙ e = ρ e p 0 e ˙ ,
which gives
d ρ e ρ = p 0 d e .
Hence, we obtain
ρ e = ρ e 0 exp ( p 0 e ) ρ e 0 ( 1 p 0 e ) ,
where ρ e 0 is the charge density when the strain vanishes.
Then, we obtain
ρ ˙ e = ρ e 0 ( 1 p 0 e ) p 0 e ˙ .
While taking into account the Equation (22), from Equation (14) we deduce that the Fourier’s law, in its generalized form, receives the form:
q i , i = [ k T , i i + k T ˙ , i i ] + M ρ e 0 ( 1 p 0 e ) p 0 e ˙ .
In the case of null heat supply, the balance energy becomes
ρ T 0 η ˙ * =    q i , i .
Taking into account Equations (9) and (23), from Equation (24), we deduce that the equation of heat conduction can be written in the form
k T , i i + k * T ˙ , i i m T 0 ϕ ˙ = ρ    C e T ¨ + β    T 0 u ¨ i , i + M ρ e 0 ( 1 p 0 e ) p 0 e ˙ .
This equation can be substitute by an approximate form
k T , i i + k * T ˙ , i i m T 0 ϕ ˙ = ρ    C e T ¨ + β T 0 u ¨ i , i + M ρ e 0 p 0 e ˙ .
As a consequence, we can obtain the stress components in a simplified form. Accordingly, from Equations (6), (10), and (11), we are led to
σ x x = λ [ u x + v y ] + 2 μ u x + λ 0 ϕ β T L ,
σ y y = λ [ u x + v y ] + 2 μ v y + λ 0 ϕ β T L ,
σ z z = λ [ u x + v y ] + λ 0 ϕ β T L ,
σ x y = ( μ + L 2 ) u y + ( μ L 2 ) v x ,
σ y x = ( μ + L 2 ) v x + ( μ L 2 ) u y , σ x z = σ y z = 0 .
The equations of motion, taking into account the Lorentz force
σ j i , j + F i = ρ u i , t t .
The Lorentz force is given by
F i = μ 0 ( J × H ) i .
The current density vector J is parallel to electric intensity vector E , thus J = ( J 1 , J 2 , 0 ) .
The Ohm’s law (5) after linearization gives
J 1 = σ 0 ( E 1 + μ 0 H 0 v ˙ ) , J 2 = σ 0 ( E 2 μ 0 H 0 u ˙ ) .
Equations (1), (4), and (34) give
h y = σ 0 ( E 1 + μ 0 H 0 v t ) + ε 0 E 1 t ,
h x = σ 0 ( E 1 μ 0 H 0 u t ) ε 0 E 2 t .
From Equations (2) and (5), we get
E 1 y E 2 x = μ 0 h t .
From Equations (33) and (34), we obtain
F 1 = σ 0 μ 0 H 0 ( E 2 μ 0 H 0 u t ) , F 2 = σ 0 μ 0 H 0 ( E 1 + μ 0 H 0 v t ) , F 3 = 0 .
From Equations (27)–(32) and (38), we get
( μ L 2 ) 2 u + ( λ + μ + L 2 ) e x + b ϕ x β T x + σ 0 μ 0 H 0 ( E 2 μ 0 H 0 u t ) = ρ 2 u t 2
( μ L 2 ) 2 v + ( λ + μ + L 2 ) e y + b ϕ y β T y σ 0 μ 0 H 0 ( E 1 + μ 0 H 0 v t ) = ρ 2 v t 2 ,
in which we used the notation e = u x + v y .
For the equation of the equilibrated forces, we obtain
ρ ψ ϕ , t t = h i , i + g .
Also, while taking into account Equations (7), (8), and (41), we are led to
α ϕ , j j λ 0 u i , i ξ 1 ϕ ω 0 ϕ , t + m T = ρ ψ ϕ ¨ .
Let us define the non-dimensional sizes
( x i , u i ) = w 1 c 1 ( x i , u i ) , ( t , t 0 ) = w 1 ( t , t 0 ) , ϕ = ψ w 1 2 c 1 2 ϕ , σ i j = σ i j μ , p 1 = p 1 μ , L = L μ ,
θ = T T 0 T 0 , h = w 1 σ 0 H 0 μ 0 c 1 2 h , E i = w 1 σ 0 H 0 2 μ 0 c 1 2 E i ,
where w 1 = ρ c e c 1 2 K , c 1 2 = λ + 2 μ ρ .
For dimensionless sizes that are defined in Equation (43), we can write the above basic equations in the following from
a 1 2 u + a 2 e x + a 3 ϕ x a 4 θ x + a 5 ( a 6 E 2 u t ) = a 7 2 u t 2 ,
a 1 2 v + a 2 e y + a 3 ϕ y a 4 θ y a 5 ( a 6 E 1 + v t ) = a 7 2 v t 2 ,
( 2 a 8 a 9 t a 10 2 t 2 ) ϕ a 11 e + a 12 θ = 0 ,
[ ( k + k w 1 t ) 2 a 13 2 t 2 ] θ + a 14 ϕ t β c 1 2 e t 2 M 0 e t = 0 ,
h y = a 15 E 1 + a 16 E 1 t + v t ,
h x = a 15 E 2 a 16 E 2 t + u t ,
E 1 y E 2 x = h t ,
by dropping the dashed, for convenience. Here, M 0 = M ρ e 0 p 0 c 1 2 w 1 T 0 , is the Peltier coefficient at T 0 and a 1 = μ L μ 2 , a 2 = λ + μ + L μ 2 , a 3 = b c 1 2 ψ w 1 2 , a 4 = β T 0 , a 5 = μ 0 2 H 0 2 σ 0 c 1 2 w 1 , a 6 = μ 0 σ 0 c 1 2 w 1 , a 7 = ρ c 1 2 , a 8 = ξ 1 c 1 2 α w 1 2 , a 9 = w 0 c 1 2 α w 1 , a 12 = m ψ T 0 α , a 13 = ρ C e c 1 2 , a 14 = m c 1 4 ψ w 1 3 , a 15 = σ 0 μ 0 c 1 2 w 1 , a 16 = μ 0 ε 0 c 1 2 , a 10 = ρ ψ c 1 2 α , a 11 = λ 0 ψ α .
From Equations (44) and (45), we obtain
[ ( a 1 + a 2 ) 2 a 7 t a 5 2 t 2 ] e + a 3 2 ϕ a 4 2 θ a 6 a 5 h t = 0 .
From Equations (48) and (49), we obtain
( 2 a 15 t a 16 2 t 2 ) h e t = 0 .

3. The Solution of the Problem

3.1. Decomposition by Normal Mode Analysis

Using the normal mode analysis, we can decompose the solution of the above physical parameters in the following form
[ e , θ , ϕ , h , σ i j ] ( x , y , t ) = [ e , θ , ϕ , h , σ i j ] ( y ) exp [ i ( a x ω t ) ] ,
in which e , ϕ , θ , h , σ i j are the amplitudes of the respective fields, i = 1 , ω is the frequency, and a is the wave number.
By taking into account Equation (53) in Equations (46), (47), (51), and (52), we are led to
( r 1 D 2 r 2 ) e + a 3 ( D 2 a 2 ) ϕ * a 4 ( D 2 a 2 ) θ * + r 3 h = 0 ,
( D 2 r 4 ) ϕ a 11 e * + a 12 θ = 0 ,
( r 5 D 2 r 6 ) θ + r 7 ϕ * + r 8 e * = 0 ,
( D 2 r 9 ) h + i ω e * = 0 ,
where r 1 = a 1 + a 2 , r 2 = r a 2 i ω a 5 a 7 ω 2 , r 3 = i ω a 6 a 5 , r 4 = a 2 + a 8 i ω a 9 a 10 ω 2 , r 5 = k i ω k w , r 6 = r 5 a 2 a 13 ω 2 , r 7 = i ω a 14 , r 8 = β ω 2 c 1 2 i ω M 0 , and r 9 = a 2 i ω a 15 a 10 ω 2 .
Eliminating (55), we get the following ordinary Equation (52) between e ( y ) ϕ ( y ) , θ ( y ) , the differential Equation satisfied by h ( y ) :
( D 8 N 1 D 6 + N 2 D 4 N 3 D 2 + N 4 ) { h ( y ) , ϕ ( y ) , e ( y ) , θ ( y ) } = 0 .
We can write the Equation (58) in a decomposed form, as follows
( D 2 k 1 2 ) ( D 2 k 2 2 ) ( D 2 k 3 2 ) ( D 2 k 4 2 ) { h ( y ) , ϕ ( y ) , e ( y ) , θ ( y ) } = 0 ,
where, k n 2 ( n = 1 , 2 , 3 , 4 ) are roots of the characteristic equation of Equation (58) and s 1 = r 5 r 9 + r 4 r 5 + r 6 , s 2 = r 6 r 9 + r 4 r 5 r 9 a 12 r 7 + r 6 r 4 , s 3 = r 4 r 6 r 9 + a 12 r 7 r 9 , s 4 = r 5 r 6 a 12 r 8 + r 6 , s 5 = a 11 r 6 r 9 a 12 r 8 r 9 , s 6 = r 8 r 9 + r 4 r 8 a 11 r 3 , s 7 = r 8 r 4 r 9 a 11 r 7 r 9 , s 8 = r 6 + r 4 r 5 , s 9 = r 4 r 6 a 12 r 7 , s 10 = 1 r 1 r 5 , N 1 = s 10 ( r 1 s 1 + r 2 r 5 a 3 a 11 r 5 + a 4 r 8 ) , N 2 = s 10 ( r 1 s 2 + r 2 s 1 a 3 a 11 r 5 a 2 a 3 s 4 + a 4 s 6 + a 4 r 8 a 2 i ω r 3 r 5 ) , N 3 = s 10 ( r 1 s 3 + r 2 s 2 a 3 s 5 a 3 s 4 a 2 + a 4 s 7 + a 4 s 6 a 2 i ω r 3 r 8 ) , and N 4 = s 10 ( r 2 s 3 a 3 s 5 a 2 + a 4 s 7 a 2 i ω r 3 r 9 ) .
The general solutions of the Equation (59), bound at y ,
e ( x , y , t ) = n = 1 4 R n exp [ k n y + i ( a x ω t ) ] ,
h ( x , y , t ) = n = 1 4 H 1 n R n exp [ k n y + i ( a x ω t ) ] ,
θ ( x , y , t ) = n = 1 4 H 2 n R n exp [ k n y + i ( a x ω t ) ] ,
ϕ ( x , y , t ) = n = 1 4 H 3 n R n exp [ k n y + i ( a x ω t ) ] ,
where
H 1 n = i ω k n 2 s 9 , H 2 n = r 8 ( k n 2 r 4 ) r 7 a 11 ( k n 2 r 4 ) ( r 5 k n 2 r 6 ) r 7 a 12 , H 3 n = a 11 a 12 H 2 n k n 2 r 4 .

3.2. Boundary Conditions

In the following, we will consider some boundary conditions on the surface of Equation y = 0 , which will help us to determine the above constants R1, R2, R3 and R4.

3.2.1. The mechanical boundary condition

The mechanical boundary condition that the bounding plane to the surface y = 0 has zero strain, so we have
e ( x , 0 , t ) = 0

3.2.2. The Boundary Restriction of Heat

We assume that the boundary surface of the body is subject to a thermal shock described by the function
θ ( x , 0 , t ) = θ 0 exp ( i ( a x ω t ) ) ,
where θ 0 is constant.

3.2.3. Voids Conditions

ϕ y = 0 .

3.2.4. The Boundary Restriction for Electromagnetic Field

On the surface of the half-space h ( x , 0 , t ) = h , we consider that the electromagnetic field intensity is a continuous function. Here, the intensity of the magnetic field in free space is h .
We now assume there is no magnetic or electric field in the free space, that is,
h ( x , 0 , t ) = h = 0 .
In order to obtain the constants R1, R2, R3 and R4, we will use the dimensionless size θ 0 = θ 0 T 0 and the expressions of the variables into the boundary restrictions imposed above. Additionally, we will use the normal mode analysis in order to obtain the system of equations
n = 1 4 R n = 0 ,
n = 1 4 H 2 n R n = θ 0 ,
n = 1 4 k n H 3 n R n = 0 ,
n = 1 4 H 1 n R n = 0 .
after suppressing the primes.
After applying the inverse of matrix method for the above equations, we get the values of the constants R n ( n = 1 , 2 , 3 , 4 ) , hence; we obtain the expressions of strain, magnetic intensity, temperature distribution and the change in the volume fraction field for the generalized thermoelastic medium with voids.

4. Special Cases

4.1. Pores Neglect

First, we will neglect the presence of the voids, that is, we have ( α = λ 0 = ξ 1 = ω 0 = m = ψ = 0 ) .
While putting ( α = λ 0 = ξ 1 = ω 0 = m = ψ = 0 ) in Equations (54)–(57), we get:
( r 1 D 2 r 2 ) e a 4 ( D 2 a 2 ) θ * + r 3 h = 0 ,
( r 5 D 2 r 6 ) θ + r 8 e * = 0 ,
( D 2 r 9 ) h + i ω e * = 0 .
Eliminating e , θ , and h * among Equations (72)–(74), we obtain the following sixth order differential Equation, which is satisfied by e ( y ) , θ ( y ) and h ( y )
[ D 6 B 1 D 4 + B 2 D 2 B 3 ] { e ( y ) , θ ( y ) , h * ( y ) } = 0 ,
where s 10 = 1 r 1 r 5 , B 2 = s 10 [ r 1 r 6 r 9 + r 2 r 5 r 9 + r 2 r 6 r 8 a 4 ( r 9 + a 2 ) + i ω r 3 r 5 ] , B 1 = s 10 ( r 1 r 9 r 5 + r 1 r 6 + r 2 r 5 r 8 a 4 ) , B 3 = s 10 ( r 2 r 5 r 9 r 8 a 4 r 9 a 2 + i ω r 3 r 6 ) .
The solutions of Equation (75) are
( e , h , θ ) ( y ) = n = 1 3 ( 1 , H 4 n , H 5 n ) R n exp ( α n y ) ,
where α n 2 ( n = 1 , 2 , 3 ) are the roots of the characteristic equation of Equation (75) and H 4 n = i ω α n 2 r 9 , H 5 n = r 8 r 5 α n 2 r 6 .
The expressions for the strain, the induced magnetic field, and the temperature field in the generalized initially stressed the electro-magneto-thermoelastic half-space solid with voids are:
( e , h , θ ) ( y ) = n = 1 3 ( 1 , H 4 n , H 5 n ) R n exp ( α n y + i ( a x ω t ) ) .
We wish to determine the above coefficients R n , ( n = 1 , 2 , 3 ) .
To this end, we will keep in mind the boundary conditions in Equations (64), (65), and (67), and we will use the method of inverse of the matrix, as following:
( R 1 R 2 R 3 ) = ( 1 1 1 H 51 H 52 H 53 H 41 H 42 H 43 ) 1 ( 0 θ 0 0 ) .

4.2. Neglecting the Initial Stress

By taking ( L = 0 ) in the governing equations, the corresponding expressions of the physical variables can be obtained without initial stress.

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 λ = 2.17 × 10 10 N/m2, μ = 3.278 × 10 10 N/m2, k = 1.7 × 10 2 W/m·deg, ρ = 1.74 × 10 3 Kg/m3, β = 2.68 × 10 6 N/m2·deg, C e = 1.04 × 10 3 J/Kg·deg, ω 1 * = 3.58 × 10 11 /s, α t = 1.78 × 10 5 N/m2 and T0 = 298 K.
The voids parameters are ψ = 1.753 × 10 15 m2, α = 3.688 × 10 5 N, ξ 1 = 1.475 × 10 10 N/m2, λ 0 = 1.13849 × 10 10 N/m2, m = 2 × 10 6 N/m2·deg and ω 0 = 0.0787 × 10 3 N/m2s.
The Magnetic field parameters are ε 0 = 10 9 / ( 36 π ) F / M , μ 0 = 4 π × 10 7 H / M , H 0 = 10 5 A / M and σ 0 = 9.36 × 10 5 Col2/Cal·cm·sec.
The comparisons were carried out for a = 0.7    m , θ 0 = 0.1    k , ω = χ 0 + i χ 1 , χ 0 = 2 rad / s , χ 1 = 0.09 rad / s , x = 0.2    m , and 0 y 4 m .
Since, we have exp ( ω t ) = [ cos ( χ 1 t ) + i sin ( χ 1 t ) ] exp ( χ 0 t ) , and for small values of time we can take ω = χ 0 , 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 [ M 0 = 0.5 , H 0 = 10 5 a n d L = 10 5 ] .
Case i: In the Figure 1, Figure 2, Figure 3 and Figure 4, we made the calculations for t = 0.2 , at x = 0.2 . The values of the deformation e, the values of the temperature θ , the electromagnetic field h , 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 L * is increasing, the value of the temperature θ is decreasing. Figure 3 shows the variation of the induced magnetic field h versus y . 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 y . 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 M 0 = 0.0 , 0.5 at H 0 = 10 5 , and t = 0.2 . 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 M 0 = 0.0 , 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 M 0 = 0.5 , 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 M 0 = 0.0 , while the line with small dashes and the dot is for results for the type III G-N theory for M 0 = 0.5 . Figure 5 is for the effect of M 0 in the case that it exists and we can see that the value of the deformation e is increasing when M 0 increases with the corresponding difference. Figure 6 is for the value of the temperature θ , which is decreasing for the parameter M 0 , which is increasing. In Figure 7, the effect of parameter M 0 exists and the value of the induced magnetic field h increases when the parameter M 0 increases. In Figure 8, the value of the voids function ϕ is increasing for the case that the parameter M 0 is increasing.
Case iii: Figure 9, Figure 10 and Figure 11 present the evolution of the physical sizes with regards to the distance y in 2D during M 0 = 0.5 , H 0 = 10 5 , L = 10 5 , 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 1.9 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 e , 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 0 y 1.7 ; while, the values are the same for two cases at y 1.7 . 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 0 y 3 ; in the case y 3 . 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 h , 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.

6. Conclusions

The results concluded from the above analysis can be summarized, as follows:
(1)
We have derived the field equations of homogeneous, isotropic, electro-magneto-thermo- porous elastic half-plane with the Thomson effect and initial stress.
(2)
The analytical solutions that are based upon normal mode analysis for the thermoelastic problem in solids have been developed and utilized.
(3)
The presence of initial stress, void parameters, and Thomson effect play significant roles in all of the physical quantities.
(4)
The value of all physical quantities converges to zero with the increase in distance y and all of the functions are continuous.
(5)
The deformation of a body depends on the nature of the applied forces and Thomson effect, as well as the type of boundary conditions.

Author Contributions

All three authors conceived the framework and structured the whole paper, checked the results of the paper and completed the revision of the article. The authors have equally contributed to the elaboration of this manuscript. All authors have read and approved the final manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the anonymous referee for the comments that helped us improve this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Green, A.E.; Naghdi, P.M. A re-examination of the basic postulates of thermo-mechanics. Proc. R. Soc. Lond. A 1991, 432, 171–194. [Google Scholar] [CrossRef]
  2. Green, A.E.; Naghdi, P.M. On undamped heat wave in elastic solids. J. Therm. Stresses 1992, 15, 253–264. [Google Scholar] [CrossRef]
  3. Green, A.E.; Naghdi, P.M. Thermoelasticity without energy dissipation. J. Elast. 1993, 31, 189–209. [Google Scholar] [CrossRef]
  4. Sharma, J.N.; Chouhan, R.N. On the problems of body forces and heat sources in thermoelasticity without energy dissipation, Indian. J. Pure Appl. Math. 1999, 30, 595–610. [Google Scholar]
  5. Liang, M.; Liu, Y.; Xiao, B.; Yang, S.; Wang, Z.; Han, H. An analytical model for the transverse permeability of gas diffusion layer with electrical double layer effects in proton exchange membrane fuel cells. Int. J. Hydrogen Energy 2018, 43, 17880–17888. [Google Scholar] [CrossRef]
  6. Xiao, B.; Zhang, X.; Wang, W.; Long, G.; Chen, H.; Kang, H.; Ren, W. Afractal model for water flow through unsaturated porous rocks. Fractals 2018, 26, 1840015. [Google Scholar] [CrossRef]
  7. Hassan, M.; Marin, M.; Alsharif, A.; Ellahi, R. Convective heat transfer flow of nanofluid in a porous medium over wavy surface. Phys. Lett. A Gen. At. Solid State Phys. 2018, 382, 2749–2753. [Google Scholar] [CrossRef]
  8. Long, G.; Liu, S.; Xu, G.; Won, S.W. A perforation-erosion model for hydraulic-fracturing applications. SPE Prod. Oper. 2018, 33, 770–783. [Google Scholar] [CrossRef]
  9. Xiao, B.; Chen, H.; Xiao, S.; Cai, J. Research on relative permeability of nanofibers with capillary pressure effect by means of fractal-monte Carlo technique. J. Nanosci. Nanotechnol. 2017, 17, 6811–6817. [Google Scholar] [CrossRef]
  10. Xiao, B.; Wang, W.; Fan, J.; Chen, H.; Hu, X.; Zhao, D.; Zhang, X.; Ren, W. Optimization of the fractal-like architecture of porous fibrous materials related to permeability, diffusivity and thermal conductivity. Fractals 2017, 25, 1750030. [Google Scholar] [CrossRef]
  11. Nunziato, J.W.; Cowin, S.C. A non-linear theory of elastic materials with voids. Arch. Rat. Mech. Anal. 1979, 72, 175–201. [Google Scholar] [CrossRef]
  12. Cowin, S.C.; Nunziato, J.W. Linear theory of elastic materials with voids. J. Elast. 1983, 13, 125–147. [Google Scholar] [CrossRef]
  13. Puri, P.; Cowin, S.C. Plane waves in linear elastic materials with voids. J. Elast. 1985, 15, 167–183. [Google Scholar] [CrossRef]
  14. Iesan, D. A theory of thermoelastic materials with voids. Acta Mech. 1986, 60, 67–89. [Google Scholar] [CrossRef]
  15. Cicco, S.D.; Diaco, M. A theory of thermoelastic material with voids without energy dissipation. J. Therm. Stress. 2002, 25, 493–503. [Google Scholar] [CrossRef]
  16. Othman, M.I.A.; Edeeb, E.R.M. The effect of rotation on thermoelastic medium with voids and temperature dependent under three theories. J. Eng. Mech. 2018, 144, 04018003. [Google Scholar] [CrossRef]
  17. Montanaro, A. On singular surfaces in isotropic linear thermoelasticity with initial stress. J. Acoust. Soc. Am. 1999, 106, 1586–1588. [Google Scholar] [CrossRef]
  18. Othman, M.I.A.; Abd-Elaziz, E.M. The effect of thermal loading due to laser pulse on generalized thermoelastic medium with voids in dual phase lag model. J. Therm. Stress. 2015, 38, 1068–1082. [Google Scholar] [CrossRef]
  19. Marin, M. Cesaro means in thermoelasticity of dipolar bodies. Acta Mechanica 1997, 122, 155–168. [Google Scholar] [CrossRef]
  20. Marin, M.; Oechsner, A. The effect of a dipolar structure on the Holder stability in Green-Naghdi thermoelasticity. Contin Mech Thermodyn 2017, 29, 1365–1374. [Google Scholar] [CrossRef]
  21. Nowinski, J. Theory of Thermoelasticity with Applications; Sijthoff and Noordhoff International Publishers: Alphen Aan Den Rijn, The Netherlands, 1978. [Google Scholar]
  22. Chadwick, P. Progress in Solid Mechanics; Hill, R., Sneddon, I.N., Eds.; North Holland Publishing Company: North Holland, Amsterdam, 1960. [Google Scholar]
  23. Hussain, F.; Ellahi, R.; Zeeshan, A. Mathematical models of electromagnetohydrodynamic multiphase flows synthesis with nanosized hafnium particles. Appl. Sci. 2018, 8, 275. [Google Scholar] [CrossRef]
  24. Bhatti, M.M.; Zeeshan, A.; Ellahi, R.; Ijaz, N. Heat and mass transfer of two-phase flow with Electric double layer effects induced due to peristaltic propulsion in the presence of transverse magnetic field. J. Mol. Liq. 2017, 230, 237–246. [Google Scholar] [CrossRef]
  25. Othman, M.I.A.; Abd-Elaziz, E.M. Plane waves in a magneto-thermoelastic solids with voids and microtemperatures due to hall current and rotation. Results Phys. 2017, 7, 4253–4263. [Google Scholar] [CrossRef]
  26. Marin, M.; Craciun, E.M.; Pop, N. Consideration on mixed initial-boundary value problems for micro-polarporous bodies. Dyn. Syst. Appl. 2016, 25, 175–196. [Google Scholar]
  27. Marin, M.; Craciun, E.M. Uniqueness results for a boundary value problem in dipolar thermoelasticity to model composite materials. Compos. Part B Eng. 2017, 126, 27–37. [Google Scholar] [CrossRef]
  28. Dhaliwal, R.S.; Singh, A. Dynamic Coupled Thermoelasticity; Hindustan Publ. Corp.: New Delhi, India, 1980. [Google Scholar]
Figure 1. The strain e distribution at M 0 = 0.5 and H 0 = 10 5 .
Figure 1. The strain e distribution at M 0 = 0.5 and H 0 = 10 5 .
Symmetry 11 00413 g001
Figure 2. The temperature θ distribution at M 0 = 0.5 and H 0 = 10 5 .
Figure 2. The temperature θ distribution at M 0 = 0.5 and H 0 = 10 5 .
Symmetry 11 00413 g002
Figure 3. The induced magnetic field h distribution at M 0 = 0.5 , and H 0 = 10 5 .
Figure 3. The induced magnetic field h distribution at M 0 = 0.5 , and H 0 = 10 5 .
Symmetry 11 00413 g003
Figure 4. The change in the volume fraction field ϕ distribution at M 0 = 0.5 and H 0 = 10 5 .
Figure 4. The change in the volume fraction field ϕ distribution at M 0 = 0.5 and H 0 = 10 5 .
Symmetry 11 00413 g004
Figure 5. The strain e distribution at H 0 = 10 5 and L * = 10 5 .
Figure 5. The strain e distribution at H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g005
Figure 6. The temperature θ distribution at H 0 = 10 5 and L * = 10 5 .
Figure 6. The temperature θ distribution at H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g006
Figure 7. The induced magnetic field h distribution at H 0 = 10 5 and L * = 10 5 .
Figure 7. The induced magnetic field h distribution at H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g007
Figure 8. The change in the volume fraction field ϕ distribution at H 0 = 10 5 and L * = 10 5 .
Figure 8. The change in the volume fraction field ϕ distribution at H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g008
Figure 9. The strain e distribution at H 0 = 10 5 , M 0 = 0.5 and L * = 10 5 .
Figure 9. The strain e distribution at H 0 = 10 5 , M 0 = 0.5 and L * = 10 5 .
Symmetry 11 00413 g009
Figure 10. The temperature θ distribution at H 0 = 10 5 , M 0 = 0.5 , and L * = 10 5 .
Figure 10. The temperature θ distribution at H 0 = 10 5 , M 0 = 0.5 , and L * = 10 5 .
Symmetry 11 00413 g010
Figure 11. The induced magnetic field h distribution at H 0 = 10 5 , M 0 = 0.5 , and L * = 10 5 .
Figure 11. The induced magnetic field h distribution at H 0 = 10 5 , M 0 = 0.5 , and L * = 10 5 .
Symmetry 11 00413 g011
Figure 12. Three-dimensional (3D) curve distribution of the strain e versus the distances at: M 0 = 0.5 , H 0 = 10 5 , and L * = 10 5 .
Figure 12. Three-dimensional (3D) curve distribution of the strain e versus the distances at: M 0 = 0.5 , H 0 = 10 5 , and L * = 10 5 .
Symmetry 11 00413 g012
Figure 13. 3D Curve distribution of the temperature θ versus the distances, at: M 0 = 0.5 , H 0 = 10 5 and L * = 10 5 .
Figure 13. 3D Curve distribution of the temperature θ versus the distances, at: M 0 = 0.5 , H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g013
Figure 14. 3D Curve distribution of the induced magnetic field h versus the distances, at: M 0 = 0.5 , H 0 = 10 5 and L * = 10 5 .
Figure 14. 3D Curve distribution of the induced magnetic field h versus the distances, at: M 0 = 0.5 , H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g014
Figure 15. 3D Curve distribution of the change in the volume fraction field ϕ versus the distances, at: M 0 = 0.5 , H 0 = 10 5 and L * = 10 5 .
Figure 15. 3D Curve distribution of the change in the volume fraction field ϕ versus the distances, at: M 0 = 0.5 , H 0 = 10 5 and L * = 10 5 .
Symmetry 11 00413 g015

Share and Cite

MDPI and ACS Style

Abd-Elaziz, E.M.; Marin, M.; Othman, M.I.A. On the Effect of Thomson and Initial Stress in a Thermo-Porous Elastic Solid under G-N Electromagnetic Theory. Symmetry 2019, 11, 413. https://doi.org/10.3390/sym11030413

AMA Style

Abd-Elaziz EM, Marin M, Othman MIA. On the Effect of Thomson and Initial Stress in a Thermo-Porous Elastic Solid under G-N Electromagnetic Theory. Symmetry. 2019; 11(3):413. https://doi.org/10.3390/sym11030413

Chicago/Turabian Style

Abd-Elaziz, Elsayed M., Marin Marin, and Mohamed I. A. Othman. 2019. "On the Effect of Thomson and Initial Stress in a Thermo-Porous Elastic Solid under G-N Electromagnetic Theory" Symmetry 11, no. 3: 413. https://doi.org/10.3390/sym11030413

APA Style

Abd-Elaziz, E. M., Marin, M., & Othman, M. I. A. (2019). On the Effect of Thomson and Initial Stress in a Thermo-Porous Elastic Solid under G-N Electromagnetic Theory. Symmetry, 11(3), 413. https://doi.org/10.3390/sym11030413

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop