The wake of a finite wall-mounted circular cylinder of diameter $D$ and height $H$ is investigated for aspect ratios $3\leq H/D \leq 7$ and boundary layer thickness of $\delta /D \approx 0.98$ using tomographic particle image velocimetry. The Reynolds number based on $D$ is $Re = 750$. The mean wake topology is related to the evolution of the periodically shed vortices, educed from a low-order representation based on proper orthogonal decomposition of the three-dimensional velocity field. The main topological features are an arch vortex, defining the recirculating base region, and a quadrupole structure consisting of two pairs of opposite-sign vorticity concentrations extending downstream behind the obstacle-free end and wall junction. The quadrupole is the time-averaged signature of shed vortices. Vortex-tilting terms in the base region act to reorient flow-normal vorticity components streamwise, resulting in the reorientation of the ends of vortices initially shed parallel to the cylinder side walls. Through the action of the vortex-stretching terms, the bent ends connect successive vortices in a continuous chain. The influence of $H/D$ on the development of the quadrupole is characterized. The results demonstrate that the quadrupole in the mean field emerges as an imprint of the shed full-loop structures. This work reconciles mean and instantaneous interpretations satisfying the solenoidal condition on the vorticity field.

It is shown that shape anisotropy and intrinsic surface slip lead to equilibrium tilt of slippery particles in a creeping simple shear flow, even for nearly shape-isotropic particles with a cross-section that is close to circular provided the Navier-slip length is sufficiently large. We study a rigid particle with an elliptical cross-section, and of infinite extent in the vorticity direction, in simple shear. A Navier-slip boundary condition is imposed on its surface. When a Navier-slip length parameter $\lambda$ is infinite, an analytical solution is derived for the Stokes flow around a particle tilting in equilibrium at an angle $(1/2)\cos ^{-1}((1-k)/(1+k))$ to the flow direction where $0 \le k \le 1$ is the ratio of the semi-minor to semi-major axes of its elliptical cross-section. A regular perturbation analysis about this analytical solution is then performed for small values of $1/\lambda$ and a numerical continuation method implemented for larger values. It is found that an equilibrium continues to exist for any anisotropic particle $k < 1$ provided $\lambda \ge \lambda _{crit}(k)$ where $\lambda _{crit}(k)$ is a critical Navier-slip length parameter determined here. As the case $k \to 1$ of a circular cross-section is approached, it is found that $\lambda _{crit}(k) \to \infty$, so the range of Navier-slip lengths allowing equilibrium tilt shrinks as shape anistropy is lost. Novel theoretical connections with equilibria for constant-pressure gas bubbles with surface tension are also pointed out.

Direct numerical simulations have been performed for heat and momentum transfer in internally heated turbulent shear flow with constant bulk mean velocity and temperature, $u_{b}$ and $\theta _{b}$, between parallel, isothermal, no-slip and permeable walls. The wall-normal transpiration velocity on the walls $y=\pm h$ is assumed to be proportional to the local pressure fluctuations, i.e. $v=\pm \beta p/\rho$ (Jiménez et al., J. Fluid Mech., vol. 442, 2001, pp. 89–117). The temperature is supposed to be a passive scalar, and the Prandtl number is set to unity. Turbulent heat and momentum transfer in permeable-channel flow for the dimensionless permeability parameter $\beta u_b=0.5$ has been found to exhibit distinct states depending on the Reynolds number $Re_b=2h u_b/\nu$. At $Re_{b}\lesssim 10^4$, the classical Blasius law of the friction coefficient and its similarity to the Stanton number, $St\approx c_{f}\sim Re_{b}^{-1/4}$, are observed, whereas at $Re_{b}\gtrsim 10^4$, the so-called ultimate scaling, $St\sim Re_b^0$ and $c_{f}\sim Re_b^0$, is found. The ultimate state is attributed to the appearance of large-scale intense spanwise rolls with the length scale of $O(h)$ arising from the Kelvin–Helmholtz type of shear-layer instability over the permeable walls. The large-scale rolls can induce large-amplitude velocity fluctuations of $O(u_b)$ as in free shear layers, so that the Taylor dissipation law $\epsilon \sim u_{b}^{3}/h$ (or equivalently $c_{f}\sim Re_b^0$) holds. In spite of strong turbulence promotion there is no flow separation, and thus large-amplitude temperature fluctuations of $O(\theta _b)$ can also be induced similarly. As a consequence, the ultimate heat transfer is achieved, i.e. a wall heat flux scales with $u_{b}\theta _{b}$ (or equivalently $St\sim Re_b^0$) independent of thermal diffusivity, although the heat transfer on the walls is dominated by thermal conduction.

