The non-equilibrium-dissipation module (nonequilibrium_dissipation.py)

Tests the dissipation magnitude law (ledger row IV/C_ε). A developed ns3d.py field (skewness ≈-0.5) is released into free decay (forcing off); along the decay it computes the integral scale L from the spherically-binned energy spectrum, u'=√2E/3, ε=2ν Z, the Taylor microscale λ and Rey_λ=u'λ/ν, and the dissipation coefficient C_ε=ε L/u'³, logging the trajectory C_ε(Rey_λ).

)} Tests the dissipation magnitude law (ledger row IV/C_ε). A developed ns3d.py field (skewness ≈-0.5) is released into free decay (forcing off); along the decay it computes the integral scale L from the spherically-binned energy spectrum, u'=√2E/3, ε=2ν Z, the Taylor microscale λ and Rey_λ=u'λ/ν, and the dissipation coefficient C_ε=ε L/u'³, logging the trajectory C_ε(Rey_λ).

)} Tests the dissipation magnitude law (ledger row IV/C_\eps). A developed ns3d.py field (skewness \approx-0.5) is released into free decay (forcing off); along the decay it computes the integral scale L from the spherically-binned energy spectrum, u'=\sqrt{2E/3}, \eps=2\nu Z, the Taylor microscale \lambda and \Rey_\lambda=u'\lambda/\nu, and the dissipation coefficient C_\eps=\eps L/u'^3, logging the trajectory C_\eps(\Rey_\lambda). In the high-\Rey_\lambda (non-equilibrium) window C_\eps follows the Vassilicos law \propto\Rey_\lambda^{-1.01\pm0.02} at N=96 (well resolved, k_{\max}\eta>1.3; R^2=0.98) and -1.14 at N=64; as the decay becomes self-similar the slope flattens to \approx-0.3 and C_\eps approaches a plateau (\approx1.0) — the \Pi_L event-RG fixed point (\S§5), i.e. the zeroth law. A checkpointed driver (decay.py, analyze_ce.py) reproduces both resolutions; captured output ships with the module.