We employ the methods of statistical mechanics to obtain closures for the balance equations of momentum and fluctuation kinetic energy that govern the ballistic motion of grains rebounding at a rigid, bumpy bed that are driven by turbulent or non-turbulent shearing fluids, in the absence of mid-trajectory collisions and fluid velocity fluctuations. We obtain semi-analytical solutions for steady and fully developed saltation over horizontal beds for the vertical profiles of particle concentration and stresses and fluid and particle velocities. These compare favourably with measurements in discrete-element numerical simulations in the wide range of conditions of Earth and other planetary environments. The predictions of the particle horizontal mass flux and its scaling with the amount of particles in the system, the properties of the carrier fluid and the intensity of the shearing also agree with numerical simulations and wind-tunnel experiments.

A distinguishing feature of the bi-stable wake is that the wake persists in either of two preferred states for a sufficiently long time. Aiming to understand the persistence mechanism, this paper numerically investigates the airwake characteristics of the Chalmers ship model (CSM) using large eddy simulation with a wall-adapting local-eddy viscosity model and is complemented by experimental testings for validations. There are two cases of interest: (i) the baseline CSM with a sharp-edged superstructure front that induces massive boundary layer separation; (ii) the front-rounded (FR) CSM with suppressed flow separation. During a characteristic time ($t^*$) of 1142 (26.5 s), the baseline case has a frequently switching wake, whereas the FR wake maintains a stable asymmetric structure with only one switch attempt. To understand the different wake behaviours, the study starts by analysing wake flow structures, vortex cores and the wake dynamics, followed by investigating the instantaneous flow physics. Results suggest that the baseline wake has a weak bi-stable pattern, whereas the FR wake behaves similarly to a reflectional symmetry breaking state of a potential bi-stable wake. The wake switching is found to be driven by the tilting of (vertical-oriented) $z$-vorticity sheets from either side of the base toward the centre. This tilting behaviour is subjected to the high-magnitude vorticity that sheds from the upstream flow separation at the front sharp edges. With the sharp edges rounded in the FR case, the upstream vorticity is mitigated and the tilting effect is significantly reduced, leading to a more stable wake structure. The reasoning provided in the paper potentially explains the persistence mechanism of the bi-stable wake.

We optimize jet mixing using large eddy simulations (LES) at a Reynolds number of $3000$. Key methodological enablers consist of Bayesian optimization, a surrogate model enhanced by deep learning and persistent data topology for physical interpretation. The mixing performance is characterized by an equivalent jet radius ($R_{eq}$) derived from the streamwise velocity in a plane located $8$ diameters downstream. The optimization is performed in a 22-dimensional actuation space that comprises most known excitations. This search space parameterizes the distributed actuation imposed on the bulk flow and at the periphery of the nozzle in the streamwise and radial directions. The momentum flux measures the energy input of the actuation. The optimization quadruples the jet radius $R_{eq}$ with a $7$-armed blooming jet after around $570$ evaluations. The control input requires $2\,\%$ momentum flux of the main flow, which is one order of magnitude lower than an ad hoc dual-mode excitation. Intriguingly, a pronounced suboptimum in the search space is associated with a double-helix jet, a new flow pattern. This jet pattern results in a mixing improvement comparable to the blooming jet. A state-of-the-art Bayesian optimization converges towards this double-helix solution. The learning is accelerated and converges to another better optimum by including a deep-learning-enhanced surrogate model trained along the optimization. Persistent data topology extracts the global and many local minima in the actuation space. These minima can be identified with flow patterns beneficial to the mixing.

The flow resistance, i.e. friction factor times Reynolds number ($\,f\,{Re}$), of longitudinal-fin heat sinks with or without clearance between a shroud and the tips of the fins is an important parameter in thermal design. This is because it dictates the caloric resistance of the heat sink, i.e. change in bulk temperature of the fluid flowing through it. When there is no clearance and the common and oft-valid assumption of negligible fin thickness is invoked, $f\,{Re}$ corresponds to simply that of a rectangular duct. However, with clearance, only numerical results are available as per the well-known study by Sparrow, Baliga and Patankar (ASME J. Heat Transfer, vol. 100, 1978). We develop analytical formulae for $f\,{Re}$ for fully developed flow with clearance. The exact solution is provided by an integral formula derived via conformal mappings. Additionally, simple formulae are derived via asymptotic expansions in three cases: (1) the fin spacing is small compared to the fin height and clearance; (2) the clearance is small compared to the fin spacing, which is small compared to the fin height; (3) the same as case (2) but valid for larger clearances. The different asymptotic formulae are compared to the exact formula, and together cover almost the entire relevant parameter range (for fin spacing and clearance) with errors of less than 15 %. A formula for the limiting case of no clearance is shown to be more accurate, for any fin spacing, than a widely used correlation from the literature.

