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

doi:10.7527/S1000-6893.2025.32342

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

本期目录 | 过刊浏览 | 高级检索

前一篇 |

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

汉芮岐¹, 徐绯¹^,²(), 陈成¹, 任战鹏³, 张惠³

^1.西北工业大学航空学院，西安 710072
^2.强度与结构完整性全国重点实验室，西安 710072 3. 中国飞机强度研究所，西安 710065

收稿日期:2025-05-30 修回日期:2025-07-30 接受日期:2025-08-18 出版日期:2025-09-08 发布日期:2025-08-28
通讯作者: 徐绯 E-mail:xufei@nwpu.edu.cn
基金资助:
国家自然科学基金(D5110220438)

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

Ruiqi HAN¹, Fei XU¹^,²(), Cheng CHEN¹, Zhanpeng REN³, Hui ZHANG³

^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

摘要：

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

关键词: 过冷水, 冰, 温度, 应力, 分子动力学模拟

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 10⁷ - 10⁸ s^-1， and the triaxial compression strain rate is 10⁵ - 10⁶ 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

中图分类号:

V219

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

Ruiqi HAN, Fei XU, Cheng CHEN, Zhanpeng REN, Hui ZHANG. Molecular dynamics simulation of stress and temperature effects on supercooled water freezing[J]. Acta Aeronautica et Astronautica Sinica, 2026, 47(4): 232342.

图/表 17

图 1

图 2

图 3

图 4

图 5

图 6

图 7

图 8

图 9

图 10

图 11

图 12

图 13

图 14

图 15

图 16

图 17

参考文献 58

[1]	WU Z L. Drop “impact” on an airfoil surface［J］. Advances in Colloid and Interface Science， 2018， 256： 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 Mathematics， 2023， 429： 115230.
[3]	CHENG X， SUN T-P， GORDILLO L. Drop impact dynamics： impact force and stress distributions［J］. Annual Review of Fluid Mechanics， 2022， 54： 57-81.
[4]	SUN T P， ÁLVAREZ-NOVOA F， ANDRADE K， et al. Stress distribution and surface shock wave of drop impact［J］. Nature Communications， 2022， 13： 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 Mechanics， 2016， 795： 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 Fluids， 2016， 28（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 Refrigeration， 2011， 34（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 Transfer， 2016， 97： 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 Science， 2021， 584： 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 Transfer， 2019， 130： 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 Transfer， 2019， 142： 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 Storage， 2022， 52： 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 Mechanics， 2018， 840： 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 America， 2020， 117（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 Fluids， 2020， 5（8）： 081601.
[17]	THIÉVENAZ V， JOSSERAND C， SÉON T. Retraction and freezing of a water film on ice［J］. Physical Review Fluids， 2020， 5（4）： 041601.
[18]	ANGELL C A. Supercooled water［J］. Annual Review of Physical Chemistry， 1983， 34： 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 Physics， 2013， 138（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］. Nature， 2019， 576（7787）： 437-441.
[22]	PEREYRA R G， SZLEIFER I， CARIGNANO M A. Temperature dependence of ice critical nucleus size［J］. The Journal of Chemical Physics， 2011， 135（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 Simulation， 2022， 48（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 Physics， 2024， 161（8）： 084504.
[25]	LUO S， WANG J， LI Z G. Homogeneous ice nucleation under shear［J］. The Journal of Physical Chemistry B， 2020， 124（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 Letters， 2017， 8（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 Physics， 2018， 148（22）： 222811.
[28]	LUO S， LI C， LI F， et al. Ice crystallization in shear flows［J］. The Journal of Physical Chemistry C， 2019， 123（34）： 21042-21049.
[29]	RICHARD D， SPECK T. The role of shear in crystallization kinetics： From suppression to enhancement［J］. Scientific Reports， 2015， 5： 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 Physics， 2015， 17（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 Communications， 2022， 271： 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 Letters， 2020， 11（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 Communications， 2019， 10： 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 Physics， 2011， 145（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 B， 2009， 113（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 B， 2015， 119（29）： 9369-9376.
[38]	MOORE E B， MOLINERO V. Ice crystallization in water’s “no-man’s land”［J］. The Journal of Chemical Physics， 2010， 132（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 Physics， 2021， 155（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］. Matter， 2021， 4（2）： 441-460.
[41]	MARTYNA G J， TOBIAS D J， KLEIN M L. Constant pressure molecular dynamics algorithms［J］. The Journal of Chemical Physics， 1994， 101（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 Matter， 2007， 19（33）： 335206.
[43]	ANGELL C A. Amorphous water［J］. Annual Review of Physical Chemistry， 2004， 55： 559-583.
[44]	MOORE E B， MOLINER O V. Structural transformation in supercooled water controls the crystallization rate of ice［J］. Nature， 2011， 479（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 Physics， 1976， 64（3）： 1106-1121.
[46]	MALKIN T L， MURRAY B J， SALZMANN C G， et al. Stacking disorder in ice I［J］. Physical Chemistry Chemical Physics， 2015， 17（1）： 60-76.
[47]	王瑞，任瑛，陈卫，等. 冰水界面动态结构的分子动力学模拟研究［J］. 化工学报， 2022， 73（3）： 1315-1323.
	WANG R， REN Y， CHEN W， et al. Molecular dynamics simulation on the dynamic structure of icing interface［J］. CIESC Journal， 2022， 73（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 Physics， 2003， 68（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 Physics， 2014， 141（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 Computation， 2025， 21（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 Matter， 2024， 20（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 Flow， 2024， 181： 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 Liquids， 2022， 365： 120116.
