Propagation of Plane Waves in a Thermally Conducting Mixture

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.


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 ISRN Applied Mathematics 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.

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 u S , u F are displacement vectors in solid and fluid phase, respectively, and T is the temperature variable.The coefficients λ S , μ S , β 0 are the thermoelastic constants for isotropic elastic solids, λ F , μ F are fluid viscosities, K * is classical Fourier constant, s is the heat generation due to velocity difference, ξ is the momentum generation coefficient due to the velocity difference, C 0 is specific heat at constant deformation, T 0 is an ambient temperature, δ − ∂ψ/∂ρ F 0 , ψ is Helmholtz's free energy function, ρ, ρ S , ρ F 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 t 1 0, Ω 1 and for G-L Green-Lindsay theory t 1 > 0 and Ω 0. The thermal relaxations t 0 and t 1 satisfy the inequality t 1 ≥ t 0 ≥ 0 for the G-L theory only.

ISRN Applied Mathematics 3
By introducing the scalar potentials q S , q F and vector potentials U S , U F through Helmholtz representation of vector field, we can write

2.4
Using 2.4 into 2.1 to 2.3 , we obtain the following five equations, that is, three coupled equations in q S , q F , T and two coupled equations in U S , U F ;

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 n q S , q F , T q S , q F , T exp ιk n • r − V t , 3.1 where q S , q F , T are the constant complex scalar wave amplitudes, ι √ −1, r is the position vector, V is the phase velocity in the direction of n, k is the wave number, and ω kV is angular frequency.Using the values of q S , q F , and T from 3.1 into 2.5 , we obtain a set of three homogeneous equations in three unknown amplitudes q S , q F , T .After eliminating these three unknowns, we get the following cubic equation where Γ V 2 and the expressions for A, B, C are given in Appendix A. The three roots V 2  1 , V 2 2 , V 2 3 of 3.2 correspond to complex speeds of three dispersive coupled longitudinal 4 ISRN Applied Mathematics 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 A S , A F are constant complex vector wave amplitudes.With the help of 3.3 , 2.6 leads to homogeneous vector equations in A S and A F .After elimination of A S and A F , the following quadratic equation in V 2 is obtained:

3.4
Equation 3.4 gives two roots V 2 4 , V 2 5 , 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.

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:

4.2
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.

Reflection from Traction-Free Thermally Insulated Surface
We consider the half-space of the medium with x-axis along the thermally insulated stressfree 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,

5.2
The appropriate potentials required to satisfy the above boundary conditions are Σ η j A j exp ιk j sin θ j x cos θ j z − ιω j t , j 4, 5 .

5.3
Here A 0 , A 1 , A 2 , A 3 , B 0 , A B 0 0 for incident coupled longitudinal I CL I wave, A 0 0 for incident coupled transverse I CT I wave, and and, V j −1 c j −1 −ιω −1 q j , where c j and q j , j 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

Numerical Results and Discussion
In this section, the complex absolutes of phase speeds and reflection coefficients are computed with the following physical constants at T 0 293 • K, t 0 0.005 s, t 1 0.006 s, β 0 0.0005, δ 0.0001, n 0.15.
The densities of solid, fluid, and mixture are related as ρ 1 − n ρ S nρ F , where n 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 5 Hz ≤ ω ≤ 50 Hz 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.
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.

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.

5 Figure 1 :
Figure 1: Geometry of the problem showing incident and reflected waves.

Figure 3 :
Figure 3: Thermal effects on coefficients of reflected coupled longitudinal waves for incident CL I wave.

Figure 4 :
Figure 4: Thermal effects on coefficients of reflected coupled transverse waves for incident CL I wave.