固体力学与飞行器总体设计

应力和温度对过冷水结冰影响的分子动力学模拟

  • 汉芮岐 ,
  • 徐绯 ,
  • 陈成 ,
  • 任战鹏 ,
  • 张惠
展开
  • 1.西北工业大学 航空学院,西安 710072
    2.强度与结构完整性全国重点实验室,西安 710072 3. 中国飞机强度研究所,西安 710065
E-mail: xufei@nwpu.edu.cn

收稿日期: 2025-05-30

  修回日期: 2025-07-30

  录用日期: 2025-08-18

  网络出版日期: 2025-08-28

基金资助

国家自然科学基金(D5110220438)

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

  • Ruiqi HAN ,
  • Fei XU ,
  • Cheng CHEN ,
  • Zhanpeng REN ,
  • Hui ZHANG
Expand
  • 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
E-mail: xufei@nwpu.edu.cn

Received date: 2025-05-30

  Revised date: 2025-07-30

  Accepted date: 2025-08-18

  Online published: 2025-08-28

Supported by

National Natural Science Foundation of China(D5110220438)

摘要

过冷大液滴碰撞结冰的物理过程复杂,涉及液滴冲击动力学与凝固过程的耦合,其中温度、应力对冰形核、生长的影响引起了广泛的关注。采用分子动力学模拟方法,针对不同温度下过冷水结冰问题,建立了均匀形核、非均匀形核模型,研究了自然结冰、冰核结冰过程。分析了结冰比例达到70%后的径向分布函数特征,初步确定当研究对象的径向分布函数与立方冰一致,且在2.7×10-10、4.4×10-10 m 2个峰值处偏差分别小于10%、20%时可判断为结冰现象;计算了不同温度下非均匀形核模型的冰核临界数,在低于250 K温度下与经典形核理论给出的立方冰核临界数吻合;引入了载荷应力研究其对结冰过程的影响。结果表明,在剪切载荷、三向挤压载荷作用下,促进冰的形核、冰核生长的应变率不同,剪应变率为107~108 s-1,三向压缩应变率为105~106 s-1。冰的形核过程主要是由于成核势垒增大引起的阻碍效应与自扩散率增大引起的促进效应之间的竞争,冰核生长过程主要归因于自扩散率增大引起的促进效应。

本文引用格式

汉芮岐 , 徐绯 , 陈成 , 任战鹏 , 张惠 . 应力和温度对过冷水结冰影响的分子动力学模拟[J]. 航空学报, 2026 , 47(4) : 232342 -232342 . DOI: 10.7527/S1000-6893.2025.32342

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.

参考文献

