Non-Newtonian Fluids on Vibrating Diaphragms: Oobleck Wave

Executive Summary & Theoretical Thesis: Non-Equilibrium Hydrodynamic Transitions

Non-Newtonian Colloidal Suspension Under Vertical Acceleration

The hydrodynamic response of a dense particulate suspension subjected to high-amplitude, periodic vertical translation exposes a fundamental failure of classical continuum mechanics. When an aqueous suspension of unmodified cornstarch (Zea mays endosperm)—colloquially termed oobleck—is bounded within a rigid receptacle and mechanically coupled to an electrodynamic shaker, the system demonstrates an inversion of the classical Navier-Stokes boundary condition. Below critical excitation thresholds, the suspension behaves as a weakly shear-thinning or generalized Newtonian fluid with an effective kinematic viscosity governed by hydrodynamic lubrication layers separating adjacent starch granules.

When the non-dimensional peak acceleration $\Gamma$ surpasses the threshold of terrestrial gravity, the hydrodynamic shear field generated by vertical oscillatory displacement forces the interstitial liquid phase out from the interparticle gaps. The resulting transition converts lubricated hydrodynamic interactions into direct frictional contact networks. This phenomenon, known as discontinuous shear-thickening (DST), causes the local effective viscosity to diverge by orders of magnitude over millisecond timescales.

Subjecting a non newtonian fluid cornstarch oobleck vibrating speaker cymatics apparatus to vertical acceleration vectors induces dynamic, non-linear stress fields. These fields rupture the smooth, harmonic wave envelopes typically observed in Newtonian capillary-gravity systems.

[ Periodic Acceleration a(t) = a_0 cos(ωt) ]
                                      |
                                      v
          [ Interparticle Lubrication Layer Breakdown: h_gap -> 0 ]
                                      |
                                      v
               [ Frictional Contact Percolation (Jammed Core) ]
                                      |
                 +--------------------+--------------------+
                 |                                         |
                 v                                         v
   [ Persistent Standing Fluid Holes ]     [ Shear-Thickening Fingering Protrusions ]
   (Negative Curvature Topologies)         (Localized Hydrodynamic Solitons)

The Anomaly of Stable Dissipative Solitons

Under specific frequency-acceleration regimes ($f \in [30, 120]\text{ Hz}$, $\Gamma > 7.5$), the system generates localized macroscopic structures that do not disperse through gravitational leveling. These formations are recognized as acoustic hydrodynamic solitons and persistent standing fluid holes.

In standard Newtonian fluids, any localized void formed on an open liquid surface collapses immediately under hydrostatic head pressure $\rho g h$. The gravitational potential energy drives lateral fluid inflow, closing the void on a timescale set by the Rayleigh-Taylor or capillary time:

$$\tau_c = \sqrt{\frac{\rho R^3}{\gamma}}$$

In a dense colloidal suspension driven past the DST threshold, localized vertical kinetic energy converts into lateral compressive strain at the void boundary.

During the upward phase of the acceleration cycle, the fluid rim experiences high strain rates that trigger local dynamic jamming. The boundary of the void hardens into an anisotropic, solid-like wall capable of supporting shear and normal stresses. During the downward trajectory, the localized inertial unloading is insufficient to allow complete stress relaxation before the subsequent upward impact restiffens the perimeter.

The void remains open indefinitely as a non-equilibrium stationary defect: a stable dissipative soliton energized by the acoustic field.

Upward Acceleration Cycle                   Downward Descent Cycle
       =========================                   ======================
        a(t) > +g (Inertial Impact)                 a(t) < -g (Free Fall / Unloading)
                  ▲                                            ▼
                  │                                            │
   Solidified Rim  Void  Solidified Rim         Viscous Rim     Void    Viscous Rim
   [ JAMMED CORE ] (  ) [ JAMMED CORE ]        [ RELAXING ]    (    )  [ RELAXING ]
         ▲                 ▲                         │                    │
         └── σ_shear > σ_c ┘                         └── τ_relax > T_osc ─┘
     (Void Resists Hydrostatic Inflow)              (Incomplete Liquid Refill)

