Source author record

M. Sbragaglia

M. 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

29works
8topics
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

29 published item(s)

preprint2020arXiv

On the impact of controlled wall roughness shape on the flow of a soft-material

We explore the impact of geometrical corrugations on the near-wall flow properties of a soft-material driven in a confined rough microchannel. By means of numerical simulations, we perform a quantitative analysis of the relation between the flow rate $Φ$ and the wall stress $σ_w$ for a number of setups, by changing both the roughness values as well as the roughness shape. Roughness suppresses the flow, with the existence of a characteristic value of $σ_w$ at which flow sets in. Just above the onset of flow, we quantitatively analyze the relation between $Φ$ and $σ_w$. While for smooth walls a linear dependency is observed, steeper behaviours are found to set in by increasing wall roughness. The variation of the steepness, in turn, depends on the shape of the wall roughness, wherein gentle steepness changes are promoted by a variable space localization of the roughness.

preprint2016arXiv

Deformation and break-up of viscoelastic droplets in confined shear flow

The deformation and break-up of Newtonian/viscoelastic droplets are studied in confined shear flow. Our numerical approach is based on a combination of Lattice-Boltzmann models (LBM) and finite difference schemes, the former used to model two immiscible fluids with variable viscous ratio, and the latter used to model the polymer dynamics. 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 quantify the droplet response by changing the polymer relaxation time $τ_P$, the maximum extensibility $L$ of the polymers, and the degree of confinement, i.e. the ratio of the droplet diameter to gap spacing. In unconfined shear flow, the effects of droplet viscoelasticity on the critical Capillary number $\mbox{Ca}_{\mbox{\tiny{cr}}}$ for break-up are moderate in all cases studied. However, in confined conditions a different behaviour is observed: the critical Capillary number of a viscoelastic droplet increases or decreases, depending on the maximum elongation of the polymers, the latter affecting the extensional viscosity of the polymeric solution. Force balance is monitored in the numerical simulations to validate the physical picture.

preprint2016arXiv

Fluctuating Multicomponent Lattice Boltzmann Model

Current implementations of fluctuating lattice Boltzmann equations (FLBE) describe single component fluids. In this paper, a model based on the continuum kinetic Boltzmann equation for describing multicomponent fluids is extended to incorporate the effects of thermal fluctuations. The thus obtained fluctuating Boltzmann equation is first linearized to apply the theory of linear fluctuations, and expressions for the noise covariances are determined by invoking the fluctuation-dissipation theorem (FDT) directly at the kinetic level. Crucial for our analysis is the projection of the Boltzmann equation onto the ortho-normal Hermite basis. By integrating in space and time the fluctuating Boltzmann equation with a discrete number of velocities, the FLBE is obtained for both ideal and non-ideal multicomponent fluids. Numerical simulations are specialized to the case where mean-field interactions are introduced on the lattice, indicating a proper thermalization of the system.

preprint2016arXiv

Hybrid Lattice Boltzmann/Finite Difference simulations of viscoelastic multicomponent flows in confined geometries

We propose numerical simulations of viscoelastic fluids based on a hybrid algorithm combining Lattice-Boltzmann models (LBM) and Finite Differences (FD) schemes, the former used to model the macroscopic hydrodynamic equations, and the latter used to model the polymer dynamics. 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). The numerical model is first benchmarked by characterizing the rheological behaviour of dilute homogeneous solutions in various configurations, including steady shear, elongational flows, transient shear and oscillatory flows. As an upgrade of complexity, we study the model in presence of non-ideal multicomponent interfaces, where immiscibility is introduced in the LBM description using the "Shan-Chen" model. The problem of a confined viscoelastic (Newtonian) droplet in a Newtonian (viscoelastic) matrix under simple shear is investigated and numerical results are compared with the predictions of various theoretical models. The proposed numerical simulations explore problems where the capabilities of LBM were never quantified before.

preprint2016arXiv

Lattice Boltzmann Simulations of Droplet formation in confined Channels with Thermocapillary flows

