Local Spherical Flow Model

• Local spherical flow model is a framework that analyzes fluid dynamics on a finite, curved domain by retaining leading-order curvature and global effects.
• It generalizes the shearing box paradigm by incorporating time-dependent contraction/expansion and anisotropic background flows to capture local instabilities and turbulence.
• The approach supports high-resolution numerical simulations and analytical rescaling methods, linking global conservation laws with localized astrophysical phenomena.

A local spherical flow model provides an analytical and computational framework for studying the dynamics of fluid flows—hydrodynamic or MHD—on spherically symmetric or nearly symmetric backgrounds, by focusing on a finite, typically small, spatial domain (local box or patch). Such models retain leading-order curvature and global-scale effects but enable the high-resolution analysis of local instabilities, turbulence, and wave propagation. They generalize the shearing box paradigm to account for spherical geometry, time-dependent contraction or expansion, and anisotropic background flows. This approach is essential for astrophysical phenomena (e.g., protostellar collapse, supernovae, planet formation, helioseismology) and complements global simulations and kinetic treatments by enabling tractable local analysis with rigorous conservation properties.

1. Mathematical Formulation and Scale Separation

The foundational mathematical structure is based on expressing the governing equations (Euler, Navier–Stokes, MHD, or Boltzmann) in global spherical coordinates, then locally transforming to a co-moving patch basis. The box follows a reference spherical shell at radius $R_0(t)$, introducing scaled local coordinates $(x, y, z)$ aligned with $(\phi, \theta, r)$ via:

• $x \approx \phi$, $y \approx \theta - \pi/2$, $z = [r - R_0(t)]/L_z(t)$
• Horizontal and vertical scale factors: $a(t) \equiv R_0(t)$, $L_z(t)$, or in normalized form $b(t) = L_z(t)/R_0(t)$

The Jacobian for the box transformation is $J(t) = a^2(t) L_z(t)$. The total fluid velocity is decomposed as $u = U_{\rm bg} \hat{r} + v$, with $U_{\rm bg}$ the spherically symmetric background and $v$ the local perturbation.

Mass conservation in the local frame reads: $$D\rho = -\rho \left[ \Delta + \nabla \cdot v \right]\,, \quad D = \partial_t + v \cdot \nabla\,, \quad \Delta = 2 U_0/a + U_{R0}$$ and momentum,

$$D(\rho v_i) + \partial_j \Big[ \rho v_i v_j + \delta_{ij} \frac{p}{L_j^2} \Big] = S_i$$

$S_i$ are geometric source terms from background flow divergence, and $L_j$ are scale factors in each direction ($a$ for $x$ and $y$, $L_z$ for $z$) (Xu et al., 2024, Lynch et al., 2023).

2. Geometric Effects, Scale Factors, and Conservation Laws

Curvature and time-dependent background contraction/expansion are encoded through varying the box metric $g = \mathrm{diag}[a^2, a^2, L_z^2]$ and its evolution. The background divergence $\Delta$ modulates the conservation laws and introduces explicit angular momentum and vorticity amplification/decay: $D(J\,\omega) = 0\,,\quad \omega = \nabla \times v$ implying that vorticity amplifies as the box contracts ($J \downarrow$) and decays under expansion.

Angular momentum conservation and symmetries manifest in the invariance of the local equations under accelerating Galilean transformations $v \to v - v_g(t)$ with $dv_g/dt + 2(U_0/a)\,v_g = 0$, furnishing a local analog to global angular momentum conservation in spheres (Lynch et al., 2023).

3. Linear and Nonlinear Solution Families; Instabilities

The local spherical flow model admits both linear and nonlinear solutions that capture the essential physics of local instabilities and turbulence:

• Zonal (horizontal) and elevator (vertical) shear flows: their amplitude scales as $a^{-2}(t)$ or $b^{-2}(t)$ owing to angular momentum conservation under collapse.
• Sound waves, inertial waves, and “frozen” (non-WKB) modes: Linear perturbations satisfy time-dependent dispersion relations with effective sound speeds $c_{s,H} = c_s/a$, $c_{s,V} = c_s/L_z$ (Xu et al., 2024).
• Self-similar rarefaction or shock solutions can be constructed for background-homentropic Euler flows, with exact control of geometric source terms (Zhang, 29 Nov 2025).
• Under conformal (super-comoving) rescaling, the small-scale, isotropic limit maps directly to standard Cartesian incompressible turbulence, yielding a rescaled Kolmogorov cascade (Lynch et al., 16 Jan 2026).