Paradigm Shift: From Linear Cymatics to Discontinuous Phase Boundaries

Linear cymatics, historically grounded in the observation of modal resonance Chladni patterns, maps the spatial distribution of standing wave nodes and antinodes across an elastic or fluid membrane. In those linear systems, the medium responds passively to the driving harmonic spectrum, and the surface topography maps directly to eigenfunctions of the Helmholtz operator:

$$\nabla^2 \psi + k^2 \psi = 0$$

Introducing a discontinuous shear-thickening suspension replaces this linear spatial distribution with a state-dependent phase boundary. At any coordinate $(x, y)$ along the horizontal aperture, the local phase state—viscous liquid versus jammed solid—is determined dynamically by the local strain tensor $\dot{\gamma}_{ij}(x, y, t)$.

The system transitions from an energy-dissipating wave medium to an active, excitable non-equilibrium system. Here, the boundaries between wave crests and troughs evolve into structural phase boundaries.

The classic Faraday wave instabilities observed on fluid interfaces are altered: rather than generating continuous capillary ripple fields, the oobleck interface undergoes dynamic spatial bifurcation. It self-organizes into persistent standing fluid holes, delaminating structural tongues, and high-aspect-ratio vertical protrusions. These non-equilibrium states demonstrate how macroscale mechanical stability can emerge from the interplay between high-frequency acoustic fields and granular lubrication transitions.

Historical Lineage & Experimental Precedents: From Faraday Crisping to Granular Bifurcations

Faraday’s 1831 Discovery of Parametric Surface Oscillations

The study of vertically oscillated fluid layers was inaugurated by Michael Faraday in his 1831 treatise presented to the Royal Society. Faraday investigated the dynamic crispations produced on liquid surfaces supported by vibrating plates, discovering that the free surface of a fluid undergoing vertical shaking destabilizes into standing waves whose primary oscillation frequency is precisely half the driving frequency:

$$f_{\text{response}} = \frac{1}{2} f_{\text{drive}}$$

This fundamental subharmonic response, now designated the Faraday instability, was derived under the assumption of an incompressible, linearly viscous Newtonian fluid. Faraday observed that light powders and fluid layers organize into specific geometric fields, yet his formulations treated the fluid’s shear modulus and viscosity as static material constants.

For more than a century and a half, fluid mechanics focused on Newtonian or weakly non-Newtonian solutions. Anomalous experimental behaviors were often dismissed as transient wall-boundary effects or contaminations of the air-liquid interfacial tension.

Twentieth-Century Rheology and Colloidal Science Precursors

During the late twentieth century, the convergence of granular physics, non-linear dynamics, and colloidal rheology revealed a wider spectrum of material behaviors. Investigators including J. F. Brady and N. J. Wagner demonstrated that dense colloidal suspensions diverge from classical Newtonian and power-law shear-thinning models.

When colloidal volume fractions approach the random close packing limit:

$$\phi_{\text{rcp}} \approx 0.64$$

the system exhibits hydrodynamically clustered states under shear.

As the applied stress surpasses a critical repulsive force stabilizing the particulate surfaces, the protective hydration or steric layers collapse. The suspension leaves the Stokesian hydrodynamic lubrication regime and enters a frictional contact regime, producing discontinuous shear-thickening (DST).

These rheological insights explained the steady-state properties of oobleck under rotational viscometry. However, they were rarely extended to calculate interfacial boundary dynamics under high-amplitude, parametric acoustic acceleration.

       Hydrodynamic Lubrication Regime             Frictional Contact Regime (DST)
      (Applied Shear Stress σ < σ_critical)       (Applied Shear Stress σ > σ_critical)
      ─────────────────────────────────────       ─────────────────────────────────────
             Starch Granule   Starch Granule               Starch Granule   Starch Granule
               (○) ~~~~~~ (○)                               (●)======(●)
                     ▲                                            ▲
                     │                                            │
        Fluid film separates surfaces.             Lubrication film broken; direct
        Low resistance; shear-thinning.            frictional contact. Viscosity diverges.

