BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Succi

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

ResearcherUnclaimed source record
physics.flu-dyncond-mat.softphysics.comp-phcond-mat.stat-mechphysics.plasm-phcond-mat.mes-hallgr-qchep-phhep-thastro-ph.HEquant-phastro-ph.COastro-ph.IMastro-ph.SRBiological Physicsnlin.CGphysics.chem-phastro-ph.GAcond-mat.dis-nncond-mat.mtrl-scicond-mat.othercond-mat.str-elhep-latnlin.CDnucl-thphysics.atom-phphysics.med-ph

Catalog footprint

What is connected

53works
27topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy dissipation in flows through curved spaces

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved space: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature-induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Poiseuille flow in curved spaces

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cooling Effect of the Richtmyer-Meshkov Instability

We provide numerical evidence that the Richtmyer-Meshkov (RM) instability contributes to the cooling of a relativistic fluid. Due to the presence of jet particles traveling throughout the medium, shock waves are generated in the form of Mach cones. The interaction of multiple shock waves can trigger the RM instability, and we have found that this process leads to a down-cooling of the relativistic fluid. To confirm the cooling effect of the instability, shock tube Richtmyer-Meshkov instability simulations are performed. Additionally, in order to provide an experimental observable of the RM instability resulting from the Mach cone interaction, we measure the two particle correlation function and highlight the effects of the interaction. The simulations have been performed with an improved version of the relativistic lattice Boltzmann model, including general equations of state and external forces.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

DSMC-LBM mapping scheme for rarefied and non-rarefied gas flows

We present the formulation of a kinetic mapping scheme between the Direct Simulation Monte Carlo (DSMC) and the Lattice Boltzmann Method (LBM) which is at the basis of the hybrid model used to couple the two methods in view of efficiently and accurately simulate isothermal flows characterized by variable rarefaction effects. Owing to the kinetic nature of the LBM, the procedure we propose ensures to accurately couple DSMC and LBM at a larger Kn number than usually done in traditional hybrid DSMC-Navier-Stokes equation models. We show the main steps of the mapping algorithm and illustrate details of the implementation. Good agreement is found between the moments of the single particle distribution function as obtained from the mapping scheme and from independent LBM or DSMC simulations at the grid nodes where the coupling is imposed. We also show results on the application of the hybrid scheme based on a simpler mapping scheme for plane Poiseuille flow at finite Kn number. Potential gains in the computational efficiency assured by the application of the coupling scheme are estimated for the same flow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Flow through randomly curved manifolds

We have found that the relation between the flow through campylotic (generically curved) media, consisting of randomly located curvature perturbations, and the average Ricci scalar of the system exhibits two distinct functional expressions (hysteresis), depending on whether the typical spatial extent of the curvature perturbation lies above or below the critical value maximizing the overall Ricci curvature. Furthermore, the flow through such systems as a function of the number of curvature perturbations presents a sublinear behavior for large concentrations due to the interference between curvature perturbations that, consequently, produces a less curved space. For the purpose of this study, we have developed and validated a lattice kinetic model capable of describing fluid flow in arbitrarily curved manifolds, which allows to deal with highly complex spaces in a very compact and efficient way.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-Order Kinetic Relaxation Schemes as High-Accuracy Poisson Solvers

We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretise the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice Boltzmann Model for Electronic Structure Simulations

Recently, a new connection between density functional theory and kinetic theory has been proposed. In particular, it was shown that the Kohn-Sham (KS) equations can be reformulated as a macroscopic limit of the steady-state solution of a suitable single-particle kinetic equation. By using a discrete version of this new formalism, the exchange and correlation energies of simple atoms and the geometrical configuration of the methane molecule were calculated accurately. Here, we discuss the main ideas behind the lattice kinetic approach to electronic structure computations, offer some considerations for prospective extensions, and also show additional numerical results, namely the geometrical configuration of the water molecule.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice Boltzmann model for resistive relativistic magnetohydrodynamics

In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to $σ^{-\frac{1}{2}}$, $σ$ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Simulator for Transport Phenomena in Fluid Flows

