Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHASE VELOCITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Broadband surface-wave transformation cloak

Guiding surface electromagnetic waves around disorder without disturbing the wave amplitude or phase is in great demand for modern photonic and plasmonic devices, but is fundamentally difficult to realize because light momentum must be conserved in a scattering event. A partial realization has been achieved by exploiting topological electromagnetic surface states, but this approach is limited to narrow-band light transmission and subject to phase disturbances in the presence of disorder. Recent advances in transformation optics apply principles of general relativity to curve the space for light, allowing one to match the momentum and phase of light around any disorder as if that disorder were not there. This feature has been exploited in the development of invisibility cloaks. An ideal invisibility cloak, however, would require the phase velocity of light being guided around the cloaked object to exceed the vacuum speed of light—a feat potentially achievable only over an extremely narrow band. In this paper, we theoretically and experimentally show that the bottlenecks encountered in previous studies can be overcome. We introduce a class of cloaks capable of remarkable broadband surface electromagnetic waves guidance around ultrasharp corners and bumps with no perceptible changes in amplitude and phase. These cloaks consist of specifically designed nonmagnetic metamaterials and achieve nearly ideal transmission efficiency over a broadband frequency range from 0 + to 6 GHz. Finally, this work provides strong support for the application of transformation optics to plasmonic circuits and could pave the way toward high-performance, large-scale integrated photonic circuits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Development of long-pulse heating & current drive actuators & operational techniques compatible with a high-Z divertor & first wall

This was a collaboration between UCSD and MIT to study the effective application of ion-cyclotron heating (ICRH) on the EAST tokamak, located in China. The original goal was for UCSD to develop a diagnostic that would allow measurement of the steady state, or DC, convection pattern that develops on magnetic field lines that attach or connect to the ICRH antenna. This diagnostic would then be used to develop techniques and approaches that minimize or even eliminate such DC convection during application of strong ICRH heating. This was thought to then indicate reduction or elimination of parasitic losses of heating power, and thus be an indicator of effective RF heating. The original plan to use high speed digital gas-puff imaging (GPI) of the antenna-edge plasma region in EAST was ultimately unsuccessful due to limitations in machine and camera operations. We then decided to attempt the same experiment on the ALCATOR C-MOD tokamak at MIT which had a similar instrument already installed. This effort was ultimately successful, and demonstrated that the underlying idea of using GPI as a diagnostic for ICRH antenna physics would, in fact, work. The two-dimensional velocity fields of the turbulent structures, which are advected by RF-induced E x B flows, are obtained via the time-delay estimation (TDE) techniques. Both the magnitude and radial extension of the radial electric field E-r were observed to increase with the toroidal magnetic field strength B and the ICRF power. The TDE estimations of RF-induced plasma potentials are consistent with previous results based on the probe measurements of poloidal phase velocity. The results suggest that effective ICRH heating with reduced impurity production is possible when the antenna/box system is designed so as to reduce the RF-induced image currents that flow in the grounded conducting antenna frame elements that surround the RF antenna current straps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Mean flows and blob velocities in scrape-off layer (SOLT) simulations of an L-mode discharge on Alcator C-Mod

Two-dimensional scrape-off layer turbulence (SOLT) code simulations are compared with an L-mode discharge on the Alcator C-Mod tokamak [M. Greenwald, et al., Phys. Plasmas 21, 110501 (2014)]. Density and temperature profiles for the simulations were obtained by smoothly fitting Thomson scattering and mirror Langmuir probe (MLP) data from the shot. Simulations differing in turbulence intensity were obtained by varying a dissipation parameter. Mean flow profiles and density fluctuation amplitudes are consistent with those measured by MLP in the experiment and with a Fourier space diagnostic designed to measure poloidal phase velocity. Blob velocities in the simulations were determined from the correlation function for density fluctuations, as in the analysis of gas-puff-imaging (GPI) blobs in the experiment. In the simulations, it was found that larger blobs moved poloidally with the ExB flow velocity, v E , in the near-SOL, while smaller fluctuations moved with the group velocity of the dominant linear (interchange) mode, v E + 1/2 v di , where v di is the ion diamagnetic drift velocity. Comparisons are made with the measured GPI correlation velocity for the discharge. The saturation mechanisms operative in the simulation of the discharge are also discussed. In conclusion, it is found that neither sheared flow nor pressure gradient modification can be excluded as saturation mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Mean flows and blob velocities in scrape-off layer (SOLT) simulations of an L-mode discharge on Alcator C-Mod

