Debris flows are a growing natural hazard as a result of climate change and population density. To effectively assess this hazard, simulating field-scale debris flows at a reasonable computational cost is crucial. We enhance existing debris flow models by rigorously deriving a series of depth-averaged shallow models with varying complexities describing the behaviour of grain–fluid flows, considering granular mass dilatancy and pore fluid pressure feedback. The most complete model includes a mixture layer with an upper fluid layer, and solves for solid and fluid velocity in the mixture and for the upper fluid velocity. Simpler models are obtained by assuming velocity equality in the mixture or single-layer descriptions with a virtual thickness. Simulations in a uniform configuration mimicking submarine landslides and debris flows reveal that these models are extremely sensitive to the rheology, the permeability (grain diameter) and initial volume fraction, parameters that are hard to measure in the field. Notably, velocity equality assumptions in the mixture hold true only for low permeability (corresponding to grain diameter $d=10^{-3}$ m). The one-layer models’ results can strongly differ from those of the complete model, for example, the mass can stop much earlier. One-layer models, however, provide a rough estimate of two-layer models when permeability is low, initial volume fraction is far from critical and the upper fluid layer is very thin. Our work with uniform settings highlights the need of developing two-layer models accounting for dilatancy and for an upper layer made either of fluid or grains.

We investigate the effects of thermal boundary conditions and Mach number on turbulence close to walls. In particular, we study the near-wall asymptotic behaviour for adiabatic and pseudo-adiabatic walls, and compare to the asymptotic behaviour recently found near isothermal cold walls (Baranwal et al. 2022. J. Fluid Mech. 933, A28). This is done by analysing a new large database of highly-resolved direct numerical simulations of turbulent channels with different wall thermal conditions and centreline Mach numbers. We observe that the asymptotic power-law behaviour of Reynolds stresses as well as heat fluxes does change with both centreline Mach number and thermal condition at the wall. Power-law exponents transition from their analytical expansion for solenoidal fields to those for non-solenoidal field as the Mach number is increased, though this transition is found to be dependent on the thermal boundary conditions. The correlation coefficients between velocity and temperature are also found to be affected by these factors. Consistent with recent proposals on universal behaviour of compressible turbulence, we find that dilatation at the wall is the key scaling parameter for these power-law exponents, providing a universal functional law that can provide a basis for general models of near-wall behaviour.

Turbulent flames in practical devices are subject to a superposition of broadband turbulence and narrowband harmonic flow oscillations. In such cases, flames have a superposition of space–time correlated wrinkles, superposed with broadband turbulent disturbances that interact nonlinearly. This paper extends our prior experimental work to characterise and quantify these flame dynamics. We extract ensemble-averaged flame edge and velocity by ensemble-averaging the instantaneous data at the same phase with respect to the forcing cycle. This paper shows that the ensemble-averaged spatio-temporal dynamics of the flame changes significantly with turbulence intensity. From a spatial viewpoint, the ensemble-averaged flame at weak turbulence intensities exhibits clear cusps and a large ratio between curvature in concave and convex regions. In contrast, at high turbulence intensities, the concave and convex parts of the ensemble-averaged flame are nearly symmetric. From a temporal viewpoint, increasing turbulence intensity monotonically suppresses higher harmonics of the forcing frequency that are manifestations of flame nonlinearities. Taken together, these both point to the interesting observation that the ensemble-averaged flame exhibits increasingly linear dynamics with increasing turbulence intensities, in contrast to its very strong nonlinear behaviours at weak turbulence intensities and juxtaposed with the increasingly nonlinear nature of its instantaneous dynamics with increasing turbulence intensity. In addition, prior studies have shown clear coherent modulation of turbulent flame speed correlated with coherent curvature modulation and that this relationship could be quantified via a ‘turbulent Markstein number’, $M_{T}$. We develop correlations for $M_{T}$ showing how it scales with turbulent and narrowband disturbance quantities, such as turbulent flame brush thickness and convective length scale.

