Calibration of lubrication force measurements by lattice Boltzmann simulations

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Many experiments explore the hydrodynamic boundary of a surface by approaching a colloidal sphere and measuring the occurring lubrication force. However, in this case many different parameters like wettability and surface roughness influence the result. In the experiment these cannot be separated easily. For a deeper understanding of such surface effects a tool is required that predicts the influence of different surface properties. Here computer simulations can help. In this paper we present lattice Boltzmann simulations of a sphere submerged in a Newtonian liquid and show that our method is able to reproduce the theoretical predictions. In order to provide high precision simulation results the influence of finite size effects has to be controlled. We study the influence of the required system size and resolution of the sphere and demonstrate that already moderate computing ressources allow to keep the error below 1%.This study is funded by DFG priority program SPP 1164

Simulations of slip flow on nanobubble-laden surfaces

On microstructured hydrophobic surfaces, geometrical patterns may lead to the appearance of a superhydrophobic state, where gas bubbles at the surface can have a strong impact on the fluid flow along such surfaces. In particular, they can strongly influence a detected slip at the surface. We present two-phase lattice Boltzmann simulations of a flow over structured surfaces with attached gas bubbles and demonstrate how the detected slip depends on the pattern geometry, the bulk pressure, or the shear rate. Since a large slip leads to reduced friction, our results allow to assist in the optimization of microchannel flows for large throughput.Comment: 22 pages, 12 figure

Lattice Boltzmann simulations of apparent slip in hydrophobic microchannels

Various experiments have found a boundary slip in hydrophobic microchannel flows, but a consistent understanding of the results is still lacking. While Molecular Dynamics (MD) simulations cannot reach the low shear rates and large system sizes of the experiments, it is often impossible to resolve the needed details with macroscopic approaches. We model the interaction between hydrophobic channel walls and a fluid by means of a multi-phase lattice Boltzmann model. Our mesoscopic approach overcomes the limitations of MD simulations and can reach the small flow velocities of known experiments. We reproduce results from experiments at small Knudsen numbers and other simulations, namely an increase of slip with increasing liquid-solid interactions, the slip being independent of the flow velocity, and a decreasing slip with increasing bulk pressure. Within our model we develop a semi-analytic approximation of the dependence of the slip on the pressure.Comment: 7 pages, 4 figure

Optimization of chaotic micromixers using finite time Lyapunov exponents

In microfluidics mixing of different fluids is a highly non-trivial task due to the absence of turbulence. The dominant process allowing mixing at low Reynolds number is therefore diffusion, thus rendering mixing in plain channels very inefficient. Recently, passive chaotic micromixers such as the staggered herringbone mixer were developed, allowing efficient mixing of fluids by repeated stretching and folding of the fluid interfaces. The optimization of the geometrical parameters of such mixer devices is often performed by time consuming and expensive trial and error experiments. We demonstrate that the application of the lattice Boltzmann method to fluid flow in highly complex mixer geometries together with standard techniques from statistical physics and dynamical systems theory can lead to a highly efficient way to optimize micromixer geometries. The strategy applies massively parallel fluid flow simulations inside a mixer, where massless and noninteracting tracer particles are introduced. By following their trajectories we can calculate finite time Lyapunov exponents in order to quantify the degree of chaotic advection inside the mixer. The current report provides a review of our results published in [1] together with additional details on the simulation methodology

Lattice Boltzmann simulations in microfluidics: probing the no-slip boundary condition in hydrophobic, rough, and surface nanobubble laden microchannels

In this contribution we review recent efforts on investigations of the effect of (apparent) boundary slip by utilizing lattice Boltzmann simulations. We demonstrate the applicability of the method to treat fundamental questions in microfluidics by investigating fluid flow in hydrophobic and rough microchannels as well as over surfaces covered by nano- or microscale gas bubbles.Comment: 11 pages, 6 figure

Quantification of the performance of chaotic micromixers on the basis of finite time Lyapunov exponents

Chaotic micromixers such as the staggered herringbone mixer developed by Stroock et al. allow efficient mixing of fluids even at low Reynolds number by repeated stretching and folding of the fluid interfaces. The ability of the fluid to mix well depends on the rate at which "chaotic advection" occurs in the mixer. An optimization of mixer geometries is a non trivial task which is often performed by time consuming and expensive trial and error experiments. In this paper an algorithm is presented that applies the concept of finite-time Lyapunov exponents to obtain a quantitative measure of the chaotic advection of the flow and hence the performance of micromixers. By performing lattice Boltzmann simulations of the flow inside a mixer geometry, introducing massless and non-interacting tracer particles and following their trajectories the finite time Lyapunov exponents can be calculated. The applicability of the method is demonstrated by a comparison of the improved geometrical structure of the staggered herringbone mixer with available literature data.Comment: 9 pages, 8 figure

Boundary Effects in Microfluidic Setups

Due to large surface to volume ratios in microfluidic setups, the roughness of channel surfaces must not be neglected since it is not any longer small compared to the length scale of the system. In addition, the wetting properties of the wall have an important influence on the flow. Even though these effects are getting more and more important for industrial and scientific applications, the knowledge about the interplay of surface roughness and hydrophobic fluid-surface interaction is still very limited because these properties cannot be decoupled easily in experiments. We investigate the problem by means of lattice Boltzmann (LB) simulations of rough microchannels with tunable fluid-wall interaction. We introduce an effective no-slip plane'' at an intermediate position between peaks and valleys of the surface and observe how the position of the wall may change due to surface roughness and hydrophobic interactions. We find that the position of the effective wall, in the case of a Gaussian distributed roughness depends linearly on the width of the distribution. Further we are able to show that roughness creates a non-linear effect on the slip length for hydrophobic boundaries