Showing papers in "Journal of Fluid Mechanics in 2020"

The flow physics of COVID-19

Abstract: Flow physics plays a key role in nearly every facet of the COVID-19 pandemic. This includes the generation and aerosolization of virus-laden respiratory droplets from a host, its airborne dispersion and deposition on surfaces, as well as the subsequent inhalation of these bioaerosols by unsuspecting recipients. Fluid dynamics is also key to preventative measures such as the use of face masks, hand washing, ventilation of indoor environments and even social distancing. This article summarizes what we know and, more importantly, what we need to learn about the science underlying these issues so that we are better prepared to tackle the next outbreak of COVID-19 or a similar disease.

Effects of ventilation on the indoor spread of COVID-19

Abstract: Although the relative importance of airborne transmission of the SARS-CoV-2 virus is controversial, increasing evidence suggests that understanding airflows is important for estimation of the risk of contracting COVID-19. The data available so far indicate that indoor transmission of the virus far outstrips outdoor transmission, possibly due to longer exposure times and the decreased turbulence levels (and therefore dispersion) found indoors. In this paper we discuss the role of building ventilation on the possible pathways of airborne particles and examine the fluid mechanics of the processes involved.

Surfactant dynamics: hidden variables controlling fluid flows

Abstract: Surfactants - molecules and particles that preferentially adsorb to fluid interfaces - play a ubiquitous role in the fluids of industry, of nature, and of life. Since most surfactants cannot be seen directly, their behavior must be inferred from their impact on observed flows, like the buoyant rise of a bubble, or the thickness of a coating film. In so doing, however, a difficulty arises: physically distinct surfactant processes can affect measurable flows in qualitatively identical ways, raising the specter of confusion or even misinterpretation. This Perspective describes, in one coherent piece, both the equilibrium properties and dynamic processes of surfactants, to better enable the fluid mechanics community to understand, interpret, and design surfactant/fluid systems. Specifically, §2 treats the equilibrium thermodynamics of surfactants at interfaces, including surface pressure, isotherms of soluble and insoluble surfactants, and surface dilatational moduli (Gibbs and Marangoni). §3 describes surfactant dynamics in fluid systems, including surfactant transport and interfacial stress boundary conditions, the competition between surface diffusion, advection, and adsorption/desorption, Marangoni stresses and flows, and surface excess rheology. §4 discusses paradigmatic problems from fluid mechanics that are impacted by surfactants, including translating drops and bubbles, surfactant adsorption to clean and oscillating interfaces; capillary wave damping, thin film dynamics, foam drainage, and the dynamics of particles and probes at surfactant-laden interfaces. Finally, §5 discusses the additional richness and complexity that frequently arise in 'real' surfactants, including phase transitions, phase coexistence, and polycrystalline phases within surfactant monolayers, and their impact on non-Newtonian surface rheology.

Prediction of turbulent heat transfer using convolutional neural networks

Abstract: With the recent rapid development of artificial intelligence (AI) and wide applications in many areas, some fundamental questions in turbulence research can be addressed, such as: ‘Can turbulence be learned by AI? If so, how and why?’ In order to provide answers to these questions, we applied deep learning to the prediction of turbulent heat transfer based only on wall information using data obtained from direct numerical simulations (DNS) of turbulent channel flow. Through this attempt, we investigated whether deep learning could help to improve our understanding of the physics of turbulent heat transfer. Under the assumption that the wall-normal local heat flux can be explicitly expressed through a multilayer nonlinear network in terms of the nearby wall-shear stresses and wall pressure fluctuations, we applied convolutional neural networks (CNNs) to predict the local heat flux. After optimizing the network hyperparameters using a grid searching method, we found that the network can predict the heat flux very accurately with a correlation coefficient of 0.980 between the DNS and the prediction by CNN for the trained Reynolds number, , and shows similar accuracy at a Reynolds number three times higher than the trained number. This result indicates that relationships between the local heat flux and the nearby inputs are quite insensitive to the Reynolds number within the tested range. In addition, observing the gradient maps of the trained network, we identified the part of the input data that is essential for the local heat flux prediction and the spatial relationship between the local heat flux and the nearby input fields. In addition to obtaining an understanding of the underlying physics, we investigated whether our model could be utilized for turbulence modelling.

Fragmentation versus Cohesion