[55]	CHEN M B， SONG W J， LIN W Y， et al. Ice nucleation in supercooled water under shear［J］. Chemical Engineering Science， 2024， 300： 120674.
[56]	GOSWAMI A， DALAL I S， SINGH J K. Universal nucleation behavior of sheared systems［J］. Physical Review Letters， 2021， 126（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 Letters， 2016， 117（13）： 135702.
[58]	HARRIS K R， NEWITT P J. Self-diffusion of water at low temperatures and high pressure［J］. Journal of Chemical & Engineering Data， 1997， 42（2）： 346-348.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

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

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

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 17

参考文献 58

相关文章 15

编辑推荐

Metrics

本文评价

[1]	齐海峰, 常士楠, 任战鹏, 杨映麟. 基于沿程冷却法的冰晶沿程参数变化模拟[J]. 航空学报, 2026, 47(4): 132227-132227.
[2]	范春浩, 胡俊山, 杨智勇, 侯富森, 陈培林, 田威. 固化残余应力和变形对钛合金损伤构件CFRP单面贴补结构拉伸性能的影响[J]. 航空学报, 2026, 47(3): 431997-431997.
[3]	唐海洋, 任哲钒, 陈宁立, 易贤, 陈志勇. 基于PLE网络的电加热防冰功率分布优化设计[J]. 航空学报, 2026, 47(2): 132179-132179.
[4]	于立, 李文豪, 张健, 张卓然, 孙窈, 胡光源. 航空油冷双凸极高压直流发电机绕组故障下的温升特性[J]. 航空学报, 2026, 47(1): 632094-632094.
[5]	刘欣乐, 姜亚楠, 辛荣提, 李庆辉, 蔡晋生. 超疏水电热复合分区防冰策略[J]. 航空学报, 2025, 46(9): 130784-130784.
[6]	于宝石, 雷勇军, 申志彬, 张大鹏. 固体发动机药柱结构固化残余应力分析与控制技术[J]. 航空学报, 2025, 46(8): 31083-031083.
[7]	刘峰, 杨森, 韦振鹏. 基于复合材料弹性周期结构和预应力蒙皮的变弦长机翼[J]. 航空学报, 2025, 46(7): 230966-230966.
[8]	冯蕴雯, 腾达, 路成, 王锐, 陈俊宇. 机理与生成对抗代理建模融合的飞机系统可靠性评估[J]. 航空学报, 2025, 46(7): 230948-230948.
[9]	易贤, 任靖豪, 赖庆仁, 刘宇, 王强. 大型飞机全尺寸多段翼结冰特性计算和试验[J]. 航空学报, 2025, 46(5): 531575-531575.
[10]	高天翔, 罗振兵, 周岩, 彭文强, 程盼. 合成双射流控制水滴轨迹特性实验[J]. 航空学报, 2025, 46(4): 130833-130833.
[11]	杨娟, 胡佳宁, 佟佳成, 张青松. 航空锂电池热失控高温喷射冲击实验[J]. 航空学报, 2025, 46(23): 430965-430965.
[12]	李秋华. 基于格林函数的热障涂层结构热力耦合场分析[J]. 航空学报, 2025, 46(22): 431764-431764.
[13]	雷凯, 吴敬涛, 邓文亮. 热力耦合作用下复合材料-金属连接结构力学性能分析[J]. 航空学报, 2025, 46(21): 532276-532276.
[14]	白春玉, 赵子豪, 刘小川, 张宇, 郭亚周. 飞机结构抗冰雹冲击特性研究进展与展望[J]. 航空学报, 2025, 46(21): 532424-532424.
[15]	汤林, 孟宣市, 刘蕊迪, 顾兴士, 易贤. 飞机翼面冰的结构建模与弹性模量预测[J]. 航空学报, 2025, 46(20): 531895-531895.