Planar and Radial Flows in Astrophysics and Graphs

• Planar and radial flows are fundamental transport phenomena characterized by 2D surface confinement and center-directed motion, with applications in astrophysics and network theory.
• In astrophysical disks, radial flows governed by mass and angular momentum conservation influence gas accretion and metallicity gradients, typically showing velocities of 0.1–1 km/s.
• In computational contexts, planar flows on graphs yield efficient algorithms via duality principles, enabling rapid solutions in flow optimization and network problems.

Planar and radial flows constitute fundamental classes of transport phenomena in both physical systems (e.g., fluid and galaxy disks) and mathematical models (e.g., network flows on planar graphs). The term "planar flow" may refer to flows confined to a two-dimensional surface, while "radial flow" specifies outward or inward transport along radii from a center. In astrophysics, planar/radial flows mediate gas transport in galactic disks and directly impact chemical evolution. In applied mathematics and computer science, flows on planar graphs benefit from algorithmic properties and dualities that yield efficient solutions. The precise analysis of such flows combines observational, analytic, and computational methods.

1. Radial and Planar Flows in Astrophysical Disks

In the context of galactic disks, radial flows are driven by the conservation of mass and angular momentum, connecting cosmological accretion to in-plane motion. The disk is treated as a thin, rotating, axisymmetric structure characterized by a time-dependent gas surface density $\Sigma_g(R,t)$. The mass continuity equation, accounting for star formation and accretion, is

$$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$

where $u_R(R,t)$ is the radial gas velocity, $\dot\Sigma_{\mathrm{acc}}$ the local accretion rate per unit area, and $\dot\Sigma_*$ the SFR surface density (Pezzulli et al., 2015, Pezzulli et al., 2016).

Angular momentum conservation imposes constraints on the decomposition of $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$ into direct accretion and radial inflow. For a mismatch parameter $$\alpha(R,t) = 1 - v_{\mathrm{acc}}/v_{\mathrm{disc}}$$ (flat rotation curve), the induced radial velocity is

$$u_R(R,t) = -\alpha(R,t)\,R\,\frac{\dot\Sigma_{\mathrm{acc}}(R,t)}{\Sigma_g(R,t)}$$

Typical values for the Milky Way (from chemical abundance constraints) have $\alpha \simeq 0.2$--$0.3$ ($$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$0--$$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$1), yielding $$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$2–$$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$3 and producing radial metallicity gradients $$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$4 to $$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$5 in alignment with observations (Pezzulli et al., 2015, Pezzulli et al., 2016).

2. Measurement and Properties of Radial Flows in Spiral Galaxies

Empirically, radial gas flows are constrained through kinematic modeling of neutral hydrogen (H I) velocity fields. Tilted-ring decomposition fits the observed data cube with parameters for center, systemic velocity, inclination, position angle, rotation velocity $$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$6, and radial velocity $$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$7. The line-of-sight velocity is modeled as

$$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$8

Three-dimensional modeling (e.g., with 3D-BAROLO) differentiates true radial motions from disk warps by fitting $$\frac{\partial \Sigma_g}{\partial t} + \frac{1}{R} \frac{\partial}{\partial R} (R\,\Sigma_g\,u_R) = \dot\Sigma_{\mathrm{acc}} - \dot\Sigma_*$$9 using dedicated weighting schemes near the minor axis (Teodoro et al., 2021).

Across 54 local spirals, $u_R(R,t)$0 is typically a few km/s (rarely exceeding 10 km/s) and lacks a systematic radial trend. Both inflows ($u_R(R,t)$1) and outflows ($u_R(R,t)$2) occur, often with sign changes within the same disk. Monte Carlo uncertainties per ring are $u_R(R,t)$3–$u_R(R,t)$4, making detections below $u_R(R,t)$5 uncertain (Teodoro et al., 2021).

3. Mass Flow Rates and Constraints on Star Formation

The net radial mass flux is given by

$u_R(R,t)$6

with a correction for helium: total neutral gas $u_R(R,t)$7. The mean inflow rate outside the optical disk ($u_R(R,t)$8) is $u_R(R,t)$9 ($\dot\Sigma_{\mathrm{acc}}$0 median), and in the far H I outskirts, $\dot\Sigma_{\mathrm{acc}}$1. These are an order of magnitude below the typical star formation rate ($\dot\Sigma_{\mathrm{acc}}$2). Only $\dot\Sigma_{\mathrm{acc}}$3–$\dot\Sigma_{\mathrm{acc}}$4 galaxies have radial inflows at or above the SFR ($\dot\Sigma_{\mathrm{acc}}$5 significance). Therefore, in-plane radial flows, as directly measured, cannot supply sufficient gas to sustain current SFRs in most local disks (Teodoro et al., 2021).