Transport phenomena still stand as one of the most challenging problems in computational physics. By exploiting the analogies between Dirac and lattice Boltzmann equations, we develop a quantum simulator based on pseudospin-boson quantum systems, which is suitable for encoding fluid dynamics transport phenomena within a lattice kinetic formalism. It is shown that both the streaming and collision processes of lattice Boltzmann dynamics can be implemented with controlled quantum operations, using a heralded quantum protocol to encode non-unitary scattering processes. The proposed simulator is amenable to realization in controlled quantum platforms, such as ion-trap quantum computers or circuit quantum electrodynamics processors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Short-lived lattice quasiparticles for strongly interacting fluids

It is shown that lattice kinetic theory based on short-lived quasiparticles proves very effective in simulating the complex dynamics of strongly interacting fluids (SIF). In particular, it is pointed out that the shear viscosity of lattice fluids is the sum of two contributions, one due to the usual interactions between particles (collision viscosity) and the other due to the interaction with the discrete lattice (propagation viscosity). Since the latter is {\it negative}, the sum may turn out to be orders of magnitude smaller than each of the two contributions separately, thus providing a mechanism to access SIF regimes at ordinary values of the collisional viscosity. This concept, as applied to quantum superfluids in one-dimensional optical lattices, is shown to reproduce shear viscosities consistent with the AdS-CFT holographic bound on the viscosity/entropy ratio. This shows that lattice kinetic theory continues to hold for strongly coupled hydrodynamic regimes where continuum kinetic theory may no longer be applicable.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formal analogy between the Dirac equation in its Majorana form and the discrete-velocity version of the Boltzmann kinetic equation

We point out a formal analogy between the Dirac equation in Majorana form and the discrete-velocity version of the Boltzmann kinetic equation. By a systematic analysis based on the theory of operator splitting, this analogy is shown to turn into a concrete and efficient computational method, providing a unified treatment of relativistic and non-relativistic quantum mechanics. This might have potentially far-reaching implications for both classical and quantum computing, because it shows that, by splitting time along the three spatial directions, quantum information (Dirac-Majorana wavefunction) propagates in space-time as a classical statistical process (Boltzmann distribution).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic Formulation of the Kohn-Sham Equations for ab initio Electronic Structure Calculations

We introduce a new approach to density functional theory based on kinetic theory, showing that the Kohn-Sham equations can be derived as a macroscopic limit of a suitable Boltzmann kinetic equation in the limit of small mean free path versus the typical scale of density gradients (Chapman-Enskog expansion). To derive the approach, we first write the Schrödinger equation as a special case of a Boltzmann equation for a gas of quasi-particles, with the potential playing the role of an external source that generates and destroys particles, so as to drive the system towards the ground state. The ions are treated as classical particles, using the Born-Oppenheimer dynamics, or by imposing concurrent evolution with the electronic orbitals. In order to provide quantitative support to our approach, we implement a discrete (lattice) model and compute, the exchange and correlation energies of simple atoms, and the geometrical configuration of the methane molecule. Excellent agreement with values in the literature is found.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lattice Boltzmann model for ultra-relativistic flows

We develop a relativistic lattice Boltzmann model capable of describing relativistic fluid dynamics at ultra-high velocities, with Lorentz factors up to $γ\sim 10$. To this purpose, we first build a new lattice kinetic scheme by expanding the Maxwell-Jüttner distribution function in an orthogonal basis of polynomials and applying an appropriate quadrature, providing the discrete versions of the relativistic Boltzmann equation and the equilibrium distribution. To achieve ultra-high velocities, we include a flux limiter scheme, and introduce the bulk viscosity by a suitable extension of the discrete relativistic Boltzmann equation. The model is validated by performing simulations of shock waves in viscous quark-gluon plasmas and comparing with existing models, finding very good agreement. To the best of our knowledge, we for the first time successfully simulate viscous shock waves in the highly relativistic regime. Moreover, we show that our model can also be used for near-inviscid flows even at very high velocities. Finally, as an astrophysical application, we simulate a relativistic shock wave, generated by, say, a supernova explosion, colliding with a massive interstellar cloud, e.g. molecular gas.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic effects on the Richtmyer-Meshkov instability