Turbulent mixing in a supercritical CO$_2$ shear layer is examined using both experimental and numerical methods. Boundary conditions are selected to focus on the rarely studied near-critical regime, where thermophysical properties vary nonlinearly with respect to temperature and pressure. Experimental results are obtained via Raman spectroscopy and shadowgraphy, while numerical results are obtained via direct numerical simulation. The shear layer growth rate is found to be 0.2. Additionally, density profiles indicate a relaxation of density gradients between the mixed fluid and heavy fluid as the flow evolves downstream, which runs counter to existing supercritical shear layer data in the literature. The computational results identify significant anisotropy in the turbulence in the shear layer, which is discussed in terms of the development of regions of high density gradient magnitude. The Reynolds-averaged enstrophy budget at various streamwise locations indicates no significant dilatational or baroclinic contribution within the shear layer.

It is generally accepted that the evolution of the deep-water surface gravity wave spectrum is governed by quartet resonant and quasi-resonant interactions. However, it has also been reported in both experimental and computational studies that non-resonant triad interactions can play a role, e.g. generation of bound waves. In this study, we investigate the effects of triad and quartet interactions on the spectral evolution, by numerically tracking the contributions from quadratic and cubic terms in the dynamical equation. In a finite time interval, we find that the contribution from triad interactions follows the trend of that from quartet resonances (with comparable magnitude) for most wavenumbers, except that it peaks at low wavenumbers with very low initial energy. This result reveals two effects of triad interactions. (1) The non-resonant triad interactions can be connected to form quartet resonant interactions (hence exhibiting the comparable trend), which is a reflection of the normal form transformation applied in wave turbulence theory of surface gravity waves. (2) The triad interactions can fill energy into the low-energy portion of the spectrum (low wavenumber part in this case) on a very fast time scale, with energy distributed in both bound and free modes at the same wavenumber. We further analyse the latter mechanism using a simple model with two initially active modes in the wavenumber domain. Analytical formulae describing the distribution of energy in free and bound modes are provided, along with numerical validations.

Thermal Marangoni effects play important roles in bubble dynamics such as bubbles generated by water electrolysis or boiling. As macroscopic bubbles often originate from nucleated nanobubbles, it is crucial to understand how thermocapillarity operates at the nanoscale. In this study, the motion of transient bulk gas nanobubbles in water driven by a vertical temperature gradient between two solid plates is investigated using molecular dynamics simulations and analytical theory. The simulation results show that due to the thermal Marangoni force, nanobubbles move towards the hot plate at a constant velocity, similar to the behaviour of macroscale gas bubbles. However, unlike macroscale gas bubbles whose thermal conductivity and viscosity are negligible compared to those of water, the thermal conductivity and viscosity of nanoscale gas bubbles are significantly increased due to their large gas density. The thermal resistance and the slip length are also found to matter at the liquid–gas interface, though they decrease with increasing gas densities. The previous thermocapillary theory for macroscale bubbles is extended by considering all these nanoscopic effects. Expressions of the Marangoni force and the drag force are derived. By balancing the Marangoni force and the drag force, the theoretical velocity of the nanobubble migration in a thermal gradient is obtained. When using the measured transport properties of liquid, gas, and their interfaces, the theoretically obtained velocity is consistent with the result of the molecular simulations. We find that the slip length is too small to have considerable effects on nanobubble motions in the current liquid–gas system.

This paper presents detailed analyses of the Reynolds stresses and their budgets in temporally evolving stratified wakes using direct numerical simulation. Ensemble averaging is employed to mitigate statistical errors in the data, and the results are presented as functions of both the transverse and vertical coordinates – at time instants across the near-wake, non-equilibrium, and quasi-two-dimensional regimes for wakes in weakly and strongly stratified environments. Key findings include the identification of dominant terms in the Reynolds stress transport equations and their spatial structures, the generation and destruction processes of the Reynolds stresses, and the energy transfer between the Reynolds stress and the mean flow. The study also clarifies the effects of the Reynolds number and the Froude number. Additionally, we assess the validity of the eddy-viscosity type models and some existing closures for the Reynolds stress model, highlighting the limitations of isotropy and return-to-isotropy hypotheses in stratified flows.

