# A Two-Dimensional Generalized Electromagnetothermoelastic Diffusion Problem for a Rotating Half-Space

Authors: Jingrui Zhang; Yanyan Li
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

In the context of the theory of generalized thermoelastic diffusion, a two-dimensional generalized electromagnetothermoelastic problem with diffusion for a rotating half-space is investigated. The rotating half-space is placed in an external magnetic field with constant intensity and its bounding surface is subjected to a thermal shock and a chemical potential shock. The problem is formulated based on finite element method and the derived finite element equations are solved directly in time domain. The nondimensional temperature, displacement, stress, chemical potential, concentration, and induced magnetic field are obtained and illustrated graphically. The results show that all the considered variables have a nonzero value only in a bounded region and vanish identically outside this region, which fully demonstrates the nature of the finite speeds of thermoelastic wave and diffusive wave.


---

## Body

## 1. Introduction

Biot [1] proposed the coupled thermoelasticity to amend a defect in uncoupled thermoelasticity that elastic deformation has no effect on temperature. However, this theory shares another defect in uncoupled thermoelasticity in that it predicts infinite speed for heat propagation, which is physically impossible. To overcome such defect, the generalized thermoelastic theories have been developed by Lord and Shulman (L-S) [2] and Green and Lindsay (G-L) [3]. Both theories can characterize the so-called second sound effect; that is, heat propagates in medium with a finite speed. The L-S theory was later extended by Dhaliwal and Sherief [4] to the case of anisotropic media. Based on these generalized thermoelastic theories, many efforts have been devoted to dealing with the generalized dynamic problems. Sherief and Dhaliwal [5] studied a one-dimensional thermal shock problem by the Laplace transform technique and its inverse transform. Dhaliwal and Rokne [6] solved a thermal shock problem of a half-space with its plane boundary either held rigidly fixed or stress-free and an approximate small-time solution was obtained by using the Laplace transform method. Sherief and Anwar [7] considered the thermoelastic problem of a homogeneous isotropic thick plate of infinite extent with heating on a part of the surface by means of state space approach together with Laplace and Fourier integral transforms and their inverse counterparts.Investigation of the propagation of electromagnetothermoelastic waves in a thermoelastic solid has attracted much attention due to its extensive potential applications in diverse fields, such as geophysics for understanding the effect of the Earth’s magnetic field on seismic waves, damping of acoustic waves in a magnetic field, and emissions of electromagnetic radiations from nuclear devices. Sharma and Chand [8] analyzed a one-dimensional transient magnetothermoelastic problem by introducing a potential function. Ezzat et al. [9] researched a two-dimensional electromagnetothermoelastic plane wave problem of a medium of perfect conductivity in terms of normal mode analysis. Sherief and Helmy [10] dealt with a two-dimensional electromagnetothermoelastic problem for a finitely conducting half-space by Laplace and Fourier transforms. Tianhu et al. [11] studied the electromagnetic-thermoelastic interactions in a semi-infinite perfectly conducting solid by hybrid Laplace transform-finite element method. Ezzat and Youssef [12] solved the problem of generalized magnetothermoelasticity in a perfectly conducting medium by means of Laplace and Fourier transform techniques. Sharma and Thakur [13] studied the effect of rotation on Rayleigh-Lamb waves in magnetothermoelastic media. Othman and Song [14] studied the effect of rotation on plane waves of generalized electromagnetothermoviscoelasticity with two relaxation times. Guan [15] studied a two-dimensional of a rotating half-space by using Laplace transform and its numerical inversion. He and Jia [16] studied a two-dimensional of a rotating half-space by using the normal mode analysis. Othman and Song [17] investigated reflection of magnetothermoelastic waves in a rotating medium. Recently, Deswal and Kalkal [18] considered a two-dimensional generalized electromagnetothermoviscoelastic problem for a half-space with diffusion whose surface is subjected to mechanical and thermal loads by introducing potential functions along with the normal modes based on G-L theory.Diffusion can be defined as the random walk of an ensemble of particles from regions of high concentration to regions of lower concentration. There is now a great deal of interest in the study of this phenomenon due to its diverse applications in geophysics and industry. In geology, diffusion principle has been applied to measuring the diffusion coefficients of various cations in minerals which are present in the Earth’s crust. In diffusion bonding, diffusion technique is used to join metallic or nonmetallic materials together. In heat treatment of metals, the surface characteristics of metals, such as wear and corrosion resistance and hardness, can be improved by carburizing through diffusion. In integrated circuit fabrication, diffusion is used to introduce dopants in controlled amounts into the semiconductor substrate. In particular, diffusion is used to form the base and emitter in bipolar transistors, integrated resistors, the source/drain regions in MOS transistors, and dope polysilicon gates in MOS transistors. In the above cases, temperature plays a vital role in the process of diffusion and it is urgent to explore the interactions among diffusion field, strain field, temperature field, and so forth. Normally, the diffusion process is modeled by what is known as Fick’s law, which does not take into consideration the mutual interplay between the introduced substance and the substrate or the effect of temperature on the interplay. Nowacki [19–21] put forward the theory of thermoelastic diffusion in which the coupled thermoelastic model was formulated and infinite speed of propagation of thermoelastic wave was predicted. Recently, Sherief et al. [22] extended this theory associating with L-S model to a generalized thermoelastic diffusion theory that predicts finite speeds of propagation for thermoelastic and diffusive waves. Following this theory, Sherief and Saleh [23] studied a one-dimensional problem of a half-space by using Laplace transform and its numerical inversion. Singh [24] analyzed the reflection problem of SV wave from free surface in an elastic solid. Aouadi [25] examined the thermoelastic diffusion problem for an infinite elastic body with a spherical cavity. Xia et al. [26] worked on the dynamic response of an infinite body with a cylindrical cavity by using finite element method.In the present work, a two-dimensional generalized electromagnetothermoelastic problem with diffusion for a rotating half-space is studied in the context of the theory of the generalized thermoelastic diffusion. The problem is formulated based on finite element method and the derived finite element equations are solved directly in time domain. The variations of the considered variables are obtained and illustrated graphically.


