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

Sauro Succi

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

ResearcherUnclaimed source record
physics.comp-phphysics.flu-dyncond-mat.softcond-mat.mes-hallcond-mat.stat-mechgr-qcquant-phnlin.CGphysics.soc-phastro-ph.IMBiological Physicscomp-gasComputer VisionDigital Librariesphysics.data-anastro-ph.COComputationComputational Engineering, Finance, and Sciencecond-mat.dis-nncond-mat.mtrl-scicond-mat.quant-gascs.CYmath.APphysics.chem-phphysics.class-phphysics.ed-phphysics.hist-phphysics.pop-phPopulations and Evolution

Catalog footprint

What is connected

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

74 published item(s)

preprint2026arXiv

Ghost-Mode Filtered Fluctuating Lattice Boltzmann Method

Fluctuating lattice Boltzmann solvers are widely employed to model mesoscopic fluid behavior in soft-matter systems, including colloidal suspensions and dilute polymer solutions. Despite their utility, these methods can lose accuracy and stability when non-hydrodynamic modes interfere with the dynamics, especially in single--relaxation-time schemes. Here, we introduce a ghost-mode filtered fluctuating lattice Boltzmann method (GMF-FLBM) for the D3Q27 lattice, obtained by selectively eliminating the propagation of the ghost deterministic content while preserving the necessary stochastic forcing. We show, over a broad range of relaxation times, that GMF-FLBM recovers the amplitudes of equilibrium fluctuations with a comparable accuracy as a fully regularized high-order formulation, while requiring only minor adjustments to the conventional BGK collision framework.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Quantum Algorithm for the Fixed-Radius Neighbor Search

Neighbor search is a computationally demanding problem, usually both time- and memory-consuming. The main problem of this kind of algorithms is the long execution time due to cache misses. In this work, we propose a quantum algorithm for the Fixed RAdius Neighbor Search problem (FRANS) based on the fixed-point version of Grover's algorithm. We propose an efficient circuit for solving the FRANS with linear query complexity with the number of particles $N$. The quantum circuit returns the list of all the neighbors' pairs within the fixed radius, together with their distance, avoiding the slow down given by cache miss. We analyzed the gate and the query complexity of the circuit. Our FRANS algorithm presents a query complexity of $\mathcal{O}(N/\sqrt{M})$, where $M$ is the number of solutions, reaching the optimal lower bound of the Grover's algorithm. We propose different implementations of the oracle, which must be chosen depending on the precise structure of the database. Among these, we present an implementation using the Chebyshev distance with depth $\mathcal{O}(q_1)$, where $2^{q_1}$ is the number of grid points used to discretize a spatial dimension. State-of-the-art algorithms for state preparation allow for a trade-off between depth and width of the circuit, with a volume (depth$\times$ width) of $\mathcal{O}(N\log(N))$. This unfavorable scaling can be brought down to $\mathcal{O}(\text{poly}(\log N))$ in case of structured datasets. We proposed a stopping criterion based on Bayes interference and tested its validity on $1D$ simulations. Finally, we accounted for the readout complexity and assessed the resilience of the model to the readout error, suggesting an error correction-free strategy to check the accuracy of the results.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

A new kind of science

We discuss whether science is in the process of being transformed from a quest for causality to a quest for correlation in light of the recent development in artificial intelligence. We observe that while a blind trust in the most seductive promises of AI is surely to be avoided, a judicious combination of computer simulation based on physical insight and the machine learning ability to explore ultra-dimensional spaces, holds potential for transformative progress in the way science is going to be pursued in the years to come.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Benchmarking YOLOv5 and YOLOv7 models with DeepSORT for droplet tracking applications

Tracking droplets in microfluidics is a challenging task. The difficulty arises in choosing a tool to analyze general microfluidic videos to infer physical quantities. The state-of-the-art object detector algorithm You Only Look Once (YOLO) and the object tracking algorithm Simple Online and Realtime Tracking with a Deep Association Metric (DeepSORT) are customizable for droplet identification and tracking. The customization includes training YOLO and DeepSORT networks to identify and track the objects of interest. We trained several YOLOv5 and YOLOv7 models and the DeepSORT network for droplet identification and tracking from microfluidic experimental videos. We compare the performance of the droplet tracking applications with YOLOv5 and YOLOv7 in terms of training time and time to analyze a given video across various hardware configurations. Despite the latest YOLOv7 being 10% faster, the real-time tracking is only achieved by lighter YOLO models on RTX 3070 Ti GPU machine due to additional significant droplet tracking costs arising from the DeepSORT algorithm. This work is a benchmark study for the YOLOv5 and YOLOv7 networks with DeepSORT in terms of the training time and inference time for a custom dataset of microfluidic droplets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Collective Behavior of Crowded Drops in Microfluidic Systems

Droplet microfluidics, in which micro-droplets serve as individual reactors, has enabled a wide range of high-throughput biochemical processes. Unlike solid wells typically used in current biochemical assays, droplets are subject to instability and can undergo breakup, especially under fast flow conditions. Although the physics of single drops has been studied extensively, the flow of crowded drops or concentrated emulsions, where droplet volume fraction exceeds 80 percent, is relatively unexplored in microfluidics. In this article and the related invited lecture from the 74th Annual Meeting of the American Physical Society's Division of Fluid Dynamics, we describe the collective behavior of drops in a concentrated emulsion by tracking the dynamics and the fate of individual drops within the emulsion. At the slow flow limit of the concentrated emulsion, we observe an unexpected order, where the velocity of individual drops in the emulsion exhibits spatiotemporal periodicity. As the flow rate increases, the emulsion transitions from a solid-like to a liquid-like material, and the spatiotemporal order in the flow is lost. At the fast flow limit, droplet breakup starts to occur. We show that droplet breakup within the emulsion follows a probability distribution, in stark contrast to the deterministic behavior in classical single-drop studies. In addition to capillary number and viscosity ratio, break-up probability is governed by a confinement factor that measures drop size relative to a characteristic channel length. The breakup probability arises from the time-varying packing configuration of the drops. Finally, we discuss recent progress in computation methods for recapitulating the flow of concentrated emulsions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

