Localization-delocalization transition for light particles in turbulence.

Small bubbles in fluids rise to the surface due to Archimede's force. Remarkably, in turbulent flows this process is severely hindered by the presence of vortex filaments, which act as moving potential wells, dynamically trapping light particles and bubbles. Quantifying the statistical weights and roles of vortex filaments in turbulence is, however, still an outstanding experimental and computational challenge due to their small scale, fast chaotic motion, and transient nature. Here we show that, under the influence of a modulated oscillatory forcing, the collective bubble behavior switches from a dynamically localized to a delocalized state. Additionally, we find that by varying the forcing frequency and amplitude, a remarkable resonant phenomenon between light particles and small-scale vortex filaments emerges, likening particle behavior to a forced damped oscillator. We discuss how these externally actuated bubbles can be used as a type of microscopic probe to investigate the space-time statistical properties of the smallest turbulence scales, allowing to quantitatively measure physical characteristics of vortex filaments. We develop a superposition model that is in excellent agreement with the simulation data of the particle dynamics which reveals the fraction of localized/delocalized particles as well as characteristics of the potential landscape induced by vortices in turbulence. Our approach paves the way for innovative ways to accurately measure turbulent properties and to the possibility to control light particles and bubble motions in turbulence with potential applications to oceanography, medical imaging, drug/gene delivery, chemical reactions, wastewater treatment, and industrial mixing.

Efficient point-based simulation of four-way coupled particles in turbulence at high number density.

In many natural and industrial applications, turbulent flows encompass some form of dispersed particles. Although this type of multiphase turbulent flow is omnipresent, its numerical modeling has proven to be a remarkably challenging problem. Models that fully resolve the particle phase are computationally very expensive, strongly limiting the number of particles that can be considered in practice. This warrants the need for efficient reduced order modeling of the complex system of particles in turbulence that can handle high number densities of particles. Here we present an efficient method for point-based simulation of particles in turbulence that are four-way coupled. In contrast with traditional one-way coupled simulations, where only the effect of the fluid phase on the particle phase is modeled, this method additionally captures the back-reaction of the particle phase on the fluid phase, as well as the interactions between particles themselves. We focus on the most challenging case of very light particles or bubbles, which show strong clustering in the high-vorticity regions of the fluid. This strong clustering poses numerical difficulties which are systematically treated in our work. Our method is valid in the limit of small particles with respect to the Kolmogorov scales of the flow and is able to handle very large number densities of particles. This methods paves the way for comprehensive studies of the collective effect of small particles in fluid turbulence for a multitude of applications.

Coexistence of Defect Morphologies in Three-Dimensional Active Nematics.

We establish how active stress globally affects the morphology of disclination lines of a three-dimensional active nematic liquid crystal under chaotic flow. Thanks to a defect detection algorithm based on the local nematic orientation, we show that activity selects a crossover length scale in between the size of small defect loops and that of long and tangled defect lines of fractal dimension 2. This length scale crossover is consistent with the scaling of the average separation between defects as a function of activity. Moreover, on the basis of numerical simulation in a 3D periodic geometry, we show the presence of a network of regular defect loops, contractible onto the 3-torus, always coexisting with wrapping defect lines. While the length of regular defects scales linearly with the emerging active length scale, it verifies an inverse quadratic dependence for wrapping defects. The shorter the active length scale, the more the defect lines wrap around the periodic boundaries, resulting in extremely long and buckled structures.

Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions.

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

Extreme statistics and extreme events in dynamical models of turbulence.

We present a study of the intermittent properties of a shell model of turbulence with statistics of ∼10^{7} eddy turn over time, achieved thanks to an implementation on a large-scale parallel GPU factory. This allows us to quantify the inertial range anomalous scaling properties of the velocity fluctuations up to the 24th-order moment. Through a careful assessment of the statistical and systematic uncertainties, we show that none of the phenomenological and theoretical models previously proposed in the literature to predict the anomalous power-law exponents in the inertial range are in agreement with our high-precision numerical measurements. We find that at asymptotically high-order moments, the anomalous exponents tend toward a linear scaling, suggesting that extreme turbulent events are dominated by one leading singularity. We found that systematic corrections to scaling induced by the infrared and ultraviolet (viscous) cutoffs are the main limitations to precision for low-order moments, while high orders are mainly affected by the finite statistical samples.. The high-fidelity numerical results reported in this work offer an ideal benchmark for the development of future theoretical models of intermittency in dynamical systems for either extreme events (high-order moments) or typical fluctuations (low-order moments). For the latter, we show that we achieve a precision in the determination of the inertial range scaling exponents of the order of one part over ten thousand (fifth significant digit), which may be considered a record for out-of-equilibrium fluid-mechanics systems and models.