Abstract: Capillarity is the familiar manifestation of the cohesion of liquids Since Laplace (Traite de mecanique celeste, vol IV, supplement au livre X: Sur l’action capillaire, 1805, pp 1–65), we know that intense attractive forces between the molecules bridge the small with the large as they shape liquid/vapour interfaces at the macroscopic scale through the concept of surface tension (menisci, drops, bubbles, puddles, liquid rise in tubes, etc …) We concentrate on situations where liquids ‘disgregate’, following the neologism of Clausius (Phil Mag, vol 24 (159), 1862, pp 81–97), meaning that they fragment by the action of deformation stresses whose intensity competes with that of cohesion forces Various examples, including explosions, blow-ups, hard and soft impacts and shears applied to liquid jets, sheets and drops are reviewed They concern applications ranging from liquid propulsion, agricultural spraying, to the formation of ocean spray, raindrops and human exhalations by violent respiratory events In spite of their diversity, the various modes of fragment production share an ultimate common phenomenology – the ligament dynamics – suggesting that the final stable droplet size distribution can be interpreted from elementary principles

Dynamics and control of premixed combustion systems based on flame transfer and describing functions

Abstract: This article describes recent progress on premixed flame dynamics interacting with acoustic waves. Expressions are derived to determine the stability of combustors with respect to thermoacoustic oscillations. The validity of these expressions is general, but they are illustrated in laminar systems. Laminar burners are commonly used to elucidate the response of premixed flames to incoming flow perturbations, highlight the role of acoustic radiation in their stability, identify modes associated with thermoacoustic intrinsic instabilities and decipher the leading mechanisms in annular systems with multiple injectors. Many industrial devices also operate in a laminar premixed mode such as, for example, domestic gas boilers and heaters equipped with matrix burners for material processing in which unconfined flames are stabilized at one extremity of the system. This article proposes a systematic approach to determine the stability of all these systems with respect to thermoacoustic oscillations by highlighting the key role of the burner impedance and the flame transfer function (FTF). This transfer function links in frequency space incoming flow perturbations to heat release rate disturbances. This concept can be used in the turbulent flame case as well. Weakly nonlinear stability analysis can also easily be conducted by replacing the FTF by a flame describing function in the expressions derived in this work. The response of premixed flames to harmonic mixture compositions and flow-rate perturbations is then revisited and the main parameters controlling the FTF are described. A theoretical framework is finally developed to reduce the system thermoacoustic sensitivity by tailoring the FTF.

Nonlinear mode decomposition with convolutional neural networks for fluid dynamics

Abstract: We present a new nonlinear mode decomposition method to visualize decomposed flow fields, named the mode decomposing convolutional neural network autoencoder (MD-CNN-AE). The proposed method is applied to a flow around a circular cylinder at the Reynolds number as a test case. The flow attributes are mapped into two modes in the latent space and then these two modes are visualized in the physical space. Because the MD-CNN-AEs with nonlinear activation functions show lower reconstruction errors than the proper orthogonal decomposition (POD), the nonlinearity contained in the activation function is considered the key to improving the capability of the model. It is found by applying POD to each field decomposed using the MD-CNN-AE with hyperbolic tangent activation such that a single nonlinear MD-CNN-AE mode contains multiple orthogonal bases, in contrast to the linear methods, i.e. POD and MD-CNN-AE with linear activation. We further assess the proposed MD-CNN-AE by applying it to a transient process of a circular cylinder wake in order to examine its capability for flows containing high-order spatial modes. The present results suggest a great potential for the nonlinear MD-CNN-AE to be used for feature extraction of flow fields in lower dimensions than POD, while retaining interpretable relationships with the conventional POD modes.

A momentum-conserving wake superposition method for wind farm power prediction

Abstract: Analytical wind turbine wake models and wake superposition methods are prevailing tools widely adopted by the wind energy community to predict the power production of wind farms. However, none of the existing wake superposition methods conserve the streamwise momentum. In this study, a novel wake superposition method capable of conserving the total momentum deficit in the streamwise direction is derived theoretically, and its performance is validated with both particle imaging velocimetry measurements and large-eddy simulation results. Detailed inter-method comparisons show that the novel wake superposition method outperforms all the existing methods by delivering an accurate prediction of the power production and the centreline wake velocity deficit, with a typical error of less than 5 % (excluding the near-wake region). Additionally, the momentum-conserving wake superposition method is extended to combine the transverse velocities induced by yawed wind turbines, and the secondary wake steering effect crucial to the power optimization in active wake control is well reproduced.