DropTrack -- automatic droplet tracking using deep learning for microfluidic applications

Deep neural networks are rapidly emerging as data analysis tools, often outperforming the conventional techniques used in complex microfluidic systems. One fundamental analysis frequently desired in microfluidic experiments is counting and tracking the droplets. Specifically, droplet tracking in dense emulsions is challenging as droplets move in tightly packed configurations. Sometimes the individual droplets in these dense clusters are hard to resolve, even for a human observer. Here, two deep learning-based cutting-edge algorithms for object detection (YOLO) and object tracking (DeepSORT) are combined into a single image analysis tool, DropTrack, to track droplets in microfluidic experiments. DropTrack analyzes input videos, extracts droplets' trajectories, and infers other observables of interest, such as droplet numbers. Training an object detector network for droplet recognition with manually annotated images is a labor-intensive task and a persistent bottleneck. This work partly resolves this problem by training object detector networks (YOLOv5) with hybrid datasets containing real and synthetic images. We present an analysis of a double emulsion experiment as a case study to measure DropTrack's performance. For our test case, the YOLO networks trained with 60% synthetic images show similar performance in droplet counting as with the one trained using 100% real images, meanwhile saving the image annotation work by 60%. DropTrack's performance is measured in terms of mean average precision (mAP), mean square error in counting the droplets, and inference speed. The fastest configuration of DropTrack runs inference at about 30 frames per second, well within the standards for real-time image analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine learning assisted droplet trajectories extraction in dense emulsions and their analysis

This work analyzes trajectories obtained by YOLO and DeepSORT algorithms of dense emulsion systems simulated by Lattice Boltzmann methods. The results indicate that the individual droplet's moving direction is influenced more by the droplets immediately behind it than the droplets in front of it. The analysis also provides hints on constraints on writing down a dynamical model of droplets for the dense emulsion in narrow channels.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predicting today's cosmological constant via the Zel'dovich-Holographic connection

This Letter proposes a solution of the Vacuum Energy and the Cosmological Constant (CC) paradox based on the Zel'dovich's ansatz, which states that the observable contribution to the vacuum energy density is given by the gravitational energy of virtual particle-antiparticle pairs, continually generated and annihilated in the vacuum state. The novelty of this work is the use of an ultraviolet cut-off length based on the Holographic Principle, which is shown to yield current values of the CC in good agreement with experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Role of Oxygen Functionalities in Graphene Oxide Architectural Laminate Subnanometer Spacing and Water Transport

Active research in nanotechnology contemplates the use of nanomaterials for engineering applications. However, a primary challenge is understanding the effects of nanomaterial properties on industrial device performance and translating unique nanoscale properties to the macroscale. One emerging example is graphene oxide (GO) membranes for separation processes. Thus, here we investigate how individual GO properties can impact layered GO characteristics and water permeability. GO chemistry and morphology were controlled with easy-to-implement photo-reduction and sonication techniques and were quantitatively correlated offering a valuable tool to speed up the characterization process. For example, one could perform chemical analysis and concurrently obtain morphology information or vice versa. Chemical GO modification allows for fine control of GO oxidation state and GO laminate nanoarchitecture enabling controlled synthesis of a GO architectural laminate (GOAL). The GOAL can be considered as the selective layer of the membrane created by interconnected sub-nanometer channels, characterized by a length and a height (i.e., GO spacing), through which water molecules permeate. Water permeability was measured for eight GOAL characterized by different GO chemistry and morphology, and indicate that GO nanochannel height dictates water transport. The simulations indicate a no-slip Darcy-like water transport regime inside the GOAL due to the presence of basal oxygen functionalities. The experimental and simulation evidence presented in this letter help create a clearer picture of water transport in GO and can be used to rationally design more effective and efficient GO membranes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stochastic jetting and dripping in confined soft granular flows

We report new dynamical modes in confined soft granular flows, such as stochastic jetting and dripping, with no counterpart in continuum viscous fluids. The new modes emerge as a result of the propagation of the chaotic behaviour of individual grains -- here, monodisperse emulsion droplets to the level of the entire system as the emulsion is focused into a narrow orifice by an external viscous flow. We observe avalanching dynamics and the formation of remarkably stable jets -- singlefile granular chains -- which occasionally break, resulting in a non-Gaussian distribution of cluster sizes. We find that the sequences of droplet rearrangements that lead to the formation of such chains resemble unfolding of cancer cell clusters in narrow capillaries, overall demonstrating that microfluidic emulsion systems could serve to model various aspects of soft granular flows, including also tissue dynamics at the meso-scale.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tricritical point in the quantum Hamiltonian mean-field model

Engineering long-range interactions in experimental platforms has been achieved with great success in a large variety of quantum systems in recent years. Inspired by this progress, we propose a generalization of the classical Hamiltonian mean-field model to fermionic particles. We study the phase diagram and thermodynamic properties of the model in the canonical ensemble for ferromagnetic interactions as a function of temperature and hopping. At zero temperature, small charge fluctuations drive the many-body system through a first order quantum phase transition from an ordered to a disordered phase at zero temperature. At higher temperatures, the fluctuation-induced phase transition remains first order initially and switches to second order only at a tricritical point. Our results offer an intriguing example of tricriticality in a quantum system with long-range couplings, which bears direct experimental relevance. The analysis is performed by exact diagonalization and mean-field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Lattice Boltzmann Method for Relativistic Rarefied Flows in (2 + 1) Dimensions

We propose an extension to recently developed Relativistic Lattice Boltzmann solvers (RLBM), which allows the simulation of flows close to the free streaming limit. Following previous works [Phys. Rev. C 98 (2018) 035201], we use product quadrature rules and select weights and nodes by separately discretising the radial and the angular components. This procedure facilitates the development of quadrature-based RLBM with increased isotropy levels, thus improving the accuracy of the method for the simulation of flows beyond the hydrodynamic regime. In order to quantify the improvement of this discretisation procedure over existing methods, we perform numerical tests of shock waves in one and two spatial dimensions in various kinetic regimes across the hydrodynamic and the free-streaming limits.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

LBcuda: a high-performance CUDA port of LBsoft for simulation of colloidal systems