We describe a mechanism that leads to the spontaneous formation of a thermohaline staircase in the high-latitude oceans. Our analysis of this mechanism is based upon a model in which uniform gradients of temperature and salinity are assumed and is applied to a simplified mean-field model of stratified turbulence. Detailed analysis employs a parametrization of turbulent diapycnal diffusivities (Bouffard & Boegman, Dyn. Atmos. Oceans, vol. 61, 2013, pp. 14–34). This parametrization is apparently unique in that it distinguishes between the diapycnal diffusivities for heat and salt on the basis of their Prandtl (Schmidt) numbers. Our model predicts that the temperature and salinity profiles will be susceptible to linear instability if the buoyancy Reynolds number lies in the range 0.18–91, and a nonlinear mean-field model simulation demonstrates that it evolves into a well-defined thermohaline staircase that matches the characteristics of those found in the high-latitude oceans. The criterion for initial instability is furthermore shown to be consistent with the observed regional variability of staircase occurrence in the Arctic Ocean as determined by the most recent observational datasets.

Rayleigh–Taylor instability of a thin liquid film overlying a passive fluid is examined when the film is attached to a periodic wavy deep corrugated wall. A reduced-order long-wave model shows that the wavy wall enhances the instability toward rupture when the interface pattern is sub-harmonic to the wall pattern. An expression that approximates the growth constant of instability is obtained for any value of wall amplitude for the special case when the wall consists of two full waves and the interface consists of a full wave. Nonlinear computations of the interface evolution show that sliding is arrested by the wavy wall if a single liquid film residing over a passive fluid is considered but not necessarily when a bilayer sandwiched by a top wavy wall and bottom flat wall is considered. In the latter case interface tracking shows that primary and secondary troughs will evolve and subsequently slide along the flat wall due to symmetry-breaking. It is further shown that this sliding motion of the interface can ultimately be arrested by the top wavy wall, depending on the holdup of the fluids. In other words, there exists a critical value of the interface position beyond which the onset of the sliding motion is observed and below which the sliding is always arrested.

The wake of a notchback Ahmed body presenting a bi-stable nature is investigated by performing wind tunnel experiments and large-eddy simulations. Attention is confined to the Reynolds number ($Re$) influence on the wake state instability within $5\times 10^{4}\leq Re \leq 25\times 10^{4}$. Experimental observations suggest a wake bi-stability with low-frequency switches under low $Re$. The wake becomes ‘tri-stable’ with the increase of $Re$ with the introduction of a new symmetric state. The higher presence of the symmetric state can be considered as a symmetrization of the wake bi-stability with an increasing $Re$. The wake symmetry under high $Re$ attributed to the highly frequent switches of the wake is extremely sensitive to small yaw angles, showing the feature of bi-stable flows. The wake asymmetry is confirmed in numerical simulations with both low and high $Re$. The wake asymmetries are indicated by the wake separation, the reattachment and the wake dynamics identified by the proper orthogonal decomposition. However, the turbulence level is found to be significantly higher with a higher $Re$. This leads to a higher possibility to break the asymmetric state, resulting in highly frequent switches showing symmetry.

Direct numerical simulation of spanwise-rotation-driven flow transitions in viscoelastic plane Couette flow from a drag-reduced inertial to a drag-enhanced elasto-inertial turbulent flow state followed by full relaminarization is reported for the first time. Specifically, this novel flow transition begins with a drag-reduced inertial turbulent flow state at a low rotation number $0\leqslant Ro \leqslant 0.1$, and then transitions to a rotation/polymer-additive-driven drag-enhanced inertial turbulent regime, $0.1\leqslant Ro \leqslant 0.3$. In turn, the flow transitions to a drag-enhanced elasto-inertial turbulent state, $0.3\leqslant Ro \leqslant 0.9$, and eventually relaminarizes at $Ro=1$. In addition, two novel rotation-dependent drag enhancement mechanisms are proposed and substantiated. (1) The formation of large-scale roll cells results in enhanced convective momentum transport along with significant polymer elongation and stress generated in the extensionally dominated flow between adjacent roll cells at $Ro\leqslant 0.2$. (2) Coriolis-force-generated turbulent vortices cause strong incoherent transport and homogenization of significant polymer stress in the bulk via their vortical circulations at $Ro=0.5 - 0.9$.