Two-dimensional scrape-off layer turbulence (SOLT) code simulations are compared with an L-mode discharge on the Alcator C-Mod tokamak [M. Greenwald, et al., Phys. Plasmas 21, 110501 (2014)]. Density and temperature profiles for the simulations were obtained by smoothly fitting Thomson scattering and mirror Langmuir probe (MLP) data from the shot. Simulations differing in turbulence intensity were obtained by varying a dissipation parameter. Mean flow profiles and density fluctuation amplitudes are consistent with those measured by MLP in the experiment and with a Fourier space diagnostic designed to measure poloidal phase velocity. Blob velocities in the simulations were determined from the correlation function for density fluctuations, as in the analysis of gas-puff-imaging (GPI) blobs in the experiment. In the simulations, it was found that larger blobs moved poloidally with the ExB flow velocity, vE , in the near-SOL, while smaller fluctuations moved with the group velocity of the dominant linear (interchange) mode, vE + 1/2 vdi, where vdi is the ion diamagnetic drift velocity. Comparisons are made with the measured GPI correlation velocity for the discharge. The saturation mechanisms operative in the simulation of the discharge are also discussed. It is found that neither sheared flow nor pressure gradient modification can be excluded as saturation mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Study of the film boiling heat transfer and two-phase flow interface behavior using image processing

Here, in this study, based on a small-scale quench test facility, the two-phase flow interface behavior during quench transients is visualized and analyzed utilizing an image processing framework. The high-fidelity experimental results obtained for two-phase flow in the current framework can support various studies both in the time domain and in the frequency domain. In particular, visualization of the data obtained from different heating surfaces under different test conditions are used to perform a full-scale transient 2-D vapor film reconstruction. The liquid-vapor interface variations in various heat transfer regimes as well as at the initial film breakup point can be directly obtained through the processed data. Moreover, the temporal variation of the interfacial wave frequency approaching quench is investigated in detail. Based on the high-resolution data obtained for the liquid-vapor interface, the detailed phase velocity and temperature profiles are obtained through theoretic analysis, based on which the film boiling heat transfer coefficient (HTC) can be determined. In addition, an improved film boiling HTC model is developed considering the effects of wall superheat, liquid subcooling temperature, vapor film thickness as well as fluid properties. The model is found to predict film boiling HTC well within 15% error.

42 ENGINEERING↗
DOE OSTI2018

Characteristics of infrasound signals from North Korean underground nuclear explosions on 2016 January 6 and September 9

This study investigates infrasound signals from North Korean underground nuclear explosions on 2016 January 6 at 01:30:01 UTC (UNEJ16) and on 2016 September 9 at 00:30:01 UTC (UNES16). Infrasound observations from eight, seismo-acoustic arrays at distances of 300–600 km from the explosions and aligned approximately east–west across South Korea are used in the analysis. Progressive multichannel correlation was used followed by analyst review to identify the infrasound arrivals with estimates of azimuth, phase velocity, celerity, F-statistic, signal-to-noise ratio, correlation and signal duration. These observations are compared with model predictions of travel path, ray turning height and relative amplitudes, using Ground-to-Space atmospheric specifications at the time of each explosion. Stratospheric conditions at the time of the UNEJ16 are favourable to the Korean arrays when all arrays detected infrasound signals, while the few detections for UNES16 are consistent with unfavourable stratospheric winds and reduced amplitude thermospheric arrivals. Infrasound locations of the explosions were estimated using the Bayesian Infrasonic Source Location method (BISL), with wind-corrected backazimuth estimates based on the ray tracing calculations. Adding wind corrections to BISL resulted in infrasound locations up to 50 per cent closer to the seismic epicentre over uncorrected locations. The UNES16 infrasound location estimate using atmospheric corrections had larger differences (~67 km difference from seismic location) than the UNEJ16 location (~11 km difference from seismic location) as a result of the reduced number of observations and poor azimuth coverage for UNES16. Based on stratospheric phases recorded at six arrays and the empirical yield-scaling relations, wind corrected amplitudes from UNEJ16 lead to equivalent infrasound source energy estimates between 0.9 and 16.1 tons of TNT detonated in the atmosphere with an average of 6.4 tons and a standard deviation of 4.6 tons.

