Phase separation in thermal systems: A lattice Boltzmann study and morphological characterization

Yanbiao Gan, Aiguo Xu, Guangcai Zhang, Yingjun Li, and Hua Li
Phys. Rev. E 84, 046715 – Published 31 October 2011

Abstract

We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: σ=σ0[(TcT)/(TcT0)]3/2, where Tc is the critical temperature and σ0 is the interfacial stress at a reference temperature T0, which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
8 More
  • Received 14 December 2010

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

©2011 American Physical Society

Authors & Affiliations

Yanbiao Gan1,2,3, Aiguo Xu2,*, Guangcai Zhang2, Yingjun Li1,*, and Hua Li2

  • 1State Key Laboratory for GeoMechanics and Deep Underground Engineering, SMCE, China University of Mining and Technology (Beijing), Beijing 100083, P.R. China
  • 2National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box 8009-26, Beijing 100088, P.R. China
  • 3North China Institute of Aerospace Engineering, Langfang 065000, P.R. China

  • *Corresponding authors: Xu_Aiguo@iapcm.ac.cn; lyj@aphy.iphy.ac.cn

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 84, Iss. 4 — October 2011

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

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
×