Unconsolidated soils typically develop a physical surface crust after wetting and drying. We reproduced this process in the laboratory by wetting with fog and simulated rain on fallow agricultural soils from 26 locations, representing 15 soil types from Pinal County, Arizona. Through correlative analyses, we found that carbonate content was a strong predictor of physical crust strength with fog (p < 0.0001, R2 = 0.48) and rain (p = 0.004, R2 = 0.30). Clay content increased crust strength (p = 0.04) but was not a useful predictor. Our results extend the current understanding of the soil crusting process by highlighting the preeminence of carbonate cementation in desert agricultural soils. Consequently, we identify carbonate as a pragmatic tool for estimating crust strength, a surrogate measure of a soil’s potential to produce fugitive dust, which can help prioritize interventions to curb airborne dust in arid lands.

Vapour bubbles produced by long-pulsed laser often have complex non-spherical shapes that reflect some characteristics of the laser beam. The transition between two commonly observed shapes, namely, a rounded pear-like shape and an elongated conical shape, is studied using a new computational model that combines compressible multiphase fluid dynamics with laser radiation and phase transition. Two laboratory experiments are simulated, in which a holmium:YAG or thulium fibre laser is used to generate bubbles of different shapes. In both cases, the predicted bubble nucleation and morphology agree reasonably well with the experimental observation. The full-field results of laser irradiance, temperature, velocity and pressure are analysed to explain bubble dynamics and energy transmission. It is found that due to the lasting energy input, the vapour bubble's dynamics is driven not only by advection, but also by the continued vaporisation at its surface. Vaporisation lasts less than $1~{\rm \mu}$s in the case of the pear-shaped bubble, compared with over $50~{\rm \mu}$s for the elongated bubble. It is thus hypothesised that the bubble's morphology is determined by competition. When the speed of advection is higher than that of vaporisation, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem, then estimated for the experiments using simulation results. The results support the hypothesis. They also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation.

We consider Euler flows on two-dimensional (2-D) periodic domain and are interested in the stability, both linear and nonlinear, of a simple equilibrium given by the 2-D Taylor–Green vortex. As the first main result, numerical evidence is provided for the fact that such flows possess unstable eigenvalues embedded in the band of the essential spectrum of the linearized operator. However, the unstable eigenfunction is discontinuous at the hyperbolic stagnation points of the base flow and its regularity is consistent with the prediction of Lin (Intl Math. Res. Not., vol. 2004, issue 41, 2004, pp. 2147–2178). This eigenfunction gives rise to an exponential transient growth with the rate given by the real part of the eigenvalue followed by passage to a nonlinear instability. As the second main result, we illustrate a fundamentally different, non-modal, growth mechanism involving a continuous family of uncorrelated functions, instead of an eigenfunction of the linearized operator. Constructed by solving a suitable partial differential equation (PDE) optimization problem, the resulting flows saturate the known estimates on the growth of the semigroup related to the essential spectrum of the linearized Euler operator as the numerical resolution is refined. These findings are contrasted with the results of earlier studies of a similar problem conducted in a slightly viscous setting where only the modal growth of instabilities was observed. This highlights the special stability properties of equilibria in inviscid flows.

Three-dimensional (3-D) clustering characteristics of large-Stokes-number sprays interacting with turbulent swirling co-flows are investigated experimentally. The astigmatic interferometric particle imaging (AIPI) technique is used for simultaneous measurement of the spray droplets position in 3-D space and their corresponding diameter. The Stokes number estimated based on the Kolmogorov time scale varies from 34 to 142. The results show that the degree of droplet clustering plateaus at about 0.4 and at large Stokes numbers. It is obtained that the mean length scale of the clusters normalized by the Kolmogorov length scale follows a power-law relation, and the mean void length scale normalized by the integral length scale plateaus at about 1.5 and at large Stokes numbers. It is shown that the ratio of the number density of the droplets residing within the clusters to the global number density increases with increasing Stokes number and is about 8 for the largest Stokes number examined in this study. The joint characteristics of cluster's normalized volume and the mean diameter of droplets residing within the clusters show that small-volume clusters accommodate droplets with a relatively broad range of diameters. However, large clusters carry droplets with the most probable diameter. The developed AIPI technique in the present study and the corresponding spray characteristics are of importance for engineering applications that aim to understand the 3-D clustering characteristics of large-Stokes-number droplets sprayed into turbulent swirling co-flows.