Based on mesoscale lattice Boltzmann simulations with the "Shan-Chen" model, we explore the influence of thermocapillarity on the break-up properties of fluid threads in a microfluidic T-junction, where a dispersed phase is injected perpendicularly into a main channel containing a continuous phase, and the latter induces periodic break-up of droplets due to the cross-flowing. Temperature effects are investigated by switching on/off both positive/negative temperature gradients along the main channel direction, thus promoting a different thread dynamics with anticipated/delayed break-up. Numerical simulations are performed at changing the flow-rates of both the continuous and dispersed phases, as well as the relative importance of viscous forces, surface tension forces and thermocapillary stresses. The range of parameters is broad enough to characterize the effects of thermocapillarity on different mechanisms of break-up in the confined T-junction, including the so-called "squeezing" and "dripping" regimes, previously identified in the literature. Some simple scaling arguments are proposed to rationalize the observed behaviour, and to provide quantitative guidelines on how to predict the droplet size after break-up.

preprint2016arXiv

Sliding drops across alternating hydrophobic and hydrophilic stripes

We perform a joint numerical and experimental study to sistematically characterize the motion of drops sliding over a periodic array of alternating hydrophobic and hydrophilic stripes with large wettability contrast, and typical width of hundreds of $μ\textrm{m}$. The fraction of the hydrophobic stripes has been varied from about 20% to 80%. The effects of the heterogeneous patterning can be described by a renormalized value of the critical Bond number, i.e. the critical dimensionless force needed to depin the drop before it starts to move. Close to the critical Bond number we observe a jerkily motion characterized by an evident stick-slip dynamics. As a result, dissipation is strongly localized in time, and the mean velocity of the drops can easily decrease by an order of magnitude compared to the sliding on homogeneous surface. Lattice Boltzmann (LB) numerical simulations are crucial for disclosing to what extent the sliding dynamics can be deduced from the computed balance of capillary, viscous and body forces at varying the Bond number, the surface composition and the liquid viscosity. Away from the critical Bond number, we characterize both experimentally and numerically the dissipation inside the droplet by studying the relation between the average velocity and the applied volume forces.

preprint2016arXiv

Thermal fluctuations of an interface near a contact line

The effect of thermal fluctuations near a contact line of a liquid interface partially wetting an impenetrable substrate is studied analytically and numerically. Promoting both the interface profile and the contact line position to random variables, we explore the equilibrium properties of the corresponding fluctuating contact line problem based on an interfacial Hamiltonian involving a "contact" binding potential. To facilitate an analytical treatment we consider the case of a one-dimensional interface. The effective boundary condition at the contact line is determined by a dimensionless parameter that encodes the relative importance of thermal energy and substrate energy at the microscopic scale. We find that this parameter controls the transition from a partially wetting to a pseudo-partial wetting state, the latter being characterized by a thin prewetting film of fixed thickness. In the partial wetting regime, instead, the profile typically approaches the substrate via an exponentially thinning prewetting film. We show that, independently of the physics at the microscopic scale, Young's angle is recovered sufficiently far from the substrate. The fluctuations of the interface and of the contact line give rise to an effective disjoining pressure, exponentially decreasing with height. Fluctuations therefore provide a regularization of the singular contact forces occurring in the corresponding deterministic problem.

preprint2016arXiv

Viscoelastic Multicomponent Fluids in confined Flow-Focusing Devices

The effects of elasticity on the break-up of liquid threads in microfluidic cross-junctions is investigated using numerical simulations based on the "lattice Boltzmann models" (LBM). Working at small Capillary numbers, we investigate the effects of non-Newtonian phases in the transition from droplet formation at the cross-junction (DCJ) and droplet formation downstream of the cross-junction (DC) (Liu & Zhang, ${\it Phys. Fluids.}$ ${\bf 23}$, 082101 (2011)). 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.

preprint2015arXiv

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

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.

preprint2015arXiv

Pair separation of magnetic elements in the quiet Sun

