We explore predictions of two models of one-dimensional capillary rise in rigid and partially saturated porous media. One is an existing one from the literature and the second is a free-boundary model based on Richards’ equation with two moving boundaries of the evolving partially saturated region. Both models involve the specification of saturation-dependent functions for local capillary pressure and permeability and connect to classical models for saturated porous media. Existing capillary-rise experiments show two notable regimes: (i) an early-time regime typically well-described by classical capillary-rise theory in a fully saturated porous media, and (ii) a long-time regime that has anomalous dynamics in which the capillary-rise height may scale with a non-classical power law in time or have more complicated dynamics. We demonstrate that the predictions of both models compare well with experimental capillary-rise data over early- and long-time regimes gathered from three independent studies in the literature. The model predictions also shed light on recent scaling laws that relate the capillary pressure and permeability of the partially saturated media to the capillary-rise height. We use these models to probe computationally observed permeability relationships to capillary-rise height. We demonstrate that a recently proposed permeability scaling for the anomalous capillary-rise regime is indeed realized and is particularly apparent in the lower portion of the partially saturated media. For our free-boundary model we also compute capillary pressure measures and show that these reveal the linear relation between the capillary pressure and capillary-rise height expected for a capillarity–gravity balance in the upper portion of the partially saturated porous media.

While the problem governing Stokes flow about a single particle that is subject to an external force is ill posed in two dimensions (the ‘Stokes paradox’), the related problem of two mutually repellent particles is well posed. Motivated by self-assembly phenomena in thin viscous membranes, we consider this problem in the limit of remote particles. Such limits are typically handled in the literature using reflection techniques, which provide successive approximations to the mutual hydrodynamic interactions. Since their starting point is a single particle in an unbounded fluid domain, these techniques are futile in the present two-dimensional problem. We show how this apparent contradiction is resolved via use of singular perturbations. We obtain a two-term approximation for the velocity acquired by circular disks, considering both rigid and free particle surfaces. We also illustrate our perturbation scheme for elliptic disks, deriving a renormalised single-particle velocity. The utility of our asymptotic scheme is illustrated in the general problem of hydrodynamic interaction between a cluster of remote disks.

Coherent small-amplitude unsteadiness of the shock wave and the separation region over a canonical double cone flow, termed in literature as oscillation-type unsteadiness, is experimentally studied at Mach 6. The double cone model is defined by three non-dimensional geometric parameters: fore- and aft-cone angles ($\theta _1$ and $\theta _2$), and ratio of the conical slant lengths ($\varLambda$). Previous studies of oscillations have been qualitative in nature, and mostly restricted to a special case of the cone model with fixed $\theta _1 = 0^\circ$ and $\theta _2 = 90^\circ$ (referred to as the spike-cylinder model), where $\varLambda$ becomes the sole governing parameter. In the present effort we investigate the self-sustained flow oscillations in the $\theta _1$-$\varLambda$ parameter space for fixed $\theta _2 = 90^\circ$ using high-speed schlieren visualisation. The experiments reveal two distinct subtypes of oscillations, characterised by the motion (or lack thereof) of the separation point on the fore-cone surface. The global time scale associated with flow oscillation is extracted using spectral proper orthogonal decomposition. The non-dimensional frequency (Strouhal number) of oscillation is seen to exhibit distinct scaling for the two oscillation subtypes. The relationship observed between the local flow properties, instability of the shear layer, and geometric constraints on the flow suggests that an aeroacoustic feedback mechanism sustains the oscillations. Based on this understanding, a simple model with no empiricism is developed for the Strouhal number. The model predictions are found to match well with experimental measurements. The model provides helpful physical insight into the nature of the self-sustained flow oscillations over a double cone at high speeds.

