Researcher profile

S. Succi

S. Succi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 published item(s)

preprint2022arXiv

Curvature dynamics and long-range effects on fluid-fluid interfaces with colloids

We investigate the dynamics of a phase-separating binary fluid, containing colloidal dumbbells anchored to the fluid-fluid interface. Extensive Lattice Boltzmann-Immersed Boundary method simulations reveal that the presence of soft dumbbells can significantly affect the curvature dynamics of the interface between phase-separating fluids, even though the coarsening dynamics is left nearly unchanged. In addition, our results show that the curvature dynamics exhibits distinct non-local effects, which might be exploited for the design of new soft mesoscale materials. We point out that the inspection of the statistical dynamics of the curvature can disclose new insights into local inhomogeneities of the binary fluid configuration, as a function of the volume fraction and aspect ratio of the dumbbells.

preprint2022arXiv

Dynamics of polydisperse multiple emulsions in microfluidic channels

Multiple emulsions are a class of soft fluid in which small drops are immersed within a larger one and stabilized over long periods of time by a surfactant. We recently showed that, if a monodisperse multiple emulsion is subject to a pressure-driven flow, a wide variety of nonequilibrium steady states emerges at late times, whose dynamics relies on a complex interplay between hydrodynamic interactions and multibody collisions among internal drops. In this work, we use lattice Boltzmann simulations to study the dynamics of polydisperse double emulsions driven by a Poiseuille flow within a microfluidic channel. Our results show that their behavior is critically affected by multiple factors, such as initial position, polydispersity index, and area fraction occupied within the emulsion. While at low area fraction inner drops may exhibit either a periodic rotational motion (at low polydispersity) or arrange into nonmotile configurations (at high polydispersity) located far from each other, at larger values of area fraction they remain in tight contact and move unidirectionally. This decisively conditions their close-range dynamics, quantitatively assessed through a time-efficiency-like factor. Simulations also unveil the key role played by the capsule, whose shape changes can favor the formation of a selected number of nonequilibrium states in which both motile and nonmotile configurations are found.

preprint2021arXiv

The vortex-driven dynamics of droplets within droplets

Understanding the fluid-structure interaction is crucial for an optimal design and manufacturing of soft mesoscale materials. Multi-core emulsions are a class of soft fluids assembled from cluster configurations of deformable oil-water double droplets (cores), often employed as building-blocks for the realisation of devices of interest in bio-technology, such as drug-delivery, tissue engineering and regenerative medicine. Here, we study the physics of multi-core emulsions flowing in microfluidic channels and report numerical evidence of a surprisingly rich variety of driven non-equilibrium states (NES), whose formation is caused by a dipolar fluid vortex triggered by the sheared structure of the flow carrier within the microchannel. The observed dynamic regimes range from long-lived NES at low core-area fraction, characterised by a planetary-like motion of the internal drops, to short-lived ones at high core-area fraction, in which a pre-chaotic motion results from multi-body collisions of inner drops, as combined with self-consistent hydrodynamic interactions. The onset of pre-chaotic behavior is marked by transitions of the cores from one vortex to another, a process that we interpret as manifestations of the system to maximize its entropy by filling voids, as they arise dynamically within the capsule.

preprint2020arXiv

Dissipative hydrodynamics of relativistic shock waves in a Quark Gluon Plasma: comparing and benchmarking alternate numerical methods

This paper presents numerical cross-comparisons and benchmark results for two different kinetic numerical methods, capable of describing relativistic dissipative fluid dynamics in a wide range of kinematic regimes, typical of relevant physics applications, such as transport phenomena in quark-gluon plasmas. We refer to relativistic lattice Boltzmann versus Montecarlo Test-Particle methods. Lacking any realistic option for accurate validation vis-a-vis experimental data, we check the consistency of our results against established simulation packages available in the literature. We successfully cross-compare the results of the two aforementioned numerical approaches for momentum integrated quantities like the hydrostatic and dynamical pressure profiles, the collective flow and the heat flux. These results corroborate the confidence on the robustness and correctness of these computational methods and on the accurate calibration of their numerical parameters with respect to the physical transport coefficients. Our numerical results are made available as supplemental material, with the aim of establishing a reference benchmark for other numerical approaches.

