🔬 Acoustic Metamaterials
Arrays of subwavelength Helmholtz resonators give a material an effective negative bulk modulus near resonance, opening acoustic band gaps and enabling negative refraction of sound.
How it Works
Each unit cell of this 1D lattice contains a small local resonator — a Helmholtz-resonator-like mass-spring element with natural (angular) frequency ω₀ = √(k/m) — weakly coupled to a background medium that would otherwise carry sound at a fixed speed. Far from ω₀ the resonators barely respond and the lattice behaves like an ordinary elastic medium. Near ω₀, however, the resonators absorb and re-radiate energy with a phase lag that can flip the sign of the material's effective, frequency-dependent bulk modulus K_eff(ω) — the Lorentz-oscillator dispersion familiar from optics, transplanted into acoustics.
Wherever K_eff(ω) is negative, the wave equation admits no propagating solution: the wavenumber becomes imaginary, and an incident wave decays exponentially with distance instead of oscillating through the lattice. This forbidden frequency range is the acoustic band gap. Because the resonators can be built far smaller than the sound wavelength they control, this locally-resonant band gap opens at frequencies that would normally need a lattice many times larger (as in ordinary Bragg-scattering phononic crystals) — the defining trick of acoustic metamaterials.
ω₀ = √(k/m) (unit-cell resonance)
Band gap: f₀ < f < √(f₀² + fₚ²)
k(ω) = ω · √(ρ₀ / K_eff(ω)) — complex wavenumber, wave ∝ e^(i(kx−ωt))
Frequently Asked Questions
What makes acoustic metamaterials different from ordinary materials?
Ordinary acoustic materials get their properties from chemistry — the atoms and bonds of the bulk medium. Acoustic metamaterials instead get their properties from engineered, subwavelength structure: periodic arrays of resonant unit cells (such as Helmholtz resonators) that respond to sound in ways no natural homogeneous material can, including effective negative bulk modulus or negative mass density near resonance.
Does a "negative" effective bulk modulus mean the material has negative volume or negative stiffness?
No. It means that over a narrow frequency band near the resonators' natural frequency, the material's dynamic, frequency-dependent response to an oscillating pressure is out of phase with the driving pressure — the effective modulus used in the wave equation flips sign in that band. The material itself is perfectly normal; it is the driven, frequency-dependent effective description that goes negative.
What is a band gap in an acoustic metamaterial?
A band gap is a range of frequencies for which no propagating wave solution exists in the periodic structure. Sound entering the material at a band-gap frequency cannot travel through it as an oscillating wave; instead its amplitude decays exponentially with distance (an evanescent wave), so the structure acts as a highly effective, often subwavelength, acoustic barrier at those frequencies.
What are Helmholtz resonators and why are they used?
A Helmholtz resonator is a cavity connected to the outside by a narrow neck; the air in the neck acts like a mass and the air in the cavity acts like a spring, giving it a single natural resonance frequency. Because a Helmholtz resonator can be made much smaller than the sound wavelength it responds to, arrays of them let you build a locally-resonant metamaterial whose unit cell is far smaller than the wavelength of the sound it controls.
What is negative refraction of sound and how does it arise?
Negative refraction occurs when the effective phase velocity of sound in a metamaterial points opposite to the direction energy actually flows (the group velocity), so an incident beam bends to the same side of the normal as the incoming beam instead of the opposite side. It arises near a resonance where the effective bulk modulus and/or effective density become negative simultaneously, producing an effective negative refractive index for sound.
What real-world applications use acoustic metamaterials?
Applications include thin locally-resonant soundproofing panels that block low frequencies far better than their mass alone would allow, acoustic lenses and superlenses for subwavelength imaging, acoustic cloaking devices that route sound waves around an object, and compact low-frequency noise barriers for machinery, HVAC ducts, and traffic noise.
How do acoustic metamaterials relate to electromagnetic metamaterials?
They are mathematical cousins. Electromagnetic metamaterials use split-ring resonators and wire arrays to create effective negative permittivity or permeability, described by Lorentz-oscillator dispersion. Acoustic metamaterials use Helmholtz resonators or local mass-spring resonators to create an effective negative bulk modulus or negative mass density, described by the same Lorentz-oscillator mathematics — only the physical fields differ.
Can a locally-resonant metamaterial block sound without adding much mass?
Yes, and this is one of their most useful properties. Ordinary mass-law sound insulation requires adding mass to block low frequencies. A locally-resonant metamaterial panel can present a very large (or negative) effective mass density near its resonators' natural frequency, blocking a target low-frequency band far more effectively than the mass-law prediction for the same physical mass — a genuine bypass of the classical mass law in a narrow band.
What limits real acoustic metamaterials in practice?
Damping and material losses narrow and shallow the band gap compared to the idealized lossless prediction. The strong effects are also inherently narrowband, centered on the resonator design frequency, and fabrication tolerances shift that frequency in real devices. Achieving broadband or multi-band performance generally requires stacking multiple resonator designs with different resonance frequencies.
About this simulation
This simulator drives a one-dimensional chain of N locally-resonant unit cells with a plane wave at frequency f and computes, cell by cell, whether that frequency propagates through the lattice or decays evanescently inside an acoustic band gap. The effective bulk modulus K_eff(ω) is evaluated directly from a Lorentz-oscillator model of the embedded Helmholtz-like resonators, and the same formula drives both the animated lattice and the transmission-vs-frequency curve beneath it.
🔬 What it shows
A live lattice of N cells whose instantaneous displacement is computed from the complex wavenumber k(f), so the driven wave either propagates smoothly across the chain or decays exponentially within a few cells — plus a transmission-vs-frequency curve, computed from the same K_eff(ω) model, that traces out the band gap analytically.
🎮 How to use
Set the unit-cell resonance f₀, coupling strength ωp (band-gap width), damping γ, driving frequency f, and cell count N, then watch the Status readout flip between "Propagating" and "Band-gap (evanescent)" as f crosses into [f₀, √(f₀²+fₚ²)]. P pauses, R resets.
💡 Did you know?
Locally-resonant sonic crystals can open a band gap at wavelengths dozens of times larger than the unit cell itself — unlike ordinary Bragg-gap phononic crystals, which need a cell size comparable to half the sound wavelength to block the same frequency.
Frequently asked questions
Why does the displayed band gap range depend only on f₀ and fₚ, not on damping?
bandGapEdges() returns [params.f0, Math.sqrt(f0*f0+fp*fp)], the zero-crossing of the undamped Lorentz form of keffOf(); damping γ only enters the imaginary part of the denominator inside keffOf(), which smooths and shallows the transition at the edges rather than moving the nominal edges themselves.
How does the Damping slider change the picture without changing the [f_low, f_high] readout?
Increasing params.gamma raises the imaginary part of the denominator inside keffOf(f), which raises Im(waveNumberOf(f)) even for frequencies outside the nominal gap — so higher damping adds real loss to the "propagating" region on the lattice and transmission curve while the analytic band-gap numbers shown in the stats stay fixed.
Why does raising the Number of unit cells N change the decay envelope but not the transmission formula itself?
transmissionAt(f) computes Math.exp(-waveNumberOf(f).im * params.N); increasing N only multiplies the same per-cell decay rate by a larger exponent, so a longer lattice attenuates a band-gap frequency far more, while keffOf() and waveNumberOf() — the physics per unit cell — never change with N.
Why is the imaginary part of k forced to be non-negative in waveNumberOf()?
waveNumberOf() takes the complex square root of ρ₀/K_eff(ω) and then applies Math.abs() to its imaginary part before scaling by ω, selecting the physically causal branch so that energy always decays away from the driven left edge of the lattice rather than spuriously growing — the correct behavior for a passive, lossy medium.
Why does increasing the Coupling strength ωp widen the band gap instead of shifting its center?
params.fp enters bandGapEdges() only inside the square root as f_high = √(f0² + fp²), so a larger fp pushes the upper edge outward while f_low stays pinned to params.f0 — exactly matching the physical picture that stronger coupling between the host medium and the local resonators widens the frequency window over which keffOf() stays negative.
Arrays of subwavelength Helmholtz resonators give a material an effective negative bulk modulus near resonance, opening acoustic band gaps and enabling negative refraction of sound.
2D · HTML5 Canvas 2D · 60 FPS target · runs fully client-side, no install