No Arabic abstract
Hydrodynamic interactions in systems comprised of self-propelled particles, such as swimming microorganisms, and passive tracers have a significant impact on the tracer dynamics compared to the equivalent dry sample. However, such interactions are often difficult to take into account in simulations due to their computational cost. Here, we perform a systematic investigation of swimmer-tracer interaction using an efficient force/counter-force based lattice-Boltzmann (LB) algorithm [J. de Graaf~textit{et al.}, J. Chem. Phys.~textbf{144}, 134106 (2016)] in order to validate its ability to capture the relevant low-Reynolds-number physics. We show that the LB algorithm reproduces far-field theoretical results well, both in a system with periodic boundary conditions and in a spherical cavity with no-slip walls, for which we derive expressions here. The force-lattice coupling of the LB algorithm leads to a smearing out of the flow field, which strongly perturbs the tracer trajectories at close swimmer-tracer separations, and we analyze how this effect can be accurately captured using a simple renormalized hydrodynamic theory. Finally, we show that care must be taken when using LB algorithms to simulate systems of self-propelled particles, since its finite momentum transport time can lead to significant deviations from theoretical predictions based on Stokes flow. These insights should prove relevant to the future study of large-scale microswimmer suspensions using these methods.
We study numerically the effect of thermal fluctuations and of variable fluid-substrate interactions on the spontaneous dewetting of thin liquid films. To this aim, we use a recently developed lattice Boltzmann method for thin liquid film flows, equipped with a properly devised stochastic term. While it is known that thermal fluctuations yield shorter rupture times, we show that this is a general feature of hydrophilic substrates, irrespective of the contact angle. The ratio between deterministic and stochastic rupture times, though, decreases with $theta$. Finally, we discuss the case of fluctuating thin film dewetting on chemically patterned substrates and its dependence on the form of the wettability gradients.
By means of lattice-Boltzmann simulations the drag force on a sphere of radius R approaching a superhydrophobic striped wall has been investigated as a function of arbitrary separation h. Superhydrophobic (perfect-slip vs. no-slip) stripes are characterized by a texture period L and a fraction of the gas area $phi$. For very large values of h/R we recover the macroscopic formulae for a sphere moving towards a hydrophilic no-slip plane. For h/R=O(1) and smaller the drag force is smaller than predicted by classical theories for hydrophilic no-slip surfaces, but larger than expected for a sphere interacting with a uniform perfectly slipping wall. At a thinner gap, $hll R$ the force reduction compared to a classical result becomes more pronounced, and is maximized by increasing $phi$. In the limit of very small separations our simulation data are in quantitative agreement with an asymptotic equation, which relates a correction to a force for superhydrophobic slip to texture parameters. In addition, we examine the flow and pressure field and observe their oscillatory character in the transverse direction in the vicinity of the wall, which reflects the influence of the heterogeneity and anisotropy of the striped texture. Finally, we investigate the lateral force on the sphere, which is detectable in case of very small separations and is maximized by stripes with $phi=0.5$.
The squirmer is a simple yet instructive model for microswimmers, which employs an effective slip velocity on the surface of a spherical swimmer to describe its self-propulsion. We solve the hydrodynamic flow problem with the lattice Boltzmann (LB) method, which is well-suited for time-dependent problems involving complex boundary conditions. Incorporating the squirmer into LB is relatively straight-forward, but requires an unexpectedly fine grid resolution to capture the physical flow fields and behaviors accurately. We demonstrate this using four basic hydrodynamic tests: Two for the far-field flow---accuracy of the hydrodynamic moments and squirmer-squirmer interactions---and two that require the near field to be accurately resolved---a squirmer confined to a tube and one scattering off a spherical obstacle---which LB is capable of doing down to the grid resolution. We find good agreement with (numerical) results obtained using other hydrodynamic solvers in the same geometries and identify a minimum required resolution to achieve this reproduction. We discuss our algorithm in the context of other hydrodynamic solvers and present an outlook on its application to multi-squirmer problems.
In recent years, the lattice Boltzmann (LB) method has been widely employed to simulate boiling phenomena [A. Markus and G. Hazi, Phys. Rev. E 83, 046705 (2011); Biferale et al., Phys. Rev. Lett. 108, 104502 (2012); Li et al., Phys. Rev. E 96, 063303 (2017); Wu et al., Int. J. Heat Mass Transfer 126, 773 (2018)]. However, a very important issue still remains open, i.e., how does boiling occur in the LB simulations? For instance, the existing LB studies showed that the boiling on a hydrophobic surface begins at a lower wall superheat than that on a hydrophilic surface, which qualitatively agrees well with experimental studies, but no one has yet explained how this phenomenon appears in the LB simulations and what happened in the simulations after changing the wettability of the heating surface. In this paper, the LB boiling mechanism is revealed by analyzing boiling on a flat surface with mixed wettability and boiling on a structured surface with homogeneous wettability. Through a theoretical analysis, we demonstrate that, when the same wall superheat is applied, in the LB boiling simulations the fluid density near the heating surface decreases faster on a hydrophobic surface than that on a hydrophilic surface. Accordingly, a lower wall superheat can induce the phase transition from liquid to vapor on a hydrophobic surface than that on a hydrophilic surface. Furthermore, a similar theoretical analysis shows that the fluid density decreases fastest at concave corners in the case of a structured surface with homogeneous wettability, which explains why vapor bubbles are nucleated at concave corners in the LB simulations of boiling on structured surfaces.
The effects of mid-range repulsion in Lattice Boltzmann models on the coalescence/breakup behaviour of single-component, non-ideal fluids are investigated. It is found that mid-range repulsive interactions allow the formation of spray-like, multi-droplet configurations, with droplet size directly related to the strength of the repulsive interaction. The simulations show that just a tiny ten-percent of mid-range repulsive pseudo-energy can boost the surface/volume ratio of the phase- separated fluid by nearly two orders of magnitude. Drawing upon a formal analogy with magnetic Ising systems, a pseudo-potential energy is defined, which is found to behave like a quasi-conserved quantity for most of the time-evolution. This offers a useful quantitative indicator of the stability of the various configurations, thus helping the task of their interpretation and classification. The present approach appears to be a promising tool for the computational modelling of complex flow phenomena, such as atomization, spray formation and micro-emulsions, break-up phenomena and possibly glassy-like systems as well.