A hybrid lattice Boltzmann and finite difference method is applied to study the influence of soluble surfactants on droplet deformation and breakup in simple shear flow. First, the influence of bulk surfactant parameters on droplet deformation in two-dimensional shear flow is investigated, and the surfactant solubility is found to influence droplet deformation by changing average interface surfactant concentration and non-uniform effects induced by non-uniform interfacial tension and Marangoni forces. In addition, the droplet deformation first increases and then decreases with Biot number, increases significantly with adsorption number $k$ and decreases with Péclet number or adsorption depth; and among the parameters, $k$ is the most influential one. Then, we consider three-dimensional shear flow and investigate the roles of surfactants on droplet deformation and breakup for different capillary numbers and viscosity ratios. Results show that in the soluble case with $k=0.429$, the droplet exhibits nearly the same deformation as in the insoluble case due to the balance between surfactant adsorption and desorption; upon increasing $k$ from 0.429 to 1, the average interface surfactant concentration is greatly enhanced, leading to significant increase in droplet deformation. The critical capillary number of droplet breakup $C{{a}_{cr}}$ is identified for varying viscosity ratios in clean, insoluble and soluble ($k=0.429$ and 1) systems. As the viscosity ratio increases, $C{{a}_{cr}}$ first decreases and then increases rapidly in all systems. The addition of surfactants always favours droplet breakup, and increasing solubility or $k$ could further reduce $C{{a}_{cr}}$ by increasing average interface surfactant concentration and local surfactant concentration near the neck during the necking stage.

In this study, the asymptotic solutions of the pressure variations induced by two trains passing each other in a tunnel are theoretically investigated. The one-dimensional inviscid compressible airflow is analysed, and two methods to obtain numerically exact solutions and $M_{H}$ expansion formulas for approximate equations are presented, where $M_{H}$ is the Mach number of the high-speed train. The pressure coefficient, corresponding to the maximum value of the magnitude of the pressure, is expressed as $|c_{p}|_{max}=|c_{p,min}|=[({R}/({1-R}))$ $(1+\alpha )^{2}+({R(1-R)}/{(1-2R)^{2}})(1-\alpha )^{2}]+O[M_{H}]$, where $c_{p,min}<0$, $\alpha =U_{L}/U_{H}$ and $U_{L}$ and $U_{H}$ denote the speeds of the low- and high-speed trains, respectively, and $R$ is the cross-sectional area ratio of the train to the tunnel. The theoretical results indicate the dependence of the speeds of the two trains on the pressure distribution and that the maximum magnitude of the asymptotic pressure for a fixed value of $M_{H}$ is obtained for $\alpha =1$ and $\alpha =0$ when $R< R_{c}$ and $R>R_{c}$, respectively, where $R_{c}$ denotes the critical blockage ratio. Because the airflow along the side of the low-speed train, induced by the low-speed train, is along the running direction of the high-speed train and reduces the relative velocity of the high-speed train as the two trains pass each other, $|c_{p}|_{max}$ for $\alpha =0$ is larger than $|c_{p}|_{max}$ for $\alpha =1$ when $R>R_{c}$. It is theoretically demonstrated that, as conventional high-speed railway systems satisfy $R< R_{c}$, a conservative pressure estimation can be established assuming $\alpha =1$.

The non-Oberbeck–Boussinesq effects on the stability of a vertical natural convection boundary layer are investigated using the linearised disturbance equations for air flows up to a temperature difference of $\Delta T=100\,{\rm K}$. Based on the linear stability results, the neutral curve is shown to be sensitive to the choice of reference temperature. When evaluated using the film temperature $T_f$, a lower film Grashof number is required to trigger the linear instability for larger $\Delta T$. The relative contributions of shear and buoyant production to the perturbation kinetic energy budget reveals that the marginally unstable modes are amplified based on different mechanisms: for lower wavenumbers at relatively small Grashof number, the instability is driven by buoyancy; whereas for higher wavenumbers and larger Grashof number, the flow becomes unstable due to a shear instability. The use of reference temperature is found to scale the shear- and buoyant-driven instabilities differently so that no single reference temperature definition would collapse the neutral curves. The linear stability result further demonstrates that at a given Grashof number a higher temperature difference would give a larger amplification rate of the perturbation, which then leads to an earlier onset of the nonlinearities when evaluated at $T_f$. Finally, by comparing the amplification rates obtained from direct numerical simulation and the linear stability results, the extent of the linear regime is determined for $\Delta T = 100\,{\rm K}$.