4. Algorithmic and Mathematical Aspects of Planar and Radial Flows

In mathematical and computational contexts, planar flows are central in the study of network flows and optimization on planar graphs. For an undirected planar graph with non-negative edge capacities, the minimum $\dot\Sigma_{\mathrm{acc}}$6-$\dot\Sigma_{\mathrm{acc}}$7 cut and the corresponding maximum flow are linked via the classical max-flow/min-cut theorem.

Planarity enables dual graph constructions: every $\dot\Sigma_{\mathrm{acc}}$8-$\dot\Sigma_{\mathrm{acc}}$9 cut corresponds to a simple cycle in the dual graph separating the dual faces of $\dot\Sigma_*$0 and $\dot\Sigma_*$1. This duality underpins efficient algorithms. For instance, minimum cut and maximum flow in undirected planar graphs can be computed in $\dot\Sigma_*$2 time, by leveraging $\dot\Sigma_*$3-divisions, dense boundary metrics, and divide-and-conquer on dual paths. Dynamic data structures achieve $\dot\Sigma_*$4 update/query times for these problems (Italiano et al., 2010).

5. Special Classes: Circular Flows in Planar Graphs

Beyond classical flow values, the concept of circular $\dot\Sigma_*$5-flows—flows constrained to take values in a discrete symmetric set—arises in the study of coloring and flows in graphs. For planar graphs, the Planar Circular Flow Conjecture posits that every $\dot\Sigma_*$6-edge-connected planar graph admits a circular $\dot\Sigma_*$7-flow. Results show (i) every 10-edge-connected planar graph admits a circular $\dot\Sigma_*$8-flow (modulo 5-orientation), and (ii) every 16-edge-connected planar graph admits a circular $\dot\Sigma_*$9-flow (modulo 7-orientation). Proofs for these results operate via modulo-orientation reformulation and combine combinatorial techniques (weight functions, reducible configurations, discharging) without requiring extensive computational assistance. By planar duality, circular $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$0-flows correspond to circular colorings of the dual graph. The existence and sharpness of such flows are tightly bound to planarity properties (Cranston et al., 2018).

6. Analytical and Physical Models of Planar Radial Flows

In fluid dynamics, planar radial flows are sharply analyzable under axisymmetry, as in the braking of a vortex on a flat plate (e.g., tornado sole or vortex chamber outflow). The governing steady-state, axisymmetric radial momentum equation reads

$$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$1

with $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$2 the mean radial speed, $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$3 the wall flow thickness, and $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$4 a skin-friction coefficient. For incompressible flows with $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$5, the solution is expressible as

$$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$6

where $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$7 and $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$8 are determined by flow parameters, and $$\dot\Sigma_{\mathrm{eff}} \equiv \partial(\Sigma_g+\Sigma_*)/\partial t$$9 is fixed by boundary conditions. Analogous closed forms exist for polytropic and isothermal flows, and cumulative force profiles can be explicitly calculated, giving insight into pressure distribution and engineering design of radial-flow devices (e.g., for steam–water cloud compression) (Budarin, 2013).

7. Planar Curvature and Geometric Flows

In geometric analysis, planar (curve) flows include the evolution of closed curves under normal velocities that may depend linearly or nonlinearly on curvature. The Andrews–Bloore flow parametrizes normal velocity by $$\alpha(R,t) = 1 - v_{\mathrm{acc}}/v_{\mathrm{disc}}$$0; special cases include Eikonal and curve-shortening flows. For star-shaped curves, radial and curvature-critical points evolve via ODEs derived from the flow PDE. Computational schemes based on polar representation and spline approximation permit precise tracking of evolving shape descriptors under various flow types (Fehér et al., 2020). This formalism is distinct from mass-transport radial flows but shares foundational 2D geometry.

---

In sum, planar and radial flows provide a unifying framework, linking the physics of disks, the theory of flows on planar graphs, advanced algorithmic techniques, and geometric curve evolution. Their mathematical structure is robust against generalizations, while empirical studies in astrophysics reveal that observed radial motions, though ubiquitous, are systematically weaker than pure secular growth requires, implying the necessity of additional, non-steady modes of mass and angular momentum transport in disk systems (Teodoro et al., 2021, Pezzulli et al., 2015, Pezzulli et al., 2016).