skip to main content


Title: Temperature-controlled spatiotemporally modulated phononic crystal for achieving nonreciprocal acoustic wave propagation
We computationally investigate a method for spatiotemporally modulating a material's elastic properties, leveraging thermal dependence of elastic moduli, with the goal of inducing nonreciprocal propagation of acoustic waves. Acoustic wave propagation in an aluminum thin film subjected to spatiotemporal boundary heating from one side and constant cooling from the other side was simulated via the finite element method. Material property modulation patterns induced by the asymmetric boundary heating are found to be non-homogenous with depth. Despite these inhomogeneities, it will be shown that such thermoelasticity can still be used to achieve nonreciprocal acoustic wave propagation.  more » « less
Award ID(s):
1640860
NSF-PAR ID:
10353185
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
The Journal of the Acoustical Society of America
Volume:
151
Issue:
6
ISSN:
0001-4966
Page Range / eLocation ID:
3669 to 3675
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The existing concepts of non-reciprocity in propagation of acoustic or elastic waves are based either on nonlinear effects, or on local circulation of linear elastic fluid that leads to red or blue Doppler shift, depending on the direction of sound wave. The same concepts exist for electromagnetic non-reciprocity, where external magnetic field may produce the effect similar to local rotation of the medium. These two concepts originate from two known methods of breaking a time-reversal symmetry (T-symmetry), that is necessary for observation of nonreciprocal wave propagation. Both concepts require additional electrical or mechanical devices to be installed with their own power sources. Here we propose to explore viscosity of fluid as a natural factor of the T-symmetry breaking through energy dissipation. We report experimental observation of the nonreciprocal transmission of ultrasound through a water-submerged phononic crystals consisting of several layers of aluminum rods arranged in a square lattice. While viscous losses break the T-symmetry, making the wave propagation thermodynamically irreversible, the transmission remains reciprocal if the scatterers are symmetrical. To generate different energy losses for opposite directions of propagation, the P-symmetry of the crystal is broken by using asymmetric scatterers. Due to asymmetry, two sound waves propagating in the opposite directions produce different distributions of velocity and pressure that leads to different local absorption. Dissipation of acoustic energy occurs mostly near the surface of the scatterers and it strongly depends on surface roughnesses. Using two phononic crystal with smooth and rough aluminum rods we demonstrate low (2-5 dB) and high (10-15 dB) level of non-reciprocity within a wide range of frequencies, 300-600 kHz. Experimental results are in agreement with numerical simulations based on the Navier-Stokes equation. This nonreciprocal linear device is very cheap, robust and does not require energy source. 
    more » « less
  2. Abstract

    Magnetoelastic coupling is considered as one of the most reliable method to induce nonreciprocity of propagation losses of microwave‐frequency surface acoustic waves (SAW) and other acoustic modes propagating in nonmagnetic‐ferromagnetic heterostructures. Here, it is demonstrated theoretically that magnetoelastic coupling can also induce phase nonreciprocity of SAW, which is necessary for the development of SAW circulators and other nonreciprocal solid‐state‐acoustic devices. In contrast to previous studies, induction of the phase nonreciprocity requires the coupling of SAW to a strongly nonreciprocal spin wave (SW), having the nonreciprocal splitting of the SW spectrum much larger than the strength of the magnetoelastic coupling, which, in turn, should be much larger than the geometric mean of the SW and SAW damping rates. In this case, the hybridized SAW in the spectral region between the magnetoelastic gaps demonstrate significant phase nonreciprocity, retaining, at the same time, propagation losses that are close to those of unhybridized SAW. Possible practical realization of nonreciprocal SAW phase shifters and SAW‐ring‐based circulators based on hybridized waves in acoustic crystal and synthetic antiferromagnetic heterostructures is discussed.

     
    more » « less
  3. Abstract

    The ability to create linear systems that manifest broadband nonreciprocal wave propagation would provide for exquisite control over acoustic signals for electronic filtering in communication and noise control. Acoustic nonreciprocity has predominately been achieved by approaches that introduce nonlinear interaction, mean-flow biasing, smart skins, and spatio-temporal parametric modulation into the system. Each approach suffers from at least one of the following drawbacks: the introduction of modulation tones, narrow band filtering, and the interruption of mean flow in fluid acoustics. We now show that an acoustic media that is non-local and active provides a new means to break reciprocity in a linear fashion without these deleterious effects. We realize this media using a distributed network of interlaced subwavelength sensor–actuator pairs with unidirectional signal transport. We exploit this new design space to create a stable metamaterial with non-even dispersion relations and electronically tunable nonreciprocal behavior over a broad range of frequencies.

     
    more » « less
  4. SUMMARY

    Tsunami generation by offshore earthquakes is a problem of scientific interest and practical relevance, and one that requires numerical modelling for data interpretation and hazard assessment. Most numerical models utilize two-step methods with one-way coupling between separate earthquake and tsunami models, based on approximations that might limit the applicability and accuracy of the resulting solution. In particular, standard methods focus exclusively on tsunami wave modelling, neglecting larger amplitude ocean acoustic and seismic waves that are superimposed on tsunami waves in the source region. In this study, we compare four earthquake-tsunami modelling methods. We identify dimensionless parameters to quantitatively approximate dominant wave modes in the earthquake-tsunami source region, highlighting how the method assumptions affect the results and discuss which methods are appropriate for various applications such as interpretation of data from offshore instruments in the source region. Most methods couple a 3-D solid earth model, which provides the seismic wavefield or at least the static elastic displacements, with a 2-D depth-averaged shallow water tsunami model. Assuming the ocean is incompressible and tsunami propagation is negligible over the earthquake duration leads to the instantaneous source method, which equates the static earthquake seafloor uplift with the initial tsunami sea surface height. For longer duration earthquakes, it is appropriate to follow the time-dependent source method, which uses time-dependent earthquake seafloor velocity as a forcing term in the tsunami mass balance. Neither method captures ocean acoustic or seismic waves, motivating more advanced methods that capture the full wavefield. The superposition method of Saito et al. solves the 3-D elastic and acoustic equations to model the seismic wavefield and response of a compressible ocean without gravity. Then, changes in sea surface height from the zero-gravity solution are used as a forcing term in a separate tsunami simulation, typically run with a shallow water solver. A superposition of the earthquake and tsunami solutions provides an approximation to the complete wavefield. This method is algorithmically a two-step method. The complete wavefield is captured in the fully coupled method, which utilizes a coupled solid Earth and compressible ocean model with gravity. The fully coupled method, recently incorporated into the 3-D open-source code SeisSol, simultaneously solves earthquake rupture, seismic waves and ocean response (including gravity). We show that the superposition method emerges as an approximation to the fully coupled method subject to often well-justified assumptions. Furthermore, using the fully coupled method, we examine how the source spectrum and ocean depth influence the expression of oceanic Rayleigh waves. Understanding the range of validity of each method, as well as its computational expense, facilitates the selection of modelling methods for the accurate assessment of earthquake and tsunami hazards and the interpretation of data from offshore instruments.

     
    more » « less
  5. SUMMARY

    In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

     
    more » « less