Inertial effects in three-dimensional spinodal decomposition of a symmetric binary fluid mixture: a lattice Boltzmann study

The late-stage demixing following spinodal decomposition of a three-dimensional symmetric binary fluid mixture is studied numerically, using a thermodynamically consistent lattice Boltzmann method. We combine results from simulations with different numerical parameters to obtain an unprecedented range of length and time scales when expressed in reduced physical units. (These are the length and time units derived from fluid density, viscosity, and interfacial tension.) Using eight large (2563) runs, the resulting composite graph of reduced domain size l against reduced time t covers 1 [lsim ] l [lsim ] 105, 10 [lsim ] t [lsim ] 108. Our data are consistent with the dynamical scaling hypothesis that l(t) is a universal scaling curve. We give the first detailed statistical analysis of fluid motion, rather than just domain evolution, in simulations of this kind, and introduce scaling plots for several quantities derived from the fluid velocity and velocity gradient fields. Using the conventional definition of Reynolds number for this problem, Reϕ = ldl/dt, we attain values approaching 350. At Reϕ [gsim ] 100 (which requires t [gsim ] 106) we find clear evidence of Furukawa's inertial scaling (l ∼ t2/3), although the crossover from the viscous regime (l ∼ t) is both broad and late (102 [lsim ] t [lsim ] 106). Though it cannot be ruled out, we find no indication that Reϕ is self-limiting (l ∼ t1/2) at late times, as recently proposed by Grant & Elder. Detailed study of the velocity fields confirms that, for our most inertial runs, the RMS ratio of nonlinear to viscous terms in the Navier–Stokes equation, R2, is of order 10, with the fluid mixture showing incipient turbulent characteristics. However, we cannot go far enough into the inertial regime to obtain a clear length separation of domain size, Taylor microscale, and Kolmogorov scale, as would be needed to test a recent ‘extended’ scaling theory of Kendon (in which R2 is self-limiting but Reϕ not). Obtaining our results has required careful steering of several numerical control parameters so as to maintain adequate algorithmic stability, efficiency and isotropy, while eliminating unwanted residual diffusion. (We argue that the latter affects some studies in the literature which report l ∼ t2/3 for t [lsim ] 104.) We analyse the various sources of error and find them just within acceptable levels (a few percent each) in most of our datasets. To bring these under significantly better control, or to go much further into the inertial regime, would require much larger computational resources and/or a breakthrough in algorithm design.

[1]  T. N. Stevenson,et al.  Fluid Mechanics , 2021, Nature.

[2]  M. Cates,et al.  Interfacial dynamics in 3D binary fluid demixing: animation studies* , 2001 .

[3]  Kendon Scaling theory of three-dimensional spinodal turbulence , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[4]  P. Coveney,et al.  Entropic lattice Boltzmann methods , 2000, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[5]  Peter Bladon Edinburgh Parallel Computing Centre,et al.  LUDWIG: A parallel Lattice-Boltzmann code for complex fluids , 2000, cond-mat/0005213.

[6]  Olvera de la Cruz M,et al.  Hydrodynamic coarsening of binary fluids , 2000, Physical review letters.

[7]  M. Cates,et al.  Inertia, coarsening and fluid motion in binary mixtures , 1999, cond-mat/9908007.

[8]  J. Yeomans,et al.  Phase separation in two-dimensional fluids: the role of noise. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[9]  M. E. Cates,et al.  3D spinodal decomposition in the inertial regime , 1999, cond-mat/9902346.

[10]  M. Grant,et al.  Spinodal Decomposition in Fluids , 1999 .

[11]  Dimitrios V. Rovas,et al.  An Improved Hydrodynamics Formulation for Multiphase Flow Lattice-Boltzmann Models , 1998 .

[12]  Leeds,et al.  Tests of dynamical scaling in three-dimensional spinodal decomposition , 1998, cond-mat/9810298.