The dynamic properties of the quiet Sun photosphere can be investigated by analyzing the pair dispersion of small-scale magnetic fields (i.e., magnetic elements). By using $25$ hr-long Hinode magnetograms at high spatial resolution ($0".3$), we tracked $68,490$ magnetic element pairs within a supergranular cell near the disk center. The computed pair separation spectrum, calculated on the whole set of particle pairs independently of their initial separation, points out what is known as a super-diffusive regime with spectral index $γ=1.55\pm0.05$, in agreement with the most recent literature, but extended to unprecedented spatial and temporal scales (from granular to supergranular). Furthermore, for the first time, we investigated here the spectrum of the mean square displacement of pairs of magnetic elements, depending on their initial separation $r_0$. We found that there is a typical initial distance above (below) which the pair separation is faster (slower) than the average. A possible physical interpretation of such a typical spatial scale is also provided.

preprint2014arXiv

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

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.

preprint2014arXiv

Internal dynamics and activated processes in Soft-Glassy materials

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.

preprint2014arXiv

Natural convection with mixed insulating and conducting boundary conditions: low and high Rayleigh numbers regimes

We investigate the stability and dynamics of natural convection in two dimensions, subject to inhomogeneous boundary conditions. In particular, we consider a Rayleigh-Bènard (RB) cell, where the horizontal top boundary contains a periodic sequence of alternating thermal insulating and conducting patches, and we study the effects of the heterogeneous pattern on the global heat exchange, both at low and high Rayleigh numbers. At low Rayleigh numbers, we determine numerically the transition from a regime characterized by the presence of small convective cells localized at the inhomogeneous boundary to the onset of bulk convective rolls spanning the entire domain. Such a transition is also controlled analytically in the limit when the boundary pattern length is small compared with the cell vertical size. At higher Rayleigh number, we use numerical simulations based on a lattice Boltzmann method to assess the impact of boundary inhomogeneities on the fully turbulent regime up to $Ra \sim 10^{10}$.

preprint2014arXiv

Plastic flow of foams and emulsions in a channel

In order to understand the flow profiles of complex fluids, a crucial issue concerns the emergence of spatial correlations among plastic rearrangements exhibiting cooperativity flow behaviour at the macroscopic level. In this paper, the rate of plastic events in a Poiseuille flow is experimentally measured on a confined foam in a Hele-Shaw geometry. The correlation with independently measured velocity profiles is quantified. To go beyond a limitation of the experiments, namely the presence of wall friction which complicates the relation between shear stress and shear rate, we compare the experiments with simulations of emulsion droplets based on the lattice-Boltzmann method, which are performed both with, and without, wall friction. Our results indicate a correlation between the localisation length of the velocity profiles and the localisation length of the number of plastic events. Finally, unprecedented results on the distribution of the orientation of plastic events show that there is a non-trivial correlation with the underlying local shear strain. These features, not previously reported for a confined foam, lend further support to the idea that cooperativity mechanisms, originally invoked for concentrated emulsions (Goyon et al. 2008), have parallels in the behaviour of other soft-glassy materials.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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)].

preprint2013arXiv

Dynamic crossover length in soft-glassy materials

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.

preprint2013arXiv

Interaction Pressure Tensor for a class of Multicomponent Lattice Boltzmann models

We present a theory to obtain the pressure tensor for a class of non-ideal multicomponent lattice Boltzmann models, thus extending the theory presented by Shan (X. Shan, Phys. Rev. E 77, 066702 (2008)) for single component fluids. We obtain the correct form of the pressure tensor directly on the lattice and the resulting equilibrium properties are shown to agree very well with those measured from numerical simulations. Results are compared with those of alternative theories.

preprint2011arXiv

Convection in multiphase flows using Lattice Boltzmann methods

We present high resolution numerical simulations of convection in multiphase flows (boiling) using a novel algorithm based on a Lattice Boltzmann method. We first validate the thermodynamical and kinematical properties of the algorithm. Then, we perform a series of 3d numerical simulations at changing the mean properties in the phase diagram and compare convection with and without phase coexistence at $Ra \sim 10^7$. We show that in presence of nucleating bubbles non-Oberbeck Boussinesq effects develops, mean temperature profile becomes asymmetric, heat-transfer and heat-transfer fluctuations are enhanced. We also show that small-scale properties of velocity and temperature fields are strongly affected by the presence of buoyant bubble leading to high non-Gaussian profiles in the bulk.

preprint2011arXiv

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

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.

preprint2011arXiv

Second order closure for stratified convection: bulk region and overshooting

