Showing 1–47 of 47 results for author: Benzi, R

Intermittent thermal convection in jammed emulsions

Authors: Francesca Pelusi, Andrea Scagliarini, Mauro Sbragaglia, Massimo Bernaschi, Roberto Benzi

Abstract: We study the process of thermal convection in jammed emulsions with a yield-stress rheology. We find that heat transfer occurs via an intermittent mechanism, whereby intense short-lived convective "heat bursts" are spaced out by long-lasting conductive periods. This behaviour is the result of a sequence of fluidization-rigidity transitions, rooted in a non-trivial interplay between emulsion yield-… ▽ More We study the process of thermal convection in jammed emulsions with a yield-stress rheology. We find that heat transfer occurs via an intermittent mechanism, whereby intense short-lived convective "heat bursts" are spaced out by long-lasting conductive periods. This behaviour is the result of a sequence of fluidization-rigidity transitions, rooted in a non-trivial interplay between emulsion yield-stress rheology and plastic activity, which we characterize via a statistical analysis of the dynamics at the droplet scale. We also show that droplets' coalescence induced during heat bursts leads to a spatially heterogeneous phase-inversion of the emulsion which eventually supports a sustained convective state. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 6 pages, 4 figures

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

Authors: Thibaut Divoux, Elisabeth Agoritsas, Stefano Aime, Catherine Barentin, Jean-Louis Barrat, Roberto Benzi, Ludovic Berthier, Dapeng Bi, Giulio Biroli, Daniel Bonn, Philippe Bourrianne, Mehdi Bouzid, Emanuela Del Gado, Hélène Delanoë-Ayari, Kasra Farain, Suzanne Fielding, Matthias Fuchs, Jasper van der Gucht, Silke Henkes, Maziyar Jalaal, Yogesh M. Joshi, Anaël Lemaître, Robert L. Leheny, Sébastien Manneville, Kirsten Martens , et al. (15 additional authors not shown)

Abstract: 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 worksh… ▽ More 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. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 21 pages, 7 figures, perspective on the workshop 'Yield stress and fluidization in brittle and ductile amorphous systems' held from January 9 to 13, 2023 at the Lorentz Center in Leiden

Lagrangian statistics of dense emulsions

Authors: Ivan Girotto, Andrea Scagliarini, Roberto Benzi, Federico Toschi

Abstract: The dynamics of dense stabilized emulsions presents a rich phenomenology including chaotic emulsification, non-Newtonian rheology and ageing dynamics at rest. Macroscopic rheology results from the complex droplet microdynamics and, in turn, droplet dynamics is influenced by macroscopic flows via the competing action of hydrodynamic and interfacial stresses, giving rise to a complex tangle of elast… ▽ More The dynamics of dense stabilized emulsions presents a rich phenomenology including chaotic emulsification, non-Newtonian rheology and ageing dynamics at rest. Macroscopic rheology results from the complex droplet microdynamics and, in turn, droplet dynamics is influenced by macroscopic flows via the competing action of hydrodynamic and interfacial stresses, giving rise to a complex tangle of elastoplastic effects, diffusion, breakups and coalescence events. This tight multiscale coupling, together with the daunting challenge of experimentally investigating droplets under flow, hindered the understanding of dense emulsions dynamics. We present results from 3D numerical simulations of dense stabilised emulsions, resolving the shape and dynamics of individual droplets, along with the macroscopic flows. We investigate droplet dispersion statistics, measuring probability density functions (PDF) of droplet displacements and velocities, changing the concentration, in the stirred and ageing regimes. We provide the first measurements ever, in concentrated emulsions, of the relative droplet-droplet separations PDF and of the droplet acceleration PDF, which becomes strongly non-Gaussian as the volume fraction is increased above the jamming point. Cooperative effects, arising when droplets are in contact, are argued to be responsible of the anomalous superdiffusive behaviour of the mean square displacement and of the pair separation at long times, in both the stirred and in the ageing regimes. This superdiffusive behaviour is reflected in a non-Gaussian pair separation PDF, whose analytical form is investigated, in the ageing regime, by means of theoretical arguments. This work paves the way to developing a connection between Lagrangian dynamics and rheology in dense stabilised emulsions. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Continuum modeling of Soft Glassy Materials under shear

Authors: Roberto Benzi, Thibaut Divoux, Catherine Barentin, Sébastien Manneville, Mauro Sbragaglia, Federico Toschi

Abstract: Soft Glassy Materials (SGM) consist in dense amorphous assemblies of colloidal particles of multiple shapes, elasticity, and interactions, which confer upon them solid-like properties at rest. They are ubiquitously encountered in modern engineering, including additive manufacturing, semi-solid flow cells, dip-coating, adhesive locomotion, where they are subjected to complex mechanical histories. S… ▽ More Soft Glassy Materials (SGM) consist in dense amorphous assemblies of colloidal particles of multiple shapes, elasticity, and interactions, which confer upon them solid-like properties at rest. They are ubiquitously encountered in modern engineering, including additive manufacturing, semi-solid flow cells, dip-coating, adhesive locomotion, where they are subjected to complex mechanical histories. Such processes often include a solid-to-liquid transition induced by large enough shear, which results in complex transient phenomena such as non-monotonic stress responses, i.e., stress overshoot, and spatially heterogeneous flows, e.g., shear-banding or brittle failure. In the present article, we propose a pedagogical introduction to a continuum model based on a spatially-resolved fluidity approach that we recently introduced to rationalize shear-induced yielding in SGMs. Our model, which relies upon non-local effects, quantitatively captures salient features associated with such complex flows, including the rate dependence of the stress overshoot, as well as transient shear-banded flows together with nontrivial scaling laws for fluidization times. This approach offers a versatile framework to account for subtle effects, such as avalanche-like phenomena, or the impact of boundary conditions, which we illustrate by including in our model the elasto-hydrodynamic slippage of soft particles compressed against solid surfaces. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 8 pages, 4 figures

Self-consistent stochastic dynamics for finite-size networks of spiking neurons

Authors: Gianni V. Vinci, Roberto Benzi, Maurizio Mattia

Abstract: Despite the huge number of neurons composing a brain network, ongoing activity of local cell assemblies composing cortical columns is intrinsically stochastic. Fluctuations in their instantaneous rate of spike firing $ν(t)$ scale with the size of the assembly and persist in isolated network, i.e., in absence of external source of noise. Although deterministic chaos due to the quenched disorder of… ▽ More Despite the huge number of neurons composing a brain network, ongoing activity of local cell assemblies composing cortical columns is intrinsically stochastic. Fluctuations in their instantaneous rate of spike firing $ν(t)$ scale with the size of the assembly and persist in isolated network, i.e., in absence of external source of noise. Although deterministic chaos due to the quenched disorder of the synaptic couplings likely underlies this seemingly stochastic dynamics, an effective theory for the network dynamics of a finite ensemble of spiking neurons is lacking. Here, we fill this gap by extending the so-called population density approach including an activity- and size-dependent stochastic source in the Fokker-Planck equation for the membrane potential density. The finite-size noise embedded in this stochastic partial derivative equation is analytically characterized leading to a self-consistent and non-perturbative description of $ν(t)$ valid for a wide class of spiking neuron networks. Its power spectra of $ν(t)$ are found in excellent agreement with those from detailed simulations both in the linear regime and across a synchronization phase transition, when a size-dependent smearing of the critical dynamics emerges. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 6 pages, 4 figures

