Live Physical & Visual Interpretation Particle Drift • 2D Plate
What Is Happening Physically

What You See in the Viewport

Active Setting Breakdown

Mathematical Physics & Continuum Mechanics Reference

1. Kirchhoff-Love Biharmonic Plate Theory (∇⁴W − k⁴W = 0)

Transverse deflection $W(x,y) \cos(\omega t)$ of a thin isotropic elastic plate of thickness $h$, Young's modulus $E$, density $\rho$, and Poisson's ratio $\nu$ is governed by the fourth-order biharmonic PDE with flexural rigidity $D$:

D ∇⁴W − ρh ω² W = 0,   D = E h³ / [ 12(1 − ν²) ]

At any arbitrary driving frequency $f$ (or across a frequency band $[f_{\min}, f_{\max}]$), the forced plate response is the Lorentzian modal admittance superposition over eigenmodes $(n, m)$:

W(x, y, f) = ∑ₙ,ₘ [ Γₙₘ² / ((f − fₙₘ)² + Γₙₘ²) ] · Φₙₘ(x, y)
2. Circular Kirchhoff Plates: Bessel Jₙ(kr) & Modified Bessel Iₙ(kr)

Because $\nabla^4 - k^4 = (\nabla^2 + k^2)(\nabla^2 - k^2)$, the radial eigenfunction on a circular plate $r \le R$ is a linear combination of ordinary Bessel functions $J_n(k_{nm}r)$ and modified Bessel functions $I_n(k_{nm}r)$:

Rₙₘ(r) = Jₙ(αₙₘ r/R) − [ Jₙ(αₙₘ) / Iₙ(αₙₘ) ] · Iₙ(αₙₘ r/R)

For a clamped edge ($W(R) = \partial_r W(R) = 0$), the dimensionless eigenvalue $\alpha_{nm}$ satisfies the exact Wronskian condition $J_n(\alpha)I'_n(\alpha) - I_n(\alpha)J'_n(\alpha) = 0$.

3. Degenerate Spherical Harmonics Yₗᵐ(θ, φ) on S²

On a spherical shell $S^2$, the Laplace-Beltrami operator $\Delta_{S^2} Y_l^m = -l(l+1)Y_l^m$ has a $(2l+1)$-dimensional degenerate eigenspace for each degree $l$. Superposing complementary orders $m \in [0, l]$ and $m' = l - m$ within the same degree $l$ preserves an exact eigenstate:

W(θ, φ) = w₂ Nₗₘ Pₗᵐ(cos θ) cos(mφ) ± w₁ Nₗ,ₗ₋ₘ Pₗˡ⁻ᵐ(cos θ) sin((l−m)φ)
4. Oₕ Cubic Acoustic Cavity & Analytical Newton-Raphson Projection

Inside a rigid cubic acoustic cavity $[-L/2, L/2]^3$ ($c_0 = 343.2\text{ m/s}$), all 6 permutations of wavevector $(n, m, k)$ share the exact Helmholtz eigenvalue $\lambda = \pi^2(n^2+m^2+k^2)/L^2$. Candidate points drawn from a 3D Halton sequence are projected onto the zero-set $W(\mathbf{x}) = 0$ using exact analytical partial derivatives:

x⁽ᵗ⁺¹⁾ = x⁽ᵗ⁾ − [ W(x⁽ᵗ⁾) / |∇W(x⁽ᵗ⁾)|² ] · ∇W(x⁽ᵗ⁾)