Low-order model for successive bifurcations of the fluidic pinball

Abstract: We propose the first least-order Galerkin model of an incompressible flow undergoing two successive supercritical bifurcations of Hopf and pitchfork type. A key enabler is a mean-field consideration exploiting the symmetry of the mean flow and the asymmetry of the fluctuation. These symmetries generalize mean-field theory, e.g. no assumption of slow growth rate is needed. The resulting five-dimensional Galerkin model successfully describes the phenomenogram of the fluidic pinball, a two-dimensional wake flow around a cluster of three equidistantly spaced cylinders. The corresponding transition scenario is shown to undergo two successive supercritical bifurcations, namely a Hopf and a pitchfork bifurcation on the way to chaos. The generalized mean-field Galerkin methodology may be employed to describe other transition scenarios.

Is vortex stretching the main cause of the turbulent energy cascade

Abstract: In three-dimensional turbulence there is on average a cascade of kinetic energy from the largest to the smallest scales of the flow. While the dominant idea is that the cascade occurs through the process of vortex stretching, evidence for this is debated. Here we show theoretically and numerically that vortex stretching is in fact not the main contributor to the average cascade. The main contributor is the self-amplification of the strain-rate field, and we provide several arguments for why its role must not be conflated with that of vortex stretching. Numerical results, however, indicate that vortex stretching plays a more important role during fluctuations of the cascade about its average behaviour. We also resolve a paradox regarding the differing role of vortex stretching on the energy cascade and energy dissipation rate dynamics.

Effect of adverse pressure gradients on turbulent wing boundary layers

Abstract: The characteristics of turbulent boundary layers (TBLs) subjected to adverse pressure gradients are analysed through well-resolved large-eddy simulations. The geometries under study are the NAC ...

Spectral proper orthogonal decomposition and resolvent analysis of near-wall coherent structures in turbulent pipe flows

Abstract: Direct numerical simulations, performed with a high-order spectral-element method, are used to study coherent structures in turbulent pipe flow at friction Reynolds numbers . The database was analysed using spectral proper orthogonal decomposition (SPOD) to identify energetically dominant coherent structures, most of which turn out to be streaks and quasi-streamwise vortices. To understand how such structures can be modelled, the linear flow responses to harmonic forcing were computed using the singular value decomposition of the resolvent operator, using the mean field as a base flow. The SPOD and resolvent analysis were calculated for several combinations of frequencies and wavenumbers, allowing the mapping out of similarities between SPOD modes and optimal responses for a wide range of relevant scales in turbulent pipe flows. In order to explore physical reasons behind the agreement between both methods, an indicator of lift-up mechanism in the resolvent analysis was introduced, activated when optimal forcing is dominated by the wall-normal and azimuthal components, and associated response corresponds to streaks of streamwise velocity. Good agreement between leading SPOD and resolvent modes is observed in a large region of parameter space. In this region, a significant gain separation is found in resolvent analysis, which may be attributed to the strong amplification associated with the lift-up mechanism, here understood as nonlinear forcing terms leading to the appearance of streamwise vortices, which in turn form high-amplitude streaks. For both Reynolds numbers, the observed concordances were generally for structures with large energy in the buffer layer. The results highlight resolvent analysis as a pertinent reduced-order model for coherent structures in wall-bounded turbulence, particularly for streamwise elongated structures corresponding to near-wall streamwise vortices and streaks.

Global instability of wing shock-buffet onset