Time-resolved particle image velocimetry measurements were performed downstream of a swept backward-facing step. The measurements allow detailed analysis of the interactions between the unsteady instabilities and the stationary crossflow vortices. Different mechanisms are identified that lead to the modulation of the different families of unsteady instabilities that occur downstream of the step. For the low-frequency spanwise-travelling mode, the modulation occurs due to a redistribution of momentum when the instability encounters regions of large spanwise shear of the wall-normal and streamwise velocity. However, the higher-frequency streamwise-travelling instabilities undergo the familiar ‘lift-up’ mechanism when they encounter the regions of large vertical velocity due to the presence of the stationary crossflow vortices. The process leading to large velocity spikes, and ultimately to a laminar breakdown to turbulence, is identified as a constructive interaction between the different unsteady instabilities, coupled with an interaction with the stationary crossflow vortices when the phases align properly.

Steady doublet-type flow takes place in a porous formation, where the log-transform $Y = \ln K$ of the spatially variable hydraulic conductivity $K$ is regarded as a stationary random field of two-point autocorrelation $\rho _Y$. A passive solute is injected at the source in the porous formation and we aim to quantify the resulting dispersion process between the two lines by means of spatial moments. The latter depend on the distance $\ell$ between the lines, the variance $\sigma ^2_Y$ of $Y$ and the (anisotropy) ratio $\lambda$ between the vertical and the horizontal integral scales of $Y$. A simple (analytical) solution to this difficult problem is obtained by adopting a few simplifying assumptions: (i) a perturbative solution, which regards $\sigma ^2_Y$ as a small parameter, of the velocity field is sought; (ii) pore-scale dispersion is neglected; and (iii) we deal with a highly anisotropic formation ($\lambda \lesssim 0.1$). We focus on the longitudinal spatial moment, as it is of most importance for the dispersion mechanism. A general expression is derived in terms of a single quadrature, which can be straightforwardly carried out once the shape of $\rho _Y$ is specified. Results permit one to grasp the main features of the dispersion processes as well as to assess the difference with similar mechanisms observed in other non-uniform flows. In particular, the dispersion in a doublet-type flow is observed to be larger than that generated by a single line. This effect is explained by noting that the advective velocity in a doublet, unlike that in source/line flows, is rapidly increasing in the far field owing to the presence there of the singularity. From the standpoint of the applications, it is shown that the solution pertaining to $\lambda \to 0$ (stratified formation) provides an upper bound for the dispersion mechanism. Such a bound can be used as a conservative limit when, in a remediation procedure, one has to select the strength as well as the distance $\ell$ of the doublet. Finally, the present study lends itself as a valuable tool for aquifer tests and to validate more involved numerical codes accounting for complex boundary conditions.

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within $6\,\%$ at all wall-normal separations for particle Reynolds numbers up to 24. The discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

