We propose a novel multiple-scale spatial marching method for flows with slow streamwise variation. The key idea is to couple the boundary region equations, which govern large-scale flow evolution, with local exact coherent structures that capture the small-scale dynamics. This framework is consistent with high-Reynolds-number asymptotic theory and offers a promising approach to constructing time-periodic finite-amplitude solutions in a broad class of spatially developing shear flows. As a first application, we consider a non-uniformly curved channel flow, assuming that a finite-amplitude travelling-wave solution of plane Poiseuille flow is sustained at the inlet. The method allows for the estimation of momentum transport and highlights the impact of the inlet condition on both the transport properties and the overall flow structure. We then consider a case with gradually decreasing curvature, starting with Dean vortices at the inlet. In this setting, small external oscillatory disturbances can give rise to subcritical self-sustained states that persist even after the curvature vanishes.

In the paper, we consider a two-dimensional free-surface flow past a single point vortex in fluid of infinite depth. The flow moves from left to right with uniform speed $c$ far upstream and is subject to the downward acceleration $g$ of gravity. A point vortex of circulation $\varGamma$ is located at depth $H$. The positive direction of circulation is counterclockwise. The flow is characterised by two dimensionless parameters which are the dimensionless vortex circulation $\gamma =\varGamma /(\textit{cH}\,)$ and the Froude number $ \textit{Fr}=c/\sqrt {gH}$. The goal of the paper is to find the solutions of the solitary wave type with one or several crests on the free surface. These solutions are waveless far downstream and have a vertical line of symmetry. We have established that for a fixed Froude number $ \textit{Fr}\le 0.8$, there exists a finite set of positive $\gamma$ for which the solutions of the solitary wave type occur.

Granular flow down an inclined plane is ubiquitous in geophysical and industrial applications. On rough inclines, the flow exhibits Bagnold’s velocity profile and follows the so-called $\mu (I)$ local rheology. On insufficiently rough or smooth inclines, however, velocity slip occurs at the bottom and a basal layer with strong agitation emerges below the bulk, which is not predicted by the local rheology. Here, we use discrete element method simulations to study detailed dynamics of the basal layer in granular flows down both smooth and rough inclines. We control the roughness via a dimensionless parameter, $R_a$, varied systematically from 0 (flat, frictional plane) to near 1 (very rough plane). Three flow regimes are identified: a slip regime ($R_a \lesssim 0.45$) where a dilated basal layer appears, a no-slip regime ($R_a \gtrsim 0.6$) and an intermediate transition regime. In the slip regime the kinematics profiles (velocity, shear rate and granular temperature) of the basal layer strongly deviate from Bagnold’s profiles. General basal slip laws are developed that express the slip velocity as a function of the local shear rate (or granular temperature), base roughness and slope angle. Moreover, the basal layer thickness is insensitive to flow conditions but depends somewhat on the interparticle coefficient of restitution. Finally, we show that the rheological properties of the basal layer do not follow the $\mu (I)$ rheology, but are captured by Bagnold’s stress scaling and an extended kinetic theory for granular flows. Our findings can help develop more predictive granular flow models in the future.

The flow past a $6:1$ prolate spheroid at a moderate pitch angle $\alpha =10^\circ$ is investigated with a focus on the turbulent wake in a high-fidelity large eddy simulation (LES) study. Two length-based Reynolds numbers, ${\textit{Re}}_L=3\times 10^4$ and $9\times 10^4$, and four Froude numbers, ${\textit{Fr}} = \infty \text{(unstratified)}, 6, 1.9 \text{ and }1$, are selected for the parametric study. Spectral proper orthogonal decomposition (SPOD) analysis of the flow reveals the leading coherent modes in the unsteady separated flow at the tail of the body. At the higher ${\textit{Re}}_L=9\times 10^4$, a high-frequency spanwise flapping of shear layers on either side of the body is observed in the separated boundary layer for all cases. The flapping does not perturb the lateral symmetry of the wake. At ${\textit{Fr}}=\infty$, a low-frequency oscillating laterally asymmetric mode, which is found in addition to the shear-layer mode, leads to a sidewise unsteady lateral load. All temporally averaged wakes at ${\textit{Re}}=9\times 10^4$ are found to be spanwise symmetric in the mean as opposed to the lower ${\textit{Re}}=3\times 10^4$, at which the ${\textit{Fr}}=\infty \text{ and }6$ wakes exhibit asymmetry. The turbulent kinetic energy (TKE) budget is compared among cases. Here, ${\textit{Fr}}=\infty$ exhibits higher production and dissipation compared with ${\textit{Fr}}=6 \text{ and }1.9$. The streamwise vortex pair in the wake induces a significant mean vertical velocity ($U_z$). Therefore, in contrast to straight-on flow, the terms involving gradients of $U_z$ matter to TKE production. Buoyancy reduces $U_z$ and also the Reynolds shear stresses involving $u^{\prime}_z$. Through this indirect mechanism, buoyancy exerts control on the wake TKE budget, albeit being small relative to production and dissipation. Buoyancy, through the baroclinic torque, is found to qualitatively affect the streamwise vorticity. In particular, the primary vortex pair is extinguished in the intermediate wake and two new vortex pairs form with opposite-sense circulation relative to the primary.