58 GEOSCIENCES↗
DOE OSTI2017

Electrically Tunable Integrated Thin-Film Magnetoelectric Resonators

Magnetoelectrics have attracted much attention for their ability to control magnetic behavior electrically and electrical behavior magnetically. This feature provides numerous benefits to electronic systems and can potentially serve as the bridge needed to integrate magnetic devices into mainstream electronics. This natural next step is pursued and thin-film integrated magnetoelectric devices are produced for radio-frequency (RF) electronics. The first fully integrated, thin-film magnetoelectric modulators for tunable RF electronics are presented. Moreover, these devices provide electric field control of magnetic permeability in order to change the phase velocity and resonance frequency of coplanar waveguides. During this study, the various thin-film material phenomena, trade-offs, and integration considerations for composite magnetoelectrics are analyzed and discussed. The fabricated devices achieve reversible tunability of the resonance frequency, characterized by a remarkable converse magnetoelectric coupling coefficient of up to 24 mG cm V -1 using just thin films. Based on this work, suggestions are given for additional optimizations of future designs that will maximize the thin-film magnetoelectric interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Dynamics of High Sound-Speed Metal Confiners Driven By Non-Ideal High-Explosive Detonation

Here, the results of 14 tests examining the behavior of aluminum (Al) conifners driven by non-ideal ANFO detonation in a cylinder test configuration are presented. In each test, the measured detonation phase velocity is slower than the aluminum sound speed. Thus, in the detonation reference frame, the ow in the Al is both shockless and subsonic. The tests involve: 3-inch inner diameter (ID) cylinders with Al wall thicknesses of 1/4, 3/8, 1/2, 1 and 2 inches; a 4-inch ID cylinder with a 1/2-inch Al wall thickness; and 6-inch ID cylinders with Al wall thicknesses of 1/2, 1 and 2 inches. The ANFO detonation velocity is seen to increase with increasing wall thickness for both the 3- and 6-inch ID tests, with no limiting velocity reached for the wall thicknesses used. The motion of the outer Al wall due to precursor elastic waves in the Al running ahead of the detonation is also measured at various axial locations along the cylinders. It is found that the magnitude of the outer wall motion due to the precursor elastic waves is small, while the associated wall motion is unsteady and decays in amplitude as the elastic disturbances move further ahead of the detonation front. The variations in the expansion history of the main outer wall motion of the cylinders are presented for increasing wall thickness at fixed ID, and for increasing cylinder inner diameter at a fixed wall thickness. Finally, we also explore the existence of a geometric similarity scaling of the wall expansion history for three geometrically scaled tests (3- and 6-inch ID cylinders with 1/4- and 1/2-inch walls respectively, 3- and 6-inch ID cylinders with 1/2- and 1-inch walls and 3- and 6-inch ID cylinders with 1- and 2-inch walls respectively). We find that the wall velocity histories for each of the three scaled tests, when plotted directly against time relative to start of main motion of the wall, are similar over a certain range of wall velocities without any geometric based rescaling in time. The range of wall velocities where the overlap occurs increases as the ratio of the wall thickness to inner diameter decreases. In conclusion, this is in contrast to ideal high explosives, where the outer wall velocity histories are only similar when the geometric scale factor (in this case a factor of 2) is applied to the wall velocity motion.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