The 2004 Texas Experiments: Discovery of Persistent Holes

The intersection of non-Newtonian rheology and parametric wave mechanics was clarified experimentally in 2004 at the University of Texas at Austin. Investigators F. S. Merkt, R. D. Deegan, D. I. Levinson, and Harry L. Swinney published observations of an anomalous hydrodynamic state: vertically shaken, dense suspensions of cornstarch and water yielded autonomous, localized circular voids that defied gravitational collapse.

By driving a 5 mm layer of dense oobleck on an electromechanical piston at driving accelerations exceeding:

$$\Gamma \approx 10$$

within a frequency band of 30 to 120 Hz, the Texas group observed that a mechanical perturbation (such as an air jet or physical probe) nucleated an isolated, persistent standing fluid hole.

Rather than closing due to hydrostatic pressure, these holes maintained an average diameter stable over millions of oscillation cycles. The perimeter of the hole rose above the mean fluid depth, forming an elevated rim that underwent subharmonic oscillations while maintaining an open interior that exposed the dry substrate.

Subsequent investigations by Ebata et al. (2009) and Falcon et al. (2007) confirmed that these voids represent localized, non-equilibrium dissipative solitons. They are sustained dynamically by a spatial balance between parametric energy injection, localized dynamic jamming, and stress relaxation.

Mathematical Formalism & Physical Mechanics: Navier-Stokes Discontinuity and Constitutive Rheology

Modified Navier-Stokes Formulations for Discontinuous Shear-Thickening (DST)

The hydrodynamic evolution of a viscous, incompressible continuum under an external acceleration field $\mathbf{g}(t)$ is governed by the generalized momentum conservation equation:

$$\rho \left( \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u} \right) = -\nabla p + \nabla \cdot \boldsymbol{\tau} + \rho \mathbf{g}_{\text{eff}}(t)$$

where $\rho$ denotes the suspension mass density, $\mathbf{u}$ is the velocity vector field, $p$ is the isotropic hydrostatic pressure, and $\boldsymbol{\tau}$ is the deviatoric stress tensor.

In a classical Newtonian system, the deviatoric stress tensor depends linearly on the symmetric strain rate tensor $\mathbf{D} = \frac{1}{2}\left(\nabla \mathbf{u} + (\nabla \mathbf{u})^T\right)$ via a constant dynamic viscosity:

$$\boldsymbol{\tau} = 2 \eta_0 \mathbf{D}$$

For dense colloidal suspensions undergoing discontinuous shear-thickening, this constitutive closure fails. The viscosity must be redefined as a dynamic, non-linear function of the second invariant of the rate-of-strain tensor:

$$\dot{\gamma} = \sqrt{2 \mathbf{D} : \mathbf{D}}$$

as well as the solid volume fraction $\phi$. The Wyart-Cates constitutive model formalizes this transition by mapping the suspension between two limiting jamming volume fractions: an unjammed, lubricated close-packing fraction $\phi_{\text{lub}} \approx 0.64$ and a lower, frictionally jammed packing fraction $\phi_{\text{fric}} \approx 0.55$.

The fraction of particle contacts that transition from lubricated to frictional regimes is governed by the local shear stress $\sigma$:

$$f(\sigma) = \exp\left(-\frac{\sigma^*}{\sigma}\right)$$

where $\sigma^*$ is the characteristic microscopic force threshold required to displace the stabilizing hydration layer.

The effective volume fraction is expressed through an interpolated jamming limit $\phi_J(f) = (1 - f)\phi_{\text{lub}} + f \phi_{\text{fric}}$, yielding the modified dynamic viscosity:

$$\eta(\dot{\gamma}, \phi) = \eta_s \left( 1 - \frac{\phi}{\phi_J(f(\sigma))} \right)^{-2}$$

As $\sigma \gg \sigma^*$ under rapid acoustic strain, $\phi_J(f) \to \phi_{\text{fric}}$. If the suspension is prepared at $\phi \approx \phi_{\text{fric}}$, the effective viscosity diverges discontinuously ($\eta \to \infty$).