Spatial population genetics with fluid flow

Authors: Roberto Benzi, David R. Nelson, Suraj Shankar, Federico Toschi, Xiaojue Zhu

Abstract: 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 m… ▽ More 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 these genetic drops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well. △ Less

Submitted 30 June, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 29 pages, 22 figures

Journal ref: Rep. Prog. Phys. 85 096601, 2022

Build up of yield stress fluids via chaotic emulsification

Authors: Ivan Girotto, Roberto Benzi, Gianluca Di Staso, Andrea Scagliarini, Sebastiano Fabio Schifano, Federico Toschi

Abstract: Stabilized dense emulsions display a rich phenomenology connecting microstructure and rheology. In this work we study how an emulsion with a finite yield stress can be built via large-scale stirring. By gradually increasing the volume fraction of the dispersed minority phase, under the constant action of a stirring force, we are able to achieve volume fractions close to $80\%$. Despite the fact th… ▽ More Stabilized dense emulsions display a rich phenomenology connecting microstructure and rheology. In this work we study how an emulsion with a finite yield stress can be built via large-scale stirring. By gradually increasing the volume fraction of the dispersed minority phase, under the constant action of a stirring force, we are able to achieve volume fractions close to $80\%$. Despite the fact that our system is highly concentrated and not yet turbulent we observe a droplet size distribution consistent with the $-10/3$ scaling, often associated to inertial range droplets breakup. We report that the polydispersity of droplet sizes correlates with the dynamics of the emulsion formation process. Additionally we quantify the visco-elastic properties of the dense emulsion finally obtained and we demonstrate the presence of a finite yield stress. The approach reported can pave the way to a quantitative understanding of the complex interplay between the dynamics of mesoscale constituents and the large scale flow properties yield-stress fluids. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Self Similar Properties of Avalanche Statistics in a Simple Turbulent Model

Authors: Roberto Benzi, Ilaria Castaldi, Federico Toschi, Jeannot Trampert

Abstract: 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 characterised by two well defined phases: a laminar phase where the kinetic energy grows linearly for a (random) time $t_w$ followed by abrupt avalanche-like energy drops of sizes… ▽ More 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 characterised by two well defined phases: a laminar phase where the kinetic energy grows linearly for a (random) time $t_w$ followed by abrupt avalanche-like energy drops of sizes $S$ due to strong intermittent fluctuations of energy dissipation. We study the probability distribution $P[t_w]$ and $P[S]$ which both exhibit a quite well defined scaling behaviour. Although $t_w$ and $S$ 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. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: 12 pages, 10 figures

Continuum modelling of shear start-up in soft glassy materials

Authors: Roberto Benzi, Thibaut Divoux, Catherine Barentin, Sébastien Manneville, Mauro Sbragaglia, Federico Toschi

Abstract: Yield stress fluids (YSFs) display a dual nature highlighted by the existence of a yield stress such that YSFs are solid below the yield stress, whereas they flow like liquids above it. Under an applied shear rate $\dotγ$, the solid-to-liquid transition is associated with a complex spatiotemporal scenario. Still, the general phenomenology reported in the literature boils down to a simple sequence… ▽ More Yield stress fluids (YSFs) display a dual nature highlighted by the existence of a yield stress such that YSFs are solid below the yield stress, whereas they flow like liquids above it. Under an applied shear rate $\dotγ$, the solid-to-liquid transition is associated with a complex spatiotemporal scenario. 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 long-lasting, which usually involves the growth of a shear band that can be only transient or that may persist at steady-state, or abrupt, in which case the solid-to-liquid transition resembles the failure of a brittle material, involving avalanches. Here 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. Our model provides a scaling for the coordinates of the stress maximum as a function of $\dotγ$, which shows excellent agreement with experimental and numerical data extracted from the literature. Moreover, our approach shows that such a scaling 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. 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 of the stress maximum with $\dotγ$. Our work offers a unified picture of shear-induced yielding in YSFs, whose complex spatiotemporal dynamics are deeply connected to non-local effects. △ Less

Submitted 1 July, 2021; v1 submitted 31 March, 2021; originally announced March 2021.

Comments: 24 pages, 19 figures

Journal ref: Phys. Rev. E 104, 034612 (2021)

Stress Overshoots in Simple Yield Stress Fluids

Authors: Roberto Benzi, Thibaut Divoux, Catherine Barentin, Sébastien Manneville, Mauro Sbragaglia, Federico Toschi