Direct numerical simulations (DNS) were carried out to investigate flow control for transition delay using steady blowing/suction strips at the wall of a flared cone at Mach 6 and zero angle of attack. For the numerical investigations of the transition control strategy, the flared cone geometry and the flow conditions of the experiments in the Boeing/Air Force Office of Scientific Research (AFOSR) Mach 6 Quiet Tunnel (BAM6QT) at Purdue University were chosen. For the DNS, transition was initiated by introducing random disturbances at the inflow of the computational domain, emulating ‘natural’ transition in wind-tunnel experiments caused by free-stream noise. In both wind-tunnel experiments and numerical simulations, streamwise ‘hot’ streaks were found on the surface of the flared cone, which are caused by a nonlinear interaction of an axisymmetric second-mode wave and a pair of oblique waves of the same frequency (‘fundamental resonance’). The objective of the flow control strategy proposed here is to delay the transition onset, and thus mitigate the negative consequences associated with the nonlinear transition stages, i.e. the development of hot streaks and large wall-pressure amplitudes that were observed in experiments and DNS. Our previous so-called ‘controlled’ transition simulations have shown that flow control using steady blowing and suction strips can lead to a significant delay of the hot streak development on the surface of the flared cone. The simulation results presented in this paper show that this flow control strategy remains effective, even in a natural transition scenario characterized by broadband disturbances.

We present a theory to describe the Nusselt number, $\operatorname {\mathit {Nu}}$, corresponding to the heat or mass flux, as a function of the Rayleigh–Darcy number, $\operatorname {\mathit {Ra}}$, the ratio of buoyant driving force over diffusive dissipation, in convective porous media flows. First, we derive exact relationships within the system for the kinetic energy and the thermal dissipation rate. Second, by segregating the thermal dissipation rate into contributions from the boundary layer and the bulk, which is inspired by the ideas of the Grossmann and Lohse theory (J. Fluid Mech., vol. 407, 2000; Phys. Rev. Lett., vol. 86, 2001), we derive the scaling relation for $\operatorname {\mathit {Nu}}$ as a function of $\operatorname {\mathit {Ra}}$ and provide a robust theoretical explanation for the empirical relations proposed in previous studies. Specifically, by incorporating the length scale of the flow structure into the theory, we demonstrate why heat or mass transport differs between two-dimensional and three-dimensional porous media convection. Our model is in excellent agreement with the data obtained from numerical simulations, affirming its validity and predictive capabilities.

Different types of neural networks have been used to solve the flow sensing problem in turbulent flows, namely to estimate velocity in wall-parallel planes from wall measurements. Generative adversarial networks (GANs) are among the most promising methodologies, due to their more accurate estimations and better perceptual quality. This work tackles this flow sensing problem in the vicinity of the wall, addressing for the first time the reconstruction of the entire three-dimensional (3-D) field with a single network, i.e. a 3-D GAN. With this methodology, a single training and prediction process overcomes the limitation presented by the former approaches based on the independent estimation of wall-parallel planes. The network is capable of estimating the 3-D flow field with a level of error at each wall-normal distance comparable to that reported from wall-parallel plane estimations and at a lower training cost in terms of computational resources. The direct full 3-D reconstruction also unveils a direct interpretation in terms of coherent structures. It is shown that the accuracy of the network depends directly on the wall footprint of each individual turbulent structure. It is observed that wall-attached structures are predicted more accurately than wall-detached ones, especially at larger distances from the wall. Among wall-attached structures, smaller sweeps are reconstructed better than small ejections, while large ejections are reconstructed better than large sweeps as a consequence of their more intense footprint.

The nonlinear dynamics of a helical vortex disturbed by a long-wave-instability mode is studied by direct numerical simulation. Vortex reconnection or self-reconnection of the helical vortex is shown to play a crucial role depending on the pitch of the helical vortex. For the larger pitch, a vortex ring is created after the vortex reconnection; it detaches from the remaining helical vortex, whose pitch is doubled. A vortex ring is also created for the smaller pitch; however, it forms a linked system with the remaining vortex. The topological constraint due to this linkage forces strong interaction between the different parts of the helical vortex, leading to turbulent transition.