We present LBcuda, a GPU accelerated version of LBsoft, our open-source MPI-based software for the simulation of multi-component colloidal flows. We describe the design principles, the optimization and the resulting performance as compared to the CPU version, using both an average cost GPU and high-end NVidia GPU cards (V100 and the latest A100). The results show a substantial acceleration for the fluid solver reaching up to 200 GLUPS (Giga Lattice Updates Per Second) on a cluster made of 512 A100 NVIDIA cards simulating a grid of eight billion lattice points. These results open attractive prospects for the computational design of new materials based on colloidal particles.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Wet to dry self-transitions in dense emulsions: from order to disorder and back

One of the most distinctive hallmarks of many-body systems far from equilibrium is the spontaneous emergence of order under conditions where disorder would be plausibly expected. Here, we report on the self-transition between ordered and disordered emulsions in divergent microfluidic channels, i.e. from monodisperse assemblies to heterogeneous polydisperse foam-like structures, and back again to ordered ones. The transition is driven by the nonlinear competition between viscous dissipation and surface tension forces as controlled by the device geometry, particularly the aperture angle of the divergent microfluidic channel. An unexpected route back to order is observed in the regime of large opening angles, where a trend towards increasing disorder would be intuitively expected.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Coupled lattice Boltzmann-Multiparticle collision method for multi-resolution hydrodynamics

In this work we discuss the coupling of two mesoscopic approaches for fluid dynamics, namely the lattice Boltzmann method (LB) and the multiparticle collision dynamics (MPCD) \cite{kapral2008multiparticle} to design a new class of flexible and efficient multiscale schemes based on a dual representation of the fluid observables. At variance with other commonly used multigrid methods, mostly oriented to high Reynolds and turbulent flows, the present approach is designed to capture the physics at the smallest scales whenever the lattice Boltzmann alone falls short of providing the correct physical information due to a lack of resolution, as it occurs for example in thin films between interacting bubbles or droplets in microfluidic crystals. The coupling strategies employed to pass the hydrodynamic information between the LB and the MPCD are detailed and the algorithm is tested against the steady isothermal Poiseuille flow problem for different coupling configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A moving-grid approach for fluid-structure interaction problems with hybrid lattice Boltzmann method

In this paper, we propose a hybrid lattice Boltzmann method (HLBM) for solving fluid-structure interaction problems. The proposed numerical approach is applied to model the flow induced by a vibrating thin lamina submerged in a viscous quiescent fluid. The hydrodynamic force exerted by the fluid on the solid body is described by means of a complex hydrodynamic function, whose real and imaginary parts are determined via parametric analysis. Numerical results are validated by comparison with those from other numerical as well as experimental works available in the literature. The proposed hybrid approach enhances the capability of lattice Boltzmann methods to solve fluid dynamic problems involving moving geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multi-component discrete Boltzmann model for nonequilibrium reactive flows

We propose a multi-component discrete Boltzmann model (DBM) for premixed, nonpremixed, or partially premixed nonequilibrium reactive flows. This model is suitable for both subsonic and supersonic flows with or without chemical reaction and/or external force. A two-dimensional sixteen-velocity model is constructed for the DBM. In the hydrodynamic limit, the DBM recovers the modified Navier-Stokes equations for reacting species in a force field. Compared to standard lattice Boltzmann models, the DBM presents not only more accurate hydrodynamic quantities, but also detailed nonequilibrium effects that are essential yet long-neglected by traditional fluid dynamics. Apart from nonequilibrium terms (viscous stress and heat flux) in conventional models, specific hydrodynamic and thermodynamic nonequilibrium quantities (high order kinetic moments and their departure from equilibrium) are dynamically obtained from the DBM in a straightforward way. Due to its generality, the developed methodology is applicable to a wide range of phenomena across many energy technologies, emissions reduction, environmental protection, mining accident prevention, chemical and process industry.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Multiresolution Mesoscale Approach for Microscale Hydrodynamics

A new class of multiscale scheme is presented for micro-hydrodynamic problems based on a dual representation of the fluid observables. The hybrid model is first tested against the classical flow between two parallel plates and then applied to a plug flow within a micrometer-sized striction and a shear flow within a microcavity. Both cases demonstrate the capability of the multiscale approach to reproduce the correct macroscopic hydrodynamics also in the presence of refined grids (one and two levels), while retaining the correct thermal fluctuations, embedded in the multiparticle collision method. This provides the first step toward a novel class of fully mesoscale hybrid approaches able to capture the physics of fluids at the micro- and nanoscales whenever a continuum representation of the fluid falls short of providing the complete physical information, due to a lack of resolution and thermal fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Combined effects of fluid type and particle shape on particles flow in microfluidic platforms

Recent numerical analyses to optimize the design of microfluidic devices for more effective entrapment or segregation of surrogate circulating tumor cells (CTCs) from healthy cells have been reported in the literature without concurrently accommodating the non-Newtonian nature of the body fluid and the non-uniform geometric shapes of the CTCs. Through a series of two-dimensional proof-of-concept simulations with increased levels of complexity (e.g., number of particles, inline obstacles), we investigated the validity of the assumptions of the Newtonian fluid behavior for pseudoplastic fluids and the circular particle shape for different-shaped particles (DSPs) in the context of microfluidics-facilitated shape-based segregation of particles. Simulations with a single DSP revealed that even in the absence of internal geometric complexities of a microfluidics channel, the aforementioned assumptions led to 0.11-0.21W (W is the channel length) errors in lateral displacements of DSPs, up to 3-20% errors in their velocities, and 3-5% errors in their travel times. When these assumptions were applied in simulations involving multiple DSPs in inertial microfluidics with inline obstacles, errors in the lateral displacements of DSPs were as high as 0.78W and in their travel times up to 23%, which led to different (un)symmetric flow and segregation patterns of DSPs. Thus, the fluid type and particle shape should be included in numerical models and experiments to assess the performance of microfluidics for targeted cell (e.g., CTCs) harvesting.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamic symmetry-breaking in mutually annihilating fluids with selective interfaces

The selective entrapment of mutually annihilating species within a phase-changing carrier fluid is explored by both analytical and numerical means. The model takes full account of the dynamic heterogeneity which arises as a result of the coupling between hydrodynamic transport, dynamic phase-transitions and chemical reactions between the participating species, in the presence of a selective droplet interface. Special attention is paid to the dynamic symmetry breaking between the mass of the two species entrapped within the expanding droplet as a function of time. It is found that selective sources are much more effective symmetry breakers than selective diffusion. The present study may be of interest for a broad variety of advection-diffusion-reaction phenomena with selective fluid interfaces, including the problem of electroweak baryogenesis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Elucidating the mechanism of step-emulsification

Three-dimensional, time-dependent direct simulations of step emulsification micro-devices highlight two essential mechanisms for droplet formation: first, the onset of an adverse pressure gradient driving a back-flow of the continuous phase from the external reservoir to the micro-channel. Second, the striction of the flowing jet which leads to its subsequent rupture. It is also shown that such a rupture is delayed and eventually suppressed by increasing the flow speed of the dispersed phase within the channel, due to the stabilising effect of dynamic pressure. This suggests a new criterion for dripping-jetting transition, based on local values of the Capillary and Weber numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enhanced computational performance of the lattice Boltzmann model for simulating micron- and submicron-size particle flows and non-Newtonian fluid flows

Significant improvements in the computational performance of the lattice-Boltzmann (LB) model, coded in FORTRAN90, were achieved through application of enhancement techniques. Applied techniques include optimization of array memory layouts, data structure simplification, random number generation outside the simulation thread(s), code parallelization via OpenMP, and intra- and inter-timestep task pipelining. Effectiveness of these optimization techniques was measured on three benchmark problems: (i) transient flow of multiple particles in a Newtonian fluid in a heterogeneous fractured porous domain, (ii) thermal fluctuation of the fluid at the sub-micron scale and the resultant Brownian motion of a particle, and (iii) non-Newtonian fluid flow in a smooth-walled channel. Application of the aforementioned optimization techniques resulted in an average 21 performance improvement, which could significantly enhance practical uses of the LB models in diverse applications, focusing on the fate and transport of nano-size or micron-size particles in non-Newtonian fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluid flow around NACA 0012 airfoil at low-Reynolds numbers with hybrid lattice Boltzmann method

We simulate the two-dimensional fluid flow around National Advisory Committee for Aeronautics (NACA) 0012 airfoil using a hybrid lattice Boltzmann method (HLBM), which combines the standard lattice Boltzmann method with an unstructured finite-volume formulation. The aim of the study is to assess the numerical performances and the robustness of the computational method. To this purpose, after providing a convergence study to estimate the overall accuracy of the method, we analyze the numerical solution for different values of the angle of attack at a Reynolds number equal to $10^3$. Subsequently, flow fields at Reynolds numbers up to $10^4$ are computed for a zero angle of attack configuration. A grid refinement scheme is applied to the uniformly spaced component of the overlapping grid system to further enhance the numerical efficiency of the model. The results demonstrate the capability of the HLBM to achieve high accuracy near solid curved walls, thus providing a viable alternative in the realm of off-lattice Boltzmann methods based on body-fitted mesh.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

LBsoft: a parallel open-source software for simulation of colloidal systems

We present LBsoft, an open-source software developed mainly to simulate the hydro-dynamics of colloidal systems based on the concurrent coupling between lattice Boltzmann methods for the fluid and discrete particle dynamics for the colloids. Such coupling has been developed before, but, to the best of our knowledge, no detailed discussion of the programming issues to be faced in order to attain efficient implementation on parallel architectures, has ever been presented to date. In this paper, we describe in detail the underlying multi-scale models, their coupling procedure, along side with a description of the relevant input variables, to facilitate third-parties usage. The code is designed to exploit parallel computing platforms, taking advantage also of the recent AVX-512 instruction set. We focus on LBsoft structure, functionality, parallel implementation, performance and availability, so as to facilitate the access to this computational tool to the research community in the field. The capabilities of LBsoft are highlighted for a number of prototypical case studies, such as pickering emulsions, bicontinuous systems, as well as an original study of the coarsening process in confined bijels under shear.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mesoscale modeling of near-contact interactions for complex flowing interfaces

We present a mesoscale kinetic model for multicomponent flows, augmented with a short range forcing term, aimed at describing the combined effect of surface tension and near-contact interactions operating at the fluid interface level. Such mesoscale approach is shown to i) accurately capture the complex dynamics of bouncing colliding droplets for different values of the main governing parameters, ii) predict quantitatively the effective viscosity of dense emulsions in micro-channels and iii) simulate the formation of the so-called soft flowing crystals in microfluidic focusers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling pattern formation in soft flowing crystals

We present a mesoscale representation of near-contact interactions between colliding droplets which permits to reach up to the scale of full microfluidic devices, where such droplets are produced. The method is demonstrated for the case of colliding droplets and the formation of soft flowing crystals in flow-focussing microfluidic devices. This model may open up the possibility of multiscale simulation of microfluidic devices for the production of new droplet/bubble-based mesoscale porous materials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multiparticle collision dynamics for fluid interfaces with near-contact interactions

We present an extension of the multiparticle collision dynamics method for flows with complex interfaces, including supramolecular near-contact interactions mimicking the effect of surfactants. The new method is demonstrated for the case of (i) short range repulsion of droplets in close contact, (ii) arrested phase separation and (iii) different pattern formation during spinodal decomposition of binary mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Particle Shape Influences Settling and Sorting Behavior in Microfluidic Domains

We present a new numerical model to simulate settling trajectories of discretized individual or a mixture of particles of different geometrical shapes in a quiescent fluid and their flow trajectories in a flowing fluid. Simulations unveiled diverse particle settling trajectories as a function of their geometrical shape and density. The effects of the surface concavity of a boomerang particle and aspect ratio of a rectangular particle on the periodicity and amplitude of oscillations in their settling trajectories were numerically captured. Use of surrogate circular particles for settling or flowing of a mixture of non-circular particles were shown to miscalculate particle velocities by a factor of 0.9-2.2 and inaccurately determine the particles' trajectories. In a microfluidic chamber with particles of different shapes and sizes, simulations showed that steady vortices do not necessarily always control particle entrapments, nor do larger particles get selectively and consistently entrapped in steady vortices. Strikingly, a change in the shape of large particles from circular to elliptical resulted in stronger entrapments of smaller circular particles, but enhanced outflows of larger particles, which could be an alternative microfluidics-based method for sorting and separation of particles of different sizes and shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Semi-Lagrangian implicit Bhatnagar-Gross-Krook collision model for the finite-volume discrete Boltzmann method

A new implicit BGK collision model using a semi-Lagrangian approach is proposed in this paper. Unlike existing models, in which the implicit BGK collision is resolved either by a temporal extrapolation or by a variable transformation, the new model removes the implicitness by tracing the particle distribution functions (PDFs) back in time along their characteristic paths during the collision process. An interpolation scheme is needed to evaluate the PDFs at the traced-back locations. By using the first-order interpolation, the resulting model allows for the straightforward replacement of ${f_α}^{eq,n+1}$ by ${f_α}^{eq,n}$ no matter where it appears. After comparing the new model with the existing models under different numerical conditions (e.g. different flux schemes and time marching schemes) and using the new model to successfully modify the variable transformation technique, three conclusions can be drawn. First, the new model can improve the accuracy by almost an order of magnitude. Second, it can slightly reduce the computational cost. Therefore, the new scheme improves accuracy without extra cost. Finally, the new model can significantly improve the $Δt/τ$ limit compared to the temporal interpolation model while having the same $Δt/τ$ limit as the variable transformation approach. The new scheme with a second-order interpolation is also developed and tested; however, that technique displays no advantage over the simple first-order interpolation approach. Both numerical and theoretical analyses are also provided to explain why the new implicit scheme with simple first-order interpolation can outperform the same scheme with second-order interpolation, as well as the existing temporal extrapolation and variable transformation schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards a self-consistent Boltzmann's kinetic model of fluid turbulence

A closure for the effective relaxation time of the Boltzmann-BGK kinetic equation for fluid turbulence is presented, based on a double-averaging procedure over both kinetic and turbulent fluctuations. The resulting effective relaxation time appears to agree with values obtained via a renormalization group treatment of the Navier-Stokes equation only at low values of $k/T$, the ratio of turbulent kinetic energy to fluid temperature. For $k/T>0.1$ the kinetic treatment delivers a significantly longer effective relaxation time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards Exascale Lattice Boltzmann computing

We discuss the state of art of Lattice Boltzmann (LB) computing, with special focus on prospective LB schemes capable of meeting the forthcoming Exascale challenge. After reviewing the basic notions of LB computing, we discuss current techniques to improve the performance of LB codes on parallel machines and illustrate selected leading-edge applications in the Petascale range. Finally, we put forward a few ideas on how to improve the communication/computation overlap in current large-scale LB simulations, as well as possible strategies towards fault-tolerant LB schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Disordered interfaces in soft fluids with suspended colloids

Computer simulations of bi-continuous two-phase fluids with intersparsed dumbbells show that, unlike rigid colloids, soft dumbbells do not lead to arrested coarsening. However, they significantly alter the curvature dynamics of the fluid-fluid interface, whose probability density distributions are shown to exhibit (i) a universal spontaneous transition from an initial broad-shape distribution towards a highly localized one and (ii) super-diffusive dynamics with long-range effects. Both features may prove useful for the design of novel families of soft porous materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamic mesh refinement for discrete models of jet electro-hydrodynamics

Nowadays, several models of unidimensional fluid jets exploit discrete element methods. In some cases, as for models aiming at describing the electrospinning nanofabrication process of polymer fibers, discrete element methods suffer a non constant resolution of the jet representation. We develop a dynamic mesh-refinement method for the numerical study of the electro-hydrodynamic behavior of charged jets using discrete element methods. To this purpose, we import ideas and techniques from the string method originally developed in the framework of free-energy landscape simulations. The mesh-refined discrete element method is demonstrated for the case of electrospinning applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mapping Reactive Flow Patterns in Monolithic Nanoporous Catalysts

The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nano-scale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domain. To address this problem, we developed a multi-scale mass transport computational framework based on the Lattice Boltzmann Method (LBM) that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical evidence of Sinai diffusion of random-mass Dirac particles

We present quantum Lattice Boltzmann simulations of the Dirac equation for quantum-relativistic particles with random mass. By choosing zero-average random mass fluctuation, the simulations show evidence of localization and ultra-slow Sinai diffusion, due to the interference of oppositely propagating branches of the quantum wavefunction which result from random sign changes of the mass around a zero-mean. The present results indicate that the quantum lattice Boltzmann scheme may offer a viable tool for the numerical simulation of quantum-relativistic transport phenomena in topological materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical solution of the non-linear Schrodinger equation using smoothed-particle hydrodynamics

We formulate a smoothed-particle hydrodynamics numerical method, traditionally used for the Euler equations for fluid dynamics in the context of astrophysical simulations, to solve the non-linear Schrodinger equation in the Madelung formulation. The probability density of the wavefunction is discretized into moving particles, whose properties are smoothed by a kernel function. The traditional fluid pressure is replaced by a quantum pressure tensor, for which a novel, robust discretization is found. We demonstrate our numerical method on a variety of numerical test problems involving the simple harmonic oscillator, Bose-Einstein condensates, collapsing singularities, and dark matter halos governed by the Gross-Pitaevskii-Poisson equation. Our method is conservative, applicable to unbounded domains, and is automatically adaptive in its resolution, making it well suited to study problems with collapsing solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Striated Populations in Disordered Environments with Advection

Growth in static and controlled environments such as a Petri dish can be used to study the spatial population dynamics of microorganisms. However, natural populations such as marine microbes experience fluid advection and often grow up in heterogeneous environments. We investigate a generalized Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation describing single species population subject to a constant flow field and quenched random spatially inhomogeneous growth rates with a fertile overall growth condition. We analytically and numerically demonstrate that the non-equilibrium steady-state population density develops a flow-driven striation pattern. The striations are highly asymmetric with a longitudinal correlation length that diverges linearly with the flow speed and a transverse correlation length that approaches a finite velocity-independent value. Linear response theory is developed to study the statistics of the steady states. Theoretical predictions show excellent agreement with the numerical steady states of the generalized FKPP equation obtained from Lattice Boltzmann simulations. These findings suggest that, although the growth disorder can be spatially uncorrelated, correlated population structures with striations emerge naturally at sufficiently strong advection.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermo-hydrodynamic non-equilibrium effects on compressible Rayleigh-Taylor instability