Abstract: Shock buffet on wings encountered in edge-of-the-envelope transonic flight remains an unresolved and disputed flow phenomenon, challenging both fundamental fluid mechanics and applied aircraft aerodynamics. Its dynamics is revealed through the interaction of spanwise shock-wave oscillations and intermittent turbulent boundary-layer separation. Resulting unsteady aerodynamic loads, and their mutual working with the flexible aircraft structure, need to be accounted for in establishing the safe flight envelope. The question of global instability leading to this flow unsteadiness is addressed herein. It is shown for the first time on an industrially relevant configuration that the dynamics of a single unstable oscillatory eigenmode plays a prominent role in near-onset shock buffet on a quasi-rigid wing. Its three-dimensional spatial structure, previously inferred both from experiment and time-marching simulation, describes a spanwise-localised pocket of shear-layer pulsation synchronised with an outboard-propagating shock oscillation. The results also suggest that the concept of a critical global shock-buffet mode commonly reported for two-dimensional aerofoils also applies to three-dimensional finite and swept wings, albeit different modes at play. Specifically, the modern wing design, NASA Common Research Model, with publicly available geometry and experimental data for code validation is studied at a free-stream Mach number of 0.85 with Reynolds number per reference chord of and varying angle of attack between 3. 5° and 4. 0° targeting the instability onset. Strouhal number at instability onset just above 3. 7° is approximately 0.39. At the same time, a band of eigenmodes shows reduced decay rate in the Strouhal-number range of 0.3 to 0.7, with additional unstable oscillatory modes appearing beyond onset. Importantly, those emerging modes seem to discretise the continuous band of medium-wavelength modes, as recently reported for infinite swept wings using stability analysis, hence generalising those findings to finite wings. Through conventional time-marching unsteady simulation it is explored how the critical linear eigenmode feeds into the nonlinearly saturated limit-cycle oscillation near instability onset. The established numerical strategy, using an iterative inner–outer Krylov approach with shift-and-invert spectral transformation and sparse iterative linear solver, to solve the arising large-scale eigenvalue problem with an industrial Reynolds-averaged Navier–Stokes flow solver means that such a practical non-canonical test case at a high-Reynolds-number condition can be investigated. The numerical findings can potentially be exploited for more effective unsteady flow analysis in future wing design and inform routes to flow control and model reduction.

Wall-Attached and wall-detached eddies in wall-bounded turbulent flows

Abstract: According to Townsend’s attached eddy hypothesis (AEH), a boundary layer flow is comprised of wall-attached eddies, but to extract the part of the flow whose statistical behaviours are well described by the AEH is not at all straightforward. The objective of this work is to extract the part of the flow that can be described by the AEH, and study the statistical behaviours of the other part, which cannot be described by the AEH. In this process, two types of eddies are identified in addition to the Kolmogorov-scale eddies, i.e. wall-attached eddies and wall-detached eddies. The statistical behaviours of the wall-attached eddies are shown to be very well described by the AEH, i.e. the eddies are wall-attached, self-similar and, importantly, Gaussian, whereas the wall-detached eddies cannot be modelled by the AEH. Specifically, a decomposition scheme is proposed following Townsend (The Structure of Turbulent Shear Flow, Cambridge University Press, 1976). We apply our decomposition scheme to three different flows, i.e. channel, boundary layer and atmospheric surface layer flows. The results are similar with only quantitative differences, suggesting possible universality in both the wall-attached eddies and the wall-detached eddies.

A physical model of turbulence cascade via vortex reconnection sequence and avalanche

Abstract: Viscous anti-parallel vortex reconnection is studied by means of direct numerical simulation for vortex Reynolds numbers ( , circulation/viscosity) up to 40 000. To suppress the inherent symmetry breaking due to the Kelvin–Helmholtz (planar jet) instability, as prevalent in prior studies, and to better explore the progression of the mechanism details, the simulation is performed by imposing symmetry and using double-precision arithmetic. We show, for the first time, the evidence of vortex reconnection cascade scenario initially proposed by Melander and Hussain (CTR Report, 1988), who suggested that the remnant threads, following the first reconnection, undergo successive reconnections in a cascade. Secondary reconnection (the details distinctly captured and visualized at a lower ) leads to the successive generation of numerous small-scale structures, including vortex rings, hairpin-like vortex packets and vortex tangles. As increases, the third and higher generations of reconnection form a turbulent cloud avalanche consisting of a tangle of fine vortices. The energy is rapidly transferred to finer scales during reconnection, and a distinct - 5/3 inertial range is observed for the kinetic energy spectrum, associated with numerous resulting fine-scale bridgelets and thread filaments. In addition, we also discover an inverse cascade at large scales through the accumulation of bridgelets. The separation distance before the first reconnection is found to scale as , which is different from the typical 1/2 scaling for classical and quantum vortex filament reconnections. Both peak enstrophy and its production rate grow with faster than the power law suggested by Hussain and Duraisamy (Phys. Fluids, vol. 23, 2011, 021701). Our simulations not only reveal the detailed mechanisms of high- reconnection, but also shed light on the physics of turbulence cascade and present the reconnection avalanche as a realistic physical model for turbulence cascade.

