Interface tracking characteristics of color-gradient lattice Boltzmann model for immiscible fluids.

We study the interface tracking characteristics of a color-gradient-based lattice Boltzmann model for immiscible flows. Investigation of the local density change in one of the fluid phases, via a Taylor series expansion of the recursive lattice Boltzmann equation, leads to the evolution equation of the order parameter that differentiates the fluids. It turns out that this interface evolution follows a conservative Allen-Cahn equation with a mobility which is independent of the fluid viscosities and surface tension. The mobility of the interface, which solely depends upon lattice speed of sound, can have a crucial effect on the physical dynamics of the interface. Further, we find that, when the equivalent lattice weights inside the segregation operator are modified, the resulting differential operators have a discretization error that is anisotropic to the leading order. As a consequence, the discretization errors in the segregation operator, which ensures a finite interface width, can act as a source of the spurious currents. These findings are supported with the help of numerical simulations.

Modeling the flow in diffuse interface methods of solidification.

Fluid dynamical equations in the presence of a diffuse solid-liquid interface are investigated via a volume averaging approach. The resulting equations exhibit the same structure as the standard Navier-Stokes equation for a Newtonian fluid with a constant viscosity, the effect of the solid phase fraction appearing in the drag force only. This considerably simplifies the use of the lattice Boltzmann method as a fluid dynamics solver in solidification simulations. Galilean invariance is also satisfied within this approach. Further, we investigate deviations between the diffuse and sharp interface flow profiles via both quasiexact numerical integration and lattice Boltzmann simulations. It emerges from these studies that the freedom in choosing the solid-liquid coupling parameter h provides a flexible way of optimizing the diffuse interface-flow simulations. Once h is adapted for a given spatial resolution, the simulated flow profiles reach an accuracy comparable to quasiexact numerical simulations.

Shear-flow-controlled mode selection in a nonlinear autocatalytic medium.

The effect of shear flow on mode selection and the length scale of patterns formed in a nonlinear autocatalytic reaction-diffusion model is investigated. We predict analytically the existence of transverse and longitudinal modes. The type of the selected mode strongly depends on the difference in the flow rates of the participating species, quantified by the differential flow parameter. Spatial structures are obtained by varying the length scale of individual modes and superposing them via the differential flow parameter. Our predictions are in line with numerical results obtained from lattice Boltzmann simulations.

Fluctuating multicomponent lattice Boltzmann model.

Current implementations of fluctuating lattice Boltzmann equations (FLBEs) describe single component fluids. In this paper, a model based on the continuum kinetic Boltzmann equation for describing multicomponent fluids is extended to incorporate the effects of thermal fluctuations. The thus obtained fluctuating Boltzmann equation is first linearized to apply the theory of linear fluctuations, and expressions for the noise covariances are determined by invoking the fluctuation-dissipation theorem directly at the kinetic level. Crucial for our analysis is the projection of the Boltzmann equation onto the orthonormal Hermite basis. By integrating in space and time the fluctuating Boltzmann equation with a discrete number of velocities, the FLBE is obtained for both ideal and nonideal multicomponent fluids. Numerical simulations are specialized to the case where mean-field interactions are introduced on the lattice, indicating a proper thermalization of the system.

Simulation of static critical phenomena in nonideal fluids with the lattice Boltzmann method.

A fluctuating nonideal fluid at its critical point is simulated with the lattice Boltzmann method. It is demonstrated that the method, employing a Ginzburg-Landau free energy functional, correctly reproduces the static critical behavior associated with the Ising universality class. A finite-size scaling analysis is applied to determine the critical exponents related to the order parameter, compressibility and specific heat. A particular focus is put on finite-size effects and issues related to the global conservation of the order parameter.

Modelling thermal fluctuations in non-ideal fluids with the lattice Boltzmann method.

Recently, we proposed a theoretical framework to include thermal fluctuations into the Lattice Boltzmann (LB) method for non-ideal fluids. Here, we apply a variant thereof to a certain class of force-based non-ideal fluid LB models. We find that ideal-gas-like noise is an exact result of the fluctuation-dissipation theorem in the hydrodynamic regime. It is shown that satisfactory equilibration of the density and fluid momentum can be obtained in a simulation over a wide range of length scales.

