ORIGINAL RESEARCH article

Front. Mech. Eng., 17 August 2021
Sec. Micro- and Nanoelectromechanical Systems
https://doi.org/10.3389/fmech.2021.704192

Boundary Reflections of Rolling Waves in Cubic Anisotropic Material

  • Department of Mechanical Engineering, University of Utah, Salt Lake City, UT, United States

Rolling waves have unconventional circular polarizations enabled by the equal-speed propagation of longitudinal and transverse waves in elastic solids. They can transport non-paraxial intrinsic (i.e. spin) mechanical angular momentum in the media. In this work, we analyze the rolling wave reflections and their effects on the non-paraxial spins in a cubic elastic half-space with an elastically supported boundary. Reflected waves from both normal and general oblique incidences are investigated. We show that, by adjusting the stiffness of the elastic boundary, we can precisely control the spin properties of the reflected waves, paving the way towards a broad category of spin manipulation techniques for bulk elastic waves.

1 Introduction

Recent advances in mechanical and acoustic metamaterials have made it possible to manipulate the general intrinsic spin angular momentum (Holanda et al., 2018; Zhu et al., 2018; An et al., 2020; Burns et al., 2020; Rückriegel et al., 2020; Shi and Yang, 2020) carried by elastic waves in the bulk (Long et al., 2018; Wang et al., 2018; Shi et al., 2019; Toftul et al., 2019; Long et al., 2020). Our previous work (Zhang et al., 2020) demonstrated that, when longitudinal and transverse elastic waves propagate at the same speed, rolling waves may emerge, carrying a spin that is orthogonal to, and in general not parallel to, the wave vector. This is only possible in anisotropic materials that satisfies certain equal-wave-speed criteria (Zhang et al., 2020). As the simplest anisotropic material characterized by only three independent elastic constants, cubic materials have been investigated in many studies (Stoneley, 1955; Borcherds, 1973). Thomas (1966) analyzed the longitudinal wave velocity at different directions using orientation components, and Mielnicki (1972) studied both longitudinal and transverse waves using approximation method. The bounds of elastic wave speeds (Zuo and Hjelmstad, 1997), surface wave polarizations (Chadwick and Wilson, 1992), and supersonic surface waves (Every, 2015) have all been the topics of an active line of research. Benefited from the fast development in additive manufacturing technology, cubic symmetric metamaterials can now provide us with a broader range of possible elastic constants (Favre et al., 2018; Tancogne-Dejean and Mohr, 2018; Lohmuller et al., 2019) even to the higher-order elasticity effects, such as the compression torsion metamaterial (Frenzel et al., 2017; Wang and Liu, 2020), hyper-stress tensor with second-gradient elasticity (Weeger, 2021), and chiral-micropolar metamaterials (Chen et al., 2021).

Although the stable propagation of rolling waves with non-paraxial spins have been demonstrated (Zhang et al., 2020), effects of boundary reflections on the wave spin angular momenta have yet to be investigated. The elastic wave reflection process is often both rich in physics and useful in applications (Schoenberg, 1980; Bedford and Drumheller, 1994; Vavrycuk and Psencík, 1998; Lee et al., 2017; Chen et al., 2019). However, studies on reflections from elastically supported boundaries of cubic materials are lacking. This further motivates us to analyze the reflection of rolling waves in a cubic material incident either normally or obliquely at such boundaries. The present work fills the gap in understanding and predicting rolling reflections and offers a new technique to manipulate the spin of elastic bulk waves.

In this paper, we first review the general governing equations for plane waves in an anisotropic medium. Next, the equal-wave-speed criteria in cubic materials are derived. Furthermore, we provide a general formulation for the spin angular momentum of bulk elastic waves. Lastly, the normal and oblique incident rolling waves are analyzed with different boundary conditions, where both spin-preserving and spin-flipping phenomena may emerge.

2 Plane Waves in Anisotropic Medium

With consideration of long-wavelength limit (i.e. quasi-static limit) only, we focus on the plane wave with wavelength much larger than any micro-structure size, which implies the dispersive behavior from geometry is negligible. The general bulk elastic wave equation without body force is,

σ=ρü,(1)

where σ is the stress tensor and u is the displacement vector.

The general solution of plane waves is

u=u0exp[i(kxωt)],(2)

where u0 denotes wave amplitude, and k the wave vector. We consider a plane wave propagating along an arbitrary direction specified by the unit vector,

k̂=l1e1+l2e2+l3e3,(3)

where l1, l2, l3 are direction cosines, and e1, e2, e3 are the standard orthonormal basis vectors of Cartesian coordinates. Then, we have

k=kk̂,(4)

where k is the wavenumber.

The time derivative of eq. 2 gives us,

tu=iωu.(5)

With the definition of a matrix L,

L=l1000l3l20l20l30l100l3l2l10.(6)

We have the spatial derivative in Voigt notation (Carcione, 2007),

=ikL.(7)

Combining the above equations, we arrive at

k2Γu=ρω2u,(k2Γijuj=ρω2ui)(8)

where

Γ=LCLT,(Γij=LiICIJLJj)(9)

is the Kelvin-Christoffel matrix.

Here, we consider the phase velocity v and its amplitude, wave speed v given as,

v=vk̂,v=ωk.(10)

Then, we obtain the governing equation for bulk elastic wave speeds in a general anisotropic material,

Γuρv2u=(Γρv2I)u=0.(11)

The expressions of elements in Γ are given below as,

Γ11=C11l12+C66l22+C55l32+2C56l2l3+2C15l1l3+2C16l1l2,Γ22=C66l12+C22l22+C44l32+2C24l2l3+2C46l1l3+2C26l1l2,Γ33=C55l12+C44l22+C33l32+2C34l2l3+2C35l1l3+2C45l1l2,Γ12=C16l12+C26l22+C45l32+(C46+C25)l2l3+(C14+C56)l1l3+(C12+C66)l1l2,Γ13=C15l12+C46l22+C35l32+(C45+C36)l2l3+(C13+C55)l1l3+(C14+C56)l1l2,Γ23=C56l12+C24l22+C34l32+(C44+C23)l2l3+(C36+C45)l1l3+(C25+C46)l1l2,(12)