The fluid undergoes an instantaneous transformation into an acoustic-inertially jammed state. The corresponding momentum balance must accommodate this rheological phase transition across moving internal boundaries.

[ Applied Acoustic Strain Rate \dot{\gamma} ]
                                                        |
                                                        v
                                    [ Shear Stress σ = η(\dot{\gamma}) \dot{\gamma} ]
                                                        |
                               +------------------------+------------------------+
                               |                                                 |
                               v                                                 v
                   Low Stress (σ < σ*)                              High Stress (σ > σ*)
                   f(σ) -> 0; φ_J -> φ_lub                          f(σ) -> 1; φ_J -> φ_fric
             =================================               =====================================
             Viscosity bounded; fluid flows freely.          Denominator (1 - φ/φ_J) -> 0; η -> ∞.
             System follows classical Navier-Stokes.         Discontinuous dynamic jamming occurs.

Parametric Mathieu Instability and Subharmonic Bifurcation

The fluid layer rests on a rigid horizontal diaphragm vibrating vertically with displacement:

$$z(t) = a \cos(\omega t)$$

In the non-inertial reference frame of the container, this translation maps to a time-modulated effective gravitational acceleration:

$$\mathbf{g}_{\text{eff}}(t) = -\left( g - a \omega^2 \cos(\omega t) \right) \hat{\mathbf{k}}$$

Linearizing the free surface perturbation $\zeta(\mathbf{x}\perp, t)$ across horizontal coordinates $\mathbf{x}\perp = (x, y)$ in terms of spatial Fourier modes $\zeta_k(t) e^{i \mathbf{k} \cdot \mathbf{x}_\perp}$ yields a generalized Mathieu equation for each wavenumber $k = |\mathbf{k}|$:

$$\frac{d^2 \zeta_k}{d t^2} + 2 \gamma_k \frac{d \zeta_k}{d t} + \left[ \Omega_k^2 - k \left( a \omega^2 \cos(\omega t) \right) \tanh(k h) \right] \zeta_k = 0$$

where $h$ is the fluid depth, $\gamma_k$ is the phenomenological wave damping rate, and $\Omega_k$ is the unforced dispersion relation for capillary-gravity waves:

$$\Omega_k^2 = \left( g k + \frac{\gamma_{\text{st}}}{\rho} k^3 \right) \tanh(k h)$$

with $\gamma_{\text{st}}$ representing interfacial surface tension.

In a Newtonian fluid, parametric resonance occurs in instability tongues within the $(a, \omega)$ parameter plane, with the strongest instability centered at the subharmonic resonance $\Omega_k \approx \frac{\omega}{2}$.

In oobleck suspensions, the damping coefficient $\gamma_k$ is not a static attenuation constant; it is an active variable driven by the local strain rate $\dot{\gamma} \sim k \dot{\zeta}_k$. When vertical acceleration drives the interface upward, the local strain rate forces the material past $\sigma^*$, shifting $\gamma_k$ to large values within fractions of a cycle.

This non-linear damping limits the growth of delocalized Fourier modes, preventing spatial cascade into classical hexagonal or square Faraday wave patterns. Instead, it confines kinetic energy to localized topological bifurcations.

Hydrodynamic Soliton Mechanics: Acoustic Radiation Pressure vs. Gravitational Collapse

The structural equilibrium of persistent standing fluid holes requires balance between two opposing force vectors integrated across the fluid cycle period $T = \frac{2\pi}{\omega}$.

The first is the inward hydrostatic driving pressure, which tends to collapse the hole via gravitational potential:

$$P_{\text{hydro}}® = \rho g (h_0 - z®)$$

The second is the outward acoustic-inertial radiation stress:

$$S_{\text{ac}} = \langle \rho u_r^2 - \tau_{rr} \rangle_T$$

where $u_r$ is the radial velocity of the suspension and $\tau_{rr}$ is the radial component of the deviatoric stress tensor.

During the downward phase of diaphragm motion ($a(t) < -g$), the fluid is in effective free fall. The normal stress drops, causing the rim to dilate laterally.