Time of flight (TOF) and energy analysis in vacuum are used in series to determine jet velocity Uj, diameter dj, electrical potential Vj and energy dissipated ΔV at the breakup point of electrified nanojets of the ionic liquid 1-Ethyl-3-methylimidazolium tris(pentafluoroethyl)trifluorophosphate (EMI-FAP) (Ignat'ev et al., J. Fluorine Chem., vol. 126, issue 8, 2008, pp.1150–1159). The full spray is periodically gated by a grid held at a high voltage Vg, and received at a collector where the measured flight times provide the distribution of drop speeds u. Varying Vg provides the bivariate distribution of drop energies ξ and velocities. The collector plate, centred with the beam axis, is divided into eight concentric rings, yielding the angular distribution of the spray current, and high resolution (u,ξ) values in the whole spray. The energies of various particles of given u are all well defined, but depend uniquely on u, even though u and ξ are in principle independent experimental variables. Slow and fast particles have energies respectively well above and below the capillary voltage Ve (1.64 kV). As previously shown by Gamero-Castaño & Hruby (J. Fluid Mech., vol. 459, 2002, pp. 245–276), this behaviour is due to the 2-stage acceleration process, first jointly in the jet for all particles, and then separately for free flying drops or ions of different mass/charge. The measured two-dimensional distributions of u and ξ provide the jet velocity Uj (~0.44 km s−1) and electrical potential Vj (1.2 kV) at the breakup point. All molecular ions originate near the breakup point rather than the meniscus neck. A measurable fraction of anomalously fast drops is observed that must come from Coulomb fissions of the main drops.

Tandem configurations of two self-propelled flexible flappers of finite span are explored by means of numerical simulations. The same sinusoidal vertical motion is imposed on the leading edge of both flappers, but with a phase shift ($\phi$). In addition, a vertical offset, $H$, is prescribed between the flappers. The configurations that emerge are characterized in terms of their hydrodynamic performance and topology. The flappers reach a stable configuration with a constant mean propulsive speed and a mean equilibrium horizontal distance. Depending on $H$ and $\phi$, two different tandem configurations are observed, namely compact and regular configurations. The performance of the upstream flapper (i.e. the leader) is virtually equal to the performance of an isolated flapper, except in the compact configuration, where the close interaction with the downstream flapper (i.e. the follower) results in higher power requirements and propulsive speed than an isolated flapper. Conversely, the follower's performance is significantly affected by the wake of the leader in both regular and compact configurations. The analysis of the flow shows that the follower's performance is influenced by the interaction with the vertical jet induced by the vortex rings shed by the leader. This interaction can be beneficial or detrimental for the follower's performance, depending on the alignment of the jet velocity with the follower's vertical motion. Finally, a qualitative prediction of the performance of a hypothetical follower is presented. The model is semi-empirical, and it uses the flow field of an isolated flapper.

Modulation of fluid temperature fluctuations by particles due to thermal interaction in homogeneous isotropic turbulence is studied. For simplicity, only thermal coupling between the fluid and particles is considered, and momentum coupling is neglected. Application of the statistical theory used in cloud turbulence research leads to the prediction that modulation of the intensity of fluid temperature fluctuations by particles is expressed as a function of the Damköhler number, which is defined as the ratio of the turbulence large-eddy turnover time to the fluid thermal relaxation time. Direct numerical simulations are conducted for two-way thermal coupling between the fluid temperature field and point particles in homogeneous isotropic turbulence. The simulation results are shown to agree well with the theoretical predictions.

We report measurements of the temperature frequency spectra $P(\,f, z, r)$, the variance $\sigma ^2(z,r)$ and the Nusselt number $Nu$ in turbulent Rayleigh–Bénard convection (RBC) over the Rayleigh number range $4\times 10^{11} \underset{\smash{\scriptscriptstyle\thicksim}} { < } Ra \underset{\smash{\scriptscriptstyle\thicksim}} { < } 5\times 10^{15}$ and for a Prandtl number $Pr \simeq ~0.8$ ($z$ is the vertical distance from the bottom plate and $r$ is the radial position). Three RBC samples with diameter $D = 1.12$ m yet different aspect ratios $\varGamma \equiv D/L = 1.00$, $0.50$ and $0.33$ ($L$ is the sample height) were used. In each sample, the results for $\sigma ^2/\varDelta ^2$ ($\varDelta$ is the applied temperature difference) in the classical state over the range $0.018 \underset{\smash{\scriptscriptstyle\thicksim}} { < } z/L \underset{\smash{\scriptscriptstyle\thicksim}} { < } 0.5$ can be collapsed onto a single curve, independent of $Ra$, by normalizing the distance $z$ by the thermal boundary layer thickness $\lambda = L/(2 Nu)$. One can derive the equation $\sigma ^2/\varDelta ^2 = c_1\times \ln (z/\lambda )+c_2+c_3(z/\lambda )^{-0.5}$ from the observed $f^{-1}$ scaling of the temperature frequency spectrum. It fits the collapsed $\sigma ^2(z/\lambda )$ data in the classical state over the large range $20 \underset{\smash{\scriptscriptstyle\thicksim}} { < } z/\lambda \underset{\smash{\scriptscriptstyle\thicksim}} { < } 10^4$. In the ultimate state ($Ra \underset{\smash{\scriptscriptstyle\thicksim}} { > } Ra^*_2$) the data can be collapsed only when an adjustable parameter $\tilde \lambda = L/(2 \widetilde {Nu})$ is used to replace $\lambda$. The values of $\widetilde {Nu}$ are larger by about 10 % than the experimentally measured $Nu$ but follow the predicted $Ra$ dependence of $Nu$ for the ultimate RBC regime. The data for both the global heat transport and the local temperature fluctuations reveal the ultimate-state transitions at $Ra^*_2(\varGamma )$. They yield $Ra^*_2 \propto \varGamma ^{-3.0}$ in the studied $\varGamma$ range.

The inevitable formation of a Mach disk at the central axis of a convergent conical shock wave may suffer from fundamental changes when the flow deviates from the axisymmetric condition. In this paper, the behaviours of near-conical shocks, which are generated by a circular ring wedge of $10^{\circ }$ at typical angles of attack (AoAs), are investigated at a free stream Mach number of 6 in a shock tunnel. To reveal the characteristics and mechanism of the flow, numerical analyses are carried out under the same conditions. The results indicate that when the flow deviates from axial symmetry, the circumferential non-uniformity is remarkably intensified as the shock converges downstream. The converging centre shifts against the inclination of the incoming flow and moves to the leeward side. For a sufficiently small AoA, the formation of a Mach disk remains similar to that in the axisymmetric case, although the Mach disk shrinks in size and is slightly flattened. As the circumferential non-uniformity of the shock increases at an AoA of approximately $3^{\circ }$, a pair of kinks separate the shock surface into two discontinuous segments with the stronger shock segment on the windward side and the weaker shock segment on the leeward side. When the AoA increases further, the shrinkage of the Mach disk continuously occurs, and the Mach disk is eventually replaced by a regular reflection. The discontinuity of a convergent shock with flattening on the separated shock segments and the insufficient strength increase during the subsequent convergence are responsible for the appearance of regular reflection.

The settling velocity of porous particles in linear stratification is affected by the diffusive exchange between interstitial and ambient water. The extent to which buoyancy and interstitial mass adaptation alters the settling velocity depends on the ratio of the diffusive and viscous time scales. We conducted schlieren experiments and lattice Boltzmann simulations for highly porous (95 %) but impermeable spheres settling in linear stratification. For a parameter range that resembles marine porous particles, ‘marine aggregates’, i.e. low Reynolds numbers ($0.05\leq \textit {Re}\leq 10$), intermediate Froude numbers ($0.1\leq \textit {Fr}\leq 100$) and Schmidt number of salt ($\textit {Sc}=700$), we observe delayed mass adaptation of the interstitial fluid due to lower-density fluid being dragged by a particle that forms a density boundary layer around the particle. The boundary layer buffers the diffusive exchange of stratifying agent with the ambient fluid, leading to an enhanced density contrast of the interstitial pore fluid. Stratification-related drag enhancement by means of additional buoyancy of dragging lighter fluid and buoyancy-induced vorticity resembles earlier findings for solid spheres. However, the exchange between density boundary layer and pore fluid substantially increases stratification drag for small $\textit {Fr}$. To estimate the effect of stratification on marine aggregates settling in the ocean, we derived scaling laws and show that small particles ($\leq$0.5 mm) experience enhanced drag which increases retention times by 10 % while larger porous particle (>0.5 mm) settling is dominated by delayed mass adaptation that diminishes settling velocity by 10 % up to almost 100 %. The derived relationships facilitate the integration of stratification-dependent settling velocities into biogeochemical models.

The present paper uses the detailed flow data produced by direct numerical simulation (DNS) of a three-dimensional, spatially developing plane free shear layer to assess several commonly used turbulence models in compressible flows. The free shear layer is generated by two parallel streams separated by a splitter plate, with a naturally developing inflow condition. The DNS is conducted using a high-order discontinuous spectral element method (DSEM) for various convective Mach numbers. The DNS results are employed to provide insights into turbulence modelling. The analyses show that with the knowledge of the Reynolds velocity fluctuations and averages, the considered strong Reynolds analogy models can accurately predict temperature fluctuations and Favre velocity averages, while the extended strong Reynolds analogy models can correctly estimate the Favre velocity fluctuations and the Favre shear stress. The pressure–dilatation correlation and dilatational dissipation models overestimate the corresponding DNS results, especially with high compressibility. The pressure–strain correlation models perform excellently for most pressure–strain correlation components, while the compressibility modification model gives poor predictions. The results of an a priori test for subgrid-scale (SGS) models are also reported. The scale similarity and gradient models, which are non-eddy viscosity models, can accurately reproduce SGS stresses in terms of structure and magnitude. The dynamic Smagorinsky model, an eddy viscosity model but based on the scale similarity concept, shows acceptable correlation coefficients between the DNS and modelled SGS stresses. Finally, the Smagorinsky model, a purely dissipative model, yields low correlation coefficients and unacceptable accumulated errors.

In the coastal ocean, interactions of waves and currents with large roughness elements, similar in size to wave orbital excursions, generate drag and dissipate energy. These boundary layer dynamics differ significantly from well-studied small-scale roughness. To address this problem, we derived spatially and phase-averaged momentum equations for combined wave–current flows over rough bottoms, including the canopy layer containing obstacles. These equations were decomposed into steady and oscillatory parts to investigate the effects of waves on currents, and currents on waves. We applied this framework to analyse large-eddy simulations of combined oscillatory and steady flows over hemisphere arrays (diameter $D$), in which current ($U_c$), wave velocity ($U_w$) and period ($T$) were varied. In the steady momentum budget, waves increase drag on the current, and this is balanced by the total stress at the canopy top. Dispersive stresses from oscillatory flow around obstacles are increasingly important as $U_w/U_c$ increases. In the oscillatory momentum budget, acceleration in the canopy is balanced by pressure gradient, added-mass and form drag forces; stress gradients are small compared to other terms. Form drag is increasingly important as the Keulegan–Carpenter number $KC=U_wT/D$ and $U_c/U_w$ increase. Decomposing the drag term illustrates that a quadratic relationship predicts the observed dependences of steady and oscillatory drag on $U_c/U_w$ and $KC$. For large roughness elements, bottom friction is well represented by a friction factor ($f_w$) defined using combined wave and current velocities in the canopy layer, which is proportional to drag coefficient and frontal area per unit plan area, and increases with $KC$ and $U_c/U_w$.

Previous analyses of the flow of low-Reynolds-number, viscous gravity currents down inclined planes are investigated further and extended. Particular emphasis is on the motion of the fluid front and tail, which previous analyses treated somewhat cavalierly. We obtain reliable, approximate, analytic solutions in these regions, the accuracies of which are satisfactorily tested against our numerical evaluations. The solutions show that the flow has several time scales determined by the inclination angle, $\alpha$. At short times, the influence of initial and boundary conditions is important and the flow is governed by both the pressure gradient and the direct action of gravity due to inclination. Later on, the areas where the boundary conditions are important shrink. This fact explains why previous solutions, being inaccurate near the front and the tail, described experimental data with high accuracy. At larger times, of the order of $\alpha ^{-5/2}$, the influence of the pressure gradient may be neglected and the fluid profile converges to the square-root shape predicted in previous works. Important extensions of our approach are also outlined.

The mechanisms of leading-edge vortex (LEV) formation and its stable attachment to revolving wings depend highly on Reynolds number ($\textit {Re}$). In this study, using numerical methods, we examined the $\textit {Re}$ dependence of LEV formation dynamics and stability on revolving wings with $\textit {Re}$ ranging from 10 to 5000. Our results show that the duration of the LEV formation period and its steady-state intensity both reduce significantly as $\textit {Re}$ decreases from 1000 to 10. Moreover, the primary mechanisms contributing to LEV stability can vary at different $\textit {Re}$ levels. At $\textit {Re} <200$, the LEV stability is mainly driven by viscous diffusion. At $200<\textit {Re} <1000$, the LEV is maintained by two distinct vortex-tilting-based mechanisms, i.e. the planetary vorticity tilting and the radial–tangential vorticity balance. At $\textit {Re}>1000$, the radial–tangential vorticity balance becomes the primary contributor to LEV stability, in addition to secondary contributions from tip-ward vorticity convection, vortex compression and planetary vorticity tilting. It is further shown that the regions of tip-ward vorticity convection and tip-ward pressure gradient almost overlap at high $\textit {Re}$. In addition, the contribution of planetary vorticity tilting in LEV stability is $\textit {Re}$-independent. This work provides novel insights into the various mechanisms, in particular those of vortex tilting, in driving the LEV formation and stability on low-$\textit {Re}$ revolving wings.