A theoretical study of the upper bound of surface elevation variance in the Phillips initial stage during wind-wave generation

The resonance mechanism in the initial of wind-wave generation proposed by Phillips is a foundation of wind-wave generation theory, but a precise theoretical quantification of wave energy growth in this initial stage has not been obtained yet after more than six decades of research. In this study, we aim to address this knowledge gap by developing an analytical approach based on a novel complex analysis method to theoretically investigate the temporal evolution of the wave energy in the Phillips initial stage. We quantitatively derive and analyse the growth behaviour of the surface wave energy and obtain an analytical solution for its upper bound. Our result highlights the crucial effects of surface tension. Because the phase velocity of gravity–capillary waves has a minimal value at a critical wavenumber, gravity–capillary waves and gravity waves (which neglect surface tension) exhibit distinct resonance curve properties and wave energy growth behaviours. For gravity waves, the resonance curve extends indefinitely; for gravity–capillary waves, it either forms a finite-length curve or does not exist, depending on the wind speed. The leading-order term of the upper-bound solution of the energy of gravity waves increases linearly over time, while for gravity–capillary waves, the term increases linearly over time under strong wind conditions but remains finite under weak wind conditions. This theoretical study provides an analytical framework for the generation of wind-waves in the Phillips initial stage, which may inspire further theoretical, numerical and experimental research.

58 GEOSCIENCES↗
DOE OSTI2017

Generalized fluid theory including non-Maxwellian kinetic effects

The results obtained by the plasma physics community for the validation and the prediction of turbulence and transport in magnetized plasmas come mainly from the use of very central processing unit (CPU)-consuming particle-in-cell or (gyro)kinetic codes which naturally include non-Maxwellian kinetic effects. To date, fluid codes are not considered to be relevant for the description of these kinetic effects. Here, after revisiting the limitations of the current fluid theory developed in the 19th century, we generalize the fluid theory including kinetic effects such as non-Maxwellian super-thermal tails with as few fluid equations as possible. The collisionless and collisional fluid closures from the nonlinear Landau Fokker–Planck collision operator are shown for an arbitrary collisionality. Indeed, the first fluid models associated with two examples of collisionless fluid closures are obtained by assuming an analytic non-Maxwellian distribution function. One of the main differences with the literature is our analytic representation of the distribution function in the velocity phase space with as few hidden variables as possible thanks to the use of non-orthogonal basis sets. These new non-Maxwellian fluid equations could initiate the next generation of fluid codes including kinetic effects and can be expanded to other scientific disciplines such as astrophysics, condensed matter or hydrodynamics. As a validation test, we perform a numerical simulation based on a minimal reduced INMDF fluid model. The result of this test is the discovery of the origin of particle and heat diffusion. The diffusion is due to the competition between a growing INMDF on short time scales due to spatial gradients and the thermalization on longer time scales. Here, the results shown here could provide the insights to break some of the unsolved puzzles of turbulence.

42 ENGINEERING↗
DOE OSTI2021

Phase contrast imaging measurements and numerical simulations of turbulent density fluctuations in gas-fuelled ECRH discharges in Wendelstein 7-X

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the E × B rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the E × B velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Correlation between Spin and Orbital Dynamics during Laser-Induced Femtosecond Demagnetization

Spin and orbital angular momenta are two intrinsic properties of an electron and are responsible for the physics of a solid. How the spin and orbital evolve with respect to each other on several hundred femtoseconds is largely unknown, but it is at the center of laser-induced ultrafast demagnetization. In this paper, we introduce a concept of the spin–orbital correlation diagram, where spin angular momentum is plotted against orbital angular momentum, much like the position-velocity phase diagram in classical mechanics. We use four sets of highly accurate time-resolved X-ray magnetic circular dichroism data to construct four correlation diagrams for iron and cobalt. To our surprise, a pattern emerges. The trace on the correlation diagram for iron is an arc, and at the end of demagnetization, it has a pronounced cusp. The correlation diagram for cobalt is different and appears more linear but with kinks. We carry out first-principles calculations with two different methods: time-dependent density functional theory (TDDFT) and time-dependent Liouville density functional theory. These two methods agree that the experimental findings for both Fe and Co are not due to experimental errors. It is the spin–orbit coupling that correlates the spin dynamics to the orbital dynamics. Microscopically, Fe and Co have different orbital occupations, which leads to distinctive correlation diagrams. We believe that this correlation diagram presents a useful tool to better understand spin and orbital dynamics on an ultrafast time scale. A brief discussion on the magnetic anisotropy energy is also provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Tuning Kinetic Inductance with Doping in Superconducting Electron Gases at the KTaO 3 (111) Interface