During the upward stroke, the base plate accelerates into the fluid layer with $\ddot{z}(t) \gg g$. The fluid responds with high vertical shear against the rigid floor, triggering the DST transition:

$$\sigma_{zr} = \eta(\dot{\gamma}) \frac{\partial u_z}{\partial r} > \sigma^*$$

The rim of the hole solidifies dynamically, converting vertical momentum into outward radial deflection. This dynamic dilatant wall acts as an inertial barrier, resisting the hydrostatic pressure of the surrounding fluid pool.

The hole remains stable because the cycle-averaged radial momentum flux balances the hydrostatic head:

$$\oint_{0}^{T} \left( P_{\text{hydro}}(R_{\text{hole}}) - \left[ \rho u_r^2 - \tau_{rr}(\dot{\gamma}) \right] \right) dt = 0$$

If this condition holds at a finite radius $R_{\text{hole}}$, the void avoids both capillary closure and unbounded expansion, stabilizing as an acoustic hydrodynamic soliton.

Empirical Evidence & Observational Data: Laboratory Mapping of Fingering Instabilities and Stable Voids

Parametric Phase Diagrams Across Acceleration-Frequency Space

Systematic laboratory mapping reveals distinct morphological regimes across the acceleration-frequency domain ($\Gamma$ vs. $f$). The following phase diagram summarizes empirical stability thresholds recorded across the parameter space:

    Acceleration (Γ)
         ▲
      16 ┼─────────────────────────────────────────────────────────────
         │  REGIME IV: High-Aspect-Ratio Fingering & Droplet Ejection
      12 ┼─────────────────────────────────┬───────────────────────────
         │                                 │ REGIME III: Persistent
       8 ┼─── REGIME II: Subharmonic       │ Standing Fluid Holes
         │    Faraday Instabilities        │ (Bistable Dynamic Solitons)
       4 ┼─────────────────────────────────┴───────────────────────────
         │  REGIME I: Quiescent Viscous Planar Interface
       0 ┴───────────┬─────────────────────────────┬───────────────────►
         0          30                            120             f (Hz)

At drive accelerations $\Gamma < 1.0$, the fluid layer remains planar and quiescent. Surface shear remains well below the critical activation threshold $\sigma^*$, and the system behaves as a standard, viscous liquid absorbing the acoustic energy through viscous dissipation.

Between $\Gamma \approx 1.0$ and $\Gamma \approx 4.0$, the system undergoes a classical subharmonic bifurcation, forming standing Faraday ripples with wavelengths governed by the effective gravity-capillary dispersion relation:

$$\lambda_F = \frac{2\pi}{k_F}$$

As acceleration increases past $\Gamma_c \approx 7.5$ within the 30–120 Hz excitation band, the uniform Faraday pattern destabilizes. The wave crests experience shear rates high enough to trigger dynamic jamming at their peaks.

This transition breaks spatial homogeneity, allowing isolated topological defects to emerge. If an initial void is introduced within this regime, it does not collapse. Instead, it enters a bistable equilibrium, forming a stable, persistent standing fluid hole.

At extreme accelerations ($\Gamma > 12.0$), the system enters a fingering and ejection regime. Here, shear-thickened structures delaminate from the bulk fluid and form high-aspect-ratio vertical protrusions.

High-Speed Optical Profilometry of Persistent Holes

High-speed profilometry (captured at frame rates $\ge 2000\text{ fps}$) paired with laser-sheet imaging clarifies the interior profile of persistent standing fluid holes.

  Laser Profilometry Cross-Section of Persistent Hole
  ==================================================
           Elevated Oscillating Rim (h_rim > h_0)
                 ┌──┐                     ┌──┐
                 │  │                     │  │        Mean Layer Depth (h_0)
  ───────────────┘  │                     │  └─────────────── 
  Bulk Fluid        │                     │        Bulk Fluid
  (Viscous Phase)   │  Exposed Dry Floor  │   (Viscous Phase)
  ══════════════════╧═════════════════════╧══════════════════ (Substrate Base)
                    │◄───  D_hole  ──────►│