Particle stress in suspensions of soft objects.

In a suspension of extended objects such as colloidal particles, capsules or vesicles, the contribution of particles to the stress is usually evaluated by first determining the stress originating from a single particle (e.g. via integrating the fluid stress over the surface of a particle) and then adding up the contributions of individual particles. While adequate for a computation of the average stress over the entire system, this approach fails to correctly reproduce the local stress. In this work, we propose and validate a variant of the method of planes which overcomes this problem. The method is particularly suited for many-body interactions arising from, for example, shear and bending rigidity of red blood cells.

Microfluidic emulsion separation-simultaneous separation and sensing by multilayer nanofilm structures.

Emulsion separation is of high relevance for filtration applications, liquid-liquid-partitioning of biomolecules like proteins and recovery of products from droplet microreactors. Selective interaction of various components of an emulsion with substrates is used to design microfluidic flow chambers for efficient separation of emulsions into their individual components. Our lab-on-a-chip device consists of an emulsion separation cell with an integrated silicon sensor chip, the latter allowing the detection of liquid motion via the field-effect signal. Thus, within our lab-on-a-chip device, emulsions can be separated while the separation process is monitored simultaneously. For emulsion separation a surface energy step gradient, namely a sharp interface between the hydrophobic and hydrophilic parts of the separation chamber, is used. The key component of the lab-on-a-chip system is a multilayer and multifunctional nanofilm structure which not only provides the surface energy step gradient for emulsion separation but also constitutes the functional parts of the field-effect transistors. The proof-of-principle was performed using a model emulsion consisting of immiscible aqueous and organic solvent components. Droplet coalescence was identified as a key aspect influencing the separation process, with quite different effects during separation on open surfaces as compared to slit geometry. For a detailed description of this observation, an analytical model was derived and lattice Boltzmann computer simulations were performed. By use of grazing incidence small angle x-ray scattering (GISAXS) interfacial nanostructures during gold nanoparticle deposition in a flow field were probed to demonstrate the potential of GISAXS for in situ investigations during flow.

Structure and flow of droplets on solid surfaces.

The structure and flow of droplets on solid surfaces is investigated with imaging and scattering techniques and compared to simulations. To access nanostructures at the liquid-solid interface advanced scattering techniques such as grazing incidence small-angle x-ray scattering (GISAXS) with micro- and nanometer-sized beams, GISAXS and in situ imaging ellipsometry and GISAXS tomography are used. Using gold nanoparticle suspensions, structures observed in the wetting area due to deposition are probed in situ during the drying of the droplets. After drying, nanostructures in the wetting area and inside the dried droplets are monitored. In addition to drying, a macroscopic movement of droplets is caused by body forces acting on an inclined substrate. The complexity of the solid surfaces is increased from simple silicon substrates to binary polymer brushes, which undergo a switching due to the liquid in the droplet. Nanostructures introduced in the polymer brush due to the movement of droplets are observed.

Stability and dynamics of droplets on patterned substrates: insights from experiments and lattice Boltzmann simulations.

The stability and dynamics of droplets on solid substrates are studied both theoretically and via experiments. Focusing on our recent achievements within the DFG-priority program 1164 (Nano- and Microfluidics), we first consider the case of (large) droplets on the so-called gradient substrates. Here the term gradient refers to both a change of wettability (chemical gradient) or topography (roughness gradient). While the motion of a droplet on a perfectly flat substrate upon the action of a chemical gradient appears to be a natural consequence of the considered situation, we show that the behavior of a droplet on a gradient of topography is less obvious. Nevertheless, if care is taken in the choice of the topographic patterns (in order to reduce hysteresis effects), a motion may be observed. Interestingly, in this case, simple scaling arguments adequately account for the dependence of the droplet velocity on the roughness gradient (Moradi et al 2010 Europhys. Lett. 89 26006). Another issue addressed in this paper is the behavior of droplets on hydrophobic substrates with a periodic arrangement of square shaped pillars. Here, it is possible to propose an analytically solvable model for the case where the droplet size becomes comparable to the roughness scale (Gross et al 2009 Europhys. Lett. 88 26002). Two important predictions of the model are highlighted here. (i) There exists a state with a finite penetration depth, distinct from the full wetting (Wenzel) and suspended (Cassie-Baxter, CB) states. (ii) Upon quasi-static evaporation, a droplet initially on the top of the pillars (CB state) undergoes a transition to this new state with a finite penetration depth but then (upon further evaporation) climbs up the pillars and goes back to the CB state again. These predictions are confirmed via independent numerical simulations. Moreover, we also address the fundamental issue of the internal droplet dynamics and the terminal center of mass velocity on a flat substrate.

