Generalized Lattice Boltzmann Method with multi-range pseudo-potential

The physical behaviour of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudo-potential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or high order isotropy. Higher order schemes to implement self-consistent forcings are rigorously computed for 2d and 3d models. The extended scenario developed in this work clarifies the theoretical foundations of the Shan-Chen methodology for the lattice Boltzmann method and enhances its applicability and flexibility to the simulation of multiphase flows to density ratios up to O(100)

A planar surface acoustic wave micropump for closed-loop microfluidics

We have designed and characterized a simple Rayleigh-surface acoustic wave-based micropump, integrated directly with a fully enclosed 3D microfluidic system, which improves significantly the pumping efficiency within a coupled fluid whilst maintaining planar integration of the micropump and microfluidics. We achieve this by exploiting the Rayleigh-scattering angle of surface acoustic waves into pressure waves on contact with overlaid fluids, by designing a microfluidic channel aligned almost co-linearly with the launched pressure waves and by minimizing energy losses by reflections from, or absorption within, the channel walls. This allows the microfluidic system to remain fully enclosed—a pre-requisite for point-of-care applications—removing sources of possible contamination, whilst achieving pump efficiencies up to several orders of magnitude higher than previously reported, at low operating powers of 0.5 W

Reaction-diffusion dynamics: confrontation between theory and experiment in a microfluidic reactor

We confront, quantitatively, the theoretical description of the reaction-diffusion of a second order reaction to experiment. The reaction at work is \ca/CaGreen, and the reactor is a T-shaped microchannel, 10 $\mu$m deep, 200 $\mu$m wide, and 2 cm long. The experimental measurements are compared with the two-dimensional numerical simulation of the reaction-diffusion equations. We find good agreement between theory and experiment. From this study, one may propose a method of measurement of various quantities, such as the kinetic rate of the reaction, in conditions yet inaccessible to conventional methods

Intermittency and the Slow Approach to Kolmogorov Scaling

From a simple path integral involving a variable volatility in the velocity differences, we obtain velocity probability density functions with exponential tails, resembling those observed in fully developed turbulence. The model yields realistic scaling exponents and structure functions satisfying extended self-similarity. But there is an additional small scale dependence for quantities in the inertial range, which is linked to a slow approach to Kolmogorov (1941) scaling occurring in the large distance limit.Comment: 10 pages, 5 figures, minor changes to mirror version to appear in PR

Mesoscopic two-phase model for describing apparent slip in micro-channel flows

The phenomenon of apparent slip in micro-channel flows is analyzed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactins. The weakly-inhomogeneous limit of this model is solved analytically. The present mesoscopic approach permits to access much larger scales than molecular dynamics, and comparable with those attained by continuum methods. However, at variance with the continuum approach, the existence of a gas layer near the wall does not need to be postulated a priori, but emerges naturally from the underlying non-ideal mesoscopic dynamics. It is therefore argued that a mesoscopic Lattice Boltzmann approach with non-ideal fluid-fluid and fluid-wall interactions might achieve an optimal compromise between physical realism and computational efficiency for the study of channel micro-flows.Comment: 5 pages, 3 figure

Numerical Simulation of Vortex Crystals and Merging in N-Point Vortex Systems with Circular Boundary

In two-dimensional (2D) inviscid incompressible flow, low background vorticity distribution accelerates intense vortices (clumps) to merge each other and to array in the symmetric pattern which is called ``vortex crystals''; they are observed in the experiments on pure electron plasma and the simulations of Euler fluid. Vortex merger is thought to be a result of negative ``temperature'' introduced by L. Onsager. Slight difference in the initial distribution from this leads to ``vortex crystals''. We study these phenomena by examining N-point vortex systems governed by the Hamilton equations of motion. First, we study a three-point vortex system without background distribution. It is known that a N-point vortex system with boundary exhibits chaotic behavior for N\geq 3. In order to investigate the properties of the phase space structure of this three-point vortex system with circular boundary, we examine the Poincar\'e plot of this system. Then we show that topology of the Poincar\'e plot of this system drastically changes when the parameters, which are concerned with the sign of ``temperature'', are varied. Next, we introduce a formula for energy spectrum of a N-point vortex system with circular boundary. Further, carrying out numerical computation, we reproduce a vortex crystal and a vortex merger in a few hundred point vortices system. We confirm that the energy of vortices is transferred from the clumps to the background in the course of vortex crystallization. In the vortex merging process, we numerically calculate the energy spectrum introduced above and confirm that it behaves as k^{-\alpha},(\alpha\approx 2.2-2.8) at the region 10^0<k<10^1 after the merging.Comment: 30 pages, 11 figures. to be published in Journal of Physical Society of Japan Vol.74 No.

Quasi-Gaussian Statistics of Hydrodynamic Turbulence in 3/4+\epsilon dimensions