Numerical simulations of thermoelectrohydrodynamic convection in a dielectric liquid inside a finite-length cylindrical annulus with a fixed temperature difference have been performed with increasing high-frequency electric tension under microgravity conditions. The electric field, coupled with the permittivity gradient, generates a dielectrophoretic buoyancy force whose non-conservative part can induce thermoelectric convection in the liquid. The liquid remains in a conductive state below a critical value of the applied electric voltage. At a critical value, a supercritical bifurcation occurs from the conductive state to a convective state made of stationary helicoidal vortices. A further increase of electric voltage leads to oscillatory helicoidal vortices and then to wavy patterns before spoke patterns dominate the convective flow. The dielectrophoretic force is shown to enhance the heat transfer from the hot to cold walls due to induced convective flows. Particularly, these results demonstrate that the dielectrophoretic buoyancy force holds promise to replace the gravitational force to induce efficient heat transfer in microgravity conditions, and they contribute to a better fundamental understanding of heat transfer in microgravity.

Recent theoretical progress using multiscale asymptotic analysis has revealed various possible regimes of stratified turbulence. Notably, buoyancy transport can either be dominated by advection or diffusion, depending on the effective Péclet number of the flow. Two types of asymptotic models have been proposed, which yield measurably different predictions for the characteristic vertical velocity and length scale of the turbulent eddies in both diffusive and non-diffusive regimes. The first, termed a ‘single-scale model’, is designed to describe flow structures having large horizontal and small vertical scales, while the second, termed a ‘multiscale model’, additionally incorporates flow features with small horizontal scales, and reduces to the single-scale model in their absence. By comparing predicted vertical velocity scaling laws with direct numerical simulation data, we show that the multiscale model correctly captures the properties of strongly stratified turbulence within regions dominated by small-scale isotropic motions, whose volume fraction decreases as the stratification increases. Meanwhile its single-scale reduction accurately describes the more orderly, layer-like, quiescent flow outside those regions.

This study explores the dynamics of flexible ribbons with an added weight $G$ at the tail in uniform flow, considering key parameters like inflow Reynolds number ($Re_u$), mass ratio ($M_t$) and aspect ratio (${A{\kern-4pt}R}$). For two-dimensional ribbons, a simplified theoretical model accurately predicts equilibrium configurations and forces. Inspired by Barois & De Langre (J. Fluid Mech., vol. 735, 2013, R2), we introduce an important control parameter ($C_G$) that effectively collapses normalized forces and angle data. Vortex-induced vibration is observed, and Strouhal number ($St$) scaling laws with $C_G$ are identified. In three-dimensional scenarios, the model effectively predicts lift, but its accuracy in predicting drag is limited to situations with small $Re_u$ values. The flow along the side edges mitigates pressure differences, thereby suppressing vibration and uplift, particularly noticeable in the case of narrow ribbons. This study offers new insights into the dynamics of flexible bodies in uniform flow.

Previous studies have shown that low-frequency vortex oscillations occur around a hemisphere–cylinder body at different angles of attack, but the underlying mechanism is still unclear. In this study, we examine the origin of the vortex oscillation using numerical simulations and global linear stability analysis. The vortex oscillation is reproduced using numerical simulations, and the oscillatory modes are computed through dynamic mode decomposition (DMD). We obtain the base flow through a selective frequency damping method, which exhibits a pair of steady leeward vortices over the body. The four unstable modes are computed using a modified Arnoldi iteration. The antisymmetric mode with a Strouhal number of 0.105 is discovered to be responsible for the alternate oscillation of the vortex pair, and the mode with a Strouhal number of 0.220 corresponds to the in-phase vortex oscillation. Their frequencies have good agreement with the modes of DMD. The other two unstable modes with higher frequencies, one antisymmetric and one symmetric, are harmonic frequencies of the above two modes. The study conclusively verifies that the vortex oscillation over a hemisphere–cylinder body originates from a global flow instability.

