Abstract
The governing equations for generalized thermoelasticity of a mixture of an elastic solid and a Newtonian fluid are formulated in the context of Lord-Shulman and Green-Lindsay theories of generalized thermoelasticity. These equations are solved to show the existence of three coupled longitudinal waves and two coupled transverse waves, which are dispersive in nature. Reflection from a thermally insulated stress-free surface is considered for incidence of coupled longitudinal wave. The speeds and reflection coefficients of plane waves are computed numerically for a particular model.
1. Introduction
The soil consists of an assemblage of particles with different sizes and shapes which form a skeleton whose voids are filled by water and air or gas. The word โsoil,โ therefore, implies a mixture of assorted mineral grains with various fluids. The first continuum theory of mixtures was proposed by Truesdell [1] and Truesdell and Toupin [2] in terms of the kinematic and thermodynamic variables associated with each constituent of the mixture. Another theory presented by Green and Naghdi [3] uses the kinematic and thermodynamic variables associated with the mixture as a whole. The classical theory of mixtures was discussed by Atkin and Craine [4] and Bowen [5]. Bedford and Drumheller [6] formulated the theories of immiscible and structured mixtures. Various other theories of mixtures were also developed by Muller [7], Dunwoody [8], Krishnaswamy and Batra [9], Iesan [10โ12], Rajagopal and Tao [13], and Pompei and Scalia [14].
Many engineering materials, as well as soils, rocks, granular materials, sand, and underground water mixtures may be modeled more realistically by means of micropolar continua. Eringen [15] developed a theory of micropolar mixture of porous media (nonreacting mixture of micropolar elastic solid and a micropolar viscous fluid at a single temperature) to include the rotational degrees of freedom. Eringen [15] also obtained the field equations of a mixture of an elastic solid and a Newtonian fluid as a special case by dropping the micropolar effects. In the present paper, the governing equations for generalized thermoelasticity of mixture of an elastic solid and a Newtonian fluid are formulated in context of Lord-Shulman [16] and Green-Lindsay [17] theories of generalized thermoelasticity. These equations are solved to show the existence of various coupled plane waves. The reflection of these plane waves is studied at thermally insulated stress-free surface and the reflection coefficients are computed and shown graphically with the angle of incidence. Wave propagation in such a media may find its applications in consolidation and oil exploration problems.
2. Equations of Motion
We consider a binary mixture of an elastic solid and a Newtonian fluid at the same temperature where no chemical reaction takes place between the two species of the mixture. According to the last section given in Eringen [15], the equations of motion in a linear isotropic thermoelastic mixture of an elastic solid and Newtonian fluid in absence of external loads are written as where are displacement vectors in solid and fluid phase, respectively, and is the temperature variable. The coefficients are the thermoelastic constants for isotropic elastic solids, are fluid viscosities, is classical Fourier constant, is the heat generation due to velocity difference, is the momentum generation coefficient due to the velocity difference, is specific heat at constant deformation, is an ambient temperature, , is Helmholtzโs free energy function, are the densities of mixture, solid, and fluid in the natural state, superposed dot indicates the temporal derivative, and other symbols have their usual meanings.
The use of symbol , in (2.3), makes these fundamental equations possible for the two different theories of the generalized thermoelasticity. For the L-S (Lord-Shulman) theory and for G-L (Green-Lindsay) theory and . The thermal relaxations and satisfy the inequality for the G-L theory only.
By introducing the scalar potentials and vector potentials through Helmholtz representation of vector field, we can write
Using (2.4) into (2.1) to (2.3), we obtain the following five equations, that is, three coupled equations in and two coupled equations in ; where, .
3. Wave Propagation
In this section, the plane wave propagation in an infinite thermally conducting mixture of elastic solid and Newtonian fluid is studied. In order to solve (2.5), we consider the following form of plane waves propagating in the positive direction of a unit vector where are the constant complex scalar wave amplitudes, , is the position vector, is the phase velocity in the direction of , k is the wave number, and is angular frequency. Using the values of , and from (3.1) into (2.5), we obtain a set of three homogeneous equations in three unknown amplitudes . After eliminating these three unknowns, we get the following cubic equation where and the expressions for , , are given in Appendix A. The three roots of (3.2) correspond to complex speeds of three dispersive coupled longitudinal waves in a thermally conducting mixture of elastic solid and Newtonian fluid for both L-S and G-L theories.
Equation (2.6) can be solved by considering the following form of vector potentials: where are constant complex vector wave amplitudes. With the help of (3.3), (2.6) leads to homogeneous vector equations in and . After elimination of and , the following quadratic equation in is obtained: Equation (3.4) gives two roots , which correspond to complex speeds of two dispersive coupled transverse waves in a thermally conducting mixture of elastic solid and Newtonian fluid. Hence, the analytical solution indicates the existence of five dispersive plane waves in two-dimensional model of a thermally conducting mixture of an elastic solid and a Newtonian fluid for both L-S and G-L theories.
4. Particular Cases
(i)In absence of fluid, (3.2) reduces to the following quadratic equation: which gives complex speeds of two coupled longitudinal waves for both L-S and G-L theories whereas (3.4) reduces to a linear equation which gives the speed of transverse wave. (ii)In absence of thermal effects, the cubic equation (3.2) reduces to the following quadratic equation: The roots of (4.2) correspond to the complex phase speeds of two coupled longitudinal waves in a mixture of an elastic solid and a Newtonian fluid.
5. Reflection from Traction-Free Thermally Insulated Surface
We consider the half-space of the medium with x-axis along the thermally insulated stress-free surface and z-axis into the medium. For incident coupled longitudinal waves or coupled transverse waves, there will be five reflected waves into the medium as shown in Figure 1. The required boundary conditions at stress-free thermally insulated surface are Here,
The appropriate potentials required to satisfy the above boundary conditions are Here and are the amplitude of incident coupled longitudinal I (CL I), reflected coupled longitudinal I (CL I), reflected coupled longitudinal II (CL II), reflected coupled longitudinal III (CL III), incident coupled transverse I (CT I), reflected coupled transverse I (CT I) and reflected coupled transverse II (CT II) waves, respectively, where for incident coupled longitudinal I (CL I) wave, for incident coupled transverse I (CT I) wave, and and the relations between angles of incidence and reflection are given as and, , where and , ( = 1,2,โฆ,5) are real phase speeds and attenuations of coupled longitudinal and transverse waves.
Using the appropriate potentials given by (5.3) with the required boundary conditions (5.1), the following nonhomogeneous system of five equations is obtained where , and are given in Appendix B.
6. Numerical Results and Discussion
In this section, the complex absolutes of phase speeds and reflection coefficients are computed with the following physical constants at ,, ,, ,, , , , , .
The densities of solid, fluid, and mixture are related as , where is porosity of the mixture.
Numerical computations of phase speeds and reflection coefficients are restricted to the particular case of Lord-Shulman theory only. The phase speeds of the coupled waves are shown graphically in Figure 2 for the range of frequency. It is observed that the phase speeds of all coupled longitudinal waves increase with the increase in frequency. In absence of thermal effects, the phase speed of CL I wave decreases at each value of frequency, whereas the phase speed of CL II increases. The wave CL III will disappear in absence of thermal parameters. The phase speeds of coupled transverse waves (CT I and CT II) remain unaffected by thermal effects.
(a)
(b)
(c)
(d)
Reflection coefficients are computed for incident CL I wave only. The reflection coefficients of coupled longitudinal waves are shown graphically in Figure 3 against the angle of incidence. The reflection coefficients of CL I and CL II waves change with the increase in angle of incidence. The comparison of solid and dotted curves shows the effects of thermal parameters on these coefficients. The reflected CL III wave will disappear in absence of thermal effects. The reflection coefficients of CT I and CT II waves also change with the increase in angle of incidence. These waves are also affected by thermal parameters as shown in Figure 4.
(a)
(b)
(c)
(a)
(b)
7. Conclusions
The generalized thermoelasticity of a mixture of an elastic solid and a Newtonian fluid is developed in accordance with Lord-Shulman and Green-Lindsay theories. The solutions of governing equations suggest that there will exist three coupled longitudinal waves and two coupled transverse waves in the present model. From numerical results, it is observed that the phase speeds of coupled waves increase with the increase in frequency. The presence of thermal parameters affects the phase speeds of coupled longitudinal waves only. Reflection coefficients of all reflected coupled waves are affected significantly by thermal disturbances in the material.
Appendices
A. Expressions for , , and
and, , , , , , .
B. Expressions for , and
and for incident CL I wave (), for incident CT I wave (), Here, the amplitude ratios , and correspond to reflected CL I wave, CL II wave, CL III wave, CT I wave, and CT II wave, respectively.