where Cij are components of the stiffness matrix in Voigt notation (Carcione, 2007).

3 Constitutive Relation of Cubic Materials

Next, we focus our discussions on cubic materials characterized by the cubic crystallographic point groups. As illustrated in Figure 1, the symmetry conditions include four axes of three-fold rotational symmetry that can be identified with the body diagonals of a cubic unit cell (Bückmann et al., 2014; Dirrenberger et al., 2013; Authier, 2003). Rather counter-intuitively, this symmetry requirement does not include any axis of four-fold (90 degree) rotational symmetry (Ashcroft and Mermin, 1976, p. 121). The stiffness matrix of cubic material is given as (Authier, 2003; Norris, 2006; Zhang et al., 2020),

C(Cubic)=C11C12C12000C11C12000C11000C4400 SymC440C44,3constants(13)

FIGURE 1
www.frontiersin.org

FIGURE 1. The symmetry requirement for cubic materials: Body diagonals (orange solid lines) of a cubic unit cell serve as the four axes of three-fold rotational symmetry. They are the sufficient and necessary conditions for the cubic crystallographic point groups. On the contrary to some common misconceptions, cubic materials do not need to have any four-fold rotational symmetry.

For waves propagating along the z-direction, the direction cosines are

l1=0,l2=0,l3=1.(14)

Then, the Kelvin-Christoffel matrix components become,

Γ11=Γ22=C44,Γ33=C11,Γ12=Γ13=Γ23=0.(15)

Therefore, the governing eq. 11 becomes,

0=(Γρv2I)u=C44ρv2000C44ρv2000C11ρv2u1u2u3.(16)

We note that u3 represents the longitudinal wave, while u1 and u2 represent the transverse waves. Therefore, from eq. 16, we obtain the longitudinal and transverse wave speeds as,

vL=C11ρ,vT=C44ρ.(17)

Thus, for waves traveling in the z-direction of a cubic material, in order to have a rolling wave with non-praxial spin polarization (Zhang et al., 2020), we need the equal-speed criterion:

C11=C44.(18)

By symmetry, for the waves propagating in x− or y−direction, the equal-speed criterion is exactly the same as eq. 18.

4 Intrinsic Spin Angular Momentum of Bulk Elastic Waves

To calculate the spin density of the bulk elastic waves, we consider the displacement field of a general plane wave u=ũexp[i(kxωt)] with

ũ=mnl.(19)

Importantly, here m, n and l are complex-valued, so that they contain the information not only about amplitudes but also about the relative phase differences among displacement components of the traveling wave. The spin angular momentum density, as a real-valued vector, can be calculated as (Berry, 2009; Long et al., 2018; Burns et al., 2020):

s=ρω2<ũ|Ŝ|ũ>=ρω2Im[ũ*×ũ],(20)

where (⋅)* denotes complex conjugation, and the spin-1 operator is a third-order tensor defined as

Ŝ=00000i0i0,00i000i00,0i0i00000.(21)

Hence, the spin density for a general traveling wave is

s=ρωImn*ll*mm*n.(22)

We note that the above formulation makes physical sense only if all displacement components propagate at the same wave speed. Particularly for cubic materials, the equal-speed criterion eq. 18 needs to be true. In this case, we can have stable propagation of spins pointing at any direction by adjusting the complex amplitudes m, n and l. This can be realized by simply changing the relative phase differences. Hence, the direction of the spin vector s can be completely independent from the direction of the wave vector k.

5 Normal Incidence and Reflection

We now investigate normal reflections of a rolling wave at a general elastic boundary. This serves as an example non-paraxial spin manipulation. Considering a rolling wave along the z-direction normally incident on a flat surface, we have the instantaneous wave displacement fields at time t = 0 as uI=ũIexp(ikz) and uR=ũRexp(ikz) with

ũI=mInIlIandũR=mRnRlR,(23)

where the superscripts, I and R, denote the incident and reflected waves, respectively.

For the displacement with small amplitude, the linear strain terms are dominant. In this case, by neglecting higher-order strain terms, the strains can be calculated from displacements by

εij=12(ui,j+uj,i),i,j=x,y,z(24)

where the comma , in subscripts denotes the derivative operation.

The reflection occurs at z = 0. So we have eikz = e−ikz = 1, and the strain vector in Voigt notation becomes

εI=εxxIεyyIεzzI2εyzI2εxzI2εxyI=00lIiknIikmIik0,εR=εxxRεyyRεzzR2εyzR2εxzR2εxyR=00lRiknRikmRik0.(25)

The cubic constitutive relations are

σxxσyyσzzσyzσxzσxy=C11C12C12000C12C11C12000C12C12C11000000C44000000C44000000C44εxxεyyεzz2εyz2εxz2εxy.(26)

It is easy to compute the stress components for incident and reflection waves, respectively:

σxxIσyyIσzzIσyzIσxzIσxyI=C12lIikC12lIikC11lIikC44nIikC44mIik0,σxxRσyyRσzzRσyzRσxzRσxyR=C12lRikC12lRikC11lRikC44nRikC44mRik0.(27)

For a elastically supported cubic half-space, the boundary conditions are

σzx0=Kxux0,σzy0=Kyuy0,σzz0=Kzuz0.(28)

where (Kx, Ky, Kz) are components of the distributed stiffness per unit area (Zhang et al., 2017) representing a general elastic foundation supporting the solid surface, and the total stress and displacement are

σzx0=σzxI+σzxR,σzy0=σzyI+σzyR,σzz0=σzzI+σzzR.(29)
ux0=uxI+uxR,uy0=uyI+uyR,uz0=uzI+uzR.(30)

Substituting the displacement and stress components into elastic boundary conditions, we get