When turbulent convection interacts with a turbulent shear flow, the cores of convective cells become aligned with the mean current, and these cells (which span the height of the domain) may interact with motions closer to the solid boundary. In this work, we use coupled Eulerian–Lagrangian direct numerical simulations of a turbulent channel flow to demonstrate that, under conditions of turbulent mixed convection, interactions between motions associated with ejections and low-speed streaks near the solid boundary and coherent superstructures in the interior of the flow interact and lead to significant vertical transport of strongly settling Lagrangian particles. We show that the primary suspension mechanism is associated with strong ejection events (canonical low-speed streaks and hairpin vortices characterised by $u'\lt 0$ and $w'\gt 0$, where $u'$ and $w'$ are the streamwise and vertical turbulent velocity fluctuations), whereas secondary suspension is strongly associated with large-scale plume structures aligned with the mean shear (characterised by $w'\gt 0$ and $\theta '\gt 0$, where $\theta$ represents temperature fluctuations). This coupling, which is absent in the limiting cases (pure channel flow or free convection) is shown to lead to a sudden increase in the interior concentration profiles as ${Ri}_\tau$, the friction Richardson number, increases, resulting in concentrations that are larger by roughly an order of magnitude at the channel midplane.

Estimation of near-wall turbulence in channel flow from outer observations is investigated using adjoint-variational data assimilation. We first consider fully resolved velocity data, starting at a distance from the wall. By enforcing the estimated flow to exactly satisfy the Navier–Stokes equations, we seek a statistically stationary turbulent state that reproduces the instantaneous outer measurements. Such an estimated state provides full access to the unknown near-wall turbulence, including the wall shear stresses and pressure. When the first observation is within 50 viscous units from the wall, the correlation coefficient between the true and estimated state exceeds 95 %. As the observations are further separated from the wall, at 90 viscous units, the accuracy of the assimilated wall stresses decreases to 40 % at the wall. This trend is nearly independent of the Reynolds number. The Fourier spectrum of the estimation error is qualitatively consistent with the coherence spectrum between the outer and the inner state variables: observed long wavelength structures in the outer flow have deeper coherence into the unobserved near-wall region and, therefore, the error is lowest at large scales. Nevertheless, the adjoint-variational approach provides a more rigorous quantification of the capacity to accurately predict the instantaneous near-wall turbulence from outer measurements. Lastly, we demonstrate the robustness of the estimation accuracy using filtered and sub-sampled outer observations.

An important feature of the dynamics of double-diffusive fluids is the spontaneous formation of thermohaline staircases, where wide regions of well-mixed fluid are separated by sharp density interfaces. Recent developments have produced a number of one-dimensional reduced models to describe the evolution of such staircases in the salt fingering regime relevant to mid-latitude oceans; however, there has been significantly less work done on layer formation in the diffusive convection regime. We aim to fill this gap by presenting a new model for staircases in diffusive convection based on a regularisation of the $\gamma$-instability (Radko 2003 J. Fluid Mech. vol. 805, 147–170), with a range of parameter values relevant to both polar oceans and astrophysical contexts. We use the results of numerical simulations to inform turbulence-closure parametrisations as a function of the horizontally averaged kinetic energy $e$, and ratio of the haline to thermal gradients $R_0^*$. These parametrisations result in a one-dimensional model that reproduces the critical value of $R_0^*$ for the layering instability, and the spatial scale of layers, for a wide range of parameter values, although there is a mismatch between the range of $R_0^*$ for layer formation in the model and observational values from polar oceans. Staircases form in the one-dimensional model, evolving gradually through layer merger events that closely resemble simulations.