[1] WU Z L. Drop “impact” on an airfoil surface[J]. Advances in Colloid and Interface Science2018256: 23-47.
[2] BELLOSTA T, BALDAN G, SIRIANNI G, et al. Lagrangian and Eulerian algorithms for water droplets in in-flight ice accretion[J]. Journal of Computational and Applied Mathematics2023429: 115230.
[3] CHENG X, SUN T-P, GORDILLO L. Drop impact dynamics: impact force and stress distributions[J]. Annual Review of Fluid Mechanics202254: 57-81.
[4] SUN T P, áLVAREZ-NOVOA F, ANDRADE K, et al. Stress distribution and surface shock wave of drop impact[J]. Nature Communications202213: 1703.
[5] THANH-VINH N, MATSUMOTO K, SHIMOYAMA I. Pressure distribution on the contact area during the impact of a droplet on a textured surface[C]∥2016 IEEE 29th International Conference on Micro Electro Mechanical Systems (MEMS). Piscataway: IEEE Press, 2016: 177-180.
[6] PHILIPPI J, LAGRéE P Y, ANTKOWIAK A. Drop impact on a solid surface: Short-time self-similarity[J]. Journal of Fluid Mechanics2016795: 96-135.
[7] ZHANG C, LIU H. Effect of drop size on the impact thermodynamics for supercooled large droplet in aircraft icing[J]. Physics of Fluids201628(6): 062107.
[8] YANG G M, GUO K H, LI N. Freezing mechanism of supercooled water droplet impinging on metal surfaces[J]. International Journal of Refrigeration201134(8): 2007-2017.
[9] JIN Z Y, WANG Z N, SUI D Y, et al. The impact and freezing processes of a water droplet on different inclined cold surfaces[J]. International Journal of Heat and Mass Transfer201697: 211-223.
[10] FANG W Z, ZHU F Q, TAO W Q, et al. How different freezing morphologies of impacting droplets form[J]. Journal of Colloid and Interface Science2021584: 403-410.
[11] WANG L P, KONG W L, WANG F X, et al. Effect of nucleation time on freezing morphology and type of a water droplet impacting onto cold substrate[J]. International Journal of Heat and Mass Transfer2019130: 831-842.
[12] SUN M M, KONG W L, WANG F X, et al. Impact freezing modes of supercooled droplets determined by both nucleation and icing evolution[J]. International Journal of Heat and Mass Transfer2019142: 118431.
[13] WANG L P, WANG F X, LU C L, et al. Nucleation in supercooled water triggered by mechanical impact: Experimental and theoretical analyses[J]. Journal of Energy Storage202252: 104755.
[14] GORDILLO L, SUN TP, CHENG X. Dynamics of drop impact on solid surfaces: evolution of impact force and self-similar spreading[J]. Journal of Fluid Mechanics2018840: 190-214.
[15] KANT P, KOLDEWEIJ R B J, HARTH K, et al. Fast-freezing kinetics inside a droplet impacting on a cold surface[J]. Proceedings of the National Academy of Sciences of the United States of America2020117(6): 2788-2794.
[16] HU M, WANG F, TAO Q, et al. Frozen patterns of impacted droplets: From conical tips to toroidal shapes[J]. Physical Review Fluids20205(8): 081601.
[17] THIéVENAZ V, JOSSERAND C, SéON T. Retraction and freezing of a water film on ice[J]. Physical Review Fluids20205(4): 041601.
[18] ANGELL C A. Supercooled water[J]. Annual Review of Physical Chemistry198334: 593-630.
[19] HOLTEN V, LIMMER D T, MOLINERO V, et al. Nature of the anomalies in the supercooled liquid state of the mW model of water[J]. The Journal of Chemical Physics2013138(17): 174501.
[20] KALIKMANOV V I. Nucleation Theory[M]. Dordrecht: Springer Netherlands, 2013.
[21] BAI G Y, GAO D, LIU Z, et al. Probing the critical nucleus size for ice formation with graphene oxide nanosheets[J]. Nature2019576(7787): 437-441.
[22] PEREYRA R G, SZLEIFER I, CARIGNANO M A. Temperature dependence of ice critical nucleus size[J]. The Journal of Chemical Physics2011135(3): 034508.
[23] PEREYRA R G, SEBASTIANELLI P, áVILA E E. Homogeneous nucleation in supercooled liquid water. Determination of ice germ size and activation energy barrier in Molecular Dynamics simulations[J]. Molecular Simulation202248(12): 1112-1121.
[24] LIN M, XIONG Z W, CAO H S. Bridging classical nucleation theory and molecular dynamics simulation for homogeneous ice nucleation[J]. The Journal of Chemical Physics2024161(8): 084504.
[25] LUO S, WANG J, LI Z G. Homogeneous ice nucleation under shear[J]. The Journal of Physical Chemistry B2020124(18): 3701-3708.
[26] ESPINOSA J R, SORIA G D, RAMIREZ J, et al. Role of salt, pressure, and water activity on homogeneous ice nucleation[J]. The Journal of Physical Chemistry Letters20178(18): 4486-4491.
[27] SORIA G D, ESPINOSA J R, RAMIREZ J, et al. A simulation study of homogeneous ice nucleation in supercooled salty water[J]. The Journal of Chemical Physics2018148(22): 222811.
[28] LUO S, LI C, LI F, et al. Ice crystallization in shear flows[J]. The Journal of Physical Chemistry C2019123(34): 21042-21049.
[29] RICHARD D, SPECK T. The role of shear in crystallization kinetics: From suppression to enhancement[J]. Scientific Reports20155: 14610.
[30] ICKES L, WELTI A, HOOSE C, et al. Classical nucleation theory of homogeneous freezing of water: thermodynamic and kinetic parameters[J]. Physical Chemistry Chemical Physics201517(8): 5514-5537.
[31] THOMPSON A P, AKTULGA H M, BERGER R, et al. LAMMPS-A flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales[J]. Computer Physics Communications2022271: 108171.
[32] FARMER T O, MARKVARDSEN A J, ROD T H, et al. Dynamical accuracy of water models on supercooling[J]. The Journal of Physical Chemistry Letters202011(18): 7469-7475.
[33] DICK T J, MADURA J D. Annual reports in computational chemistry[M]. Pittsburgh: Elsevier, 2005: 59-74.
[34] CHAN H, CHERUKARA M J, NARAYANAN B, et al. Machine learning coarse grained models for water[J]. Nature Communications201910: 379.
[35] HUJO W, SHADRACK JABES B, RANA V K, et al. The rise and fall of anomalies in tetrahedral liquids[J]. Journal of Statistical Physics2011145(2): 293-312.
[36] MOLINERO V, MOORE E B. Water modeled as an intermediate element between carbon and silicon[J]. The Journal of Physical Chemistry B2009113(13): 4008-4016.
[37] NGUYEN A H, MOLINERO V. Identification of clathrate hydrates, hexagonal ice, cubic ice, and liquid water in simulations: The CHILL+ algorithm[J]. The Journal of Physical Chemistry B2015119(29): 9369-9376.
[38] MOORE E B, MOLINERO V. Ice crystallization in water’s “no-man’s land”[J]. The Journal of Chemical Physics2010132(24): 244504.
[39] ALVAREZ F, ARBE A, COLMENERO J. Unraveling the coherent dynamic structure factor of liquid water at the mesoscale by molecular dynamics simulations[J]. The Journal of Chemical Physics2021155(24): 244509.
[40] SOUZA J B Jr, SCHLEDER G R, BETTINI J, et al. Pair distribution function obtained from electron diffraction: an advanced real-space structural characterization tool[J]. Matter20214(2): 441-460.
[41] MARTYNA G J, TOBIAS D J, KLEIN M L. Constant pressure molecular dynamics algorithms[J]. The Journal of Chemical Physics1994101(5): 4177-4189.
[42] SOPER A K. Joint structure refinement of X-ray and neutron diffraction data on disordered materials[J]. Journal of Physics: Condensed Matter200719(33): 335206.
[43] ANGELL C A. Amorphous water[J]. Annual Review of Physical Chemistry200455: 559-583.
[44] MOORE E B, MOLINER O V. Structural transformation in supercooled water controls the crystallization rate of ice[J]. Nature2011479(7374): 506-508.
[45] NARTEN A H, VENKATESH C G, RICE S A. Diffraction pattern and structure of amorphous solid water at 10 and 77?°K[J]. The Journal of Chemical Physics197664(3): 1106-1121.
[46] MALKIN T L, MURRAY B J, SALZMANN C G, et al. Stacking disorder in ice I[J]. Physical Chemistry Chemical Physics201517(1): 60-76.
[47] 王瑞, 任瑛, 陈卫, 等. 冰水界面动态结构的分子动力学模拟研究[J]. 化工学报202273(3): 1315-1323.
  WANG R, REN Y, CHEN W, et al. Molecular dynamics simulation on the dynamic structure of icing interface[J]. CIESC Journal202273(3): 1315-1323 (in Chinese).