Profilometry measurements demonstrate that the interior of the hole remains dry, exposing the rigid floor of the container. The fluid rim bounding the void does not rest at the quiescent layer depth $h_0$; it forms an elevated wall ($h_{\text{rim}} > h_0$) that undergoes an oscillation phase-locked to the subharmonic frequency $\frac{\omega}{2}$.

The cross-sectional dynamic profile reveals that the rim experiences maximum shear during the upward stroke of the transducer. As the base plate moves upward:

$$\ddot{z}(t) > 0$$

the fluid pushed into the rim undergoes high vertical shear rates:

$$\dot{\gamma}_{zr} \approx \frac{\Delta u_z}{\Delta r} \sim \mathcal{O}(10^3)\text{ s}^{-1}$$

This strain rate drives the particulate suspension beyond its critical packing fraction, jamming the rim into a semi-rigid ring.

During the downward stroke, the effective gravitational unloading allows partial surface relaxation. However, the suspension’s relaxation time constant $\tau_{\text{relax}}$ is longer than the half-period of oscillation:

$$\tau_{\text{relax}} > \frac{\pi}{\omega}$$

The rim retains residual mechanical integrity until the next cycle’s compression restiffens the structure.

Shear Thickening Fingering Phenomena: Dynamic Protrusions and Delamination

At high accelerations ($\Gamma > 10.0$), the free surface exhibits shear thickening fingering phenomena. These structures are not chaotic droplet splashes. Splash droplets in Newtonian fluids detach when fluid velocity exceeds the capillary escape velocity:

$$v_{\text{splash}} > \sqrt{\frac{2\pi \gamma_{\text{st}}}{\rho \lambda}}$$

forming spherical droplets driven by surface tension minimization.

In contrast, oobleck fingers emerge as coherent, high-aspect-ratio vertical columns that oscillate, bend, and sustain shear stresses without detaching.

       Newtonian Splash Protrusion               Non-Newtonian Oobleck Finger
      (Inertia-Capillary Instability)             (Dynamic Jamming Instability)
      ───────────────────────────────             ─────────────────────────────
                    ( ) Droplet Delamination                  ▲   Rigidified Tip
                     │                                        │   (Jammed Granular State)
                     ▼ Capillary Necking                     [ ]
                    ( )                                      [ ]  Coherent Column
                     │                                       [ ]  (Maintains Vertical Structure)
                  ┌─────┐                                    [ ]
                  │     │ Free Surface                     ┌─┴─┴─┐
             ─────┘     └─────                        ─────┘     └─────

These fingers are generated by local normal stress differences. Under high cyclic acceleration, starch granules experience kinematic confinement. Lateral strain fluctuations create local dilatant zones that solidify faster than surrounding fluid domains.

The fluid beneath these jammed zones pushes them upward, generating vertical columns. High-speed recordings reveal that these columns act like flexible, solid rods during upward acceleration, resisting lateral flow and deformation.

Only when the driving acceleration is interrupted does the dynamic yield stress drop back to zero:

$$\sigma_y \to 0$$

allowing the fingers to melt back into the bulk liquid pool over several hundred milliseconds.

Metaphysical Implications & Unified Synthesis: Self-Organized Morphogenesis and Universal Cymatic Archetypes

Dissipative Solitons as Archetypes of Discrete Matter Condensation

The emergence of persistent standing fluid holes and shear-thickened projections in an oscillating fluid layer demonstrates how continuous energy inputs can produce stable, localized, particle-like forms.

In classical field theories, matter particles are often modeled as localized, non-linear excitations—solitons—derived from continuous underlying fields. The non-Newtonian fluid under vertical vibration provides a macroscopic fluid-mechanical analog for this transition.

The system receives spatially uniform, non-localized mechanical energy from the vibrating diaphragm:

$$\mathbf{E}{\text{in}} = \oint \mathbf{F}{\text{speaker}}(t) \cdot d\mathbf{z}(t)$$

Without requiring localized external boundary constraints, the medium’s non-linear constitutive response generates stable, individual structures.

