For turbulent bubbly flows, multi-phase simulations resolving both the liquid and bubbles are prohibitively expensive in the context of different natural phenomena. One example is breaking waves, where bubbles strongly influence wave impact loads, acoustic emissions and atmospheric-ocean transfer, but detailed simulations in all but the simplest settings are infeasible. An alternative approach is to resolve only large scales, and model small-scale bubbles adopting sub-resolution closures. Here, we introduce a large eddy simulation smoothed particle hydrodynamics (SPH) scheme for simulations of bubbly flows. The continuous liquid phase is resolved with a semi-implicit isothermally compressible SPH framework. This is coupled with a discrete Lagrangian bubble model. Bubbles and liquid interact via exchanges of volume and momentum, through turbulent closures, bubble breakup and entrainment, and free-surface interaction models. By representing bubbles as individual particles, they can be tracked over their lifetimes, allowing closure models for sub-resolution fluctuations, bubble deformation, breakup and free-surface interaction in integral form, accounting for the finite time scales over which these events occur. We investigate two flows: bubble plumes and breaking waves, and find close quantitative agreement with published experimental and numerical data. In particular, for plunging breaking waves, our framework accurately predicts the Hinze scale, bubble size distribution, and growth rate of the entrained bubble population. This is the first coupling of an SPH framework with a discrete bubble model, with potential for cost-effective simulations of wave–structure interactions and more accurate predictions of wave impact loads.

Microfluidic flow focusing is a versatile method for the production of monodisperse microbubbles for biomedical applications involving ultrasound. Existing studies propose several theoretical models to predict bubble size and production rate as a function of the liquid and gas flow rate. Yet, they typically do not include physical fluid parameters such as density, viscosity and surface tension. Here, we present an exhaustive experimental and numerical investigation of the influence of physical properties of the gas and liquid, and of the channel geometry on bubble size and production rate. We find a particularly strong effect of (i) gas density on the production rate and (ii) liquid viscosity on the bubble size. We further discuss our findings within the context of existing theoretical models to reflect on gaps in our current understanding of the fluid mechanics of bubble formation by flow focusing.

Aerobreakup of drops is a fundamental two-phase flow problem that is essential to many spray applications. A parametric numerical study was performed by varying the gas stream velocity, focusing on the regime of moderate Weber numbers, in which the drop deforms to a forward bag. When the bag is unstable, it inflates and disintegrates into small droplets. Detailed numerical simulations were conducted using the volume-of-fluid method on an adaptive octree mesh to investigate the aerobreakup dynamics. Grid-refinement studies show that converged three-dimensional simulation results for drop deformation and bag formation are achieved by the refinement level equivalent to 512 cells across the initial drop diameter. To resolve the thin liquid sheet when the bag inflates, the mesh is refined further to 2048 cells across the initial drop diameter. The simulation results for the drop length and radius were validated against previous experiments, and good agreement was achieved. The high-resolution results of drop morphological evolution were used to identify the different phases in the aerobreakup process, and to characterize the distinct flow features and dominant mechanisms in each phase. In the early time, the drop deformation and velocity are independent of the Weber number, and a new internal-flow deformation model, which respects this asymptotic limit, has been developed. The pressure and velocity fields around the drop were shown to better understand the internal flow and interfacial instability that dictate the drop deformation. Finally, the impact of drop deformation on the drop dynamics was discussed.

An experiment on convective flows induced by the dielectrophoretic force was performed under the microgravity condition provided during a sounding rocket flight. The dielectrophoretic force possesses a non-conservative term that can be seen as resulting from an electric gravity. That gravity can be responsible for an electric Rayleigh–Bénard convection between a hot inner cylinder and a cold outer cylinder when an electric field is applied in the radial direction. Four cells with independent temperature and electric field controls allowed the investigation of eight different values of the electric Rayleigh number relatively close to the onset of the thermo-electric instability. A linear stability analysis is performed to predict the stability threshold and the evolution of the growth rate of the instability. The three-dimensional structure of the flow is captured by simultaneous particle image velocimetry and by shadowgraphy. The amplitude of the instability modes and the time evolution of the flow is analysed, and various methods are proposed to extrapolate the experimental critical value of the electric Rayleigh number for the onset of convection. The measured critical electric Rayleigh number is in agreement with the prediction of the linear stability theory. The comparison of the new experimental results with previous ones from parabolic flight campaigns highlights the importance of long-term microgravity for the achievement of thermal convection at low values of the control parameters.

