We investigate numerically, by a hybrid lattice Boltzmann method, the morphology and the dynamics of an emulsion made of a polar active gel, contractile or extensile, and an isotropic passive fluid. We focus on the case of a highly off-symmetric ratio between the active and passive components. In absence of any activity we observe an hexatic-ordered droplets phase, with some defects in the layout. We study how the morphology of the system is affected by activity both in the contractile and extensile case. In the extensile case a small amount of activity favors the elimination of defects in the array of droplets, while at higher activities, first aster-like rotating droplets appear, and then a disordered pattern occurs. In the contractile case, at sufficiently high values of activity, elongated structures are formed. Energy and enstrophy behavior mark the transitions between the different regimes.

We propose the use of a mathematical model in order to denoise X-ray twodimensional patterns. The model, which makes use of a generalized diffusion equation whose diffusion constant depends on the image gradients, enables to obtain an efficient reduction of pattern noise as witnessed by the computed peak of signal to noise ratio. The corresponding MATLAB code is made available.

Aiming to study the bubble cavitation problem in quiescent and sheared liquids, a third-order isothermal lattice Boltzmann model that describes a two-dimensional (2D) fluid obeying the van der Waals equation of state, is introduced. The evolution equations for the distribution functions in this off-lattice model with 16 velocities are solved using the corner-transport-upwind (CTU) numerical scheme on large square lattices (up to 6144 x 6144 nodes). The numerical viscosity and the regularization of the model are discussed for first- and second-order CTU schemes finding that the latter choice allows to obtain a very accurate phase diagram of a nonideal fluid. In a quiescent liquid, the present model allows us to recover the solution of the 2D Rayleigh-Plesset equation for a growing vapor bubble. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation, and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient D and the capillary number Ca is found at small Ca but with a different factor than in equilibrium liquids. A nonlinear regime is observed for Ca greater than or similar to 0.2.

We numerically study the behavior of self-propelled liquid droplets whose motion is triggered by a Marangoni-like flow. This latter is generated by variations of surfactant concentration which affect the droplet surface tension promoting its motion. In the present paper a model for droplets with a third amphiphilic component is adopted. The dynamics is described by Navier-Stokes and convection-diffusion equations, solved by the lattice Boltzmann method coupled with finite-difference schemes. We focus on two cases. First, the study of self-propulsion of an isolated droplet is carried on and, then, the interaction of two self-propelled droplets is investigated. In both cases, when the surfactant migrates towards the interface, a quadrupolar vortex of the velocity field forms inside the droplet and causes the motion. A weaker dipolar field emerges instead when the surfactant is mainly diluted in the bulk. The dynamics of two interacting droplets is more complex and strongly depends on their reciprocal distance. If, in a head-on collision, droplets are close enough, the velocity field initially attracts them until a motionless steady state is achieved. If the droplets are vertically shifted, the hydrodynamic field leads to an initial reciprocal attraction followed by a scattering along opposite directions. This hydrodynamic interaction acts on a separation of some droplet radii otherwise it becomes negligible and droplets motion is only driven by the Marangoni effect. Finally, if one of the droplets is passive, this latter is generally advected by the fluid flow generated by the active one.

Vesicles are involved in a vast variety of transport processes in living organisms. Additionally, they serve as a model for the dynamics of cell suspensions. Predicting the rheological properties of their suspensions is still an open question, as even the interaction of pairs is yet to be fully understood. Here we analyse the effect of a single vesicle, undergoing tank-treading motion, on its surrounding shear flow by studying the induced disturbance field delta(V) over right arrow, the difference between the velocity field in its presence and absence. The comparison between experiments and numerical simulations reveals an impressive agreement. Tracking ridges in the disturbance field magnitude landscape, we identify the principal directions along which the velocity difference field is analysed in the vesicle vicinity. The disturbance magnitude is found to be significant up to about 4 vesicle radii and can be described by a power law decay with the distance d from the vesicle parallel to delta(V) over right arrow parallel to proportional to d(-3/2). This is consistent with previous experimental results on the separation distance between two interacting vesicles under similar conditions, for which their dynamics is altered. This is an indication of vesicles long-range effect via the disturbance field and calls for the proper incorporation of long-range hydrodynamic interactions when attempting to derive rheological properties of vesicle suspensions. Copyright (C) EPLA, 2016

