Three dimensional hysdrodynamic lattice-gas simulations of binary immiscible and ternary amphiphilic flow through porous media

We report the results of a study of multiphase flow in porous media. A Darcy's law for steady multiphase flow was investigated for both binary and ternary amphiphilic flow. Linear flux-forcing relationships satisfying Onsager reciprocity were shown to be a good approximation of the simulation data. The dependence of the relative permeability coefficients on water saturation was investigated and showed good qualitative agreement with experimental data. Non-steady state invasion flows were investigated, with particular interest in the asymptotic residual oil saturation. The addition of surfactant to the invasive fluid was shown to significantly reduce the residual oil saturation.Comment: To appear in Phys. Rev.

Lattice Boltzmann Thermohydrodynamics

We introduce a lattice Boltzmann computational scheme capable of modeling thermohydrodynamic flows of monatomic gases. The parallel nature of this approach provides a numerically efficient alternative to traditional methods of computational fluid dynamics. The scheme uses a small number of discrete velocity states and a linear, single-time-relaxation collision operator. Numerical simulations in two dimensions agree well with exact solutions for adiabatic sound propagation and Couette flow with heat transfer.Comment: 11 pages, Physical Review E: Rapid Communications, in pres

Lattice Boltzmann Simulation of Non-Ideal Fluids

A lattice Boltzmann scheme able to model the hydrodynamics of phase separation and two-phase flow is described. Thermodynamic consistency is ensured by introducing a non-ideal pressure tensor directly into the collision operator. We also show how an external chemical potential can be used to supplement standard boundary conditions in order to investigate the effect of wetting on phase separation and fluid flow in confined geometries. The approach has the additional advantage of reducing many of the unphysical discretisation problems common to previous lattice Boltzmann methods.Comment: 11 pages, revtex, 4 Postscript figures, uuencode

Hydrodynamic Spinodal Decomposition: Growth Kinetics and Scaling Functions