C44mIikC44mRik=Kx(mI+mR),(31)
C44nIikC44nRik=Ky(nI+nR),(32)
C11lIikC11lRik=Kz(lI+lR).(33)

Solving the above system of equations, we obtain

mR=C44ikKxC44ik+KxmI,(34)
nR=C44ikKyC44ik+KynI,(35)
lR=C11ikKzC11ik+KzlI.(36)

Because of the equal-speed requirement of rolling wave (C11 = C44), the amplitude of reflection wave becomes,

mR=RxmI,nR=RynI,lR=RzlI,(37)

with the generalized amplitude ratio of normal reflection given as

Rj=C11ikKjC11ik+Kj,j=x,y,z.(38)

We can also write in the form of a reflection matrix (Zoeppritz, 1919; Ursin and Tjäland, 1996):

mRnRlR=Rx000Ry000RzmInIlI.(39)

with the reflection matrix R̂=diag(Rx,Ry,Rz). This complex-valued non-dimensional parameter Rj is the key connection between the incident and reflected waves and deserves further analysis. We note that |Rj| = 1 is consistent with the fact that all wave energy is reflected. Hence, we have

Rj=eiϕwithtanϕ=2Kj/C11k1(Kj/C11k)2.(40)

Here, the phase angle ϕ of Rj corresponds to the phase change during the reflection, and it dependents on the boundary-bulk stiffness ratio Kj/C11 and the wave number k, as shown in Figure 2. Therefore, by adjusting the elastic stiffness at the boundary, one can manipulate the spin of the reflected waves.

FIGURE 2
www.frontiersin.org

FIGURE 2. The complex-valued amplitude ratio of normal reflection, Rj, is determined by the boundary-bulk stiffness ratio, Kj/C11, and wave number k.

Although all properties of both the bulk and the reflection surface are assumed to be independent of the incidence wave frequency, here the reflection phase change can be frequency-dependent, as the angular wave number k appears in the ratio. The emergence of frequency dependency can be intuitively explained by the role of wavelength during reflection. We note from eq. 25 that k first appears in the strain calculations since, for a given wave displacement amplitude, the strains in the propagation direction, ɛzj, are actually inversely proportional to the wavelength: Longer wavelength gives rise to a smaller strain and vice versa. Consequently, the stresses, σzj, and the force acting on the boundary springs are wavelength-dependent as well. If the boundary is supported by an elastic foundation with finite stiffness per area, Kj, we have the following: At the low-frequency and long-wavelength limit, the force per area acting on the boundary, |σzj|∝ C11k → 0. So, the boundary hardly moves, and the incidence wave effectively “sees” a rigid surface; At the high-frequency and short-wavelength limit, the force per area acting on the boundary, |σzj|∝ C11k. So, the boundary moves a lot, and the incidence wave effectively “sees” a free surface.

Next, we consider some special cases and focus on the in-xz-plane rolling waves with nI = nR = 0. We note that the elastic boundary condition degenerates into traction-free boundary condition (Neumann type) when the boundary stiffness Kj → 0. Then, the amplitudes of reflected wave become

mR=mI,lR=lI.(41)

Similarly, the elastic boundary condition degenerates into the fixed rigid boundary condition (Dirichlet type) when the stiffness Kj. Then, the amplitudes of reflected wave become

mR=mI,lR=lI.(42)

In both cases above, by the definition of spin given in eq. 22, we can conclude sR = sI, as shown in Figures 3B,C. So, the spin is unaffected by the reflection process.

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) Rolling wave reflections with normal incidence at a general elastic boundary. In all sub-plots, the incidence wave carries a spin angular momentum, syI=1. Here sI and sR indicate the spin density of incidence wave and reflection wave, respectively. They can be calculated using the eq. 22. The reflection process is (B) spin-preserving at the fixed rigid boundary, (C) spin-preserving at the free boundary, (D) spin-flipping at the fixed-free boundary, and (E) spin-flipping at the free-fixed boundary.

In addition, for the fixed-free hybrid boundary (Kx and Kz → 0, Figure 3D), we have

mR=mI,lR=lIsR=sI.(43)

Similarly, for the free-fixed hybrid boundary (Kx → 0 and Kz, Figure 3E), we have

mR=mI,lR=lIsR=sI.(44)

Thus, both hybrid boundaries will flip the spin for any incident rolling wave.

Summarizing different cases for the normal incidence, we note that, for Kx → 0 and Kz → 0, the boundary becomes traction-free (Neumann type) and we have ũR=ũI with no phase change. For Kx and Kz, the boundary becomes rigid (Dirichlet type) and we have an out-of-phase reflected wave with ũR=ũI. In both cases, we have sR = sI, so the spin is preserved. In contrast, for hybrid boundaries (Kx, Kz → 0) and (Kx → 0, Kz), we get sR = −sI, so the spin is flipped due to the difference in phase change between the longitudinal and transverse displacement components during the reflection process.

Figure 3 illustrates the results using the example of an incident rolling wave carrying a non-paraxial spin of syI=1. The reflection process is spin-preserving in both rigid and free boundaries, while being spin-flipping for both hybrid free-fixed and hybrid fixed-free boundary conditions.

6 The Oblique Incidence and Reflection at 45 degree

Next we investigate the scenario in which the boundary plane is oblique to the incident rolling wave at the 45° angle. We focus on the incident rolling wave propagating towards the positive z−direction,

uI=mInIlIeikIz(45)

As shown in Figure 4A, the incidence angle here is αI = 45°. Importantly, we note that, due to material anisotropy, the reflection angle(s) will not necessarily be equal to the incidence angle. In fact, a single incidence wave, in general, may result in multiple reflected waves towards different directions (Achenbach, 2012; Graff, 2012). Hence, the reflection angle(s) need to be determined by a detailed analysis on boundary conditions and material properties. In this case, we can write the reflection wave as a general form in xz-plane,

uR=mRnRlReikR[cos(θ)x+sin(θ)z](46)