Pattern formation by turbulent cascades.

Fully developed turbulence is a universal and scale-invariant chaotic state characterized by an energy cascade from large to small scales at which the cascade is eventually arrested by dissipation1-6. Here we show how to harness these seemingly structureless turbulent cascades to generate patterns. Pattern formation entails a process of wavelength selection, which can usually be traced to the linear instability of a homogeneous state7. By contrast, the mechanism we propose here is fully nonlinear. It is triggered by the non-dissipative arrest of turbulent cascades: energy piles up at an intermediate scale, which is neither the system size nor the smallest scales at which energy is usually dissipated. Using a combination of theory and large-scale simulations, we show that the tunable wavelength of these cascade-induced patterns can be set by a non-dissipative transport coefficient called odd viscosity, ubiquitous in chiral fluids ranging from bioactive to quantum systems8-12. Odd viscosity, which acts as a scale-dependent Coriolis-like force, leads to a two-dimensionalization of the flow at small scales, in contrast with rotating fluids in which a two-dimensionalization occurs at large scales4. Apart from odd viscosity fluids, we discuss how cascade-induced patterns can arise in natural systems, including atmospheric flows13-19, stellar plasma such as the solar wind20-22, or the pulverization and coagulation of objects or droplets in which mass rather than energy cascades23-25.

Stochastic fluctuations of diluted pedestrian dynamics along curved paths.

As we walk towards our destinations, our trajectories are constantly influenced by the presence of obstacles and infrastructural elements; even in the absence of crowding our paths are often curved. Since the early 2000s pedestrian dynamics have been extensively studied, aiming at quantitative models with both fundamental and technological relevance. Walking kinematics along straight paths have been experimentally investigated and quantitatively modeled in the diluted limit (i.e., in absence of pedestrian-pedestrian interactions). It is natural to expect that models for straight paths may be an accurate approximations of the dynamics even for paths with curvature radii much larger than the size of a single person. Conversely, as paths curvature increase one may expect larger and larger deviations. As no clear experimental consensus has been reached yet in the literature, here we accurately and systematically investigate the effect of paths curvature on diluted pedestrian dynamics. Thanks to a extensive and highly accurate set of real-life measurements campaign, we derive a Langevin-like social-force model quantitatively compatible with both averages and fluctuations of the walking dynamics. Leveraging on the differential geometric notion of covariant derivative, we generalize previous work by some of the authors, effectively casting a Langevin social-force model for the straight walking dynamics in a curved geometric setting. We deem this the necessary first step to understand and model the more general and ubiquitous case of pedestrians following curved paths in the presence of crowd traffic.

Toward learning Lattice Boltzmann collision operators.

In this work, we explore the possibility of learning from data collision operators for the Lattice Boltzmann Method using a deep learning approach. We compare a hierarchy of designs of the neural network (NN) collision operator and evaluate the performance of the resulting LBM method in reproducing time dynamics of several canonical flows. In the current study, as a first attempt to address the learning problem, the data were generated by a single relaxation time BGK operator. We demonstrate that vanilla NN architecture has very limited accuracy. On the other hand, by embedding physical properties, such as conservation laws and symmetries, it is possible to dramatically increase the accuracy by several orders of magnitude and correctly reproduce the short and long time dynamics of standard fluid flows.

Spatial population genetics with fluid flow.

The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2D), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1D and 2D flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how thesegeneticdrops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.

Self-similar properties of avalanche statistics in a simple turbulent model.

In this paper, we consider a simplified model of turbulence for large Reynolds numbers driven by a constant power energy input on large scales. In the statistical stationary regime, the behaviour of the kinetic energy is characterized by two well-defined phases: a laminar phase where the kinetic energy grows linearly for a (random) time [Formula: see text] followed by abrupt avalanche-like energy drops of sizes [Formula: see text] due to strong intermittent fluctuations of energy dissipation. We study the probability distribution [Formula: see text] and [Formula: see text] which both exhibit a quite well-defined scaling behaviour. Although [Formula: see text] and [Formula: see text] are not statistically correlated, we suggest and numerically checked that their scaling properties are related based on a simple, but non-trivial, scaling argument. We propose that the same approach can be used for other systems showing avalanche-like behaviour such as amorphous solids and seismic events. This article is part of the theme issue 'Scaling the turbulence edifice (part 1)'.