The statistics of 2-dimensional turbulence exhibit a riddle: the scaling exponents in the regime of inverse energy cascade agree with the K41 theory of turbulence far from equilibrium, but the probability distribution functions are close to Gaussian like in equilibrium. The skewness \C S \equiv S_3(R)/S^{3/2}_2(R) was measured as \C S_{\text{exp}}\approx 0.03. This contradiction is lifted by understanding that 2-dimensional turbulence is not far from a situation with equi-partition of enstrophy, which exist as true thermodynamic equilibrium with K41 exponents in space dimension of $d=4/3$. We evaluate theoretically the skewness \C S(d) in dimensions ${4/3}\le d\le 2$, show that \C S(d)=0 at $d=4/3$, and that it remains as small as \C S_{\text{exp}} in 2-dimensions.Comment: PRL, submitted, REVTeX 4, 4 page

Mesoscopic modeling of a two-phase flow in the presence of boundaries: the Contact Angle

We present a mesoscopic model, based on the Boltzmann Equation, for the interaction between a solid wall and a non-ideal fluid. We present an analytic derivation of the contact angle in terms of the surface tension between the liquid-gas, the liquid-solid and the gas-solid phases. We study the dependency of the contact angle on the two free parameters of the model, which determine the interaction between the fluid and the boundaries, i.e. the equivalent of the wall density and of the wall-fluid potential in Molecular Dynamics studies. We compare the analytical results obtained in the hydrodynamical limit for the density profile and for the surface tension expression with the numerical simulations. We compare also our two-phase approach with some exact results for a pure hydrodynamical incompressible fluid based on Navier-Stokes equations with boundary conditions made up of alternating slip and no-slip strips. Finally, we show how to overcome some theoretical limitations connected with a discretized Boltzmann scheme and we discuss the equivalence between the surface tension defined in terms of the mechanical equilibrium and in terms of the Maxwell construction.Comment: 29 pages, 12 figure

Studying Flow Close to an Interface by Total Internal Reflection Fluorescence Cross Correlation Spectroscopy: Quantitative Data Analysis

Total Internal Reflection Fluorescence Cross Correlation Spectroscopy (TIR-FCCS) has recently (S. Yordanov et al., Optics Express 17, 21149 (2009)) been established as an experimental method to probe hydrodynamic flows near surfaces, on length scales of tens of nanometers. Its main advantage is that fluorescence only occurs for tracer particles close to the surface, thus resulting in high sensitivity. However, the measured correlation functions only provide rather indirect information about the flow parameters of interest, such as the shear rate and the slip length. In the present paper, we show how to combine detailed and fairly realistic theoretical modeling of the phenomena by Brownian Dynamics simulations with accurate measurements of the correlation functions, in order to establish a quantitative method to retrieve the flow properties from the experiments. Firstly, Brownian Dynamics is used to sample highly accurate correlation functions for a fixed set of model parameters. Secondly, these parameters are varied systematically by means of an importance-sampling Monte Carlo procedure in order to fit the experiments. This provides the optimum parameter values together with their statistical error bars. The approach is well suited for massively parallel computers, which allows us to do the data analysis within moderate computing times. The method is applied to flow near a hydrophilic surface, where the slip length is observed to be smaller than 10nm, and, within the limitations of the experiments and the model, indistinguishable from zero.Comment: 18 pages, 12 figure

Mean- Field Approximation and a Small Parameter in Turbulence Theory

Numerical and physical experiments on two-dimensional (2d) turbulence show that the differences of transverse components of velocity field are well described by a gaussian statistics and Kolmogorov scaling exponents. In this case the dissipation fluctuations are irrelevant in the limit of small viscosity. In general, one can assume existence of critical space-dimensionality $d=d_{c}$, at which the energy flux and all odd-order moments of velocity difference change sign and the dissipation fluctuations become dynamically unimportant. At $d<d_{c}$ the flow can be described by the ``mean-field theory'', leading to the observed gaussian statistics and Kolmogorov scaling of transverse velocity differences. It is shown that in the vicinity of $d=d_{c}$ the ratio of the relaxation and translation characteristic times decreases to zero, thus giving rise to a small parameter of the theory. The expressions for pressure and dissipation contributions to the exact equation for the generating function of transverse velocity differences are derived in the vicinity of $d=d_{c}$. The resulting equation describes experimental data on two-dimensional turbulence and demonstrate onset of intermittency as $d-d_{c}>0$ and $r/L\to 0$ in three-dimensional flows in close agreement with experimental data. In addition, some new exact relations between correlation functions of velocity differences are derived. It is also predicted that the single-point pdf of transverse velocity difference in developing as well as in the large-scale stabilized two-dimensional turbulence is a gaussian.Comment: 25 pages, 1 figur