Theoretical and numerical analysis of the relativistic effects on the Richtmyer-Meshkov (RM) instability reveals new and potentially very useful effects. We find that, in contrast with the non- relativistic case, the growth rate of the RM instability depends strongly on the equation of state of the fluid, opening up the possibility to infer equations of state from experimental observations of the RM instability. As opposed to the non-relativistic case, we also discover that, above a critical value of the fluid velocity, the growth rate of the instability counter-intuitively decreases due to the Lorentz's factor, and vanishes in the ultrarelativistic limit, as the speed of the particles approaches the speed of light. Both effects might prove very useful for leading-edge applications, such as the study of the equation of state of quark-gluon matter, and the design of fast ignition inertial confinement fusion (ICF) schemes. We perform a linear stability analysis to characterize the instability, for an arbitrary equation of state, and implement numerical simulations to study the instability in the non-linear regime, using the equation of state of an ideal gas. Furthermore, based on the numerical results, we propose a general expression that characterizes the long term evolution of the instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic Lattice Boltzmann Model with Improved Dissipation

We develop a relativistic lattice Boltzmann (LB) model, providing a more accurate description of dissipative phenomena in relativistic hydrodynamics than previously available with existing LB schemes. The procedure applies to the ultra-relativistic regime, in which the kinetic energy (temperature) far exceeds the rest mass energy, although the extension to massive particles and/or low temperatures is conceptually straightforward. In order to improve the description of dissipative effects, the Maxwell-Juettner distribution is expanded in a basis of orthonormal polynomials, so as to correctly recover the third order moment of the distribution function. In addition, a time dilatation is also applied, in order to preserve the compatibility of the scheme with a cartesian cubic lattice. To the purpose of comparing the present LB model with previous ones, the time transformation is also applied to a lattice model which recovers terms up to second order, namely up to energy-momentum tensor. The approach is validated through quantitative comparison between the second and third order schemes with BAMPS (the solution of the full relativistic Boltzmann equation), for moderately high viscosity and velocities, and also with previous LB models in the literature. Excellent agreement with BAMPS and more accurate results than previous relativistic lattice Boltzmann models are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ultrarelativistic Transport Coefficients in Two Dimensions

We compute the shear and bulk viscosities, as well as the thermal conductivity of an ultrarelativistic fluid obeying the relativistic Boltzmann equation in 2+1 space-time dimensions. The relativistic Boltzmann equation is taken in the single relaxation time approximation, based on two approaches, the first, due to Marle and using the Eckart decomposition, and the second, proposed by Anderson and Witting and using the Landau-Lifshitz decomposition. In both cases, the local equilibrium is given by a Maxwell-Juettner distribution. It is shown that, apart from slightly different numerical prefactors, the two models lead to a different dependence of the transport coefficients on the fluid temperature, quadratic and linear, for the case of Marle and Anderson-Witting, respectively. However, by modifying the Marle model according to the prescriptions given in Ref.[1], it is found that the temperature dependence becomes the same as for the Anderson-Witting model.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A fully relativistic lattice Boltzmann algorithm

Starting from the Maxwell-Juettner equilibrium distribution, we develop a relativistic lattice Boltzmann (LB) algorithm capable of handling ultrarelativistic systems with flat, but expanding, spacetimes. The algorithm is validated through simulations of quark-gluon plasma, yielding excellent agreement with hydrodynamic simulations. The present scheme opens the possibility of transferring the recognized computational advantages of lattice kinetic theory to the context of both weakly and ultra-relativistic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analogy between turbulence and quantum gravity: beyond Kolmogorov's 1941 theory

