SLING Parameters in Turbulent Flows

• SLING parameters are a set of physical, dynamical, and dimensionless controls that define the onset and intensity of caustic (sling) events in particle-laden turbulent flows.
• They establish critical thresholds—such as the critical radius and Stokes number—that predict when particle trajectories decouple, leading to strong local density amplifications.
• Quantitative analyses employ metrics like vortex circulation, particle relaxation time, and flow gradients to model clustering, collision kernels, and subsequent microphysical effects.

A variety of research domains employ the term “SLING parameters,” but in the context of caustic (sling) formation in particle-laden turbulent or vortex-dominated flows, the phrase denotes the precise set of physical, dynamical, and dimensionless parameters governing the onset, occurrence, and intensity of caustic singularities (sling events) in the velocity and density fields of inertial particles. These phenomena are fundamental in understanding clustering, collision rates, and the microphysical mechanisms behind processes such as cloud droplet growth and rain initiation in atmospheric sciences. The parameters define the phase space in which fluid-inertia decoupling causes particle trajectories to cross, resulting in multi-valued velocity fields and strong local density amplifications.

1. Mathematical Definition and Physical Meaning of Sling Parameters

Sling (caustic) formation is characterized by the finite-time crossing of inertial particle trajectories, typically resulting from extreme centrifugal ejection or strong straining in turbulent or vortex-dominated flows. For small, heavy particles under Stokes drag (Maxey–Riley regime),

$$\frac{d\mathbf{r}}{dt} = \mathbf{v}, \quad \frac{d\mathbf{v}}{dt} = \frac{\mathbf{u}(\mathbf{r}, t) - \mathbf{v}}{\tau}$$

where $\mathbf{r}$ is position, $\mathbf{v}$ is particle velocity, $\mathbf{u}$ is local fluid velocity, and $\tau$ is the relaxation (Stokes) time which scales as $\tau \propto \rho_p a^2$ (particle density times radius squared). The key dimensionless group is the Stokes number: $St = \frac{\tau}{T_{flow}}$ with $T_{flow}$ a characteristic flow timescale. In the vortex context, an additional critical radius emerges: $r_{cr} \simeq 0.5 \sqrt{\Gamma \tau}$ with $\Gamma$ the vortex circulation. Only particles initially positioned within a $O(1)$-thick annulus around $r_{cr}$ experience sling events, and these are responsible for dominant local density and collision rate enhancements (Ravichandran et al., 2015).

2. Fundamental Sling Parameter List and Dimensionless Control Groups

The principal parameters governing sling/caustic formation are:

Symbol	Physical Meaning	Typical Unit/Interpretation
$\tau$	Particle relaxation time	Time (s)
$\Gamma$	Vortex circulation	Area/time ($\mathrm{m}^2/\mathrm{s}$)
$r$	Radial distance from vortex center	Length (m)
$St$	(Global) Stokes number	Dimensionless ($\tau/T_{flow}$)
$St_l$	Local Stokes number ($\tau \sigma$)	Dimensionless (velocity-gradient weighted)
$r_{cr}$	Caustic formation critical radius	Length (m), $r_{cr} \sim 0.5 \sqrt{\Gamma \tau}$
$\Delta r_{\mathrm{ann}}$	“Sling” annulus width	$0.4 r_{cr}$ (m)
$\rho_{p,\max}/\rho_{p,0}$	Max density amplification	$10^2$–$10^3$ (dimensionless ratio)

In turbulent suspensions, additional controlling groups appear:

Symbol	Physical Meaning	Typical Unit/Interpretation
$St = \tau_p/\tau_\eta$	Stokes number (DNS/turbulence)	$\tau_p$ = relaxation time, $\tau_\eta$ = Kolmogorov time
$Re$	Reynolds number	Dimensionless
$Q, R$	Flow-gradient invariants (second/third)	$Q = -\frac{1}{2}\mathrm{tr}A^2$, $R = -\frac{1}{3}\mathrm{tr}A^3$
$St_{eff}$	Filter-effective Stokes number (LES)	$\tau_p/\bar{\tau}_\eta$
$F$	Mean sling (caustic) frequency	$\tau_K^{-1}$ units

Particle inertia, fluid strain and vorticity fluctuations, and intermittent events in the tails of $Q,R$ distributions all set the frequency and intensity of sling events in high-$Re$ flows (Codispoti et al., 2024).

3. Quantitative Criteria, Thresholds, and Onset Regimes

Sling events (caustics) occur only when the relevant parameters surpass well-defined thresholds.

• Vortex-centric flows: Only particles with initial position $r_0 < r_{cr} \simeq 0.5 \sqrt{\Gamma \tau}$ can generate sling caustics, with strongest effects in the annulus $0.8 r_{cr} \leq r_0 \leq 1.2 r_{cr}$ (Ravichandran et al., 2015). At this radius the local Stokes number $St_{cr} = \Gamma \tau / r_{cr}^2 \sim O(1)$, marking the balance between inertia and local orbital timescale.
• Turbulent suspensions: The onset Stokes number for detectable sling/caustic events is $St \gtrsim 0.1$–$0.2$ (experiment (Bewley et al., 2013), simulation (Voßkuhle et al., 2013)), with an exponentially sharp transition in caustic frequency for $St<0.3$ and a $\sim St^{1/2}$ scaling for collision rates at high $St$.
• Mathematical criterion (turbulence): In high-$Re$ turbulence, caustic (sling) formation along a trajectory is triggered when the smallest eigenvalue $a(t)$ of the fluid-velocity gradient drops below $-1/(4\tau)$ for a long enough interval such that $\sqrt{\Delta} D > 1$ with $\Delta = -(a_{\min} + 1/(4\tau))$ and $D$ the excursion duration. This quantifies both the depth and persistence requirement for the "potential well" leading to finite-time gradient blow-up (Bätge et al., 2022).