In this study we propose a novel data-driven reduced-order model for complex dynamics, including nonlinear, multi-attractor, multi-frequency and multiscale behaviours. The starting point is a fully automatable cluster-based network model (CNM) (Li et al., J. Fluid Mech., vol. 906, 2021, A21) that kinematically coarse grains the state with clusters and dynamically predicts the transitions in a network model. In the proposed dynamics-augmented CNM (dCNM) the prediction error is reduced with trajectory-based clustering using the same number of centroids. The dCNM is first exemplified for the Lorenz system and then demonstrated for the three-dimensional sphere wake featuring periodic, quasi-periodic and chaotic flow regimes. For both plants, the dCNM significantly outperforms the CNM in resolving the multi-frequency and multiscale dynamics. This increased prediction accuracy is obtained by stratification of the state space aligned with the direction of the trajectories. Thus, the dCNM has numerous potential applications to a large spectrum of shear flows, even for complex dynamics.

In 2017, Brosseau & Vlahovska (Phys. Rev. Lett, vol. 119, no. 3, 2017, p. 034501) found that, in a strong electric field, a weakly conductive, low-viscosity droplet immersed in a highly conductive, high-viscosity medium formed a lens shape, and liquid rings continuously detached from its equatorial plane and subsequently broke up into satellite droplets. This fascinating multiphase electrohydrodynamic (EHD) phenomenon is known as droplet equatorial streaming. In this paper, based on the unified lattice Boltzmann method framework proposed by Luo et al. (Phil. Trans. R. Soc. A Math. Phys. Engng Sci, vol. 379, no. 2208, 2021, p. 20200397), a novel lattice Boltzmann (LB) model is constructed for multiphase EHD by coupling the Allen–Cahn type of multiphase LB model and two new LB equations to solve the Poisson equation of the electric field and the conservation equation of the surface charge. Using the proposed LB model, we successfully reproduced, for the first time, the complete process of droplet equatorial streaming, including the continuous ejection and breakup of liquid rings on the equatorial plane. In addition, it is found that, under conditions of high electric field strength or significant electrical conductivity contrast, droplets exhibit fingering equatorial streaming that was unknown before. A power-law relationship is discovered for droplet total charge evolution and a theoretical model is then proposed to describe the droplet radius and height over time. The breakup of liquid rings is found to be dominated by capillary instability, while the breakup of liquid fingers is governed by the end-pinching mechanism. Finally, a phase diagram is constructed for fingering equatorial streaming and ring equatorial streaming, and a criterion equation is established for the phase boundary.

The evolution of the water-entry cavity affects the impact load and the motion of the body. This paper adopts the Eulerian finite element method for multiphase flow for simulations of the high-speed water-entry process. The accuracy and convergence of the numerical method are verified by comparing it with the experimental data and the results of the transient cavity dynamics theory. Based on the results, the representative characteristics of the cavity are discussed from the perspective of the cavity cross-section. It is found that the asymmetry of the cavity expansion and contraction durations is related to the motion of the free surface and the closure of the cavity. The uplift of the free surface suppresses cavity expansion, while the jet generated from free surface closure accelerates cavity contraction. The duration of the contraction of the cavity near the free surface is shorter than the expansion duration due to the change in the velocity distribution caused by the free surface motion. The necking phenomenon during deep closure leads to an increase in the internal pressure of the cavity, prolonging cavity contraction near the deep closure area. This work provides new insights into the cavity dynamics in high-speed water entry.

Plane turbulent wall jets are traditionally considered to be composed of a turbulent boundary layer (TBL) topped by a half-free jet. However, certain peculiar features, such as counter-gradient momentum flux occurring below velocity maximum in experiments and numerical simulations, suggest a different structure of turbulence therein. Here, we hypothesize that turbulence in wall jets has two distinct structural modes, wall mode scaling on wall variables and free-jet mode scaling on jet variables. To investigate this hypothesis, experimental data from our wall jet facility are acquired using single hot-wire anemometry and two-dimensional particle image velocimetry at three nozzle Reynolds numbers 10 244, 15 742 and 21 228. Particle image velocimetry measurements with four side-by-side cameras capture the longest field of view studied so far in wall jets. Direct spatial spectra of these fields reveal modal spectral contributions to variances of velocity fluctuations, Reynolds shear stress, shear force, turbulence production, velocity fluctuation triple products and turbulent transport. The free-jet mode has wavelengths scaling on the jet length scale ${z_{T}}$, and contains two dominant submodes with wavelengths $5{z_{T}}$ and $2.5{z_{T}}$. The region of flow above the velocity maximum shows the presence of the outer jet mode whereas the region below it shows robust bimodal behaviour attributed to both wall and inner jet modes. Counter-gradient momentum flux is effected by the outer jet mode intruding into the region below velocity maximum. These findings support the hypothesis of wall and free-jet structural modes, and indicate that the region below velocity maximum could be much complex than a conventional TBL.