This paper investigates the transient characteristics of uniform momentum zones (UMZs) in a rapidly accelerating turbulent pipe flow using direct numerical simulation datasets starting from an initial friction Reynolds number ($Re_{\tau 0}) = 500$ up to a final friction Reynolds number ($Re_{\tau 1}) = 670$. Instantaneous UMZs are identified following the identification methodology proposed by Adrian et al. (2000 J. Fluid Mech. vol. 422, pp. 1–54). The present results reveal that, as the flow rapidly accelerates, the average number of UMZs drops. However, as the flow recovers, it is regained. This result is complemented by the temporal evolution of the average number of internal shear layers. The temporal evolution of UMZs reveals that UMZs sustain their hierarchical flow arrangement with slower zones near the wall and faster zones away from the wall throughout the rapid turbulent flow acceleration. The results show that UMZs speed up during the inertial and pre-transition phases, and progressively slow down during the transition and core-relaxation stages. It is also revealed that UMZs near the wall respond first to flow instability and show earlier signs of recovery based on UMZ kinematic results. Finally, the dominant quadrant behaviour of Reynolds shear stress within UMZs has been investigated. It is found that, prior to the flow excursion, the UMZs nearest to the wall are always $Q2$ dominated, while the rest of the UMZs are always $Q4$ dominated. This behaviour is detected to not change during and after the flow excursion, suggesting that this is a characteristic behaviour of UMZs in accelerating turbulent wall-bounded flows.

For smectic C* (SmC*) liquid crystals, configured in a bookshelf-type geometry between two horizontal parallel plates, with the bottom plate fixed and the top plate free to move, it is known from experiment that pumping can occur when an electric field is applied, i.e. an upward movement of the top plate through mechanical vibrations when the electric field is suddenly reversed. In this paper we revisit an earlier mathematical model for fast electric field reversal by removing an assumption made there on the velocity field; instead, we arrive at a time-dependent, two-dimensional squeeze-film model, which can ultimately be formulated in terms of a highly nonlinear integro-differential equation. Subsequent analysis leads to an unexpected solvability condition involving the five SmC* viscosity coefficients regarding the existence and uniqueness of solutions. Furthermore, we find that, when solutions do exist, they imply that the plate can move down as well as up, with the final resting position turning out to be dependent on the initial conditions; this is in stark contrast to the results of the earlier model.

Permanent gravity waves propagating in deep water, spanning amplitudes from infinitesimal to their theoretical limiting values, remain a classical yet challenging problem due to its inherent nonlinear complexities. Traditional analytical and numerical methods encounter substantial difficulties near the limiting wave condition due to singularities at sharp wave crests. In this study, we propose a novel hybrid framework combining the homotopy analysis method (HAM) with machine learning (ML) to efficiently compute convergent series solutions of Stokes waves in deep water for arbitrary wave amplitudes from small to theoretical limiting values, which show excellent agreement with established benchmarks. We introduce a neural network trained using only 20 representative cases whose series solution are given by means of HAM, which can rapidly predict series solutions across arbitrary steepness levels, substantially improving computational efficiency. Additionally, we develop a neural network to gain the inverse mapping from the conformal coordinates $(\theta , r)$ to the physical coordinates $(x,y)$, facilitating explicit and intuitive representations of series solutions in physical plane. This HAM–ML hybrid framework represents a powerful and efficient approach to compute convergent series in a whole range of physical parameters for water waves with arbitrary wave height including even limiting waves. In this way we establish a new paradigm to quickly obtain convergent series solutions of complex nonlinear systems for a whole range of physical parameters, thereby significantly broadening the scope of series solutions that can be easily gained by means of HAM even for highly nonlinear problems in science and engineering.