where θ is the angle of wave vector with respect to the x-axis. It follows that the reflection angle can be expressed as αR = θ + 45°, θ ∈ (−45°, + 45°).

FIGURE 4
www.frontiersin.org

FIGURE 4. (A) Illustration of rolling wave reflections at 45° incidence. In all sub-plots, the incidence wave carries a spin angular momentum, syI=+1. Here sI and sR indicate the spin density of incidence wave and reflection wave, respectively. They can be calculated using the eq. 22. The reflection process is (B) spin-preserving at the fixed rigid boundary, (C) spin-preserving at the free boundary, (D) spin-flipping at the fixed-free boundary, and (E) spin-flipping at the free-fixed boundary. Note that the free boundary in (C) has the specially effect of giving rise to an additional linearly polarized reflected wave. Although the spin is still exactly preserved, the overall reflected wave has an elliptical rolling polarization.

Next, we look into the general conditions for rolling waves propagating in the xz-plane. The normalized wave vector, as given in eq. 3, is now characterized by cosine directions l2 = 0 and l12+l32=1. The element expressions of matrix Γ become,

Γ11=C11l12+C44l32,Γ22=C44,Γ33=C44l12+C11l32,Γ12=0,Γ13=(C12+C44)l1l3,Γ23=0.(47)

Therefore, the governing eq. 11 becomes,

0=(Γρv2I)u=C11l12+C44l32ρv20(C12+C44)l1l30C44ρv20(C12+C44)l1l30C44l12+C11l32ρv2u1u2u3.(48)

Combining with the equal speed criterion in eq. 18, i.e., C11 = C44, we obtain

0=(Γρv2I)u=C11ρv20(C12+C11)l1l30C11ρv20(C12+C11)l1l30C11ρv2u1u2u3.(49)

This shows that the in-plane waves (u1, u3) are in general coupled together, and the out-of-plane wave u2 is decoupled from the two in-plane polarizations.

For non-trivial solution to exist for in-plane wave velocities in eq. 49, the determinant must vanish. This leads to the following

ρ2v42C11ρv2+C112(C12+C11)2l12l32=0.(50)

This is a quadratic polynomial with respect to v2. It has the discriminant as Δ=4ρ2(C12+C11)2l12l320, and the roots are

v2=C11±|(C12+C11)l1l3|ρ.(51)

For any linear material to be statically stable, it needs a positive strain energy function, and hence a positive-definite stiffness matrix (Ting, 1996), which leads to the following conditions:

C11>0andC11<C12<C11(52)

Then, in order to determine the reflection angle(s) in this case, we apply Snell’s law for reflection

ωIvIsin(αI)=ωRvRsin(αR).(53)

Here, the incidence angle is αI = 45°. Combining it with ωI = ωR and vI=C11/ρ, we arrive at

12C11=sin(αR)C11±|C12+C11|l1l3=sin(θ+45°)C11±|C12+C11|l1l3(54)

with the cosine directions,

l1=cos(θ),l3=sin(θ).(55)

Further simplification gives us,

sin(θ)cos(θ)(2C11±|C11+C12|)=0.(56)

For θ ∈ (−45°, + 45°) with the positive-definite conditions of eq. 52, the only solution to above equation is

θ=0(57)

which implies αR = 45°. Therefore, there can be one reflection wave vector only, and the reflected wave will propagate along x-direction at the speed vR = vI. Hence, we can rewrite the reflection wave as,

uR=mRnRlReikRxwithkR=kI.(58)

For brevity, denoting the wave number as kI = kR = k, we obtain

uI=mInIlIeikz,uR=mRnRlReikx.(59)

The corresponding strain wave amplitudes are

εI=εxxIεyyIεzzI2εyzI2εxzI2εxyI=00lIiknIikmIik0,εR=εxxRεyyRεzzR2εyzR2εxzR2εxyR=mRik000lRiknRik.(60)

And the stress wave amplitudes are

σxxIσyyIσzzIσyzIσxzIσxyI=C12lIikC12lIikC11lIikC11nIikC11mIik0,σxxRσyyRσzzRσyzRσxzRσxyR=C11mRikC12mRikC12mRik0C11lRikC11nRik.(61)

To analyze the normal and transverse elastic boundary condition, we introduce a local coordinate system x′ = Ψx specified by the transformation matrix,

Ψ=cos(φ)0sin(φ)010sin(φ)0cos(φ)=22101020101(62)

where the φ = 45° as shown in Figure 4.

The general elastic boundary conditions in the local coordinate system are

σzx0=Kxux0,σzy0=Kyuy0,σzz0=Kzuz0.(63)

where (Kx, Ky, Kz) are components of distributed stiffness per unit area, and

σzx0=σzxI+σzxR,σzy0=σzyI+σzyR,σzz0=σzzI+σzzR.(64)
ux0=uxI+uxR,uy0=uyI+uyR,uz0=uzI+uzR.(65)

are components of total stress and total displacement on the boundary surface.

We next perform coordinate transforms using

σ=ΨσΨT,u=Ψu.(66)

In the local coordinates, x′, we obtain

σI=ik2(C11+C12)lI+2C11mI2C11nI(C11C12)lI2C11nI2C12lI2C11nI(C11C12)lI2C11nI(C11+C12)lI2C11mI,(67)
uI=22lI+mI2nIlImI.(68)

Similarly, the local stress and displacement of reflection wave are,

σR=ik2(C11+C12)mR+2C11lR2C11nR(C12C11)mR2C11nR2C12mR2C11nR(C12C11)mR2C11nR(C11+C12)mR2C11lR,(69)
uR=22lR+mR2nRlRmR.(70)

Substituting the stress and displacement components eqs  67 –eqs  70 into the boundary conditions eqs  64,65, we obtain

ik2[(C11C12)lI+(C12C11)mR]=Kx22(lI+mI+lR+mR),(71)
ik22[C11nIC11nR]=Ky(nI+nR),(72)
ik2[(C11+C12)lI2C11mI+(C11+C12)mR2C11lR]=Kz22(lImI+lRmR).(73)