The effects of compressibility on Rayleigh-Taylor instability (RTI) are investigated by inspecting the interplay between thermodynamic and hydrodynamic non-equilibrium phenomena (TNE, HNE, respectively) via a discrete Boltzmann model (DBM). Two effective approaches are presented, one tracking the evolution of the \emph{local} TNE effects and the other focussing on the evolution of the mean temperature of the fluid, to track the complex interfaces separating the bubble and the spike regions of the flow. It is found that, both the compressibility effects and the \emph{global} TNE intensity show opposite trends in the initial and the later stages of the RTI. Compressibility delays the initial stage of RTI and accelerates the later stage. Meanwhile, the TNE characteristics are generally enhanced by the compressibility, especially in the later stage. The global or mean thermodynamic non-equilibrium indicators provide physical criteria to discriminate between the two stages of the RTI.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Three-Dimensional Model for Electrospinning Processes in Controlled Gas Counterflow

We study the effects of a controlled gas flow on the dynamics of electrified jets in the electrospinning process. The main idea is to model the air drag effects of the gas flow by using a non-linear Langevin-like approach. The model is employed to investigate the dynamics of electrified polymer jets at different conditions of air drag force, showing that a controlled gas counterflow can lead to a decrease of the average diameter of electrospun fibers, and potentially to an improvement of the quality of electrospun products. We probe the influence of air drag effects on the bending instabilities of the jet and on its angular fluctuations during the process. The insights provided by this study might prove useful for the design of future electrospinning experiments and polymer nanofiber materials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

JETSPIN: a specific-purpose open-source software for simulations of nanofiber electrospinning

We present the open-source computer program JETSPIN, specifically designed to simulate the electrospinning process of nanofibers. Its capabilities are shown with proper reference to the underlying model, as well as a description of the relevant input variables and associated test-case simulations. The various interactions included in the electrospinning model implemented in JETSPIN are discussed in detail. The code is designed to exploit different computational architectures, from single to parallel processor workstations. This paper provides an overview of JETSPIN, focusing primarily on its structure, parallel implementations, functionality, performance, and availability.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-linear Langevin model for the early-stage dynamics of electrospinning jets

We present a non-linear Langevin model to investigate the early-stage dynamics of electrified polymer jets in electrospinning experiments. In particular, we study the effects of air drag force on the uniaxial elongation of the charged jet, right after ejection from the nozzle. Numerical simulations show that the elongation of the jet filament close to the injection point is significantly affected by the non-linear drag exerted by the surrounding air. These result provide useful insights for the optimal design of current and future electrospinning experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Lattice Boltzmann is a quantum walk

Numerical methods for the 1-D Dirac equation based on operator splitting and on the quantum lattice Boltzmann (QLB) schemes are reviewed. It is shown that these discretizations fall within the class of quantum walks, i.e. discrete maps for complex fields, whose continuum limit delivers Dirac-like relativistic quantum wave equations. The correspondence between the quantum walk dynamics and these numerical schemes is given explicitly, allowing a connection between quantum computations, numerical analysis and lattice Boltzmann methods. The QLB method is then extended to the Dirac equation in curved spaces and it is demonstrated that the quantum walk structure is preserved. Finally, it is argued that the existence of this link between the discretized Dirac equation and quantum walks may be employed to simulate relativistic quantum dynamics on quantum computers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rayleigh-Bénard Instability in Graphene

Motivated by the observation that electrons in graphene, in the hydrodynamic regime of transport, can be treated as a two-dimensional ultra-relativistic gas with very low shear viscosity, we examine the existence of the Rayleigh-Bénard instability in a massless electron-hole plasma. Firstly, we perform a linear stability analysis, derive the leading contributions to the relativistic Rayleigh number, and calculate the critical value above which the instability develops. By replacing typical values for graphene, such as thermal conductivity, shear viscosity, temperature, and sample sizes, we find that the instability might be experimentally observed in the near future. Additionally, we have performed simulations for vanishing reduced chemical potential and compare the measured critical Rayleigh number with the theoretical prediction, finding good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Semi-Spectral Method for the Wigner equation

We propose a numerical method to solve the Wigner equation in quantum systems of spinless, non-relativistic particles. The method uses a spectral decomposition into $L^2(\mathbb{R}^d)$ basis functions in momentum-space to obtain a system of first-order advection-reaction equations. The resulting equations are solved by splitting the reaction and advection steps so as to allow the combination of numerical techniques from quantum mechanics and computational fluid dynamics by identifying the skew-hermitian reaction matrix as a generator of unitary rotations. The method is validated for the case of particles subject to a one-dimensional (an-)harmonic potential using finite-differences for the advection part. Thereby, we verify the second order of convergence and observe non-classical behavior in the evolution of the Wigner function.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sub-ms dynamics of the instability onset of electrospinning

Electrospun polymer jets are imaged for the first time at an ultra-high rate of 10,000 frames per second, investigating the process dynamics, and the instability propagation velocity and displacement in space. The polymer concentration, applied voltage bias and needle-collector distance are systematically varied, and their influence on the instability propagation velocity and on the jet angular fluctuations analyzed. This allows us to unveil the instability formation and cycling behavior, and its exponential growth at the onset, exhibiting radial growth rates of the order of 10^3 s^-1. Allowing the conformation and evolution of polymeric solutions to be studied in depth, high-speed imaging at sub-ms scale shows a significant potential for improving the fundamental knowledge of electrified jets, leading to obtain finely controllable bending and solution stretching in electrospinning, and consequently better designed nanofibers morphologies and structures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tailoring boundary geometry to optimize heat transport in turbulent convection

