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

Nélido González-Segredo, Maziar Nekovee, and Peter V. Coveney
Phys. Rev. E 67, 046304 – Published 21 April 2003
PDFExport Citation

Abstract

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. The model is derivable from a continuous-time Boltzmann-BGK equation in the presence of an intercomponent body force. We find the average domain size grows with time as tγ, where γ increases in the range 0.545±0.014<γ<0.717±0.002, consistent with a crossover between diffusive t1/3 and hydrodynamic viscous, t1.0, behavior. We find good collapse onto a single scaling function, yet the domain growth exponents differ from previous results for similar values of the unique characteristic length L0 and time T0 that can be constructed out of the fluid’s parameters. This rebuts claims of universality for the dynamical scaling hypothesis. For Re=2.7 and small wave numbers q we also find a q2q4 crossover in the scaled structure function, which disappears when the dynamical scaling reasonably improves at later stages (Re=37). This excludes noise as the cause for a q2 behavior, as analytically derived from Yeung and proposed by Appert et al. and Love et al. on the basis of their lattice-gas simulations. We also observe exponential temporal growth of the structure function during the initial stages of the dynamics and for wave numbers less than a threshold value, in accordance with the diffusive Cahn-Hilliard Model B. However, this exponential growth is also present in regimes proscribed by that model. There is no evidence that regions of parameter space for which the scheme is numerically stable become unstable as the simulations proceed, in agreement with finite-difference relaxational models and in contradistinction with an unconditionally unstable lattice-BGK free-energy model previously reported. Those numerical instabilities that do arise in this model are the result of large intercomponent forces which turn the equilibrium distribution negative.

  • Received 28 March 2002

DOI:https://doi.org/10.1103/PhysRevE.67.046304

©2003 American Physical Society

Authors & Affiliations

Nélido González-Segredo*

  • Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom

Maziar Nekovee

  • Complexity Research Group, BT Laboratories, Martlesham Heath, Ipswich, Suffolk IP5 3RE, United Kingdom

Peter V. Coveney

  • Centre for Computational Science, Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom

  • *Also at Grup de Física Estadística, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Spain. Email address: n.gonzalez-segredo@ucl.ac.uk
  • Email address: maziar.nekovee@bt.com
  • Email address: p.v.coveney@ucl.ac.uk

References (Subscription Required)

Click to Expand
Issue

Vol. 67, Iss. 4 — April 2003

Reuse & Permissions
Access Options

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review E

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×