We examine the effects of hydrodynamics on the late stage kinetics in spinodal decomposition. From computer simulations of a lattice Boltzmann scheme we observe, for critical quenches, that single phase domains grow asymptotically like $t^{\alpha}$, with $\alpha \approx .66$ in two dimensions and $\alpha \approx 1.0$ in three dimensions, both in excellent agreement with theoretical predictions.Comment: 12 pages, latex, Physical Review B Rapid Communication (in press

Multi-component lattice-Boltzmann model with interparticle interaction

A previously proposed [X. Shan and H. Chen, Phys. Rev. E {\bf 47}, 1815, (1993)] lattice Boltzmann model for simulating fluids with multiple components and interparticle forces is described in detail. Macroscopic equations governing the motion of each component are derived by using Chapman-Enskog method. The mutual diffusivity in a binary mixture is calculated analytically and confirmed by numerical simulation. The diffusivity is generally a function of the concentrations of the two components but independent of the fluid velocity so that the diffusion is Galilean invariant. The analytically calculated shear kinematic viscosity of this model is also confirmed numerically.Comment: 18 pages, compressed and uuencoded postscript fil

Domain Growth, Wetting and Scaling in Porous Media

The lattice Boltzmann (LB) method is used to study the kinetics of domain growth of a binary fluid in a number of geometries modeling porous media. Unlike the traditional methods which solve the Cahn-Hilliard equation, the LB method correctly simulates fluid properties, phase segregation, interface dynamics and wetting. Our results, based on lattice sizes of up to $4096\times 4096$, do not show evidence to indicate the breakdown of late stage dynamical scaling, and suggest that confinement of the fluid is the key to the slow kinetics observed. Randomness of the pore structure appears unnecessary.Comment: 13 pages, latex, submitted to PR

Simulation of Flow of Mixtures Through Anisotropic Porous Media using a Lattice Boltzmann Model

We propose a description for transient penetration simulations of miscible and immiscible fluid mixtures into anisotropic porous media, using the lattice Boltzmann (LB) method. Our model incorporates hydrodynamic flow, diffusion, surface tension, and the possibility for global and local viscosity variations to consider various types of hardening fluids. The miscible mixture consists of two fluids, one governed by the hydrodynamic equations and one by diffusion equations. We validate our model on standard problems like Poiseuille flow, the collision of a drop with an impermeable, hydrophobic interface and the deformation of the fluid due to surface tension forces. To demonstrate the applicability to complex geometries, we simulate the invasion process of mixtures into wood spruce samples.Comment: Submitted to EPJ

Fluctuations of elastic interfaces in fluids: Theory and simulation

We study the dynamics of elastic interfaces-membranes-immersed in thermally excited fluids. The work contains three components: the development of a numerical method, a purely theoretical approach, and numerical simulation. In developing a numerical method, we first discuss the dynamical coupling between the interface and the surrounding fluids. An argument is then presented that generalizes the single-relaxation time lattice-Boltzmann method for the simulation of hydrodynamic interfaces to include the elastic properties of the boundary. The implementation of the new method is outlined and it is tested by simulating the static behavior of spherical bubbles and the dynamics of bending waves. By means of the fluctuation-dissipation theorem we recover analytically the equilibrium frequency power spectrum of thermally fluctuating membranes and the correlation function of the excitations. Also, the non-equilibrium scaling properties of the membrane roughening are deduced, leading us to formulate a scaling law describing the interface growth, W^2(L,T)=L^3 g[t/L^(5/2)], where W, L and T are the width of the interface, the linear size of the system and the temperature respectively, and g is a scaling function. Finally, the phenomenology of thermally fluctuating membranes is simulated and the frequency power spectrum is recovered, confirming the decay of the correlation function of the fluctuations. As a further numerical study of fluctuating elastic interfaces, the non-equilibrium regime is reproduced by initializing the system as an interface immersed in thermally pre-excited fluids.Comment: 15 pages, 11 figure

Three-dimensional lattice-Boltzmann simulations of critical spinodal decomposition in binary immiscible fluids

We use a modified Shan-Chen, noiseless lattice-BGK model for binary immiscible, incompressible, athermal fluids in three dimensions to simulate the coarsening of domains following a deep quench below the spinodal point from a symmetric and homogeneous mixture into a two-phase configuration. We find the average domain size growing with time as $t^\gamma$, where $\gamma$ increases in the range $0.545 < \gamma < 0.717$, consistent with a crossover between diffusive $t^{1/3}$ and hydrodynamic viscous, $t^{1.0}$, behaviour. We find good collapse onto a single scaling function, yet the domain growth exponents differ from others' works' for similar values of the unique characteristic length and time that can be constructed out of the fluid's parameters. This rebuts claims of universality for the dynamical scaling hypothesis. At early times, we also find a crossover from $q^2$ to $q^4$ in the scaled structure function, which disappears when the dynamical scaling reasonably improves at later times. This excludes noise as the cause for a $q^2$ behaviour, as proposed by others. We also observe exponential temporal growth of the structure function during the initial stages of the dynamics and for wavenumbers less than a threshold value.Comment: 45 pages, 18 figures. Accepted for publication in Physical Review

Numerical study of wetting transitions on biomimetic surfaces using a lattice Boltzmann approach with large density ratio

The hydrophobicity of natural surfaces have drawn much attention of scientific communities in recent years. By mimicking natural surfaces, the manufactured biomimetic hydrophobic surfaces have been widely applied to green technologies such as self-cleaning surfaces. Although the theories for wetting and hydrophobicity have been developed, the mechanism of wetting transitions between heterogeneous wetting state and homogeneous wetting state is still not fully clarified. As understanding of wetting transitions is crucial for manufacturing a biomimetic superhydrophobic surface, more fundamental discussions in this area should be carried out. In the present work the wetting transitions are numerically studied using a phase field lattice Boltzmann approach with large density ratio, which should be helpful in understanding the mechanism of wetting transitions. The dynamic wetting transition processes between Cassie-Baxter state and Wenzel state are presented, and the energy barrier and the gravity effect on transition are discussed. It is found that the two wetting transition processes are irreversible for specific inherent contact angles and have different transition routes, the energy barrier exists on an ideally patterned surface and the gravity can be crucial to overcome the energy barrier and trigger the transition