BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Sauro Succi

Sauro Succi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phcond-mat.softquant-phComputer Visioncond-mat.stat-mechcond-mat.mes-hallcs.CYgr-qcphysics.ed-phphysics.hist-phphysics.soc-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
36works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

36 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