4. Practical Numerical Implementations

The local spherical box framework features in several high-precision computational environments:

• Athena++ “collapsing box” coordinates, with explicit rescaling of grid and eigenvalues, operator-splitting for analytic source-term integration, and conservation to machine precision (Xu et al., 2024).
• Vielbein-based Lattice Boltzmann methods for flows on the 2-sphere, encoding geometric forces through orthonormal frames, and using Hermite quadrature for velocity-space accuracy and robust pole boundary treatments (Ambrus et al., 22 Apr 2025).
• Discrete Boltzmann Models (DBMs) in spherical symmetry employ a D3V26 velocity scheme, with geometric “force” terms arising from curvature, leading to faithful recovery of Navier–Stokes dynamics and measurement of thermodynamic nonequilibrium via higher moments (Xu et al., 2018).

Characteristic features include periodic boundary conditions; explicit time-dependent metric factors in the divergence and flux terms; anisotropic sound speed entering Riemann solvers; and modifications to CFL conditions.

5. Analytical and Rescaling Approaches

Local spherical flow models can be recast via time-dependent variable rescaling (“super-comoving” coordinates), enabling mapping of solutions between spherical and Cartesian backgrounds. The rescaled variables $\tilde{\rho}, \tilde{\mathbf{v}}, \tilde{P}$, along with new time $\tau = \int^t v_0(t')\,dt'$, transform the equations such that: $$\tilde{D}\tilde{\rho} = -\tilde{\rho}\nabla \cdot \tilde{\mathbf{v}}\,, \quad \tilde{\rho}\tilde{D}\tilde{\mathbf{v}} = -\nabla \tilde{P}$$ where $\tilde{D}$ is the Lagrangian derivative in rescaled (comoving) variables. In the isotropic monatomic case, pressure source terms vanish and the physical evolution is governed by transformed Cartesian dynamics, with physically observable growth rates for instabilities $\gamma_{\rm phys}(t) = \sigma a^{-2}(t)$ (Lynch et al., 16 Jan 2026).

6. Helioseismic and Astrophysical Applications

Local spherical flow modeling is central for interpreting solar interior flows via time-distance helioseismology. The formalism relies on finite-frequency sensitivity kernels expanded in vector spherical harmonics: $$K_i(r,\theta,\phi; x_1,x_2) = \sum_{\ell m \alpha} A^{(i,\alpha)}_{\ell m}(r;x_1,x_2) P^{\alpha}_{\ell m}(\theta,\phi)$$ and line-of-sight projections, enabling efficient inversion for flow fields from travel-time shift data. Spherical symmetry and Wigner D-rotations allow for significant computational savings (Bhattacharya, 2020).

In star and planet formation contexts, local spherical flow boxes enable highly resolved studies of turbulent fragmentation, dust-gas interactions, hydrodynamic and MHD instabilities during collapse/expansion, and can be generalized to include self-gravity and magnetic fields (Xu et al., 2024, Lynch et al., 2023).

7. Limitations and Extensions

Local spherical flow models presuppose the box size is small compared to the global radius, which constrains their domain of validity to regimes where higher-order curvature, large-scale gradients, or strong departures from spherical symmetry are negligible. Neglect of explicit viscosity, magnetic fields, and self-gravity is common but extensions are straightforward. Key assumptions may limit the ability to capture slow-growing instabilities in rapidly expanding or contracting flows because of time-dependent suppression of physical growth rates (Lynch et al., 16 Jan 2026).

The approach is not restricted to hydrodynamics; it extends to MHD (with appropriate vector and induction equations), kinetic theory (via local Boltzmann/DBM or lattice Boltzmann formalisms), and real-time inviscid flows/aeroacoustics on curved surfaces, including coupling to acoustic emission models by direct computation of pressure-force time derivatives (Sinha et al., 22 Jan 2026).

---

In summary, the local spherical flow model constitutes a unifying framework for high-resolution, curvature-aware local analysis of hydrodynamic and MHD processes on spherically symmetric or close-to-symmetric backgrounds, with deep connections to global conservation laws, rescaling symmetries, turbulence theory, and numerical implementations spanning Boltzmann, Eulerian fluid, and lattice-based methods (Lynch et al., 2023, Xu et al., 2024, Ambrus et al., 22 Apr 2025, Lynch et al., 16 Jan 2026, Xu et al., 2018, Bhattacharya, 2020, Sinha et al., 22 Jan 2026, Zhang, 29 Nov 2025).