Data-driven decomposition of the streamwise turbulence kinetic energy in boundary layers. Part 1. Energy spectra

Abstract: In wall-bounded turbulence, a multitude of coexisting turbulence structures form the streamwise velocity energy spectrum from the viscosity- to the inertia-dominated range of scales. Definite scaling-trends for streamwise spectra have remained empirically elusive, although a prominent school of thought stems from the works of Perry & Abell (J. Fluid Mech., vol. 79, 1977, pp. 785–799) and Perry et al. (J. Fluid Mech., vol. 165, 1986, pp. 163–199), which were greatly inspired by the attached-eddy hypothesis of Townsend (The Structure of Turbulent Shear Flow, Cambridge University Press, 1976). In this paper, we re-examine the turbulence kinetic energy of the streamwise velocity component in the context of the spectral decompositions of Perry and co-workers. Two universal spectral filters are derived via spectral coherence analysis of two-point velocity signals, spanning a Reynolds-number range .

Resolvent-based estimation of space-time flow statistics

Abstract: We develop a method to estimate space–time flow statistics from a limited set of known data. While previous work has focused on modelling spatial or temporal statistics independently, space–time statistics carry fundamental information about the physics and coherent motions of the flow and provide a starting point for low-order modelling and flow control efforts. The method is derived using a statistical interpretation of resolvent analysis. The central idea of our approach is to use known data to infer the statistics of the nonlinear terms that constitute a forcing on the linearized Navier–Stokes equations, which in turn imply values for the remaining unknown flow statistics through application of the resolvent operator. Rather than making an a priori assumption that the flow is dominated by the leading singular mode of the resolvent operator, as in some previous approaches, our method allows the known input data to select the most relevant portions of the resolvent operator for describing the data, making it well suited for high-rank turbulent flows. We demonstrate the predictive capabilities of the method, which we call resolvent-based estimation, using two examples: the Ginzburg–Landau equation, which serves as a convenient model for a convectively unstable flow, and a turbulent channel flow at low Reynolds number.

Competitive evaporation of multiple sessile droplets

Abstract: An asymptotic model is derived for the competitive diffusion-limited evaporation of multiple thin sessile droplets under the assumption that the droplets are well separated. Exact solutions of the model are obtained for a pair of and for a polygonal array of identical droplets, and the model is found to perform well even outside its formal range of validity, up to and including the limit of touching droplets. The shielding effect of droplets on each other is demonstrated, and the model is used to investigate the effect of this shielding on droplet evolutions and lifetimes, as well as on the coffee-ring effect. The theoretical predictions of the model are found to be in good agreement with recent experimental results for seven relatively closely spaced droplets, suggesting that the model could be a useful tool for studying a wide range of other droplet configurations.

Extreme wave statistics of long-crested irregular waves over a shoal

Abstract: We report laboratory experiments of long-crested irregular water surface waves propagating over a shoal. For a sufficiently shallow shoal we find that the surface elevation can have a local maximum of skewness and kurtosis above the shallower part of the shoal close to the edge on the incoming side, and a local minimum of skewness over the downward slope on the lee side of the shoal. We find that the horizontal fluid velocity can have a local maximum and minimum of skewness at the same locations as those for the surface elevation. However, the kurtosis of the horizontal fluid velocity can have a local maximum over the downward slope on the lee side of the shoal, different from the location of the maximum of kurtosis of the surface elevation.

Numerical simulations of flow past three circular cylinders in equilateral-triangular arrangements

Abstract: Flow past three identical circular cylinders is numerically investigated using the immersed boundary method. The cylinders are arranged in an equilateral-triangle configuration with one cylinder placed upstream and the other two side-by-side downstream. The focus is on the effect of the spacing ratio , Reynolds number and three-dimensionality on the flow structures, hydrodynamic forces and Strouhal numbers, where is the cylinder centre-to-centre spacing and is the cylinder diameter. The fluid dynamics involved is highly sensitive to both and , leading to nine distinct flow structures, namely single bluff-body flow, deflected flow, flip-flopping flow, steady symmetric flow, steady asymmetric flow, hybrid flow, anti-phase flow, in-phase flow and fully developed in-phase co-shedding flow. The time-mean drag and lift of each cylinder are more sensitive to than while fluctuating forces are less sensitive to than . The three-dimensionality of the flow affects the development of the wake patterns, changing the ranges of different flow structures. A diagram of flow regimes, together with the contours of hydrodynamic forces, in the space, is given, providing physical insights into the complex interactions of the three cylinders.

