Source author record

Mauro Sbragaglia

Mauro Sbragaglia appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2023arXiv

A reduced model for droplet dynamics in shear flows at finite capillary numbers

We propose an extension of the Maffettone-Minale (MM) model to predict the dynamics of an ellipsoidal droplet in a shear flow. The parameters of the MM model are traditionally retrieved in the framework of the perturbation theory for small deformations, i.e., small capillary numbers ($\mbox{Ca} \ll 1$) applied to Stokes equations. In this work, we take a novel route, in that we determine the model parameters at finite capillary numbers ($\mbox{Ca}\sim {\cal O}(1)$) without relying on perturbation theory results, while retaining a realistic representation in creeping time and steady deformation attained by the droplet for different realizations of the viscosity ratio $λ$ between the inner and the outer fluids. This extended MM (EMM) model hinges on an independent characterization of the process of droplet deformation via numerical simulations of Stokes equations employing the Immersed Boundary - Lattice Boltzmann (IB-LB) numerical techniques. Issues on droplet breakup are also addressed and discussed within the EMM model.

preprint2022arXiv

Lattice Boltzmann simulations of Plasma Wakefield Acceleration

We explore a novel simulation route for Plasma Wakefield Acceleration (PWFA) by using the computational method known as the Lattice Boltzmann Method (LBM). LBM is based on a discretization of the continuum kinetic theory while assuring the convergence towards hydrodynamics for coarse-grained fields (i.e., density, velocity, etc.). LBM is an established numerical analysis tool in computational fluid dynamics, able to efficiently bridge between kinetic theory and hydrodynamics, but its application in the context of PWFA has never been investigated so far. This paper aims at filling this gap. Results of LBM simulations for PWFA are discussed and compared with those of a code (Architect) implementing a Cold Fluid (CF) model for the plasma. In the hydrodynamic framework, we discuss the importance of regularization effects related to diffusion properties intrinsic of the LBM, allowing to go beyond the traditional CF approximations. Issues on computational efficiency are also addressed.

preprint2022arXiv

Mesoscale Modelling of the Tolman Length in Multi-component Systems

In this paper we analyze the curvature corrections to the surface tension in the context of the Shan-Chen (SC) multi-component Lattice Boltzmann method (LBM). We demonstrate that the same techniques recently applied in the context of the Shan-Chen multi-phase model can be applied to multi-component mixtures. We implement, as a new application, the calculation of the surface of tension radius $R_s$ through the minimization of the generalized surface tension $σ[R]$. In turn we are able to estimate the Tolman length, i.e. the first order coefficient of the curvature expansion of the surface tension $σ(R)$, as well as the higher order corrections, i.e. the curvature- and the Gaussian-rigidity coefficients. The SC multi-component model allows to model both fully-symmetric as well as asymmetric interactions among the components. By performing an extensive set of simulations we present a first example of tunable Tolman length in the mesoscopic model, being zero for symmetric interactions and different from zero otherwise. This result paves the way for controlling such interface properties which are paramount in presence of thermal fluctuations. All reported results can be independently reproduced through the "idea.deploy" framework available at https://github.com/lullimat/idea.deploy.

preprint2021arXiv

A Mesoscale Perspective on the Tolman Length

We demonstrate that the multi-phase Shan-Chen lattice Boltzmann method (LBM) yields a curvature dependent surface tension $σ$ as computed from three-dimensional hydrostatic droplets/bubbles simulations. Such curvature dependence is routinely characterized, at first order, by the so-called {\it Tolman length} $δ$. LBM allows to precisely compute $σ$ at the surface of tension $R_s$ and determine the Tolman length from the coefficient of the first order correction. The corresponding values of $δ$ display universality for different equations of state, following a power-law scaling near the critical temperature. The Tolman length has been studied so far mainly via computationally demanding molecular dynamics (MD) simulations or by means of density functional theory (DFT) approaches playing a pivotal role in extending Classical Nucleation Theory. The present results open a new hydrodynamic-compliant mesoscale arena, in which the fundamental role of the Tolman length, alongside real-world applications to cavitation phenomena, can be effectively tackled. All the results can be independently reproduced through the "idea.deploy" framework.

preprint2021arXiv

TLBfind: a Thermal Lattice Boltzmann code for concentrated emulsions with FINite-size Droplets

In this paper, we present TLBfind, a GPU code for simulating the hydrodynamics of droplets along with a dynamic temperature field. TLBfind hinges on a two-dimensional multi-component lattice Boltzmann (LB) model simulating a concentrated emulsion with finite-size droplets evolving in a thermal convective state, just above the transition from conduction to convection. The droplet concentration of the emulsion system is tenable, and at the core of the code lies the possibility to measure a large number of physical observables characterising the flow and droplets. Furthermore, TLBfind includes a parallel implementation on GPU of the Delaunay triangulation useful for the detection of droplets' plastic rearrangements, and several types of boundary conditions, supporting simulations of channels with structured rough walls.

preprint2020arXiv

Lattice Boltzmann simulations of droplet breakup in confined and time-dependent flows

We study droplet dynamics and breakup in generic time-dependent flows via a multicomponent Lattice Boltzmann algorithm, with emphasis on flow start up conditions. We first study droplet breakup in a confined oscillatory shear flow via two different protocols. In one set up, we start from an initially spherical droplet and turn on the flow abruptly ("shock method"); in the other protocol, we start from an initially spherical droplet as well, but we progressively increase the amplitude of the flow, by allowing the droplet to relax to the steady state for each increase in amplitude, before increasing the flow amplitude again ("relaxation method"). The two protocols are shown to produce substantially different breakup scenarios. The mismatch between these two protocols is also studied for variations in the flow topology, the degree of confinement and the inertia of the fluid. All results point to the fact that under extreme conditions of confinement the relaxation protocols can drive the droplets into metastable states, which break only for very intense flow amplitudes, but their stability is prone to external perturbations, such as an oscillatory driving force.

preprint2019arXiv

Unified theoretical and experimental view on transient shear banding

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.

preprint2016arXiv

A Lattice Boltzmann study of the effects of viscoelasticity on droplet formation in microfluidic cross-junctions

Based on mesoscale lattice Boltzmann (LB) numerical simulations, we investigate the effects of viscoelasticity on the break-up of liquid threads in microfluidic cross-junctions, where droplets are formed by focusing a liquid thread of a dispersed (d) phase into another co-flowing continuous (c) immiscible phase. Working at small Capillary numbers, we investigate the effects of non-Newtonian phases in the transition from droplet formation at the cross-junction (DCJ) to droplet formation downstream of the cross-junction (DC) (Liu $\&$ Zhang, ${\it Phys. ~Fluids.}$ ${\bf 23}$, 082101 (2011)). We will analyze cases with ${\it Droplet ~Viscoelasticity}$ (DV), where viscoelastic properties are confined in the dispersed phase, as well as cases with ${\it Matrix ~Viscoelasticity}$ (MV), where viscoelastic properties are confined in the continuous phase. Moderate flow-rate ratios $Q \approx {\cal O}(1)$ of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, where viscoelasticity is found to influence the break-up point of the threads, which moves closer to the cross-junction and stabilizes. This is attributed to an increase of the polymer feedback stress forming in the corner flows, where the side channels of the device meet the main channel. Quantitative predictions on the break-up point of the threads are provided as a function of the Deborah number, i.e. the dimensionless number measuring the importance of viscoelasticity with respect to Capillary forces.

preprint2016arXiv

Deformation and break-up of viscoelastic droplets Using Lattice Boltzmann Models

We investigate the break-up of Newtonian/viscoelastic droplets in a viscoelastic/Newtonian matrix under the hydrodynamic conditions of a confined shear flow. Our numerical approach is based on a combination of Lattice-Boltzmann models (LBM) and Finite Difference (FD) schemes. LBM are used to model two immiscible fluids with variable viscosity ratio (i.e. the ratio of the droplet to matrix viscosity); FD schemes are used to model viscoelasticity, and the kinetics of the polymers is introduced using constitutive equations for viscoelastic fluids with finitely extensible non-linear elastic dumbbells with Peterlin's closure (FENE-P). We study both strongly and weakly confined cases to highlight the role of matrix and droplet viscoelasticity in changing the droplet dynamics after the startup of a shear flow. Simulations provide easy access to quantities such as droplet deformation and orientation and will be used to quantitatively predict the critical Capillary number at which the droplet breaks, the latter being strongly correlated to the formation of multiple neckings at break-up. This study complements our previous investigation on the role of droplet viscoelasticity (A. Gupta \& M. Sbragaglia, {\it Phys. Rev. E} {\bf 90}, 023305 (2014)), and is here further extended to the case of matrix viscoelasticity.

