The present investigation is concerned with the propagation of waves at an imperfect boundary of heat conducting elastic solid and micropolar fluid media. The amplitude ratios of various reflected and transmitted waves are obtained in a closed form due to incidence of longitudinal wave (P-wave), thermal wave (T-wave), and transverse wave (SV-wave). The variation of various amplitude ratios with angle of incidence is obtained for normal force stiffness, transverse force stiffness, thermal contact conductance, and perfect bonding. Numerical results are shown graphically to depict the effect of stiffness and thermal relaxation times on resulting quantities. Some particular cases are also deduced in the present investigation.
1. Introduction
The fluids in which coupling between the spin of each fluid particle and the microscopic velocity field is taken into account are termed as micropolar fluids. They represent fluids consisting of rigid, randomly oriented, or spherical particles suspended in a viscous medium, where the deformation of fluid particles is ignored. The theory of microfluids was introduced by Eringen [1]. A microfluid possesses three gyration vector fields in addition to its classical translatory degrees of freedom. As a subclass of these fluids, Eringen introduced the micropolar fluids [2] to describe the physical systems, which do not fall in the realm of viscous fluids. In micropolar fluids the local fluid elements were allowed to undergo only rigid rotations without stretch. These fluids support couple stress, the body couples, and asymmetric stress tensor and possess a rotational field, which is independent of the velocity of fluid. A large class of fluids such as anisotropic fluids, liquid crystals with rigid molecules, magnetic fluids, cloud with dust, muddy fluids, biological fluids, and dirty fluids (dusty air, snow) can be modeled more realistically as micropolar fluids. The importance of micropolar fluids in industrial applications has motivated many researchers to extend the study in numerous ways to explain various physical effects.
Ariman et al. [3, 4] studied microcontinuum fluid mechanics and its applications and Říha [5] investigated the theory of heat-conducting micropolar fluid with microtemperature. Eringen and Kafadar [6] discussed polar field theories, Brulin [7] studied linear micropolar media, Gorla [8] investigated combined forced and free convection in the boundary layer flow of a micropolar fluid on a continuous moving vertical cylinder, and Eringen [9] studied the theory of thermo-microstretch fluids and bubbly liquid. Aydemir and Venart [10] investigated flow of a thermomicropolar fluid with stretch. Ciarletta [11] established two sorts of spatial decay estimates for describing the spatial behaviour of the solutions for the flow of a heat-conducting micropolar fluid in a semi-infinite cylindrical pipe. The lateral surface is taken to be thermal insulated and the adherence of the fluid to the lateral boundary is assumed. A time-dependent velocity and angular velocity profile are prescribed together with the temperature field at the finite end of the pipe. Hsia and Cheng [12] studied longitudinal plane waves propagation in elastic micropolar porous media. Propagation of transverse waves in elastic micropolar porous semispaces is discussed by Hsia et al. [13].
Different researchers studied the problems of reflection and transmission of plane waves at an interface of micropolar elastic half-spaces (Tomar and Gogna [14–16] and Kumar et al. [17, 18]). Singh and Tomar [19] discussed the longitudinal waves at an interface of micropolar fluid/micropolar solid half-spaces. Sun et al. [20] studied propagation characteristics of longitudinal displacement wave in micropolar fluid with micropolar elastic plate. Xu et al. [21] discussed reflection and transmission characteristics of coupled wave through micropolar elastic solid interlayer in micropolar fluid. Dang et al. [22] investigated propagation characteristics of coupled wave through micropolar fluid interlayer in micropolar elastic solid. Fu and Wei [23] investigated wave propagation through the imperfect interface between two micropolar solids. Sharma and Marin [24] studied reflection and transmission of waves from imperfect boundary between two heat conducting micropolar thermoelastic solids. Vishwakarma et al. [25] studied torsional wave propagation in an earth’s crustal layer under the influence of imperfect interface.
The purpose of the present study is to study the problem of reflection and transmission of plane waves at an imperfect interface between heat conducting elastic solid and micropolar fluid media. Effects of stiffness and thermal relaxation times on the amplitude ratios for incidence of various plane waves, for example, longitudinal wave (P-wave), thermal wave (T-wave), and transverse wave (SV-wave), are depicted numerically and shown graphically with angle of incidence.
2. Basic Equations
The field equations in a homogeneous and isotropic elastic medium in the context of generalized theories of thermoelasticity, without body forces and heat sources, are given by Lord and Shulman [26] and Green and Lindsay [27] as
(1)(λ+2μ)∇(∇·u→)-μ(∇×(∇×u→))-ν(1+τ1∂∂t)∇T=ρ∂2u→∂t2,K*ΔT=ρc*(∂T∂t+τ0∂2T∂t2)+νT0(∂∂t+η0τ0∂2∂t2)(∇·u→),
and the constitutive relations are
(2)tij=λur,rδij+μ(ui,j+uj,i)-ν(T+τ1∂T∂t)δij,hhhhhhhhhhhhhhhhhhhhhhhhhi,j,r=1,2,3,
where ν=(3λ+2μ)αT and the meaning of symbols is defined in the list at the end of the paper. The thermal relaxation times τ0 and τ1 satisfy the inequalities τ0≥τ1≥0 for G-L theory only and Δ=∂2/∂x12+∂2/∂x32 is the Laplacian operator. For Lord Shulman (L-S) theory η0=1, τ1=0 and for Green-Lindsay (G-L) theory η0=0, τ1>0.
Following Ciarletta [11], the field equations and the constitutive relations for heat conducting micropolar fluids without body forces, body couples, and heat sources are given by
(3)D1v→+(λf+μf)∇(∇·v→)+Kf(∇×Ψ→)-b∇Tf-c0∇ϕ*f=0,D2Ψ→+(αf+βf)∇(∇·Ψ→)+Kf(∇×v→)=0,(4)K1*ΔTf-bT0f(∇·v→)=ρfaT0f∂Tf∂t,(5)ρf∂ϕ*f∂t=∇·v→,
where
(6)D1=(μf+Kf)Δ-ρf∂∂t,D2=γfΔ-I∂∂t-2Kf,
such that superscript f denotes physical quantities and material constants related to the fluid and the constitutive relations are
(7)tijf=-pδij+σijf,p=bTf+c0ϕ*f,σijf=λfγrrδij+(μf+Kf)γij+μfγji,mijf=αfυrrδij+βfυji+γfυij,
where γij=vj,i+εjirΨr, υij=Ψj,i, b=(3λf+2μf+Kf)αTf and symbols are defined in the list at the end of the paper.
3. Formulation of the Problem
An imperfect interface of a homogeneous, isotropic generalized thermoelastic half-space (medium M1) in contact with heat conducting micropolar fluid half-space (medium M2) is considered. The rectangular Cartesian coordinate system Ox1x2x3 having origin on the surface x3=0 seperating the two media is taken. Let us take the x1-axis along the interface between two half-spaces, namely, M1 (0<x3<∞) and M2 (-∞<x3<0), in such a way that x3-axis is pointing vertically downward into the medium M1. The geometry of the problem is shown in Figure 1.
Geometry of the problem.
For two dimensional problem in x1x3-plane, we take the displacement vector u→, velocity vector v→, and microrotation velocity vector Ψ→ as
(8)u→=(u1(x1,x3),0,u3(x1,x3)),v→=(v1(x1,x3),0,v3(x1,x3)),Ψ→=(0,Ψ2(x1,x3),0).
The following nondimensional quantities are defined:
(9)xi′=ω*xic1,ui′=ρω*c1νT0ui,vi′=ρc1νT0vi,ψ2′=ρc12ω*νT0ψ2,(t′,τ0′,τ1′)=(ω*t,ω*τ0,ω*τ1),(T′,Tf′)=(TT0,TfT0),(tij′,tijf′)=1νT0(tij,tijf),mijf′=ω*c1νT0mijf,ϕ*f′=ρϕ*f,Kn′=c1νT0Kn,Kt′=c1νT0Kt,Kθ′=1νc1Kθ,
where ω*=ρc*c12/K*, c12=(λ+2μ)/ρ.
The displacement components u1, u3 and velocity components v1, v3 are related to the potential functions ϕ, ϕf and ψ, ψf in dimensionless form as
(10)(u1,v1)=(∂∂x1(ϕ,ϕf)-∂∂x3(ψ,ψf)),(u3,v3)=(∂∂x3(ϕ,ϕf)+∂∂x1(ψ,ψf)).
Using (10) in (1) and (3)–(5) and with the aid of (8) and (9) (after suppressing the primes), we obtain
(11)∇2ϕ-(1+τ1∂∂t)T-∂2ϕ∂t2=0,∇2ψ-a1∂2ψ∂t2=0,∇2T=(1+τ0∂∂t)∂T∂t+ε1(∂∂t+η0τ0∂2∂t2)∇2ϕ,∇2ϕf-b1Tf-b2ϕ*f-b3∂ϕf∂t=0,∇2ψf+b4Ψ2-b5∂ψf∂t=0,∇2Tf-b9∇2ϕf-b10∂Tf∂t=0,∇2Ψ2-b6∇2ψf-b7Ψ2-b8∂Ψ2∂t=0,b11∇2ϕf-∂∂tϕ*f=0,
where
(12)a1=ρc12μ,ε1=ν2T0K*ω*ρ,b1=bρc12(λf+2μf+Kf)ω*ν,b2=c0c12(λf+2μf+Kf)ω*νT0,b3=ρfc12(λf+2μf+Kf)ω*,b4=Kfμf+Kf,b5=ρfc12(μf+Kf)ω*,b6=Kfc12γfω*2,b7=2b6,b8=Ic12γfω*,b9=bνT0fK1*ρω*,b10=ρfaT0fc12K1*ω*,b11=νT0ρfc12.
4. Boundary Conditions
The boundary conditions at the interface x3=0 are defined as
(13)t33f=Kn(∂u3∂t-v3),t31f=Kt(∂u1∂t-v1),K1*∂Tf∂x3=Kθ(T-Tf),t33=t33f,t31=t31f,m32f=0,K*∂T∂x3=K1*∂Tf∂x3,
where Kn, Kt, and Kθ are the normal force stiffness, transverse force stiffness, and thermal contact conductance coefficients of unit layer thickness having dimensions Nsec/m3, Nsec/m3, and N/msecK.
5. Reflection and Transmission
The longitudinal wave (P-wave) or thermal wave (T-wave) or transverse wave (SV-wave) propagating through the medium M1 which is designated as the region x3>0 is considered. The plane wave is taken to be incident at the plane x3=0 with its direction of propagation with angle θ0 normal to the surface. Each incident wave corresponds to reflected P-wave, T-wave, and SV-wave in medium M1 and transmitted longitudinal wave (L-wave), thermal wave (T-wave), and transverse longitudinal wave coupled with transverse microrotational wave (C-I and C-II waves) in medium M2 as shown in Figure 1.
In order to solve (11), we assume the solutions of the system of the form
(14){ϕ,T,ψ,ϕf,ϕ*f,Tf,ψf,Ψ2}={ϕ¯,T¯,ψ¯,ϕf¯,ϕ*f¯,Tf,¯ψf¯,Ψ2¯}eι{k(x1sinθ-x3cosθ)-ωt},
where k is the wave number, ω is the angular frequency, and ϕ¯, T¯, ψ¯, ϕf¯, ϕ*f¯, Tf¯, ψf¯, Ψ2¯ are arbitrary constants.
Making use of (14) in (11), we obtain
(15)V4+D1V2+E1=0,(16)V4+D2V2+E2=0,V4+D3V2+E3=0,
where
(17)D1=-1τ00[1+(1-ιτ1ω)ε1(ιω+η0τ0)]-1,E1=1τ00,D2=ιωb3(1-ιb2b11ω)+ιωb10(1+ιb1b9ωb3),E2=-ω2b3b10(1-ιb2b11ω),D3=ιωb5+ιω[(1-ιb4b6/ωb5)(b8+(ι/ω)b7)],E3=-ω2(b8+(ι/ω)b7)b5,V2=ω2k2,τ00=(ιω+τ0).
Here V1, V2 are the velocities of P-wave and T-wave in medium M1 and these are roots of (15) and V3=1/a1 is the velocity of SV-wave in medium M1. V1¯, V2¯, V3¯, and V4¯ are the velocities of transmitted longitudinal wave (L-wave), thermal wave (T-wave), and transverse longitudinal wave coupled with microrotational wave (C-I and C-II) in medium M2. These are roots of (16).
In view of (14), the appropriate solutions of (11) for medium M1 and medium M2 are taken as follows.
Medium M1 is as follows:
(18){ϕ,T}=∑i=12{1,fi}[S0ieι{ki(x1sinθ0i-x3cosθ0i)-ωit}+Pi],ψ=S03eι{k3(x1sinθ03-x3cosθ03)-ω3t}+S3eι{k3(x1sinθ03+x3cosθ03)-ω3t},
where
(19)fi=ε1ωi2(ι/ωi+η0τ0)-1/Vi2+(ι/ωi+τ0)-ιε1ωi(ι/ωi+η0τ0)(ι/ωi+τ1),Pi=Sieι{ki(x1sinθi+x3cosθi)-ωit}.
Medium M2 is as follows:
(20){ϕf,Tf,ϕ*f}=∑i=12{1,fi¯,gi¯}Si¯eι{ki¯(x1sinθi¯-x3cosθi¯)-ωi¯t},(21){ψf,Ψ2}=∑j=34{1,fj¯}Sj¯eι{kj¯(x1sinθj¯-x3cosθj¯)-ωj¯t},
where
(22)fi¯=-b3b9b1b9(ι/ωi)+ιωi(1-ιb2b11/ωi)(1/Vi¯2-b10(ι/ωi)),fj¯=(b6b5(ι/ωi))-1/Vj¯2-b7/ωj2+b8(ι/ωi)+b4b6/ωj2,gi¯=-ιb11ωiVi¯2,
and S0i, S03 are the amplitudes of incident (P-wave, T-wave) and SV-wave, respectively. Si and S3 are the amplitudes of reflected (P-wave, T-wave) and SV-wave and Si¯, Sj¯ are the amplitudes of transmitted longitudinal wave, thermal wave, and transverse longitudinal wave coupled with transverse microrotational wave, respectively.
Equation (20) represents the relation between (ϕf and Tf) and (ϕf and ϕf*).
Snell’s law is given by
(23)sinθ0V0=sinθ1V1=sinθ2V2=sinθ3V3=sinθ1¯V1¯=sinθ2¯V2¯=sinθ3¯V3¯=sinθ4¯V4¯,
where
(24)k1V1=k2V2=k3V3=k1¯V1¯=k2¯V2¯=k3¯V3¯=k4¯V4¯=ω,hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhatx3=0.
Making use of (18)–(21) in the boundary conditions (13) and with the help of (2), (7), (9), (10), (23), and (24), we obtain a system of seven nonhomogeneous equations which can be written as
(25)∑j=17aijZj=Yi;(i=1,2,3,4,5,6,7),
where the values of aij are given in the Appendix.
For incident P-wave,
(26)A*=S01,S02=S03=0,Y1=a11,Y2=-a21,Y3=-a31,Y4=-a41,Y5=a51,Y6=0,Y7=a71.
For incident T-wave,
(27)A*=S02,S01=S03=0,Y1=a12,Y2=-a22,Y3=-a32,Y4=-a42,Y5=a52,Y6=0,Y7=a72.
For incident SV-wave,
(28)A*=S03,S01=S02=0,Y1=-a13,Y2=a23,Y3=a33=0,Y4=-a43,Y5=-a53,Y6=0,Y7=a73=0,Z1=S1A*,Z2=S2A*,Z3=S3A*,Z4=S1¯A*,Z5=S2¯A*,Z6=S3¯A*,Z7=S4¯A*,
where Z1, Z2, and Z3 are the amplitude ratios of reflected P-wave, T-wave, and SV-wave in medium M1 and Z4, Z5, Z6, and Z7 are the amplitude ratios of transmitted longitudinal wave (L-wave), thermal wave (T-wave), and transverse longitudinal wave coupled with transverse microrotational wave (C-I and C-II waves) in medium M2.
6. Particular Cases6.1. Case I: Normal Force Stiffness
Kn≠0, Kt→∞, and Kθ→∞ in (25) yield the resulting quantities for normal force stiffness and lead a system of seven nonhomogeneous equations given by (A.1) with the changed values of aij as
(29)a2i=-ω2V0sinθ0,a23=ω2V31-V32V02sin2θ0,a24=ιωV0sinθ0,a25=ιωV0sinθ0,a26=ιωV3¯1-V32¯V02sin2θ0,a27=ιωV4¯1-V42¯V02sin2θ0,a3i=fi,a33=0,a34=-f1¯,a35=-f2¯,a36=a37=0.
6.2. Case II: Transverse Force Stiffness
Kt≠0, Kn→∞, and Kθ→∞ in (25) provide the case of transverse force stiffness giving a system of seven nonhomogeneous equations given by (A.1) with the changed values of aij as
(30)a1i=-ω2Vi1-Vi2V02sin2θ0,a13=-ω2V0sinθ0,a14=-{ιωV1¯(1-V12¯V02sin2θ0)},a15=-{ιωV2¯(1-V22¯V02sin2θ0)},a16=ιωV0sinθ0,a17=ιωV0sinθ0,a3i=fi,a33=0,a34=-f1¯,a35=-f2¯,a36=a37=0.
6.3. Case III: Thermal Contact Conductance
Taking Kθ≠0, Kt→∞, and Kn→∞ in (25) corresponds to the case of thermal contact conductance and yields a system of seven nonhomogeneous equations as given by (A.1) with the changed values of aij as
(31)a1i=-ω2Vi1-Vi2V02sin2θ0,a13=-ω2V0sinθ0,a14=-{ιωV1¯(1-V12¯V02sin2θ0)},a15=-{ιωV2¯(1-V22¯V02sin2θ0)},a16=a17=ιωV0sinθ0,a2i=-ω2V0sinθ0,a23=ω2V31-V32V02sin2θ0,a24=ιωV0sinθ0,a25=ιωV0sinθ0,a26=ιωV3¯1-V32¯V02sin2θ0,a27=ιωV4¯1-V42¯V02sin2θ0.
6.4. Case IV: Perfect Bonding
If we take Kn→∞, Kt→∞, and Kθ→∞ in (25), we obtain the case of perfect bonding and yield a system of seven nonhomogeneous equations as given by (A.1) with the changed values of aij as
(32)a1i=-ω2Vi1-Vi2V02sin2θ0,a13=-ω2V0sinθ0,a14=-{ιωV1¯(1-V12¯V02sin2θ0)},a15=-{ιωV2¯(1-V22¯V02sin2θ0)},a16=a17=ιωV0sinθ0,a2i=-ω2V0sinθ0,a23=ω2V31-V32V02sin2θ0,a24=ιωV0sinθ0,a25=ιωV0sinθ0,a26=ιωV3¯1-V32¯V02sin2θ0,a27=ιωV4¯1-V42¯V02sin2θ0,a3i=fi,a33=0,a34=-f1¯,a35=-f2¯,a36=a37=0.
6.5. Subcases
If micropolar heat conducting fluid medium is absent, we obtain the amplitude ratios at the free surface of thermoelastic solid with one relaxation time and two relaxation times.
These results are similar to those obtained by A. N. Sinha and S. B. Sinha [28] for one relaxation time (L-S theory) and Sinha and Elsibai [29] for two relaxation times (G-L theory).
In the absence of upper medium M2, we obtain the amplitude ratios at the free surface of thermoelastic solid for CT-theory (η0=τ0=τ1=0).
These results are similar to those obtained by Deresiewicz [30] for CT-theory.
7. Special Cases
If η0=1, τ1=0 in (25), (A.1), (29), (30), and (31) then we obtain the corresponding amplitude ratios at an interface of thermoelastic solid with one relaxation time and heat conducting micropolar fluid half-space for normal force stiffness, transverse force stiffness, and thermal contact conductance.
If η0=0, τ1>0 in (25), (A.1), (29), (30), and (31) then we obtain the corresponding amplitude ratios at an interface of thermoelastic solid with two relaxation times and heat conducting micropolar fluid half-space for normal force stiffness, transverse force stiffness, and thermal contact conductance.
8. Numerical Results and Discussion
The following values of relevant parameters for both the half-spaces for numerical computations are taken.
Following Singh and Tomar [19], the values of elastic constants for medium M1 are taken as
(33)λ=0.209730×1010Nm-2,μ=0.91822×109Nm-2,ρ=0.0034×103Kgm-3,
and thermal parameters are taken from Dhaliwal and Singh [31]:
(34)ν=0.268×107Nm-2K-1,c*=1.04×103NmKg-1K-1,K*=1.7×102Nsec-1K-1,T0=0.298K,τ0=0.613×10-12sec,τ1=0.813×10-12sec,ω=1.
Following Singh and Tomar [19], the values of micropolar constants for medium M2 are taken as
(35)λf=1.5×108Nsecm-2,μf=0.03×108Nsecm-2,Kf=0.000223×108Nsecm-2,γf=0.0000222Nsec,ρf=0.8×103Kgm-3,I=0.00400×10-16Nsec2m-2.
Thermal parameters for the medium M2 are taken as of comparable magnitude:
(36)T0f=0.196K,K1*=0.89×102Nsec-1K-1,c0=0.005×1011Nsec2m-6,a=1.5×105m2sec-2K-2,b=1.6×105Nm-2K-1.
The values of amplitude ratios have been computed at different angles of incidence.
In Figures 2–22, for L-S theory, we represent the solid line for stiffness (ST1), small dashes line for normal force stiffness (NS1), medium dashes line for transverse force stiffness (TS1), and dash dot dash line for thermal contact conductance (TCS1). For G-L theory, we represent the dash double dot dash line for stiffness (ST2), solid line with center symbol “plus” for normal force stiffness (NS2), solid line with center symbol “diamond” for transverse force stiffness (TS2), and solid line with center symbol “cross” for thermal contact conductance (TCS2).
Variation of |Z1| with angle of incidence (P-wave).
Variation of |Z2| with angle of incidence (P-wave).
Variation of |Z3| with angle of incidence (P-wave).
Variation of |Z4| with angle of incidence (P-wave).
Variation of |Z5| with angle of incidence (P-wave).
Variation of |Z6| with angle of incidence (P-wave).
Variation of |Z7| with angle of incidence (P-wave).
Variation of |Z1| with angle of incidence (P-wave).
Variation of |Z2| with angle of incidence (P-wave).
Variation of |Z3| with angle of incidence (T-wave).
Variation of |Z4| with angle of incidence (T-wave).
Variation of |Z5| with angle of incidence (T-wave).
Variation of |Z6| with angle of incidence (T-wave).
Variation of |Z7| with angle of incidence for T-wave.
Variation of |Z1| with angle of incidence (SV-wave).
Variation of |Z2| with angle of incidence (SV-wave).
Variation of |Z3| with angle of incidence (SV-wave).
Variation of |Z4| with angle of incidence (SV-wave).
Variation of |Z5| with angle of incidence (SV-wave).
Variation of |Z6| with angle of incidence (SV-wave).
Variation of |Z7| with angle of incidence (SV-wave).
8.1. P-Wave Incident
Variations of amplitude ratios |Zi|, 1≤i≤7, with the angle of incidence θ0, for incident P-wave are shown in Figures 2–8.
Figure 2 shows that the values of |Z1| for NS1 and NS2 increase in the whole range. The values for ST1, ST2, transverse force stiffness, and thermal contact conductance decrease in the whole range, except near the grazing incidence, where the values get increased. The values for ST1 remain more than the values for TS1, TS2, ST2, TCS1, and TCS2 in the whole range.
From Figure 3 it is evident that the values of |Z2| for all the stiffnesses, except ST1 and ST2, decrease in the whole range. The values of |Z2| for NS1 and NS2 are magnified by multiplying by 10.
Figure 4 shows that the values of |Z3| for all the stiffnesses for L-S theory and G-L theory first increase up to intermediate range and then decrease with the increase in θ0 and the values for ST2 are greater than the values for ST1, NS1, NS2, TS1, TS2, TCS1, and TCS2 in the whole range. The values of |Z2| for NS1 and NS2 are magnified by multiplying by 10.
From Figure 5 it is noticed that the values of |Z4| for all the stiffnesses start with maximum value at normal incidence and then decrease to attain minimum value at grazing incidence. The values of amplitude ratios for ST1 and NS1 are more than the values for ST2 and NS2, respectively, in the whole range. There is slight difference in the values of TCS1 and TCS2 in the whole range. The values of |Z4| for ST1, ST2, NS1, and NS2 are magnified by multiplying by 102 and the values for TS1, TS2, TCS1, and TCS2 are magnified by multiplying by 10.
Figure 6 shows that the values of |Z5| for all the boundary stiffnesses decrease with increase in θ and the values of amplitude ratio for TS1 are greater than the values for all other boundary stiffnesses in the whole range that shows the effect of transverse force stiffness. The values of |Z5| for all the stiffnesses are magnified by multiplying by 102.
From Figure 7 it is evident that the amplitude of |Z6| for normal force stiffness increases in the range 0°<θ0<48° and then decreases in the further range. The values for transverse force stiffness increase in the range 0°<θ0<59° and for thermal contact conductance increase in the range 0°<θ0<56° and then decrease. The values of |Z6| for all the stiffnesses, except TCS1 and TCS2, are magnified by multiplying by 103, while the values for TCS1 and TCS2 are magnified by multiplying by 102.
Figure 8 depicts that the values of |Z7| for ST1, ST2, TS1 and TS2 increase in the range 0°<θ0<56° and decrease in the remaining range. The values for ST2 and TS2 are greater than the values for ST1 and TS2 respectively in the whole range. The values for normal force stiffness for G-L theory remain more than the values for L-S theory. The values of |Z7| for ST1, ST2 are magnified by multiplying by 108, the values for TCS1, TCS2 are magnified by multiplying by 106 and the values for NS1, NS2, TS1 and TS2 are magnified by multiplying by 107.
8.2. T-Wave Incident
The values of amplitude ratio |Z1| for transverse force stiffness and thermal contact conductance decrease in the whole range with slight increase in the initial range. The values for ST1 and ST2 increase in the range 0°<θ0<48° and then decrease. These variations have been shown in Figure 9. The values for TS1, TS2, TCS1 and TCS2 are reduced by dividing by 10.
The values of amplitude ratio |Z2| for NS1 and NS2 increase in the whole range, while the values for TS1, TS2, ST1, ST2, TCS1, TCS2 first decrease and then increase to attain maximum value at grazing incidence. These variations are shown in Figure 10.
From Figure 11, it is noticed that the values of |Z3| for normal force stiffness, transverse force stiffness, and thermal contact conductance for G-L theory are greater than the corresponding values for L-S theory. It is seen that the values for ST1 are greater than the values for ST2 in the range 0°<θ0<38° and, in the further range, the values for ST2 are more.
Figure 12 depicts that the values of amplitude ratios |Z4| for transverse force stiffness and thermal contact conductance oscillate up to intermediate range and then decrease in the further range to attain minimum value at grazing incidence. The values for ST1, ST2, NS1, and NS2 decrease in the whole range. The values of |Z4| for ST1 and ST2 are magnified by multiplying by a factor of 102 and NS1, NS2, TS1, TS2, TCS1, and TCS2 are magnified by a factor of 10.
Figure 13 shows that the behavior of variation of |Z5| for all the boundary stiffnesses is similar to that of |Z4|, but the magnitude of variation is different. The values of |Z5| for ST1 and ST2 are magnified by multiplying by a factor of 102 and NS1, NS2, TS1, TS2, TCS1, and TCS2 are magnified by a factor of 10.
From Figure 14 it is evident that the values of |Z6| for all the boundary stiffnesses increase to attain maximum value and then decrease up to grazing incidence. The values for NS1 and TS1 are greater than the values for NS2 and TS2, respectively, that reveals the thermal relaxation time effect. The values of |Z5| for ST1 and ST2 are magnified by multiplying by a factor of 103 and NS1, NS2, TS1, TS2, TCS1, and TCS2 are magnified by a factor of 102.
Figure 15 depicts that the values of |Z7| for ST1 and ST2 increase in the interval 0°<θ0<39° and then decrease in the further range. It is noticed that there is slight difference in the values of amplitude ratio for L-S and G-L theory. The values of |Z7| for ST1 and ST2 are magnified by multiplying by a factor of 108 and NS1, NS2, TS1, TS2, TCS1, and TCS2 are magnified by a factor of 106.
8.3. SV Wave Incident
Variations of amplitude ratios |Zi|, 1≤i≤7, with the angle of incidence θ0, for incident SV-wave are shown in Figures 16–22.
Figure 16 depicts that the values of |Z1| for transverse force stiffness and thermal contact conductance increase to attain peak values in the range 15°<θ0<25° and then decrease in the further range. The values for ST1 and ST2 increase in the range 0°<θ0<35° and 45°<θ0<66°, respectively, and decrease in the remaining range. The values for normal force stiffness increase from normal incidence to attain maximum value in the range 45°<θ0<55° and then decrease up to grazing incidence. The values for TS1, TS2, TCS1, and TCS2 are reduced by dividing by 10.
Figure 17 shows that values of amplitude ratio |Z2| for all the boundary stiffnesses follow oscillatory pattern in the whole range. The maximum value is attained by TCS1 in the range 15°<θ0<25°. It is seen that the values for L-S theory are greater than the values for G-L theory in the whole range.
Figure 18 shows that the values of |Z3| for transverse force stiffness increase from normal incidence to grazing incidence and attain peak value in the interval 15°<θ0<25°. The values for normal stiffness and thermal contact conductance decrease in the intervals 0°<θ0<56° and 0°<θ0<17°, respectively, and then increase in the further range.
Figure 19 shows that the values of |Z4| for all the boundary stiffnesses oscillate in the whole range. The values for ST1 are greater than the values for ST2 in the whole range, except some intermediate range. The values of |Z4| for ST1 and ST2 are magnified by multiplying by a factor of 103 and NS1, NS2, TS1, and TS2 by a factor of 102 and TCS1 and TCS2 are magnified by a factor of 10.
Figure 20 depicts that the behavior of variation of |Z5| for all the boundary stiffnesses is similar to that of |Z4|, but the magnitude of variation is different. The values of |Z5| for ST1, ST2 are magnified by multiplying by a factor of 103 and NS1, NS2, TCS1, TCS2, TS1, and TS2 are magnified by a factor of 102.
From Figure 21, it is noted that the values of |Z6| for ST1 and ST2 decrease in the range 0°<θ0<45°and then increase with increase in angle of incidence. The values of amplitude ratio for transverse force stiffness and thermal contact conductance decrease in the whole range, except the ranges 16°<θ0<24° and 17°<θ0<23°, respectively. The amplitude ratio for normal force stiffness attains maximum value at normal incidence. The values of |Z6| for ST1, ST2 are magnified by multiplying by a factor of 103 and NS1, NS2, TCS1, TCS2, TS1, and TS2 are magnified by a factor of 102.
Figure 22 depicts that the values of |Z7| for all the boundary stiffnesses attain maximum value at the normal incidence and then decrease with oscillation to attain minimum value at the grazing incidence. The values of |Z7| for ST1 and ST2 are magnified by multiplying by a factor of 108 and NS1, NS2, TS1, TS2, TCS1, and TCS2 are magnified by a factor of 106.
9. Conclusion
Reflection and transmission at an interface between heat conducting elastic solid and micropolar fluid media are discussed in the present paper. Effect of normal force stiffness, transverse force stiffness, thermal contact conductance, and thermal relaxation times is observed on the amplitude ratios for incidence of various plane waves (P-wave, T-wave, and SV-wave). When P-wave is incident, it is noticed that the values of amplitude ratio for transverse force stiffness for transmitted T-wave are greater than all the other boundary stiffnesses. When plane wave (SV-wave) is incident, the trend of variation of amplitude ratio for transmitted transverse wave coupled with transverse microrotational wave, that is, C-I and C-II waves, is similar, but magnitude of oscillation is different. The values of amplitude ratio of transmitted LD-wave and T-wave for L-S theory are greater than the value for G-L theory (when T-wave is incident). The model considered is one of the more realistic forms of earth models and it may be of interest for experimental seismologists in exploration of valuable materials such as minerals and crystal metals.
Appendix
Consider the following:
(A.1)a1i=-c1Knω2Vi1-Vi2V02sin2θ0,a13=-c1Knω2V0sinθ0,a14=-{(1-V12¯V02sin2θ0)[d2f+d3f(1-V12¯V02sin2θ0)]ω2V12¯+d1fff1¯+g1¯d1f+ιc1KnωV1¯(1-V12¯V02sin2θ0)},a15=-{(1-V22¯V02sin2θ0)[d2f+d3f(1-V22¯V02sin2θ0)]ω2V22¯+d1fff1¯+g2¯d1f+ιc1KnωV2¯(1-V22¯V02sin2θ0)},a16=d3fω2V3¯V0sinθ01-V32¯V02sin2θ0+ιc1KnωV0sinθ0,a17=d3fω2V4¯V0sinθ01-V42¯V02sin2θ0+ιc1KnωV0sinθ0,a2i=-c2Ktω2V0sinθ0,a23=c2Ktω2V31-V32V02sin2θ0,a24=(2d4f+d5f)ω2V1¯V0sinθ01-V12¯V02sin2θ0+ιc2KtωV0sinθ0,a25=(2d4f+d5f)ω2V2¯V0sinθ01-V22¯V02sin2θ0+ιc2KtωV0sinθ0,a26=d4fω2V32¯(1-2V32¯V02sin2θ0)+d5f[1-V32¯V02sin2θ0-f3¯]+ιc2KtωV3¯1-V32¯V02sin2θ0,a27=d4fω2V42¯(1-2V42¯V02sin2θ0)+d5f[1-V42¯V02sin2θ0-f4¯]+ιc2KtωV4¯1-V42¯V02sin2θ0,a3i=c4Kθfi,a33=0,a34=[ιωV1¯1-V12¯V02sin2θ0-c4Kθ]f1¯,a35=[ιωV2¯1-V22¯V02sin2θ0-c4Kθ]f2¯,a36=a37=0,a4i=(d1+d2(1-Vi2V02sin2θ0))ω2Vi2+(1-τ1ιω)fi,a43=d2ω2V3V0sinθ01-V32V02sin2θ0,a44=-[(d2f+d3f(1-V12¯V02sin2θ0))ω2V12¯+(d1fff1¯+d1fg1¯)],a45=-[(d2f+d3f(1-V22¯V02sin2θ0))ω2V22¯+(d1fff1¯+d1fg1¯)],a46=d3fω2V3¯V0sinθ01-V32¯V02sin2θ0,a47=d3fω2V4¯V0sinθ01-V42¯V02sin2θ0,a51=-(2d4)ω2V1V0sinθ01-V12V02sin2θ0,a52=-(2d4)ω2V2V0sinθ01-V22V02sin2θ0,a53=d4ω2V32(1-2V32V02sin2θ0),a54=-(2d4f+d5f)ω2V1¯V0sinθ01-V12¯V02sin2θ0,a55=-(2d4f+d5f)ω2V2¯V0sinθ01-V22¯V02sin2θ0,a56=-[d5f(ω2V32¯(1-V32¯V02sin2θ0)-f3¯)d4fω2V32¯(1-2V32¯V02sin2θ0)+d5f(ω2V32¯(1-V32¯V02sin2θ0)-f3¯)],a57=-[(ω2V42¯(1-V42¯V02sin2θ0)-f4¯)d4fω2V42¯(1-2V42¯V02sin2θ0)+d5f(ω2V42¯(1-V42¯V02sin2θ0)-f4¯)],a6i=a63=a64=a65=0,a66=ιωV3¯f3¯1-V32¯V02sin2θ0,a67=ιωV4¯f4¯1-V42¯V02sin2θ0,a7i=ιωVi1-Vi2V02sin2θ0fi,a73=0,a74=ιp2ωV4¯1-V12¯V02sin2θ0f1¯,a75=ιp2ωV5¯1-V22¯V02sin2θ0f2¯,a76=a77=0,(i=1,2),d1=λρc12,d2=2μρc12,d3=d22,d1ff=bν,d1f=c0ρνT0,d2f=λfω*ρc12,d3f=(2μf+Kf)ω*ρc12,d4f=μfω*ρc12,d5f=Kfω*ρc12,p2=K1*K*,c1=c2=νT0ρc12,c4=νc12ω*K1*.
Symbolsλ, μ:
Lame’s constants
tij:
Components of the stress tensor
u→:
Displacement vector
ρ:
Density
K*:
Thermal conductivity
c*:
Specific heat at constant strain
T0:
Uniform temperature
T:
Temperature change
αT:
Coefficient of linear thermal expansion
δij:
Kronecker delta
εijr:
Alternating symbol
λf, μf, Kf,αf, βf, γf, c0:
Material constants of the fluid
σijf:
Components of stress tensor in the fluid
mijf:
Components of couple stress tensor in the fluid
v→:
Velocity vector
Ψ→:
Microrotation velocity vector
ρf:
Density
I:
Scalar constant with the dimension of moment of inertia of unit mass
p:
Pressure
K1*:
Thermal conductivity
aT0f:
Specific heat at constant strain
T0f:
Absolute temperature
Tf:
Temperature change
ϕ*f:
Variation in specific volume
αTf:
Coefficient of linear thermal expansion.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
Thanks are due to Punjab Technical University, Kapurthala, for providing research facilities and enrolling one of the authors (Mandeep Kaur) as a research scholar to Ph.D. Programme.
EringenA. C.Simple microfluids1964220521710.1016/0020-7225(64)90005-9MR0169468ZBL0136.450032-s2.0-49749223867EringenA. C.Theory of micropolar fluids196616118MR0204005ArimanT.TurkM. A.SylvesterN. D.Microcontinuum fluid mechanics-A review197311890593010.1016/0020-7225(73)90038-4ZBL0259.760012-s2.0-0015656660ArimanT.SylvesterN. D.TurkM. A.Applications of microcontinuum fluid mechanics197412427329310.1016/0020-7225(74)90059-72-s2.0-0016046496ŘíhaP.On the theory of heat-conducting micropolar fluids with microtemperatures1975231-21810.1007/BF01177664ZBL0326.760062-s2.0-0016587526EringenA. C.KafadarC. B.EringenA. C.Polar field theories19764New York, NY, USAAcademic PressBrulinO.1982edited by O. Brulin and R. K. T. HsiehSingaporeWorld ScientificGorlaR. S. R.Combined forced and free convection in the boundary layer flow of a micropolar fluid on a continuous moving vertical cylinder1989271778610.1016/0020-7225(89)90169-92-s2.0-0024862390EringenA. C.Theory of thermo-microstretch fluids and bubbly liquids199028213314310.1016/0020-7225(90)90063-OMR1040726ZBL0709.760152-s2.0-0025555720AydemirN. U.VenartJ. E. S.Flow of a thermomicropolar fluid with stretch199028121211122210.1016/0020-7225(90)90069-U2-s2.0-0025558675CiarlettaM.Spatial decay estimates for heat-conducting micropolar fluids200139665566810.1016/S0020-7225(00)00056-2MR1815358ZBL1210.760142-s2.0-0035314294HsiaS.ChengJ.Longitudinal plane wave propagation in elastic-micropolar porous media20064531743174810.1143/JJAP.45.17432-s2.0-33644916155HsiaS. Y.ChiuS. M.SuC. C.ChenT. H.Propagation of transverse waves in elastic-micropolar porous semispaces200746117399740510.1143/JJAP.46.73992-s2.0-35948950277TomarS. K.GognaM. L.Reflection and refraction of a longitudinal microrotational wave at an interface between two micropolar elastic solids in welded contact199230111637164610.1016/0020-7225(92)90132-ZZBL0825.731542-s2.0-0026954139TomarS. K.GognaM. L.Reflection and refraction of longitudinal wave at an interface between two micropolar elastic solids in welded contact199597282283010.1121/1.4130862-s2.0-0028795350TomarS. K.GognaM. L.Reflection and refraction of coupled transverse and micro-rotational waves at an interface between two different micropolar elastic media in welded contact199533448549610.1016/0020-7225(94)00077-8ZBL0899.730922-s2.0-0029267738KumarR.SharmaN.RamP.Reflection and transmission of micropolar elastic waves at an imperfect boundary200841153610.1163/1573611087834703882-s2.0-38649140669KumarR.SharmaN.RamP.Interfacial imperfection on reflection and transmission of plane waves in anisotropic micropolar media200849330531210.1016/j.tafmec.2008.02.0072-s2.0-42649105792SinghD.TomarS. K.Longitudinal waves at a micropolar fluid/solid interface200845122524410.1016/j.ijsolstr.2007.07.0152-s2.0-43049112499SunQ.XuH.LiangB.Propagation characteristics of longitudinal displacement wave in micropolar fluid with micropolar elastic plate201169492392710.4028/www.scientific.net/MSF.694.9232-s2.0-80052056374XuH. Y.DangS. Y.SunQ. Y.LiangB.Reflection and transmission characteristics of coupled wave through micropolar elastic solid interlayer in micropolar fluid201326347354DangS. Y.XuH. Y.SunQ. Y.LiangB.Propagation characteristics of coupled wave through micropolar fluid interlayer in micropolar elastic solid201372172973210.4028/www.scientific.net/AMR.721.7292-s2.0-84882966388FuC.WeiP. J.The wave propagation through the imperfect interface between two micropolar solids201380341942210.4028/www.scientific.net/AMR.803.419SharmaK.MarinM.Reflection and transmission of waves from imperfect boundary between two heat conducting micropolar thermoelastic solids201422151175VishwakarmaS. K.GuptaS.VermaA. K.Torsional wave propagation in Earth's crustal layer under the influence of imperfect interface201420335536910.1177/1077546312461029MR3179357LordH. W.ShulmanY.A generalized dynamical theory of thermoelasticity196715529930910.1016/0022-5096(67)90024-5ZBL0156.227022-s2.0-49949149766GreenA. E.LindsayK. A.Thermoelasticity1972211710.1007/BF000456892-s2.0-0010974262SinhaA. N.SinhaS. B.Reflection of thermoelastic waves at a solid half-space with thermal relaxation19742222372442-s2.0-0016314215SinhaS. B.ElsibaiK. A.Reflection of thermoelastic waves at a solid half-space with two relaxation times199619874976210.1080/014957396089462052-s2.0-0030286755DeresiewiczH.Effect of boundaries on waves in a thermoelastic solid: reflexion of plane waves from a plane boundary1960816417210.1016/0022-5096(60)90035-1MR0118105ZBL0122.194062-s2.0-0001472688DhaliwalR. S.SinghA.1980New Delhi, IndiaHindustan Publication Corporation