Acoustic Metamaterials: Negative Bulk Modulus and Subwavelength Sound Attenuation

Acoustic Metamaterials: Negative Bulk Modulus and Subwavelength Sound Attenuation

Traditional acoustic noise barriers are bound by the classical Mass Law: mitigating low-frequency sound ($<200\text{ Hz}$, with long acoustic wavelengths $\lambda > 1.7\text{ m}$) requires prohibitively heavy, thick concrete or lead partitions ($>100\text{ kg/m}^2$). Sound attenuation scales as $6\text{ dB}$ per doubling of mass density.

Acoustic Metamaterials bypass the mass law through engineered subwavelength periodic unit cells, exhibiting effective negative mass density ($\rho_{\text{eff}} < 0$) or negative bulk modulus ($K_{\text{eff}} < 0$).

1. Mechanics of Negative Elastic Moduli

The constitutive wave equation in fluid acoustics relates pressure $p$ and particle velocity $\mathbf{v}$:

$$\nabla p = -\rho \frac{\partial \mathbf{v}}{\partial t}, \quad \nabla \cdot \mathbf{v} = -\frac{1}{K} \frac{\partial p}{\partial t}$$

Where sound velocity is $c = \sqrt{K / \rho}$.

• Locally Resonant Membranes (Negative Effective Density $\rho_{\text{eff}} < 0$): When driven above their structural resonance frequency, an array of elastic membranes with central attached masses oscillates out-of-phase with the incident acoustic pressure field. The system accelerates opposite to the acoustic pressure gradient, reflecting sound waves completely.
• Helmholtz Resonator Arrays (Negative Effective Modulus $K_{\text{eff}} < 0$): Cavities connected to the air channel via narrow necks act as acoustic inductors and capacitors. Near resonance, volume flow into the resonator is in anti-phase with acoustic pressure, creating an opaque subwavelength acoustic bandgap.

       Membrane Acoustic Metamaterial Unit Cell
       +---------------------------------------------+
       |             Rigid Frame Housing             |
       |     +---------------------------------+     |
       |     |        Elastic Membrane         |     |
       |     |        (Silicone / Latex)       |     |
       |     |             +-------+           |     |
       |     |             | Mass  |           |     |
       |     |             +-------+           |     |
       |     |                                 |     |
       |     +---------------------------------+     |
       |                                             |
       +---------------------------------------------+

2. Transmission Loss Formulations

The acoustic power transmission coefficient $T(\omega)$ across a unit cell of thickness $d$ with negative parameters demonstrates exponential evanescence:

$$k = \frac{\omega}{c_{\text{eff}}} = \omega \sqrt{\frac{\rho_{\text{eff}}}{K_{\text{eff}}}}$$

When only one parameter is negative, $k$ becomes purely imaginary ($k = j \alpha$), causing the acoustic wave amplitude to decay exponentially: $p(x) = p_0 e^{-\alpha x}$ without acoustic energy transmission.

import numpy as np

def compute_helmholtz_effective_modulus(freq_range, f_res=120.0, F_geom=0.25):
    '''
    Computes effective bulk modulus dispersion across Helmholtz cavity arrays.
    K_eff becomes negative in narrow band immediately above resonance.
    '''
    omega = 2 * np.pi * freq_range
    omega_res = 2 * np.pi * f_res
    gamma = 0.05 * omega_res # Damping factor

    # Dynamic response function
    k_ratio = 1.0 - (F_geom * (omega**2)) / (omega**2 - omega_res**2 - 1j * gamma * omega)
    return k_ratio

Acoustic metamaterials enable lightweight, sub-centimeter acoustic baffles for aircraft passenger cabins, MRI machine enclosures, and whisper-quiet ventilation ducts.