Lattice Boltzmann study of pattern formation in reaction-diffusion systems.

Pattern formation in reaction-diffusion systems is of great importance in surface micropatterning [Grzybowski et al., Soft Matter 1, 114 (2005)], self-organization of cellular micro-organisms [Schulz et al., Annu. Rev. Microbiol. 55, 105 (2001)], and in developmental biology [Barkai et al., FEBS Journal 276, 1196 (2009)]. In this work, we apply the lattice Boltzmann method to study pattern formation in reaction-diffusion systems. As a first methodological step, we consider the case of a single species undergoing transformation reaction and diffusion. In this case, we perform a third-order Chapman-Enskog multiscale expansion and study the dependence of the lattice Boltzmann truncation error on the diffusion coefficient and the reaction rate. These findings are in good agreement with numerical simulations. Furthermore, taking the Gray-Scott model as a prominent example, we provide evidence for the maturity of the lattice Boltzmann method in studying pattern formation in nonlinear reaction-diffusion systems. For this purpose, we perform linear stability analysis of the Gray-Scott model and determine the relevant parameter range for pattern formation. Lattice Boltzmann simulations allow us not only to test the validity of the linear stability phase diagram including Turing and Hopf instabilities, but also permit going beyond the linear stability regime, where large perturbations give rise to interesting dynamical behavior such as the so-called self-replicating spots. We also show that the length scale of the patterns may be tuned by rescaling all relevant diffusion coefficients in the system with the same factor while leaving all the reaction constants unchanged.

Second-order convergence of the deviatoric stress tensor in the standard Bhatnagar-Gross-Krook lattice Boltzmann method.

It is shown numerically that the deviatoric stress tensor is second-order accurate in the bulk Bhatnagar-Gross-Krook lattice Boltzmann (LB) method. In an earlier work [T. Krüger, Phys. Rev. E 79, 46704 (2009)], we have already predicted the second-order convergence. However, numerical simulations using a duct flow were not fully in line with this prediction. In particular, the convergence rate of the stress tensor was observed to depend on the LB boundary condition. In the present paper, we examine a pure bulk system, the decaying Taylor-Green vortex flow. Our prediction on the second-order accuracy of the stress tensor is unambiguously evidenced via these studies.

Acoustic driven flow and lattice Boltzmann simulations to study cell adhesion in biofunctionalized mu-fluidic channels with complex geometry.

Accurately mimicking the complexity of microvascular systems calls for a technology which can accommodate particularly small sample volumes while retaining a large degree of freedom in channel geometry and keeping the price considerably low to allow for high throughput experiments. Here, we demonstrate that the use of surface acoustic wave driven microfluidics systems successfully allows the study of the interrelation between melanoma cell adhesion, the matrix protein collagen type I, the blood clotting factor von Willebrand factor (vWF), and microfluidic channel geometry. The versatility of the tool presented enables us to examine cell adhesion under flow in straight and bifurcated microfluidic channels in the presence of different protein coatings. We show that the addition of vWF tremendously increases (up to tenfold) the adhesion of melanoma cells even under fairly low shear flow conditions. This effect is altered in the presence of bifurcated channels demonstrating the importance of an elaborate hydrodynamic analysis to differentiate between physical and biological effects. Therefore, computer simulations have been performed along with the experiments to reveal the entire flow profile in the channel. We conclude that a combination of theory and experiment will lead to a consistent explanation of cell adhesion, and will optimize the potential of microfluidic experiments to further unravel the relation between blood clotting factors, cell adhesion molecules, cancer cell spreading, and the hydrodynamic conditions in our microcirculatory system.

Thermal fluctuations in the lattice Boltzmann method for nonideal fluids.