Magneto-gravity-precessional instability, which results from the excitation of resonant magneto-inertia-gravity (MIG) waves by a background shear generic to precessional flows, is addressed here. Two simple background precession flows, that of Kerswell (1993 Geophys. Astrophys. Fluid Dyn. vol 72, no. 1–4, pp. 107–144), and that of Mahalov (1993 Phys. Fluids A: Fluid Dyn. vol. 5, no. 4, pp. 891–900), are considered. We analytically perform an asymptotic analysis to order ${ O}(\varepsilon ),$ where $\varepsilon$ denotes the Poincaré number, i.e. the precession parameter, and determine the maximum growth rate of the destabilizing subharmonic resonances of MIG waves: that between two fast modes, that between two slow modes and that between a fast mode and a slow mode (mixed modes). The domains of the $(K_0 B_0/\varOmega _0, N/\varOmega _0)\hbox{-}$plane for which this instability operates are identified, where $1/K_0$ denotes a characteristic length scale, $B_0$ is the unperturbed Alfvén velocity, $\varOmega _0$ is the rotation rate and $N$ denotes the Brunt–Väisälä frequency. We demonstrate that the $N\rightarrow 0$ limit is, in fact, singular (discontinuous). At large $K_0B_0/\varOmega _0,$ stable stratification acts to suppress the destabilizing resonance between two fast modes as well as that between two slow modes, whereas it revives the destabilizing resonance between a fast mode and a slow mode provided $N\lt \varOmega _0,$ because, without stratification, the maximal growth rate of this instability approaches zero as $K_0B_0/\varOmega _0\rightarrow +\infty .$ This would be relevant for the generation of the mean electromotive force, and hence, the $\alpha \hbox{-}$effect in helical magnetized precessional flows under weak stable stratification. Diffusive effects on the instability is considered in the simple case where the magnetic and thermal Prandtl numbers are both equal to one.

We examine the linear stability of a shear flow driven by wind stress at the free surface and rotation at the lower boundary, mimicking oceanic flows influenced by surface winds and the Earth’s rotation. The linearised eigenvalue problem is solved using the Chebyshev spectral collocation method and a long-wave asymptotic analysis. Our results reveal new long-wave instability modes that emerge for non-zero rotational Reynolds numbers. It is observed that the most unstable mode, characterised by the lowest critical parameters, corresponds to long-wave spanwise disturbances with vanishing streamwise wavenumber. The asymptotic analysis, which shows excellent agreement with numerical results, analytically confirms the existence of this instability. Thus, the present study demonstrates the hitherto unreported combined influence of wind stress and the Earth’s rotation on ocean dynamics.

Two-way diffusion equations arising in kinetic problems relating to electron scattering and in Brownian particle dynamics present singularities absent from conventional diffusion equations. Although calculations by Stein & Bernstein, and Fisch & Kruskal have revealed the formation of entry and exit slope discontinuities at the critical points where the velocity changes sign, the analytical structure of these discontinuities remains unclear. Here we fill this gap via a local similarity variable analysis, illustrated through the two-way diffusion equation $y \partial n/\partial x=\partial ^2 n/\partial y^2$ in $-1 \leq y \leq 1$; $0 \leq x \leq L$, with $n(x,\pm 1)=0$ with various entry conditions $n(0,y)_{y\gt 0}$, and the exit condition $n(L,y)_{y\lt 0}=0$. The similarity variable $\eta =y/x^{1/3}$ permits the analytical characterization of the entry discontinuity, except for constants determined by matching with numerical solutions obtained with two numerical schemes: separation of variables following the construction of Beals, or finite-difference discretization of the transient partial differential equation, which converges in time to a solution almost identical to the separation of variables solution. Although the slope discontinuity depends markedly on the initial condition $n(0,y)_{y\gt 0}$, a simple general similarity solution structure emerges empirically, always involving a spontaneous singular contribution $C |y|^{1/2}$ at $x=0,y\lt 0$. Slow convergence of both numerical solutions near $\{x,y\}=\{0,0\}$ is attributed to the poor eigenfunction representation of the ever-present singular solution component $|y|^{1/2}$. The similarity approach applies equally to other two-way diffusion equations when the coefficient of $\partial n/\partial x$ changes sign linearly with $y$. It can also be extended to situations where this coefficient is discontinuous at the critical points.

