Submesoscale fronts with large horizontal buoyancy gradients and $O(1)$ Rossby numbers are common in the upper ocean. These fronts are associated with large vertical transport and are hotspots for biological activity. Submesoscale fronts are susceptible to symmetric instability (SI) – a form of stratified inertial instability which can occur when the potential vorticity is of the opposite sign to the Coriolis parameter. Here, we use a weakly nonlinear stability analysis to study SI in an idealised frontal zone with a uniform horizontal buoyancy gradient in thermal wind balance. We find that the structure and energetics of SI strongly depend on the front strength, defined as the ratio of the horizontal buoyancy gradient to the square of the Coriolis frequency. Vertically bounded non-hydrostatic SI modes can grow by extracting potential or kinetic energy from the balanced front and the relative importance of these energy reservoirs depends on the front strength and vertical stratification. We describe two limiting behaviours as ‘slantwise convection’ and ‘slantwise inertial instability’ where the largest energy source is the buoyancy flux and geostrophic shear production, respectively. The growing linear SI modes eventually break down through a secondary shear instability, and in the process transport considerable geostrophic momentum. The resulting breakdown of thermal wind balance generates vertically sheared inertial oscillations and we estimate the amplitude of these oscillations from the stability analysis. We finally discuss broader implications of these results in the context of current parameterisations of SI.

This paper presents an experimental study on how both variable solid volume fractions and aspect ratios (length/width) of a centre-channel rectangular porous patch under aligned configuration of rigid and emergent stems impact the flow behaviour and wake structure. This study forms an essential extension to the existing fundamentals and knowledge on this topic. Through rigorous experimental tests by velocity measurement and dye visualization, the aspect ratio, rarely addressed before, is confirmed to play a critical role. Vortex street, unable to be triggered under a low solid volume fraction, however, can be generated by elongating the patch (increasing the aspect ratio). The key reason is that patch elongation promotes the generation of the wake vortex street by producing a relatively high transverse velocity gradient in the wake region. Meanwhile, Kelvin–Helmholtz vortex streets are triggered along the two patch lateral edges, re-increasing the in-patch velocity and imposing contributions to the wake vortex streets generation. By scaling the characteristic velocity (at the wake vortex initiation position and patch trailing edge) and solid volume fraction with the patch aspect ratio, three non-dimensional threshold maps can be established to express the combined effects of the solid volume fraction and aspect ratio on the initiation of the wake vortex street. They could be alternatively used for theoretical analysis and implementation on wake formation and structure subject to parameter availability.

Engineering flow systems operating under low pressures and/or at the micro/nano scale generally include a physically adsorbed gas layer next to the surface. In this paper, we develop a scattering kernel that accounts for the effect of adsorption, arising from van der Waals interactions, on the dynamics of molecules impinging on solid smooth surfaces. In the limit of low bulk density, surface adsorption becomes negligible and the scattering kernel recovers consistently the Cercignani–Lampis model, which best describes molecular collisions with a clean, smooth surface. In the limit of high bulk density, a dense adsorbed molecular layer forms next to the surface and its presence is picked up by the Maxwell model with complete diffuse reflection, which better captures the multiple collisions suffered by molecules. A weight coefficient based on the Langmuir adsorption isotherm is incorporated into the modelling to handle the transition between these two limiting conditions of low and high densities. The proposed model is validated against high-fidelity molecular dynamics simulations that are performed for a variety of gas–surface combinations and adsorbed molecular layers with different densities. It is shown that the proposed model very well captures the scattering patterns of beams of gas molecules at different velocities impinging on surfaces, as well as momentum and energy accommodation coefficients in the entire range of explored conditions.