Then, we can solve the system of eqs 71–73 for the complex-valued amplitude components of the reflected wave:

mR=Kx(κ+Kz)Kz+2C11ik2(κKx)lI+2KxKzmIKx(Kz+κ+)+Kz+2C11ik2(κ+Kx),(74)
nR=ik2C11Kyik2C11+KynI,(75)
lR=2(κ+κKxKz)lI+Kz2C11ik2(κ+Kx)Kx(κ++Kz)mIKx(Kz+κ+)+Kz+2C11ik2(κ+Kx),(76)

where

κ+=ik2(C11+C12),κ=ik2(C11C12).(77)

We can also write in the form of a reflection matrix (Zoeppritz, 1919; Ursin and Tjäland, 1996).

mRnRlR=Rmm0Rml0Rnn0Rlm0RllmInIlI.(78)

where

Rmm=2KxKzDRml=KzκKxκ++2C11ik(κKx)DRnn=ik2C11Kyik2C11+KyRlm=KzκKxκ+2C11ik(κ+Kx)DRll=2(κ+κKxKz)D(79)

with denominator D=Kx(Kz+κ+)+(Kz+2C11ik2)(κ+Kx). We note that RmlRlm, which means that the complex-valued reflection matrix is not Hermitian. A closer scrutiny reveals that the amplitudes of Rml and Rlm are always the same. The two complex-valued reflective coefficients only differ in their phase angles.

We also note that the term Rnn for the out-of-xz-plane wave is independent and very similar to the normal incidence scenario, as shown in eqs 38,40 and Figure  2. All other terms in eq. 79 are illustrated in Figure  5 to show how they change with the wave number, k, and the boundary stiffness ratio, Kz/Kx. We observe in Figure  5B that both Rml and Rlm shows a phase jump of 180°. This is because of two reasons: 1) The small wave number k = 0.001 leads to the fact that both Rml and Rlm are almost purely imaginary; and 2) At the boundary stiffness ratio of Kz/Kx = 3, values transition from a negative imaginary number to a positive imaginary number. This results in the 180° phase jump.

FIGURE 5
www.frontiersin.org

FIGURE 5. The amplitude and phase of reflection matrix elements under different wave numbers k at Cr = C12/C11 = 0. (A,B) Amplitude and phase at k = 0.001, (C,D) Amplitude and phase at k = 1, (E,F) Amplitude and phase at k = 1,000.

Clearly, the out-of-xz-plane wave nR only depends on nI:

nR=nIforKy=0,(80)
nR=nIforKy=.(81)

The four extreme cases for in-xz-plane wave polarizations are:

Case I (Figure 4B), Kx, Kz:

mR=mI,lR=lI,sR=sI.(82)

Case II (Figure 4C), Kx → 0, Kz → 0:

mR=lI,lR=mI+C11+C12C11lI=mI+1+CrmRsR=sI(83)

with Cr = C12/C11. Here, the second term in lR is associated with lI, and mR = lI. So, this term does not have any contribution to the circularly rolling polarization, and it gives rise to an additional linearly polarized wave. Consequently, the spin is still exactly preserved, but the overall reflected wave has a elliptically rolling polarization.

Case III (Figure 4D), Kx, Kz → 0:

mR=lI,lR=mIsR=sI.(84)

Case IV (Figure 4E), Kx → 0, Kz:

mR=lI,lR=mIsR=sI.(85)

As illustrated in Figure  4C, the results in eq. 83 for Case II are special and remarkable. The reflected wave from the free boundary is different from all other cases. However, the reflected spin density is indeed still the same to that of the incidence wave. The reason is that the additional term in lR does not contribute to the circularly rolling polarization. To further investigate this especially interesting case, we calculate the polarization of reflected waves with different ratios of Cr = C12/C11 and plot the result in Figure 6. Note that the ratios of Cr = −1 in Figure 6A and Cr = 1 in Figure 6E are actually extreme scenarios not attainable in statically stable materials due the positive definite conditions listed in eq. 52.

FIGURE 6
www.frontiersin.org

FIGURE 6. In the oblique reflection at 45°, elliptically rolling polarizations of the reflected waves from the free boundary may appear due to different material stiffness constants C11 and C12. An incident spin which is the same of that shown (A) is assumed throughout. (A) also represents the purely circularly rolling polarization for Cr = −1. (B–F) show the elliptically rolling polarizations for the cases of Cr = −0.5, Cr = 0, Cr = 0.5 and Cr = 1, respectively. Note that all polarizations are scaled for the purpose of clearly demonstrating differences in elliptical eccentricity.

7 The Oblique Incidence and Reflection at Arbitrary Angle

Next, we look into the situation with an arbitrary incidence angle. From Snell’s law eq. 53, we can obtain

sin(αI)C11=sin(αR)C11±(C12+C11)l1l3(86)

where l1 = cos(θ), l3 = sin(θ). Noticed from Figure 4A that the incidence angle is identical to the rotation angle, i.e., αI = φ and then the reflection angle αR = π/2 − φ + θ. The reflection angle αR is the solution of the transcendental equation depending on the incidence angle αI and the material constant ratio Cr = C12/C11:

1±(1+Cr)sin(αR+αIπ/2)cos(αR+αIπ/2)sin(αI)=sin(αR)(87)

where “±” gives rise to two different reflected wave modes. In general, those two modes have two different reflection angles.

Next, we numerically solve eq. 87 for the reflection angle αR and illustrate the results in Figure 7. The points P1 and P2 in the plots indicate the cases of the normal incidence and the 45° incidence, respectively. Both cases have been thoroughly discussed in previous sections. The next critical point P3 is where the reflection angle αRπ/2 for Mode 1. Plugging αR = π/2 into eq. 87, we obtain:

(1+Cr)2sin6(α3I)+sin2(α3I)1=0.(88)

FIGURE 7
www.frontiersin.org