Abstract: Soft glassy materials such as mayonnaise, wet clays, or dense microgels display under external shear a solid-to-liquid transition. Such a shear-induced transition is often associated with a non-monotonic 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 $σ_\text{M}$ an… ▽ More Soft glassy materials such as mayonnaise, wet clays, or dense microgels display under external shear a solid-to-liquid transition. Such a shear-induced transition is often associated with a non-monotonic 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 $σ_\text{M}$ and strain $γ_\text{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 start-up of shear from the local flow behavior to the global rheological observables. △ Less

Submitted 13 September, 2021; v1 submitted 31 March, 2021; originally announced March 2021.

Comments: 5 pages, 4 figures, supplemental materials with 2 figures

Journal ref: Phys. Rev. Lett. 127, 148003 (2021)

Rayleigh-Bénard convection of a model emulsion: anomalous heat-flux fluctuations and finite-size droplets effects

Authors: Francesca Pelusi, Mauro Sbragaglia, Roberto Benzi, Andrea Scagliarini, Massimo Bernaschi, Sauro Succi

Abstract: We present mesoscale numerical simulations of Rayleigh-Bénard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration Phi_0 is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from… ▽ More We present mesoscale numerical simulations of Rayleigh-Bénard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration Phi_0 is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from conductive to convective states, where it is known that a homogeneous Newtonian system exhibits a steady flow and a time-independent heat flux. In marked contrast, emulsions exhibit a non-steady dynamics with fluctuations in the heat flux. In this paper, we aim at the characterisation of such non-steady dynamics via detailed studies on the time-averaged heat flux and its fluctuations. To understand the time-averaged heat flux, we propose a side-by-side comparison between the emulsion system and a single-phase (SP) system, whose viscosity is constructed from the shear rheology of the emulsion. We show that such local closure works well only when a suitable degree of coarse-graining (at the droplet scale) is introduced in the local viscosity. To delve deeper into the fluctuations in the heat flux, we propose a side-by-side comparison between a Newtonian emulsion and a non-Newtonian emulsion, at fixed time-averaged heat flux. This comparison elucidates that finite-size droplets and the non-Newtonian rheology cooperate to trigger enhanced heat-flux fluctuations at the droplet scales. These enhanced fluctuations are rooted in the emergence of space correlations among distant droplets, which we highlight via direct measurements of the droplets displacement and the characterisation of the associated correlation function. The observed findings offer insights on heat transfer properties for confined systems possessing finite-size constituents. △ Less

Submitted 19 February, 2021; v1 submitted 9 January, 2021; originally announced January 2021.

Comments: 15 pages, 22 figures. arXiv admin note: substantial text overlap with arXiv:2010.02139, arXiv:2010.02138

Report number: 17, 3709 - 3721

Journal ref: Soft Matter 2021

Convective heat transfer of a model emulsion at the droplet scale

Authors: Francesca Pelusi, Mauro Sbragaglia, Roberto Benzi, Andrea Scagliarini, Massimo Bernaschi, Sauro Succi

Abstract: We numerically study the Rayleigh-Bénard (RB) convection in two-dimensional model emulsions confined between two parallel walls at fixed temperatures. The systems under study are heterogeneous, with finite-size droplets dispersed in a continuous phase. The droplet concentration is chosen to explore the convective heat transfer of both Newtonian (low droplet concentration) and non-Newtonian (high d… ▽ More We numerically study the Rayleigh-Bénard (RB) convection in two-dimensional model emulsions confined between two parallel walls at fixed temperatures. The systems under study are heterogeneous, with finite-size droplets dispersed in a continuous phase. The droplet concentration is chosen to explore the convective heat transfer of both Newtonian (low droplet concentration) and non-Newtonian (high droplet concentration) emulsions, the latter exhibiting shear-thinning rheology, with a noticeable increase of viscosity at low shear rates. It is well known that the transition to convection of a homogeneous Newtonian system is accompanied by the onset of steady flow and time-independent heat flux; in marked contrast, the heterogeneity of emulsions brings in an additional and previously unexplored phenomenology. As a matter of fact, when the droplet concentration increases, we observe that the heat transfer process is mediated by a non-steady flow, with neat heat-flux fluctuations, obeying a non-Gaussian statistics. The observed findings are ascribed to the emergence of space correlations among distant droplets, which we highlight via direct measurements of the droplets displacement and the characterisation of the associated correlation functions. △ Less

Submitted 16 April, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: This contribution was merged with the arXiv:2010.02138. The interested reader is referred to the merged contribution arXiv:2101.03374

Effective modelling of the Rayleigh-Bénard convection of concentrated emulsions with finite-size droplets

Authors: Francesca Pelusi, Mauro Sbragaglia, Roberto Benzi, Andrea Scagliarini, Massimo Bernaschi, Sauro Succi

Abstract: We present mesoscale numerical simulations of Rayleigh-Bénard convection in a two-dimensional concentrated emulsion, confined between two parallel walls, heated from below and cooled from above, under the effect of buoyancy forces. The systems under study comprise finite-size droplets, whose concentration $Φ_0$ is varied, ranging from the dilute limit up to the point where the emulsion starts to b… ▽ More We present mesoscale numerical simulations of Rayleigh-Bénard convection in a two-dimensional concentrated emulsion, confined between two parallel walls, heated from below and cooled from above, under the effect of buoyancy forces. The systems under study comprise finite-size droplets, whose concentration $Φ_0$ is varied, ranging from the dilute limit up to the point where the emulsion starts to be packed and exhibits non-Newtonian rheology. We focus on the characterisation of the convective heat transfer properties close to the transition from conductive to convective states. The convective flow is confined and heterogeneous, which causes the emulsion to exhibit concentration heterogeneities in space $φ_0(y)$, depending on the location in the wall-to-wall direction ($y$). With the aim of assessing quantitatively the heat transfer efficiency of such heterogeneous systems, we resort to a side-by-side comparison between the concentrated emulsion system and a single-phase (SP) system, whose local viscosity $η^{\mbox{SP}}(y)$ is suitably constructed from the shear rheology of the emulsion. Such comparison highlights that a suitable degree $Λ$ of coarse-graining needs to be introduced in the local viscosity $η_Λ^{\mbox{SP}}(y)$, in order for the single-phase system to attain the same heat transfer efficiency of the emulsion. Specifically, it is shown that a quantitative matching between the two systems is possible whenever the coarse-graining is performed over a scale of the order of the droplet size. △ Less

Submitted 16 April, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: This contribution was merged with the arXiv:2010.02139. The interested reader is referred to the merged contribution arXiv:2101.03374

Deep learning velocity signals allows to quantify turbulence intensity

Authors: Alessandro Corbetta, Vlado Menkovski, Roberto Benzi, Federico Toschi

Abstract: Turbulence, the ubiquitous and chaotic state of fluid motions, is characterized by strong and statistically non-trivial fluctuations of the velocity field, over a wide range of length- and time-scales, and it can be quantitatively described only in terms of statistical averages. Strong non-stationarities hinder the possibility to achieve statistical convergence, making it impossible to define the… ▽ More Turbulence, the ubiquitous and chaotic state of fluid motions, is characterized by strong and statistically non-trivial fluctuations of the velocity field, over a wide range of length- and time-scales, and it can be quantitatively described only in terms of statistical averages. Strong non-stationarities hinder the possibility to achieve statistical convergence, making it impossible to define the turbulence intensity and, in particular, its basic dimensionless estimator, the Reynolds number. Here we show that by employing Deep Neural Networks (DNN) 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 rate of convergence of the central limit theorem, and provide, for the same statistical sample, an error at least $100$ times larger. Our findings open up new perspectives in the possibility to quantitatively define and, therefore, study highly non-stationary turbulent flows as ordinarily found in nature as well as in industrial processes. △ Less

Submitted 14 November, 2019; v1 submitted 13 November, 2019; originally announced November 2019.

Journal ref: Science Advances 7, eaba7281, 2021

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

Authors: Wolfram Möbius, Francesca Tesser, Kim M. J. Alards, Roberto Benzi, David R. Nelson, Federico Toschi

Abstract: 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 investig… ▽ More 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 speedup, which we characterise as function of hotspot strength and density. Our findings emphasise the importance of taking the dimensionality of the environment into account. △ Less

Submitted 20 September, 2021; v1 submitted 11 October, 2019; originally announced October 2019.

Comments: Authors Wolfram Möbius and Francesca Tesser contributed equally

Unified theoretical and experimental view on transient shear banding

Authors: Roberto Benzi, Thibaut Divoux, Catherine Barentin, Sébastien Manneville, Mauro Sbragaglia, Federico Toschi

Abstract: Dense emulsions, colloidal gels, microgels, and foams all display a solid-like behavior at rest characterized by a yield stress, above which the material flows like a liquid. Such a fluidization transition often consists of long-lasting transient flows that involve shear-banded velocity profiles. The characteristic time for full fluidization, $τ_\text{f}$, has been reported to decay as a power-law… ▽ More Dense emulsions, colloidal gels, microgels, and foams all display a solid-like behavior at rest characterized by a yield stress, above which the material flows like a liquid. Such a fluidization transition often consists of long-lasting transient flows that involve shear-banded velocity profiles. The characteristic time for full fluidization, $τ_\text{f}$, has been reported to decay as a power-law of the shear rate $\dot γ$ and of the shear stress $σ$ with respective exponents $α$ and $β$. Strikingly, the ratio of these exponents was empirically observed to coincide with the exponent of the Herschel-Bulkley law that describes the steady-state flow behavior of these complex fluids. Here we introduce a continuum model, based on the minimization of a "free energy", that captures quantitatively all the salient features associated with such \textit{transient} shear-banding. More generally, our results provide a unified theoretical framework for describing the yielding transition and the steady-state flow properties of yield stress fluids. △ Less

Submitted 2 October, 2019; v1 submitted 20 July, 2019; originally announced July 2019.

Comments: 5 pages, 4 figures - supplemental 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 123, 248001 (2019)

Interevent time distributions of avalanche dynamics

Authors: Pinaki Kumar, Evangelos Korkolis, Roberto Benzi, Dmitry Denisov, Andre Niemeijer, Peter Schall, Federico Toschi, Jeannot Trampert

Abstract: Physical systems characterized by stick-slip dynamics often display avalanches. Regardless of the diversity of their microscopic structure, these systems are governed by a power-law distribution of avalanche size and duration. Here we focus on the interevent times between avalanches and show that, unlike their distributions of size and duration, the interevent time distributions are able to distin… ▽ More Physical systems characterized by stick-slip dynamics often display avalanches. Regardless of the diversity of their microscopic structure, these systems are governed by a power-law distribution of avalanche size and duration. Here we focus on the interevent times between avalanches and show that, unlike their distributions of size and duration, the interevent time distributions are able to distinguish different mechanical states of the system, characterized by different volume fractions or confining pressures. We use experiments on granular systems and numerical simulations of emulsions to show that systems having the same probability distribution for avalanche size and duration can have different interevent time distributions. Remarkably, for large packing ratios, these interevent time distributions look similar to those for earthquakes and are indirect evidence of large space-time correlations in the system. Our results therefore indicate that interevent time statistics are more informative to characterize the dynamics of avalanches. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 14 pages, 11 figures

Avalanche statistics during coarsening dynamics

Authors: Francesca Pelusi, Mauro Sbragaglia, Roberto Benzi

Abstract: We study the coarsening dynamics of a two dimensional system via lattice Boltzmann numerical simulations. The system under consideration is a biphasic system consisting of domains of a dispersed phase closely packed together in a continuous phase and separated by thin interfaces. Such system is elastic and typically out of equilibrium. The equilibrium state is attained via the coarsening dynamics,… ▽ More We study the coarsening dynamics of a two dimensional system via lattice Boltzmann numerical simulations. The system under consideration is a biphasic system consisting of domains of a dispersed phase closely packed together in a continuous phase and separated by thin interfaces. Such system is elastic and typically out of equilibrium. The equilibrium state is attained via the coarsening dynamics, wherein the dispersed phase slowly diffuses through the interfaces, causing domains to change in size and eventually rearrange abruptly. The effect of rearrangements is propagated throughout the system via the intrinsic elastic interactions and may cause rearrangements elsewhere, resulting in intermittent bursts of activity and avalanche behaviour. Here we aim at quantitatively characterizing the corresponding avalanche statistics (i.e. size, duration, inter-avalanche time). Despite the coarsening dynamics is triggered by an internal driving mechanism, we find quantitative indications that such avalanche statistics displays scaling-laws very similar to those observed in the response of disordered materials to external loads. △ Less

Submitted 5 April, 2019; originally announced April 2019.

Report number: 45, 4518

Journal ref: Soft Matter 2019

Polymers in Fluid Flows

Authors: Roberto Benzi, Emily S. C. Ching

Abstract: The interaction of flexible polymers with fluid flows leads to a number of intriguing phenomena observed in laboratory experiments, namely drag reduction, elastic turbulence and heat transport modification in natural convection, and is one of the most challenging subjects in soft matter physics. In this paper we review our present knowledge on the subject. Our present knowledge is mostly based on… ▽ More The interaction of flexible polymers with fluid flows leads to a number of intriguing phenomena observed in laboratory experiments, namely drag reduction, elastic turbulence and heat transport modification in natural convection, and is one of the most challenging subjects in soft matter physics. In this paper we review our present knowledge on the subject. Our present knowledge is mostly based on direct numerical simulations performed in the last twenty years, which have successfully explained, at least qualitatively, most of the experimental results. Our goal is to disentangle as much as possible the basic mechanisms acting in the system in order to capture the basic features underlying different theoretical approaches and explanations. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Journal ref: Annu. Rev. Condens. Matter Phys. 2018. 9:163-181

Fixation Probabilities in Weakly Compressible Fluid Flows

Authors: Abigail Plummer, Roberto Benzi, David R. Nelson, Federico Toschi

Abstract: Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we co… ▽ More Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we couple an off-lattice agent-based simulation of two populations in 1D to a weakly compressible velocity field--first a sine wave and then a shell model of turbulence. We find for both cases that even in a regime where the overall population structure is approximately unaltered, the flow can significantly diminish the effect of a selective advantage on fixation probabilities. We understand this effect in terms of the enhanced survival of organisms born at sources in the flow and the influence of Fisher genetic waves. △ Less

Submitted 23 January, 2019; v1 submitted 21 August, 2018; originally announced August 2018.

Comments: 8 pages, 6 figures; SI: 9 pages, 3 figures

Journal ref: PNAS, vol. 116 no. 2, 373-378 (2019)

Metastability at the Yield-Stress Transition in Soft Glasses

Authors: Matteo Lulli, Roberto Benzi, Mauro Sbragaglia

Abstract: We study the solid-to-liquid transition in a two-dimensional fully periodic soft-glassy model with an imposed spatially heterogeneous stress. The model we consider consists of droplets of a dispersed phase jammed together in a continuous phase. When the peak value of the stress gets close to the yield stress of the material, we find that the whole system intermittently tunnels to a metastable "flu… ▽ More We study the solid-to-liquid transition in a two-dimensional fully periodic soft-glassy model with an imposed spatially heterogeneous stress. The model we consider consists of droplets of a dispersed phase jammed together in a continuous phase. When the peak value of the stress gets close to the yield stress of the material, we find that the whole system intermittently tunnels to a metastable "fluidized" state, which relaxes back to a metastable "solid" state by means of an elastic-wave dissipation. This macroscopic scenario is studied through the microscopic displacement field of the droplets, whose time statistics displays a remarkable bimodality. Metastability is rooted in the existence, in a given stress range, of two distinct stable rheological branches as well as long-range correlations (e.g., large dynamic heterogeneity) developed in the system. Finally, we show that a similar behavior holds for a pressure-driven flow, thus suggesting possible experimental tests. △ Less

Submitted 5 June, 2018; v1 submitted 2 October, 2017; originally announced October 2017.

Comments: 13 pages, 11 figures

Journal ref: Phys. Rev. X 8, 021031 (2018)

Cooperativity flows and Shear-Bandings: a statistical field theory approach

Authors: R. Benzi, M. Sbragaglia, M. Bernaschi, S. Succi, F. Toschi

Abstract: Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin \& A. Ajdari ({\em Phys. Rev. Lett.} {\bf 103}, 036001 (2009)), we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity), are intimately related to the eme… ▽ More Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin \& A. Ajdari ({\em Phys. Rev. Lett.} {\bf 103}, 036001 (2009)), we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity), are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. Compactons coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material. △ Less