Modelling double emulsion formation in planar flow-focusing microchannels

Abstract: Double emulsion formation in a hierarchical flow-focusing channel is systematically investigated using a free-energy ternary lattice Boltzmann model. A three-dimensional formation regime diagram is constructed based on the capillary numbers of the inner ( ), middle ( ) and outer ( ) phase fluids. The results show that the formation diagram can be classified into periodic two-step region, periodic one-step region, and non-periodic region. By varying and in the two-step formation region, different morphologies are obtained, including the regular double emulsions, decussate regimes with one or two alternate empty droplets, and structures with multiple inner droplets contained in the continuous middle phase thread. Bidisperse behaviours are also frequently encountered in the two-step formation region. In the periodic one-step formation region, scaling laws are proposed for the double emulsion size and for the size ratio between the inner droplet and the overall double emulsion. Furthermore, we show that the interfacial tension ratio can greatly change the morphologies of the obtained emulsion droplets, and the channel geometry plays an important role in changing the formation regimes and the double emulsion sizes. In particular, narrowing the side inlets or the distance between the two side inlets promotes the conversion from the two-step formation regime to the one-step formation regime.

Noise reduction mechanisms of an open-cell metal-foam trailing edge

Abstract: Open-cell porous materials have been reported as a promising concept for mitigating turbulent boundary-layer trailing-edge noise. This manuscript examines the aeroacoustics of a porous trailing edge to study its noise reduction mechanisms. Numerical investigations have been carried out for a NACA 0018 aerofoil with three different types of trailing edge: a baseline solid trailing edge, a fully porous trailing edge and a blocked-porous variant in which a solid core is added at the symmetry plane. The latter prevents flow interaction between the two sides of the aerofoil. Flow-field solutions are obtained by solving the explicit, transient and compressible lattice-Boltzmann equation, while the Ffowcs-Williams and Hawkings acoustic analogy has been used to compute far-field noise. The porous material is modelled using an equivalent fluid region governed by Darcy's law, in which the properties of a Ni-Cr-Al open-cell metal foam are applied. The simulation results are validated against reference data from experiments. The regular porous trailing edge reduces noise substantially, particularly at low frequency, whereas the blocked variant retains similar noise characteristics as the solid one. By employing a beamforming technique, the dominant source is found at the trailing edge for the solid and blocked trailing edges, while for the fully porous one, the dominant source is located near the solid-porous junction. The analysis of the scattered sound suggests that the permeability of the porous trailing edge allows for acoustic scattering along the porous medium surface that promotes destructive interference, and in turn, attenuates far-field noise intensity. The spectra and spanwise coherence of surface pressure fluctuations at the trailing edge are hardly affected by the presence of the porous material, which are found to be insufficient to justify the noise reduction. The flow field inside the porous medium is also examined to explain the differences between the fully porous and blocked-porous trailing edges. While the mean velocity components are similar for both, substantial difference is found for the velocity fluctuations. The impedance of the porous medium is computed as the ratio of velocity and pressure fluctuations. Unlike the blocked variant, the impedance in the fully porous trailing edge gradually decreases along the downstream direction, which leads to the distributed noise scattering along the porous medium surface. Additionally, the scattering efficiency at the actual trailing edge location is reduced due to the smaller impedance discontinuity.

Velocity distributions, dispersion and stretching in three-dimensional porous media

Abstract: Using index matching and particle tracking, we measure the three-dimensional velocity field in an isotropic porous medium composed of randomly packed solid spheres. This high-resolution experimental dataset provides new insights into the dynamics of dispersion and stretching in porous media. Dynamic-range velocity measurements indicate that the distribution of the velocity magnitude, ). While such a distribution should lead to a persistent anomalous dispersion process for advected non-diffusive point particles, we show that the dispersion of non-diffusive tracers nonetheless becomes Fickian beyond a time set by the smallest effective velocity of the tracers. We derive expressions for the onset time of the Fickian regime and the longitudinal and transverse dispersion coefficients as a function of the velocity field properties. The experimental velocity field is also used to study, by numerical advection, the stretching histories of fluid material lines. The mean and the variance of the line elongations are found to grow exponentially in time and the distribution of elongation is log-normal. These results confirm the chaotic nature of advection within three-dimensional porous media. By providing the laws of dispersion and stretching, the present study opens the way to a complete description of mixing in porous media.