## 2. Basic Equations

The linear electrodynamic equations of slowly moving medium for a homogeneous and perfectly conducting elastic solid are given by Maxwell’s equations as follows:(1)
∇
×
h
=
J
+
ε
0
E
˙
,
(2)
∇
×
E
=
-
μ
0
h
˙
,
(3)
E
=
-
μ
0
(
u
˙
×
H
)
,
(4)
∇
·
h
=
0
,
where H is the applied external magnetic field intensity vector, h is the induced magnetic field vector, E is the induced electric field vector, J is the current density vector, u is the displacement vector, μ
0 and ε
0 are the magnetic permeability and electric permeability, respectively, and ∇ is Hamilton’s operator.In the absence of body force and inner heat source, the generalized electromagnetothermoelastic diffusive governing equations based on the generalized thermoelastic diffusion theory put forth by Sherief et al. [22] can be written as
(5)
σ
i
j
=
2
μ
ε
i
j
+
δ
i
j
(
λ
1
ε
k
k
-
γ
1
θ
-
γ
2
P
)
,
(6)
C
=
γ
2
ε
k
k
+
d
θ
+
n
P
,
(7)
ρ
S
=
γ
1
ε
k
k
+
d
2
θ
+
d
P
,
(8)
σ
i
j
,
j
+
F
i
=
ρ
{
u
¨
i
+
[
Ω
×
(
Ω
×
u
)
]
i
+
(
2
Ω
×
u
˙
)
i
}
,
(9)
F
i
=
μ
0
(
J
×
H
)
i
,
(10)
q
i
+
τ
0
q
˙
i
=
-
κ
i
j
θ
,
j
,
(11)
q
i
,
i
=
-
ρ
T
0
S
˙
,
(12)
η
i
+
τ
η
˙
i
=
-
D
i
j
P
,
j
,
(13)
η
i
,
i
=
-
C
˙
,
(14)
ε
i
j
=
1
2
(
u
i
,
j
+
u
j
,
i
)
,
where
(15)
θ
=
T
-
T
0
,
γ
1
=
β
1
+
a
b
β
2
,
λ
1
=
λ
-
β
2
2
b
γ
2
=
β
2
b
,
d
2
=
ρ
C
E
T
0
+
a
2
b
,
d
=
a
b
,
n
=
1
b
,
β
1
=
(
3
λ
+
2
μ
)
α
t
,
β
2
=
(
3
λ
+
2
μ
)
α
c
.
In the above equations, a superimposed dot denotes the derivative with respect to time, a comma followed by a suffix denotes material derivative, and the summation convention is used. σ
i
j are the components of the stress tensor, ε
i
j are the components of the strain tensor, u
i are the components of displacement vector, κ
i
j are the coefficients of thermal conductivity, S is the entropy density, F
i are the components of Lorentz force, Ω is the angular velocity, η
i is the flow of the diffusing mass vector, q
i are the components of heat flux vector, τ
0 is the thermal relaxation time, τ is the diffusion relaxation time, T is the absolute temperature, T
0 is the initial reference temperature, ρ is the mass density, C
E is the specific heat at constant strain, α
t is the coefficient of linear thermal expansion, α
c is the coefficient of linear diffusion expansion, D
i
j are the coefficients of diffusion, C is the concentration of diffusive material, λ
,
μ are Lame’s constants, P is the chemical potential, “a” is a measure of thermodiffusion effect, and “b” is a measure of diffusive effect.We consider the problem of a homogeneous, isotropic, and perfectly conducting thermoelastic rotating half-space(
x
≥
0
). A magnetic field with constant intensity H
=
(
0,0
,
H
0
) acts parallel to the bounding surface (taken as the direction of the z-axis). At the same time, angular velocity Ω
=
(
0,0
,
Ω
) in half-space goes around the z-axis. Considering rotating effect, the equation of motion is included in the centripetal acceleration related to time and Coriolis acceleration Ω
×
(
Ω
×
u
) item 2
Ω
×
u
˙. The surface of the half-space is subjected at time t
=
0 to a thermal shock and a chemical potential that are functions of y and t. Thus, all the variables will be functions of time t and coordinates x and y. Due to the application of H, this results in an induced magnetic field h and an induced electric field E in the half-space when it undergoes deformation.The displacement components have the form(16)
u
x
=
u
(
x
,
y
,
t
)
,
u
y
=
υ
(
x
,
y
,
t
)
,
u
z
=
0
.
From (14) and (16), we obtain
(17)
e
x
x
=
∂
u
∂
x
,
e
y
y
=
∂
υ
∂
y
,
e
x
y
=
1
2
(
∂
u
∂
y
+
∂
υ
∂
x
)
,
e
x
z
=
e
y
z
=
e
z
z
=
0
.
From (5), the stress components for a homogeneous isotropic solid are given by(18a)
σ
x
x
=
2
μ
∂
u
∂
x
+
λ
1
e
-
γ
1
θ
-
γ
2
P
,
(18b)
σ
y
y
=
2
μ
∂
υ
∂
y
+
λ
1
e
-
γ
1
θ
-
γ
2
P
,
(18c)
σ
x
y
=
μ
(
∂
u
∂
y
+
∂
υ
∂
x
)
,wheree is the cubical dilatation; it takes the expression e
=
e
x
x
+
e
y
y
=
u
,
x
+
υ
,
y.From (1)–(3), we obtain
(19)
E
=
μ
0
H
0
(
-
υ
˙
,
u
˙
,
0
)
,
(20)
h
=
-
H
0
(
0,0
,
e
)
,
(21)
J
=
(
-
H
0
e
,
y
+
ε
0
μ
0
H
0
υ
¨
,
H
0
e
,
x
-
ε
0
μ
0
H
0
u
¨
,
0
)
.
From (9), we get(22a)
F
x
=
μ
0
H
0
2
(
∂
e
∂
x
-
ε
0
μ
0
∂
2
u
∂
t
2
)
,
(22b)
F
y
=
μ
0
H
0
2
(
∂
e
∂
y
-
ε
0
μ
0
∂
2
υ
∂
t
2
)
,
(22c)
F
z
=
0
.It can be noted from (19)–(22c) that the induced electric field, the induced magnetic field, and the Lorentz force are functions of the components of displacement, which implies that the generalized electromagnetothermoelastic problems with diffusion can then be treated as a generalized thermoelastic one with diffusion. Once the components of displacement are obtained, the induced electric field and the induced magnetic field can be calculated from (19) and (20), respectively.Generally speaking, for generalized multifield problems, the involved physical fields, such as electromagnetic field, temperature field, strain field, and diffusion field, would couple with each other, which makes the governing equations of such problems usually too complex to get the solutions by analytical method, so that numerical methods would be powerful tools to solve such problems. One feasible way can be the integral transform techniques. By means of this method, the partial differential governing equations can be converted into ordinary differential equations and solved in transform domain. By applying inverse transform, the solutions of the problem in time domain can then be obtained. However, this method encounters loss of precision believed to be caused by discretization error and truncation error introduced inevitably in the process of numerical inverse Laplace and Fourier transforms, which leads to identifying heat wave front and prediction of the second sound effect, the result of the [27], is thus not so obvious as that the result of [5, 6] showing clear step in the temperature field for small time. An alternative choice to such problems is the hybrid Laplace transform-finite element method presented by Chen and Weng [28, 29]. The same as depicted above, this method also encounters loss of precision and the step of the temperature in the heat wave front is not obvious either. Therefore, the applicability of the integral transform techniques as well as the hybrid Laplace transform-finite element method to generalized thermoelastic problems is limited. To avoid the defects of the above methods, we are inspired to formulate our problem by finite element method and directly solve the derived nonlinear finite element equations in time domain as reported by Tian et al. in [30] in which the obtained results show that this method can achieve a high calculation precision. The diffusion problem done by Xia et al. [26] was just solved by using this method recently.