The convection velocity in high-Reynolds-number pipe flow was investigated using two-point correlations obtained from two laser Doppler velocimetry systems. The Reynolds number ranged from ${\textit{Re}}_{{\tau}}=3000$ to 20 800, and profiles were obtained from $y/R=0.002$ up to the pipe centre, where $R$ is the pipe radius. This study examines the scaling behaviour of convection velocity profiles derived from raw velocity signals, and the convection characteristics of very large-scale motions (VLSMs) and large-scale motions extracted via scale-separated or time-resolved velocity signals. The profiles show that convection velocities from raw signals exceed the local mean velocity near the wall and gradually approach it toward the centre. These profiles can be scaled using inner variables, namely $y^+$ and $\Delta x^+$, where $\Delta x^+$ represents the measurement distance. Scale-separated convection velocities for VLSM-scale structures – defined as those larger than $5R$ – were higher than the unfiltered values and remained nearly constant up to $y^+ \leq 2000$ at ${\textit{Re}}_{{\tau}} \approx 20\,000$. In this constant region, the convection velocity of VLSMs scaled well with the bulk velocity $U_{\textit{b}}$, taking values of approximately $0.85U_{\textit{b}}$. Furthermore, analysis of the time-resolved data highlights that, when applying Taylor’s frozen turbulence hypothesis, it is essential to consider both the scale dependence and the temporal fluctuations of the convection velocity, which reflect the underlying spatio-temporal dynamics of the flow structures. The present study provides valuable data for discussions on converting frequency-domain measurements into wavenumber space using Taylor’s hypothesis.

Surface roughness is often present in flight systems travelling at high speeds, but its interaction with compressible turbulence is not well understood. Using direct numerical simulations, we study how prism-shaped roughness influences supersonic turbulent boundary layers at a free-stream Mach number $M_\infty =2$. The dataset includes four simulations featuring cubic- and diamond-shaped elements in aligned and staggered configurations. All cases have an initial smooth region where a fully turbulent boundary layer transitions to a rough wall with positively skewed roughness elements relative to the smooth-wall zero plane. This causes a sudden boundary layer growth at the smooth-to-rough transition, generating an oblique shock wave. Individual roughness elements downstream do not generate shock or expansion waves, as they do not protrude into the supersonic region. For cubical elements, the staggered arrangement increases drag and produces more pronounced boundary layer growth than the aligned case. Rotating the cubes along their vertical axis further enhances these effects, yielding the highest drag. Interestingly, diamond-shaped elements in a staggered arrangement exhibit a dynamics similar to aligned cubes, producing lower drag than other cases. We explain the relative drag induced by each roughness shape by examining viscous and pressure drag components separately. The analysis reveals that, for staggered diamonds, the flow skims more easily over roughness, drastically reducing recirculation in troughs and gaps. In other cases, wake interactions are more prominent, causing spikes of highly positive and negative skin friction, a feature often neglected in reduced-order model formulations.

Droplet vaporisation can exhibit distinct shrinkage kinetic laws depending on the experimental set-up and ambient conditions. In this work, we present a unified approach that combines experiment and theory to identify true shrinkage kinetics across a broad range of droplet vaporisation processes extending beyond the classical D2-law – particularly under realistic conditions involving support fibres or/and inevitable convective effects. Experimentally, we assume a power law $D^n= D_0^n- \textit{Kt}$, where K is the vaporisation rate constant, and re-express it as $(D/D_0)^n = 1 - t/t_{\textit{life}}$ in terms of the normalised droplet diameter $D/D_0$ and time t$ / $tlife relative to the droplet’s initial diameter D0 and lifetime tlife. Taking D as the diameter of a volume-equivalent sphere, the exponent n can be reliably extracted from the slope of the log–log plot of $( 1 - t/t_{\textit{life}})$ against $D/D_0$. The robustness of this method is demonstrated by re-confirming the D2-law for pure fuel droplet evaporation and validating the $D^{3/2} $-law for droplet evaporation under forced convection. We further apply this method to droplet combustion, revealing a significant departure from the D2-law with n $=$ 2.56 ± 0.20–2.65 ± 0.17 across various liquid fuels, unaffected by the presence of support fibres. An even more pronounced departure, with n approaching 3, is observed in droplet combustion within a continuous flame sustained by an auxiliary burner. Theoretically, we develop a more general theory to describe these droplet combustion processes, showing that the observed positive departures mainly result from flame-driven buoyant convection with 2.33 < n < 3, capturing well the experimental data. The same theoretical framework can also account for the negative departures in convection-driven vaporisation processes without flame, thereby providing a unified interpretation for the fundamental distinctions between flame-driven and non-flame-driven droplet vaporisation processes. The present study not only identifies distinct shrinkage power laws that emerge from complexities in these processes, but also reveals the central role of an inherent length scale – arising from underlying convective mechanisms – in shaping the true shrinkage kinetics that lead to violations of the D2-law.