The evolution of the flow structure around an impulsively stopped sphere is investigated in an incompressible viscous fluid under a transverse magnetic field. The study focuses on the wake structure and drag force over the range of Reynolds numbers $60 \leqslant {\textit{Re}}_{\!D} \leqslant 300$ and $ {\textit{Re}}_{\!D}=1000$, with the interaction parameters $0 \leqslant N \leqslant 10$, where $N$ characterises the strength of the magnetic field. The wake is fully developed before the impulsive stop, after which it moves downstream and interacts with the sphere under the influence of a transverse magnetic field. The complex flow structures are characterised by skin friction lines on the downstream side of the sphere and categorised into five regimes in the $\{N, {\textit{Re}}_{\!D}\}$ phase diagram based on nearly 200 cases. The drag force generally decays over time following the impulsive stop. A drag decomposition model based on the vorticity diffusion scale is proposed, attributing the drag decay to three components: the original Stokes contribution, an inertia correction at high Reynolds numbers and a magnetohydrodynamic (MHD) correction, where the inertia and MHD effects both contribute a temporal power-law decay with an exponent of $-1/6$. Temporal scaling laws of the drag decay are derived by coupling these three different effects, considering flow structures at short and long time scales, as well as the dependence on ${\textit{Re}}_{\!D}$ and $N$. The prediction results are consistent with present simulations. Furthermore, the proposed drag decomposition model is successfully extended to complex vortex flow past a sphere at ${\textit{Re}}_{\!D}=1000$, to an anisotropic ellipsoidal particle and to different magnetic field orientations.

A combined experimental and direct numerical simulation (DNS) investigation is undertaken to study the laminar boundary-layer (BL) flow adjacent to a melting vertical ice face at two far-field water salinities ($S_\infty =0$ and 34 ‰) and a range of far-field temperatures ($T_\infty$). Wall-normal distributions of vertical velocity and temperature within the BL are measured by a modified molecular tagging velocimetry and thermometry technique. Experimental data match with DNS only when a nonlinear equation of state (EoS) for density is used rather than a linear EoS. For all $S_\infty =0$, i.e. freshwater cases, the flow remains uni-directional, although the flow reverses direction at $T_\infty =4^{\,\circ} \text{C}$. A bi-directional flow, however, exists for $S_\infty =$ 34 g kg−1, where an inner salinity-driven upward flow of fresher water is accompanied by a downward-flowing temperature-driven outer flow. Although the contribution of temperature to density relative to salinity is small $({\approx}1/40)$, the thermal BL region is larger owing to higher diffusivity. This results in increased total buoyancy force when the buoyancy is integrated across the BL, which combined with effects of wall shear stress on salinity BL and a freer thermal BL growth reveals that buoyancy from temperature contributes almost equally to the overall flow. Melt rates ($V$) also show differing features in uni- and bi-directional flows. The uni-directional flows exhibit the standard scaling of increasing velocity magnitude and BL thickness, and decreasing $V$ with distance along the flow direction. Such scalings are not followed in the bi-directional flows. These show a more uniform $V$ with height, which is attributed to the counteracting effects of an upward-growing salinity BL and a downward-growing temperature BL, combined with the necessity of maintaining salinity and temperature flux balance at the ice–water interface.

Numerical simulations of turbulent flows at realistic Reynolds numbers generally rely on filtering out small scales from the Navier–Stokes equations and modelling their impact through the subgrid-scale stress tensor ${\tau }_{\textit{ij}}$. Traditional models approximate ${\tau }_{\textit{ij}}$ solely as a function of the filtered velocity gradient, leading to deterministic subgrid-scale closures. However, small-scale fluctuations can locally exhibit instantaneous values whose deviation from the mean can have a significant influence on the flow dynamics. In this work, we investigate these effects by employing direct numerical simulations combined with Gaussian filtering to quantify subgrid-scale effects and evaluating the local energy flux in both space and time. The mean performance of the canonical Clark model is assessed by conditioning the energy flux distributions on the invariants of the filtered velocity gradient tensor, $Q$ and $R$. The Clark model captures to a good degree the mean energy flux. However, the fluctuations around these mean values for given ($Q,R$) are of the order of the mean, displaying fat-tailed distributions. To be more precise, we examine the joint distributions of true energy flux and the predictions from both the Clark and the Smagorinsky models. This approach mirrors the strategy adopted in early stochastic subgrid-scale models. Clear non-Gaussian characteristics emerge from the obtained distributions, particularly through the appearance of heavy tails. The mean, the variance, the skewness and the flatness of these distributions are quantified. Our results emphasise that fluctuations are an integral component of the small-scale feedback onto the large-scale dynamics and should be incorporated into subgrid-scale modelling through an appropriate stochastic framework.

In this paper, we perform a Floquet-based linear stability analysis of the centrifugal parametric resonance phenomenon in a Taylor–Couette system subjected to a time-quasiperiodic forcing where both the inner and outer cylinders are oscillating with the same amplitude and different angular velocities given respectively by $\varOmega _0 \cos (\omega _1t)$ and $\varOmega _0 \cos (\omega _2t)$. In this context, the frequencies $\omega _1$ and $\omega _2$ are incommensurate, where the ratio $\omega _2/\omega _1$ is irrational. Taking into account non-axisymmetric disturbances, a new set of partial differential equations is derived and solved using the spectral method along with the Runge–Kutta numerical scheme. The obtained results in this framework show that this forcing triggers new and numerous reversing and non-reversing Taylor vortex flows arising via either synchronous or period-doubling bifurcations. A rich and complex dynamics is found owing to strong mode competition between these modes that alters significantly the topology of the marginal stability curves. The latter exhibit a multitude of small and condensed parabolas, giving rise to several codimension-two bifurcation points, discontinuities and cusp points in the stability diagrams. Furthermore, a proper tuning of the frequency ratio leads to a significant control of both the instability threshold and the axisymmetric nature of the primary bifurcation. Moreover, using a local quasi-steady analysis when the cylinders are slowly oscillating, intermittent instabilities are detected, characterised by spike-like behaviour in the stability diagrams with several successive growths, dampings and periods of quietness. In this limit case, the inner cylinder drive becomes the responsible forcing of the Taylor vortices’ formation where the calculated critical instability parameters correspond to those of the inner oscillating cylinder case with fixed outer cylinder. The potentially unstable regions between the cylinders are determined on the basis of the Rayleigh discriminant, where an excellent agreement with the linear stability analysis results is pointed out.