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 numerically investigate the rheological response of a noncoalescing multiple emulsion under a symmetric shear flow. We find that the dynamics significantly depends on the magnitude of the shear rate and on the number of the encapsulated droplets, two key parameters whose control is fundamental to accurately select the resulting nonequiibrium steady states. The double emulsion, for instance, attains a static steady state in which the external droplet stretches under flow and achieves an elliptical shape (closely resembling the one observed in a sheared isolated fluid droplet), while the internal one remains essentially unaffected. Novel nonequiibrium steady states arise in a multiple emulsion. Under low/moderate shear rates, for instance, the encapsulated droplets display a nontrivial planetarylike motion that considerably affects the shape of the external droplet. Some features of this dynamic behavior are partially captured by the Taylor deformation parameter and the stress tensor. Besides a theoretical interest on its own, our results can potentially stimulate further experiments, as most of the predictions could be tested in the lab by monitoring droplets' shapes and position over time. Published under license by AIP Publishing.

The study of the underlying physics of soft flowing materials depends heavily on numerical simulations, due to the complex structure of the governing equations reflecting the competition of concurrent mechanisms acting at widely disparate scales in space and time. A full-scale computational modelling remains a formidable challenge since it amounts to simultaneously handling six or more spatial decades in space and twice as many in time. Coarse-grained methods often provide a viable strategy to significantly mitigate this issue, through the implementation of mesoscale supramolecular forces designed to capture the essential physics at a fraction of the computational cost of a full-detail description. Here, we review some recent advances in the design of a lattice Boltzmann mesoscale approach for soft flowing materials, inclusive of near-contact interactions (NCIs) between dynamic interfaces, as they occur in high packing-fraction soft flowing crystals. The method proves capable of capturing several aspects of the rheology of soft flowing crystals, namely, (i) a 3/2 power-law dependence of the dispersed phase flow rate on the applied pressure gradient, (ii) the structural transition between an ex-two and ex-one (bamboo) configurations with the associated drop of the flow rate, (iii) the onset of interfacial waves once NCI is sufficiently intense. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

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.

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.

Bicontinuous interfacially jammed emulsion gels ("bijels") represent a new class of soft materials made of a densely packed monolayer of solid particles sequestered at the interface of a bicontinuous fluid. Their mechanical properties are relevant to many applications, such as catalysis, energy conversion, soft robotics, and scaffolds for tissue engineering. While their stationary bulk properties have been covered in depth, much less is known about their behavior in the presence of an external shear. In this paper, we numerically study the dynamics of a bijel confined within a three-dimensional rectangular domain and subject to a symmetric shear flow sufficiently intense to break the material. Extensive numerical simulations reveal that the shear flow generally promotes the detachment of a sizable amount of particles from the fluid interface and their accumulation in the bulk. Fluid interfaces undergo large stretching and deformations along the flow direction, an effect that reduces their capability of entrapping particles. These results are supported by a series of quantitative indicators such as (i) curvature of the fluid interface, (ii) spatial distribution of the colloidal particles, and (iii) fluid flow structure within the microchannel. (c) 2020 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

We provide a brief survey of our current developments in the simulation-based design of novel families of mesoscale porous materials using computational kinetic theory. Prospective applications on exascale computers are also briefly discussed and commented on, with reference to two specific examples of soft mesoscale materials: microfluid crystals and bi-continuous jels.

We numerically study the dynamic behavior under a symmetric shear flow of selected examples of concentrated phase emulsions with multicore morphology confined within a microfluidic channel. A variety of nonequilibrium steady states is reported. Under low shear rates, the emulsion is found to exhibit a solidlike behavior, in which cores display a periodic planetarylike motion with approximately equal angular velocity. At higher shear rates, two steady states emerge, one in which all inner cores align along the flow and become essentially motionless and a further one in which some cores accumulate near the outer interface and produce a dynamical elliptical-shaped ring chain, reminiscent of a treadmillinglike structure, while others occupy the center of the emulsion. A quantitative description in terms of the (i) motion of the cores, (ii) rate of deformation of the emulsion, and (iii) structure of the fluid flow within the channel is also provided.

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.

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 largescale LB simulations, as well as possible strategies towards fault-tolerant LB schemes. (C) 2019 Published by Elsevier Ltd.

Fully three-dimensional, time-dependent, direct simulations of the non-ideal Navier-Stokes equations for a two-component fluid shed light into the mechanism which inhibits droplet breakup in step emulsifiers below a critical threshold of the width-to-height (w/h) ratio of the microfluidic nozzle. Below w/h similar to 2.6, the simulations provide evidence of a smooth topological transition of the fluid from the confined rectangular channel geometry to an isotropic (spherical) expansion of the fluid downstream the nozzle step. Above such threshold, the transition from the inner to the outer space involves a series of dynamical rearrangements which keep the free surface in mechanical balance. Such rearrangements also induce a backflow of the ambient fluid which, in turn, leads to jet pinching and ultimately to its rupture, namely, droplet formation. The simulations show remarkable agreement with the experimental value of the threshold, which is found around w/h similar to 2.56. Published under license by AIP Publishing.

This article is part of the theme issue 'Multiscale modelling, simulation and computing: from the desktop to the exascale'. We outline the main ideas behind the numerical modelling of soft flowing crystals, paying special attention to their application to microfluidic devices for the design of novel mesoscale porous materials.

Computer simulations of bi-continuous two-phase fluids with interspersed 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 (observed even in the absence of colloids) 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.

We present a mesoscale representation of near-contact interactions between colliding droplets which permits one 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-focusing microfluidic devices. This model may open up the possibility of multiscale simulation of microfluidic devices for the production of new droplet and bubble-based mesoscale porous materials.

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

The dynamic interaction of complex fluid interfaces is highly sensitive to near-contact interactions occurring at the scale of ten of nanometers. Such interactions are difficult to analyze because they couple self-consistently to the dynamic morphology of the evolving interface, as well as to the hydrodynamics of the interstitial fluid film. In this work, we show that, above a given magnitude threshold, near-contact interactions trigger nontrivial microvorticity patterns, which in turn affect the effective near-contact interactions, giving rise to persistent fluctuating ripples at the fluid interface. In such a regime, near-contact interactions may significantly affect the macroscopic arrangement of emulsion configurations, such as those arising in soft-flowing microfluidic crystals.

We propose a mesoscopic model of binary fluid mixtures with tunable viscosity ratio based on a two-range pseudopotential lattice Boltzmann method, for the simulation of soft flowing systems. In addition to the short-range repulsive interaction between species in the classical single-range model, a competing mechanism between the short-range attractive and midrange repulsive interactions is imposed within each species. Besides extending the range of attainable surface tension as compared with the single-range model, the proposed scheme is also shown to achieve a positive disjoining pressure, independently of the viscosity ratio. The latter property is crucial for many microfluidic applications involving a collection of disperse droplets with a different viscosity from that of the continuum phase. As a preliminary application, the relative effective viscosity of a pressure-driven emulsion in a planar channel is computed.

Three-dimensional, time-dependent direct simulations of step emulsification microdevices highlight two essential mechanisms for droplet formation: first, the onset of an adverse pressure gradient driving a backflow of the continuous phase from the external reservoir to the microchannel, and second, the stnction 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 stabilizing effect of dynamic pressure. This suggests a new criterion for dripping-jetting transition, based on local values of the capillary and Weber numbers.