Simple arguments based on the general properties of quantum fluctuations have been recently shown to imply that quantum fluctuations of spacetime obey the same scaling laws of the velocity fluctuations in a homogeneous incompressible turbulent flow, as described by Kolmogorov 1941 (K41) scaling theory. Less noted, however, is the fact that this analogy rules out the possibility of a fractal quantum spacetime, in contradiction with growing evidence in quantum gravity research. In this Note, we show that the notion of a fractal quantum spacetime can be restored by extending the analogy between turbulence and quantum gravity beyond the realm of K41 theory. In particular, it is shown that compatibility of a fractal quantum-space time with the recent Horava-Lifshitz scenario for quantum gravity, implies singular quantum wavefunctions. Finally, we propose an operational procedure, based on Extended Self-Similarity techniques, to inspect the (multi)-scaling properties of quantum gravitational fluctuations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lattice Boltzmann models for non-ideal fluids with arrested phase-separation

The effects of mid-range repulsion in Lattice Boltzmann models on the coalescence/breakup behaviour of single-component, non-ideal fluids are investigated. It is found that mid-range repulsive interactions allow the formation of spray-like, multi-droplet configurations, with droplet size directly related to the strength of the repulsive interaction. The simulations show that just a tiny ten-percent of mid-range repulsive pseudo-energy can boost the surface/volume ratio of the phase- separated fluid by nearly two orders of magnitude. Drawing upon a formal analogy with magnetic Ising systems, a pseudo-potential energy is defined, which is found to behave like a quasi-conserved quantity for most of the time-evolution. This offers a useful quantitative indicator of the stability of the various configurations, thus helping the task of their interpretation and classification. The present approach appears to be a promising tool for the computational modelling of complex flow phenomena, such as atomization, spray formation and micro-emulsions, break-up phenomena and possibly glassy-like systems as well.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nanoflows through disordered media: a joint Lattice Boltzmann and Molecular Dynamics investigation

We investigate nanoflows through dilute disordered media by means of joint lattice Boltzmann (LB) and molecular dynamics (MD) simulations -- when the size of the obstacles is comparable to the size of the flowing particles -- for randomly located spheres and for a correlated particle-gel. In both cases at sufficiently low solid fraction, $Φ<0.01$, LB and MD provide similar values of the permeability. However, for $Φ> 0.01$, MD shows that molecular size effects lead to a decrease of the permeability, as compared to the Navier-Stokes predictions. For gels, the simulations highlights a surplus of permeability, which can be accommodated within a rescaling of the effective radius of the gel monomers.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hydrokinetic simulations of nanoscopic precursor films in rough channels

We report on simulations of capillary filling of high-wetting fluids in nano-channels with and without obstacles. We use atomistic (molecular dynamics) and hydrokinetic (lattice-Boltzmann) approaches which point out clear evidence of the formation of thin precursor films, moving ahead of the main capillary front. The dynamics of the precursor films is found to obey a square-root law as the main capillary front, z^2(t) ~ t, although with a larger prefactor, which we find to take the same value for the different geometries (2D-3D) under inspection. The two methods show a quantitative agreement which indicates that the formation and propagation of thin precursors can be handled at a mesoscopic/hydrokinetic level. This can be considered as a validation of the Lattice-Boltzmann (LB) method and opens the possibility of using hydrokinetic methods to explore space-time scales and complex geometries of direct experimental relevance. Then, LB approach is used to study the fluid behaviour in a nano-channel when the precursor film encounters a square obstacle. A complete parametric analysis is performed which suggests that thin-film precursors may have an important influence on the efficiency of nanochannel-coating strategies.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lattice Boltzmann scheme for relativistic fluids

A Lattice Boltzmann formulation for relativistic fluids is presented and numerically verified through quantitative comparison with recent hydrodynamic simulations of relativistic shock-wave propagation in viscous quark-gluon plasmas. This formulation opens up the possibility of exporting the main advantages of Lattice Boltzmann methods to the relativistic context, which seems particularly useful for the simulation of relativistic fluids in complicated geometries.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Lattice formulation of the Fokker-Planck equation

A lattice version of the Fokker-Planck equation (FPE), accounting for dissipative interactions, not resolved on the molecular scale, is introduced. The lattice FPE is applied to the study of electrorheological transport of a one-dimensional charged fluid, and found to yield quantitative agreement with a recent analytical solution. Future extensions, including inelastic ion-ion collisions, are also outlined.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access