Over the last decade, the Lattice Boltzmann method has found major scope for the simulation of a large spectrum of problems in soft matter, from multiphase and multi-component microfluidic flows, to foams, emulsions, colloidal flows, to name but a few. Crucial to many such applications is the role of supramolecular interactions which occur whenever mesoscale structures, such as bubbles or droplets, come in close contact, say of the order of tens of nanometers. Regardless of their specific physico-chemical origin, such near-contact interactions are vital to preserve the coherence of the mesoscale structures against coalescence phenomena promoted by capillarity and surface tension, hence the need of including them in Lattice Boltzmann schemes. Strictly speaking, this entails a complex multiscale problem, covering about six spatial decades, from centimeters down to tens of nanometers, and almost twice as many in time. Such a multiscale problem can hardly be taken by a single computational method, hence the need for coarse-grained models for the near-contact interactions. In this review, we shall discuss such coarse-grained models and illustrate their application to a variety of soft flowing matter problems, such as soft flowing crystals, strongly confined dense emulsions, flowing hierarchical emulsions, soft granular flows, as well as the transmigration of active droplets across constrictions. Finally, we conclude with a few considerations on future developments in the direction of quantum-nanofluidics, machine learning, and quantum computing for soft flows applications.

Droplet microfluidics has emerged as highly relevant technology in diverse fields such as nanomaterials synthesis, photonics, drug delivery, regenerative medicine, food science, cosmetics, and agriculture. While significant progress has been made in understanding the fundamental mechanisms underlying droplet generation in microchannels and in fabricating devices to produce droplets with varied functionality and high throughput, challenges persist along two important directions. On one side, the generalization of numerical results obtained by computational fluid dynamics would be important to deepen the comprehension of complex physical phenomena in droplet microfluidics, as well as the capability of predicting the device behavior. Conversely, truly three-dimensional architectures would enhance microfluidic platforms in terms of tailoring and enhancing droplet and flow properties. Recent advancements in artificial intelligence (AI) and additive manufacturing (AM) promise unequaled opportunities for simulating fluid behavior, precisely tracking individual droplets, and exploring innovative device designs. This review provides a comprehensive overview of recent progress in applying AI and AM to droplet microfluidics. The basic physical properties of multiphase flows and mechanisms for droplet production are discussed, and the current fabrication methods of related devices are introduced, together with their applications. Delving into the use of AI and AM technologies in droplet microfluidics, topics covered include AI-assisted simulations of droplet behavior, real-time tracking of droplets within microfluidic systems, and AM-fabrication of three-dimensional systems. The synergistic combination of AI and AM is expected to deepen the understanding of complex fluid dynamics and active matter behavior, expediting the transition toward fully digital microfluidic systems.

In this work an optimized multicomponent lattice Boltzmann (LB) model is deployed to simulate axisymmetric turbulent jets of a fluid evolving in a quiescent, immiscible environment over a wide range of dynamic regimes. The implementation of the multicomponent LB code achieves peak performance on graphic processing units (GPUs) with a significant reduction of the memory footprint, retains the algorithmic simplicity inherent to standard LB computing, and, being based on a high-order extension of the thread-safe LB algorithm, allows us to perform stable simulations at vanishingly low viscosities. The proposed approach opens attractive prospects for high-performance computing simulations of realistic turbulent flows with interfaces on GPU-based architectures.

We present a mathematical framework to include quantum interfacial interactions, provided by Keldysh nonequilibrium quantum transport formalism, bottom-up coupled to a nanoscale lattice Boltzmann method. As an applicative scenario, we simulate a two-dimensional water flow between two parallel solid plates hosting electrons and phonons in the solid bottom wall. The corresponding tool may prove useful for the computational design of quantum-engineered nanofluidic devices, showing its capability to explore the effects of the interfacial quantum transport phenomena at scales of experimental relevance.

In this work, we present accLB, a high-performance Fortran-based lattice Boltzmann (LB) solver tailored to multiphase turbulent flows on multi-GPU architectures. The code couples a conservative phase-field formulation of the Allen–Cahn equation with a thread-safe, regularized LB method to capture complex interface dynamics. Designed from the ground up for HPC environments, accLB employs MPI for distributed memory parallelism and OpenACC for GPU acceleration, achieving excellent portability and scalability on leading pre-exascale systems such as Leonardo and LUMI. Benchmark tests demonstrate strong and weak scaling efficiencies on multiple GPUs. Physical validation includes direct numerical simulations of homogeneous isotropic turbulence (HIT). Further, we examine bubble-laden HIT and observe a transition to a -3 energy scaling, as in experiments and theoretical predictions, due to bubble-induced dissipation, along with enhanced small-scale intermittency. These results highlight accLB as a robust and scalable platform for the simulation of multiphase turbulence in extreme computational regimes.

We introduce a regularized fluctuating lattice Boltzmann model (Reg-FLBM) for the D3Q27 lattice, which incorporates thermal fluctuations through Hermite-based projections to ensure compliance with the fluctuation-dissipation theorem. By leveraging the recursive regularization framework, the model achieves thermodynamic consistency for both hydrodynamic and ghost modes. Compared to the conventional single-relaxation-time BGK-FLBM, the Reg-FLBM provides improved stability and a more accurate description of thermal fluctuations. The implementation is optimized for large-scale parallel simulations on graphics processing unit-accelerated architectures, enabling systematic investigation of fluctuation-driven phenomena in mesoscale and nanoscale fluid systems.