The persistent hole operates as a discrete macroscopic entity: it possesses a defined geometric boundary, retains its identity through mechanical disturbances, and can interact with or repel neighboring holes.

This transformation shows how homogeneous continuous energy fields can spontaneously condense into discrete material structures when mediated by a threshold-dependent constitutive law.

+-----------------------------------------------------------------------------------+
|               CONTINUOUS FIELD INPUT: a(t) = a_0 cos(ωt)                          |
|         (Uniform, Delocalized Acoustic Energy Across Entire Container Floor)      |
+-----------------------------------------------------------------------------------+
                                         │
                                         ▼
+-----------------------------------------------------------------------------------+
|                  NON-LINEAR THRESHOLD: σ_shear > σ_critical                       |
|               (Local Viscous Dissipation -> Frictional Jamming)                   |
+-----------------------------------------------------------------------------------+
                                         │
                                         ▼
+-----------------------------------------------------------------------------------+
|            DISCRETE TOPOLOGICAL CONDENSATION: DISSIPATIVE SOLITONS                 |
|             (Localized Persistent Holes, Self-Reinforcing Protrusions)            |
+-----------------------------------------------------------------------------------+

Continuous Field to Discontinuous Form: Non-Linear Morphogenesis

This hydrodynamic phase behavior offers insights into morphogenesis across physical scales. In biological development, cellular matrices transition between fluid-like behavior (yielding to tissue migration) and solid-like behavior (maintaining structural integrity) through jamming transitions controlled by cellular packing density and cortical tension.

The non-Newtonian cymatic system illustrates this morphogenetic transition in an accessible physical medium.

The system does not follow linear superposition. While linear acoustic fields generate distributed cymatic patterns defined by continuous sinusoidal functions, the non-linear colloidal interface generates discontinuous geometries characterized by sharp, moving phase boundaries:

$$\mathcal{H}(\mathbf{x}) \in {\text{Fluid, Jammed}}$$

The transition from a continuous oscillating field to a defined geometric form is mediated by the material’s internal dissipation thresholds.

The acoustic diaphragm does not imprint a specific shape into the fluid like a mold. Instead, it provides a dynamic scalar potential field.

The fluid’s own internal constraints (interparticle friction, steric repulsion, and dilatancy) translate this continuous field into discrete, stable physical structures. These principles of acoustic field manipulation also appear in systems utilizing acoustic levitation mechanics.

Resonant Geometry as the Fundamental Organizing Operator of Matter

The dynamic behaviors observed on the vibrating diaphragm illustrate how geometric form can emerge from the interplay between acoustic energy and non-linear material responses.

Physical form is revealed not as a static arrangement of inert matter, but as a dynamic balance between continuous energy fluxes and internal material dissipation thresholds.

The persistent hole, the oscillating rim, and the shear-thickened fingering projection are all stabilized by continuous energy throughput:

$$\frac{d E_{\text{internal}}}{dt} = \dot{E}{\text{acoustic}} - \dot{E}{\text{viscous}} = 0$$

These stable configurations mirror broader patterns across physics, from subatomic wave-particle dualities to dynamic plasma equilibria and dielectric polarization dynamics.

When driven beyond linear response regimes, continuous mechanical oscillations produce discrete, localized macroscopic structures. The oobleck wave demonstrates that material structure can be generated dynamically: geometric form emerges as the natural physical resolution when cyclic stress fields interact with non-linear material thresholds.

Frequently Asked Questions: Nonlinear Acoustics and Colloidal Rheology

Hydrodynamic Stability of Negative Curvature Voids

How does a persistent standing fluid hole prevent gravitational hydrostatic collapse without an internal physical support?

A persistent standing fluid hole resists hydrostatic collapse through cycle-averaged acoustic-inertial momentum transfer coupled with localized dynamic jamming.

In a static fluid, a surface depression with negative curvature experiences an unbalanced inward pressure:

$$P_{\text{hydro}} = \rho g h$$

In the shaken non-Newtonian suspension, the vertical oscillation of the base plate introduces a high-shear boundary condition along the hole’s perimeter during the upward acceleration phase ($\ddot{z} > 0$).