Fluctuations in pedestrian dynamics routing choices.

Routing choices of walking pedestrians in geometrically complex environments are regulated by the interplay of a multitude of factors such as local crowding, (estimated) time to destination, and (perceived) comfort. As individual choices combine, macroscopic traffic flow patterns emerge. Understanding the physical mechanisms yielding macroscopic traffic distributions in environments with complex geometries is an outstanding scientific challenge, with implications in the design and management of crowded pedestrian facilities. In this work, we analyze, by means of extensive real-life pedestrian tracking data, unidirectional flow dynamics in an asymmetric setting, as a prototype for many common complex geometries. Our environment is composed of a main walkway and a slightly longer detour. Our measurements have been collected during a dedicated high-accuracy pedestrian tracking campaign held in Eindhoven (The Netherlands). We show that the dynamics can be quantitatively modeled by introducing a collective discomfort function, and that fluctuations on the behavior of single individuals are crucial to correctly recover the global statistical behavior. Notably, the observed traffic split substantially departs from an optimal, transport-wise, partition, as the global pedestrian throughput is not maximized.

A lattice Boltzmann study of particle settling in a fluctuating multicomponent fluid under confinement.

We present mesoscale numerical simulations based on the coupling of the fluctuating lattice Boltzmann method for multicomponent systems with a wetted finite-size particle model. This newly coupled methodologies are used to study the motion of a spherical particle driven by a constant body force in a confined channel with a fixed square cross section. The channel is filled with a mixture of two liquids under the effect of thermal fluctuations. After some validations steps in the absence of fluctuations, we study the fluctuations in the particle's velocity at changing thermal energy, applied force, particle size, and particle wettability. The importance of fluctuations with respect to the mean settling velocity is quantitatively assessed, especially in comparison with unconfined situations. Results show that the expected effects of confinement are very well captured by the numerical simulations, wherein the confinement strongly enhances the importance of velocity fluctuations, which can be one order of magnitude larger than what expected in unconfined domains. The observed findings underscore the versatility of the proposed methodology in highlighting the effects of confinement on the motion of particles in the presence of thermal fluctuations.

Asymmetric invasion in anisotropic porous media.

We report and discuss, by means of pore-scale numerical simulations, the possibility of achieving a directional-dependent two-phase flow behavior during the process of invasion of a viscous fluid into anisotropic porous media with controlled design. By customising the pore-scale morphology and heterogeneities with the adoption of anisotropic triangular pillars distributed with quenched disorder, we observe a substantially different invasion dynamics according to the direction of fluid injection relative to the medium orientation, that is depending if the triangular pillars have their apex oriented (flow aligned) or opposed (flow opposing) to the main flow direction. Three flow regimes can be observed: (i) for low values of the ratio between the macroscopic pressure drop and the characteristic pore-scale capillary threshold, i.e., for Δp_{0}/p_{c}≤1, the fluid invasion dynamics is strongly impeded and the viscous fluid is unable to reach the outlet of the medium, irrespective of the direction of injection; (ii) for intermediate values, 1<Δp_{0}/p_{c}≤2, the viscous fluid reaches the outlet only when the triangular pillars are flow-opposing oriented; (iii) for larger values, i.e., for Δp_{0}/p_{c}>2, the outlet is again reached irrespective of the direction of injection. The porous medium anisotropy induces a lower effective resistance when the pillars are flow-opposing oriented, suppressing front roughening and capillary fingering. We thus argue that the invasion process occurs as long as the pressure drop is larger then the macroscopic capillary pressure determined by the front roughness, which in the case of flow-opposing pillars is halved. We present a simple approximated model, based on Darcy's assumptions, that links the macroscopic effective permeability with the directional-dependent front roughening, to predict the asymmetric invasion dynamics. This peculiar behavior opens up the possibility of fabrication of porous capillary valves to control the flow along certain specific directions.

Strong noise limit for population dynamics in incompressible advection.