preprint2020arXiv

Micro-vorticity fluctuations affect the structure of thin fluid films

The dynamic interaction of complex fluid interfaces is highly sensitive to near-contact interactions occurring at the scale of ten of nanometers. Such interactions are difficult to analyse because they couple self-consistently to the dynamic morphology of the evolving interface, as well as to the hydrodynamics of the interstitial fluid film. In this work, we show that, above a given magnitude threshold, near-contact interactions trigger non-trivial micro-vorticity patterns, which in turn affect the effective near-contact interactions, giving rise to persistent fluctuating ripples at the fluid interface. In such regime, near-contact interactions may significantly affect the macroscopic arrangement of emulsion configurations, such as those arising in soft-flowing microfluidic crystals.

preprint2020arXiv

Modelling drug delivery from multiple emulsions

We present a mechanistic model of drug release from a multiple emulsion into an external surrounding fluid. We consider a single multi-layer droplet where the drug kinetics are described by a pure diffusive process through different liquid shells. The multi-layer problem is described by a system of diffusion equations coupled via interlayer conditions imposing continuity of drug concentration and flux. Mass resistance is imposed at the outer boundary through the application of a surfactant at the external surface of the droplet. The two-dimensional problem is solved numerically by finite volume discretization. Concentration profiles and drug release curves are presented for three typical round-shaped (circle, ellipse and bullet) droplets and the dependency of the solution on the mass transfer coefficient at the surface analyzed. The main result shows a reduced release time for an increased elongation of the droplets.

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.

preprint2020arXiv

Probing bulk viscosity in relativistic flows

We derive an analytical connection between kinetic relaxation rate and bulk viscosity of a relativistic fluid in d spatial dimensions, all the way from the ultra-relativistic down to the near non-relativistic regime. Our derivation is based on both Chapman-Enskog asymptotic expansion and Grad's method of moments. We validate our theoretical results against a benchmark flow, providing further evidence of the correctness of the Chapman-Enskog approach; we define the range of validity of this approach and provide evidence of mounting departures at increasing Knudsen number. Finally, we present numerical simulations of transport processes in quark gluon plasmas, with special focus on the effects of bulk viscosity which might prove amenable to future experimental verification.

preprint2020arXiv

Relativistic Lattice Boltzmann Methods: Theory and Applications

We present a systematic account of recent developments of the relativistic Lattice Boltzmann method (RLBM) for dissipative hydrodynamics. We describe in full detail a unified, compact and dimension-independent procedure to design relativistic LB schemes capable of bridging the gap between the ultra-relativistic regime, $k_{\rm B} T \gg mc^2$, and the non-relativistic one, $k_{\rm B} T \ll mc^2$. We further develop a systematic derivation of the transport coefficients as a function of the kinetic relaxation time in $d=1,2,3$ spatial dimensions. The latter step allows to establish a quantitative bridge between the parameters of the kinetic model and the macroscopic transport coefficients. This leads to accurate calibrations of simulation parameters and is also relevant at the theoretical level, as it provides neat numerical evidence of the correctness of the Chapman-Enskog procedure. We present an extended set of validation tests, in which simulation results based on the RLBMs are compared with existing analytic or semi-analytic results in the mildly-relativistic ($k_{\rm B} T \sim mc^2$) regime for the case of shock propagations in quark-gluon plasmas and laminar electronic flows in ultra-clean graphene samples. It is hoped and expected that the material collected in this paper may allow the interested readers to reproduce the present results and generate new applications of the RLBM scheme.

preprint2019arXiv

Mathematical modelling of drug delivery from pH-responsive nanocontainers