We introduce thermal fluctuations in the lattice Boltzmann method for nonideal fluids. A fluctuation-dissipation theorem is derived within the Langevin framework and applied to a specific lattice Boltzmann model that approximates the linearized fluctuating Navier-Stokes equations for fluids based on square-gradient free-energy functionals. The obtained thermal noise is shown to ensure equilibration of all degrees of freedom in a simulation to high accuracy. Furthermore, we demonstrate that satisfactory results for most practical applications of fluctuating hydrodynamics can already be achieved using thermal noise derived in the long-wavelength limit.

Effect of aspect ratio on transverse diffusive broadening: a lattice Boltzmann study.

We study scaling laws characterizing the interdiffusive zone between two miscible fluids flowing side by side in a Y-shape laminar micromixer using the lattice Boltzmann method. The lattice Boltzmann method solves the coupled three-dimensional (3D) hydrodynamics and mass transfer equations and incorporates intrinsic features of 3D flows related to this problem. We observe the different power-law regimes occurring at the center of the channel and close to the top/bottom wall. The extent of the interdiffusive zone scales as the square root of the axial distance at the center of the channel. At the top/bottom wall, we find an exponent 1/3 at early stages of mixing as observed in the experiments of Ismagilov [Appl. Phys. Lett. 76, 2376 (2000)]. At a larger distance from the entrance, the scaling exponent close to the walls changes to 1/2 [J.-B. Salmon and A. Adjari, J. Appl. Phys. 101, 074902 (2007)]. Here, we focus on the effect of the finite aspect ratio on diffusive broadening. Interestingly, we find the same scaling laws regardless of the channel's aspect ratio. However, the point at which the exponent 1/3 characterizing the broadening at the top/bottom wall reverts to the normal diffusive behavior downstream strongly depends on the aspect ratio. We propose an interpretation of this observation in terms of the shear rate at the side walls. A criterion for the range of aspect ratios with non-negligible effect on diffusive broadening is also provided.

Profile blunting and flow blockage in a yield-stress fluid: a molecular dynamics study.

The flow of a simple glass forming system (a 80:20 binary Lennard-Jones mixture) through a planar channel is studied via molecular dynamics simulations. The flow is driven by an external body force similar to gravity. Previous studies show that the model exhibits both a static [F. Varnik, J. Chem. Phys. 120, 2788 (2004)] and a dynamic [F. Varnik and O. Henrich, Phys. Rev. B 73, 174209 (2006)] yield stress in the glassy phase. These observations are corroborated by the present work, where we investigate how the presence of a yield stress may affect the system behavior in a Poiseuille-type flow geometry. In particular, we observe a blunted velocity profile across the channel: A relatively wide region in the channel center flows with a constant velocity (zero shear rate) followed by a nonlinear change of the shear rate as the walls are approached. The observed velocity gradients are compared to those obtained from the knowledge of the shear stress across the channel and the flow curves (stress versus shear rate), the latter being determined in our previous simulations of homogeneous shear flow. Furthermore, using the value of the (dynamic) yield stress known from previous simulations, we estimate the threshold body force for a complete arrest of the flow. Indeed, a blockage is observed as the imposed force falls below this threshold value. Small but finite shear rates are observed at stresses above the dynamic but below the static yield stress. We discuss the possible role of the stick-slip-like motion for this observation.

Structural relaxation and rheological response of a driven amorphous system.

The interplay between the structural relaxation and the rheological response of a simple amorphous system {a 80:20 binary Lennard-Jones mixture [W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 (1994)]} is studied via molecular dynamics simulations. In the quiescent state, the model is well known for its sluggish dynamics and a two step relaxation of correlation functions at low temperatures. An ideal glass transition temperature of Tc=0.435 has been identified in the previous studies via the analysis of the system's dynamics in the framework of the mode coupling theory of the glass transition [W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 (1995)]. In the present work, we focus on the question whether a signature of this ideal glass transition can also be found in the case where the system's dynamics is driven by a shear motion. Indeed, the following distinction in the structural relaxation is found: In the supercooled state, the structural relaxation is dominated by the shear at relatively high shear rates gamma, whereas at sufficiently low gamma the (shear-independent) equilibrium relaxation is recovered. In contrast to this, the structural relaxation of a glass is always driven by shear. This distinct behavior of the correlation functions is also reflected in the rheological response. In the supercooled state, the shear viscosity eta decreases with increasing shear rate (shear thinning) at high shear rates, but then converges toward a constant as the gamma is decreased below a (temperature-dependent) threshold value. Below Tc, on the other hand, the shear viscosity grows as eta proportional, etax 1/gamma, suggesting a divergence at gamma=0. Thus, within the accessible observation time window, a transition toward a nonergodic state seems to occur in the driven glass as the driving force approaches zero. As to the flow curves (stress versus shear rate), a plateau forms at low shear rates in the glassy phase. A consequence of this stress plateau for Poiseuille-type flows is demonstrated.