Submitted 14 October, 2015; v1 submitted 27 July, 2015; originally announced July 2015.

Comments: 33 pages, 10 figures

Homogeneous and Isotropic Turbulence: a short survey on recent developments

Authors: R. Benzi, L. Biferale

Abstract: We present a detailed review of some of the most recent developments on Eulerian and Lagrangian turbulence in homogeneous and isotropic statistics. In particular, we review phenomenological and numerical results concerning the issue of universality with respect to the large scale forcing and the viscous dissipative physics. We discuss the state-of-the-art of numerical versus experimental compariso… ▽ More We present a detailed review of some of the most recent developments on Eulerian and Lagrangian turbulence in homogeneous and isotropic statistics. In particular, we review phenomenological and numerical results concerning the issue of universality with respect to the large scale forcing and the viscous dissipative physics. We discuss the state-of-the-art of numerical versus experimental comparisons and we discuss the dicotomy between phenomenology based on coherent structures or on statistical approaches. A detailed discussion of finite Reynolds effects is also presented. △ Less

Submitted 30 September, 2015; v1 submitted 4 January, 2015; originally announced January 2015.

Comments: based on the talk presented by R. Benzi at DSFD 2-14. postprint version, published online on 6 July 2015 J. Stat. Phys

Internal dynamics and activated processes in Soft-Glassy materials