The effect of perforation on the wake of a thin flat plate placed normal to the free stream at Reynolds number ($Re$) 250 (based on plate width $d$, and inflow velocity $U_o$) is studied by means of direct numerical simulation. The perforated plate of length $6d$ consist of six equidistant square holes of varying sizes corresponding to porosity $\beta$ (ratio of open area to total plate area) of 0 %, 4 %, 9 %, 12.25 %, 16 %, 20.25 % and 25 %. It is observed that the bleed or jet flow through perforations pushes the shear layer interaction farther downstream with increasing $\beta$. This causes a monotonic decrease in the drag coefficient with increasing porosity, and a sharp fall seeming to begin at $\beta \approx 4\,\%$. On the other hand, the Strouhal number increases with $\beta$ up to 16 % (at $\beta =16\,\%$, loss of flow three-dimensionality leads to a ‘quasi-laminar’ state of flow). This is followed by a sharp fall in the Strouhal number at $\beta \approx 20\,\%$. The behaviour of the large-scale vortical structures in the far wake is influenced by the near-wake behaviour of the bleed flow, where the local $Re$ based on the perforation hole size determines the overall flow three-dimensionality. It is also observed that the jet or bleed flow undergoes meandering instability when pitch separation is equivalent to the hole size (at $\beta =25\,\%$). The low-$Re$ turbulent flow for a non-perforated plate is altered to a transitional state by the presence of perforation. The streamwise vortex pairs (secondary instabilities) become fairly organized as $\beta$ is increased from 0 % to 16 %. The secondary instability at $\beta =16\,\%$ appears similar to mode-B with wavelength ${\approx }1d$. On the contrary, the secondary instability at $\beta =25\,\%$ appears similar to mode-A with a wavelength of ${\approx }2d$.

We study the propagation of water waves in a three-dimensional device alternating open canals and resonant canals with subwavelength resonances. The dispersion of water waves in such a medium is obtained by analysing the full three-dimensional problem and combining Bloch–Floquet analysis with an asymptotic technique. We obtain the closed forms of the dispersions for resonant canals containing one or two resonators, which depend on only two functions associated with symmetric and antisymmetric modes, and on a geometric parameter analogous to the hopping parameter in topological systems. The analysis of the complete band structure reveals frequency ranges alternating between elliptical and hyperbolic dispersions; in particular, the hyperbolic regime gives rise to a negative effective water depth with a consequent negative refraction. Throughout the course of our study, our theoretical results are validated by comparison with numerical calculations of the full three-dimensional problem.

A new two-equation model for gravity-driven liquid film flow based on the long-wave expansion has been derived. The novelty of the model consists in using a base velocity profile combining parabolic (Ruyer-Quil & Manneville, Eur. Phys. J. B, vol. 15, issue 2, 2000, pp. 357–369) and ellipse (Usha et al., Phys. Fluids, vol. 32, issue 1, 2020, 013603) profile functions in the wall-normal coordinate. The dependence on a free parameter $A$ related to the eccentricity of an ellipse serves as an adjustable parameter. The resulting models are consistent at $O(\varepsilon )$ for inertia terms and at $O(\varepsilon ^2)$ for viscous diffusion effects, and predict accurately the primary instability. Appropriate tuning of the adjustable parameter helps to recover accurate predictions for the asymptotic wave celerity of nonlinear solitary waves. Further, the model is shown to capture the closed separation vortices that can form underneath the troughs of precursory capillary ripples.

Advective trapping occurs when solute enters low velocity zones in heterogeneous porous media. Classical local modelling approaches combine the impact of slow advection and diffusion into a hydrodynamic dispersion coefficient and many temporally non-local approaches lump these mechanisms into a single memory function. This joint treatment makes parameterization difficult and thus prediction of large-scale transport a challenge. Here, we investigate the mechanisms of advective trapping and their impact on transport in media composed of a high conductivity background and isolated low permeability inclusions. Breakthrough curves show that effective transport changes from a streamtube-like behaviour to genuine random trapping as the degree of disorder of the inclusion arrangement increases. We upscale this behaviour using a Lagrangian view point, in which idealized solute particles transition over a fixed distance at random advection times combined with Poissonian advective trapping events. We discuss the mathematical formulation of the upscaled model in the continuous time random walk and mobile–immobile mass transfer frameworks, and derive a model for large-scale solute non-Fickian dispersion. These findings give new insight into transport in highly heterogeneous media.