Liquid droplet dynamics are widely used in biological and engineering applications, which contain complex interfacial instabilities and pattern formation such as droplet merging, splitting and transport. This paper studies a class of mean field control formulations for these droplet dynamics, which can be used to control and manipulate droplets in applications. We first formulate the droplet dynamics as gradient flows of free energies in modified optimal transport metrics with nonlinear mobilities. We then design an optimal control problem for these gradient flows. As an example, a lubrication equation for a thin volatile liquid film laden with an active suspension is developed, with control achieved through its activity field. Lastly, we apply the primal–dual hybrid gradient algorithm with high-order finite-element methods to simulate the proposed mean field control problems. Numerical examples, including droplet formation, bead-up/spreading, transport, and merging/splitting on a two-dimensional spatial domain, demonstrate the effectiveness of the proposed mean field control mechanism.

The depth-integrated horizontal momentum equations and continuity equation are employed to develop a new model. The vertical velocity and pressure can be expressed exactly in terms of horizontal velocities and free-surface elevation, which are the only unknowns in the model. Dividing the water column into elements and approximating horizontal velocities using linear shape function in each element, a set of model equations for horizontal velocities at element nodes is derived by adopting the weighted residual method. These model equations can be applied for transient or steady free-surface flows by prescribing appropriate lateral boundary conditions and initial conditions. Here, only the wave–current–bathymetry interaction problems are investigated. Theoretical analyses are conducted to examine various linear wave properties of the new models, which outperform the Green–Naghdi-type models for the range of water depth to wavelength ratios and the Boussinesq-type models as they are capable of simulating vertically sheared currents. One-dimensional horizontal numerical models, using a finite-difference method, are applied to a wide range of wave–current–bathymetry problems. Numerical validations are performed for nonlinear Stokes wave and bichromatic wave group propagation in deep water, sideband instability, regular wave transformation over a submerged shoal and focusing wave group interacting with linearly sheared currents in deep water. Very good agreements are observed between numerical results and laboratory data. Lastly, numerical experiments of wave shoaling from deep to shallow water are conducted to further demonstrate the capability of the new model.

The research on elasto-inertial turbulence (EIT), a new type of turbulent flow, has reached the stage of identifying the minimal flow unit (MFU). On this issue, direct numerical simulations of FENE-P fluid flow in two-dimensional channels with variable sizes are conducted in this study. We demonstrate with the increase of channel length that the simulated flow experiences several different flow patterns, and there exists an MFU for EIT to be self-sustained. At Weissenberg number ($Wi$) higher than the one required to excite EIT, when the channel length is relatively small, a steady arrowhead regime (SAR) flow structure and a laminar-like friction coefficient is achieved. However, as the channel length increases, the flow can fully develop into EIT characterized with high flow drag. Close to the size of the MFU, the simulated flow behaves intermittently between the SAR state with low drag and EIT state with high drag. The flow falling back to ‘laminar flow’ is caused by the insufficient channel size below the MFU. Furthermore, we give the relationship between the value of the MFU and the effective $Wi$, and explain its physical reasons. Moreover, the intermittent flow regime obtained based on the MFU gives us an opportunity to look into the origin and exciting process of EIT. Through capturing the onset process of EIT, we observed that EIT originates from the sheet-like extension structure located near the wall, which is maybe related to the wall mode rather than the centre mode. The fracture and regeneration of this sheet-like structure is the key mechanism for the self-sustaining of EIT.

