2026-01-19 | | Total: 6
We present a differential geometric framework for the motion of a non-Brownian particle in the presence of fixed obstacles in a quiescent fluid, in the deterministic Stokesian regime. While the Helmholtz Minimum Dissipation Theorem suggests that the hydrodynamic resistance tensor $R_{ij}$ acts as the natural Riemannian metric of the fluid domain, we demonstrate that particle trajectories driven by constant external forces are \emph{not} geodesics of this pure resistance metric. Instead, they experience a geometric drift perpendicular to the geodesic path due to the manifold's curvature. To reconcile this, we introduce a unified geometric formalism, proving that physical trajectories are geodesics of a conformally scaled metric, $\tilde{g}_{ij} = \mathcal{D}(\mathbf{x})R_{ij}$, where $\mathcal{D}$ is the local power dissipation. This framework establishes that the affine parameter along the trajectory corresponds to the cumulative energy dissipated. We apply this theory to the scattering of a spherical particle by a fixed obstacle, showing that the previously derived trajectory of the particle is recovered as a direct consequence of the curvature of this dissipation-scaled manifold.
Analytical solutions for the yaw response of tail fins for small wind turbines, and wind vanes for wind direction measurement, are derived for any planform and any release angle $γ_0$. This extends current linear models limited to small $|γ_0|$ and low aspect ratio planforms. The equation studied here is the minimal form of the general second order equation for the yaw angle, $γ$, derived by Hammam and Wood (2023). The nonlinear damping is controlled by a small parameter that depends on the vortex flow coefficient, $K_v$, which is absent from all linear models. The minimal equation is analysed using perturbation techniques. A truncated series solution from the Krylov-Bogoliubov-Mitropolskii averaging method compares favourably with a numerical solution apart from some small deviations at large time. Another form of averaging due to Beecham and Titchener (1971) yields a compact solution in terms of the rate of amplitude decay, and the rate of change of phase angle. This allows the identification of an equivalent linear system with equivalent frequency and damping ratio. Two limiting analytic solutions for small and large $|γ_0|$ are obtained. The former is used to identify the model parameters from experimental data. Both approximate solutions showed that high $K_v$ is important for fast decay of yaw amplitude for tail fins at high $|γ_0|$. High aspect ratios for wind vanes would reduce the nonlinearity to minimize yaw error. Linear response that is independent of $K_v$ occurs whenever $\sin{(πγ_0)\approx πγ_0}$. Further, the low angle analytical solution allows an exact identification of the nonlinearity which could be used to extend the modelling of wind vanes to high $γ$.
Self-propelled phoretic swimmers are generally studied in the laminar flow regime, where their low speed renders inertial effects negligible and trajectories highly predictable. This research tackles the challenge of propulsion in the inertial regime, at moderate Reynolds numbers (100 < Re < 200), where fluid dynamics becomes non-linear. By using a chemically driven macroscopic hydrogel this work demonstrates, through experiments and modeling, the existence of stable resonant states under confined geometry: as the swimmer circles, it interacts with its own lasting chemical wake. This chemical self-feedback creates a complex, stable motion characterized by both universal exponential speed decay with superimposed significant periodic speed oscillations. Furthermore, a critical threshold speed is identified, where the system abruptly transitions from the resonant oscillatory regime to a stochastic stop & go behavior. These findings provide a fundamental understanding of how chemical fields, hydrodynamic inertia, and confinement couple to determine the motion properties of high-speed active matter.
Transonic airfoil buffet is a class of shock-wave/boundary-layer interaction (SBLI) known to exhibit self-sustained two-dimensional (2D) chordwise shock wave oscillations (Strouhal number St=0.05-0.1), and three-dimensional (3D) spanwise-modulated flow separation/reattachment (St=0.2-0.4). Due to computational cost, scale-resolving simulations of span-periodic configurations to date have been limited to narrow airfoils, insufficient to accommodate the 3D buffet cell instability reported in low-fidelity simulations and experiments. In this work, implicit large-eddy simulations (ILES) and modal analysis are performed on infinite wings up to AR=3 with sweep effects for the first time. Two flow conditions are examined, corresponding to minimally and largely separated mean flow at the shock location. For the minimally separated case, the shock dynamics remain essentially spanwise-uniform (quasi-2D), with only weak and intermittent separation cells confined to the trailing-edge region and exhibiting negligible interaction with the shock. In contrast, increased mean separation leads to the emergence of pronounced 3D buffet cells with a characteristic spanwise wavelength: 1-1.5c. Spectral proper orthogonal decomposition reveals that a stationary low-frequency 3D separation mode previously identified on unswept wings (St=0.02) becomes a spanwise travelling mode as sweep is imposed, shifting monotonically to intermediate frequencies (St=0.06-0.35). The 2D shock mode is largely insensitive to sweep, whereas the frequency and energy content of the 3D mode increase with sweep while its wavelength remains unchanged. The results demonstrate that transonic buffet arises from the superposition of distinct but coupled 2D shock motion and separation-driven 3D instabilities, with mean flow separation at the shock identified as a necessary condition for dominant 3D buffet dynamics to emerge.
Particulate matter in a fluid injected into a porous reservoir impairs its permeability spatio-temporally due to pore clogging. As particle volume fraction increases near the pore throats, inter-particle contact mechanics determine their jamming and subsequent pore clogging behavior. During contact of particles submerged in a fluid, in addition to sliding friction, a rolling resistance develops due to a several micromechanical and hydrodynamic factors. A coefficient of rolling friction is often used as a lumped parameter to characterize particle rigidity, particle shape, lubrication and fluid mediated resistance, however its direct influence on the clogging behavior is not well studied in literature. We study the effect of rolling resistance on the clogging behavior of a dense suspension at pore scale using direct numerical simulations (DNS). A discrete element method (DEM) library is developed and coupled with an open-source immersed boundary method (IBM) based solver to perform pore and particle resolved simulations. Several 3D validations are presented for the DEM library and the DEM-IBM coupling and the effect of rolling resistance on clogging at a pore entry is studied.
The time dynamics of flagellar and ciliary beating is often neglected in theories of microswimmers, with the most common models prescribing a time-constant actuation of the surrounding fluid. By explicitly introducing a metachronal wave, coarse-grained to a sinusoidal surface slip velocity, we show that a spatial resonance between the metachronal wave and the corrugation of a confining cylindrical channel enables a ciliate to swim even when it cannot move forward in a bulk fluid. Using lubrication theory, we reduce the problem to the Adler equation that reveals an oscillatory and ballistic swimming regime. Interestingly, a ciliate can even reverse its swimming direction in a corrugated channel compared to the bulk fluid.