We used coplanar waveguide resonators fabricated from the two-dimensional superconductor formed at KTaO 3 (111) interfaces to characterize the superconductor’s sheet kinetic inductance L K/□ . Upon varying the carrier density, in samples with T c ranging from 0.62 to 1.81 K, we can tune L K/□ between 1.88 and 7.42 nH/□ at the lowest temperatures, exceeding the highest values reported for granular aluminum films. The temperature dependence of L K/□ is consistent with a superconducting gap without nodes. The high L K/□ of superconducting KTaO 3 (111) interfacial electron gases combined with the high dielectric constant of KTaO 3 results in resonators with exceedingly low phase velocities (“slow light”) and small mode volumes. In addition, superconducting KTaO 3 (111) interfacial electron gases are robust in magnetic fields well above a Tesla. In conclusion, these features can enable unique device applications in superconducting electronics and quantum information science.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Self-Assembled Epitaxial Au–Oxide Vertically Aligned Nanocomposites for Nanoscale Metamaterials

Metamaterials made of nanoscale inclusions or artificial unit cells exhibit exotic optical properties that do not exist in natural materials. Promising applications, such as super-resolution imaging, cloaking, hyperbolic propagation, and ultrafast phase velocities have been demonstrated based on mostly micrometer-scale metamaterials and few nanoscale metamaterials. To date, most metamaterials are created using costly and tedious fabrication techniques with limited paths toward reliable large-scale fabrication. In this work, we demonstrate the one-step direct growth of self-assembled epitaxial metal–oxide nanocomposites as a drastically different approach to fabricating large-area nanostructured metamaterials. Using pulsed laser deposition, we fabricated nanocomposite films with vertically aligned gold (Au) nanopillars (~20 nm in diameter) embedded in various oxide matrices with high epitaxial quality. Strong, broad absorption features in the measured absorbance spectrum are clear signatures of plasmon resonances of Au nanopillars. By tuning their densities on selected substrates, anisotropic optical properties are demonstrated via angular dependent and polarization resolved reflectivity measurements and reproduced by full-wave simulations and effective medium theory. Our model predicts exotic properties, such as zero permittivity responses and topological transitions. In conclusion, our studies suggest that these self-assembled metal–oxide nanostructures provide an exciting new material platform to control and enhance optical response at nanometer scales.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling Frequency–Independent Q Viscoacoustic Wave Propagation in Heterogeneous Media

Quantifying the attenuation of seismic waves propagating in the Earth interior is critical to study the subsurface structure. Previous studies have proposed fractional anelastic wave equations to model the frequency–independent Q seismic wave propagation. Such wave equations involve fractional derivatives that pose computational challenges for the numerical schemes in terms of accuracy and efficiency when dealing with heterogeneous Earth media. To tackle these challenges, here we derive a new viscoacoustic wave equation, where the power terms of the fractional Laplacian operators are spatially independent, thus accurate and efficient methods (e.g., the Fourier pseudospectral method) can be adopted. Our derivation enables the resultant equation to capture both amplitude and phase signatures of the anelastic wave propagation by matching the complex wave numbers for all the frequencies of interest. We verify the derivation by comparing the dispersion curves of both the attenuation factor and the phase velocity produced by the new wave equation with their theoretical values as well as the Pierre Shale in situ measurements. Following that, we use a synthetic attenuating gas chimney model to demonstrate the attenuation effects on seismic waveforms and then construct the Q–compensated reverse time migration to undo these effects for seismic image enhancement. Finally, we find that our forward modeling results can characterize the spatiotemporal attenuation effects revealed in the Frio–II CO2 injection time–lapse seismic monitoring data. Here, we expect this proposed equation to be useful to quantify the attenuation in seismic data to push the resolution limits of seismic imaging and inversion.

