This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.
Introduction
Thermoelasticity describes the behavior of elastic bodies under the influence of non-uniform temperature fields. It represents, therefore, a generalization of the theory of elasticity. The constitutive equations, i.e., the equations characterizing the particular material, are the temperature dependent and include an additional relation connecting the heat flux in the body with the local temperature gradient. This relation was known in its simplest form as Fourier’s law.
Biot [1] develop the coupled theory of thermo-elasticity to deal with defeat of the uncoupled theory that mechanical cause has no effect on the temperature field. The theory of generalized thermoelasticity with one relaxation time was introduced by Lord and Shulman [2] which in their replace Fourier’s law with Maxwell-Cattaneo law of heat conduction. Muller [3] first introduced the theory of generalized thermo-elasticity with two relaxation times. Later Green and Naghdi [4], [5] and [6] have proposed three models, labeled as types I, II, and type III. When they are linearized, type I is the same as the classical heat equation whereas the linearized versions of type-II (thermo-elasticity without energy dissipation) and type-III (thermoelasticity with energy dissipation) theories permit the propagation of thermal waves at finite speed.
Tzou [7] introduced two-phase-lag models to both the heat flux vector and the temperature gradient. According to this model, classical Fourier’s $\begin{array}{}
{\rm{law}}\vec q = - K\nabla \vec \theta
\end{array}$ has been replaced by $\begin{array}{}
\vec q(p,t + {\tau _q}) = - K\vec \nabla \theta (p,t + {\tau _T})
\end{array}$, where the temperature gradient $\begin{array}{}
\vec \nabla \theta
\end{array}$ at a point P of a material at time t + τT corresponds to the heat flux vector $\begin{array}{}
\vec q
\end{array}$ at the same point in time t + τq. Here K is the thermal conductivity of the microstructural interactions and is called the phase-lag of the temperature gradient. The other delay time τq is interpreted as the relaxation time due to the fast transient effects of thermal inertial and is called the phase-lag of the heat flux. The case τq = τT = 0 corresponds to classical Fourier’s law. If τq = τ and τT = 0, Tzou refers to the model as single-phase-lag model. Recently the Three-phase-lag, heat conduction equation in which the Fourier law of heat conduction is replaced by an approximation to a modification of the Fourier law with the introduction of three different phases-lags for the heat flux vector, the temperature gradient and the thermal displacement gradient [8]. The stability of the three-phase-lag, the heat conduction is discussed [9]. M. I. Othman and S. M. Said [10] have studied the 2D problem of magneto-thermoelasticity fiber-reinforced medium under temperature dependent properties with three-phase-lag model.
Youssef [11] developed a new theory of generalized thermoelasticity by taking into account the theory of heat conduction in deformable bodies, which depends on two distinct temperatures, the conductive temperature ϕ and the thermodynamic temperatureT. The two-temperature theory involves a material parameter b > 0. The limit b → 0 implies that T → φ and hence the classical theory can be recovered from two temperature theory. Sarhan Y. Atwa [12] has studied the generalized magneto-thermoelasticity with two temperature and initial stress under Green-Naghdi theory.
Diffusion is the process by which atoms move in a material. Atoms are able to move throughout solids because they are not stationary, but execute rapid, small-amplitude vibrations about their equilibrium positions, such vibrations increase with temperature and at any temperature a very small fraction of the atoms has sufficient amplitude to move from one atomic position to an adjacent one. Diffusion can be defined as the mass flow process in which atoms change their positions relative to neighbors in a given phase under the influence of thermal and a gradient, or stress gradient.
Thermodiffusion in the solids is one of a transport process that has great practical importance. Most of research associated with the presence of concentration and temperature gradients has been made with metals and alloys. The first critical review was published in the work of Oriani [13]. Sherif, et al. [14] developed the theory of generalized thermoelastic diffusion that predicts finite speeds of propagation for thermoelastic and diffusive waves. Othman et al. [15] have studied the effect of diffusion on the two-dimensional problem of generalized thermoelasticity with Green-Naghdi theory. Othman et al. [16] have showed the effect of fractional parameter on plane waves of generalized magneto-thermelastic diffusion with reference temperature-dependent elastic medium. Sarhan Y. Atwa [17] has studied the effect of fractional parameter on plane waves of generalized thermoelastic diffusion with temperature-dependent elastic medium.
Basic equations
The governing equations for an isotropic, homogeneous elastic solid with generalized thermoelastic diffusion in the absence of body forces using three-phase-lag model are:
Where, the list of symbols is given in the nomenclature.
Formulation of the problem
We consider an isotropic, homogeneous elastic solid with generalized thermo-elastic diffusion. All quantities are considered as functions of the time variable t and of the coordinates x and y. We consider the normal source acting on the plane surface of generalized thermoelastic half-space under the effect of two temperature, we assume
$$\begin{array}{}
\vec{u}=(u,v,0)
\end{array}$$
The equation of motion in the absence of body force
Where ω is the complex time constant (frequency), i is the imaginary unit, k is the wave number in the x-direction and u*, v*, φ*, Ψ*, q*, T*, C* are the amplitudes of the functions u, v, φ, Ψ, q, T, C.
Using equation (18), equations (13)-(16) become respectively
In this section, the general solutions for displacement, stresses, temperature field, concentration distribution and chemical potential, will be used to yield the response of a half-space subject to a fixed load acting with uniform constant p(x). The final solution in the original domain (x,y,t) is obtained numerically.
The boundary conditions on the plane surface y=0 are
Where p is a function in x and t. Substituting the expressions of the variables considered into the above boundary conditions, we can obtain the following equations satisfied by the parrameters
Solving Eqs. (47)-(50), we get the parameters of Mn (n=1, 2, 3)
Particular case
Neglecting diffusion effect, by taking C = a = b1 = β1 = 0, we obtain the expressions for displacement components, stresses and temperature field in the generalized thermoelastic medium are obtained as
$$\begin{array}{}
\sigma_{xx}=pe^{i(\omega t+kx)}, \frac{\partial \phi}{\partial y}=0, \sigma_{xy}=0 at y=0
\end{array}$$
Substituting from the expressions of the variables considered into the above boundary conditions, we can obtain the following equations satisfied by the parameters
With an aim to illustrate the problem, we will present some numerical results using copper as the thermoelastic material for which we take the following values of the different physical constants.
The numerical result was used to show the distribution of the displacement component u the stress components σyy, σxy the temperature distribution φ, the chemical potential P, and the concentration distribution C on the vertical axis with y at x=0.5, for two theories: Three phase-lag (TPL) model and (G-N III) theory with and without diffusion (W & WD respectively). All curves reaching to zero with the increasing of y. The diffusion effect is appearing on the curves of u, σyy and σxy as in the case of diffusion the curves of the two theories are somewhere far from each and u only start from different values other while in case of absence of diffusion are close from each other, while the effect on φ curve is appearing on the starting values of curves only as in presence of diffusion the two theories start from negative value while in absence of it (TPL) curve start from negative, while (G-N III) start from positive.
Fig. 1 shows that the distribution of u with the two theories in the presence and absence of diffusion starts from positive values. In the presence of diffusion u increases to a maximum value in the range 0<y<0.4 and then decreases to a minimum value in the range 0.5<y<1.4, while in the absence of diffusion u increases to a maximum value in the range 0<y<0.6 and then decreases to a minimum value in the range 0.9<y<1.4.
Fig. 1
Distribution the displacement component u with and without diffusion.
Fig. 2 shows the distribution of the temperature φ with the two theories in the presence and absence of diffusion start from the negative value except the curve belongs to (G-N III) theory without diffusion G-N (WD). The curves belong to TPL theory (W & WD) and G-N (WD) are smooth curves and they fastly reach to zero, while the curve belongs to G-N (D) acts as a wave has a minimum value in the range 0<y<0.9 and a maximum value in the range 0.9<y<2 and then become to decreasing wave until vanishing.
Fig. 2
Distribution the temperature function φ with and without diffusion.
Fig. 3 shows that the stress distribution σyy with the two theories in the presence and absence of diffusion starts from a positive value. In the case of diffusion presence they firstly decrease to a minimum value in the range 0.2<y<1, and then it increases to its maximum value in the range 1<y<1.9 while in the absence of it the minimum value is in the range 0.4<y<1, and the maximum value in the range 1.2<y<1.9.
Fig. 3
Distribution the stress function σyy with and without diffusion.
Fig. 4 shows that the stress distribution σxy in the two theories in the presence and absence of diffusion starts from zero and it satisfy the boundary condition. In presence of diffusion the two theories decreases to a minimum value in the range 0.2<y<0.9, and then it increase to its maximum value in the range 0.9<y<1.9, while in absence of diffusion the two theories decreases to a minimum value in the range 0.3<y<0.9, and then it increase to its maximum value in the range 1.2<y<1.9.
Fig. 4
Distribution the stress function σxy with and without diffusion.
Fig. 5 shows that the chemical potential P in the two theories in the presence of diffusion starts from positive values. P increases to a maximum value in the range 0<y<0.6 and then decreases to a minimum value in the range 0.7<y<1.2.
Fig. 5
Distribution the chemical potential p with and without diffusion.
Fig. 6 shows that the concentration distribution C in the two theories in the presence of diffusion starts from positive values, in TPL theory the curve is smooth and directly reaches to zero, while in (G-N III) theory, it acts as a wave with maximum value in the range 0.1<y<1.2 and a minimum value in the range 1.5<y<2.9.
Fig. 6
Distribution the concentration C with and without diffusion.
3D curves 7-12 are giving 3D surface curves for the physical quantities, under (TPL) model with diffusion. These figures are important to study the dependence of these physical quantities on the vertical component of the distance. All physical quantities are satisfied the boundary conditions.
Conclusions
The analytical solutions based upon normal mode analysis of the thermoelastic problem in solids have been developed.
The method that was used in the present article is applicable to a wide range of problems in thermodynamics and thermoelasticity.
The comparison of the two theories of thermoelasticity, (G-N) and (TPL) theories is carried out.
The physical quantities are very depending on the vertical distance and horizontal distance.
Fig. 7
Distribution the displacement component u against both components of distance in the TPL model with diffusion.
Fig. 8
Distribution the temperature component fagainst both components of distance in the TPL model with diffusion.