Flow past superhydrophobic surfaces with cosine variation in local slip length

Anisotropic super-hydrophobic surfaces have the potential to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent work has focused mostly on cases of super-hydrophobic stripes. Here, we analyze a relevant situation of cosine variation of the local slip length. We derive approximate formulae for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that are in good agreement with the exact numerical solution and lattice-Bolzmann simulations for any surface slip fraction. The cosine texture can provide a very large effective (forward) slip, but it was found to be less efficient in generating a transverse flow as compared to super-hydrophobic stripes.Comment: 8 pages, 6 figure

Random-roughness hydrodynamic boundary conditions

We report results of lattice Boltzmann simulations of a high-speed drainage of liquid films squeezed between a smooth sphere and a randomly rough plane. A significant decrease in the hydrodynamic resistance force as compared with that predicted for two smooth surfaces is observed. However, this force reduction does not represent slippage. The computed force is exactly the same as that between equivalent smooth surfaces obeying no-slip boundary conditions, but located at an intermediate position between peaks and valleys of asperities. The shift in hydrodynamic thickness is shown to depend on the height and density of roughness elements. Our results do not support some previous experimental conclusions on very large and shear-dependent boundary slip for similar systems.Comment: 4 pages, 4 figure

Implementation of on-site velocity boundary conditions for D3Q19 lattice Boltzmann

On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.Comment: 13 pages, 4 figures, revised versio

Lattice-Boltzmann simulations of the drag force on a sphere approaching a superhydrophobic striped plane

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, $h\ll 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$.Comment: 9 pages, 7 figure

Quantitative analysis of numerical estimates for the permeability of porous media from lattice-Boltzmann simulations

During the last decade, lattice-Boltzmann (LB) simulations have been improved to become an efficient tool for determining the permeability of porous media samples. However, well known improvements of the original algorithm are often not implemented. These include for example multirelaxation time schemes or improved boundary conditions, as well as different possibilities to impose a pressure gradient. This paper shows that a significant difference of the calculated permeabilities can be found unless one uses a carefully selected setup. We present a detailed discussion of possible simulation setups and quantitative studies of the influence of simulation parameters. We illustrate our results by applying the algorithm to a Fontainebleau sandstone and by comparing our benchmark studies to other numerical permeability measurements in the literature.Comment: 14 pages, 11 figure

Simulation of Claylike Colloids

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.Comment: 20 pages, 14 figure

Steering in computational science: mesoscale modelling and simulation

This paper outlines the benefits of computational steering for high performance computing applications. Lattice-Boltzmann mesoscale fluid simulations of binary and ternary amphiphilic fluids in two and three dimensions are used to illustrate the substantial improvements which computational steering offers in terms of resource efficiency and time to discover new physics. We discuss details of our current steering implementations and describe their future outlook with the advent of computational grids.Comment: 40 pages, 11 figures. Accepted for publication in Contemporary Physic

The interplay between shell effects and electron correlations in quantum dots

We use the Path Integral Monte Carlo method to investigate the interplay between shell effects and electron correlations in single quantum dots with up to 12 electrons. By use of an energy estimator based on the hypervirial theorem of Hirschfelder we study the energy contributions of different interaction terms in detail. We discuss under which conditions the total spin of the electrons is given by Hund's rule, and the temperature dependence of the crystallization effects.Comment: 6 pages, 4 figure

The soft and the hard pomerons in hadron elastic scattering at small t

We consider simple-pole descriptions of soft elastic scattering for pp, pbar p, pi+ p, pi- p, K+ p and K- p. We work at t and s small enough for rescatterings to be neglected, and allow for the presence of a hard pomeron. After building and discussing an exhaustive dataset, we show that simple poles provide an excellent description of the data in the region - 0.5 GeV^2 < t < -0.1 GeV^2, 6 GeV<sqrt(s)< 63 GeV. We show that new form factors have to be used, and get information on the trajectories of the soft and hard pomerons.Comment: 27 pages, 9 figures, LaTeX. A few typos fixed, and references correcte

Capillary filling with wall corrugations] Capillary filling in microchannels with wall corrugations: A comparative study of the Concus-Finn criterion by continuum, kinetic and atomistic approaches

We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods - Computational Fluid Dynamics (CFD), Lattice-Boltzmann Equations (LBE) and Molecular Dynamics (MD). The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles theta > 45. While for theta = 30, theta = 40 (no flow) and theta = 60 (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 the numerical experiments provide different results. Whilst pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle heigth.Comment: 25 pages, 8 figures, Langmuir in pres