Orientation dynamics of a settling spheroid in simple shear flow: bifurcations and stochastic alignment

Anubhab Roy; Himanshu Mishra

arxiv: 2604.13676 · v1 · submitted 2026-04-15 · ⚛️ physics.flu-dyn

Orientation dynamics of a settling spheroid in simple shear flow: bifurcations and stochastic alignment

Himanshu Mishra , Anubhab Roy This is my paper

Pith reviewed 2026-05-10 12:25 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn

keywords settling spheroidsimple shear floworientation dynamicsSNIC bifurcationstochastic alignmentJeffery orbitsPeclet number

0 comments

The pith

A single effective parameter R switches settling spheroids in shear flow from continuous rotation to fixed alignment via a saddle-node bifurcation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that when gravity lies in the shear plane the competition between Jeffery shear torques and settling torques collapses the azimuthal dynamics to overdamped motion inside a tilted periodic potential whose shape is fixed by one combined parameter R. At R equals one a saddle-node bifurcation on an invariant circle destroys the periodic orbit, after which the particle approaches a stable fixed orientation. The rotation period diverges as the inverse square root of one minus R when the transition is approached from below. In the stochastic setting, noise produces Kramers-type phase slips whose frequency drops exponentially with the Peclet number once barriers appear, in contrast to the barrier-free Jeffery case where noise simply diffuses along the orbit.

Core claim

For gravity in the shear plane the azimuthal angle obeys overdamped dynamics in a tilted periodic potential controlled by the single parameter R that combines particle shape anisotropy with settling strength. A saddle-node bifurcation on an invariant circle occurs exactly at R=1 and annihilates the rotational orbit, leaving a stable equilibrium whose basin enlarges for R greater than one. The period of the orbit diverges as (1-R)^{-1/2}. When gravity is parallel to the vorticity axis the attractor remains a periodic orbit for all values of settling strength.

What carries the argument

Overdamped motion in a tilted periodic potential controlled by the single effective parameter R

If this is right

When gravity aligns with the vorticity axis the particle continues to rotate for any settling strength.
With thermal noise, phase slips across potential barriers occur at a rate exponentially sensitive to the Peclet number once R exceeds one.
Orientation moments admit asymptotic expressions in both the small-Pe and large-Pe limits that match numerical solutions of the Fokker-Planck equation.
Langevin simulations confirm intermittent dynamics in which phase slips become rarer as barrier height or Peclet number increases.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Small changes in particle aspect ratio or density ratio could produce abrupt switches between tumbling and aligned sedimentation in practical flows.
The bifurcation structure suggests that suspension rheology or particle separation processes may exhibit sharp transitions controlled by the same effective parameter R.
The analysis could be extended by relaxing the assumption that gravity lies exactly in the shear plane to map the full three-dimensional bifurcation diagram.

Load-bearing premise

The full orientation dynamics reduce accurately to one-dimensional overdamped motion in a tilted periodic potential governed by only the single combined parameter R when gravity lies in the shear plane.

What would settle it

Experimental tracking of the average rotation period of a single settling spheroid while its aspect ratio or density is varied so that R passes through one, checking whether the period grows proportionally to (1-R)^{-1/2}.

Figures

Figures reproduced from arXiv: 2604.13676 by Anubhab Roy, Himanshu Mishra.