preprint2016arXiv

Effects of viscoelasticity on droplet dynamics and break-up in microfluidic T-Junctions: a lattice Boltzmann study

The effects of viscoelasticity on the dynamics and break-up of fluid threads in microfluidic T-junctions are investigated using numerical simulations of dilute polymer solutions at changing the Capillary number ($\mbox {Ca}$), i.e. at changing the balance between the viscous forces and the surface tension at the interface, up to $\mbox{Ca} \approx 3 \times 10^{-2}$. A Navier-Stokes (NS) description of the solvent based on the lattice Boltzmann models (LBM) is here coupled to constitutive equations for finite extensible non-linear elastic dumbbells with the closure proposed by Peterlin (FENE-P model). We present the results of three-dimensional simulations in a range of $\mbox{Ca}$ which is broad enough to characterize all the three characteristic mechanisms of breakup in the confined T-junction, i.e. ${\it squeezing}$, ${\it dripping}$ and ${\it jetting}$ regimes. The various model parameters of the FENE-P constitutive equations, including the polymer relaxation time $τ_P$ and the finite extensibility parameter $L^2$, are changed to provide quantitative details on how the dynamics and break-up properties are affected by viscoelasticity. We will analyze cases with ${\it Droplet ~Viscoelasticity}$ (DV), where viscoelastic properties are confined in the dispersed (d) phase, as well as cases with ${\it Matrix ~Viscoelasticity}$ (MV), where viscoelastic properties are confined in the continuous (c) phase. Moderate flow-rate ratios $Q \approx {\cal O}(1)$ of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, as the flow driving the break-up process upstream of the emerging thread can be sensibly perturbed by the polymer stresses.

preprint2016arXiv

Fluidisation and plastic activity in a model soft-glassy material flowing in micro-channels with rough walls

By means of mesoscopic numerical simulations of a model soft-glassy material, we investigate the role of boundary roughness on the flow behaviour of the material, probing the bulk/wall and global/local rheologies. We show that the roughness reduces the wall slip induced by wettability properties and acts as a source of fluidisation for the material. A direct inspection of the plastic events suggests that their rate of occurrence grows with the fluidity field, reconciling our simulations with kinetic elasto-plastic descriptions of jammed materials. Notwithstanding, we observe qualitative and quantitative differences in the scaling, depending on the distance from the rough wall and on the imposed shear. The impact of roughness on the orientational statistics is also studied.

preprint2015arXiv

Sliding droplets of Xanthan solutions: a joint experimental and numerical study

We have investigated the sliding of droplets made of solutions of Xanthan, a stiff rodlike polysaccharide exhibiting a non-newtonian behavior, notably characterized by a shear-rate dependence of the viscosity. These experimental results are quantitatively compared with those of newtonian fluids (water). The impact of the non-newtonian behavior on the sliding process was shown through the relation between the average dimensionless velocity (i.e. the Capillary number) and the dimensionless volume forces (i.e. the Bond number). To this aim, it is needed to define operative strategies to compute the Capillary number for the shear thinning fluids and compare with the corresponding newtonian case. Results from experiments were complemented with lattice Boltzmann numerical simulations of sliding droplets, aimed to disentangle the influence that non-newtonian flow properties have on the sliding.

preprint2015arXiv

The importance of chemical potential in the determination of water slip in nanochannels

We investigate the slip properties of water confined in graphite-like nano-channels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle -- but a different microscopic surface potential -- do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.

preprint2014arXiv

Mesoscale Structures at Complex Fluid-Fluid Interfaces: a Novel Lattice Boltzmann / Molecular Dynamics Coupling