[48] TANAKA H. Possible resolution of the Kauzmann paradox in supercooled liquids[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics200368(1 Pt 1): 011505.
[49] ESPINOSA J R, SANZ E, VALERIANI C, et al. Homogeneous ice nucleation evaluated for several water models[J]. The Journal of Chemical Physics2014141(18): 18C529.
[50] LI M D, HUANG Y P, XIA Y J, et al. Effective nucleation size for ice crystallization[J]. Journal of Chemical Theory and Computation202521(4): 1990-1996.
[51] European Union Aviation Safety Agency. Easy access rules for large aeroplanes: CS-25 [S]. Cologne: European Union Aviation Safety Agency, 2023:1341-1351.
[52] SHAH P, DRISCOLL M M. Drop impact dynamics of complex fluids: A review[J]. Soft Matter202420(25): 4839-4858.
[53] GONZALEZ-AVILA S R, ZENG Q Y, OHL C D. Pressure and wall shear stress from high-speed droplet impact[J]. International Journal of Multiphase Flow2024181: 104981.
[54] KE Q, GONG X T, LIAO S W, et al. Effects of thermostats/barostats on physical properties of liquids by molecular dynamics simulations[J]. Journal of Molecular Liquids2022365: 120116.
[55] CHEN M B, SONG W J, LIN W Y, et al. Ice nucleation in supercooled water under shear[J]. Chemical Engineering Science2024300: 120674.
[56] GOSWAMI A, DALAL I S, SINGH J K. Universal nucleation behavior of sheared systems[J]. Physical Review Letters2021126(19): 195702.
[57] ESPINOSA J R, ZARAGOZA A, ROSALES-PELAEZ P, et al. Interfacial free energy as the key to the pressure-induced deceleration of ice nucleation[J]. Physical Review Letters2016117(13): 135702.
[58] HARRIS K R, NEWITT P J. Self-diffusion of water at low temperatures and high pressure[J]. Journal of Chemical & Engineering Data199742(2): 346-348.
文章导航

/