Authors: R. Benzi, M. Sbragaglia, A. Scagliarini, P. Perlekar, M. Bernaschi, S. Succi, F. Toschi

Abstract: Plastic rearrangements play a crucial role in the characterization of soft-glassy materials, such as emulsions and foams. Based on numerical simulations of soft-glassy systems, we study the dynamics of plastic rearrangements at the hydrodynamic scales where thermal fluctuations can be neglected. Plastic rearrangements require an energy input, which can be either provided by external sources, or ma… ▽ More Plastic rearrangements play a crucial role in the characterization of soft-glassy materials, such as emulsions and foams. Based on numerical simulations of soft-glassy systems, we study the dynamics of plastic rearrangements at the hydrodynamic scales where thermal fluctuations can be neglected. Plastic rearrangements require an energy input, which can be either provided by external sources, or made available through time evolution in the coarsening dynamics, in which the total interfacial area decreases as a consequence of the slow evolution of the dispersed phase from smaller to large droplets/bubbles. We first demonstrate that our hydrodynamic model can quantitatively reproduce such coarsening dynamics. Then, considering periodically oscillating strains, we characterize the number of plastic rearrangements as a function of the external energy-supply, and show that they can be regarded as activated processes induced by a suitable "noise" effect. Here we use the word noise in a broad sense, referring to the internal non-equilibrium dynamics triggered by spatial random heterogeneities and coarsening. Finally, by exploring the interplay between the internal characteristic time-scale of the coarsening dynamics and the external time-scale associated with the imposed oscillating strain, we show that the system exhibits the phenomenon of stochastic resonance, thereby providing further credit to the mechanical activation scenario. △ Less

Submitted 24 October, 2014; v1 submitted 10 March, 2014; originally announced March 2014.

Comments: 21 Pages, 9 figures

The emergence of supramolecular forces from lattice kinetic models of non ideal fluids: applications to the rheology of soft glassy meterials

Authors: M. Sbragaglia, R. Benzi, M. Bernaschi S. Succi

Abstract: A systematic study for a single-specie lattice Boltzmann model with frustrated-short range attractive and mid/long-range repulsive-interactions is presented. The equilibrium analysis is performed along the guidelines proposed by [X. Shan, Phys. Rev. E 77, 066702 (2008)] and allows us to determine the surface tension density and the resulting disjoining pressure developing in a thin film when two i… ▽ More A systematic study for a single-specie lattice Boltzmann model with frustrated-short range attractive and mid/long-range repulsive-interactions is presented. The equilibrium analysis is performed along the guidelines proposed by [X. Shan, Phys. Rev. E 77, 066702 (2008)] and allows us to determine the surface tension density and the resulting disjoining pressure developing in a thin film when two interfaces overlap. Numerical simulations of confined flows are then performed with a multicomponent model and are successfully tested against the recent suggestion by Bocquet and coworkers on the existence of a cooperative length underlying the non-local rheology of highly confined soft-glassy materials [Goyon et al., Nature 454, 84-87 (2008); Soft Matter 6, 2668-2678 (2010)]. △ Less

Submitted 24 February, 2014; originally announced February 2014.

Journal ref: Soft Matter 8, 10773-10782 (2012)

Direct evidence of plastic events and dynamic heterogeneities in soft-glasses

Authors: R. Benzi, M. Sbragaglia, P. Perlekar, M. Bernaschi, S. Succi, F. Toschi

Abstract: By using fluid-kinetic simulations of confined and concentrated emulsion droplets, we investigate the nature of space non-homogeneity in soft-glassy dynamics and provide quantitative measurements of the statistical features of plastic events in the proximity of the yield-stress threshold. Above the yield stress, our results show the existence of a finite stress correlation scale, which can be mapp… ▽ More By using fluid-kinetic simulations of confined and concentrated emulsion droplets, we investigate the nature of space non-homogeneity in soft-glassy dynamics and provide quantitative measurements of the statistical features of plastic events in the proximity of the yield-stress threshold. Above the yield stress, our results show the existence of a finite stress correlation scale, which can be mapped directly onto the {\it cooperativity scale}, recently introduced in the literature to capture non-local effects in the soft-glassy dynamics. In this regime, the emergence of a separate boundary (wall) rheology with higher fluidity than the bulk, is highlighted in terms of near-wall spontaneous segregation of plastic events. Near the yield stress, where the cooperative scale cannot be estimated with sufficient accuracy, the system shows a clear increase of the stress correlation scale, whereas plastic events exhibit intermittent clustering in time, with no preferential spatial location. A quantitative measurement of the space-time correlation associated with the motion of the interface of the droplets is key to spot the long-range amorphous order at the yield stress threshold. △ Less

Submitted 18 February, 2014; originally announced February 2014.

Selective advantage of diffusing faster

Authors: Simone Pigolotti, Roberto Benzi

Abstract: We study a stochastic spatial model of biological competition in which two species have the same birth and death rates, but different diffusion constants. In the absence of this difference, the model can be considered as an off-lattice version of the Voter model and presents similar coarsening properties. We show that even a relative difference in diffusivity on the order of a few percent may lead… ▽ More We study a stochastic spatial model of biological competition in which two species have the same birth and death rates, but different diffusion constants. In the absence of this difference, the model can be considered as an off-lattice version of the Voter model and presents similar coarsening properties. We show that even a relative difference in diffusivity on the order of a few percent may lead to a strong bias in the coarsening process favoring the more agile species. We theoretically quantify this selective advantage and present analytical formulas for the average growth of the fastest species and its fixation probability. △ Less

Submitted 16 May, 2014; v1 submitted 31 July, 2013; originally announced July 2013.

Comments: 8 pages, 5 figures (Main Text + Supplementary Information). Accepted version

Journal ref: Phys. Rev. Lett. 112, 188102 (2014)

Rheological properties of soft-glassy flows from hydro-kinetic simulations

Authors: R. Benzi, M. Bernaschi, M. Sbragaglia, S. Succi

