Japan Geoscience Union Meeting 2015

Presentation information

Poster

Symbol S (Solid Earth Sciences) » S-VC Volcanology

[S-VC46] Dynamics of volcanic eruptions and igneous activities

Mon. May 25, 2015 6:15 PM - 7:30 PM Convention Hall (2F)

Convener:*Tomofumi Kozono(Department of Geophysics, Graduate School of Science, Tohoku University), Yujiro Suzuki(Earthquake Research Institute, The University of Tokyo), Satoshi Okumura(Division of Earth and Planetary Materials Science, Department of Earth Science, Graduate School of Science, Tohoku University)

6:15 PM - 7:30 PM

[SVC46-P06] Direct molecular dynamics simulations of homogeneous bubble nucleation and improvements of classical theory

*Kyoko TANAKA1, Hidekazu TANAKA1, Juerg DIEMAND2, Raymond ANGELIL2 (1.Institute of Low Temperature Science, Hokkaido University, 2.University of Zurich)

Keywords:bubble nucleation, liquid to vapor transition, phase transition, molecular dynamics simulation

Bubble nucleation in liquid is a liquid-to-vapor transition phenomenon and plays an important role in vulcanism. Studies of homogeneous liquid-vapour nucleation typically use the classical formula (CNT) for the bubble nucleation rate. However, the applicability of the CNT is not well understood.

Numerical techniques such as molecular dynamics and Monte-Carlo simulations are powerful methods to resolve details of the nucleation process and provide useful test cases for nucleation models. Typically, these simulations show large deviations from the CNT predictions. Most of the simulations for bubble nucleation in the literature use around 105 or fewer atoms, making it difficult to measure nucleation rates directly.

Recently, we presented large-scale, micro-canonical molecular dynamics simulations of homogeneous bubble nucleation with 5 * 108 Lennard-Jones atoms, and succeeded to directly measure nucleation rates in the range of 1021-25cm-3s-1 for argon by resolving bubble nucleation events in the steady state nucleation phase [1,2]. The unprecedented size of the simulated volumes allows us to resolve the nucleation and growth of many bubbles per run in simple direct micro-canonical (NVE) simulations while the ambient pressure and temperature remain almost perfectly constant.

We find bubble nucleation rates which are lower than in most of the previous, smaller simulations. It is widely believed that classical nucleation theory (CNT) generally underestimates bubble nucleation rates by very large factors. However, our measured rates are within two orders of magnitude of CNT predictions - only at very low temperatures the CNT underestimates the nucleation rate significantly.

We also derive an improved classical formula for the homogeneous bubble nucleation rate, where we revise the prefactor in the nucleation rate and compare it with the widely used classical nucleation theory (CNT) [3]. Our large-scale molecular dynamics simulations and laboratory experiments for argon bubble nucleation enable us to precisely test our theoretical models. The improved formula including the Tolman correction with a small positive Tolman length leads to good agreement with both MD simulations and laboratory experiments.

[1] J. Diemand, R. Angelil, K. K. Tanaka, and H. Tanaka, Pys. Rev E 90, 052407 (2014)
[2] R. Angelil, J. Diemand, K. K. Tanaka, and H. Tanaka, Pys. Rev E 90, 063301 (2014)
[3] K. K. Tanaka, H. Tanaka, R. Angelil, and J. Diemand, submitted