FIGURE 7. The reflection angle αR of two modes with respect to the incidence angle αI. (A) Cr = 0.9, (B) Cr = 0.1, (C) Cr = 0, (D) Cr = −0.1, (E) Cr = −0.149452, (F) Cr = −0.9.

Thus, we arrive at the analytical solution of the above equation,

α3I=arcsin12(1+Cr)2+Δ3+12(1+Cr)2Δ3(89)

with

Δ=14(1+Cr)4+127(1+Cr)6>0.(90)

We know that, with Δ > 0, the cubic function always has one real root and two complex conjugate roots. In this case, we can only obtain one critical angle α3I. Therefore, we conclude that, for incidence angle larger than the critical angle α3I, the reflection wave of Mode 1 will no longer exist.

Here, from eq. 52 we have the constraint −1 < Cr = C12/C11 < 1. Hence, we investigate all values of Cr within this range and find the following four representative scenarios for Mode 2 of the reflected waves:

• For 0 < Cr < 1, as shown in Figures  7A,B, we have one local maximum (point P4) and no local minimum for Mode 2. The reflection angle αR will increase together with the incidence angle αI until it reaches to the maximum at point P4. Then, the anomalous trend appears, and the reflection angle αR will decrease if we further increase the incidence angle αI. Consequently, for Mode 2, we may have two different incidence angles with the same reflection angle.

• For −0.149 < Cr ≤ 0, as shown in Figures  7C,D, we have one local maximum (point P4) and one local minimum (point P5). Remarkably, as shown in the inset of Figure 7D, we may have three different incidence angles with the same reflection angle for Mode 2. Also, we note that, with the decreasing value of Cr, point P4 and point P5 move closer to each other.

• For Cr = −0.149, as shown in Figure 7E, we have no local maximum, no local minimum, but a saddle point (point P6). This saddle point P6 is the result of point P4 and point P5 merging together.

• For −1 < Cr < −0.149452 as shown in Figure 7F, we have a monotonic curve for Mode 2. The reflection angle αR always increases with the increasing incidence angle αI. And this is the normal phenomenon generally understood.

In the first two scenarios listed above, we have the anomalous trend phenomenon: The reflection angle αR does not have a monotonic relationship with the incidence angle αI. This may provide us with a new technique to control the direction of reflected waves.

Furthermore, we can take the derivative of both sides of eq. 87 with respect to the incidence angle αI and obtain the changing rate of reflection angle αR of Mode 2,

dαRdαI=11+Cr2sin(2αR+2αI)sin(2αI)(1+Cr)sin2(αI)cos(2αR+2αI)sin(2αR)+(1+Cr)sin2(αI)cos(2αR+2αI)(91)

We note that this derivative of αR with respect to the αI is the root of another transcendental equation and rarely has closed-form solution. Thus, we solve it numerically and plot the results in Figure 8. The same four represented scenarios discussed above can be easily categorized here as different number of roots of eq. 91. The results shown in Figures 7, 8 are consistent with each other.

FIGURE 8
www.frontiersin.org

FIGURE 8. The derivative of reflection angle αR with respect to the incidence angle αI. (A) Whole range of Cr, (B) Zoom-in range of Cr with two zero roots.

8 Conclusion

Cubic materials are anisotropic elastic materials that satisfy the symmetry requirements illustrated in Figure 1, and their stiffness matrix is given by eq. 13. If the equal-speed criteria for longitudinal and transverse waves, as given in eq. 18, also holds in a cubic material, then we can have stable propagation of any spin angular momentum carried by bulk elastic waves. In this paper, we focus on rolling waves, which have a spin vector orthogonal to its wave vector, in cubic materials. We investigate the reflections of rolling waves incident either normally or obliquely at an elastically supported flat boundary. Behaviors of spin-preserving, spin-flipping, and general arbitrary phase shifts are shown by adjusting the boundary conditions. This work provides a new approach to manipulating bulk elastic wave spins. Several other key findings are summarized below:

1) For the normal incidence of rolling waves, the reflection is spin-preserving at both the fixed rigid boundary and the free boundary, while being spin-flipping at both the fixed-free and free-fixed boundaries. Most generally, we can even introduce the arbitrary phase shifts between reflection and incidence waves by adjusting the elastic boundary conditions.

2) For oblique incidence at 45°, we prove that, contrary to the general anisotropic case, there will be only 1 reflected wave with a unique reflection angle of 45°.

3) With the 45° incidence, the reflection is spin-preserving at a fixed rigid boundary, while being spin-flipping at both fixed-free and free-fixed boundaries. The free boundary is a uniquely special case here: While the reflection is still spin-preserving, depending on the cubic material constant ratio, C12/C11, it will also introduce an additional linearly polarized part. Hence, an incidence wave with circularly rolling polarization will result in a reflected wave with elliptically rolling polarization.

4) For an arbitrary incidence angle, we showed that each incidence wave will in general result in two reflected wave modes with different reflection angles. One of the mode (Mode 1 in Figure 7) exists only when the incidence angle is less than a critical value α3I. Also, The other mode (Mode 2 in Figure 7) may show anomalous trend.

5) We developed a new way to manipulate the elastic wave spin. Further, the present study can be extended to more general cases, i.e., arbitrary incidence angle for reflection and refraction problem, the reflection matrix properties of arbitrary incidence angle, visco-elastically supported boundary condition, or breaking the limit of statically stable materials by using metamaterials.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

PZ is the primary contributor to this study and this manuscript. PW is the secondary contributor to this study and this manuscript.

Funding

This work was supported by start-up funds of the Dept. Mechanical Engineering at Univ. Utah.

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.

Acknowledgments

The support and resources from the Center for High Performance Computing at Univ. Utah are gratefully acknowledged.

References

Achenbach, J. (2012). Wave Propagation in Elastic Solids. Amsterdam: Elsevier.

An, K., Litvinenko, A., Kohno, R., Fuad, A., Naletov, V., Vila, L., et al. (2020). Coherent Long-Range Transfer of Angular Momentum Between Magnon Kittel Modes by Phonons. Phys. Rev. B 101, 060407. doi:10.1103/physrevb.101.060407