The effects of surfactants on a mechanically generated plunging breaker are studied experimentally in a laboratory wave tank. Waves are generated using a dispersively focused wave packet with a characteristic wavelength of $\lambda _0 = 1.18$ m. Experiments are performed with two sets of surfactant solutions. In the first set, increasing amounts of the soluble surfactant Triton X-100 are mixed into the tank water, while in the second set filtered tap water is left undisturbed in the tank for wait times ranging from 15 min to 21 h. Increasing Triton X-100 concentrations and longer wait times lead to surfactant-induced changes in the dynamic properties of the free surface in the tank. It is found that low surface concentrations of surfactants can dramatically change the wave breaking process by changing the shape of the jet and breaking up the entrained air cavity at the time of jet impact. Direct numerical simulations (DNS) of plunging breakers with constant surface tension are used to show that there is significant compression of the free surface near the plunging jet tip and dilatation elsewhere. To explore the effect of this compression/dilatation, the surface tension isotherm is measured in all experimental cases. The effects of surfactants on the plunging jet are shown to be primarily controlled by the surface tension gradient ($\Delta \mathcal {E}$) while the ambient surface tension of the undisturbed wave tank ($\sigma _0$) plays a secondary role.

Buoyancy-driven turbulent convection leads to a fully compressible flow with a prominent top-down asymmetry of first- and second-order statistics when the adiabatic equilibrium profiles of temperature, density and pressure change very strongly across the convection layer. The growth of this asymmetry and the formation of an increasingly thicker stabilized sublayer with a slightly negative mean convective heat flux $J_c(z)$ at the top of the convection zone is reported here by a series of highly resolved three-dimensional direct numerical simulations beyond the Oberbeck–Boussinesq and anelastic limits for dimensionless dissipation numbers, $0.1 \le D\le 0.8$, at fixed Rayleigh number $Ra=10^6$ and superadiabaticity $\epsilon =0.1$. The highly stratified compressible convection regime appears for $D > D_{crit}\approx 0.65$, when density fluctuations collapse to those of pressure; it is characterized by an up to nearly 50 % reduced global turbulent heat transfer and a sparse network of focused thin and sheet-like thermal plumes falling through the top sublayer deep into the bulk.

In impermeable media, a hydraulic fracture can continue expanding even without additional fluid injection if its volume exceeds the limiting volume of a hydrostatically loaded radial fracture. This limit depends on the mechanical properties of the surrounding solid and the density contrast between the fluid and the solid. We show that two dimensionless numbers characterize self-sustained fracture growth. The first is a buoyancy factor that compares the total released volume to the volume of a hydrostatically loaded radial fracture to determine whether buoyant growth occurs. The second number is the dimensionless viscosity of a radial fracture when buoyant effects become of order one. Notably, this dimensionless viscosity depends on the rate at which the fluid volume is released, indicating that both the total volume and release history impact self-sustained buoyant growth. We identify six well-defined propagation histories based on these two dimensionless numbers. Their growth evolves between distinct limiting regimes of radial and buoyant propagation, resulting in different fracture shapes. Notably, our findings reveal two growth rates depending on the dominant energy dissipation mechanism (viscous flow versus fracture creation) in the fracture head. For finite values of material toughness, the toughness-dominated limit represents a late-time solution for all fractures in growth rate and head shape (possibly reached only at a very late time). The viscosity-dominated limit can appear at intermediate times. Our three-dimensional simulations confirm the predicted scalings. This contribution highlights the importance of the entire propagation and release history for accurate analysis of buoyant hydraulic fractures.

We investigate the coupling effect of buoyancy and shear based on an annular centrifugal Rayleigh–Bénard convection (ACRBC) system in which two cylinders rotate with an angular velocity difference. Direct numerical simulations are performed in a Rayleigh number range $10^6\leq Ra\leq 10^8$, at fixed Prandtl number $Pr=4.3$, inverse Rossby number $Ro^{-1}=20$, and radius ratio $\eta =0.5$. The shear, represented by the non-dimensional rotational speed difference $\varOmega$, varies from $0$ to $10$, corresponding to an ACRBC without shear and a radially heated Taylor–Couette flow with only the inner cylinder rotating, respectively. A stable regime is found in the middle part of the interval for $\varOmega$, and divides the whole parameter space into three regimes: buoyancy-dominated, stable and shear-dominated. Clear boundaries between the regimes are given by linear stability analysis, meaning the marginal state of the flow. In the buoyancy-dominated regime, the flow is a quasi-two-dimensional flow on the $r\varphi$ plane; as shear increases, both the growth rate of instability and the heat transfer are depressed. In the shear-dominated regime, the flow is mainly on the $rz$ plane. The shear is so strong that the temperature acts as a passive scalar, and the heat transfer is greatly enhanced. The study shows that shear can stabilize buoyancy-driven convection, makes a detailed analysis of the flow characteristics in different regimes, and reveals the complex coupling mechanism of shear and buoyancy, which may have implications for fundamental studies and industrial designs.

