We present a numerical method for simulations of spinodal decomposition of liquid-vapor systems. The results are in excellent agreement with theoretical predictions for all expected time regimes from the initial growth of “homophase fluctuations” up to the inertial hydrodynamics regime. The numerical approach follows a modern formulation of the smoothed particle hydrodynamics method with a van der Waals equation of state and thermal conduction. The dynamics and thermal evolution of instantaneously temperature-quenched systems are investigated. Therefore, we introduce a simple scaling thermostat that allows thermal fluctuations at a constant predicted mean temperature. We find that the initial stage spinodal decomposition is strongly affected by the temperature field. The separated phases react on density changes with a change in temperature. Although, the thermal conduction acts very slowly, thermal deviations are eventually compensated. The domain growth in the late stage of demixing is found to be rather unaffected by thermal fluctuations. We observe a transition from the Lifshitz-Slyozov growth rate with $1/3$ exponent to the inertial hydrodynamics regime with a rate of $2/3$, only excepted from simulations near the critical point where the liquid droplets are observed to nucleate directly in a spherical shape. The transition between the growth regimes is found to occur earlier for higher initial temperatures. We explain this time dependency with the phase interfaces that become more diffuse and overlap with approaching the critical point. A prolonging behavior of the demixing process is observed and also expected to depend on temperature. It is further found that the observations can excellently explain the growth behavior for pure nonisothermal simulations that are performed without thermostat.

• Received 22 September 2014
• Revised 5 February 2015

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