导航

Acta Aeronautica et Astronautica Sinica ›› 2026, Vol. 47 ›› Issue (4): 232342.doi: 10.7527/S1000-6893.2025.32342

• Solid Mechanics and Vehicle Conceptual Design • Previous Articles    

Molecular dynamics simulation of stress and temperature effects on supercooled water freezing

Ruiqi HAN1, Fei XU1,2(), Cheng CHEN1, Zhanpeng REN3, Hui ZHANG3   

  1. 1.School of Aeronautics,Northwestern Polytechnical University,Xi’an 710072,China
    2.National Key Laboratory of Strength and Structural Integrity,Xi’an 710072,China
    3.Aircraft Strength Research Institute of China,Xi’an 710065,China
  • Received:2025-05-30 Revised:2025-07-30 Accepted:2025-08-18 Online:2025-09-08 Published:2025-08-28
  • Contact: Fei XU E-mail:xufei@nwpu.edu.cn
  • Supported by:
    National Natural Science Foundation of China(D5110220438)

Abstract:

The physical process of supercooled large droplet impingement freezing is complex, involving the coupling of droplet impact dynamics and solidification process, among which the influence of temperature and stress on ice nucleus formation and growth has received widespread attention. The molecular dynamics simulation method is employed to establish homogeneous nucleation and heterogeneous nucleation models for the freezing problem of supercooled water at different temperatures, and the natural freezing and ice nucleation freezing processes are studied. First, the radial distribution function characteristics are analyzed after the freezing fraction reached 70%. It is preliminarily determined that when the radial distribution function of the research object is consistent with cubic ice and the deviations at the two peaks at 2.7×10-10 m and 4.4×10-10 m are less than 10% and 20% respectively, it can be used as a basis for judging the freezing phenomenon. Then, the critical number of ice nuclei for heterogeneous nucleation models at different temperatures is calculated, which is consistent with the critical number of cubic ice nuclei given by the classical nucleation theory at temperatures below 250 K. On this basis, load stress is introduced to study its influence on the icing process. The results of this paper show that under the action of different load, the strain rates promoting nucleation and growth of ice are different. The shear strain rate is 107 - 108 s-1, and the triaxial compression strain rate is 105 - 106 s-1. The nucleation process of ice is mainly due to the competition between the hindering effect caused by the increase of nucleation barrier and the promoting effect caused by the increase of self-diffusion rate, and the growth process of ice nuclei is mainly attributed to the promoting effect caused by the increase of self-diffusion rate.

Key words: supercooled water, ice, temperature, stress, molecular dynamics simulation

CLC Number: