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 behavior 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 - single-file 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 mesoscale.

In this perspective we take stock of the current state of the art of computational models for droplets microfluidics and we suggest some strategies which may open the way to the full-scale simulation of microfluidic phenomena with interfaces, from near-contact interactions to the device operational lengths.

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. Program summary: Program Title: LBcuda CPC Library link to program files: https://doi.org/10.17632/v6fvmzpcrn.1 Developer's repository link: https://github.com/copmat/LBcuda Licensing provisions: 3-Clause BSD License Programming language: CUDA Fortran Nature of problem: Hydro-dynamics of colloidal multi-component systems and Pickering emulsions. Solution method: Lattice-Boltzmann method solving the Navier-Stokes equations for the fluid dynamics within an Eulerian description. Particle solver describing colloidal particles within a Lagrangian representation coupled to the fluid solver. The numerical solution of the coupling algorithm includes the back reaction effects for each force terms according to a fluid-particle multi-scale paradigm.

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 due to inherently small droplets moving 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 [you only look once (YOLO)] and object tracking (DeepSORT) are combined into a single image analysis tool, DropTrack, to track droplets in the microfluidic experiments. DropTrack analyzes input microfluidic experimental 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. In this work, this problem is partly resolved by training many object detector networks (YOLOv5) with several 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 network trained by combining 40% real images and 60% synthetic images yields the best accuracy in droplet detection and droplet counting in real experimental videos. Also, this strategy reduces labor-intensive image annotation work by 60%. DropTrack's performance is measured in terms of mean average precision of droplet detection, mean squared error in counting the droplets, and image analysis speed for inferring droplets' trajectories. The fastest configuration of DropTrack can detect and track the droplets at approximately 30 frames per second, well within the standards for a real-time image analysis.

Photocurable polymers are used ubiquitously in 3D printing, coatings, adhesives, and composite fillers. In the present work, the free radical polymerization of photocurable compounds is studied using reactive classical molecular dynamics combined with a dynamical approach of the nonequilibrium molecular dynamics (D-NEMD). Different concentrations of radicals and reaction velocities are considered. The mechanical properties of the polymer resulting from 1,6-hexanediol dimethacrylate systems are characterized in terms of viscosity, diffusion constant, and activation energy, whereas the topological ones through the number of cycles (polymer loops) and cyclomatic complexity. Effects like volume shrinkage and delaying of the gel point for increasing monomer concentration are also predicted, as well as the stress-strain curve and Young's modulus. Combining ab initio, reactive molecular dynamics, and the D-NEMD method might lead to a novel and powerful tool to describe photopolymerization processes and to original routes to optimize additive manufacturing methods relying on photosensitive macromolecular systems.

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

In this paper, we numerically investigate the breakup dynamics of droplets in an emulsion flowing in a tapered microchannel with a narrow constriction. The mesoscale approach for multicomponent fluids with near contact interactions is shown to capture the deformation and breakup dynamics of droplets interacting within the constriction, in agreement with experimental evidence. In addition, it permits us to investigate in detail the hydrodynamic phenomena occurring during breakup stages. Finally, a suitable deformation parameter is introduced and analyzed to characterize the state of deformation of the system by inspecting pairs of interacting droplets flowing in the narrow channel.

The phase separation of a two-dimensional active binary mixture is studied under the action of an applied shear through numerical simulations. It is highlighted how the strength of the external flow modifies the initial shape of growing domains. The activity is responsible for the formation of isolated droplets which affect both the coarsening dynamics and the morphology of the system. The characteristic dimensions of domains along the flow and the shear direction are modulated in time by oscillations whose amplitudes are reduced when the activity increases. This induces a broadening of the distribution functions of domain lengths with respect to the passive case due to the presence of dispersed droplets of different sizes.

A recently proposed mesoscale approach for the simulation of multicomponent flows with near-contact interactions is employed to investigate the early stage formation and clustering statistics of soft flowing crystals in microfluidic channels. Specifically, we first demonstrate the ability of the aforementioned mesoscale model to accurately reproduce main mechanisms leading to the formation of two basic droplet patterns (triangular and hexagonal), in close agreement with experimental evidence. Next, we quantitatively evaluate the device-scale clustering efficiency of the crystal formation process by introducing a new orientational order parameter, based on the Delaunay triangulation and Voronoi diagrams analysis of the droplet patterns. The mesoscale computational approach employed in this work proves to be an efficient tool to shed new light on the complex dynamics of dense emulsions, from short-scale thin-film hydrodynamics, all the way up to global structure formation and statistics of the resulting droplets ensembles.