This study presents a high-order, thread-safe version of the lattice Boltzmann method, incorporating an interface-capturing equation, based on the conservative Allen-Cahn equation, to simulate incompressible two-component systems with high-density and viscosity contrasts. The method utilizes a recently proposed thread-safe implementation optimized for shared-memory architectures, and it is employed to reproduce the dynamics of droplets and bubbles in several test cases with results in agreement with experiments and other numerical simulations from the literature. The proposed approach offers promising opportunities for high-performance computing simulations of realistic fluid systems with high-density and viscosity contrasts for advanced applications in environmental, atmospheric, and meteorological flows, all the way down to microfluidic and biological systems, particularly on graphic processing unit-based architectures.

Accurate prediction of rarefied gas dynamics is crucial for optimizing flows through microelectromechanical systems, air filtration devices, and shale gas extraction. Traditional methods, such as discrete velocity and direct simulation Monte Carlo (DSMC), demand intensive memory and computation, especially for microflows in non-convex domains. Recently, physics-informed neural networks (PINNs) emerged as a meshless and adaptable alternative for solving non-linear partial differential equations. We trained a PINN using a limited number of DSMC-generated rarefied gas microflows in the transition regime (0.1<3), incorporating continuity and Cauchy momentum exchange equations in the loss function. The PINN achieved under 2 % error on these residuals and effectively filtered DSMC's intrinsic statistical noise. Predictions remained strong for a tested flow field with Kn=0.7, and showed limited extrapolation performance on a flow field with Kn=5 with a local overshoot of about 20 %, while maintaining physical consistency. Notably, each DSMC field required ∼20 hours on 4 graphics processing units (GPU), while the PINN training took <2 hours on one GPU, with evaluations under 2 seconds.

In microfluidic systems, droplets undergo intricate deformations as they traverse flow-focusing junctions, posing a challenging task for accurate measurement, especially during short transit times. This study investigates the physical behavior of droplets within dense emulsions in diverse microchannel geometries, specifically focusing on the impact of varying opening angles within the primary channel and injection rates of fluid components. Employing a sophisticated droplet tracking tool based on deep-learning techniques, we analyze multiple frames from flow-focusing experiments to quantitatively characterize droplet deformation in terms of ratio between maximum width and height and propensity to form liquid with hexagonal spatial arrangement. Our findings reveal the existence of an optimal opening angle where shape deformations are minimal and hexagonal arrangement is maximal. Variations of fluid injection rates are also found to affect size and packing fraction of the emulsion in the exit channel. This paper offers insight into deformations, size, and structure of fluid emulsions relative to microchannel geometry and other flow-related parameters captured through machine learning, with potential implications for the design of microchips utilized in cellular transport and tissue engineering applications.

The shape of liquid droplets in air plays an important role in the aerodynamic behavior and combustion dynamics of miniaturized propulsion systems such as microsatellites and small drones. Their precise manipulation can yield optimal efficiency in such systems. It is desired to have a minimal representation of droplet shapes using as few parameters as possible to automate shape manipulation using self-learning algorithms, such as reinforcement learning. In this paper, we use a neural compression algorithm to represent, with only two parameters, elliptical and bullet-shaped droplets initially represented with 200 points (400 real numbers) at the droplet boundary. The mapping of many to two points is achieved in two stages. Initially, a Fourier series is formulated to approximate the contour of the droplet. Subsequently, the coefficients of this Fourier series are condensed to lower dimensions utilizing a neural network with a bottleneck architecture. Finally, 5000 synthetically generated droplet shapes were used to train the neural network. With a two-real-number representation, the recovered droplet shapes had excellent overlap with the original ones, with a mean square error of ∼10−3 . Hence, this method compresses the droplet contour to merely two numerical parameters via a fully reversible process, a crucial feature for rendering learning algorithms computationally tractable.

We introduce a two-step, fully reversible process designed to project the outer shape of a generic droplet onto a lower-dimensional space. The initial step involves representing the droplet's shape as a Fourier series. Subsequently, the Fourier coefficients are reduced to lower-dimensional vectors by using autoencoder models. The exploitation of the domain knowledge of the droplet shapes allows us to map generic droplet shapes to just two-dimensional (2D) space in contrast to previous direct methods involving autoencoders that could map it on minimum eight-dimensional (8D) space. This six-dimensional (6D) reduction in the dimensionality of the droplet's description opens new possibilities for applications, such as automated droplet generation via reinforcement learning, the analysis of droplet shape evolution dynamics, and the prediction of droplet breakup. Our findings underscore the benefits of incorporating domain knowledge into autoencoder models, highlighting the potential for increased accuracy in various other scientific disciplines.

We present a highly optimized thread-safe lattice Boltzmann model in which the non-equilibrium part of the distribution function is locally reconstructed via recursivity of Hermite polynomials. Such a procedure allows the explicit incorporation of non-equilibrium moments of the distribution up to the order supported by the lattice. Thus, the proposed approach increases accuracy and stability at low viscosities without compromising performance and amenability to parallelization with respect to standard lattice Boltzmann models. The high-order thread-safe lattice Boltzmann is tested on two types of turbulent flows, namely, the turbulent channel flow at R e τ = 180 and the axisymmetric turbulent jet at Re = 7000; it delivers results in excellent agreement with reference data [direct numerical simulations (DNS), theory, and experiments] and (a) achieves peak performance [ ∼ 5 × 10 12 floating point operations (FLOP) per second and an arithmetic intensity of ∼ 7 FLOP / byte on a single graphic processing unit] by significantly reducing the memory footprint, (b) retains the algorithmic simplicity of standard lattice Boltzmann computing, and (c) allows to perform stable simulations at vanishingly low viscosities. Our findings open attractive prospects for high-performance simulations of realistic turbulent flows on GPU-based architectures. Such expectations are confirmed by excellent agreement among lattice Boltzmann, experimental, and DNS reference data.

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.

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.