On the splashing of high-speed drops impacting a dry surface

Abstract: When a drop impacts a dry surface at high velocity, it atomises into secondary droplets. These small droplets are generated by one of two types of splashes: either by a prompt splash from the spreading rim at the surface or by a thin corona splash, which levitates from the surface. This study investigates the splashing mechanisms experimentally using multiple high-resolution cameras and characterises the outcome of both splashing types at high Weber and Reynolds numbers. We demonstrate that the prompt splash is well described by the Rayleigh–Taylor instability of the rapidly advancing liquid lamella and determine the boundaries defining this splashing regime, which allows us to distinguish the prompt from the corona splash. Furthermore, we provide an expression to estimate the elapsed time during which the secondary droplets are generated, which is then implemented in the theory of Riboux & Gordillo ( Phys. Rev. Lett. , vol. 113 (2), 2014, 024507). This theoretical approach together with detailed quantification of the splashing outcome allows us to completely predict the outcome of both splashing types, which includes the mean size, velocity and total ejected volume of the secondary droplets. The detailed model proposed here can be indeed used to understand, characterise and predict more accurately the underlying physics in several applications.

Data-driven discovery of governing equations for fluid dynamics based on molecular simulation

Abstract: The discovery of governing equations from data is revolutionizing the development of some research fields, where the scientific data are abundant but the well-characterized quantitative descriptions are probably scarce. In this work, we propose to combine the direct simulation Monte Carlo (DSMC) method, which is a popular molecular simulation tool for gas flows, and machine learning to discover the governing equations for fluid dynamics. The DSMC method does not assume any macroscopic governing equations a priori but just relies on the model of molecular interactions at the microscopic level. The data generated by DSMC are utilized to derive the underlying governing equations using a sparse regression method proposed recently. We demonstrate that this strategy is capable of deriving a variety of equations in fluid dynamics, such as the momentum equation, diffusion equation, Fokker–Planck equation and vorticity transport equation. The data-driven discovery not only provides the right forms of the governing equations, but also determines accurate values of the transport coefficients such as viscosity and diffusivity. This work proves that data-driven discovery combined with molecular simulations is a promising and alternative method to derive governing equations in fluid dynamics, and it is expected to pave a new way to establish the governing equations of non-equilibrium flows and complex fluids.

From zonal flow to convection rolls in Rayleigh-Bénard convection with free-slip plates

Abstract: Rayleigh–Benard (RB) convection with free-slip plates and horizontally periodic boundary conditions is investigated using direct numerical simulations. Two configurations are considered, one is two-dimensional (2-D) RB convection and the other one three-dimensional (3-D) RB convection with a rotating axis parallel to the plate, which for strong rotation mimics 2-D RB convection. For the 2-D simulations, we explore the parameter range of Rayleigh numbers from to and Prandtl numbers from to . The effect of the width-to-height aspect ratio is investigated for . We show that zonal flow, which was observed, for example, by Goluskin et al. (J. Fluid. Mech., vol. 759, 2014, pp. 360–385) for , is only stable when is smaller than a critical value, which depends on and . The regime in which only zonal flow can exist is called the first regime in this study. With increasing , we find a second regime in which both zonal flow and different convection roll states can be statistically stable. For even larger , in a third regime, only convection roll states are statistically stable and zonal flow is not sustained. How many convection rolls form (or in other words, what the mean aspect ratio of an individual roll is), depends on the initial conditions and on and . For instance, for and , the aspect ratio of an individual, statistically stable convection roll can vary in a large range between and . A convection roll with a large aspect ratio of , or more generally already with , can be seen as ‘localized’ zonal flow, and indeed carries over various properties of the global zonal flow. For the 3-D simulations, we fix and , and compare the flow for and . We first show that with increasing rotation rate both the flow structures and global quantities like the Nusselt number and the Reynolds number increasingly behave like in the 2-D case. We then demonstrate that with increasing aspect ratio , zonal flow, which was observed for small by von Hardenberg et al. (Phys. Rev. Lett., vol. 15, 2015, 134501), completely disappears for . For such large , only convection roll states are statistically stable. In-between, here for medium aspect ratio , the convection roll state and the zonal flow state are both statistically stable. What state is taken depends on the initial conditions, similarly as we found for the 2-D case.

Transfer of mass and momentum at rough and porous surfaces