A recently proposed mesoscale approach for the simulation of multicomponent flows with near-contact interactions is employed to investigate the early stage formation and clustering statistics of soft flowing crystals in microfluidic channels. Specifically, we first demonstrate the ability of the aforementioned mesoscale model to accurately reproduce main mechanisms leading to the formation of two basic droplet patterns (triangular and hexagonal), in close agreement with experimental evidence. Next, we quantitatively evaluate the device-scale clustering efficiency of the crystal formation process by introducing a new orientational order parameter, based on the Delaunay triangulation and Voronoi diagrams analysis of the droplet patterns. The mesoscale computational approach employed in this work proves to be an efficient tool to shed new light on the complex dynamics of dense emulsions, from short-scale thin-film hydrodynamics, all the way up to global structure formation and statistics of the resulting droplets ensembles.

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

This work presents a microscale approach for simulating the dielectrophoresis assembly of polarizable particles under an external electric field. The model is shown to capture interesting dynamical and topological features, such as the formation of chains of particles and their incipient aggregation into hierarchical structures. A quantitative characterization in terms of the number and size of these structures is also discussed. This computational model could represent a viable numerical tool to study the mechanical properties of particle-based hierarchical materials and suggest new strategies for enhancing their design and manufacture. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

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

We present a deep learning-based object detection and object tracking algorithm to study droplet motion in dense microfluidic emulsions. The deep learning procedure is shown to correctly predict the droplets' shape and track their motion at competitive rates as compared to standard clustering algorithms, even in the presence of significant deformations. The deep learning technique and tool developed in this work could be used for the general study of the dynamics of biological agents in fluid systems, such as moving cells and self-propelled microorganisms in complex biological flows. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

We numerically study the dynamics of a polydisperse double emulsion under a symmetric shear flow. We show that both dispersity and shear rate crucially affect the behavior of the innermost drops and of the surrounding shell. While at low/moderate values of shear rates, the inner drops rotate periodically around a common center of mass triggered by the fluid vortex formed within the emulsion generally regardless of their polydispersity; at higher values, such dynamics occurs only at increasing polydispersity, since monodisperse drops are found to align along the shear flow and become approximately motionless at late times. Our simulations also suggest that increasing polydispersity favors close-range contacts among cores and persistent collisions, while hindering shape deformations of the external droplet. A quantitative evaluation of these effects is also provided.

We numerically study the translocation dynamics of double emulsion drops with multiple close-packed inner droplets within constrictions. Such liquid architectures, which we refer to as HIPdEs (high-internal phase double emulsions), consist of a ternary fluid, in which monodisperse droplets are encapsulated within a larger drop in turn immersed in a bulk fluid. Extensive two-dimensional lattice Boltzmann simulations show that if the area fraction of the internal drops is close to the packing fraction limit of hard spheres and the height of the channel is much smaller than the typical size of the emulsion, the crossing yields permanent shape deformations persistent over long periods of time. Morphological changes and rheological response are quantitatively assessed in terms of the structure of the velocity field, circularity of the emulsion, and rate of energy dissipated by viscous forces. Our results may be used to improve the design of soft mesoscale porous materials, which employ HIPdEs as templates for tissue engineering applications.

The state-of-the-art deep learning-based object recognition YOLO algorithm and object tracking DeepSORT algorithm are combined to analyze digital images from fluid dynamic simulations of multi-core emulsions and soft flowing crystals and to track moving droplets within these complex flows. The YOLO network was trained to recognize the droplets with synthetically prepared data, thereby bypassing the labor-intensive data acquisition process. In both applications, the trained YOLO + DeepSORT procedure performs with high accuracy on the real data from the fluid simulations, with low error levels in the inferred trajectories of the droplets and independently computed ground truth. Moreover, using commonly used desktop GPUs, the developed application is capable of analyzing data at speeds that exceed the typical image acquisition rates of digital cameras (30 fps), opening the interesting prospect of realizing a low-cost and practical tool to study systems with many moving objects, mostly but not exclusively, biological ones. Besides its practical applications, the procedure presented here marks the first step towards the automatic extraction of effective equations of motion of many-body soft flowing systems.

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

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. 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) [20] to design a new class of flexible and efficient multiscale schemes based on a dual representation of the fluid observables.

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