Abstract: Based on numerical simulations of a lattice kinetic model for soft-glassy materials, we characterize the global rheology of a dense emulsion-like system, under three representative load conditions: Couette flow, time-oscillating Strain and Kolmogorov flow. It is found that in all cases the rheology is described by a Herschel-Bulkley (HB) relation, $σ= σ_{Y} + A S^β$, with the yield stress $σ_{Y}$… ▽ More Based on numerical simulations of a lattice kinetic model for soft-glassy materials, we characterize the global rheology of a dense emulsion-like system, under three representative load conditions: Couette flow, time-oscillating Strain and Kolmogorov flow. It is found that in all cases the rheology is described by a Herschel-Bulkley (HB) relation, $σ= σ_{Y} + A S^β$, with the yield stress $σ_{Y}$ largely independent of the loading scenario. A proper rescaling of the HB parameters permits to describe heterogeneous flows with space-dependent stresses, based on the notion of cooperativity, as recently proposed to characterize the degree of non-locality of stress relaxation phenomena in soft-glassy materials. △ Less

Submitted 28 February, 2014; v1 submitted 12 May, 2013; originally announced May 2013.

Comments: 6 pages, 12 Figures

Journal ref: Europhys. Lett. 104, 48006 (2013)

Growth, competition and cooperation in spatial population genetics

Authors: Simone Pigolotti, Roberto Benzi, Prasad Perlekar Mogens H. Jensen, Federico Toschi, David R. Nelson

Abstract: We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we st… ▽ More We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we study the behavior of fixation probabilities, fixation times, and heterozygosity, in a neutral setting and in cases where the two species can compete or cooperate. By concluding with examples in which individuals are transported by fluid flows, we argue that this model is a natural choice to describe competition in marine environments. △ Less

Submitted 18 December, 2012; v1 submitted 24 August, 2012; originally announced August 2012.

Comments: 29 pages, 14 figures; revised version including a section with results in the presence of fluid flows

Journal ref: Theoretical Population Biology 84, 72-86 (2013)

Dynamic crossover length in soft-glassy materials

Authors: R. Benzi, M. Bernaschi, M. Sbragaglia, S. Succi

Abstract: Based upon mesoscale simulations of binary mixtures with very low surface tension and positive disjoining pressure (frustration), we measure the correlation length of the stress field within the flowing mixture, as a function of the frequency of the applied load. Two scaling regimes are clearly identified, with a sharp crossover between the liquid and solid regions, with exponent -1/2 and -1, resp… ▽ More Based upon mesoscale simulations of binary mixtures with very low surface tension and positive disjoining pressure (frustration), we measure the correlation length of the stress field within the flowing mixture, as a function of the frequency of the applied load. Two scaling regimes are clearly identified, with a sharp crossover between the liquid and solid regions, with exponent -1/2 and -1, respectively. The crossover correlation length is shown to be closely connected to the typical scale of the coherent excitations of the density field, whose collective dynamics is shown to be responsible for non-trivial rheological effects including the emergence of yield stress. △ Less

Submitted 21 January, 2013; v1 submitted 1 August, 2012; originally announced August 2012.

Comments: This paper has been withdrawn by the authors because we are working on a new paper that provides a better description of the same concept

Population dynamics in compressible flows

Authors: Roberto Benzi, Mogens H. Jensen, David R. Nelson, Prasad Perlekar, Simone Pigolotti, Federico Toschi

Abstract: Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible tu… ▽ More Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another △ Less

Submitted 28 March, 2012; originally announced March 2012.

Comments: 16 pages, talk delivered at the Geilo Winter School 2011

Population genetics in compressible flows

Authors: Simone Pigolotti, Roberto Benzi, Mogens H. Jensen, David R. Nelson

Abstract: We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We then study three examples of compressibl… ▽ More We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We then study three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. Despite localization on convergence zones, one species goes extinct much more rapidly than in well-mixed populations. For a weak harmonic potential, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single boundary, whose location depends on the fitness advantage. △ Less

Submitted 17 June, 2011; originally announced June 2011.

Comments: 10 pages, 5 figures, submitted

Journal ref: Phys. Rev. Lett. 108, 128102, 2012

Phase-field model of long-time glass-like relaxation in binary fluid mixtures

Authors: R. Benzi, M. Sbragaglia, M. Bernaschi, S. Succi

Abstract: We present a new phase-field model for binary fluids exhibiting typical signatures of self-glassiness, such as long-time relaxation, ageing and long-term dynamical arrest. The present model allows the cost of building an interface to become locally zero, while preserving global positivity of the overall surface tension. An important consequence of this property, which we prove analytically, is the… ▽ More We present a new phase-field model for binary fluids exhibiting typical signatures of self-glassiness, such as long-time relaxation, ageing and long-term dynamical arrest. The present model allows the cost of building an interface to become locally zero, while preserving global positivity of the overall surface tension. An important consequence of this property, which we prove analytically, is the emergence of compact configurations of fluid density. Owing to their finite-size support, these "compactons" can be arbitrarily superposed, thereby providing a direct link between the ruggedness of the free-energy landscape and morphological complexity in configurational space. The analytical picture is supported by numerical simulations of the proposed phase-field equation. △ Less

Submitted 27 February, 2011; v1 submitted 3 November, 2010; originally announced November 2010.

Comments: 5 Pages, 6 Figures

Shear Banding from lattice kinetic models with competing interactions

Authors: Roberto Benzi, Mauro Sbragaglia, Massimo Bernaschi, Sauro Succi

Abstract: Soft Glassy Materials, Non Linear Rheology, Lattice Kinetic models, frustrated phase separation} We present numerical simulations based on a Boltzmann kinetic model with competing interactions, aimed at characterizating the rheological properties of soft-glassy materials. The lattice kinetic model is shown to reproduce typical signatures of driven soft-glassy flows in confined geometries, such as… ▽ More Soft Glassy Materials, Non Linear Rheology, Lattice Kinetic models, frustrated phase separation} We present numerical simulations based on a Boltzmann kinetic model with competing interactions, aimed at characterizating the rheological properties of soft-glassy materials. The lattice kinetic model is shown to reproduce typical signatures of driven soft-glassy flows in confined geometries, such as Herschel-Bulkley rheology, shear-banding and histeresys. This lends further credit to the present lattice kinetic model as a valuable tool for the theoretical/computational investigation of the rheology of driven soft-glassy materials under confinement. △ Less

Submitted 28 September, 2010; originally announced September 2010.

Comments: 8 Pages, 5 Figures

Herschel-Bulkley rheology from lattice kinetic theory of soft-glassy materials

Authors: R. Benzi, M. Bernaschi, M. Sbragaglia, S. Succi

Abstract: We provide a clear evidence that a two species mesoscopic Lattice Boltzmann (LB) model with competing short-range attractive and mid-range repulsive interactions supports emergent Herschel-Bulkley (HB) rheology, i.e. a power-law dependence of the shear-stress as a function of the strain rate, beyond a given yield-stress threshold. This kinetic formulation supports a seamless transition from flowi… ▽ More We provide a clear evidence that a two species mesoscopic Lattice Boltzmann (LB) model with competing short-range attractive and mid-range repulsive interactions supports emergent Herschel-Bulkley (HB) rheology, i.e. a power-law dependence of the shear-stress as a function of the strain rate, beyond a given yield-stress threshold. This kinetic formulation supports a seamless transition from flowing to non-flowing behaviour, through a smooth tuning of the parameters governing the mesoscopic interactions between the two species. The present model may become a valuable computational tool for the investigation of the rheology of soft-glassy materials on scales of experimental interest. △ Less

Submitted 28 April, 2010; originally announced April 2010.

Comments: 5 figures