Wave impact on solid structures is a well-studied phenomenon, but almost exclusively for the case that the impacting liquid (e.g. water) is surrounded by a non-condensable gas (such as air). In this study we turn to wave impact in a boiling liquid, a liquid that is in thermal equilibrium with its own vapour, which is of key relevance to the transport of cryogenic liquids, such as liquified natural gas and liquid hydrogen in the near future. More specifically, we use the Atmosphere facility at MARIN, NL, to prepare water/water vapour systems at different temperatures along the vapour curve. Here, we perform wave impact experiments by generating a soliton in a flume contained within the autoclave of the facility. A bathymetry profile interacts with the soliton, leading to a breaking wave that impacts onto a vertical wall, where we measure the pressures occurring during impact by means of $100$ embedded pressure sensors. In boiling liquids, we report wave impact pressures that are up to two orders of magnitude larger than those measured in comparable water–air experiments. We trace these pressures back to the collapse of the entrapped vapour pocket, which we semi-quantitatively describe using a simplified hemicylindrical vapour bubble model, which is in good agreement with the experimental findings. Finally, this allows us to predict the relevance of our findings for the transport of cryogenic liquids in huge overseas carriers where wave impact due to sloshing is the dominant cause of hydrodynamic load of containment systems in cargo tanks.

Three-dimensional wake forcing is applied to a profiled blunt trailing edge body from synthetic jet arrays distributed symmetrically on both sides of the body. The effect on the wake is experimentally studied at Reynolds numbers based on body thickness, $d$, of $2500 \leqslant Re_d=u_\infty d/\nu \leqslant 5000$ in the turbulent wake regime. The exits of the synthetic jets are rectangular slots and are oriented spanwise to the cross-flow with a uniform spacing of $2.4d$. The forcing causes spanwise variations in the separated shear layers, leading to the von Kármán vortices tilting and forming coherent streamwise vortex loops. This reorientation of the wake vorticity is associated with the attenuation of the vortex street and drag reduction, consistent with previous studies of spanwise perturbations to wakes. The effect of forcing amplitude on the drag and wake structure is examined. It is found that the mean shedding frequency is constant across the span in all cases, indicating that the forced wake has a periodic organised structure. The greatest drag reduction of approximately 25 % is achieved when the vortical structures emitted by the jets penetrate up the edges of the boundary layers of the body, which occurs at velocity ratios (defined from the mean jet exit velocity during expulsion) of about 3 when $Re_d=2500$ and about 2 when $Re_d=5000$. This study presents evidence that the forcing effectiveness is maximised when the vortex street is most tilted into the streamwise direction.

Surface gravity waves induce a drift on objects floating on the water’s surface. This study presents laboratory experiments investigating the drift of large two-dimensional floating objects on deep-water, unidirectional, regular waves, with wave steepness ranging from 0.04 to 0.31 (0.04 $\lt k{a_w}\lt$ 0.31, where $k$ is the wavenumber and $a_w$ the wave amplitude). The objects were carefully designed to have a rectangular cross-section with a constant aspect ratio; their size varied from 2.6 $\%$ to 27 $\%$ of the incident wavelength. We observed Lagrangian behaviour for small objects. Small and large objects exhibited fundamentally different drift behaviour at high compared with low wave steepness, with a regime shift observed at a certain size and wave steepness. The scaling of object drift with steepness depends on the relative size of the object. For small objects, drift scales with steepness squared, whereas drift becomes a linear function of steepness as the object size increases. For objects that are relatively large but smaller than 13–16% of a wavelength (low to high steepness), we provide experimental evidence supporting the mechanisms of drift enhancement recently identified by Xiao et al. (J. Fluid Mech., vol. 980, 2024, p. A27) and termed the ‘diffraction-modified Stokes drift’. This enhanced drift behaviour, compared with the theoretical Stokes drift for infinitely small fluid parcels, is attributed to changes in the objects’ oscillatory motion and local wave amplitude distribution (standing wave pattern) due to the presence of the object. In the case of larger objects, similar to Harms (J. Waterw. Port Coast. Ocean Eng., vol. 113(6), 1987, pp. 606–622), we relate the critical size at which drift is maximised to their vertical bobbing motion. We determine the domain of validity for both Stokes drift and the diffraction-modified Stokes drift model of Xiao et al. (J. Fluid Mech., vol. 980, 2024, A27) in terms of relative size and wave steepness and propose an empirical parametrisation based on our experimental data.

