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

Vivien M. Kendon, Michael E. Cates, Ignacio Pagonabarraga, J.C. Desplat, Peter Bladon

Research output: Contribution to journalArticle

219 Citations (Scopus)

Abstract

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 [less, similar] l [less, similar] 105, 10 [less, similar] t [less, similar] 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φ [greater, similar] 100 (which requires t [greater, similar] 106) we find clear evidence of Furukawa's inertial scaling (l [similar] t2/3), although the crossover from the viscous regime (l [similar] t) is both broad and late (102 [less, similar] t [less, similar] 106). Though it cannot be ruled out, we find no indication that Reφ is self-limiting (l [similar] 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 [similar] t2/3 for t [less, similar] 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.
LanguageEnglish
Pages147-203
Number of pages56
JournalJournal of Fluid Mechanics
Volume440
DOIs
Publication statusPublished - 13 Aug 2001

Fingerprint

Spinodal decomposition
binary fluids
decomposition
scaling
Fluids
fluids
numerical control
isotropy
microbalances
statistical analysis
Navier-Stokes equation
Navier Stokes equations
Surface tension
Reynolds number
resources
Statistical methods
crossovers
interfacial tension
indication
simulation

Keywords

  • inertial effects
  • spinodal decomposition
  • symmetric binary fluid mixture
  • lattice Boltzmann
  • Navier-Stokes

Cite this

Kendon, Vivien M. ; Cates, Michael E. ; Pagonabarraga, Ignacio ; Desplat, J.C. ; Bladon, Peter. / Inertial effects in three-dimensional spinodal decomposition of a symmetric binary fluid mixture : a lattice Boltzmann study. In: Journal of Fluid Mechanics. 2001 ; Vol. 440. pp. 147-203.
@article{ef2a3bffa6724e9eaee47975bc97e087,
title = "Inertial effects in three-dimensional spinodal decomposition of a symmetric binary fluid mixture: a lattice Boltzmann study",
abstract = "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 [less, similar] l [less, similar] 105, 10 [less, similar] t [less, similar] 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φ [greater, similar] 100 (which requires t [greater, similar] 106) we find clear evidence of Furukawa's inertial scaling (l [similar] t2/3), although the crossover from the viscous regime (l [similar] t) is both broad and late (102 [less, similar] t [less, similar] 106). Though it cannot be ruled out, we find no indication that Reφ is self-limiting (l [similar] 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 [similar] t2/3 for t [less, similar] 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.",
keywords = "inertial effects, spinodal decomposition, symmetric binary fluid mixture, lattice Boltzmann, Navier-Stokes",
author = "Kendon, {Vivien M.} and Cates, {Michael E.} and Ignacio Pagonabarraga and J.C. Desplat and Peter Bladon",
year = "2001",
month = "8",
day = "13",
doi = "10.1017/S0022112001004682",
language = "English",
volume = "440",
pages = "147--203",
journal = "Journal of Fluid Mechanics",
issn = "0022-1120",
publisher = "Cambridge University Press",

}

Inertial effects in three-dimensional spinodal decomposition of a symmetric binary fluid mixture : a lattice Boltzmann study. / Kendon, Vivien M.; Cates, Michael E.; Pagonabarraga, Ignacio; Desplat, J.C.; Bladon, Peter.

In: Journal of Fluid Mechanics, Vol. 440, 13.08.2001, p. 147-203.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Inertial effects in three-dimensional spinodal decomposition of a symmetric binary fluid mixture

T2 - Journal of Fluid Mechanics

AU - Kendon, Vivien M.

AU - Cates, Michael E.

AU - Pagonabarraga, Ignacio

AU - Desplat, J.C.

AU - Bladon, Peter

PY - 2001/8/13

Y1 - 2001/8/13

N2 - 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 [less, similar] l [less, similar] 105, 10 [less, similar] t [less, similar] 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φ [greater, similar] 100 (which requires t [greater, similar] 106) we find clear evidence of Furukawa's inertial scaling (l [similar] t2/3), although the crossover from the viscous regime (l [similar] t) is both broad and late (102 [less, similar] t [less, similar] 106). Though it cannot be ruled out, we find no indication that Reφ is self-limiting (l [similar] 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 [similar] t2/3 for t [less, similar] 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.

AB - 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 [less, similar] l [less, similar] 105, 10 [less, similar] t [less, similar] 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φ [greater, similar] 100 (which requires t [greater, similar] 106) we find clear evidence of Furukawa's inertial scaling (l [similar] t2/3), although the crossover from the viscous regime (l [similar] t) is both broad and late (102 [less, similar] t [less, similar] 106). Though it cannot be ruled out, we find no indication that Reφ is self-limiting (l [similar] 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 [similar] t2/3 for t [less, similar] 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.

KW - inertial effects

KW - spinodal decomposition

KW - symmetric binary fluid mixture

KW - lattice Boltzmann

KW - Navier-Stokes

U2 - 10.1017/S0022112001004682

DO - 10.1017/S0022112001004682

M3 - Article

VL - 440

SP - 147

EP - 203

JO - Journal of Fluid Mechanics

JF - Journal of Fluid Mechanics

SN - 0022-1120

ER -