Mesoscopic lattice Boltzmann modeling of soft-glassy systems: theory and simulations

Authors: R. Benzi, M. Sbragaglia, S. Succi, M. Bernaschi, S. Chibbaro

Abstract: A multi-component lattice Boltzmann model recently introduced (R. Benzi et al. Phys. Rev. Lett 102, 026002 (2009)) to describe some dynamical behaviors of soft-flowing materials is theoretically analyzed. Equilibrium and transport properties are derived within the framework of a continuum free-energy formulation, and checked against numerical simulations. Due to the competition between short-ran… ▽ More A multi-component lattice Boltzmann model recently introduced (R. Benzi et al. Phys. Rev. Lett 102, 026002 (2009)) to describe some dynamical behaviors of soft-flowing materials is theoretically analyzed. Equilibrium and transport properties are derived within the framework of a continuum free-energy formulation, and checked against numerical simulations. Due to the competition between short-range inter-species repulsion and mid-range intra-species attraction, the model is shown to give rise to a very rich configurational dynamics of the density field, exhibiting numerous features of soft-flowing materials, such as long-time relaxation due to caging effects, enhanced viscosity and structural arrest, ageing under moderate shear and shear-thinning flow above a critical shear threshold. △ Less

Submitted 30 July, 2009; v1 submitted 19 June, 2009; originally announced June 2009.

Comments: 25 pages, 17 figures, submitted to Journal of chemical physics;

Journal ref: J. Chem. Phys. 131 104903 (2009)

Magnetic reversals in a simple model of MHD

Authors: Roberto Benzi, Jean-Francois Pinton

Abstract: We study a simple magnetohydrodynamical approach in which hydrodynamics and MHD turbulence are coupled in a shell model, with given dynamo constrains in the large scales. We consider the case of a low Prandtl number fluid for which the inertial range of the velocity field is much wider than that of the magnetic field. Random reversals of the magnetic field are observed and it shown that the magn… ▽ More We study a simple magnetohydrodynamical approach in which hydrodynamics and MHD turbulence are coupled in a shell model, with given dynamo constrains in the large scales. We consider the case of a low Prandtl number fluid for which the inertial range of the velocity field is much wider than that of the magnetic field. Random reversals of the magnetic field are observed and it shown that the magnetic field has a non trivial evolution linked to the nature of the hydrodynamics turbulence. △ Less

Submitted 2 June, 2009; originally announced June 2009.

Comments: 4 pages, submitted to PRL

Fisher equation with turbulence in one dimension

Authors: Roberto Benzi, David R. Nelson

Abstract: We investigate the dynamics of the Fisher equation for the spreading of micro-organisms in one dimension subject to both turbulent convection and diffusion. We show that for strong enough turbulence, bacteria, for example, track in a quasilocalized fashion (with remakably long persistance times) sinks in the turbulent field. An important consequence is a large reduction in the carrying capacity… ▽ More We investigate the dynamics of the Fisher equation for the spreading of micro-organisms in one dimension subject to both turbulent convection and diffusion. We show that for strong enough turbulence, bacteria, for example, track in a quasilocalized fashion (with remakably long persistance times) sinks in the turbulent field. An important consequence is a large reduction in the carrying capacity of the fluid medium. We determine analytically the regimes where this quasi-localized behavior occurs and test our predictions by numerical simulations. △ Less

Submitted 2 February, 2009; originally announced February 2009.

Comments: 15 pages

Mesoscopic lattice Boltzmann modeling of flowing soft systems

Authors: Roberto Benzi, Sergio Chibbaro, Sauro Succi

Abstract: A mesoscopic multi-component lattice Boltzmann model with short-range repulsion between different species and short/mid-ranged attractive/repulsive interactions between like-molecules is introduced. The interplay between these composite interactions gives rise to a rich configurational dynamics of the density field, exhibiting many features of disordered liquid dispersions (micro-emulsions) and… ▽ More A mesoscopic multi-component lattice Boltzmann model with short-range repulsion between different species and short/mid-ranged attractive/repulsive interactions between like-molecules is introduced. The interplay between these composite interactions gives rise to a rich configurational dynamics of the density field, exhibiting many features of disordered liquid dispersions (micro-emulsions) and soft-glassy materials, such as long-time relaxation due to caging effects, anomalous enhanced viscosity, ageing effects under moderate shear and flow above a critical shear rate. △ Less

Submitted 28 July, 2008; originally announced July 2008.

Comments: 4 pages, 4 figures

Journal ref: Physical Review Letters 102, 026002 (2009)

Stochastic Resonance: from climate to biology

Authors: Roberto Benzi

Abstract: In this paper I will review some basic aspects of the mechanism of stochastic resonance. Stochastic resonance was first introduced as a possible mechanism to explain long term climatic variation. Since then, there have been many applications of stochastic resonance in physical and biological systems. I will show that in complex system, stochastic resonance can substantially change as a function… ▽ More In this paper I will review some basic aspects of the mechanism of stochastic resonance. Stochastic resonance was first introduced as a possible mechanism to explain long term climatic variation. Since then, there have been many applications of stochastic resonance in physical and biological systems. I will show that in complex system, stochastic resonance can substantially change as a function of the ``system complexity''. Also, I will briefly mention how to apply stochastic resonance for the case of Brownian motors. △ Less

Submitted 5 February, 2007; originally announced February 2007.

Report number: Richardson Lecture, April 2006

Generalized Lattice Boltzmann Method with multi-range pseudo-potential

Authors: M. Sbragaglia, R. Benzi, L. Biferale, S. Succi, K. Sugiyama, F. Toschi

Abstract: The physical behaviour of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudo-potential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or h… ▽ More The physical behaviour of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudo-potential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or high order isotropy. Higher order schemes to implement self-consistent forcings are rigorously computed for 2d and 3d models. The extended scenario developed in this work clarifies the theoretical foundations of the Shan-Chen methodology for the lattice Boltzmann method and enhances its applicability and flexibility to the simulation of multiphase flows to density ratios up to O(100). △ Less

Submitted 7 November, 2006; originally announced November 2006.

Journal ref: Phys. Rev. E 75, 026702 (2007)

A mesoscopic approach for multi-phase flows in nano-corrugated channels

Authors: R. Benzi, L. Biferale, M. Sbragaglia, S. Succi, F. Toschi

Abstract: An approach based on a lattice version of the Boltzmann kinetic equation for describing multi-phase flows in nano- and micro-corrugated devices is proposed. We specialize it to describe the wetting/dewetting transition of fluids in presence of nanoscopic grooves etched on the boundaries. This approach permits to retain the essential supra-molecular details of fluid-solid interactions without sur… ▽ More An approach based on a lattice version of the Boltzmann kinetic equation for describing multi-phase flows in nano- and micro-corrugated devices is proposed. We specialize it to describe the wetting/dewetting transition of fluids in presence of nanoscopic grooves etched on the boundaries. This approach permits to retain the essential supra-molecular details of fluid-solid interactions without surrendering -actually boosting- the computational efficiency of continuum methods. The mesoscopic method is first validated quantitatively against Molecular Dynamics (MD) results of Cottin-Bizonne et al. [Nature Mater. 2, 237 (2003)] and then applied to more complex situations which are hardly accessible to MD simulations. The resulting analysis confirms that surface roughness and capillary effects may conspire to promote a counter-intuitive but significant reduction of the flow drag with substantial enhancement in the mass flow rates and slip-lengths in the micrometric range for highly hydrophobic surfaces. △ Less