Drug delivery systems represent a promising strategy to treat cancer and to overcome the side effects of chemotherapy. In particular, polymeric nanocontainers have attracted major interest because of their structural and morphological advantages and the variety of polymers that can be used, allowing the synthesis of materials capable of responding to the biochemical alterations of the tumour microenvironment. While experimental methodologies can provide much insight, the generation of experimental data across a wide parameter space is usually prohibitively time consuming and/or expensive. To better understand the influence of varying design parameters on the drug release profile and drug kinetics involved, appropriately-designed mathematical models are of great benefit. Here, we developed a novel mathematical model to describe drug transport within, and release from, a hollow nanocontainer consisting of a core and a pH-responsive polymeric shell. The two-layer mathematical model fully accounts for drug dissolution, diffusion and interaction with polymer. We generated experimental drug release profiles using daunorubicin and [Cu(TPMA)(Phenantroline)](ClO_4)_2 as model drugs, for which the nanocontainers exhibited excellent encapsulation ability. The in vitro drug release behaviour was studied under different conditions, where the system proved capable of responding to the selected pH stimuli by releasing a larger amount of drug in an acidic than in the physiological environments. By comparing the results of the mathematical model with our experimental data, we were able to identify the model parameter values that best-fit the data and demonstrate that the model is capable of describing the phenomena at hand. The proposed methodology can be used to describe and predict the release profiles for a variety of drug delivery systems.

preprint2019arXiv

Novel non-equilibrium steady states in multiple emulsions

We numerically investigate the rheological response of a non-coalescing multiple emulsion under a symmetric shear flow. We find that the dynamics significantly depends on the magnitude of the shear rate and on the number of the encapsulated droplets, two key parameters whose control is fundamental to accurately select the resulting non-equilibrium steady states. The double emulsion, for instance, attains a static steady state in which the external droplet stretches under flow and achieves an elliptical shape (closely resembling the one observed in a sheared isolated fluid droplet), while the internal one remains essentially unaffected. Novel non-equilibrium steady states arise in a multiple emulsion. Under a low/moderate shear rates, for instance, the encapsulated droplets display a non-trivial planetary-like motion that considerably affects the shape of the external droplet. Some features of this dynamic behavior are partially captured by the Taylor deformation parameter and the stress tensor. Besides a theoretical interest on its own, our results can potentially stimulate further experiments, as most of the predictions could be tested in the lab by monitoring droplets shapes and position over time.

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

Hydrodynamic Model for Conductivity in Graphene

Based on the recently developed picture of an electronic ideal relativistic fluid at the Dirac point, we present an analytical model for the conductivity in graphene that is able to describe the linear dependence on the carrier density and the existence of a minimum conductivity. The model treats impurities as submerged rigid obstacles, forming a disordered medium through which graphene electrons flow, in close analogy with classical fluid dynamics. To describe the minimum conductivity, we take into account the additional carrier density induced by the impurities in the sample. The model, which predicts the conductivity as a function of the impurity fraction of the sample, is supported by extensive simulations for different values of ${\cal E}$, the dimensionless strength of the electric field, and provides excellent agreement with experimental data.

preprint2012arXiv

Klein Tunneling in the presence of random impurities

In this paper, we study Klein tunneling in random media. To this purpose, we simulate the propagation of a relativistic Gaussian wavepacket through a graphene sample with randomly distributed potential barriers (impurities). The simulations, based on a relativistic quantum lattice Boltzmann method, permit to compute the transmission coefficient across the sample, thereby providing an estimate for the conductivity as a function of impurity concentration and strength of the potentials. It is found that the conductivity loss due to impurities is significantly higher for wave-packets of massive particles, as compared to massless ones. A general expression for the loss of conductivity as a function of the impurity percentage is presented and successfully compared with the Kozeny-Carman law for disordered media in classical fluid dynamics.

preprint2012arXiv

Pre-Turbulent Regimes in Graphene Flows