The coupled effects of thermodynamic and hydrodynamic instabilities are studied during viscous fingering (VF). We introduced a modified Cahn–Hilliard phase-field model in conjunction with the Korteweg force in the classical VF model and derived consistent governing equations. The free energy of the partially miscible system is described using a modified Flory–Huggins model, which allows us to investigate the temporal evolution of spatial inhomogeneities. The mass flux in the Cahn–Hilliard equations is modified according to modern diffusion theory. The governing equations have been solved through an in-house model implementation using the COMSOL multiphysics software. We successfully demonstrated the transition from the finger-like structures to the droplet formation during spinodal decomposition as demonstrated experimentally in the literature. Our results are also in agreement with earlier numerical results obtained using a classical Landau type mixing energy. We further systematically studied the effects of the Margules parameter (interaction parameter) and the gradient parameter, which is associated to the thermodynamic length scale and capillary number on the VF. Aysmmetric features of the binary mixture are also investigated showing a stronger thermodynamic effect on the system with increasing phase separation and, hence, droplet formation.

An important problem in passive scalar transport is to parametrize the effect of a fluctuating component of the flow, in order to overcome a limited resolution. A local effective diffusivity is one such parametrization, and over the years there have been many different suggestions for ‘closures’ that relate the advective flux to gradients of the mean concentration. Souza et al. (J. Fluid Mech., 2023, in press) introduce a stochastic framework where the local effective diffusivity is replaced by an exact effective diffusivity operator. By computing this operator for various examples, they quantify deviations from the local approximation, which can suggest areas of improvement and novel closure models.

Theory and computations are applied to assess the hydrodynamic permeability of cavity- doped hydrogels, central to a variety range of contemporary technological applications. Direct volume-averaging is undertaken in a two-dimensional, Brinkman-hydrodynamic context to test an ensemble-averaging methodology recently proposed for the ion permeability of such media. In two dimensions, the ensemble-averaging integral furnishes a pre-factor $2$ linking the pressure dipole strength of a single inclusion in an unbounded continuum to the effective hydrodynamic permeability of a composite with small inclusion area fraction. The factor is verified by direct computations for dilute simple-square arrays of (cylindrical) inclusions. At area fractions up to the close-packing limit, computations address the hydrodynamic interactions. The theory is shown to accurately predict the effective hydrodynamic permeability of physically relevant composites (Brinkman length of the continuous phase $\ell$ smaller than the inclusion radius $a$) for area fractions $\phi \lesssim {\rm \pi}/ 9 \approx 0.3$. Computations for random ensembles demonstrate that the dilute theory may be extended to higher area fractions by drawing on Rayleigh's self-consistent approximation when the continuous-phase permeability places the continuous-phase flow well into the Darcy regime ($a / \ell \gtrsim 10$). Computations also demonstrate, similarly to Rayleigh theories for scalar diffusion, that microstructural order has a very weak influence on the effective permeability when $\phi \lesssim {\rm \pi}/ 9$ with $\ell / a \lesssim 1$ (Darcy hydrodynamic interactions). Finally, a cursory examination is undertaken of the fluid velocity and its fluctuations arising from shear-viscosity heterogeneity in media with perfectly uniform permeability $\ell ^2$.

Pressure-gradient-driven flows through sinusoidal channels have been studied. The analysis was carried out up to the formation of secondary nonlinear states and spanned a range of low and moderate Reynolds numbers. Direct numerical simulations were used to identify and determine the properties of steady as well as non-stationary, two-dimensional (2-D) and three-dimensional secondary flows. Our results indicate the existence of several distinct solution types. Two-dimensional, stationary flows with periodicity determined by the corrugation represent the first type. The second type is associated with the appearance of 2-D oscillatory flows arising from the onset of unstable travelling waves. Such oscillatory solutions are generally out of phase with the wall corrugation but could be in phase in special cases determined by the ratio of the critical disturbance wavelength and the channel corrugation wavelength. Consequently, several distinct types of time-dependent solutions are possible. The third type of solution results from the centrifugal effect caused by wall curvature and leads to three-dimensionalization of the flow through the onset of stationary streamwise vortices. Finally, various states resulting from the interaction of different solution types are possible. We examine those states and present a bifurcation diagram illustrating the formation of some of them. The results presented in this paper might help with the development of small-scale flow measurement and detection devices operating at low and moderate Reynolds numbers, as well as in the use of wall topographies for the intensification of mixing in flows with moderate, subturbulent Reynolds numbers.