Submitted 4 May, 2006; originally announced May 2006.

Journal ref: Phys. Rev. Lett. 97, 204503 (2006)

Mesoscopic modeling of a two-phase flow in the presence of boundaries: the Contact Angle

Authors: R. Benzi, L. Biferale, M. Sbragaglia, S. Succi, F. Toschi

Abstract: We present a mesoscopic model, based on the Boltzmann Equation, for the interaction between a solid wall and a non-ideal fluid. We present an analytic derivation of the contact angle in terms of the surface tension between the liquid-gas, the liquid-solid and the gas-solid phases. We study the dependency of the contact angle on the two free parameters of the model, which determine the interactio… ▽ More We present a mesoscopic model, based on the Boltzmann Equation, for the interaction between a solid wall and a non-ideal fluid. We present an analytic derivation of the contact angle in terms of the surface tension between the liquid-gas, the liquid-solid and the gas-solid phases. We study the dependency of the contact angle on the two free parameters of the model, which determine the interaction between the fluid and the boundaries, i.e. the equivalent of the wall density and of the wall-fluid potential in Molecular Dynamics studies. We compare the analytical results obtained in the hydrodynamical limit for the density profile and for the surface tension expression with the numerical simulations. We compare also our two-phase approach with some exact results for a pure hydrodynamical incompressible fluid based on Navier-Stokes equations with boundary conditions made up of alternating slip and no-slip strips. Finally, we show how to overcome some theoretical limitations connected with a discretized Boltzmann scheme and we discuss the equivalence between the surface tension defined in terms of the mechanical equilibrium and in terms of the Maxwell construction. △ Less

Submitted 3 February, 2006; originally announced February 2006.

Comments: 29 pages, 12 figures

Journal ref: Phys. Rev. E 74, 021509 (2006)

On the Anomalous Scaling Exponents in Nonlinear Models of Turbulence

Authors: Luiza Angheluta, Roberto Benzi, Luca Biferale, Itamar Procaccia, Federico Toschi

Abstract: We propose a new approach to the old-standing problem of the anomaly of the scaling exponents of nonlinear models of turbulence. We achieve this by constructing, for any given nonlinear model, a linear model of passive advection of an auxiliary field whose anomalous scaling exponents are the same as the scaling exponents of the nonlinear problem. The statistics of the auxiliary linear model are… ▽ More We propose a new approach to the old-standing problem of the anomaly of the scaling exponents of nonlinear models of turbulence. We achieve this by constructing, for any given nonlinear model, a linear model of passive advection of an auxiliary field whose anomalous scaling exponents are the same as the scaling exponents of the nonlinear problem. The statistics of the auxiliary linear model are dominated by `Statistically Preserved Structures' which are associated with exact conservation laws. The latter can be used for example to determine the value of the anomalous scaling exponent of the second order structure function. The approach is equally applicable to shell models and to the Navier-Stokes equations. △ Less

Submitted 4 May, 2006; v1 submitted 1 February, 2006; originally announced February 2006.

Comments: revised version with new data on Navier-Stokes eqs

Mesoscopic two-phase model for describing apparent slip in micro-channel flows

Authors: R. Benzi, L. Biferale, M. Sbragaglia, S. Succi, F. Toschi

Abstract: The phenomenon of apparent slip in micro-channel flows is analyzed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactins. The weakly-inhomogeneous limit of this model is solved analytically. The present mesoscopic approach permits to access much larger scales than molecular dynamics, and comparable with those attained by continuum… ▽ More The phenomenon of apparent slip in micro-channel flows is analyzed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactins. The weakly-inhomogeneous limit of this model is solved analytically. The present mesoscopic approach permits to access much larger scales than molecular dynamics, and comparable with those attained by continuum methods. However, at variance with the continuum approach, the existence of a gas layer near the wall does not need to be postulated a priori, but emerges naturally from the underlying non-ideal mesoscopic dynamics. It is therefore argued that a mesoscopic Lattice Boltzmann approach with non-ideal fluid-fluid and fluid-wall interactions might achieve an optimal compromise between physical realism and computational efficiency for the study of channel micro-flows. △ Less

Submitted 5 July, 2005; originally announced July 2005.

Comments: 5 pages, 3 figures

Journal ref: Europhys. Lett. 74, 651 (2006)

Mesoscopic modeling of heterogeneous boundary conditions for microchannel flows

Authors: R. Benzi, L. Biferale, M. Sbragaglia, S. Succi, F. Toschi

Abstract: We present a mesoscopic model of the fluid-wall interactions for flows in microchannel geometries. We define a suitable implementation of the boundary conditions for a discrete version of the Boltzmann equations describing a wall-bounded single phase fluid. We distinguish different slippage properties on the surface by introducing a slip function, defining the local degree of slip for mesoscopic… ▽ More We present a mesoscopic model of the fluid-wall interactions for flows in microchannel geometries. We define a suitable implementation of the boundary conditions for a discrete version of the Boltzmann equations describing a wall-bounded single phase fluid. We distinguish different slippage properties on the surface by introducing a slip function, defining the local degree of slip for mesoscopic molecules at the boundaries. The slip function plays the role of a renormalizing factor which incorporates, with some degree of arbitrariness, the microscopic effects on the mesoscopic description. We discuss the mesoscopic slip properties in terms of slip length, slip velocity, pressure drop reduction (drag reduction), and mass flow rate in microchannels as a function of the degree of slippage and of its spatial distribution and localization, the latter parameter mimicking the degree of roughness of the ultra-hydrophobic material in real experiments. We also discuss the increment of the slip length in the transition regime, i.e. at O(1) Knudsen numbers. Finally, we compare our results with Molecular Dynamics investigations of the dependency of the slip length on the mean channel pressure and local slip properties (Cottin-Bizonne et al. 2004) and with the experimental dependency of the pressure drop reduction on the percentage of hydrophobic material deposited on the surface -- Ou et al. (2004). △ Less

Submitted 25 February, 2005; originally announced February 2005.

Comments: 21 pages, 10 figures

Journal ref: J. Fluid Mech. 548, 257 (2006)

A simple model for drag reduction

Authors: Roberto Benzi, Itamar Procaccia

Abstract: Direct Numerical Simulations established that the FENE-P model of viscoelastic flows exhibits the phenomenon of turbulent drag reduction which is caused in experiments by dilute polymeric additives. To gain analytic understanding of the phenomenon we introduce in this Letter a simple 1-dimensional model of the FENE-P equations. We demonstrate drag reduction in the simple model, and explain analy… ▽ More Direct Numerical Simulations established that the FENE-P model of viscoelastic flows exhibits the phenomenon of turbulent drag reduction which is caused in experiments by dilute polymeric additives. To gain analytic understanding of the phenomenon we introduce in this Letter a simple 1-dimensional model of the FENE-P equations. We demonstrate drag reduction in the simple model, and explain analytically the main observations which include (i) reduction of velocity gradients for fixed throughput and (ii) increase of throughput for fixed dissipation. △ Less

Submitted 23 October, 2002; originally announced October 2002.

Comments: submitted to PRL