Understanding settling motion of coral grains is important in terms of protection of coral reef systems and resotoration of the associated ecosystems. In this paper, a series of laboratory experiments was conducted to investigate the settling motion, using optical microscopy to measure shape parameters of coral grains and the particle-filtering-based object tracking to reconstruct the three-dimensional trajectory. Three characteristic descent regimes, namely, tumbling, chaotic and fluttering, are classified based on the three-dimensional trajectory, the spiral radius variation and the velocity spectrum. It is demonstrated that if one randomly picks up one coral grain, then the probabilities of occurrence of the three regimes are approximately $26\,\%$, $42\,\%$ and $32\,\%$, respectively. We have shown that first, the dimensionless settling velocity generally increases with the non-dimensional diameter and Corey shape factor and second, the drag coefficient generally decreases with the Reynolds number and Corey shape factor. Based on this, the applicability of existing models on predicting settling velocity and drag coefficient for coral grains is demonstrated further. Finally, we have proposed extended models for predicting the settling velocity. This study contributes to better understanding of settling motion and improves our predictive capacity of settling velocity for coral grains with complex geometry.

Saturated flow film boiling on a sphere has been numerically studied in this work for both vertical and horizontal flow configurations. The simulations were performed using a numerical methodology developed by the authors for boiling flows on three-dimensional unstructured meshes. For interface capturing, the coupled level set and volume of fluid method is used. The interface evolution, vapour wake dynamics and heat transfer have been thoroughly investigated by varying the saturated liquid flow velocity, sphere diameter and wall superheat. The relative importance of both the buoyancy and the inertial forces is described in terms of the Froude number $(Fr)$. The vapour bubble evolves periodically at low $Fr$ values, while a stable vapour column develops at high $Fr$ values. The interface evolution pattern obtained in the present work is in good agreement with the results of experimental studies available in the literature. For all the values of $Fr$, a stable vapour column develops for a large-diameter sphere and releases vapour bubbles of varying sizes. Furthermore, for a large-diameter sphere, surface capillary waves are observed at the interface, similar to the observations of some of the experimental studies available in the literature. The flow in the liquid and vapour wakes appears to be strongly coupled. The heat transfer in the present work is estimated using the spatially and temporally averaged Nusselt numbers. Finally, an fast Fourier transform analysis of the space-averaged Nusselt number reveals a strong interaction among the different forces.

We investigate the phenomena of crater formation and gas release caused by projectile impact on underwater beds, which occurs in many natural, geophysical and industrial applications. The bed in our experiment is constructed of hydrophobic particles, which trap a substantial amount of air in the pores of the bed. In contrast to dry beds, the air–water interface in a submerged bed generates a granular skin that provides rigidity to the medium by producing skin over the bulk. The projectile's energy is used to reorganize the grains, which causes the skin to crack, allowing the trapped air to escape. The morphology of the craters as a function of impact energy in submerged beds exhibits different scaling laws than what is known for dry beds. This phenomenon is attributed to the contact line motion on the hydrophobic fractal-like surface of submerged grains. The volume of the gas released is a function of multiple factors, chiefly the velocity of the projectile, depth of the bed and depth of the water column.

Small finite-size particles suspended in fluid flow through an enclosed curved duct can focus to points or periodic orbits in the two-dimensional duct cross-section. This particle focusing is due to a balance between inertial lift forces arising from axial flow and drag forces arising from cross-sectional vortices. The inertial particle focusing phenomenon has been exploited in various industrial and medical applications to passively separate particles by size using purely hydrodynamic effects. A fixed size particle in a circular duct with a uniform rectangular cross-section can have a variety of particle attractors, such as stable nodes/spirals or limit cycles, depending on the radius of curvature of the duct. Bifurcations occur at different radii of curvature, such as pitchfork, saddle-node and saddle-node infinite period (SNIPER), which result in variations in the location, number and nature of these particle attractors. By using a quasi-steady approximation, we extend the theoretical model of Harding et al. (J. Fluid Mech., vol. 875, 2019, pp. 1–43) developed for the particle dynamics in circular ducts to spiral duct geometries with slowly varying curvature, and numerically explore the particle dynamics within. Bifurcations of particle attractors with respect to radius of curvature can be traversed within spiral ducts and give rise to a rich nonlinear particle dynamics and various types of tipping phenomena, such as bifurcation-induced tipping (B-tipping), rate-induced tipping (R-tipping) and a combination of both, which we explore in detail. We discuss implications of these unsteady dynamical behaviours for particle separation and propose novel mechanisms to separate particles by size in a non-equilibrium manner.