Vortex-induced vibrations and galloping of an elastically mounted square cylinder are investigated for cylinder mass ratio m* = 2–50, damping ratio ζ = 0–1.0, mass-damping ratio m*ζ = 0–50 and flow reduced velocity Ur = 1–80. We home in on the effects of m*, ζ, m*ζ, $({m^\ast } + m_{a\textrm{0}}^\ast )\zeta$ and $({m^\ast } + m_{ae}^\ast )\zeta$ on the critical reduced velocity Urc marking the onset of galloping, where $m_{a\textrm{0}}^\ast $ is the quiescent-fluid added mass ratio and $m_{ae}^\ast $ is the effective added mass ratio. Vibration responses, forces, vibration frequencies and added mass ratios are studied and discussed. The different branches of vortex-induced vibrations have different dependencies of $m_{ae}^\ast $ on Ur. The $m_{ae}^\ast $ in the initial branch is positive and drops rapidly with Ur, but that in the lower branch is negative and declines gently. In the galloping regime, $m_{ae}^\ast $ jumps from negative to positive at the onset of galloping, declining slightly with increasing Ur. Our results and prediction equations show that when ζ = 0, Urc is independent of m* for m* ≥ 5, albeit slightly higher for m* = 3. The latter is ascribed to mode competition. When ζ > 0, Urc linearly increases with increasing ζ. Detailed analysis substantiates that m*ζ or $({m^\ast } + m_{a\textrm{0}}^\ast )\zeta$ does not serve as the unique criterion to predict the galloping occurrence. Here, we propose a new combined mass-damping parameter $({m^\ast } + m_{ae}^\ast )\zeta$ in the relationship between galloping onsets and structural properties, which successfully scales all data of Urc at different m* and ζ values.

Dense mixtures of particles of varying size tend to segregate based on size during flow. Granular size segregation impacts many industrial and geophysical processes, but the development of coupled, continuum models capable of predicting the evolution of segregation dynamics and flow fields in dense granular media across different geometries remains a challenge. One reason is because size segregation stems from two driving forces: pressure gradients and shear-strain-rate gradients. Another reason is the challenge of integrating segregation models with rheological constitutive equations for dense granular flow. In this paper we develop a continuum model that accounts for pressure-gradient-driven and shear-strain-rate-gradient-driven segregation, coupled to rheological modelling of a dense granular medium across the quasi-static and dense inertial flow regimes. To calibrate and test the continuum model, we perform discrete element method (DEM) simulations of dense flow of bidisperse granular systems in two flow geometries in which both segregation driving forces are present: inclined plane flow and planar shear flow with gravity. Steady-state DEM data from inclined plane flow is used to determine the dimensionless material parameters in the pressure-gradient-driven segregation model for both spheres and disks. Then, predictions of the continuum model are tested against DEM data across different cases of inclined plane flow and planar shear flow with gravity, while varying parameters such as the size of the flow geometry, the flow speed and the initial conditions. We find that it is crucial to account for both driving forces to capture segregation dynamics across both flow geometries with a single set of parameters.

The orientational dynamics of a spherical magnetic particle in linear shear flow subjected to an oscillating magnetic field in the flow plane is analysed in the viscous limit. The shear is in the $X$–$Y$ plane, the magnetic field is in the $X$ direction and the vorticity is perpendicular to the flow in the $Z$ direction. The relevant dimensionless groups are $\omega ^\ast$, the ratio of the frequency of the magnetic field and the strain rate, and $\varSigma$, the ratio of the magnetic and hydrodynamic torques. As $\varSigma$ is decreased, there is a transition from in-plane rotation, where the rotation is in the flow ($X$–$Y$) plane, to out-of-plane rotation, where the orientation vector is not necessarily in the $X$–$Y$ plane and the dynamics depends on the initial orientation. The particle rotation is phase-locked for in-plane rotation with discrete odd rotation number (number of rotations in one period of magnetic field oscillation), while the orbits are quasi-periodic with non-integer rotation number for out-of-plane rotation. For $\varSigma \gg 1$, regions of odd rotation number $n_o$ are bound by the lines $8 (n_o-1) \varSigma \omega ^\ast = 1$ and $8 (n_o+1) \varSigma \omega ^\ast = 1$, and there are discontinuous changes in the rotation number and mean and root-mean-square torque at these lines. For $\varSigma \ll 1$, the domains of in-plane rotation of finite width in the $\omega ^\ast$–$\varSigma$ plane extend into downward cusps at $\omega ^\ast = {1}/{2 n_o}$. The orbits are quasi-periodic between these domains, where the rotation is out of plane.