## 3. Finite Element Formulations

Rewrite (5), (6), and (7) in matrix form as follows:
(23)
{
σ
}
=
[
C
0
]
{
ε
}
-
{
γ
1
}
θ
-
{
γ
2
}
P
,
{
C
}
=
{
γ
2
}
T
{
ε
}
+
d
θ
+
n
P
,
ρ
S
=
{
γ
1
}
T
{
ε
}
+
c
1
θ
+
d
P
.
The generalized heat conduction law and Fick’s law of mass diffusion can be written in matrix form as
(24)
{
q
}
+
τ
0
{
q
˙
}
=
-
[
κ
]
{
θ
′
}
,
{
C
˙
}
+
τ
{
C
¨
}
=
[
D
]
{
P
′′
}
,
where θ
′
=
θ
,
i, P
′′
=
P
,
i
i.According to finite element method, the half-space can be divided into elements and nodal points and any variable considered within an element can be approximated by the values of nodal points together with shape functions. To this end, we introduce two sets of shape functions to approximate the displacement, the temperature, and the chemical potential on the element level:(25)
{
u
}
=
[
N
1
e
]
{
u
e
}
,
θ
=
{
N
2
e
}
T
{
θ
e
}
,
P
=
{
N
2
e
}
T
{
P
e
}
,
where {
u
e
}, {
θ
e
}, and {
P
e
} are the vectors of nodal displacement, temperature, and chemical potential, respectively. [
N
1
e
] and {
N
2
e
} are shape functions; they are
(26)
[
N
1
e
]
=
[
N
1
0
N
2
0
⋯
N
n
0
0
N
1
0
N
2
0
⋯
N
n
]
,
{
N
2
e
}
T
=
{
N
1
N
2
⋯
N
n
}
,
where n denotes the number of nodes in the grid.In terms ofε
i
j
=
(
u
i
,
j
+
u
j
,
i
)
/
2, θ
′
=
θ
,
i, and P
′
=
P
,
i, it yields
(27)
{
ε
}
=
[
B
1
]
{
u
e
}
,
{
θ
′
}
=
[
B
2
]
{
θ
e
}
,
{
P
′
}
=
[
B
2
]
{
P
e
}
,
where [
B
1
] is the strain matrix. In view of the coordinates x and y, [
B
1
] and [
B
2
] are
(28)
[
B
1
]
=
[
∂
N
1
∂
x
0
∂
N
2
∂
x
0
⋯
∂
N
n
∂
x
0
0
∂
N
1
∂
y
0
∂
N
2
∂
y
⋯
0
∂
N
n
∂
y
∂
N
1
∂
y
∂
N
1
∂
x
∂
N
2
∂
y
∂
N
2
∂
x
⋯
∂
N
n
∂
y
∂
N
n
∂
x
]
,
[
B
2
]
=
[
∂
N
1
∂
x
∂
N
2
∂
x
⋯
∂
N
n
∂
x
∂
N
1
∂
y
∂
N
2
∂
y
⋯
∂
N
n
∂
y
]
.
The variational form of (27) is
(29)
δ
{
ε
}
=
[
B
1
]
δ
{
u
e
}
,
δ
{
θ
′
}
=
[
B
2
]
δ
{
θ
e
}
,
δ
{
P
′
}
=
[
B
2
]
δ
{
P
e
}
.In the absence of body force and inner heat source, considering the Lorentz forceF
i, the virtual displacement principle of the generalized electromagnetothermoelastic problems with diffusion can be formulated as
(30)
∫
V
[
δ
{
ε
}
T
{
σ
}
-
ρ
T
0
(
S
˙
+
τ
0
S
¨
)
δ
{
θ
}
+
δ
{
θ
′
}
T
{
q
+
τ
0
q
˙
}
+
δ
{
P
′
}
T
(
η
+
τ
η
˙
)
-
δ
{
P
}
T
(
C
˙
+
τ
C
¨
)
]
d
V
=
∫
V
δ
{
u
}
T
(
{
F
}
-
ρ
[
{
u
¨
}
+
{
Ω
×
(
Ω
×
u
)
}
m
m
m
m
m
m
m
m
m
m
m
l
l
+
{
2
Ω
×
u
˙
}
]
)
d
V
+
∫
A
σ
δ
{
u
}
T
{
T
-
}
d
A
+
∫
A
q
δ
θ
q
-
d
A
+
∫
A
η
δ
P
{
η
-
}
d
A
,
where {
T
-
} represents the traction vector, q
- the heat flux vector, and {
η
-
} the mass flux vector, and the variables with a superimposed bar mean that they are given on surface. A
σ represents the area of the stress tensor, A
q represents the area of the heat flux vector, and A
η represents the area of the mass flux vector.Substituting (23)–(25), (27), and (29) into (30), we obtain
(31)
∫
V
δ
{
ε
}
T
{
σ
}
d
V
=
∫
V
(
δ
{
u
e
}
T
[
B
1
]
T
)
[
[
C
0
]
[
B
1
]
{
u
e
}
-
{
γ
1
}
{
N
2
e
}
T
{
θ
e
}
m
m
m
m
m
m
m
l
m
m
m
m
m
-
{
γ
2
}
{
N
2
e
}
T
{
P
e
}
]
d
V
=
δ
{
u
e
}
T
∫
V
[
B
1
]
T
(
[
[
C
0
]
[
B
1
]
{
u
e
}
-
{
γ
1
}
{
N
2
e
}
T
{
θ
e
}
m
m
m
m
m
m
m
l
m
m
m
m
m
-
{
γ
2
}
{
N
2
e
}
T
{
P
e
}

[
C
0
]
[
B
1
]
]

[
[
C
0
]
[
B
1
]
{
u
e
}
-
{
γ
1
}
{
N
2
e
}
T
{
θ
e
}
)
d
V
=
δ
{
u
e
}
T
(
[
K
m
m
0
e
]
{
u
e
}
-
[
K
m
θ
e
]
{
θ
e
}
-
[
K
m
c
e
]
{
P
e
}
)
-
∫
V
ρ
T
0
(
S
˙
+
τ
0
S
¨
)
δ
{
θ
}
d
V
=
-
∫
V
δ
{
θ
e
}
{
N
2
e
}
T
ρ
T
0
(
S
˙
+
τ
0
S
¨
)
d
V
=
-
∫
V
δ
{
θ
e
}
{
N
2
e
}
T
T
0
[
(
{
γ
1
}
T
[
B
1
]
{
u
˙
e
}
+
c
1
{
N
2
e
}
T
{
θ
˙
e
}
m
m
m
m
m
m
m
m
m
m
m
m
m
l
l
+
d
{
N
2
e
}
T
{
P
˙
e
}
)
+
τ
0
m
m
m
m
m
m
m
m
m
m
m
m
m
l
l
×
(
{
γ
1
}
T
[
B
1
]
{
u
¨
e
}
+
c
1
{
N
2
e
}
T
{
θ
¨
e
}
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
l
+
d
{
N
2
e
}
T
{
P
¨
e
}
)
]
d
V
=
-
δ
{
θ
e
}
(
[
C
θ
m
e
]
{
u
˙
e
}
+
[
C
θ
θ
e
]
{
θ
˙
e
}
+
[
C
θ
c
e
]
{
P
˙
e
}
m
m
m
m
m
m
l
l
+
[
M
θ
m
e
]
{
u
¨
e
}
+
[
M
θ
θ
e
]
{
θ
¨
e
}
+
[
M
θ
c
e
]
{
P
¨
e
}
)
,
∫
V
δ
{
θ
′
}
T
{
q
+
τ
0
q
˙
}
d
V
=
∫
V
(
[
B
2
]
δ
{
θ
e
}
)
T
(
-
[
κ
]
{
θ
′
}
)
d
V
m
m
m
m
m
m
m
m
m
m
m
l
l
=
-
∫
V
δ
{
θ
e
}
[
B
2
]
T
[
κ
]
[
B
2
]
{
θ
e
}
d
V
m
m
m
m
m
m
m
m
m
m
m
l
l
=
-
δ
{
θ
e
}
[
K
θ
θ
e
]
{
θ
e
}
,
-
∫
V
δ
{
P
}
T
(
C
˙
+
τ
C
¨
)
d
V
=
-
∫
V
δ
{
P
}
T
[
(
{
γ
1
}
T
[
B
1
]
{
u
˙
e
}
+
c
1
θ
˙
+
d
P
˙
)
m
m
m
m
m
m
m
m
m
l
+
τ
(
{
γ
1
}
T
[
B
1
]
{
u
¨
e
}
+
c
1
θ
¨
+
d
P
¨
)
]
d
V
=
-
∫
V
δ
{
P
e
}
T
{
N
2
e
}
m
m
m
m
m
×
[
(
{
γ
1
}
T
[
B
1
]
{
u
˙
e
}
+
c
1
{
N
2
e
}
T
{
θ
˙
e
}
m
m
m
m
m
m
m
l
+
d
{
N
2
e
}
T
{
P
˙
e
}
)
m
m
m
m
m
m
l
l
l
l
+
τ
(
{
γ
1
}
T
[
B
1
]
{
u
¨
e
}
+
c
1
{
N
2
e
}
T
{
θ
¨
e
}
m
m
m
m
m
m
m
m
m
m
+
d
{
N
2
e
}
T
{
P
¨
e
}
)
]
d
V
=
-
δ
{
P
e
}
T
(
[
C
c
m
e
]
{
u
˙
e
}
+
[
C
c
θ
e
]
{
θ
˙
e
}
+
[
C
c
c
e
]
{
P
˙
e
}
m
m
m
m
m
m
m
l
l
+
[
M
c
m
e
]
{
u
¨
e
}
+
[
M
c
θ
e
]
{
θ
¨
e
}
+
[
M
c
c
e
]
{
P
¨
e
}
)
∫
V
δ
{
u
}
T
(
{
F
}
-
ρ
[
{
u
¨
}
+
{
Ω
×
(
Ω
×
u
)
}
m
m
m
m
m
m
m
m
m
l
+
{
2
Ω
×
u
˙
}
]
{
F
}
-
ρ
[
{
u
¨
}
+
{
Ω
×
(
Ω
×
u
)
}
)
d
V
=
∫
V
δ
{
u
e
}
T
[
N
1
e
]
T
m
m
m
m
l
l
×
(
+
2
[
0
-
1
-
1
0
]
[
N
1
e
]
{
u
˙
e
}
}
u
0
H
0
2
[
M
]
{
u
e
}
-
ε
0
u
0
2
H
0
2
[
N
1
e
]
{
u
˙
e
}
m
m
m
m
m
m
m
l
-
ρ
{
+
2
[
0
-
1
-
1
0
]
[
N
1
e
]
{
u
˙
e
}
}

[
N
1
e
]
{
u
¨
e
}
-
Ω
2
[
N
1
e
]
{
u
e
}
m
m
m
m
m
m
m
m
m
m
l
+
2
[
0
-
1
-
1
0
]
[
N
1
e
]
{
u
˙
e
}
}
)
d
V
=
δ
{
u
e
}
T
(
[
K
m
m
1
e
]
{
u
e
}
-
[
M
m
m
1
e
]
{
u
¨
e
}
+
[
K
m
m
2
e
]
{
u
e
}
m
m
m
m
m
m
l
l
-
[
M
m
m
2
e
]
{
u
¨
e
}
+
[
C
m
m
e
]
{
u
˙
e
}
)
,
∫
A
σ
δ
{
u
}
T
{
T
-
}
d
A
=
δ
{
u
e
}
T
∫
A
σ
[
N
1
e
]
T
{
T
-
}
d
A
m
m
m
m
m
m
m
m
m
=
δ
{
u
e
}
T
{
T
m
e
}
,
∫
A
q
δ
θ
q
-
d
A
=
δ
{
θ
e
}
T
∫
A
q
{
N
2
e
}
T
q
-
d
A
m
m
m
m
m
l
l
=
δ
{
θ
e
}
T
{
T
θ
e
}
,
∫
A
η
δ
P
{
η
-
}
d
A
=
δ
{
P
e
}
T
∫
A
η
{
N
2
e
}
T
{
η
-
}
d
A
m
m
m
m
m
m
l
l
l
=
δ
{
P
e
}
T
{
T
p
e
}
.
From (31), we arrive at
(32)
[
M
m
m
e
0
0
M
θ
m
e
M
θ
θ
e
M
θ
c
e
M
c
m
e
M
c
θ
e
M
c
c
e
]
{
u
¨
e
θ
¨
e
P
¨
e
}
+
[
C
m
m
e
0
0
C
θ
m
e
C
θ
θ
e
C
θ
c
e
C
c
m
e
C
c
θ
e
C
c
c
e
]
{
u
˙
e
θ
˙
e
P
˙
e
}
+
[
K
m
m
e
-
K
m
θ
e
-
K
m
c
e
0
K
θ
θ
e
0
0
0
K
c
c
e
]
{
u
e
θ
e
P
e
}
=
{
T
m
e
-
T
θ
e
-
T
p
e
}
,
where
(33)
M
c
m
e
=
τ
∫
V
{
N
2
e
}
{
γ
1
}
T
[
B
1
]
d
V
,
M
θ
θ
e
=
∫
V
{
N
2
e
}
T
T
0
τ
0
c
1
{
N
2
e
}
T
d
V
,
M
c
θ
e
=
τ
∫
V
{
N
2
e
}
c
1
{
N
2
e
}
T
d
V
,
M
θ
c
e
=
∫
V
{
N
2
e
}
T
T
0
τ
0
d
{
N
2
e
}
T
d
V
,
M
c
c
e
=
τ
∫
V
{
N
2
e
}
d
{
N
2
e
}
T
d
V
,
C
θ
m
e
=
∫
V
{
N
2
e
}
T
T
0
{
γ
1
}
T
[
B
1
]
d
V
,
C
c
m
e
=
∫
V
{
N
2
e
}
{
γ
1
}
T
[
B
1
]
d
V
,
C
θ
θ
e
=
∫
V
{
N
2
e
}
T
T
0
c
1
{
N
2
e
}
T
d
V
,
C
c
θ
e
=
∫
V
{
N
2
e
}
c
1
{
N
2
e
}
T
d
V
,
C
θ
c
e
=
∫
V
{
N
2
e
}
T
T
0
d
{
N
2
e
}
T
d
V
,
C
c
c
e
=
∫
V
{
N
2
e
}
d
{
N
2
e
}
T
d
V
,
M
θ
m
e
=
∫
V
{
N
2
e
}
T
T
0
τ
0
{
γ
1
}
T
[
B
1
]
d
V
,
K
m
θ
e
=
∫
V
[
B
1
]
T
{
γ
1
}
{
N
2
e
}
T
d
V
,
K
θ
θ
e
=
∫
V
[
B
2
]
T
[
κ
]
[
B
2
]
d
V
,
K
m
c
e
=
∫
V
[
B
1
]
T
{
γ
2
}
{
N
2
e
}
T
d
V
,
K
c
c
e
=
∫
V
[
B
2
]
T
[
D
]
[
B
2
]
d
V
,
T
m
e
=
∫
A
σ
[
N
1
e
]
T
{
T
-
}
d
A
,
T
θ
e
=
∫
A
q
{
N
2
e
}
q
-
d
A
,
T
p
e
=
∫
A
η
{
N
2
e
}
T
{
η
-
}
d
A
,
C
m
m
θ
=
2
∫
V
[
N
1
e
]
T
ρ
[
0
-
1
-
1
0
]
[
N
1
e
]
d
V
,
K
m
m
e
=
K
m
m
0
e
-
K
m
m
1
e
-
K
m
m
2
e
=
∫
V
(
[
B
1
]
T
[
C
0
]
[
B
1
]
-
[
N
1
e
]
T
μ
0
H
0
2
[
D
1
]
-
[
N
1
e
]
T
ρ
Ω
2
[
N
1
e
]
)
d
V
,
M
m
m
e
=
∫
V
[
N
1
e
]
T
(
ε
0
μ
0
2
H
0
2
[
N
1
e
]
+
ρ
[
N
1
e
]
)
d
V
,
[
D
1
]
=
[
∂
2
N
1
∂
x
2
∂
2
N
1
∂
x
∂
y
⋯
∂
2
N
n
∂
x
2
∂
2
N
n
∂
x
∂
y
∂
2
N
1
∂
x
∂
y
∂
2
N
1
∂
y
2
⋯
∂
2
N
n
∂
x
∂
y
∂
2
N
n
∂
y
2
]
.
Once the initial conditions and the boundary conditions are specified, the finite element equation in (32) can be solved directly in time domain. In the process of numerical calculation and finite element solution, the space domain and time domain are discrete. In the calculation, because the surface of O
A is subjected to a thermal shock and a chemical potential shock, this part of the unit is divided in a more detailed way; the entire model is divided into 1535 units and 3222 nodes; similarly, the initial time step is set to t
=
3
×
1
0
-
7; the variable threshold is set at t
=
1
×
1
0
-
7, which ensures the accuracy and convergence and also saves a lot of calculated time.