A deep understanding of the dynamics and rheology of suspensions of vesicles, cells, and capsules is relevant for different applications, ranging from soft glasses to blood flow [1]. I will present the study of suspensions of fluid vesicles by a combination of molecular dynamics and mesoscale hydrodynamics simulations (multi-particle collision dynamics) in two dimensions [2], pointing out the big potential of the numerical method to address problems in soft matter. The flow behavior is studied as a function of the shear rate, the volume fraction of vesicles, and the viscosity ratio between inside and outside fluids. Results are obtained for the interactions of two vesicles, the intrinsic viscosity of the suspension, and the cell-free layer near the walls [3-4]. [1] D. Barthes-Biesel, Annu. Rev. Fluid Mech. 48, 25 (2016) [2] R. Finken, A. Lamura, U. Seifert, and G. Gompper, Eur. Phys. J. E 25, 309 (2008) [3] A. Lamura and G. Gompper, EPL 102, 28004 (2013) [4] A. Lamura and G. Gompper, Procedia IUTAM 16, 3 (2015)

A deep understanding of the dynamics and rheology of suspensions of vesicles, cells, and capsules is relevant for different applications, ranging from soft glasses to blood flow [1]. I will present the study of suspensions of fluid vesicles by a combination of molecular dynamics and mesoscale hydrodynamics simulations (multi-particle collision dynamics) in two dimensions [2], pointing out the big potential of the numerical method to address problems in soft matter. The flow behavior is studied as a function of the shear rate, the volume fraction of vesicles, and the viscosity ratio between inside and outside fluids. Results are obtained for the interactions of two vesicles, the intrinsic viscosity of the suspension, and the cell-free layer near the walls [3-5]. [1] D. Barthes-Biesel, Annu. Rev. Fluid Mech. 48, 25 (2016) [2] R. Finken, A. Lamura, U. Seifert, and G. Gompper, Eur. Phys. J. E 25, 309 (2008) [3] A. Lamura and G. Gompper, EPL 102, 28004 (2013) [4] A. Lamura and G. Gompper, Procedia IUTAM 16, 3 (2015) [5] E. Afik, A. Lamura, and V. Steinberg, EPL 113, 38003 (2016)

Vesicles are involved in a vast variety of transport processes in living organisms. Additionally, they serve as a model for the dynamics of cell suspensions. Predicting the rheological properties of their suspensions is still an open question, as even the interaction of pairs is yet to be fully understood. Here we analyse the effect of a single vesicle, undergoing tank-treading motion, on its surrounding shear flow by studying the induced disturbance field $\delta \vec{V}$, the difference between the velocity field in its presence and absence. The comparison between experiments and numerical simulations reveals an impressive agreement. Tracking ridges in the disturbance field magnitude landscape, we identify the principal directions along which the velocity difference field is analysed in the vesicle vicinity. The disturbance magnitude is found to be significant up to about 4 vesicles radii and can be described by a power law decay with the distance $d$ from the vesicle $ \| \delta \vec{V} \| \propto d^{-3/2}$. This is consistent with previous experimental results on the separation distance between two interacting vesicles under similar conditions, for which their dynamics is altered. This is an indication of vesicles long-range effect via the disturbance field and calls for the proper incorporation of long-range hydrodynamic interactions when attempting to derive rheological properties of vesicle suspensions.

Numerical simulations of vesicle suspensions are performed in two dimensions to study their dynamical and rheological properties. An hybrid method is adopted, which combines a mesoscopic approach for the solvent with a curvature-elasticity model for the membrane. Shear flow is induced by two counter-sliding parallel walls, which generate a linear flow profile. The flow behavior is studied for various vesicle concentrations and viscosity ratios between the internal and the external fluid. Both the intrinsic viscosity and the thickness of depletion layers near the walls are found to increase with increasing viscosity ratio.

The rheology and dynamics of suspensions of fluid vesicles is investigated by a combination of molecular dynamics and mesoscale hydrodynamics simulations in two dimensions. The vesicle suspension is confined between two no-slip shearing walls. The flow behavior is studied as a function of the shear rate, the volume fraction of vesicles, and the viscosity ratio between inside and outside fluids. Results are obtained for the interactions of two vesicles, the intrinsic viscosity of the suspension, and the cell-free layer near the walls.

Cavitation in a liquid moving past a constraint is numerically investigated by means of a free-energy lattice Boltzmann simulation based on the van der Waals equation of state. The fluid is streamed past an obstacle and, depending on the pressure drop between inlet and outlet, vapor formation underneath the corner of the sack-wall is observed. The circumstances of cavitation formation are investigated and it is found that the local bulk pressure and mean stress are insufficient to explain the phenomenon. Results obtained in this study strongly suggest that the viscous stress, interfacial contributions to the local pressure, and the Laplace pressure are relevant to the opening of a vapor cavity. This can be described by a generalization of Joseph's criterion that includes these contributions. A macroscopic investigation, closer in spirit to engineering practice, measuring mass flow rate behavior and discharge coefficient was also performed. Simulations were carried out by fixing the total upstream pressure and varying the static downstream pressure for different kinematic viscosities. As theoretically predicted, mass flow rate increases linearly with the square root of the pressure drop. However, when cavitation occurs the mass flow growth rate is reduced and eventually it collapses into a choked flow state. Reduction of the mass flow growth rate coincides with a smaller discharge coefficient. Therefore, in the cavitating regime, as theoretically predicted and experimentally verified, the discharge coefficient grows with the Nurick cavitation number. On the other hand, in the non-cavitating regime the discharge coefficient grows with the Reynolds number due to the reduction of the boundary layer thickness.