The dynamics of a shock-induced separation unit generated by a 20$^\circ$ sharp fin placed on a cylindrical surface in a Mach 2.5 flow was investigated. Specifically, the present work investigated the mechanisms that govern the mid-frequency range of separation shock unsteadiness in the fin shock wave–boundary layer interaction (SBLI) unit. Two-dimensional pressure fields were obtained over the cylinder surface spanning the entire fin SBLI unit using high-bandwidth pressure-sensitive paint at 40 kHz imaging rate that allowed probing the low- through mid-frequency ranges of the separation shock unsteadiness. The mean pressure field showed a progressive weakening of the separation shock with downstream distance, which is an artifact of the three-dimensional relief offered by the curved mounting surface. The root-mean-square (r.m.s.) pressure field exhibited a banded structure with elevated $p_{r.m.s.}$ levels beneath the intermittent region, separation vortex and adjacent to the fin root. The power spectral density (PSD) of the surface pressure fluctuations obtained beneath the intermittent region revealed that the separation shock oscillations exhibited the mid-frequency content over the majority of its length. Interestingly, neither the PSD nor the length of the intermittent region varied noticeably with downstream distance, revealing a constant separation shock foot velocity along the entire SBLI. The pressure fluctuation PSD beneath the separation vortex also exhibited the broadband peak at the mid-frequency range of the separation shock motions over the majority of its length within the measurement domain. By contrast, the region adjacent to the fin root exhibited pressure oscillations at a substantially lower frequency compared with the separation shock and the separation vortex. Two-point coherence and cross-correlation analysis provided unique insights into the critical sources and mechanisms that drive the separation shock unsteadiness. The separation vortex and separation shock dynamics were found to be driven by a combination of convecting perturbations that originated from the vicinity of the fin leading edge and the local interactions of the separated flow with the incoming boundary layer. The boundary layer locally strengthened or weakened the convecting pressure perturbations depending on the local momentum fluctuations within the boundary layer.

We report direct numerical simulations results of the rough-wall channel, focusing on roughness with high $k_{rms}/k_a$ statistics but small to negative $Sk$ statistics, and we study the implications of this new dataset on rough-wall modelling. Here, $k_{rms}$ is the root mean square, $k_a$ is the first-order moment of roughness height, and $Sk$ is the skewness. The effects of packing density, skewness and arrangement of roughness elements on mean streamwise velocity, equivalent roughness height ($z_0$) and Reynolds and dispersive stresses have been studied. We demonstrate that two-point correlation lengths of roughness height statistics play an important role in characterizing rough surfaces with identical moments of roughness height but different arrangements of roughness elements. Analysis of the present as well as historical data suggests that the task of rough-wall modelling is to identify geometric parameters that distinguish the rough surfaces within the calibration dataset. We demonstrate a novel feature selection procedure to determine these parameters. Further, since there is no finite set of roughness statistics that distinguish between all rough surfaces, we argue that obtaining a universal rough-wall model for making equivalent sand-grain roughness ($k_s$) predictions would be challenging, and that each rough-wall model would have its applicable range. This motivates the development of group-based rough-wall models. The applicability of multi-variate polynomial regression and feedforward neural networks for building such group-based rough-wall models using the selected features has been shown.

We study the sedimentation of U-shaped circular disks in the Stokes limit of vanishing inertia. We simulate the flow past such disks using a finite-element-based solution of the three-dimensional Stokes equations, accounting for the integrable singularities that develop along their edges. We show that the purely vertical sedimentation of such disks in their upright (upside-down) U orientation is unstable to perturbations about their pitching (rolling) axes. The instability is found to depend only weakly on the size of the container in which the disks sediment, allowing us to analyse their behaviour based on the resistance matrix which governs the evolution of the disk's six rigid-body degrees of freedom in an unbounded fluid. We show that the governing equations can be reduced to two ordinary differential equations which describe the disk's inclination against the direction of gravity. A phase-plane analysis, the results of which are in good agreement with experiments, reveals that the two instabilities generally cause the disk to sediment along complex spiral trajectories while it alternates between pitching- and rolling-dominated motions. The chirality of the trajectories is set by the initial conditions rather than the (non-chiral) shape of the disk. For certain initial orientations, the disk retains its inclination and sediments along a perfectly helical path. The observed behaviour is fundamentally different from that displayed by flat circular disks which sediment without any reorientation. We therefore study the effect of variations in the disk's curvature to show how in the limit of vanishing curvature the behaviour of a flat disk is recovered.