## 4. Numerical Results and Discussions

The schematic of the considered half-space as well as the applied loads on its bounding surface is shown in Figure1(a). The bounding surface is assumed to be traction-free, and the thermal shock and the chemical potential shock applied on the bounding Surface have, respectively, the following form:
(34)
θ
=
θ
0
H
(
t
)
H
(
L
-
|
y
|
)
,
P
=
P
0
H
(
t
)
H
(
L
-
|
y
|
)
,
where H
(
·
) is the Heaviside unit step function and θ
0 and P
0 are constants.Schematic of the rotating half-space.
(a)
(b)Assume that the rotating half-space is initially at rest, so that the initial conditions are(35)
u
=
υ
=
θ
=
P
=
0
at
t
=
0
,
u
˙
=
υ
˙
=
θ
˙
=
P
˙
=
0
at
t
=
0
.
Due to the symmetries of geometrical shape and boundary conditions, the problem can be treated as a plane strain problem and only half of the half-space needs to be considered. The model for simulation is shown in Figure 1(b), where O
A
B
C outlines the region for implementing the simulation and O
D represents the region within which the thermal shock and the chemical potential shock are applied.The half-space is taken to be copper material and the material properties are(36)
λ
=
7.76
×
1
0
10
kg
/

(
m
s
2
)
,
μ
=
3.86
×
1
0
10
kg
/