**Figure 2.** Figure 2: Trajectories of the orientation vector on the unit sphere when the gravity vector [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Figure 3: Trajectories of the orientation vector on the unit sphere for both prolate and [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

**Figure 4.** Figure 4: Bifurcation diagram in 𝜙 ∗ − K plane for both prolate and oblate spheroids. Here, solid lines correspond to the stable fixed points, and dashed lines indicate the unstable fixed points. The results are shown for the case when 𝛼 = 𝛽 = 𝜋/2. on the rotating orbit, arresting the rotation. Near the bifurcation, the tumbling period diverges as 𝑇 ∼ G−1 1 ∼ (B2+K2F 2 𝑝 −1) −1/2 , which is the universal scaling of … view at source ↗

**Figure 5.** Figure 5: Trajectories of one end of the orientation vector on the unit sphere. Here, [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗

**Figure 6.** Figure 6: Phase portrait diagrams in the 𝜃 − 𝜙 plane for the case in which gravity is aligned with the flow–gradient direction 𝛼 = 𝛽 = 𝜋/2 for a prolate spheroid. The blue lines are the trajectories, and the blue dot indicates the initial orientation. The Bretherton constant is fixed, B = 0.5. (a) K = 0.1, (b) K = 0.2, and (c) K = 1. positions at K = 0.06 and for fixed 𝜃 ∗ . In particular, in the third panel, two st… view at source ↗

**Figure 7.** Figure 7: Phase portrait diagrams in the 𝜃 − 𝜙 plane for the case in which gravity is aligned with the flow–gradient direction 𝛼 = 𝛽 = 𝜋/2 for an oblate spheroid. The blue lines are the trajectories, and the blue dot indicates the initial orientation. The Bretherton constant is fixed, B = −0.5. (a) K = 0.1, (b) K = 0.06, and (c) K = 1. the magnitude of K is large, since linear stability analysis may not capture its … view at source ↗

**Figure 8.** Figure 8: Phase portrait diagrams when the dynamics is dominated by an inertial torque. [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗

**Figure 9.** Figure 9: Phase portrait diagrams in 𝜃 − 𝜙 plane when gravity is aligned in a flow direction 𝛼 = 𝜋/2 and 𝛽 = 0. A Bretherton constant is fixed at |B| = 0.5. Left-hand side figures are shown for a prolate spheroid, and right-hand side figures correspond to an oblate spheroid. Trajectories are shown in a blue line, with a blue dot representing an initial condition. (a) K = 0.2 (Prolate spheroid), (b) K = 0.06 (Oblate … view at source ↗

**Figure 10.** Figure 10: Contour plot of the normalized orientation distribution function [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗

**Figure 11.** Figure 11: (a) Variation of the orientation moment h𝑝 2 𝑧 i with Pe, and (b) variation of h𝑝 2 𝑥 𝑝 2 𝑦 i with Pe, for the case in which the gravity vector is aligned with the vorticity vector (𝛼 = 0). Insets show a comparison between the asymptotic solution (solid lines) and numerical results (dash–dotted lines). All results are presented for B = 0.9. h𝑝 2 𝑧 i ≈ 1 3 + 2K𝑝 45 Pe + (4K2 𝑝 − 3B 2 ) 1890 Pe2 , and (4.26… view at source ↗

**Figure 12.** Figure 12: Contour plot of the normalized orientation distribution function [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗

**Figure 13.** Figure 13: (a) Variation of the orientation moment h𝑝 2 𝑧 i with Pe, and (b) variation of h𝑝 2 𝑥 𝑝 2 𝑦 i with Pe, for the case in which the gravity vector is aligned with the flow direction (𝛼 = 𝜋/2 and 𝛽 = 0). Insets show a comparison between the asymptotic solution (solid lines) and numerical results (dash–dotted lines). All results are presented for B = 0.9. h𝑝 2 𝑥 𝑝 2 𝑦 i ≈ 1 315 (21 + K𝑝Pe). (4.35b) In figure 1… view at source ↗

**Figure 14.** Figure 14: Azimuthal trajectories 𝜙(𝑡) from Langevin simulations of the flow-aligned configuration (𝛼 = 𝜋/2, 𝛽 = 0, B = 0.9). Columns correspond to increasing barrier height (|K𝑝 | = 0.75, 1, 2); rows to increasing Péclet number (Pe = 10, 50, 100). Phase slips (sudden 𝜋-jumps between potential wells) become progressively rarer as Pe or |K𝑝 | increases, consistent with the Kramers prediction (4.56). particle dwells n… view at source ↗

read the original abstract

We investigate the orientation dynamics of a settling spheroid in simple shear flow, combining a deterministic dynamical-systems analysis with a stochastic Fokker-Planck treatment. The dynamics is governed by the competition between the Jeffery torque from the background shear and the inertial torque from settling. For configurations in which gravity lies in the shear plane, the azimuthal dynamics reduces to overdamped motion in a tilted periodic potential controlled by a single effective parameter $\mathcal{R}$ that combines the particle shape anisotropy and the settling strength. A saddle-node bifurcation on an invariant circle (SNIC) at $\mathcal{R}=1$ governs the transition from sustained rotational motion to steady equilibrium, with the rotation period diverging as $(1-\mathcal{R})^{-1/2}$. When gravity is parallel to the vorticity axis, the attractor is a periodic orbit for all settling strengths. The stochastic analysis reveals that noise plays a fundamentally different role depending on whether settling-induced potential barriers are present: in the classical Jeffery problem it diffuses over the orbit constant, whereas with settling it drives Kramers-type phase slips whose rate is exponentially sensitive to the P\'{e}clet number, defined as the ratio of diffusive to convective time scales. Langevin simulations confirm the predicted intermittent dynamics, with phase slips becoming progressively rarer as the barrier height or P\'{e}clet number increases. Asymptotic results in both the small- and large-$\mathrm{Pe}$ limits, together with numerical solutions of the Fokker-Planck equation at arbitrary $\mathrm{Pe}$, quantify the orientation moments across all regimes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper reduces settling-spheroid azimuthal motion to overdamped dynamics on a tilted periodic potential, yielding a SNIC at R=1 and Kramers-type phase slips under noise.

read the letter

The main thing to know is that when gravity lies in the shear plane the torque balance collapses to a single effective parameter R that tilts the periodic potential. This produces a saddle-node on an invariant circle at R=1, with the rotation period diverging as (1-R)^{-1/2}, and it switches the role of noise from simple diffusion along the orbit to exponentially rare Kramers escapes whose rate depends on barrier height and Pe. The authors work out the small-Pe and large-Pe asymptotics, solve the Fokker-Planck equation numerically, and confirm the intermittent slips with Langevin runs. That combination is the concrete advance over the classical Jeffery literature they cite. The deterministic reduction and bifurcation scaling are standard once the potential is written down, but the clean mapping and the explicit contrast in stochastic behavior are useful for this geometry. The numerics appear internally consistent with the analysis, and there are no fitted constants or circular definitions. The main limitation is the restriction to gravity exactly in the shear plane for the 1D picture; the parallel-to-vorticity case is noted to remain periodic but intermediate angles are not explored in detail. The exponential sensitivity to Pe would be stronger with more reported checks on rare-event sampling, though the abstract indicates they performed the runs. This is aimed at researchers working on low-Re suspensions, rheology, and sedimentation who already know Jeffery orbits and want a mechanistic account of settling plus noise. A reader in that group will get a clear, testable picture of the transition and the intermittent dynamics. It deserves peer review because the core claims follow directly from the torque balance and standard dynamical-systems tools, and the numerics are in principle reproducible.

Referee Report

1 major / 2 minor

Summary. The manuscript analyzes the orientation dynamics of a settling spheroid in simple shear flow by combining deterministic dynamical-systems methods with stochastic Fokker-Planck and Langevin treatments. When gravity lies in the shear plane, the azimuthal dynamics reduces exactly to overdamped motion in a tilted periodic potential controlled by a single effective parameter R that combines shape anisotropy and settling strength. This produces a saddle-node bifurcation on an invariant circle (SNIC) at R=1, with the rotation period diverging as (1-R)^{-1/2}. When gravity is parallel to the vorticity axis, the attractor remains a periodic orbit for all settling strengths. The stochastic analysis distinguishes noise effects: diffusion over the orbit constant in the classical Jeffery case versus Kramers-type phase slips (exponentially sensitive to the Péclet number) when settling barriers are present. Asymptotics for small and large Pe, together with numerical Fokker-Planck solutions and Langevin confirmation, quantify the orientation moments across regimes.

Significance. If the reduction holds, the work provides a clean analytical mapping of a 3D orientation problem onto a standard 1D bifurcation with parameter-free scaling, plus a clear distinction in the role of noise between Jeffery orbits and settling-induced potentials. This is a strength for modeling applications in sedimentation, particle transport, and suspensions. The combination of exact reduction, textbook SNIC analysis, Kramers escape rates, and cross-validation by asymptotics plus numerics adds rigor and falsifiability.

major comments (1)

[Deterministic analysis section (likely §2 or §3)] The reduction of the full 3D torque balance (Jeffery plus inertial settling) to the 1D overdamped equation on the tilted circle is load-bearing for the SNIC claim and the subsequent stochastic analysis. The manuscript should state the explicit form of the effective potential and the definition of R (including how it combines the shape factor and the dimensionless settling velocity) in the deterministic section, together with the conditions under which the reduction is exact rather than approximate.

minor comments (2)

[Introduction and stochastic analysis] The Péclet number is defined as the ratio of diffusive to convective timescales; this definition should appear explicitly in the introduction or methods before its use in the stochastic sections.
[Numerical results] Figure captions for the Langevin simulations should include the number of realizations, integration timestep, and how phase slips are detected to allow direct comparison with the Kramers-rate predictions.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive summary of the work and for the constructive major comment. We address it directly below and have prepared a revised manuscript accordingly.

read point-by-point responses

Referee: [Deterministic analysis section (likely §2 or §3)] The reduction of the full 3D torque balance (Jeffery plus inertial settling) to the 1D overdamped equation on the tilted circle is load-bearing for the SNIC claim and the subsequent stochastic analysis. The manuscript should state the explicit form of the effective potential and the definition of R (including how it combines the shape factor and the dimensionless settling velocity) in the deterministic section, together with the conditions under which the reduction is exact rather than approximate.

Authors: We agree that an explicit statement of the effective potential and the definition of the control parameter strengthens the presentation of the load-bearing reduction. In the revised deterministic section we now give the explicit form of the tilted periodic potential that governs the azimuthal dynamics, together with the precise definition of the effective parameter (denoted R in the referee's notation and script-R in the manuscript) as the ratio that combines the particle shape anisotropy factor with the dimensionless settling velocity. We also state the conditions under which the reduction from the full 3D torque balance is exact: namely, when the gravity vector lies in the shear plane, the polar-angle equation admits an exact equatorial fixed point that decouples the azimuthal motion, yielding the 1D overdamped dynamics on the circle. These additions are placed immediately before the SNIC analysis so that the subsequent stochastic treatment rests on a fully transparent foundation. No other changes to the results or conclusions are required. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained from torque balance

full rationale

The paper derives the effective parameter R directly from the torque balance between Jeffery shear torque and inertial settling torque. The reduction of azimuthal dynamics to overdamped motion in a tilted periodic potential is an exact consequence of the governing equations when gravity lies in the shear plane. The SNIC bifurcation at R=1 and the period divergence scaling (1-R)^{-1/2} are then standard analytic results from 1D dynamical systems applied to that potential; they are not fitted or redefined in terms of themselves. No load-bearing self-citations, ansatz smuggling, or renaming of known results occur. The stochastic Kramers analysis follows downstream from the deterministic potential without circular dependence on the target claims. The construction remains independent of the bifurcation outcome.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on standard low-Reynolds-number hydrodynamics for rigid spheroids; no new entities are postulated and no parameters are fitted to data.

axioms (1)

domain assumption Low-Reynolds-number flow with additive Jeffery and settling torques on a rigid spheroid
Invoked throughout the deterministic reduction and stochastic treatment.

pith-pipeline@v0.9.0 · 5586 in / 1483 out tokens · 52200 ms · 2026-05-10T12:25:37.130442+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

[1]

ANAND, P.AND RAY, S. S. & SUBRAMANIAN, G. 2020 Orientation dynamics of sedimenting anisotropic particles in turbulence. Physical Review Letters 125 (3), 034501. ASOKAN, K., RAMAMOHAN, T.R. & KUMARAN, V. 2002 A novel approach to computing the orientation moments of spheroids in simple shear flow at arbitrary péclet number. Physics of Fluids 14 (1), 75–84. ...

work page 2020
[2]

CHEN, S. B. & JIANG, L. 1999 Orientation distribution in a dilute suspension of fibers subject to simple shear flow. Physics of Fluids 11 (10), 2878–2890. CHEN, S. B. & KOCH, D. L. 1996 Rheology of dilute suspensions of charged fibers. Physics of Fluids 8 (11), 2792–2807. CHERTKOV, M., KOLOKOLOV, I., LEBEDEV, V. & TURITSYN, K. 2005 Polymer statistics in a...

work page 1999
[3]

FREIDLIN, M. I. & WENTZELL, A. D. 1998 Random Perturbations of Dynamical Systems, 2nd edn. Springer. GAVZE, E. & SHAPIRO, M. 1996 Sedimentation of spheroidal particles in a vertical shear flow near a wall. Journal of Aerosol Science 27, S585–S586. GERASHCHENKO, S. & STEINBERG, V. 2006 Statistics of tumbling of a single polymer molecule in shear flow. Phys...

work page arXiv 1998

[1] [1]

ANAND, P.AND RAY, S. S. & SUBRAMANIAN, G. 2020 Orientation dynamics of sedimenting anisotropic particles in turbulence. Physical Review Letters 125 (3), 034501. ASOKAN, K., RAMAMOHAN, T.R. & KUMARAN, V. 2002 A novel approach to computing the orientation moments of spheroids in simple shear flow at arbitrary péclet number. Physics of Fluids 14 (1), 75–84. ...

work page 2020

[2] [2]

CHEN, S. B. & JIANG, L. 1999 Orientation distribution in a dilute suspension of fibers subject to simple shear flow. Physics of Fluids 11 (10), 2878–2890. CHEN, S. B. & KOCH, D. L. 1996 Rheology of dilute suspensions of charged fibers. Physics of Fluids 8 (11), 2792–2807. CHERTKOV, M., KOLOKOLOV, I., LEBEDEV, V. & TURITSYN, K. 2005 Polymer statistics in a...

work page 1999

[3] [3]

FREIDLIN, M. I. & WENTZELL, A. D. 1998 Random Perturbations of Dynamical Systems, 2nd edn. Springer. GAVZE, E. & SHAPIRO, M. 1996 Sedimentation of spheroidal particles in a vertical shear flow near a wall. Journal of Aerosol Science 27, S585–S586. GERASHCHENKO, S. & STEINBERG, V. 2006 Statistics of tumbling of a single polymer molecule in shear flow. Phys...

work page arXiv 1998