Abstract
We study the reflection and transmission coefficients of plane waves incident at an interface between two isotropic thermoelastic half spaces and compare them with those of the elastic case. The models include the classical-Biot (B) and extended Lord-Shulman (LS) theories, and predict reflected and transmitted fast-compressional (P), thermal (T) and shear (S) waves. The coefficients are formulated in terms of incidence and inhomogeneity angles, medium properties and potential functions. We consider different incident wave types and inhomogeneity angles to analyze the magnitude, phase and energy ratio of the plane waves, and perform a comparison with the isothermal (elastic) theory. The thermoelastic and elastic models predict different energy partitions between the P and S modes, satisfying the conservation of energy. The LS model exhibits higher T-wave thermal attenuation with increasing inhomogeneity angle at high frequencies, accordingly predicting more interference energy. The angle affects the energy partitions, particularly at the critical angle and near grazing incidence for an incident P wave, which satisfies the conservation of energy. Beyond the critical angle, the energy flux perpendicular to the interface of the isothermal model vanishes, while it is significant in the thermoelastic case. The T-wave magnitudes increase when the thermal conductivity (relaxation time) increases.
1 Introduction
The thermoelasticity theory couples the fields of elastic deformation and temperature and has been largely studied in several fields during the past decades, such as mechanics of materials, ultrasonics and to a much lesser extent in exploration geophysics (Zener, 1938; Savage, 1966; Armstrong, 1984; Veres et al., 2013; Wang and Li, 2013; Carcione et al., 2019). It is also of interest in geothermal applications (Buijze et al., 2019) and generally in the analysis of deep reservoirs, where temperature effects are important, mainly in relation to seismic-reflection technology, whose physics is based on the reflection and transmission (R/T) of waves.
Biot (1956) introduced the classical thermoelasticity theory, hereafter B theory, based on the parabolic Fourier heat conduction law, where the P or T waves propagate with unrealistic infinite phase speeds at high frequencies (Deresiewicz, 1957; Rudgers, 1990, Figures 1–3). Lord and Shulman (1967) introduced a relaxation term into the thermoelasticity equations to obtain finite wave velocities (hereafter LS theory), which led to the Maxwell-Cattaneo-Vernotte hyperbolic heat equation (Maxwell, 1867; Vernotte, 1948; Cattaneo, 1958). The involved relaxation time represents the time lag from imposing the temperature disturbance to the establishment of the steady-state regime. The LS theory predicts two compressional waves (P and T) and a shear wave having similar characteristics to the waves of poroelasticity (Biot, 1956), where the existence of the T wave has been verified experimentally (Ackerman et al., 1966; Jackson et al., 1970; McNelly et al., 1970) and numerically (Carcione et al., 2019; Wang Z.-W. et al., 2020). Green and Lindsay (1972) developed a similar thermoelasticity theory by introducing additional relaxation times. Generalized approaches, based on fractional derivatives, have been developed by Kumar et al. (2013) and Hobiny and Abbas (2020).
Research on reflection and transmission phenomena in elastic, anelastic and poroelastic media (e.g., Pilant, 1979) involves numerous approaches (Rokhlin et al., 1986; Denneman et al., 2002; Carcione, 2014; Wang E. et al., 2020). Here, we consider the theory of thermoelasticity, which is more general and provides realistic results. The reflection of thermoelastic waves at the free surface of an elastic half-space, based on the generalized Green-Lindsay theory, has been studied by Sinha and Elsibai (1996) and a similar problem has been attacked by Sharma et al. (2003) and Zenkour et al. (2013). Kumar and Sarthi (2006) solved the R/T problem, but ignoring energy dissipation, and Singh and Chakraborty (2013) assumed an initial stress at a solid-liquid interface, while Sharma (2018) considered a poroelastic/elastic interface. More recently, Sarkar and Mondal (2020) considered a stress-free and thermally insulated surface on the basis of the modified Green-Lindsay theory, but their formulation neglects the presence of inhomogeneous plane waves, violating Snell’s law, which is inappropriate but appears in many papers (e.g., Sinha and Elsibai, 1997; Wei et al., 2016; Sarkar et al., 2020). Wang et al. (2021) studied the scattering coefficients at a free surface in the framework of the thermo-poroelasticity theory.
First, we consider the B and LS theories, based on the Helmholtz potential function decomposition law, and obtain the respective inhomogeneous plane-wave solutions and dispersion relations. Then, using the boundary conditions and Snell’s law, we obtain the R/T coefficients for incident P and S waves at an interface between two thermoelastic media and compare them with the elastic case to illustrate the difference between the two theories and the influence of the inhomogeneity angle. Moreover, we verify the conservation of energy and discuss the variations of the coefficients as a function of the thermal conductivity and relaxation time.
2 Thermoelasticity
2.1 Governing Equations
The equations of thermoelasticity describe the relation between the stress-deformation and the temperature fields. We consider the generalized (LS) theory proposed by Lord and Shulman (1967). Let us define by ui, i = x, y, z the displacement components and by T the increment of temperature field above the reference absolute temperature T0 for the state of zero stress and strain. In a linear isotropic medium, combining the stress-strain and strain-displacement relations with the momentum conservation equation (Carcione et al., 2019), we obtain the displacement equation of motion and the law of heat conduction:where ρ is the mass density, is the thermal conductivity, c is the specific heat of the unit volume in the absence of deformation, τ is the relaxation time, q is the external heat source, a dot above a variable denotes time differentiation and the Einstein implicit summation is assumed. The thermal modulus is , where λ and μ are the Lamé constants and is the linear thermal expansion coefficient. More details are given in Appendix A. In the B theory, τ = 0, and the heat equation is parabolic (diffusion-like), but it is hyperbolic (wave-like) in the LS theory.
2.2 Plane-Wave Solution
Let us consider that the displacement vector u can be described by a Helmholtz decomposition of the two potential functions:Substituting Eq. 2 into Eq. 1, we obtainwhere ∇2 is the Laplacian operator and the plane-wave versions of the potential and temperature arewhere Aϕ and AT are the amplitudes, ω is the angular frequency, t is the time variable, k is the complex wavenumber vector for the compressional waves (Chadwick, 1960; Jiao et al., 2019), x is the spatial vector and i2 = −1, wherewhere Re (⋅) and Im (⋅) denote real and imaginary parts, respectively, γ is the inhomogeneity angle, κ and α are the magnitudes of the real wavenumber and attenuation vectors, respectively, and the directions of these vectors arewhere the hat denotes a unit vector, θ is the angle between the real wavenumber vector and a line perpendicular to the interface. Thus, we can setWhen γ is zero, the wave is homogeneouswhere vc is the complex velocity.
Substituting the plane waves Eq. 7 into Eqs. 31 and 32, using (5) we obtainThe condition det(H) = 0 yields the dispersion equationwhereEq. 10 has the solutionsThere are two P-wave wavenumbers, a fast P-wave (minus sign) and a T-wave (plus sign). We have (Carcione, 2014, Eq. 3.34)For an inhomogeneous wave, the phase velocity and attenuation factor are (Carcione, 2014; Carcione et al., 2019)For a homogeneous wave with γ = 0°, these quantities, expressed in terms of the real and imaginary parts of the complex velocity in Eq. 8, areand the attenuation coefficient is(Deresiewicz, 1957).
Similarly, considering an S plane wavewhere Aψ is the amplitude and replacing Eq. 17 into (3)3, we haveso that Im (kS) = 0, and the corresponding wavenumber and phase velocity arerespectively. The S wave is not affected by the thermal effects in (homogeneous) thermoelastic media.
3 Reflection and Transmission Coefficients
We consider the 2D case in the (x, z) plane, an interface defined by z = 0, the incidence media I (z > 0) and transmission media II (z < 0). The properties of medium II are represented with a prime superscript. The incident wave generates six waves, namely, three reflected and three transmitted, illustrated in Figure 1 where the subscripts 0, 1, 2 and 3 represent the incident, P, T and S waves, respectively.
FIGURE 1
3.1 Potential Functions
Consider an incident P wave, where the superscripts i, r and t represent incident, reflected and transmitted, respectively. Then, the potential functions arewherewhere Am and are wave amplitudes, pm and are horizontal wavenumbers, qm and are vertical wavenumbers, and the subscripts 0, m = 1, m = 2 and m = 3 correspond to the incident, P, T and S waves, respectively. We consider Snell’s law which establishes the continuity of the horizontal wavenumbers (e.g., Carcione, 2014; Wang E. et al., 2020)wherewhere θ0 is the incidence angle and κ0 and α0 are given by Eq. 13 with k = k0 and γ = γ0 (see Figure 1). Generalized Snell’s law (Borcherdt, 2009, Eq. (5.2.20)) isBecause the phase velocities of the incident P (or T) and reflected P (or T) waves are the same, namely θ0 = θ1 (or θ0 = θ2), we get γ0 = γ1 (or γ0 = γ2) from Eqs. 12, 13. From Eqs 18 and Sharma (2018), we obtainThe reflection angle θ2 (or θ1) and transmission angles and satisfyandwhere q0 is the vertical wavenumber of the incident wave, given bywhere DR and DI denote the real and imaginary parts of the complex quantity D, pv denotes the principal value and the calculations of the vertical wavenumber are similar for the reflected (qm) and transmitted waves. The reflected P (or T) wave is homogeneous if and only if the incident P (or T) wave is homogeneous. If θc is a critical angle for the transmitted P wave, namely , the corresponding inhomogeneity angle iswhere k0 and are the wavenumbers of the incident and transmitted P-waves. If the incident P wave is homogeneous (γ0 = 0°) and not normally incident, using Eqs. 26, 29 the waves are homogeneous if and only if (Borcherdt, 2009)where the subscripts 1 and 2 correspond to the transmitted P and T waves, respectively. Then, substituting Eq. 20 into (3)1, we obtainwhere
On the other hand, for the incident S wave,and, similarly,the corresponding expressions are given in Eqs. 21, 32. From Eq. 25, the incident S wave is inhomogeneous, and so are the reflected and transmitted S waves. The properties of the reflected and transmitted longitudinal waves generated by an incident S wave are analogous to those of the incident P wave in Eqs. 26, 30.
3.2 Zoeppritz Equations
At z = 0, the boundary conditions to be satisfied are the continuity of temperature, heat flux, and normal and tangential displacements and stresses (Ignaczak and Ostoja-Starzewski, 2010), i.e.,
Substituting Eqs. 2, 20 into the boundary conditions (35), we get the Knott equations (Knott, 1899) for the incident P wavewhere
Similarly, we obtain the equations for the incident S wave, where A0 is its amplitude in Eq. 36, and
Using the relations between displacement amplitudes and their ratios, we obtainwhere k0, km and (m = 1, 2, 3) represent the complex wavenumbers of incident, reflected and transmitted waves respectively. and define the reflection and transmission amplitudes, while ϑm and are the relative phase angles.
4 Energy-Flow Balance
Let us consider the balance of energy flux between the incident wave and the reflected and transmitted waves at a surface element of unit area where the energy flux is the scalar product of the surface traction and particle velocity. The time-averaged energy flux is (Carcione, 2014, Eq. 3.106)where denotes a temporal average over a period, the star is the complex conjugate, and each component is the sum of the components of the respective waves. The energy flux perpendicular to the interface is continuous across the interface, because of the continuity of stresses and displacements. Following Sharma (2018) and Wang E. et al. (2020), we consider the energy partition in the z-direction (perpendicular to the interface plane)In the following expressions, we have omitted the subscripts z for simplicity. Denoting by the energy flux of the incidence medium, we obtainwherewhere the diagonal element a = b = 0 of this matrix corresponds to the energy flux of the incident wave, and a = b = 1, 2, 3 to the reflected P, T, and S waves, whereas the off-diagonal elements are the interference fluxes between the incident and reflected waves.
In the transmission medium, the energy flux can be express aswhereThe diagonal elements are the energy fluxes of the transmitted P, T, and S waves. By scaling the energy flux with that of the incident wave , the relative energy ratios areAccording to Eqs. 46 and 43 can be rewritten asEq. 45 is
The energy-balance equation in the z-direction at the interface is (Carcione, 2014, Eq. 6.116; Wang E. et al., 2020)where the sum of the vertical energy ratios result from the interaction between the incident wave and the three reflected waves, as well as interactions among the three reflected waves, isSimilarly, the corresponding vertical interference energy ratio for the transmitted waves isLet us denote . The energy-flow balance verifies the reflection and transmission coefficients.
5 Examples
The effects of thermoelasticity and the inhomogeneity angle are illustrated by considering the magnitudes, phase angles and energy of the reflection and transmission coefficients. We assume the following reference properties of the incidence and transmission media:Incidence medium: ρ = 2054 kg/m3, c = 960 kg/(m s2 K), = 10.5 m kg/(s3 K), = 0.33 × 10–5 K−1, T0 = 300 K, v0 = 2256.5 m/s, vS = 1302 m/s.Transmission medium: ρ′ = 2600 kg/m3, c′ = 960 kg/(m s2 K), = 10.5 m kg/(s3 K), = 0.33 × 10–5 K−1, = 300 K, = 3636.9 m/s, = 2100 m/s.
v0 and vS are the elastic longitudinal and transverse wave velocities, consistent with Eqs 114 and 19. The relaxation time is(Rudgers, 1990), equal to 2.15 and 0.83 ns for the incidence and transmission media, respectively. The corresponding velocities and densities used here are those of Pilant (1979), whose scattering coefficients in the isothermal case are represented in his Figures 12-3 and 12-4, while the thermoelasticity properties are taken from Carcione et al. (2019).
Figures 2, 3 show the phase velocities and attenuation coefficients of the P and T waves in the incidence and transmission media, as a function of frequency, respectively. The solid, dashed and dotted solid lines correspond to inhomogeneity angles of 0°, 40° and 80°, and the superscripts LS and B represent the Lord-Shulman and classical-Biot theories, respectively. We can see that the low-frequency behavior is almost the same, the T wave is strongly dissipative at low frequencies (Figures 2D, 3D), and that the P wave has a relaxation peak (Figures 2B, 3B) caused by the thermal effects at a relaxation frequency of approximatelyThe LS model predicts a finite high-frequency limit velocity and a lower velocity dispersion of the T wave, which is wave-like at high frequencies. Both the P and T waves have abnormal velocity behaviors at and beyond the inflexion point in the B case (velocity decreases) and T waves have infinite speed at high frequencies (Figures 2C, 3C). The thermal relaxation hardly affects the wave propagation at relatively low frequencies. It can be seen that inflexion points in the P-wave velocity occur at a frequency of approximately 100 MHz.
FIGURE 2
FIGURE 3
5.1 P-Wave Incidence
Figures 4, 5 show the absolute values of the reflection and transmission coefficients and the corresponding phase angles as a function of the incidence angle in the case of an incident P wave with a frequency of 100 MHz. The elastic case is also indicated. If the two media are elastic, there is a critical angle at approximately 40°. If the incident wave travels along the z-direction (vertical), an incident P wave will generate reflected and transmitted longitudinal waves without conversion to the shear modes. The incident P wave is homogeneous when γ0 = 0°, and the behavior of the curves is close to that of the elastic case. We can see a pseudo critical angle at around 33° in the LS curves. Note that there is a converted S wave when the propagation angle is 0° and γ0 is non-zero as shown in Figures 4B,D, the effects on the magnitudes are relatively strong at the critical angle, and the phases of the reflected P wave are reversed at large incidence angles.
FIGURE 4
FIGURE 5
The energy ratios are illustrated in Figure 6. Here, we plot the square root of the relative energy ratios to show the reflection and transmission, which is closer to the seismic response. The elastic phase and square root of the energy ratios are presented in Pilant (1979, Figure 12-3) and our results agree. The inhomogeneity angle affects the magnitudes of the reflected and transmitted waves (Figure 4) and consequently the energy partitions. In the elastic case, the energy flux perpendicular to the interface of the transmitted P wave (Figure 6C) vanish at the critical angle, because the wave propagates along the interface, carrying no energy flux in the vertical direction. The influence of γ0 on the reflected P wave (Figure 6A) is evident at the critical angle, while the effect on the other waves is maximum near grazing incidence. The sum of all the energy ratios is −1, which implies energy conservation at the interface as shown in Figure 7.
FIGURE 6
FIGURE 7
5.2 S-Wave Incidence
Next, we consider an incident S wave. Figures 8, 9 show the absolute values of the reflection and transmission coefficients and phase angles as a function of the incidence angle, and Figure 10 the corresponding energy ratios. Because the S-wave wavenumber kS is real (see Eq. 19), the inhomogeneity angle has no effect (see Eq. 25) as can be seen in Figures 8–11 (the three curves overlap). The energy of the reflected and transmitted P waves, at approximately 35° and 20° in Figures 10A,C, vanishes in the elastic (lossless) case, where the wave propagates along the interface and, as for the P wave, the plane wave is evanescent. The curves in Figure 10 are identical to those of Pilant (1979, Figure 12-4). At the S-wave critical angle (38°), the energy ratio of the reflected S wave is 1 (Figure 10B), whereas the energy of other waves is zero, since the wave propagating along the interface carries no energy vertically. In the LS case, the behavior is analogous to the elastic case. However, there is transmission for all incidence angles and non-zero interference energy. Moreover, significant energy conversion occurs, as shown in Figure 11A, i.e., the LS model predicts more interference energy than the elastic one. The energy conservation at the interface is satisfied, as can be seen in Figure 11B. The same results have been illustrated at the free surface of double-porosity and two double-porosity media (Sharma, 2013; Wang E. et al., 2020).
FIGURE 8
FIGURE 9
FIGURE 10
FIGURE 11
5.3 Effect of the Thermal Conductivity (Relaxation Time)
Next, we study the effect of the relaxation time (see Eqs. 52, 53. Increasing , we obtain a higher τ and a lower fr (Carcione et al., 2019; Wang Z.-W. et al., 2020). The thermal conductivities of the incidence medium are assumed to be 0.5 m kg/(s3 K), 10.5 m kg/(s3 K) and 20.5 m kg/(s3 K), respectively, and that of the transmission medium is 10.5 m kg/(s3 K). Figures 12–15 show the effects of on the amplitudes and energy ratios. We observe that it affects more the T wave conversion in the case of an incident P wave (see Figures 12B,C). Higher enhances the T waves at the expense of the S waves for an incident S wave as shown in Figures 14, 15, where the T wave is weaker attenuated, according to Eqs. 52, 53 and curves of Figures 2D, 3D.
FIGURE 12
FIGURE 13
FIGURE 14
FIGURE 15
6 Conclusion
We obtain the reflection and transmission coefficients at an interface separating two thermoelastic media, whose properties are based on the Lord-Shulman theory. Comparing these coefficients with those of the elastic (lossless) case, shows that critical angles, amplitudes and energy ratios, including interference fields, are affected by the presence of the thermal wave. Since the presence of this wave makes the media anelastic, the propagation and attenuation vectors of the incident wave do not necessarily coincide (inhomogeneous plane wave). The effect of the inhomogeneity angle is more noticeable at the critical angle and near grazing incidence angle for an incident P wave, and has no effect for an incident S wave, since the wavenumber of this wave is real. Then, we analyze the influence of the thermal conductivity (relaxation time) of the incidence media, which shows that this property affects the converted thermal wave more than the others, as expected. The calculations have been verified with the conservation of energy. This study is aimed to further improve the understanding of the behavior of the internal structure of the Earth, where the temperature effects are important.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
WH, L-YF, and JC contributed to conception and design of the study. WH wrote the first draft of the manuscript. WH, L-YF, JC, and TH wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
The research is supported by the National Natural Science Foundation of China (Grant No. 41821002), 111 project “Deep-Superdeep Oil & Gas Geophysical Exploration” (B18055) and Innovation fund project for graduate students of China University of Petroleum (East China) (No. CXJJ-2022-17).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2022.850331/full#supplementary-material
References
1
AckermanC. C.BertmanB.FairbankH. A.GuyerR. A. (1966). Second Sound in Solid Helium. Phys. Rev. Lett.16 (18), 789–791. 10.1103/physrevlett.16.789
2
ArmstrongB. H. (1984). Models for Thermoelastic Attenuation of Waves in Heterogeneous Solids. Geophysics49 (7), 1032–1040. 10.1190/1.1441718
3
BiotM. A. (1956). Thermoelasticity and Irreversible Thermodynamics. J. Appl. Phys.27 (3), 240–253. 10.1063/1.1722351
4
BorcherdtR. D. (2009). Viscoelastic Waves in Layered Media. Cambridge University Press.
5
BuijzeL.van BijsterveldtL.CremerH.PaapB.VeldkampH.WassingB. B. T.et al (2019). Review of Induced Seismicity in Geothermal Systems Worldwide and Implications for Geothermal Systems in the Netherlands. Neth. J. Geosciences98, 1–27. 10.1017/njg.2019.6
6
CarcioneJ. M. (2014). Wave fields in Real media. Theory and Numerical Simulation of Wave Propagation in Anisotropic, Anelastic and Porous media. 3rd ed. Elsevier.
7
CarcioneJ. M.WangZ.-W.LingW.SalustiE.BaJ.FuL.-Y. (2019). Simulation of Wave Propagation in Linear Thermoelastic media. Geophysics84 (1), T1–T11. 10.1190/geo2018-0448.1
8
CattaneoC. (1958). Sur une forme de l’équation de la chaleur éliminant paradoxe d’une propagation instantantée. Comptes rendus de l’Académie des Sci.247, 431–433.
9
ChadwickP. (1960). Thermoelasticity. The Dynamical Theory. Prog. Solid Mech.1, 263–328. 10.1136/bmj.1.5165.56
10
DennemanA. I. M.DrijkoningenG. G.SmeuldersD. M. J.WapenaarK. (2002). Reflection and Transmission of Waves at a Fluid/porous‐medium Interface. Geophysics67 (1), 282–291. 10.1190/1.1451800
11
DeresiewiczH. (1957). Plane Waves in a Thermoelastic Solid. The J. Acoust. Soc. America29 (2), 204–209. 10.1121/1.1908832
12
GreenA. E.LindsayK. A. (1972). Thermoelasticity. J. Elasticity2 (1), 1–7. 10.1007/bf00045689
13
HobinyA. D.AbbasI. A. (2020). Fractional Order Thermoelastic Wave Assessment in a Two-Dimension Medium with Voids. Geomechanics Eng.21 (1), 85–93.
14
IgnaczakJ.Ostoja-StarzewskiM. (2010). Thermoelasticity with Finite Wave Speeds. Oxford, United Kingdom: Oxford University Press.
15
JacksonH. E.WalkerC. T.McNellyT. F. (1970). Second Sound in NaF. Phys. Rev. Lett.25 (1), 26–28. 10.1103/physrevlett.25.26
16
JiaoF.WeiP.ZhouX.ZhouY. (2019). The Dispersion and Attenuation of the Multi-Physical fields Coupled Waves in a Piezoelectric Semiconductor. Ultrasonics92, 68–78. 10.1016/j.ultras.2018.09.009
17
KnottC. G. (1899). III. Reflexion and Refraction of Elastic Waves, with Seismological Applications. Lond. Edinb. Dublin Philosophical Mag. J. Sci.48, 64–97. 10.1080/14786449908621305
18
KumarR.SarthiP. (2006). Reflection and Refraction of Thermoelastic Plane Waves at an Interface between Two Thermoelastic media without Energy Dissipation. Arch. Mech.58 (2), 155–185.
19
KumarR.GuptaV.AbbasI. A. (2013). Plane Deformation Due to thermal Source in Fractional Order Thermoelastic media. Jnl Comp. Theo Nano10 (10), 2520–2525. 10.1166/jctn.2013.3241
20
LordH. W.ShulmanY. (1967). A Generalized Dynamical Theory of Thermoelasticity. J. Mech. Phys. Sol.15 (5), 299–309. 10.1016/0022-5096(67)90024-5
21
MaxwellJ. C. (1867). On the Dynamical Theory of Gases. Philosophical Trans. R. Soc. Lond.157, 49–88.
22
McNellyT. F.RogersS. J.ChanninD. J.RollefsonR. J.GoubauW. M.SchmidtG. E.et al (1970). Heat Pulses in NaF: Onset of Second Sound. Phys. Rev. Lett.24 (3), 100–102. 10.1103/physrevlett.24.100
23
PilantW. L. (1979). Elastic Waves in the Earth. New York: Elsevier Science Publishing Company.
24
RokhlinS. I.BollandT. K.AdlerL. (1986). Reflection and Refraction of Elastic Waves on a Plane Interface between Two Generally Anisotropic media. J. Acoust. Soc. America79 (4), 906–918. 10.1121/1.393764
25
RudgersA. J. (1990). Analysis of Thermoacoustic Wave Propagation in Elastic media. J. Acoust. Soc. America88 (2), 1078–1094. 10.1121/1.399856
26
SarkarN.MondalS. (2020). Thermoelastic Plane Waves under the Modified green–lindsay Model with Two-Temperature Formulation. J. Appl. Mathematics Mech./Z. für Angew. Mathematik Mechanik100 (11), 1–17. 10.1002/zamm.201900267
27
SarkarN.Abo-DahabS. M.MondalS. (2020). Reflection of Magneto-Thermoelastic Waves at a Solid Half-Space under Modified Green-Lindsay Model with Two Temperatures. J. Therm. Stresses43 (9), 1083–1099. 10.1080/01495739.2020.1768991
28
SavageJ. C. (1966). Thermoelastic Attenuation of Elastic Waves by Cracks. J. Geophys. Res.71 (16), 3929–3938. 10.1029/jz071i016p03929
29
SharmaJ. N.KumarV.ChandD. (2003). Reflection of Generalized Thermoelastic Waves from the Boundary of a Half-Space. J. Therm. Stresses26 (10), 925–942. 10.1080/01495730306342
30
SharmaM. D. (2013). Effect of Local Fluid Flow on Reflection of Plane Elastic Waves at the Boundary of a Double-Porosity Medium. Adv. Water Resour.61, 62–73. 10.1016/j.advwatres.2013.09.001
31
SharmaM. D. (2018). Reflection-refraction of Attenuated Waves at the Interface between a Thermo-Poroelastic Medium and a Thermoelastic Medium. Waves in Random and Complex Media28 (3), 570–587. 10.1080/17455030.2017.1370154
32
SinghM. C.ChakrabortyN. (2013). Reflection and Refraction of P-, SV- and thermal Wave, at an Initially Stressed Solid–Liquid Interface in Generalized Thermoelasticity. Appl. Math. Model.37 (1-2), 463–475. 10.1016/j.apm.2012.03.008
33
SinhaS. B.ElsibaiK. A. (1997). Reflection and Refraction of Thermoelastic Waves at an Interface of Two Semi-infinite media with Two Relaxation Times. J. Therm. Stresses20 (2), 129–145. 10.1080/01495739708956095
34
SinhaS. B.ElsibaiK. A. (1996). Reflection of Thermoelastic Waves at a Solid Half-Space with Two Relaxation Times. J. Therm. Stresses19 (8), 749–762. 10.1080/01495739608946205
35
VeresI. A.BererT.BurgholzerP. (2013). Numerical Modeling of Thermoelastic Generation of Ultrasound by Laser Irradiation in the Coupled Thermoelasticity. Ultrasonics53 (1), 141–149. 10.1016/j.ultras.2012.05.001
36
VernotteP. (1948). Théorie continue et théorie moléculaire des phénomènes thermocinétiques. Comptes rendus de l’Académie des Sci.227, 43–44.
37
WangB. L.LiJ. E. (2013). Hyperbolic Heat Conduction and Associated Transient thermal Fracture for a Piezoelectric Material Layer. Int. J. Sol. Structures50 (9), 1415–1424. 10.1016/j.ijsolstr.2013.01.013
38
WangE.CarcioneJ. M.YuanY.BaJ. (2021). Reflection of Inhomogeneous Plane Waves at the Surface of a Thermo-Poroelastic Medium. Geophys. J. Int.224 (3), 1621–1639.
39
WangE.CarcioneJ. M.BaJ.LiuY. (2020a). Reflection and Transmission of Plane Elastic Waves at an Interface between Two Double-Porosity media: Effect of Local Fluid Flow. Surv. Geophys.41 (2), 283–322. 10.1007/s10712-019-09572-6
40
WangZ.-W.FuL.-Y.WeiJ.HouW.BaJ.CarcioneJ. M. (2020b). On the Green Function of the Lord-Shulman Thermoelasticity Equations. Geophys. J. Int.220 (1), 393–403. 10.1093/gji/ggz453
41
WeiW.ZhengR.LiuG.TaoH. (2016). Reflection and Refraction of P Wave at the Interface between Thermoelastic and Porous Thermoelastic Medium. Transp Porous Med.113 (1), 1–27. 10.1007/s11242-016-0659-1
42
ZenerC. (1938). Internal Friction in Solids II. General Theory of Thermoelastic Internal Friction. Phys. Rev.53 (1), 90–99. 10.1103/physrev.53.90
43
ZenkourA. M.MashatD. S.AbouelregalA. E. (2013). The Effect of Dual-Phase-Lag Model on Reflection of Thermoelastic Waves in a Solid Half Space with Variable Material Properties. Acta Mechanica Solida Sinica26 (6), 659–670. 10.1016/s0894-9166(14)60009-4
Appendix A
Equations of Momentum Conservation and Heat Conduction
The strain (ϵ)-displacement (u) and constitutive relations (Biot, 1956) areandrespectively, where ϵ = ϵii, σij are the stress components, δij are the Kronecker components and fij are external stress forces, and the Einstein summation is asssumed.
The components of the equation of momentum conservation arewhere fi are the components of external body forces. Substituting the constitutive relations (A.2) into Eq. A.3 and using (A.1), in the absence of external body and stress forces, we obtainOn the other hand, the law of heat conduction iswhere T is the tempeature field and q a heat source.
Summary
Keywords
thermoelasticity and elasticity, reflection and transmission coefficients, energy partitions, inhomogeneous plane waves, attenuation angle
Citation
Hou W, Fu L-Y, Carcione JM and Han T (2022) Reflection and Transmission of Inhomogeneous Plane Waves in Thermoelastic Media. Front. Earth Sci. 10:850331. doi: 10.3389/feart.2022.850331
Received
07 January 2022
Accepted
29 March 2022
Published
26 April 2022
Volume
10 - 2022
Edited by
Jonas D. De Basabe, Center for Scientific Research and Higher Education in Ensenada (CICESE), Mexico
Reviewed by
Nantu Sarkar, University of Calcutta, India
Yaser Kiani, Shahrekord University, Iran
Updates
Copyright
© 2022 Hou, Fu, Carcione and Han.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Li-Yun Fu, lfu@upc.edu.cn; José M. Carcione, jose.carcione@gmail.com
This article was submitted to Solid Earth Geophysics, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.