(
m
s
2
)
,
ρ
=
8954
kg
/
m
3
,
κ
=
386
W
/

(
m
K
)
,
D
=
8.5
×
1
0
-
9
kg
s
/
m
3
,
α
c
=
1.98
×
1
0
-
4
m
3
/
kg
,
α
t
=
1.78
×
1
0
-
5
K
-
1
,
C
E
=
383.1
J
/

(
kg
K
)
,
a
=
1.2
×
1
0
4
m
2
/
(
s
2
K
)
,
b
=
9
×
1
0
5
m
5
/
(
kg
s
2
)
.
To simplify the simulation, we introduce the following nondimensional variables:
(37)
x
*
=
c
1
η
1
x
,
y
*
=
c
1
η
1
y
,
u
*
=
c
1
η
1
u
,
υ
*
=
c
1
η
1
υ
,
t
*
=
c
1
2
η
1
t
,
τ
*
=
c
1
2
η
1
τ
,
τ
0
*
=
c
1
2
η
1
τ
0
,
θ
*
=
β
1
θ
λ
+
2
μ
,
C
*
=
β
2
C
λ
+
2
μ
,
σ
i
j
*
=
σ
i
j
λ
+
2
μ
,
P
*
=
P
β
2
,
h
*
=
h
H
0
,
η
1
=
ρ
C
E
κ
,
c
1
2
=
λ
+
2
μ
ρ
,
Ω
*
=
Ω
c
1
2
η
,
i
,
j
=
1,2
.In calculation, we specifyτ
0
=
0.02, τ
=
0.2, T
0
=
293
K, θ
0
=
1, P
0
=
1, Ω
=
0.01, and L
=
O
D
=
0.2; the dimensions along x-axis and y-axis are O
A
=
3.0 and O
C
=
3.0, respectively.The calculations are carried out for three values of nondimensional times, namely,t
=
0.05, t
=
0.1, and t
=
0.15. The nondimensional temperature, displacement, stress, chemical potential, concentration, and induced magnetic field are illustrated in Figures 2–12, respectively, dropping the asterisk at the upper right corner of the nondimensional variables for convenience. It should be pointed out that wave reflection from any edges is excluded in the simulation.Figure 2
Nondimensional temperature distribution alongO
A.Figure 3
Nondimensional temperature distribution alongO
C.Figure 4
Nondimensional horizontal displacement distribution alongO
A.Figure 5
Nondimensional horizontal displacement distribution alongO
C.Figure 6
Nondimensional vertical displacement distribution alongO
C.Figure 7
Nondimensional stress distribution alongO
A.Figure 8
Nondimensional chemical potential distribution alongO
A.Figure 9
Nondimensional concentration distribution alongO
A.Figure 10
Nondimensional chemical potential distribution alongO
C.Figure 11
Nondimensional concentration distribution alongO
C.Figure 12
Nondimensional induced magnetic field distribution alongO
A.Figures2 and 3 show the distributions of the nondimensional temperature along O
A and O
C, respectively. In Figure 2, when the time t is given, the distance of the heat propagation in the x direction should be x
=
v
h
t, where v
h is nondimensional heat wave velocity. When τ
=
0.02, we can achieve v
h
=
0.07. The heat propagation in the x direction at the time t
=
0.5, t
=
0.1 is x
=
0.35, x
=
0.7 respectively. From Figure 2 a distinct temperature step on thermal wave front distribution on O
A can be readily seen, but it becomes indistinct along with the passage of time. In Figure 3, within 0
≤
y
≤
0.2, the temperature keeps constant all along, which is consistent with the thermal boundary condition along O
D. As shown in Figures 2 and 3, the temperature increases with the passage of time.Figure4 shows the distributions of the nondimensional horizontal displacement along O
A. Due to the thermal shock, the parts of the half-space near the bounding surface expand toward the unconstrained direction, thus yielding negative displacement. It can be found that different parts on O
A undergo different deformations. Some undergo expansion and some undergo compression, while the rest remain undisturbed, resulting in the displacement shift from negative to positive gradually. With the passage of time, the expansion parts enlarge and move inside dynamically, making the negative-positive region transform dynamically. We should be aware that the vertical displacement for O
A is always zero because of the symmetries.Figures5 and 6 show the distributions of the nondimensional horizontal and vertical displacements on O
C, respectively. As seen from Figure 5, the horizontal displacement on O
C is negative, which means that O
C undergoes thermal expansion deformation and moves toward the unconstrained direction. As shown in Figure 6, the vertical displacement along O
C is positive. This can be interpreted as follows. Due to the symmetries, the vertical displacement of point “o” is always zero. This implies that point “o” can not be allowed to move up and down, which prevents all the other points on O
C from moving downwards, thus leading to positive displacement. The vertical displacement on O
C firstly goes up and then goes down. It can also be seen from Figures 5 and 6 that the magnitudes of the displacement increase with the passage of time.Figure7 shows the distributions of nondimensional stress σ
x
x along O
A. Due to the symmetries, the other two components of stress, namely, σ
y
y and σ
x
y, are always zero along O
A. It can be observed that σ
x
x is all negative, which is known as compressive stress.Figures8, 9, 10, and 11 show the distributions of nondimensional chemical potential and concentration along O
A and O
C, respectively. It can be observed that the speed of diffusive wave is larger than that of thermoelastic wave, which can be deduced by comparing the distance of diffusive wave traversing across the half-space with that of thermoelastic wave along the same direction for a given time.Figure12 shows the distributions of the nondimensional induced magnetic field. Due to the mutual interaction between the applied external magnetic field and the elastic deformation, this results in an induced magnetic field in the half-space. As shown in Figure 12, the magnitude of h increases with the passage of time.It can be readily seen from Figures2–12 that all the considered variables have a nonzero value only in a bounded region and the value vanishes outside this region, which is totally dominated by the nature of the finite speeds of thermoelastic wave and diffusive wave.From Figures13, 14, 15, 16, 17, 18, 19, 20, 21, and 22, it can be readily seen that rotation acts to decrease the magnitude of the real part of displacement, stress, and induced magnetic field and not to affect the magnitude of temperature, chemical potential, and concentration.Figure 13
Nondimensional induced magnetic field distribution alongO
A.Figure 14
Nondimensional temperature distribution alongO
A.Figure 15
Nondimensional temperature distribution alongO
C.Figure 16
Nondimensional horizontal displacement distribution alongO
A.Figure 17
Nondimensional horizontal displacement distribution alongO
C.Figure 18
Nondimensional stress distribution alongO
A.Figure 19
Nondimensional chemical potential distribution alongO
A.Figure 20
Nondimensional chemical potential distribution alongO
C.Figure 21
Nondimensional concentration distribution alongO
A.Figure 22
Nondimensional concentration distribution alongO
C.


## 5. Concluding Remarks

A two-dimensional generalized electromagnetothermoelastic problem with diffusion for a rotating half-space is studied in the context of the theory of generalized thermoelastic diffusion by means of finite element method. The nondimensional temperature, displacement, stress, chemical potential, concentration, and induced magnetic field are obtained. The results show that(
1
) all the considered variables have a nonzero value only in a bounded region and vanish identically outside this region, which is governed by the nature of the finite speeds of thermoelastic wave and diffusive wave; (
2
) the speed of diffusive wave is larger than that of thermoelastic wave; (
3
) rotation acts to decrease the magnitude of the real part of displacement, stress, and induced magnetic field and not to affect the magnitude of temperature, chemical potential, and concentration.



---

*Source: 964218-2014-04-09.xml*