25 ENERGY STORAGE↗
DOE OSTI2022

Ambient Noise Attenuation Tomography Reveals an Asymmetric Damage Zone Across San Jacinto Fault Near Anza, California

Here, we perform seismic attenuation tomography of the shallow structure for the San Jacinto Fault in the Ramona Reservation of southern California. The study uses ambient seismic noise recorded by a linear array of 65 3-C sensors across the fault. We extract amplitude decay information, with uncertainty, from noise interferometry functions. To account for strong heterogeneities in the complex shallow fault zone structure, we apply a frequency-dependent amplitude correction for focusing/defocusing effects using inverted phase velocity maps obtained by solving the transport equation. We then estimate the attenuation structure based on the linear station-triplet method for both Love and Rayleigh waves. The attenuation tomography indicates strong attenuation correlated with known San Jacinto Fault surface traces. The Love wave attenuation tomography reveals an asymmetric damage zone that exists primarily on the fault side with faster seismic velocity, consistent with earthquake ruptures on a fault bimaterial interface with preferred propagation direction to the northwest.

58 GEOSCIENCES↗
DOE OSTI2025

Structure of the Ecuadorian Upper Plate From a Joint Seismic‐Gravity Inversion

The Ecuadorian portion of the South American subduction zone presents an interesting case study in the structure and complex evolution of an upper plate. There are outstanding questions about its tectonic history, composition, and magmatic processes. While previous studies have employed ambient noise tomography to image the Ecuadorian upper plate, surface wave inversions alone often lack sensitivity at relevant shallow depths. This limitation can be overcome with an independent, complementary data set, such as gravity. We have jointly inverted Rayleigh wave phase velocities and Bouguer gravity anomalies to provide a more detailed seismic velocity model of the Ecuadorian upper plate. Our joint inversion has yielded several key improvements from previous models. First, we observe much shallower slow velocities beneath major basins (the Manabí, Progreso, and Gulf of Guayaquil), better aligning with expected basin structure. Second, we identify a high-velocity block beneath the entire forearc, corresponding to the Piñon Terrane, with velocities suggesting the presence of ultramafic material. Third, we highlight a new narrow swath of slow velocities beneath the Ecuadorian Andes, which closely follows the active volcanoes along the Eastern Cordillera. The extent of these slow velocities coincides with the termination of active arc volcanism and the predicted location of the subducted Carnegie Ridge. The predicted compositions for the mid to lower crust in the region preclude a purely compositional explanation for these velocities, suggesting that some level of partial melt is necessary.

Birkey, Andrew [Univ. of Delaware, Newark, DE (Uni↗
DOE OSTI2017

Observation of acoustic Dirac-like cone and double zero refractive index

Zero index materials where sound propagates without phase variation, holds a great potential for wavefront and dispersion engineering. Recently explored electromagnetic double zero index metamaterials consist of periodic scatterers whose refractive index is significantly larger than that of the surrounding medium. This requirement is fundamentally challenging for airborne acoustics because the sound speed (inversely proportional to the refractive index) in air is among the slowest. Here, we report the first experimental realization of an impedance matched acoustic double zero refractive index metamaterial induced by a Dirac-like cone at the Brillouin zone centre. This is achieved in a two-dimensional waveguide with periodically varying air channel that modulates the effective phase velocity of a high-order waveguide mode. Using such a zero-index medium, we demonstrated acoustic wave collimation emitted from a point source. For the first time, we experimentally confirm the existence of the Dirac-like cone at the Brillouin zone centre.

36 MATERIALS SCIENCE↗