Abstract: The surface texture of materials plays a critical role in wettability, turbulence and transport phenomena. In order to design surfaces for these applications, it is desirable to characterise non-smooth and porous materials by their ability to exchange mass and momentum with flowing fluids. While the underlying physics of the tangential (slip) velocity at a fluid–solid interface is well understood, the importance and treatment of normal (transpiration) velocity and normal stress is unclear. We show that, when the slip velocity varies at an interface above the texture, a non-zero transpiration velocity arises from mass conservation. The ability of a given surface texture to accommodate a normal velocity of this kind is quantified by a transpiration length . We further demonstrate that normal momentum transfer gives rise to a pressure jump. For a porous material, the pressure jump can be characterised by so-called resistance coefficients . By solving five Stokes problems, the introduced measures of slip, transpiration and resistance can be determined for any anisotropic non-smooth surface consisting of regularly repeating geometric patterns. The proposed conditions are a subset of the effective boundary conditions derived from formal multi-scale expansion. We validate and demonstrate the physical significance of the effective conditions on two canonical problems – a lid-driven cavity and a turbulent channel flow, both with non-smooth bottom surfaces.

Single- and double-helix vortex breakdown as two dominant global modes in turbulent swirling jet flow

Abstract: In this paper, we study the shape and dynamics of helical coherent structures found in the flow field of an annular swirling jet undergoing vortex breakdown. The flow field is studied by means of time-resolved tomographic particle image velocimetry measurements. The obtained flow fields are analysed using both classic and spectral proper orthogonal decomposition. Despite the simple geometrical set-up of the annular jet, the flow field is very complex. Two distinct large-scale helical flow structures are identified: a single and a double helix, both co-rotating with the swirl direction, and it is revealed that these structures are not higher harmonics of each other. The structures have a relatively low energy content which makes it hard to separate them from other dynamics of the flow field, notably turbulent motions. Because of this, classic proper orthogonal decomposition fails to identify both structures properly. Spectral proper orthogonal decomposition, on the other hand, allows them to be identified accurately when the filter size is set at around eight times the precession period. The precession frequencies of the single and double helices correspond to Strouhal numbers of 0.273 and mode has been identified as a second harmonic of the first mode, while this study identifies both as two independent global modes. Secondly, this work shows that the simultaneous occurrence of multiple helical global modes is possible within a turbulent flow and their shapes and frequencies are very well predicted by mean field stability analysis. The latter finding is of general interest as it applies to a wide class of fluid problems dominated by multiple oscillatory structures.

Lift-up, Kelvin–Helmholtz and Orr mechanisms in turbulent jets

Abstract: Three amplification mechanisms present in turbulent jets, namely lift-up, Kelvin–Helmholtz and Orr, are characterized via global resolvent analysis and spectral proper orthogonal decomposition (SPOD) over a range of Mach numbers. The lift-up mechanism was recently identified in turbulent jets via local analysis by Nogueira et al. (J. Fluid Mech., vol. 873, 2019, pp. 211–237) at low Strouhal number (. Finally, SPOD and resolvent analyses of localized regions show that the lift-up mechanism is present throughout the jet, with a dominant azimuthal wavenumber inversely proportional to streamwise distance from the nozzle, with streaks of azimuthal wavenumber exceeding five near the nozzle, and wavenumbers one and two most energetic far downstream of the potential core.

Mechanics of bluff body drag reduction during transient near-wake reversals

Abstract: A combination of passive and active methods is used to manipulate the symmetry of the turbulent wake of an Ahmed body. Dedicated experiments to study the wake reversals occurring between symmetry-breaking states are performed. We show how transient symmetric states explored during the reversals differ from static symmetry-breaking states in the flow organization they induce. In asymmetric states, a strong interaction and coupling between the recirculating flow from one side and the shear layer from opposite side triggers shear-layer instabilities and their amplification. The resulting large-scale flow engulfment in the recirculation region increases the recirculation intensity and thus increases base drag. By contrast, during the wake reversals the disorganization of the recirculating flow leads to a transient symmetric state with prevented shear-layer interaction and triggering mechanism with a concomitant drag reduction of compared to symmetry-breaking states. Finally, recent experimental results for unbalanced wakes and methodologies addressing the control of wake asymmetries for drag reduction are discussed and interpreted based on the present findings. This leads us to provide new insights in control methods for wake asymmetries.