[发明专利]基于分子动力学模拟微观结构表面防结冰性能的评价方法

说明书

本发明公开了一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，该方法借助分子动力学理论和仿真计算软件，建立不同形状微结构表面与水接触的异质形核模型，避开计算量需求极大的形核率计算，对水的结冰温度，结冰时间进行仿真计算，基于仿真结果评价微结构形状对于水异质形核结冰难易程度的影响，从而从分子层面解释结冰形核机理，减小实际实验试错法的消耗，可为各种期望利用表面微观结构改性达到防结冰效果的实验探索，工艺优化与产品开发提供理论支持。

技术领域

本发明属于材料表面性能评价技术领域，涉及一种基于分子动力学模拟微观结构表面防结冰性能的评价方法。本发明所提供的分子动力学模拟评定方法意在对高效防结冰表面的微观结构设计提供理论指导。

背景技术

结冰是一种常见的自然现象，但是在某些特殊工业条件下，结冰会带来严重的安全危害，例如发动机叶片，机械轴承，传输电线上结冰会导致十分严重的事故，如何避免特殊工况下的结冰现象是当今颇为重要的研究课题。

材料表面微观结构对于冷凝结冰过程有着非常重要的影响，主要体现在两方面：一方面，水在基体表面结冰需要先润湿基体表面，表面特殊微观结构能够将空气“封锁”在结构内部，这些空气能够有效地隔离开水和基体，使得固液接触面积减小，接触角增大，影响水在基体上的润湿过程，从而抑制结冰形核的发生。另一方面，微观结构自身的角度与曲率等几何因素也会影响水在微观结构上的形核结冰过程。较大的曲率，特殊的角度等因素都会抑制水分子的形核结冰，判别并设计有效抑制形核结冰的微观结构对于防结冰技术领域的探索至关重要。

现有实验方法可以通过试错法来构建大量不同的微观结构，逐一检测它们的防结冰性能，从而得到一些规律。但试错法需要大量的时间成本和材料成本，并且无法穷尽各种微观结构形貌。而且实验无法捕捉到结冰形核发生时的相关情况，实验观察到结冰现象时，形核结冰其实早开始。因此，需要寻找新的理论研究方法。

分子动力学方法由于其是从分子层面观察分析问题，更加接近事物本质，易于捕捉事物的微小变化，而且不需要耗费实际的资源，逐渐成为一种重要的研究方法。由于形核是一个随机过程，虽然分子动力学技术可以通过统计学方法计算形核率，但需要耗费大量的时间与计算资源。

发明内容

本发明的目的在于提供一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，以期对防结冰表面微观结构的设计提供理论指导。

技术方案：

一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，包括以下步骤：

1)基于分子动力学仿真软件LAMMPS建立基体表面不同形状微结构与水的接触模型；采用三维周期盒子，首先在XY轴上建立无限长、厚度一定的平面基体，然后在平面基体上构建规则阵列微结构，控制微结构上下面的形状与尺寸得到特定形状的微结构，同时定义微结构高度与间距；微结构创建后，固定基体与微结构的原子坐标；在微结构上生成水分子，保证水分子与基体距离保持在截断半径内，并进行能量最小化设置；

2)基于分子动力学仿真软件LAMMPS对步骤1)中的模型进行初始化准备；具体步骤为：首先赋予原子质量与粒子初始速度，采用SW粗粒化势函数表征水与基体原子之间的相互作用，将系统处于室温，一个大气压下进行数纳秒平衡弛豫，系综选择为NPT系综，积分步长可自行选定，使系统达到平衡，密度合理；

3)基于分子动力学仿真软件LAMMPS对步骤2)得到的平衡态模型进行降温动力学弛豫；在一个大气压下，设定温度从298k均匀变化到200k或以下温度，降温速度可因精确度而定，每隔固定步长设定一个状态点，利用dump命令输出各个状态点的水分子坐标；

4)基于分子动力学仿真软件LAMMPS对步骤2)得到的平衡态模型进行低温平衡态动力学弛豫；在一个大气压下，设定保温温度，进行平衡态动力学弛豫；每隔固定步长设定一个状态点，利用dump命令输出各个状态点的水分子坐标；