[13]  D. Mccomb,et al.  Investigation of renormalization group methods for the numerical simulation of isotropic turbulence , 1998, physics/9806028.

[14]  G. Gompper,et al.  Lattice-Boltzmann model of amphiphilic systems , 1998, cond-mat/9803097.

[15]  Paul M. Chaikin,et al.  Long-range correlations in sedimentation , 1997 .

[16]  Gerard Cambon,et al.  Declic: A facility to investigate critical fluids behavior, chemical reactivity in supercritical water and directional solidification in transparent media , 1997 .

[17]  Sidney R. Nagel,et al.  Breakdown of scaling in droplet fission at high Reynolds number , 1997 .

[18]  J. Lebowitz,et al.  Spinodal Decomposition in Binary Gases , 1997, patt-sol/9703006.

[19]  Yeomans,et al.  Lattice Boltzmann simulations of liquid-gas and binary fluid systems. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[20]  Laradji,et al.  Molecular Dynamics Simulation of Spinodal Decomposition in Three-Dimensional Binary Fluids. , 1996, Physical review letters.

[21]  U. Frisch Turbulence: The Legacy of A. N. Kolmogorov , 1996 .

[22]  Chen,et al.  Spinodal decomposition in fluids: Diffusive, viscous, and inertial regimes. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[23]  Stéphane Zaleski,et al.  Phase separation in a three-dimensional, two-phase, hydrodynamic lattice gas , 1995 .

[24]  T. Hashimoto,et al.  Late stage spinodal decomposition in binary critical fluids: scaling function obtained over a wide q-space of 4 orders of magnitude , 1994 .

[25]  A. Ladd Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 1. Theoretical foundation , 1993, Journal of Fluid Mechanics.

[26]  Chen,et al.  Hydrodynamic spinodal decomposition: Growth kinetics and scaling functions. , 1993, Physical review. B, Condensed matter.

[27]  A. Bray Theory of phase-ordering kinetics , 1993, cond-mat/9501089.

[28]  Puri,et al.  Temporally linear domain growth in the segregation of binary fluids. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[29]  Ma,et al.  Dynamics of phase separation of binary fluids. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[30]  Kuwahara,et al.  Spinodal decomposition in a critical isobutyric acid and water mixture. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[31]  Koga,et al.  Spinodal decomposition in binary fluids: Effects of hydrodynamic interactions. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[32]  J. Brasseur,et al.  THE RESPONSE OF ISOTROPIC TURBULENCE TO ISOTROPIC AND ANISOTROPIC FORCING AT THE LARGE SCALES , 1991 .

[33]  Shinozaki,et al.  Asymptotic form factor for spinodal decomposition in three-space. , 1991, Physical review letters.

[34]  R. Benzi,et al.  Lattice Gas Dynamics with Enhanced Collisions , 1989 .

[35]  Stephen B. Pope,et al.  An examination of forcing in direct numerical simulations of turbulence , 1988 .

[36]  Perrot,et al.  Spinodal decomposition patterns in an isodensity critical binary fluid: Direct-visualization and light-scattering analyses. , 1987, Physical review. A, General physics.

[37]  Furukawa Effect of inertia on droplet growth in a fluid. , 1985, Physical review. A, General physics.

[38]  Eric D. Siggia,et al.  Late stages of spinodal decomposition in binary mixtures , 1979 .

[39]  Simon Ian Jury,et al.  Computer simulation of complex fluids using dissipative particle dynamics , 1999 .

[40]  M. Cates,et al.  Spinodal decomposition in the viscous hydrodynamic regime , 1999 .

[41]  D. Frenkel,et al.  Algebraic Decay of Velocity Fluctuations in a Confined Fluid , 1997 .

[42]  S. H. A. Chen,et al.  Small-angle light scattering in microemulsions (spinodal decomposition) , 1993 .

[43]  C. Brooks Computer simulation of liquids , 1989 .