By tailoring the geometry of the upper boundary in turbulent Rayleigh-Bénard convection we manipulate the boundary layer -- interior flow interaction, and examine the heat transport using the Lattice Boltzmann method. For fixed amplitude and varying boundary wavelength $λ$, we find that the exponent $β$ in the Nusselt-Rayleigh scaling relation, $Nu-1 \propto Ra^β$, is maximized at $λ\equiv λ_{\text{max}} \approx (2 π)^{-1}$, but decays to the planar value in both the large ($λ\gg λ_{\text{max}}$) and small ($λ\ll λ_{\text{max}}$) wavelength limits. The changes in the exponent originate in the nature of the coupling between the boundary layer and the interior flow. We present a simple scaling argument embodying this coupling, which describes the maximal convective heat flux.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accelerated kinetic over-relaxation for dynamic energy minimization in lattice many-body problems

A kinetic over-relaxation minimizer, based on a lattice version of the Boltzmann equation is presented. The method is validated against standard Metropolis Monte Carlo, and proves very effective in attaining (global) minima of classical pair potentials, involving solid body rotations. Linear scaling of the computational time to minimum with the system size is demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dark energy from cosmological fluids obeying a Shan-Chen non-ideal equation of state

We consider a Friedmann-Robertson-Walker universe with a fluid source obeying a non-ideal equation of state with "asymptotic freedom," namely ideal gas behavior (pressure changes directly proportional to density changes) both at low and high density regimes, following a fluid dynamical model due to Shan and Chen. It is shown that, starting from an ordinary energy density component, such fluids naturally evolve towards a universe with a substantial "dark energy" component at the present time, with no need of invoking any cosmological constant. Moreover, we introduce a quantitative indicator of darkness abundance, which provides a consistent picture of the actual matter-energy content of the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of non-linear rheology on the electrospinning process: a model study

We develop an analytical bead-spring model to investigate the role of non-linear rheology on the dynamics of electrified jets in the early stage of the electrospinning process. Qualitative arguments, parameter studies as well as numerical simulations, show that the elongation of the charged jet filament is significantly reduced in the presence of a non-zero yield stress. This may have beneficial implications for the optimal design of future electrospinning experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Friction forces in cosmological models

We investigate the dynamics of test particles undergoing friction forces in a Friedmann-Robertson-Walker (FRW) spacetime. The interaction with the background fluid is modeled by introducing a Poynting-Robertson-like friction force in the equations of motion, leading to measurable (at least in principle) deviations of the particle trajectories from geodesic motion. The effect on the peculiar velocities of the particles is investigated for various equations of state of the background fluid and different standard cosmological models. The friction force is found to have major effects on particle motion in closed FRW universes, where it turns the time-asymptotic value (approaching the recollapse) of the peculiar particle velocity from ultra-relativistic (close to light speed) to a co-moving one, i.e., zero peculiar speed. On the other hand, for open or flat universes the effect of the friction is not so significant, because the time-asymptotic peculiar particle speed is largely non-relativistic also in the geodesic case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Particle scattering by a test fluid on a Schwarzschild spacetime: the equation of state matters

The motion of a massive test particle in a Schwarzschild spacetime surrounded by a perfect fluid with equation of state $p_0= w ρ_0$ is investigated. Deviations from geodesic motion are analyzed as a function of the parameter $w$, ranging from $w=1$ which corresponds to the case of massive free scalar fields, down into the so-called "phantom" energy, with $w<-1$. It is found that the interaction with the fluid distribution leads to capture (escape) of the particle trajectory in the case $1+w>0$ ($<0$), respectively. Based on this result, it is argued that inspection of the trajectories of test particles in the vicinity of a Schwarzschild black hole may offer a new means of gaining insights into the nature of cosmic matter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scalar field inflation and Shan-Chen fluid models

A scalar field equivalent to a non-ideal "dark energy fluid" obeying a Shan-Chen-like equation of state is used as the background source of a flat Friedmann-Robertson-Walker cosmological spacetime to describe the inflationary epoch of our universe. Within the slow-roll approximation, a number of interesting features are presented, including the possibility to fulfill current observational constraints as well as a graceful exit mechanism from the inflationary epoch.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ultrathin fibres from electrospinning experiments under driven fast-oscillating perturbations

The effects of a driven fast-oscillating spinneret on the bending instability of electrified jets, leading to the formation of spiral structures in electrospinning experiments with charged polymers, are explored by means of extensive computer simulations. It is found that the morphology of the spirals can be placed in direct correspondence with the oscillation frequency and amplitude. In particular, by increasing the oscillation amplitude and frequency, thinner fibres can be extracted by the same polymer material, thereby opening design scenarios in electrospinning experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Z-index: A geometric representation of productivity and impact which accounts for information in the entire rank-citation profile

We present a simple generalization of Hirsch's h-index, Z = \sqrt{h^{2}+C}/\sqrt{5}, where C is the total number of citations. Z is aimed at correcting the potentially excessive penalty made by h on a scientist's highly cited papers, because for the majority of scientists analyzed, we find the excess citation fraction (C-h^{2})/C to be distributed closely around the value 0.75, meaning that 75 percent of the author's impact is neglected. Additionally, Z is less sensitive to local changes in a scientist's citation profile, namely perturbations which increase h while only marginally affecting C. Using real career data for 476 physicists careers and 488 biologist careers, we analyze both the distribution of $Z$ and the rank stability of Z with respect to the Hirsch index h and the Egghe index g. We analyze careers distributed across a wide range of total impact, including top-cited physicists and biologists for benchmark comparison. In practice, the Z-index requires the same information needed to calculate h and could be effortlessly incorporated within career profile databases, such as Google Scholar and ResearcherID. Because Z incorporates information from the entire publication profile while being more robust than h and g to local perturbations, we argue that Z is better suited for ranking comparisons in academic decision-making scenarios comprising a large number of scientists.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Isotropic discrete Laplacian operators from lattice hydrodynamics