This vertical velocity gradient:

$$\dot{\gamma}_{zr} \approx \frac{\partial u_z}{\partial r}$$

forces adjacent cornstarch granules into direct frictional contact, triggering discontinuous shear-thickening (DST). The effective viscosity of the rim increases discontinuously, transforming the boundary fluid into a jammed, solid-like wall that sustains shear stress.

Additionally, the upward acceleration injects an outward radial momentum flux:

$$\langle \rho u_r^2 \rangle$$

away from the void center.

During the downward stroke, normal stresses diminish, but the structural relaxation time of the dense colloidal network ($\tau_{\text{relax}}$) is longer than the duration of the downward phase.

The rim does not fully relax into a low-viscosity liquid state before the subsequent upward stroke restiffens the boundary. The hole remains stable because the integrated outward acoustic radiation stress and dynamic wall stiffness balance the inward hydrostatic pressure over each cycle.

Parametric Resonance Frequencies vs. Driving Harmonic Ratios

At what harmonic frequency do the persistent holes, rim waves, and fingering projections oscillate relative to the driving frequency $f$?

The interfacial structures in vertically shaken oobleck oscillate primarily at the subharmonic frequency:

$$f_{\text{response}} = \frac{1}{2} f_{\text{drive}}$$

matching the fundamental resonance condition of the parametric Mathieu instability:

$$\left(\frac{\omega}{2}\right)$$

High-speed profilometry confirms that an elevated hole rim reaches its maximum vertical extension once every two cycles of the underlying diaphragm.

However, because the fluid’s discontinuous shear-thickening introduces sharp, step-like changes in local viscosity, the system’s temporal response departs from a pure sinusoid. The dynamic contact network transitions rapidly:

$$f(\sigma) \to 1$$

introducing high-frequency odd and even harmonics:

$$f_{\text{harmonic}} = n \left(\frac{f_{\text{drive}}}{2}\right), \quad n \in {1, 2, 3, 4, \dots}$$

into the local stress and velocity spectra.

While the fundamental spatial envelope oscillates subharmonically at $\frac{f}{2}$, the instantaneous acceleration of the fingering projections contains higher-frequency harmonic modes generated by the rapid impact of the jamming front during the upward stroke.

Differentiation from Classical Chladni Patterns and Capillary Instabilities

What distinguishes these non-Newtonian standing wave phenomena from classical Chladni plate figures and standard Faraday capillary waves?

The phenomena observed in vertically driven oobleck differ from classical Chladni figures and Newtonian Faraday waves in three primary physical respects:

1. Active Continuum versus Kinematic Sorting: Classical Chladni patterns rely on dry, non-interacting granular particles (such as quartz sand) moving kinematically across an elastic plate. The sand grains drift passively into nodal lines where vibrational amplitude is minimal ($\nabla \cdot \mathbf{u}_{\text{plate}} \approx 0$). In contrast, non-Newtonian fluid phenomena occur within a continuous colloidal suspension; the observed structures (such as persistent holes and fingers) are active topological deformations of the fluid interface itself, not passive particulate accumulations.
2. Phase Boundary Generation versus Sinusoidal Perturbation: In standard Faraday waves on Newtonian fluids (e.g., water), the entire surface forms continuous, periodic capillary-gravity wave arrays whose spatial wavelengths match the linear dispersion relation. In dense cornstarch suspensions, the fluid exhibits discontinuous shear-thickening (DST), which suppresses continuous periodic wave trains across the container. The interface bifurcates into distinct, coexisting thermodynamic and mechanical phases: unjammed, flowing liquid regions and jammed, solid-like structural boundaries.
3. Topological Defects and Energy Localization: In Newtonian capillary wave systems, standing perturbations dissipate through viscous damping once the driving amplitude is lowered. In the non-Newtonian suspension, persistent standing fluid holes act as localized dissipative solitons. They maintain their spatial identity, resist diffusion across the container, and sustain localized negative-curvature profiles that cannot exist in linear, Newtonian cymatic systems.