We present a second-order accurate method to include arbitrary distributions of force densities in the lattice Boltzmann formulation of hydrodynamics. Our method may be used to represent singular force densities arising either from momentum-conserving internal forces or from external forces which do not conserve momentum. We validate our method with several examples involving point forces and find excellent agreement with analytical results. A minimal model for dilute sedimenting particles is presented using the method which promises a substantial gain in computational efficiency.
Liquid capillary-bridge formation between solid particles has a critical influence on the rheological properties of granular materials and, in particular, on the efficiency of fluidized bed reactors. The available analytical and semi-analytical metho
ds have inherent limitations, and often do not cover important aspects, like the presence of non-axisymmetric bridges. Here, we conduct numerical simulations of the capillary bridge formation between equally and unequally-sized solid particles using the lattice Boltzmann method, and provide an assessment of the accuracy of different families of analytical models. We find that some of the models taken into account are shown to perform better than others. However, all of them fail to predict the capillary force for contact angles larger than $pi/2$, where a repulsive capillary force attempts to push the solid particle outwards to minimize the surface energy, especially at a small separation distance.
Capillary bridges can form between colloids immersed in a two phase fluid, e.g., in a binary liquid mixture, if the surface of the colloids prefers the species other than the one favored in the bulk liquid. Here, we study the formation of liquid brid
ges induced by confining colloids to a slit, with the slit walls having a preference opposite to the one of the colloid surface. Using mean field theory, we show that there is a line of first-order phase transitions between the bridge and the no-bridge states, which ends at a critical point. By decreasing the slit width, this critical point is shifted towards smaller separations between the colloids. However, at very small separations, and far from criticality, we observe only a minor influence of the slit width on the location of the transition. Monte Carlo simulations of the Ising model, which mimics incompressible binary liquid mixtures, confirm the occurrence of the bridging transitions, as manifested by the appearance of bistable regions where both the bridge and the no-bridge configurations are (meta)stable. Interestingly, we find no bistability in the case of small colloids, but we observe a sharpening of the transition when the colloid size increases. In addition, we demonstrate that the capillary force acting between the colloids can depend sensitively on the slit width, and varies drastically with temperature, thus achieving strengths orders of magnitude higher than at criticality of the fluid.
We study the normal and lateral effective critical Casimir forces acting on a spherical colloid immersed in a critical binary solvent and close to a chemically structured substrate with alternating adsorption preference. We calculate the universal sc
aling function for the corresponding potential and compare our results with recent experimental data [Soyka F., Zvyagolskaya O., Hertlein C., Helden L., and Bechinger C., Phys. Rev. Lett., 101, 208301 (2008)]. The experimental potentials are properly captured by our predictions only by accounting for geometrical details of the substrate pattern for which, according to our theory, critical Casimir forces turn out to be a sensitive probe.
Run-and-tumble (RNT) motion is a prominent locomotion strategy employed by many living microorganisms. It is characterized by straight swimming intervals (runs), which are interrupted by sudden reorientation events (tumbles). In contrast, directional
changes of synthetic microswimmers (active particles, APs) are caused by rotational diffusion, which is superimposed with their translational motion and thus leads to rather continuous and slow particle reorientations. Here we demonstrate that active particles can also perform a swimming motion where translational and orientational changes are disentangled, similar to RNT. In our system, such motion is realized by a viscoelastic solvent and a periodic modulation of the self-propulsion velocity. Experimentally, this is achieved using light-activated Janus colloids, which are illuminated by a time-dependent laser field. We observe a strong enhancement of the effective translational and rotational motion when the modulation time is comparable to the relaxation time of the viscoelastic fluid. Our findings are explained by the relaxation of the elastic stress, which builds up during the self-propulsion, and is suddenly released when the activity is turned off. In addition to a better understanding of active motion in viscoelastic surroundings, our results may suggest novel steering strategies for synthetic microswimmers in complex environments.
When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. We performed Monte Carlo simulations to study the phase diagram of such rod-polymer mixture. The
colloidal rods were modelled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while overlap of polymers and rods is forbidden. Large amounts of depletant cause phase separation of the mixture. We estimated the phase boundaries of isotropic-isotropic coexistence both, in the bulk and in confinement. To determine the phase boundaries we applied the grand canonical ensemble using successive umbrella sampling [J. Chem. Phys. 120, 10925 (2004)], and we performed a finite-size scaling analysis to estimate the location of the critical point. The results are compared with predictions of the free volume theory developed by Lekkerkerker and Stroobants [Nuovo Cimento D 16, 949 (1994)]. We also give estimates for the interfacial tension between the coexisting isotropic phases and analyse its power-law behaviour on approach of the critical point.