An important parameter characterising the synchronisation of turbulent flows is the threshold coupling wavenumber. This study investigates the relationship between the threshold coupling wavenumber and the leading Lyapunov vector using large eddy simulations and the SABRA model. Various subgrid-scale stress models, Reynolds numbers and different coupling methods are examined. A new scaling relation is identified for the leading Lyapunov exponents in large eddy simulations, showing that they approximate those of filtered direct numerical simulations. This interpretation provides a physical basis for results related to the Lyapunov exponents of large eddy simulations, including those related to synchronisation. Synchronisation experiments show that the peak wavenumber of the energy spectrum of the leading Lyapunov vector coincides with the threshold coupling wavenumber, in large eddy simulations of box turbulence with standard Smagorinsky or dynamic mixed models as well as in the SABRA model, replicating results from direct numerical simulations of box turbulence. Although the dynamic Smagorinsky model exhibits different behaviour, the totality of the results suggests that the relationship is an intrinsic property of a certain class of chaotic systems. We also confirm that conditional Lyapunov exponents characterise the synchronisation process in indirectly coupled systems as they do in directly coupled ones, with their values insensitive to the nature of the master flow. These findings advance the understanding of the role of the Lyapunov vector in the synchronisation of turbulence.

The linear stability of miscible displacement for radial source flow at infinite Péclet number in a Hele-Shaw cell is calculated theoretically. The axisymmetric self-similar flow is shown to be unstable to viscous fingering if the viscosity ratio $m$ between ambient and injected fluids exceeds $3/2$, and to be stable if $m\lt {3/2}$. If $1\lt m\lt {3/2}$, then small disturbances decay at rates between $t^{-3/4}$ and $t^{-1}$ (the exact range depending on $m$) relative to the $t^{1/2}$ radius of the axisymmetric base-state similarity solution; if $m\lt 1$, then they decay faster than $t^{-1}$. Asymptotic analysis confirms these results and gives physical insight into various features of the numerically determined relationship between the growth rate and the azimuthal wavenumber and viscosity ratio.

The evolution of turbulent spots in a flat plate boundary layer is examined using time-resolved tomographic particle image velocimetry (Tomo-PIV) experiments and direct numerical simulation (DNS). The characteristics of flow structures are examined using timelines and material surfaces. Both the numerical and experimental results reveal a notable behaviour in the developmental process of turbulent spots: the development of low-speed streaks at the spanwise edges of turbulent spots, followed by the subsequent formation of hairpin vortices. The behaviour of these low-speed streaks is further investigated using timelines and material surfaces generated for a series of regions and development times. The results indicate that these low-speed streaks exhibit characteristic wave behaviour. The low-speed streaks are observed to lift up as three-dimensional (3-D) waves, with high-shear layers forming at the interface of these waves. These induced high-shear layers become unstable and evolve into vortices, which contribute to the expansion of the turbulent spot. These findings show the significant role of 3-D waves in the development of turbulent spots, supporting the hypothesis that 3-D waves serve as initiators of vortices at the bounding surface of a turbulent spot.

Using clean numerical simulation (CNS) in which artificial numerical noise is negligible over a finite, sufficiently long interval of time, we provide evidence, for the first time, that artificial numerical noise in direct numerical simulation (DNS) of turbulence is approximately equivalent to thermal fluctuation and/or stochastic environmental noise. This confers physical significance on the artificial numerical noise of DNS of the Navier–Stokes equations. As a result, DNS on a fine mesh should correspond to turbulence under small internal/external physical disturbance, whereas DNS on a sparse mesh corresponds to turbulent flow under large physical disturbance. The key point is that all of them have physical meanings and so are correct in terms of their deterministic physics, even if their statistics are quite different. This is illustrated herein. Our paper provides a positive viewpoint regarding the presence of artificial numerical noise in DNS.