Genetic diversity is at the basis of the evolution process of populations and it is responsible for the populations' degree of fitness to a particular ecosystem. In marine environments many factors play a role in determining the dynamics of a population, including the amount of nutrients, the temperature, and many other stressing factors. An important and yet rather unexplored challenge is to figure out the role of individuals' dispersion, due to flow advection, on population genetics. In this paper we focus on two populations, one of which has a slight selective advantage, advanced by an incompressible two-dimensional flow. In particular, we want to understand how this advective flow can modify the dynamics of the advantageous allele. We generalize, through a theoretical analysis, previous evidence according to which the fixation probability is independent of diffusivity, showing that this is also independent of fluid advection. These findings may have important implications in the understanding of the dynamics of a population of microorganism, such as plankton or bacteria, in marine environments under the influence of (turbulent) currents.

Continuum modeling of shear startup in soft glassy materials.

Yield stress fluids (YSFs) display a dual nature highlighted by the existence of a critical stress σ_{y} such that YSFs are solid for stresses σ imposed below σ_{y}, whereas they flow like liquids for σ>σ_{y}. Under an applied shear rate γ[over ̇], the solid-to-liquid transition is associated with a complex spatiotemporal scenario that depends on the microscopic details of the system, on the boundary conditions, and on the system size. Still, the general phenomenology reported in the literature boils down to a simple sequence that can be divided into a short-time response characterized by the so-called "stress overshoot," followed by stress relaxation towards a steady state. Such relaxation can be either (1) long-lasting, which usually involves the growth of a shear band that can be only transient or that may persist at steady state or (2) abrupt, in which case the solid-to-liquid transition resembles the failure of a brittle material, involving avalanches. In the present paper, we use a continuum model based on a spatially resolved fluidity approach to rationalize the complete scenario associated with the shear-induced yielding of YSFs. A key feature of our model is to provide a scaling for the coordinates of the stress overshoot, i.e., stress σ_{M} and strain γ_{M} as a function of γ[over ̇], which shows good agreement with experimental and numerical data extracted from the literature. Moreover, our approach shows that the power-law scaling σ_{M}(γ[over ̇]) is intimately linked to the growth dynamics of a fluidized boundary layer in the vicinity of the moving boundary. Yet such scaling is independent of the fate of that layer, and of the long-term behavior of the YSF, i.e., whether the steady-state flow profile is homogeneous or shear-banded. Finally, when including the presence of "long-range" correlations, we show that our model displays a ductile to brittle transition, i.e., the stress overshoot reduces into a sharp stress drop associated with avalanches, which impacts the scaling σ_{M}(γ[over ̇]). This generalized model nicely captures subtle avalanche-like features of the transient shear banding dynamics reported in experiments. Our work offers a unified picture of shear-induced yielding in YSFs, whose complex spatiotemporal dynamics are deeply connected to nonlocal effects.

The collective effect of finite-sized inhomogeneities on the spatial spread of populations in two dimensions.

The dynamics of a population expanding into unoccupied habitat has been primarily studied for situations in which growth and dispersal parameters are uniform in space or vary in one dimension. Here, we study the influence of finite-sized individual inhomogeneities and their collective effect on front speed if randomly placed in a two-dimensional habitat. We use an individual-based model to investigate the front dynamics for a region in which dispersal or growth of individuals is reduced to zero (obstacles) or increased above the background (hotspots), respectively. In a regime where front dynamics is determined by a local front speed only, a principle of least time can be employed to predict front speed and shape. The resulting analytical solutions motivate an event-based algorithm illustrating the effects of several obstacles or hotspots. We finally apply the principle of least time to large heterogeneous environments by solving the Eikonal equation numerically. Obstacles lead to a slow-down that is dominated by the number density and width of obstacles, but not by their precise shape. Hotspots result in a speed-up, which we characterize as function of hotspot strength and density. Our findings emphasize the importance of taking the dimensionality of the environment into account.

Stress Overshoots in Simple Yield Stress Fluids.