Complex fluid-fluid interfaces featuring mesoscale structures with adsorbed particles are key components of newly designed materials which are continuously enriching the field of soft matter. Simulation tools which are able to cope with the different scales characterizing these systems are fundamental requirements for efficient theoretical investigations. In this paper we present a novel simulation method, based on the approach of Ahlrichs and Dünweg [Ahlrichs and Dünweg, Int. J. Mod. Phys. C, 1998, 9, 1429], that couples the "Shan-Chen" multicomponent Lattice Boltzmann technique to off-lattice molecular dynamics to simulate efficiently complex fluid-fluid interfaces. We demonstrate how this approach can be used to study a wide class of challenging problems. Several examples are given, with an accent on bicontinuous phases formation in polyelectrolyte solutions and ferrofluid emulsions. We also show that the introduction of solvation free energies in the particle-fluid interaction unveils the hidden, multiscale nature of the particle-fluid coupling, allowing to treat symmetrically (and interchangeably) the on-lattice and off-lattice components of the system.

preprint2014arXiv

Mesoscopic simulation study of wall roughness effects in micro-channel flows of dense emulsions

We study the Poiseuille flow of a soft-glassy material above the jamming point, where the material flows like a complex fluid with Herschel- Bulkley rheology. Microscopic plastic rearrangements and the emergence of their spatial correlations induce cooperativity flow behavior whose effect is pronounced in presence of confinement. With the help of lattice Boltzmann numerical simulations of confined dense emulsions, we explore the role of geometrical roughness in providing activation of plastic events close to the boundaries. We probe also the spatial configuration of the fluidity field, a continuum quantity which can be related to the rate of plastic events, thereby allowing us to establish a link between the mesoscopic plastic dynamics of the jammed material and the macroscopic flow behaviour.

preprint2014arXiv

Non-locality and viscous drag effects on the shear localisation in soft-glassy materials

We study the Couette flow of a quasi-2d soft-glassy material in a Hele-Shaw geometry. The material is chosen to be above the jamming point, where a yield stress $σ_Y$ emerges, below which the material deforms elastically and above which it flows like a complex fluid according to a Herschel-Bulkley (HB) rheology. Simultaneously, the effect of the confining plates is modelled as an effective linear friction law, while the walls aside the Hele-Shaw cell are sufficiently close to each other to allow visible cooperativity effects in the velocity profiles (Goyon et al., Nature 454, 84-87 (2008)). The effects of cooperativity are parametrized with a steady-state diffusion-relaxation equation for the fluidity field $f = \dotγ/σ$, defined as the ratio between shear rate $\dotγ$ and shear stress $σ$. For particular rheological flow-curves (Bingham fluids), the problem is tackled analytically: we explore the two regimes $σ\gg σ_Y$ and $σ\approx σ_Y$ and quantify the effect of the extra localisation induced by the wall friction. Other rheo-thinning fluids are explored with the help of numerical simulations based on lattice Boltzmann models, revealing a robustness of the analytical findings. Synergies and comparisons with other existing works in the literature (Barry et al., Phil. Mag. Lett. 91, 432-440 (2011)) are also discussed.

preprint2014arXiv

Regularization of the slip length divergence in water nanoflows by inhomogeneities at the Angstrom scale

We performed non-equilibrium Molecular Dynamics simulations of water flow in nano-channels with the aim of discriminating {\it static} from {\it dynamic} contributions of the solid surface to the slip length of the molecular flow. We show that the regularization of the slip length divergence at high shear rates, formerly attributed to the wall dynamics, is controlled instead by its static properties. Surprisingly, we find that atomic displacements at the Angstrom scale are sufficient to remove the divergence of the slip length and realize the no-slip condition. Since surface thermal fluctuations at room temperature are enough to generate these displacements, we argue that the no-slip condition for water can be achieved also for ideal surfaces, which do not present any surface roughness.

preprint2013arXiv

Stick-Slip Sliding of Water Drops on Chemically Heterogeneous Surfaces