Depinning of liquid droplets on substrates by flow of a surrounding immiscible fluid is central to applications such as cross-flow microemulsification, oil recovery and waste cleanup. Surface roughness, either natural or engineered, can cause droplet pinning, so it is of both fundamental and practical interest to determine the flow strength of the surrounding fluid required for droplet depinning on rough substrates. Here, we develop a lubrication-theory-based model for droplet depinning on a substrate with topographical defects by flow of a surrounding immiscible fluid. The droplet and surrounding fluid are in a rectangular channel, a pressure gradient is imposed to drive flow and the defects are modelled as Gaussian-shaped bumps. Using a precursor-film/disjoining-pressure approach to capture contact-line motion, a nonlinear evolution equation is derived describing the droplet thickness as a function of distance along the channel and time. Numerical solutions of the evolution equation are used to investigate how the critical pressure gradient for droplet depinning depends on the viscosity ratio, surface wettability and droplet volume. Simple analytical models are able to account for many of the features observed in the numerical simulations. The influence of defect height is also investigated, and it is found that, when the maximum defect slope is larger than the receding contact angle of the droplet, smaller residual droplets are left behind at the defect after the original droplet depins and slides away. The model presented here yields considerably more information than commonly used models based on simple force balances, and provides a framework that can readily be extended to study more complicated situations involving chemical heterogeneity and three-dimensional effects.

Wall temperature has a significant effect on shock wave/turbulent boundary layer interactions (STBLIs) and has become a non-negligible factor in the design process of hypersonic vehicles. In this paper, direct numerical simulations are conducted to investigate the wall temperature effects on STBLIs over a 34° compression ramp at Mach number 6. Three values of the wall-to-recovery-temperature ratio (0.50, 0.75 and 1.0) are considered in the simulations. The results show that the size of the separation bubble declines significantly as the wall temperature decreases. This is because the momentum profile of the boundary layer becomes fuller with wall cooling, which means the near-wall fluid has a greater momentum to suppress flow separation. An equation based on the free-interaction theory is proposed to predict the distributions of the wall pressure upstream of the corner at different wall temperatures. The prediction results are generally consistent with the simulation results (Reynolds number Reτ ranges from 160 to 675). In addition, the low-frequency unsteadiness is studied through the weighted power spectral density of the wall pressure and the correlation between the upstream and downstream. The results indicate that the low-frequency motion of the separation shock is mainly driven by the downstream mechanism and that wall cooling can significantly suppress the low-frequency unsteadiness, including the strength and streamwise range of the low-frequency motions.

This study proposes a novel super-resolution (or SR) framework for generating high-resolution turbulent boundary layer (TBL) flow from low-resolution inputs. The framework combines a super-resolution generative adversarial neural network (SRGAN) with down-sampling modules (DMs), integrating the residual of the continuity equation into the loss function. The DMs selectively filter out components with excessive energy dissipation in low-resolution fields prior to the super-resolution process. The framework iteratively applies the SRGAN and DM procedure to fully capture the energy cascade of multi-scale flow structures, collectively termed the SRGAN-based energy cascade reconstruction framework (EC-SRGAN). Despite being trained solely on turbulent channel flow data (via ‘zero-shot transfer’), EC-SRGAN exhibits remarkable generalization in predicting TBL small-scale velocity fields, accurately reproducing wavenumber spectra compared to direct numerical simulation (DNS) results. Furthermore, a super-resolution core is trained at a specific super-resolution ratio. By leveraging this pretrained super-resolution core, EC-SRGAN efficiently reconstructs TBL fields at multiple super-resolution ratios from various levels of low-resolution inputs, showcasing strong flexibility. By learning turbulent scale invariance, EC-SRGAN demonstrates robustness across different TBL datasets. These results underscore the potential of EC-SRGAN for generating and predicting wall turbulence with high flexibility, offering promising applications in addressing diverse TBL-related challenges.

Resonant standing waves excited on the water surface in a deep narrow rectangular cavity by a fully immersed cylinder harmonically oscillating in the vertical direction are studied theoretically and experimentally. The effect of the finite wavemaker size is considered in the framework of the potential two-dimensional flow theory. Nonlinearities and weak dissipation at solid surfaces are accounted for. The spatio-temporal structure of the waves in the presence of detuning between the forcing and the natural frequency of the system is analysed. The variation of the surface shape in space and time studied in experiments supports the assumption of two-dimensional flow. The finite size of the wavemaker causes a downshift of the effective resonant frequency of the cavity; this effect is enhanced by the nonlinearity. For small amplitude waves, the surface elevation evolution in time is decomposed into the sum of the time-periodic function, corresponding to the forcing frequency, and its second harmonic; the shape of the wavenumber spectra of these components depends on the forcing frequency. For larger wave amplitudes, additional peaks in the frequency spectrum appear. The theoretical predictions are compared with the experimental results.