Shock–droplet interaction in the early stage involves intricate wave structures. Investigating this phenomenon is inherently challenging due to the fine spatial and temporal scales involved. Past research has suggested that the occurrence of cavitation, marked by a negative peak pressure, is linked to the focus of the reflected expansion wave. In this study, a high-fidelity compressible numerical approach is utilized to replicate the initial phase of shock–droplet interactions. The location of the negative peak pressure is meticulously documented and compared with experimental measurement and numerical results. Results indicate a strong alignment between the negative peak pressure positions identified through numerical simulations and the focal points identified in theoretical models for low gas–liquid wave velocity ratios. However, this alignment is notably disrupted when dealing with higher gas–liquid wave velocity ratios. Further enhancements are made to the theoretical model, enabling a more precise depiction of internal wave structures and focus points, particularly under conditions of high gas–liquid wave velocity ratios. The study delves into the various factors influencing internal pressure fluctuations within the liquid droplet, categorizing them into four phases: the shock wave effect, relaxation effect, fluctuation effect, and expansion wave effect. Analysing the pressure decrease portion reveals that while the converging of the reflected expansion wave leads to a substantial pressure drop, it accounts for only a fraction of the total pressure variation. Consequently, any model predicting negative peak pressure positions must comprehensively consider all contributing factors.

In particle-laden turbulent wall flows, lift forces can influence the near-wall turbulence. This has been observed recently in particle-resolved simulations, which, however, are too expensive to be used in upscaled models. Instead, point-particle simulations have been the method of choice to simulate the dynamics of these flows during the last decades. While this approach is simpler, cheaper and physically sound for small inertial particles in turbulence, some issues remain. In the present work, we address challenges associated with lift force modelling in turbulent wall flows and the impact of lift forces in the near-wall flow. We performed direct numerical simulations of small inertial point particles in turbulent channel flow for fixed Stokes number and mass loading while varying the particle size. Our results show that the particle dynamics in the buffer region, causing the apparent particle-to-fluid slip velocity to vanish, raises major challenges for modelling lift forces accurately. While our results confirm that lift forces have little influence on particle dynamics for sufficiently small particle sizes, for inner-scaled diameters of order one and beyond, lift forces become quite important near the wall. The different particle dynamics under lift forces results in the modulation of streamwise momentum transport in the near-wall region. We analyse this lift-induced turbulence modulation for different lift force models, and the results indicate that realistic models are critical for particle-modelled simulations to correctly predict turbulence modulation by particles in the near-wall region.

Seafloor roughness profoundly influences the pattern and dynamics of large-scale oceanic flows. However, these kilometre-scale topographic patterns are unresolved by global numerical Earth system models and will remain subgrid for the foreseeable future. To properly represent the effects of small-scale bathymetry in analytical and coarse-resolution numerical models, we develop the stratified ‘sandpaper’ theory of flow–topography interaction. This model, which is based on the multilayer shallow-water framework, extends its barotropic antecedent to stratified flows. The proposed theory is successfully tested on the configuration representing the interaction of a zonal current with a corrugated cross-flow ridge.

We study the spreading of Newtonian viscous (aqueous glycerin solution) and viscoelastic (aqueous polymer solution) drops on solid substrates with different wettabilities. For drops of the same zero-shear viscosity, we find in the early stages of spreading that viscoelastic drops (i) spread faster and (ii) their contact radius shows a different power law vs time than Newtonian drops. We argue that the effect of viscoelasticity is only observable for experimental time scales of the order of or larger than the internal relaxation time of the viscoelastic polymer solution. We attribute this behaviour to the shear thinning of the viscoelastic polymer solution. When approaching the contact line, the shear rate increases and the steady-state viscosity of the viscoelastic drop is lower than that of the Newtonian drop. We support our experimental findings with a simple (first-order) perturbation model that qualitatively agrees with our findings.