CrossRef Full Text | Google Scholar

Ashcroft, N. W., and Mermin, N. D. (1976). Solid State Physics. New York London): Holt, Rinehart and Winston.

Authier, A. (2003). International Tables for Crystallography: Volume D: Physical Properties of Crystals. Wiley Online Library.

Bedford, A., and Drumheller, D. (1994). Elastic Wave Propagation. Amsterdam: John Wileg g Sons, 151–165.

Berry, M. V. (2009). Optical Currents. J. Opt. A: Pure Appl. Opt. 11, 094001. doi:10.1088/1464-4258/11/9/094001

CrossRef Full Text | Google Scholar

Borcherds, P. H. (1973). The Scattering of Light by Elastic Waves in Cubic Crystals. Optica Acta Int. J. Opt. 20, 147–159. doi:10.1080/713818743

CrossRef Full Text | Google Scholar

Bückmann, T., Schittny, R., Thiel, M., Kadic, M., Milton, G. W., and Wegener, M. (2014). On Three-Dimensional Dilational Elastic Metamaterials. New J. Phys. 16, 033032. doi:10.1088/1367-2630/16/3/033032

CrossRef Full Text | Google Scholar

Burns, L., Bliokh, K. Y., Nori, F., and Dressel, J. (2020). Acoustic versus Electromagnetic Field Theory: Scalar, Vector, Spinor Representations and the Emergence of Acoustic Spin. New J. Phys. 22, 053050. doi:10.1088/1367-2630/ab7f91

CrossRef Full Text | Google Scholar

Carcione, J. M. (2007). Wave fields in Real media: Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic media. Amsterdam: Elsevier.

Chadwick, P., and Wilson, N. (1992). The Behaviour of Elastic Surface Waves Polarized in a Plane of Material Symmetry. Iii. Orthorhombic and Cubic media. Proc. R. Soc. Lond. A. 438, 225–247. doi:10.1098/rspa.1992.0104

CrossRef Full Text | Google Scholar

Chen, Y., Kadic, M., and Wegener, M. (2021). Chiral Triclinic Metamaterial Crystals Supporting Isotropic Acoustical Activity and Isotropic Chiral Phonons. Proc. R. Soc. A. 477, 20200764. doi:10.1098/rspa.2020.0764

CrossRef Full Text | Google Scholar

Chen, Y., Li, X., Nassar, H., Norris, A. N., Daraio, C., and Huang, G. (2019). Nonreciprocal Wave Propagation in a Continuum-Based Metamaterial with Space-Time Modulated Resonators. Phys. Rev. Appl. 11, 064052. doi:10.1103/physrevapplied.11.064052

CrossRef Full Text | Google Scholar

Dirrenberger, J., Forest, S., and Jeulin, D. (2013). Effective Elastic Properties of Auxetic Microstructures: Anisotropy and Structural Applications. Int. J. Mech. Mater. Des. 9, 21–33. doi:10.1007/s10999-012-9192-8

CrossRef Full Text | Google Scholar

Every, A. G. (2015). Supersonic Surface Acoustic Waves on the 001 and 110 Surfaces of Cubic Crystals. The J. Acoust. Soc. America 138, 2937–2944. doi:10.1121/1.4934557

CrossRef Full Text | Google Scholar

Favre, J., Lohmuller, P., Piotrowski, B., Kenzari, S., Laheurte, P., and Meraghni, F. (2018). A Continuous Crystallographic Approach to Generate Cubic Lattices and its Effect on Relative Stiffness of Architectured Materials. Additive Manufacturing 21, 359–368. doi:10.1016/j.addma.2018.02.020

CrossRef Full Text | Google Scholar

Frenzel, T., Kadic, M., and Wegener, M. (2017). Three-Dimensional Mechanical Metamaterials with a Twist. Science 358, 1072–1074. doi:10.1126/science.aao4640

PubMed Abstract | CrossRef Full Text | Google Scholar

Graff, K. F. (2012). Wave Motion in Elastic Solids. Amsterdam: Courier Corporation.

Holanda, J., Maior, D. S., Azevedo, A., and Rezende, S. M. (2018). Detecting the Phonon Spin in Magnon-Phonon Conversion Experiments. Nat. Phys 14, 500–506. doi:10.1038/s41567-018-0079-y

CrossRef Full Text | Google Scholar

Lee, H. J., Lee, J. R., Moon, S. H., Je, T. J., Jeon, E. C., Kim, K., et al. (2017). Off-Centered Double-Slit Metamaterial for Elastic Wave Polarization Anomaly. Sci. Rep. 7, 1–13. doi:10.1038/s41598-017-15746-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Lohmuller, P., Favre, J., Kenzari, S., Piotrowski, B., Peltier, L., and Laheurte, P. (2019). Architectural Effect on 3d Elastic Properties and Anisotropy of Cubic Lattice Structures. Mater. Des. 182, 108059. doi:10.1016/j.matdes.2019.108059

CrossRef Full Text | Google Scholar

Long, Y., Ge, H., Zhang, D., Xu, X., Ren, J., Lu, M.-H., et al. (2020). Symmetry Selective Directionality in Near-Field Acoustics. Natl. Sci. Rev. 7, 1024–1035. doi:10.1093/nsr/nwaa040

CrossRef Full Text | Google Scholar

Long, Y., Ren, J., and Chen, H. (2018). Intrinsic Spin of Elastic Waves. Proc. Natl. Acad. Sci. USA 115, 9951–9955. doi:10.1073/pnas.1808534115

PubMed Abstract | CrossRef Full Text | Google Scholar

Mielnicki, J. (1972). Elastic Waves in {100}, {110}, and {111} Planes of Cubic Crystals. IEEE Trans. Son. Ultrason. 19, 15–17. doi:10.1109/t-su.1972.29636

CrossRef Full Text | Google Scholar