4. Role in Clustering, Density Amplification, and Collision Rates

Sling parameters set the scale and localization of intense density fluctuations and, consequently, collision rates.

• Density amplification: In vortex flows, particles started in the sling annulus reach peak density $\sim 10^2$–$10^3$ above the background, as measured by Osiptsov’s method or equivalent kinematic tools (Ravichandran et al., 2015).
• Collision kernel decomposition: Turbulent suspensions admit an explicit splitting: $K = K_{cluster} + K_{sling}$, with the "sling" contribution arising exclusively from caustic-induced relative velocities. For $St \gg 1$, $K_{sling} \propto a^2 u_\eta St^{1/2}$ (Codispoti et al., 21 Jul 2025, Voßkuhle et al., 2013).
• Dominance regime: Beyond $St_c \sim 0.3$–$0.5$, $K_{sling}$ dominates, rendering the collision rate’s scaling with inertia and size $(a^2 St^{1/2})$ diagnostic of the sling mechanism (Voßkuhle et al., 2013). The rate of density and velocity singularities in the sling region is statistical, set by excursion probabilities in the $Q,R$-PDF tails and by the mean frequency of successful caustic excursions (Codispoti et al., 2024, Bätge et al., 2022).

5. Statistical Distributions and Universal Scaling Laws

• Probability density function of local Stokes number: $P(St_l = \tau_p \sigma)$ exhibits heavy, power-law negative tails for $St_l < -1$ in experiments and simulations, with exponents $\alpha \sim 3$–$4$ as global $St$ increases. The negative tail encodes the frequency of strong compressive events necessary for caustic initiation (Bewley et al., 2013).
• Universal collapse: Rescaling velocity-gradient evolution in terms of $St_l$ collapses experimental data across $St$ to the master curve $$\tau_p^2 d\sigma_0/dt = -(\tau_p \sigma_0 + (\tau_p \sigma_0)^2)$$, connecting the linear and quadratic growth regimes (Bewley et al., 2013).
• Excursion statistics (turbulence): The overall caustic/slingshot frequency $F$ along a trajectory is $F=f \mathbb{P}(\sqrt{\Delta} D>1\,|\, a<-1/4\tau)$. This combines the rare-event rate $f$ and the conditional likelihood of threshold crossing, producing an explicit $F(\mathrm{St}, Re)$ relation for high-$Re$ turbulence (Bätge et al., 2022).

6. Generalization and Application to Extended Systems

The formalism of sling parameters (critical radii or timescales, threshold Stokes numbers, annular regions, statistical properties of $St_l$ and gradient excursions) underpins modeling approaches in a variety of systems:

• Cloud microphysics: Quantitative prediction of droplet collision rates and growth bottlenecks depends on integrating the measured or calculated $K_{sling}$, with turbulence parameters ($Re$, $\epsilon$) and polydispersity factored in (Codispoti et al., 21 Jul 2025, Bätge et al., 2022).
• Filtered/LES turbulence: Introduction of the filter-effective Stokes number, $St_{eff} = \tau_p/\bar{\tau}_\eta$, accounts for mesh coarsening and LES subgrid approximation, rescaling the critical thresholds and rates to the filtered flow (Codispoti et al., 2024).
• Material ejection and singularities in elasticity: An analogous parameter regime emerges in elastica-sling mechanics, where stability boundaries, critical geometric parameters ($\alpha_1, \alpha_2, d$), and stiffness $B$ set the onset of kinematic ejection (Cazzolli et al., 2024).

7. Tabular Summary of Sling Parameters

Parameter	Definition / Formula	Physical role
$\tau$	Particle relaxation time	Sets inertia, controls $St$
$\Gamma$	Vortex circulation	Sets orbital timescale, $r_{cr}$
$St$	$\tau/T_{flow}$ (vortex) or $\tau_p/\tau_\eta$ (turb.)	Dimensionless inertia
$r_{cr}$	$0.5\sqrt{\Gamma \tau}$	Critical caustic radius
$St_l$	$\tau_p \sigma$	Local inertia-gradient product, controls singularity formation
$Re$	Reynolds number	Controls flow intermittency
$Q, R$	$-\frac{1}{2}\text{tr}A^2$, $-\frac{1}{3}\text{tr}A^3$	Map distribution of straining events
$F$	Sling event frequency	Sets caustic-driven collision rate
$g(2a)$	RDF at contact	Inertial clustering amplification
$K_{sling}$	Caustic-induced kernel [e.g., $a^2 u_\eta St^{1/2}$]	Sets high-inertia collision scaling

• "Caustics and clustering in the vicinity of a vortex" (Ravichandran et al., 2015)
• "Dissecting inertial clustering and sling dynamics in high-Reynolds number particle-laden turbulence" (Codispoti et al., 2024)
• "Observation of the sling effect" (Bewley et al., 2013)
• "Prevalence of the sling effect for enhancing collision rates in turbulent suspensions" (Voßkuhle et al., 2013)
• "Quantitative prediction of sling events in turbulence at high Reynolds numbers" (Bätge et al., 2022)
• "Polydisperse collision kernels in droplet-laden turbulence with implications for rain formation" (Codispoti et al., 21 Jul 2025)