The density distribution around a sphere descending in a salt-stratified fluid is measured by the laser-induced fluorescence (LIF) method. The corresponding velocity distribution is measured by particle image velocimetry (PIV), and numerical simulation is also performed to supplement the observations by LIF and PIV. In steady flow, LIF observes a thin and vertically long structure which corresponds to a buoyant jet. The bell-shaped structure, which appears under strong stratification and moderate Reynolds number (Froude number $Fr \lesssim 3$, Reynolds number $50 \lesssim Re \lesssim 500$), is also identified. The measured density distributions in the salinity boundary layer and in the jet agree with the numerical simulations which use the Schmidt number of the fluorescent dye ($Sc \sim 2000$). The initially unsteady process of the jet formation is also investigated. Under weak stratification, the LIF shows an initial development of an axisymmetric rear vortex as observed in homogeneous fluids. However, as time proceeds and the effect of stratification becomes significant, the vortex shrinks and disappears, while the jet extends vertically upward. Under strong stratification, a thin jet develops without generating a rear vortex, since the effect of stratification becomes significant in a short time before the vortex is generated.

We investigate how inhomogeneity influences the $k^{-5/3}$ inertial range scaling of turbulent kinetic energy spectra (with $k$ the wavenumber). For weak statistical inhomogeneity, the energy spectrum can be described as an equilibrium spectrum plus a perturbation. Theoretical arguments suggest that this latter contribution scales as $k^{-7/3}$. This prediction is assessed using direct numerical simulations of three-dimensional Kolmogorov flow.

A model is constructed to describe the flow field and arbitrary deformation of a drop or vesicle that contains and is embedded in an electrolyte solution, where the flow and deformation are caused by an applied electric field. The applied field produces an electrokinetic flow, which is set up on the charge-up time scale $\tau _{*c}=\lambda _{*} a_{*}/D_{*}$, where $\lambda _{*}$ is the Debye screening length, $a_{*}$ is the inclusion length scale and $D_{*}$ is an ion diffusion constant. The model is based on the Poisson–Nernst–Planck and Stokes equations. These are reduced or simplified by forming the limit of strong electrolytes, for which dissolved salts are completely ionised in solution, together with the limit of thin Debye layers. Debye layers of opposite polarity form on either side of the drop interface or vesicle membrane, together forming an electrical double layer. Two formulations of the model are given. One utilises an integral equation for the velocity field on the interface or membrane surface together with a pair of integral equations for the electrostatic potential on the outer faces of the double layer. The other utilises a form of the stress-balance boundary condition that incorporates the double layer structure into relations between the dependent variables on the layers’ outer faces. This constitutes an interfacial boundary condition that drives an otherwise unforced Stokes flow outside the double layer. For both formulations relations derived from the transport of ions in each Debye layer give additional boundary conditions for the potential and ion concentrations outside the double layer.

We present an experimental investigation of motion onset in simple yield stress fluids. In this context, motion onset refers to the transition from the motionless steady state to a steady flow, as well as the development of motion in a fluid initially at rest. We consider the natural convection of carbopol microgels in a square cavity with differentially heated sidewalls. We use particle image velocimetry and thermometry to reveal the evolution of both temperature and velocity fields. It is a hallmark of yield stress fluids that a critical ratio of the yield stress and buoyancy stresses exists above which the steady state is motionless. We observe this critical behaviour in our experiments. Contrary to the theoretical predictions, however, systematic motion is evident at the onset of all experiments, even when the steady state is motionless. Above the critical limit, extremely slow motion is observed immediately after the onset of the experiment. This is followed by very slow decay to rest, reminiscent of creep behaviour. Below the critical limit, the initial slow dynamics is followed by flow development patterns similar to theoretical predictions based on the Bingham model. We show that motion onset in carbopol microgels is dominated by subyield motion and fluidization, key processes that are not captured by viscoplastic models.