MATHEMATICAL MODELING OF NON-EQUILIBRIUM PHASE TRANSITION IN RAPIDLY HEATED THIN LIQUID FILM

In the framework of molecular dynamics the behavior of thin liquid film during its rapid homogeneous heating with rate 2-100 K/ ps is studied. The initial film thickness through the z axis is 48 nm. The overall number of particles is 96 000 and the length of calculation area is 268 nm. In the film plane periodic boundary conditions were used with the main dimension 8x8 nm. The calculation results are presented in a form of distributions of temperature, density and particle velocity in the calculation area averaged over the film plane and also in the form of two-dimension particle density distributions in the zx plane in the range of heating time from 0 to 800 ps. The obtained results suggest the existence of four different regimes of film behavior depending on the heating rate: quasi-stationary regime with surface evaporation, explosive (volume) boiling, spinodal decomposition and overcritical expansion.