A study of the static yield stress in a binary Lennard-Jones glass.

The stress-strain relations and the yield behavior of a model glass (a 80:20 binary Lennard-Jones mixture) is studied by means of molecular dynamics simulations. In a previous paper it was shown that, at temperatures below the glass transition temperature, Tg, the model exhibits shear banding under imposed shear. It was also suggested that this behavior is closely related to the existence of a (static) yield stress (under applied stress, the system does not flow until the stress sigma exceeds a threshold value sigmay). A thorough analysis of the static yield stress is presented via simulations under imposed stress. Furthermore, using steady shear simulations, the effect of physical aging, shear rate and temperature on the stress-strain relation is investigated. In particular, we find that the stress at the yield point (the "peak"-value of the stress-strain curve) exhibits a logarithmic dependence both on the imposed shear rate and on the "age" of the system in qualitative agreement with experiments on amorphous polymers, and on metallic glasses. In addition to the very observation of the yield stress which is an important feature seen in experiments on complex systems like pastes, dense colloidal suspensions and foams, further links between our model and soft glassy materials are found. An example is the existence of hysteresis loops in the system response to a varying imposed stress. Finally, we measure the static yield stress for our model and study its dependence on temperature. We find that for temperatures far below the mode coupling critical temperature of the model (Tc = 0.435 in Lennard-Jones units), sigmay decreases slowly upon heating followed by a stronger decrease as Tc is approached. We discuss the reliability of results on the static yield stress and give a criterion for its validity in terms of the time scales relevant to the problem.

Shear localization in a model glass.

Using molecular dynamics simulations, we show that a simple model of a glassy material exhibits the shear localization phenomenon observed in many complex fluids. At low shear rates, the system separates into a fluidized shear band and an unsheared part. The two bands are characterized by a very different dynamics probed by a local intermediate scattering function. Furthermore, a stick-slip motion is observed at very small shear rates. Our results, which open the possibility of exploring complex rheological behavior using simulations, are compared to recent experiments on various soft glasses.

Confinement effects on the slow dynamics of a supercooled polymer melt: Rouse modes and the incoherent scattering function.

Results of large-scale molecular-dynamics simulations of a supercooled polymer film are presented (F. Varnik, J. Baschnagel, K. Binder, J. Phys. IV 10, 239 (2000)). The dynamic and static properties of the system are studied for a wide range of film thicknesses (from 3 to about 55 times the bulk radius of gyration) and temperatures (from the normal liquid state to the supercooled region). The system is confined between two completely smooth and purely repulsive walls. Motivated by the previous results on the enhancement of the local relaxation dynamics due to the confinement (F. Varnik, J. Baschnagel, K. Binder, Eur. Phys. J. E 8, 175 (2002); Phys. Rev. E. 65, 021507 (2002)), we now study the effect of the walls on the dynamics of the Rouse modes. It has been reported from Monte Carlo studies of the Bond Fluctuation Model (BFM) that, contrary to the enhancement of the "cage dynamics" (exemplified by a faster relaxation of the incoherent scattering function in the film), Rouse modes exhibit a slower relaxation in the confined system (C. Mischler, J. Baschnagel, K. Binder, Adv. Colloid Interface Sci. 94, 197 (2001)). However, we do not observe such a discrepancy for our continuum model: At a given temperature, the relaxation of a given Rouse mode is faster in the film than in the bulk in accordance with the acceleration of the dynamics around the cage.