Norris, A. N. (2006). Poisson's Ratio in Cubic Materials. Proc. R. Soc. A. 462, 3385–3405. doi:10.1098/rspa.2006.1726

CrossRef Full Text | Google Scholar

Rückriegel, A., Streib, S., Bauer, G. E. W., and Duine, R. A. (2020). Angular Momentum Conservation and Phonon Spin in Magnetic Insulators. Phys. Rev. B 101, 104402. doi:10.1103/physrevb.101.104402

CrossRef Full Text | Google Scholar

Schoenberg, M. (1980). Elastic Wave Behavior Across Linear Slip Interfaces. J. Acoust. Soc. America 68, 1516–1521. doi:10.1121/1.385077

CrossRef Full Text | Google Scholar

Shi, C., Zhao, R., Long, Y., Yang, S., Wang, Y., Chen, H., et al. (2019). Observation of Acoustic Spin. Natl. Sci. Rev. 6, 707–712. doi:10.1093/nsr/nwz059

CrossRef Full Text | Google Scholar

Shi, X., and Yang, J. (2020). Spin-1 Weyl Point and Surface Arc State in a Chiral Phononic Crystal. Phys. Rev. B 101, 214309. doi:10.1103/physrevb.101.214309

CrossRef Full Text | Google Scholar

Stoneley, R. (1955). The Propagation of Surface Elastic Waves in a Cubic Crystal. Proc. R. Soc. Lond. A. 232, 447–458. doi:10.1098/rspa.1955.0230

CrossRef Full Text | Google Scholar

Tancogne-Dejean, T., and Mohr, D. (2018). Elastically-Isotropic Elementary Cubic Lattices Composed of Tailored Hollow Beams. Extreme Mech. Lett. 22, 13–18. doi:10.1016/j.eml.2018.04.005

CrossRef Full Text | Google Scholar

Thomas, T. Y. (1966). Longitudinal Waves in Cubic Crystals. Proc. Natl. Acad. Sci. 55, 1011–1014. doi:10.1073/pnas.55.5.1011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ting, T. C. T. (1996). Positive Definiteness of Anisotropic Elastic Constants. Maths. Mech. Sol. 1, 301–314. doi:10.1177/108128659600100302

CrossRef Full Text | Google Scholar

Toftul, I. D., Bliokh, K. Y., Petrov, M. I., and Nori, F. (2019). Acoustic Radiation Force and Torque on Small Particles as Measures of the Canonical Momentum and Spin Densities. Phys. Rev. Lett. 123, 183901. doi:10.1103/physrevlett.123.183901

PubMed Abstract | CrossRef Full Text | Google Scholar

Ursin, B., and Tjåland, E. (1996). The Information Content of the Elastic Reflection Matrix. Geophys. J. Int. 125, 214–228. doi:10.1111/j.1365-246x.1996.tb06547.x

CrossRef Full Text | Google Scholar

Vavrycuk, V., and Psencík, I. (1998). Pp-Wave Reflection Coefficients in Weakly Anisotropic Elastic media. Geophysics 63, 2129–2141.

Google Scholar

Wang, L., and Liu, H.-T. (2020). 3D Compression-Torsion Cubic Mechanical Metamaterial with Double Inclined Rods. Extreme Mech. Lett. 37, 100706. doi:10.1016/j.eml.2020.100706

CrossRef Full Text | Google Scholar

Wang, S., Ma, G., and Chan, C. T. (2018). Topological Transport of Sound Mediated by Spin-Redirection Geometric Phase. Sci. Adv. 4, eaaq1475. doi:10.1126/sciadv.aaq1475

PubMed Abstract | CrossRef Full Text | Google Scholar

Weeger, O. (2021). Numerical Homogenization of Second Gradient, Linear Elastic Constitutive Models for Cubic 3d Beam-Lattice Metamaterials. Int. J. Sol. Structures 224, 111037. doi:10.1016/j.ijsolstr.2021.03.024

CrossRef Full Text | Google Scholar

Zhang, P., Kern, C., Sun, S., Weitz, D. A., and Wang, P. (2020). Rolling Waves with Non-Paraxial Phonon Spins. arXiv preprint arXiv:2009.13014

Google Scholar

Zhang, P., Wei, P., and Li, Y. (2017). Reflection of Longitudinal Displacement Wave at the Viscoelastically Supported Boundary of Micropolar Half-Space. Meccanica 52, 1641–1654. doi:10.1007/s11012-016-0514-z

CrossRef Full Text | Google Scholar

Zhu, H., Yi, J., Li, M.-Y., Xiao, J., Zhang, L., Yang, C.-W., et al. (2018). Observation of Chiral Phonons. Science 359, 579–582. doi:10.1126/science.aar2711

PubMed Abstract | CrossRef Full Text | Google Scholar

Zoeppritz, K. (1919). Erdbebenwellen viii b, uber reflexion and durchgang seismischer wellen durch unstetigkeisflachen. Gottinger Nachr 1, 66–84.

Google Scholar

Zuo, Q. H., and Hjelmstad, K. D. (1997). Bounds and Approximations for Elastic Wave Speeds in Cubic Crystals. J. Acoust. Soc. America 101, 3415–3420. doi:10.1121/1.418381

CrossRef Full Text | Google Scholar

Keywords: spin, elastic wave, rolling wave, cubic material, metamaterial, wave reflection, free boundary, fixed boundary

Citation: Zhang P and Wang P (2021) Boundary Reflections of Rolling Waves in Cubic Anisotropic Material. Front. Mech. Eng 7:704192. doi: 10.3389/fmech.2021.704192

Received: 01 May 2021; Accepted: 16 July 2021;
Published: 17 August 2021.

Edited by:

Chengzhi Shi, Georgia Institute of Technology, United States

Reviewed by:

Mostafa Habibi, Duy Tan University, Vietnam
Fan Ye, University of Massachusetts Amherst, United States
Chen Shen, Rowan University, United States

Copyright © 2021 Zhang and Wang. 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: Pai Wang, pai.wang@utah.edu

Download