The parameterization of small-scale turbulent fluctuations in convective systems and in the presence of strong stratification is a key issue for many applied problems in oceanography, atmospheric science and planetology. In the presence of stratification, one needs to cope with bulk turbulent fluctuations and with inversion regions, where temperature, density -or both- develop highly non-linear mean profiles due to the interactions between the turbulent boundary layer and the unmixed -stable- flow above/below it. We present a second order closure able to cope simultaneously with both bulk and boundary layer regions, and we test it against high-resolution state-of-the-art 2D numerical simulations in a convective and stratified belt for values of the Rayleigh number, up to Ra = 10^9. Data are taken from a Rayleigh-Taylor system confined by the existence of an adiabatic gradient.

preprint2010arXiv

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

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.

preprint2010arXiv

High resolution numerical study of Rayleigh-Taylor turbulence using a thermal lattice Boltzmann scheme

We present results of a high resolution numerical study of two dimensional (2d) Rayleigh-Taylor turbulence using a recently proposed thermal lattice Boltzmann method (LBT). The goal of our study is both methodological and physical. We assess merits and limitations concerning small- and large-scale resolution/accuracy of the adopted integration scheme. We discuss quantitatively the requirements needed to keep the method stable and precise enough to simulate stratified and unstratified flows driven by thermal active fluctuations at high Rayleigh and high Reynolds numbers. We present data with spatial resolution up to 4096 x 10000 grid points and Rayleigh number up to Ra ~ 10^11 . The statistical quality of the data allows us to investigate velocity and temperature fluctuations, scale-by-scale, over roughly four decades. We present a detailed quantitative analysis of scaling laws in the viscous, inertial and integral range, supporting the existence of a Bolgiano-like inertial scaling, as expected in 2d systems. We also discuss the presence of small/large intermittent deviation to the scaling of velocity/temperature fluctuations and the Rayleigh dependency of gradients flatness.

preprint2010arXiv

Numerical simulations of compressible Rayleigh-Taylor turbulence in stratified fluids

We present results from numerical simulations of Rayleigh-Taylor turbulence, performed using a recently proposed lattice Boltzmann method able to describe consistently a thermal compressible flow subject to an external forcing. The method allowed us to study the system both in the nearly-Boussinesq and strongly compressible regimes. Moreover, we show that when the stratification is important, the presence of the adiabatic gradient causes the arrest of the mixing process.

preprint2010arXiv

Volumetric formulation of lattice Boltzmann models with energy conservation

We analyze a volumetric formulation of lattice Boltzmann for compressible thermal fluid flows. The velocity set is chosen with the desired accuracy, based on the Gauss-Hermite quadrature procedure, and tested against controlled problems in bounded and unbounded fluids. The method allows the simulation of thermohydrodyamical problems without the need to preserve the exact space-filling nature of the velocity set, but still ensuring the exact conservation laws for density, momentum and energy. Issues related to boundary condition problems and improvements based on grid refinement are also investigated.

preprint2009arXiv

Continuum Free-Energy formulation for a class of Lattice Boltzmann multiphase models

It is shown that the Shan-Chen (SC) model for non-ideal lattice fluids can be made compliant with a pseudo free-energy principle by simple addition of a gradient force, whose expression is uniquely specified in terms of the fluid density. This additional term is numerically shown to provide fairly negligible effects on the system evolution during phase-separation. To the best of our knowledge, these important properties of the SC model were not noted before. The approach developed in the present work is based on a continuum analysis: further extensions, more in line with a discrete lattice theory (X. Shan, {\it Phys Rev E}, {\bf 77} 066702 (2008)) can be envisaged for the future.

preprint2009arXiv

Lattice Boltzmann method with self-consistent thermo-hydrodynamic equilibria

Lattice kinetic equations incorporating the effects of external/internal force fields via a shift of the local fields in the local equilibria, are placed within the framework of continuum kinetic theory. The mathematical treatment reveals that, in order to be consistent with the correct thermo-hydrodynamical description, temperature must also be shifted, besides momentum. New perspectives for the formulation of thermo-hydrodynamic lattice kinetic models of non-ideal fluids are then envisaged. It is also shown that on the lattice, the definition of the macroscopic temperature requires the inclusion of new terms directly related to discrete effects. The theoretical treatment is tested against a controlled case with a non ideal equation of state.

preprint2009arXiv

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

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.

preprint2005arXiv

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

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.