The bubble cavitation problem in quiescent and sheared liquids is investigated using a third-order isothermal lattice Boltzmann (LB) model that describes a two-dimensional ($2D$) fluid obeying the van der Waals equation of state. The LB model has 16 off-lattice velocities and is based on the Gauss-Hermite quadrature method. The evolution equations for the distribution functions in this model are solved using the corner transport upwind numerical scheme on large square lattices (up to $4096 \times 4096$ nodes). In a quiescent liquid, the computer simulation results are in good agreement to the $2D$ Rayleigh-Plesset equation. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient $D$ and the capillary number $Ca$ is found at small $Ca$ but with a different factor than in equilibrium liquids. A non-linear regime is observed for $Ca \gtrsim 0.3$.

Cavitation in a liquid moving past a constraint is numerically investigated by means of a free-energy lattice Boltzmann simulation based on the van der Waals equation of state. The fluid is streamed past an obstacle, and depending on the pressure drop between inlet and outlet, vapor formation underneath the corner of the sack-wall is observed. The circumstances of cavitation formation are investigated and it is found that the local bulk pressure and mean stress are insufficient to explain the phenomenon. Results obtained in this study strongly suggest that the viscous stress, interfacial contributions to the local pressure, and the Laplace pressure are relevant to the opening of a vapor cavity. This can be described by a generalization of Joseph's criterion that includes these contributions. A macroscopic investigation measuring mass flow rate behavior and discharge coefficient was also performed. As theoretically predicted, mass flow rate increases linearly with the square root of the pressure drop. However, when cavitation occurs, the mass flow growth rate is reduced and eventually it collapses into a choked flow state. In the cavitating regime, as theoretically predicted and experimentally verified, the discharge coefficient grows with the Nurick cavitation number. (C) 2015 AIP Publishing LLC.

In this paper the phase behavior and pattern formation in a sheared nonideal fluid under a periodic potential is studied. An isothermal two-dimensional formulation of a lattice Boltzmann scheme for a liquid-vapor system with the van der Waals equation of state is presented and validated. Shear is applied by moving walls and the periodic potential varies along the flow direction. A region of the parameter space, where in the absence of flow a striped phase with oscillating density is stable, will be considered. At low shear rates the periodic patterns are preserved and slightly distorted by the flow. At high shear rates the striped phase loses its stability and traveling waves on the interface between the liquid and vapor regions are observed. These waves spread over the whole system with wavelength only depending on the length of the system. Velocity field patterns, characterized by a single vortex, will also be shown.

We analyze the dynamics of a two-dimensional system of interacting active dumbbells. We characterize the mean-square displacement, linear response function, and deviation from the equilibrium fluctuation-dissipation theorem as a function of activity strength, packing fraction, and temperature for parameters such that the system is in its homogeneous phase. While the diffusion constant in the last diffusive regime naturally increases with activity and decreases with packing fraction, we exhibit an intriguing nonmonotonic dependence on the activity of the ratio between the finite-density and the single-particle diffusion constants. At fixed packing fraction, the time-integrated linear response function depends nonmonotonically on activity strength. The effective temperature extracted from the ratio between the integrated linear response and the mean-square displacement in the last diffusive regime is always higher than the ambient temperature, increases with increasing activity, and, for small active force, monotonically increases with density while for sufficiently high activity it first increases and next decreases with the packing fraction. We ascribe this peculiar effect to the existence of finite-size clusters for sufficiently high activity and density at the fixed (low) temperatures at which we worked. The crossover occurs at lower activity or density the lower the external temperature. The finite-density effective temperature is higher (lower) than the single dumbbell one below (above) a crossover value of the Péclet number.

A systems of self-propelled dumbbells interacting by a Weeks-Chandler-Anderson potential is considered. At sufficiently low temperatures the system phase separates into a dense phase and a gas-like phase. The kinetics of the cluster formation and the growth law for the average cluster size are analyzed.

The dynamics and rheology of suspensions of fluid vesicles or red blood cells is investigated by a combination of molecular dynamics and mesoscale hydrodynamics simulations in two dimensions. The vesicle suspension is confined between two no-slip walls, which are driven externally to generate a shear flow with shear rate (gamma) over dot. The flow behavior is studied as a function of (gamma) over dot, the volume fraction of vesicles, and the viscosity contrast between inside and outside fluids. Results are obtained for the encounter and interactions of two vesicles, the intrinsic viscosity of the suspension, and the cell-free layer near the walls.