Particle stabilized emulsions are ubiquitous in the food and cosmetics industry, but our understanding of the influence of microscopic fluid-particle and particle-particle interactions on the macroscopic rheology is still limited. In this paper we present a simulation algorithm based on a multicomponent lattice Boltzmann model to describe the solvents combined with a molecular dynamics solver for the description of the solved particles. It is shown that the model allows a wide variation of fluid properties and arbitrary contact angles on the particle surfaces. We demonstrate its applicability by studying the transition from a "bicontinuous interfacially jammed emulsion gel" (bijel) to a "Pickering emulsion" in dependence on the contact angle, the particle concentration, and the ratio of the solvents.

Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

Abstract. We investigate properties of dense suspensions and sediments of small spherical silt particles by means of a combined Molecular Dynamics (MD) and Stochastic Rotation Dynamics (SRD) simulation. We include van der Waals and effective electrostatic interactions between the colloidal particles, as well as Brownian motion and hydrodynamic interactions which are calculated in the SRD-part. We present the simulation technique and first results. We have measured velocity distributions, diffusion coefficients, sedimentation velocity, spatial correlation functions and we have explored the phase diagram depending on the parameters of the potentials and on the volume fraction.

On hydrophobic surfaces, roughness may lead to a transition to a superhydrophobic state, where gas bubbles at the surface can have a strong impact on a detected slip. We present two-phase lattice Boltzmann simulations of a Couette flow over structured surfaces with attached gas bubbles. Even though the bubbles add slippery surfaces to the channel, they can cause negative slip to appear due to the increased roughness. The simulation method used allows the bubbles to deform due to viscous stresses. We find a decrease of the detected slip with increasing shear rate which is in contrast to some recent experimental results implicating that bubble deformation cannot account for these experiments. Possible applications of bubble surfaces in microfluidic devices are discussed.

We present a classification scheme for phase transitions in finite systems like atomic and molecular clusters based on the Lee-Yang zeros in the complex temperature plane. In the limit of infinite particle numbers the scheme reduces to the Ehrenfest definition of phase transitions and gives the right critical indices. We apply this classification scheme to Bose-Einstein condensates in a harmonic trap as an example of a higher order phase transitions in a finite system and to small Ar clusters.

Partitioning of red blood cells (RBCs) at the level of bifurcations in the microcirculatory system affects many physiological functions yet it remains poorly understood. We address this problem by using T-shaped microfluidic bifurcations as a model. Our computer simulations and in vitro experiments reveal that the hematocrit (ϕ0) partition depends strongly on RBC deformability, as long as ϕ0<20% (within the normal range in microcirculation), and can even lead to complete deprivation of RBCs in a child branch. Furthermore, we discover a deviation from the Zweifach-Fung effect which states that the child branch with lower flow rate recruits less RBCs than the higher flow rate child branch. At small enough ϕ0, we get the inverse scenario, and the hematocrit in the lower flow rate child branch is even higher than in the parent vessel. We explain this result by an intricate up-stream RBC organization and we highlight the extreme dependence of RBC transport on geometrical and cell mechanical properties. These parameters can lead to unexpected behaviors with consequences on the microcirculatory function and oxygen delivery in healthy and pathological conditions.

We present three-dimensional numerical simulations, employing the well-established lattice Boltzmann method, and investigate similarities and differences between surfactants and nanoparticles as additives at a fluid-fluid interface. We report on their respective effects on the surface tension of such an interface. Next, we subject a fluid droplet to shear and explore the deformation properties of the droplet, its inclination angle relative to the shear flow, the dynamics of the particles at the interface, and the possibility of breakup. Particles are seen not to affect the surface tension of the interface, although they do change the overall interfacial free energy. The particles do not remain homogeneously distributed over the interface, but form clusters in preferred regions that are stable for as long as the shear is applied. However, although the overall structure remains stable, individual nanoparticles roam the droplet interface, with a frequency of revolution that is highest in the middle of the droplet interface, normal to the shear flow, and increases with capillary number. We recover Taylor's law for small deformation of droplets when surfactant or particles are added to the droplet interface. The effect of surfactant is captured in the capillary number, but the inertia of adsorbed massive particles increases deformation at higher capillary number and eventually leads to easier breakup of the droplet.Comment: 17 pages, 17 figures. The figure quality was reduced to fulfill arXiv's file size restriction

Dynamics of a single vesicle under shear flow between two parallel plates is studied in two-dimensions using lattice-Boltzmann simulations. We first present how we adapted the lattice-Boltzmann method to simulate vesicle dynamics, using an approach known from the immersed boundary method. The fluid flow is computed on an Eulerian regular fixed mesh while the location of the vesicle membrane is tracked by a Lagrangian moving mesh. As benchmarking tests, the known vesicle equilibrium shapes in a fluid at rest are found and the dynamical behavior of a vesicle under simple shear flow is being reproduced. Further, we focus on investigating the effect of the confinement on the dynamics, a question that has received little attention so far. In particular, we study how the vesicle steady inclination angle in the tank-treading regime depends on the degree of confinement. The influence of the confinement on the effective viscosity of the composite fluid is also analyzed. At a given reduced volume (the swelling degree) of a vesicle we find that both the inclination angle, and the membrane tank-treading velocity decrease with increasing confinement. At sufficiently large degree of confinement the tank-treading velocity exhibits a nonmonotonous dependence on the reduced volume and the effective viscosity shows a nonlinear behavior.