We show that discrete schemes developed for lattice hydrodynamics provide an elegant and physically transparent way of deriving Laplacians with isotropic discretisation error. Isotropy is guaranteed whenever the Laplacian weights follow from the discrete Maxwell-Boltzmann equilibrium since these are, by construction, isotropic on the lattice. We also point out that stencils using as few as 15 points in three dimensions, generate isotropic Laplacians. These computationally efficient Laplacians can be used in cell-dynamical and hybrid lattice Boltzmann simulations, in favor of popular anisotropic Laplacians, which make use of larger stencils. The method can be extended to provide discretisations of higher order and for other differential operators, such the gradient, divergence and curl.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice operators from discrete hydrodynamics

We present a general scheme to derive lattice differential operators from the discrete velocities and associated Maxwell-Boltzmann distributions used in lattice hydrodynamics. Such discretizations offer built-in isotropy and recursive techniques to increase the convergence order. This provides a simple and elegant procedure to derive isotropic and accurate discretizations of differential operators, which are expected to apply across a broad range of problems in computational physics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lattice Boltzmann Method for Electromagnetic Wave Propagation

We present a new Lattice Boltzmann (LB) formulation to solve the Maxwell equations for electromagnetic (EM) waves propagating in a heterogeneous medium. By using a pseudo-vector discrete Boltzmann distribution, the scheme is shown to reproduce the continuum Maxwell equations. The technique compares well with a pseudo-spectral method at solving for two-dimensional wave propagation in a heterogeneous medium, which by design contains substantial contrasts in the refractive index. The extension to three dimensions follows naturally and, owing to the recognized efficiency of LB schemes for parallel computation in irregular geometries, it gives a powerful method to numerically simulate a wide range of problems involving EM wave propagation in complex media.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistical regularities in the rank-citation profile of scientists

Recent "science of science" research shows that scientific impact measures for journals and individual articles have quantifiable regularities across both time and discipline. However, little is known about the scientific impact distribution at the scale of an individual scientist. We analyze the aggregate scientific production and impact of individual careers using the rank-citation profile c_{i}(r) of 200 distinguished professors and 100 assistant professors. For the entire range of paper rank r, we fit each c_{i}(r) to a common distribution function that is parameterized by two scaling exponents. Since two scientists with equivalent Hirsch h-index can have significantly different c_{i}(r) profiles, our results demonstrate the utility of the β_{i} scaling parameter in conjunction with h_{i} for quantifying individual publication impact. We show that the total number of citations C_{i} tallied from a scientist's N_{i} papers scales as C_{i} \sim h_{i}^{1+β_{i}}. Such statistical regularities in the input-output patterns of scientists can be used as benchmarks for theoretical models of career progress.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards a mesoscopic model of water-like fluids with hydrodynamic interactions

We present a mesoscopic lattice model for non-ideal fluid flows with directional interactions, mimicking the effects of hydrogen-bonds in water. The model supports a rich and complex structural dynamics of the orientational order parameter, and exhibits the formation of disordered domains whose size and shape depend on the relative strength of directional order and thermal diffusivity. By letting the directional forces carry an inverse density dependence, the model is able to display a correlation between ordered domains and low density regions, reflecting the idea of water as a denser liquid in the disordered state than in the ordered one.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiple-Relaxation-Time Lattice Boltzmann Approach to Compressible Flows with Flexible Specific-Heat Ratio and Prandtl Number

A new multiple-relaxation-time lattice Boltzmann scheme for compressible flows with arbitrary specific heat ratio and Prandtl number is presented. In the new scheme, which is based on a two-dimensional 16-discrete-velocity model, the moment space and the corresponding transformation matrix are constructed according to the seven-moment relations associated with the local equilibrium distribution function. In the continuum limit, the model recovers the compressible Navier-Stokes equations with flexible specific-heat ratio and Prandtl number. Numerical experiments show that compressible flows with strong shocks can be simulated by the present model up to Mach numbers $Ma \sim 5$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Shear Banding from lattice kinetic models with competing interactions

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical simulation of conformational variability in biopolymer translocation through wide nanopores

Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores. The simulations provide clear evidence of folding quantization in the translocation process as the biopolymers undertake multi-folded configurations, characterized by a well-defined integer number of folds. Accordingly, the translocation time is no longer represented by a single-exponent power law dependence on the length, as it is the case for single-file translocation through narrow pores. The folding quantization increases with the biopolymer length, while the rate of translocated beads at each time step is linearly correlated to the number of resident beads in the pore. Finally, analysis of the statistics over the translocation work unravels the importance of the hydrodynamic interactions in the process.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Mesoscopic lattice Boltzmann modeling of flowing soft systems

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

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Multiscale modeling of biopolymer translocation through a nanopore

We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Langevin molecular dynamics (MD) with a mesoscopic lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. Both dynamical and statistical aspects of the translocation process were investigated, by simulating polymers of various initial configurations and lengths. For a representative molecule size, we explore the effects of important parameters that enter in the simulation, paying particular attention to the strength of the molecule-solvent coupling and of the external electric field which drives the translocation process. Finally, we explore the connection between the generic polymers modeled in the simulation and DNA, for which interesting recent experimental results are available.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Lattice Boltzmann Equation: Failure or Success?

The lattice Boltzmann equation (LBE) is a microscopically-inspired method designed to solve macroscopic fluid dynamics problems. As a such, it lives at the interface between the microscopic (molecular) and macroscopic (continuum) worlds, hopefully capturing the best of the two. In this paper we shall discuss whether or not, after almost a decade since its inception, LBE has lived up to the initial expectations. Open problems and future research developments are also briefly outlined.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Turbulent channel flow simulations using a coarse-grained extension of the Lattice Boltzmann method

A coarse-grained version of the Lattice Boltzmann (LB) method is developed with the intent of enhancing its geometrical flexibility so as to be able to tackle a wider class of flows of engineering interest. To this purpose, the original uniform LB technique is combined with standard finite-volume techniques based upon a blend of piecewise constant and piecewise linear interpolation schemes. A series of validation tests for the three dimensional channel flow with one-dimensional (cross-channel) statistical behaviour are presented. The main conclusion is that, although the method does indeed mark a significant stride forward with respect to the original uniform LB scheme, better interpolation schemes should be developed before the coarse-grain LB can become fully competitive with modern CFD schemes.

0 citations0 reviewsNo code linkedOpen access