Soft glassy materials such as mayonnaise, wet clays, or dense microgels display a solid-to-liquid transition under external shear. Such a shear-induced transition is often associated with a nonmonotonic stress response in the form of a stress maximum referred to as "stress overshoot." This ubiquitous phenomenon is characterized by the coordinates of the maximum in terms of stress σ_{M} and strain γ_{M} that both increase as weak power laws of the applied shear rate. Here we rationalize such power-law scalings using a continuum model that predicts two different regimes in the limit of low and high applied shear rates. The corresponding exponents are directly linked to the steady-state rheology and are both associated with the nucleation and growth dynamics of a fluidized region. Our work offers a consistent framework for predicting the transient response of soft glassy materials upon startup of shear from the local flow behavior to the global rheological observables.

Lattice Boltzmann method investigation of a reactive electro-kinetic flow in porous media: towards a phenomenological model.

A model based on the Lattice Boltzmann method is developed to study the flow of reactive electro-kinetic fluids in porous media. The momentum, concentration and electric/potential fields are simulated via the Navier-Stokes, advection-diffusion/Nernst-Planck and Poisson equations, respectively. With this model, the total density and velocity fields, the concentration of reactants and reaction products, including neutral and ionized species, the electric potential and the interaction forces between the fields can be studied, and thus we provide an insight into the interplay between chemistry, flow and the geometry of the porous medium. The results show that the conversion efficiency of the reaction can be strongly influenced by the fluid velocity, reactant concentration and by porosity of the porous medium. The fluid velocity determines how long the reactants stay in the reaction areas, the reactant concentration controls the amount of the reaction material and with different dielectric constant, the porous medium can distort the electric field differently. All these factors make the reaction conversion efficiency display a non-trivial and non-monotonic behaviour as a function of the flow and reaction parameters. To better illustrate the dependence of the reaction conversion efficiency on the control parameters, based on the input from a number of numerical investigations, we developed a phenomenological model of the reactor. This model is capable of capturing the main features of the causal relationship between the performance of the reactor and the main test parameters. Using this model, one could optimize the choice of reaction and flow parameters in order to improve the performance of the reactor and achieve higher production rates. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Impact of the prequench state of binary fluid mixtures on surface-directed spinodal decomposition.

Using lattice Boltzmann simulations we investigate the impact of the amplitude of concentration fluctuations in binary fluid mixtures prior to demixing when in contact with a surface that is preferentially wet by one of the components. We find a bicontinuous structure near the surface for an initial, prequench state of the mixture close to the critical point where the amplitude of concentration fluctuations is large. In contrast, if the initial state of the mixture is not near the critical point and concentration fluctuations are relatively weak, then the morphology is not bicontinuous but remains layered until the very late stages of coarsening. In both cases, it is the morphology of a depletion layer rich in the nonpreferred component that dictates the growth exponent of the thickness of the fluid layer that is in direct contact with the substrate. In the early stages of demixing, we find a growth exponent consistent with a value of 1/4 for a prequench state away from the critical point, which is different from the usual diffusive scaling exponent of 1/3 that we recover for a prequench state close to the critical point. We attribute this to the structure of a depletion layer that is penetrated by tubes of the preferred fluid, connecting the wetting layer to the bulk fluid even in the early stages if the initial state is characterized by concentration fluctuations that are large in amplitude. Furthermore, we find that in the late stages of demixing the flow through these tubes results in significant in-plane concentration variations near the substrate, leading to dropletlike structures with a concentration lower than the average concentration in the wetting layer. This causes a deceleration in the growth of the wetting layer in the very late stages of the demixing. Irrespective of the prequench state of the mixture, the late stages of the demixing process produce the same growth law for the layer thickness, with a scaling exponent of unity usually associated with the impact of hydrodynamic flow fields.

Deep learning velocity signals allow quantifying turbulence intensity.

Turbulence, the ubiquitous and chaotic state of fluid motions, is characterized by strong and statistically nontrivial fluctuations of the velocity field, and it can be quantitatively described only in terms of statistical averages. Strong nonstationarities impede statistical convergence, precluding quantifying turbulence, for example, in terms of turbulence intensity or Reynolds number. Here, we show that by using deep neural networks, we can accurately estimate the Reynolds number within 15% accuracy, from a statistical sample as small as two large-scale eddy turnover times. In contrast, physics-based statistical estimators are limited by the convergence rate of the central limit theorem and provide, for the same statistical sample, at least a hundredfold larger error. Our findings open up previously unexplored perspectives and the possibility to quantitatively define and, therefore, study highly nonstationary turbulent flows as ordinarily found in nature and in industrial processes.