We provide numerical evidence that electronic pre-turbulent phenomena in graphene could be observed, under current experimental conditions, through detectable current fluctuations, echoing the detachment of vortices past localized micron-sized impurities. Vortex generation, due to micron sized constriction, is also explored with special focus on the effects of relativistic corrections to the normal Navier-Stokes equations. These corrections are found to cause a delay in the stability breakout of the fluid as well as a small shift in the vortex shedding frequency. Finally, a relation between the Strouhal number, a dimensionless measure of the vortex shedding frequency, and the Reynolds number is provided under conditions of interest for future experiments.

preprint2012arXiv

Single to Double Hump Transition in the Equilibrium Distribution Function of Relativistic Particles

We unveil a transition from single peaked to bimodal velocity distribution in a relativistic fluid under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broadening of the bell-shaped distribution is observed. Such transition results from the interplay between the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell-Jüttner distributions, all exhibiting the same qualitative behavior. We characterize the nature of the transition in the framework of critical phenomena and show that it is either continuous or discontinuous, depending on the group velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on two-dimensions, for which a possible experiment in graphene, based on the measurement of the Johnson-Nyquist noise, is proposed.

preprint2012arXiv

Taylor-Couette Instability in General Manifolds: A Lattice Kinetic Approach

We present a new lattice kinetic method to simulate fluid dynamics in curvilinear geometries. A suitable discrete Boltzmann equation is solved in contravariant coordinates, and the equilibrium distribution function is obtained by a Hermite polynomials expansion of the Maxwell-Boltzmann distribution, expressed in terms of the contravariant coordinates and the metric tensor. To validate the model, we calculate the critical Reynolds number for the onset of the Taylor-Couette instability between two concentric cylinders, obtaining excellent agreement with the theory. In order to extend this study to more general geometries, we also calculate the critical Reynolds number for the case of two concentric spheres, finding good agreement with experimental data. In the case of two concentric tori, we have found that the critical Reynolds is about 10% larger than the respective value for the two concentric cylinders.

preprint2011arXiv

Comment on the paper Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng and Wei Zhang, Numerics of the lattice Boltzmann method: Effects of collision models on the lattice Boltzmann simulations, Physical Review E 83, 056710 (2011)

Critical comments on the entropic lattice Boltzmann equation (ELBE), by Li-Shi Luo, Wei Liao, Xingwang Chen, Yan Peng and Wei Zhang, Physical Review E 83, 056710 (2011), are based on simulations which make use of a model that, despite being called ELBE by the authors, is in fact fully equivalent to the standard lattice Bhatnagar-Gross-Krook equation. As a result, the conclusion of Luo et al on ELBE is circular, hence devoid of scientific bearing.

preprint2011arXiv

On the relativistic Lattice Boltzmann method for quark-gluon plasma simulations

In this paper, we investigate the recently developed lattice Boltzmann model for relativistic hydrodynamics. To this purpose, we perform simulations of shock waves in quark-gluon plasma in the low and high viscosities regime, using three different computational models, the relativistic lattice Boltzmann (RLB), the Boltzmann Approach Multi-Parton Scattering (BAMPS), and the viscous sharp and smooth transport algorithm (vSHASTA). From the results, we conclude that the RLB model departs from BAMPS in the case of high speeds and high temperature(viscosities), the departure being due to the fact that the RLB is based on a quadratic approximation of the Maxwell-Jüttner distribution, which is only valid for sufficiently low temperature and velocity. Furthermore, we have investigated the influence of the lattice speed on the results, and shown that inclusion of quadratic terms in the equilibrium distribution improves the stability of the method within its domain of applicability. Finally, we assess the viability of the RLB model in the various parameter regimes relevant to ultra-relativistic fluid dynamics.

preprint2010arXiv

Derivation of the Lattice Boltzmann Model for Relativistic Hydrodynamics

A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic problems, namely shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast-wave on massive interstellar clouds. Close to second order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time-steps, are reported.