We study the onset of spontaneous dynamics in the follower force model of an active filament, wherein a slender elastic filament in a viscous liquid is clamped normal to a wall at one end and subjected to a tangential compressive force at the other. Clarke et al. (Phys. Rev. Fluids, vol. 9, 2024, 073101) recently conducted a thorough investigation of this model using methods of computational dynamical systems; inter alia, they showed that the filament first loses stability via a supercritical double-Hopf bifurcation, with periodic ‘planar-beating’ states (unstable) and ‘whirling’ states (stable) simultaneously emerging at the critical follower-force value. We complement their numerical study by carrying out a weakly nonlinear analysis close to this unconventional bifurcation, using the method of multiple scales. The main outcome is an ‘amplitude equation’ governing the slow modulation of small-magnitude oscillations of the filament in that regime. Analysis of this reduced-order model provides insights into the onset of spontaneous dynamics, including the creation of the nonlinear whirling states from particular superpositions of linear planar-beating modes as well as the selection of whirling over planar beating in three-dimensional scenarios.

We perform simulations of a two-fluid–structure interaction problem involving liquid–gas flow past a fully submerged stationary circular cylinder. Interactions between the liquid–gas interface with finite surface tension and flow disturbances arising from the cylinder induce a variety of interfacial phenomena and wake structures. We map different interface regimes in a parameter space defined by the Bond number $Bo \in [100, 5000]$ and the submergence depth $h/D \in [1, 2.5]$ of the cylinder while keeping the Reynolds (Re) and Weber (We) numbers fixed at 150 and 1000, respectively. The emerging interface features are classified into three distinct regimes: interfacial waves generated by Strouhal vortices, the entrainment of multi-scale gas bubbles and the reduced deformation state. In the interfacial wave regime, we demonstrate that the frequency of transverse interface fluctuations at a specific streamwise location is identical to the vortex shedding frequency. Additionally, the wavelength of interfacial waves is determined by the size of vortex pairs consisting of alternating Strouhal vortices. In the gas entrainment regime at $ Bo = 1000$, our bubble-size distributions reveal that the entrained bubbles have sizes ranging from one to two orders of magnitude smaller than the cylinder. These multi-scale bubbles are formed primarily through plunging and surfing breakers at $h/D = 2.5$. In contrast, at $h/D = 1$, smaller bubbles initially emerge from the breakup of a gas finger. Over time, some of these bubbles grow in size through coalescence cascades. The influence of $ Re \in [50, 150]$ and $ We \in [700, 1100]$ on gas entrainment is quantified in terms of mean bubble size and count. Lastly, we demonstrate how the deformability of the liquid–gas interface drives the hydrodynamic lift force acting on the cylinder. The net downward lift materializes only in the gas entrainment and reduced deformation regimes due to the broken symmetry of the front stagnation point. While our study focuses on two-dimensional simulations, we also provide insights into the three-dimensional gas entrainment mechanism for one of the extreme cases at $h/D = 1$.

We conduct direct numerical simulations (DNS) to investigate the attenuation of turbulence in a periodic cube due to the addition of prolate spheroidal solid particles. Even with a dilute volume fraction of $O(10^{-2})$, particles can drastically attenuate the turbulence. Our DNS show that the turbulent kinetic energy reduces more significantly when the particles’ Stokes number is larger, size is smaller or aspect ratio is larger. We can explain these results based on the formula proposed by Oka and Goto (2022 J. Fluid Mech. 949, A45), which relates the turbulence attenuation rate to the energy dissipation rate $\epsilon _p$ around particles. More precisely, under the condition that the volume fraction of particles is fixed, $\epsilon _p$ is larger when the Stokes number and, therefore, the relative velocity between fluid and particles are larger, the particle size is smaller or the aspect ratio is larger. These results also imply that the rotation of the anisotropic particles plays only a limited role in the attenuation of turbulence when the Stokes number of particles is sufficiently large, because the main cause of the attenuation is the relative translational velocity between fluid and particles.