We present a comprehensive study of water drops sliding down chemically heterogeneous surfaces formed by a periodic pattern of alternating hydrophobic and hydrophilic stripes. Drops are found to undergo a stick-slip motion whose average speed is an order of magnitude smaller than that measured on a homogeneous surface having the same static contact angle. This motion is the result of the periodic deformations of the drop interface when crossing the stripes. Numerical simulations confirm this view and are used to elucidate the principles underlying the experimental observations.

preprint2012arXiv

Unified framework for a side-by-side comparison of different multicomponent algorithms: lattice Boltzmann vs. phase field model

Lattice Boltzmann Models (LBM) and Phase Field Models (PFM) are two of the most widespread approaches for the numerical study of multicomponent fluid systems. Both methods have been successfully employed by several authors but, despite their popularity, still remains unclear how to properly compare them and how they perform on the same problem. Here we present a unified framework for the direct (one-to-one) comparison of the multicomponent LBM against the PFM. We provide analytical guidelines on how to compare the Shan-Chen (SC) lattice Boltzmann model for non-ideal multicomponent fluids with a corresponding free energy (FE) lattice Boltzmann model. Then, in order to properly compare the LBM vs. the PFM, we propose a new formulation for the free energy of the Cahn-Hilliard/Navier-Stokes equations. Finally, the LBM model is numerically compared with the corresponding phase field model solved by means of a pseudo-spectral algorithm. This work constitute a first attempt to set the basis for a quantitative comparison between different algorithms for multicomponent fluids. We limit our scope to the few of the most common variants of the two most widespread methodologies, namely the lattice Boltzmann model (SC and FE variants) and the Phase Field Model.

preprint2011arXiv

Droplet size distribution in homogeneous isotropic turbulence

We study the physics of droplet breakup in a statistically stationary homogeneous and isotropic turbulent flow by means of high resolution numerical investigations based on the multicomponent lattice Boltzmann method. We verified the validity of the criterion proposed by Hinze (1955) for droplet breakup and we measured the full probability distribution function (pdf) of droplets radii at different Reynolds numbers and for different volume fraction. By means of a Lagrangian tracking we could follow individual droplets along their trajectories, define a local Weber number based on the velocity gradients and study its cross-correlation with droplet deformation.

preprint2010arXiv

Droplet breakup in homogeneous and isotropic turbulence

This fluid dynamics video shows the breakup of a droplet in a stationary homogeneous and isotropic turbulent flow. We consider droplets with the same density of the transporting fluid. The droplets and the fluid are numerically modelled by means of a multicompo- nent Lattice-Boltzmann method. The turbulent fluid is maintained through a large scale stirring force and the radius of stable droplets, for the parameters in our simulation, is larger than the Kolmogorov scale. Events of droplet deformation, break-up and aggregation are clearly visible from the movie. With the present database droplet evo- lution can be studied from both an Eulerian and Lagrangian point of view. The Kolmogorov-Hinze criteria for droplets break-up can be tested also by means of simulations with different viscosity contrast between the two components.

preprint2010arXiv

Lattice Boltzmann Methods for thermal flows: continuum limit and applications to compressible Rayleigh-Taylor systems

We compute the continuum thermo-hydrodynamical limit of a new formulation of lattice kinetic equations for thermal compressible flows, recently proposed in [Sbragaglia et al., J. Fluid Mech. 628 299 (2009)]. We show that the hydrodynamical manifold is given by the correct compressible Fourier- Navier-Stokes equations for a perfect fluid. We validate the numerical algorithm by means of exact results for transition to convection in Rayleigh-Bénard compressible systems and against direct comparison with finite-difference schemes. The method is stable and reliable up to temperature jumps between top and bottom walls of the order of 50% the averaged bulk temperature. We use this method to study Rayleigh-Taylor instability for compressible stratified flows and we determine the growth of the mixing layer at changing Atwood numbers up to At ~ 0.4. We highlight the role played by the adiabatic gradient in stopping the mixing layer growth in presence of high stratification and we quantify the asymmetric growth rate for spikes and bubbles for two dimensional Rayleigh- Taylor systems with resolution up to Lx \times Lz = 1664 \times 4400 and with Rayleigh numbers up to Ra ~ 2 \times 10^10.

preprint2010arXiv

Shear Banding from lattice kinetic models with competing interactions

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.