搜索
(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号(43)申请公布日(21)申请号201910046416.X(22)申请日2019.01.18(71)申请人南京航空航天大学地址211106江苏省南京市秦淮区御道街29号(72)发明人陶杰 许杨江山 沈一洲 江家威 卢阳 (74)专利代理机构北京恒创益佳知识产权代理事务所(普通合伙)11556代理人付金豹(51)Int.Cl.G06F17/50(2006.01)(54)发明名称一种基于分子动力学模拟微观结构表面防结冰性能的评价方法(57)摘要本发明公开了一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，该方法借助分子动力学理论和仿真计算软件，建立不同形状微结构表面与水接触的异质形核模型，避开计算量需求极大的形核率计算，对水的结冰温度，结冰时间进行仿真计算，基于仿真结果评价微结构形状对于水异质形核结冰难易程度的影响，从而从分子层面解释结冰形核机理，减小实际实验试错法的消耗，可为各种期望利用表面微观结构改性达到防结冰效果的实验探索，工艺优化与产品开发提供理论支持。权利要求书2页说明书7页附图2页CN109711094A2019.05.03CN109711094A1.一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，其特征在于，包括以下步骤：1)基于分子动力学仿真软件LAMMPS建立基体表面不同形状微结构与水的接触模型；采用三维周期盒子，首先在XY轴上建立无限长、厚度一定的平面基体，然后在平面基体上构建规则阵列微结构，控制微结构上下面的形状与尺寸得到特定形状的微结构，同时定义微结构高度与间距；微结构创建后，固定基体与微结构的原子坐标；在微结构上生成水分子，保证水分子与基体距离保持在截断半径内，并进行能量最小化设置；2)基于分子动力学仿真软件LAMMPS对步骤1)中的模型进行初始化准备；具体步骤为：首先赋予原子质量与粒子初始速度，采用SW粗粒化势函数表征水与基体原子之间的相互作用，将系统处于室温，一个大气压下进行数纳秒平衡弛豫，系综选择为NPT系综，积分步长可自行选定，使系统达到平衡，密度合理；3)基于分子动力学仿真软件LAMMPS对步骤2)得到的平衡态模型进行降温动力学弛豫；在一个大气压下，设定温度从298k均匀变化到200k或以下温度，降温速度可因精确度而定，每隔固定步长设定一个状态点，利用dump命令输出各个状态点的水分子坐标；4)基于分子动力学仿真软件LAMMPS对步骤2)得到的平衡态模型进行低温平衡态动力学弛豫；在一个大气压下，设定保温温度，进行平衡态动力学弛豫；每隔固定步长设定一个状态点，利用dump命令输出各个状态点的水分子坐标；5)利用脚本程序读取降温动力学过程得出的数据结果，判定各个分子的结冰情况，并输出每个状态点的结冰分子数量；得到数据后利用origin软件作图，得到结冰粒子数与温度变化之间的关系图，取结冰分子数量开始明显上升的温度作为结冰温度；6)读取低温平衡态动力学弛豫得出的数据结果，利用Q6判据判定各个分子的结冰情况，并输出每个状态点的结冰分子数量；得到数据后利用origin软件作图，得到结冰分子数与时间变化之间的关系图，取结冰分子数量开始明显上升的温度作为结冰时间；7)利用步骤5)、6)得到的结冰温度与结冰时间数据，评价微结构形状抑制结冰性能的优劣。2.根据权利要求1所述的评价方法，其特征在于：在步骤1)中，建立不同微结构与水接触的模型：微结构为周期性阵列结构，控制微观结构形状与尺寸，微观结构高度和单个微观结构之间的距离根据需求定义，高度定为1nm-4nm，间距定为1nm-4nm；微结构的形状是圆台，立方体，棱台之一。3.根据权利要求1所述的评价方法，其特征在于：在步骤2)中，参数初始化设置，基体元素种类可以改变，SW势函数适用于描述粗粒化水分子与基体的相互作用，具体表现形式为：权 利 要 求 书1/2页2CN109711094A2式中ε为能量参数，表征原子之间的相互作用强度；σ为距离参数，表征截断半径；θ为水中三个原子所成的角度；r表示两原子之间的距离；其余参数为无量纲的常数。其中Φ2描述二体势能量，Φ3描述三体势能量，系统总能量E等于二体势能量与三体势能量的加和，对于微结构与水接触的异质形核问题，水分子之间的相互作用由三体势描述，水与基体原子间的相互作用由二体势描述，由于基体原子固定，则不用设定基体原子间的作用势，以减少计算量；由于基体原子不同，水分子与基体原子的势函数值也会有所不同，具体可通过赋予不同的ε值改变基体的元素种类。4.根据权利要求1所述的评价方法，其特征在于：在步骤3)中，降温速度，降温范围与输出时间步长均可自行设定。降温速度一般在0.1k/ns-1k/ns之间，降温范围从室温到200k附近任意值。5.根据权利要求1所述的评价方法，其特征在于：在步骤4)中，设定的平衡温度可以设定为冰点以下的任意温度，优选选择230k，保温时间需要随保温温度改变，温度较高时，保温时间需要延长，温度较低时保温时间缩短。6.根据权利要求1所述的评价方法，其特征在于：在步骤5)中，结冰温度为结冰分子数目结束震荡，开始明显上升的温度。7.根据权利要求1所述的评价方法，其特征在于：在步骤6)中，结冰时间为结冰分子数目结束震荡，开始明显上升的时间。8.根据权利要求1所述的评价方法，其特征在于：在步骤7)中，结冰温度越低，结冰时间越长，微结构对于结冰的抑制能力越优。权 利 要 求 书2/2页3CN109711094A3一种基于分子动力学模拟微观结构表面防结冰性能的评价方法技术领域[0001]本发明属于材料表面性能评价技术领域，涉及一种基于分子动力学模拟微观结构表面防结冰性能的评价方法。本发明所提供的分子动力学模拟评定方法意在对高效防结冰表面的微观结构设计提供理论指导。背景技术[0002]结冰是一种常见的自然现象，但是在某些特殊工业条件下，结冰会带来严重的安全危害，例如发动机叶片，机械轴承，传输电线上结冰会导致十分严重的事故，如何避免特殊工况下的结冰现象是当今颇为重要的研究课题。[0003]材料表面微观结构对于冷凝结冰过程有着非常重要的影响，主要体现在两方面：一方面，水在基体表面结冰需要先润湿基体表面，表面特殊微观结构能够将空气“封锁”在结构内部，这些空气能够有效地隔离开水和基体，使得固液接触面积减小，接触角增大，影响水在基体上的润湿过程，从而抑制结冰形核的发生。另一方面，微观结构自身的角度与曲率等几何因素也会影响水在微观结构上的形核结冰过程。较大的曲率，特殊的角度等因素都会抑制水分子的形核结冰，判别并设计有效抑制形核结冰的微观结构对于防结冰技术领域的探索至关重要。[0004]现有实验方法可以通过试错法来构建大量不同的微观结构，逐一检测它们的防结冰性能，从而得到一些规律。但试错法需要大量的时间成本和材料成本，并且无法穷尽各种微观结构形貌。而且实验无法捕捉到结冰形核发生时的相关情况，实验观察到结冰现象时，形核结冰其实早开始。因此，需要寻找新的理论研究方法。[0005]分子动力学方法由于其是从分子层面观察分析问题，更加接近事物本质，易于捕捉事物的微小变化，而且不需要耗费实际的资源，逐渐成为一种重要的研究方法。由于形核是一个随机过程，虽然分子动力学技术可以通过统计学方法计算形核率，但需要耗费大量的时间与计算资源。发明内容[0006]本发明的目的在于提供一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，以期对防结冰表面微观结构的设计提供理论指导。[0007]技术方案：[0008]一种基于分子动力学模拟微观结构表面防结冰性能的评价方法，包括以下步骤：[0009]1)基于分子动力学仿真软件LAMMPS建立基体表面不同形状微结构与水的接触模型；采用三维周期盒子，首先在XY轴上建立无限长、厚度一定的平面基体，然后在平面基体上构建规则阵列微结构，控制微结构上下面的形状与尺寸得到特定形状的微结构，同时定义微结构高度与间距；微结构创建后，固定基体与微结构的原子坐标；在微结构上生成水分子，保证水分子与基体距离保持在截断半径内，并进行能量最小化设置；说 明 书1/7页4CN109711094A4[0010]2)基于分子动力学仿真软件LAMMPS对步骤1)中的模型进行初始化准备；具体步骤为：首先赋予原子质量与粒子初始速度，采用SW粗粒化势函数表征水与基体原子之间的相互作用，将系统处于室温，一个大气压下进行数纳秒平衡弛豫，系综选择为NPT系综，积分步长可自行选定，使系统达到平衡，密度合理；[0011]3)基于分子动力学仿真软件LAMMPS对步骤2)得到的平衡态模型进行降温动力学弛豫；在一个大气压下，设定温度从298k均匀变化到200k或以下温度，降温速度可因精确度而定，每隔固定步长设定一个状态点，利用dump命令输出各个状态点的水分子坐标；[0012]4)基于分子动力学仿真软件LAMMPS对步骤2)得到的平衡态模型进行低温平衡态动力学弛豫；在一个大气压下，设定保温温度，进行平衡态动力学弛豫；每隔固定步长设定一个状态点，利用dump命令输出各个状态点的水分子坐标；[0013]5)利用脚本程序读取降温动力学过程得出的数据结果，判定各个分子的结冰情况，并输出每个状态点的结冰分子数量；得到数据后利用origin软件作图，得到结冰粒子数与温度变化之间的关系图，取结冰分子数量开始明显上升的温度作为结冰温度；[0014]6)读取低温平衡态动力学弛豫得出的数据结果，利用Q6判据判定各个分子的结冰情况，并输出每个状态点的结冰分子数量；得到数据后利用origin软件作图，得到结冰分子数与时间变化之间的关系图，取结冰分子数量开始明显上升的温度作为结冰时间；[0015]7)利用步骤5)、6)得到的结冰温度与结冰时间数据，评价微结构形状抑制结冰性能的优劣。[0016]所述的评价方法，在步骤1)中，建立不同微结构与水接触的模型：微结构为周期性阵列结构，控制微观结构形状与尺寸，微观结构高度和单个微观结构之间的距离根据需求定义，高度定为1nm-4nm，间距定为1nm-4nm；微结构的形状是圆台，立方体，棱台之一。[0017]所述的评价方法，在步骤2)中，参数初始化设置，基体元素种类可以改变，SW势函数适用于描述粗粒化水分子与基体的相互作用，具体表现形式为：[0018][0019][0020][0021]式中ε为能量参数，表征原子之间的相互作用强度；σ为距离参数，表征截断半径；θ为水中三个原子所成的角度；r表示两原子之间的距离；其余参数为无量纲的常数。其中Φ2描述二体势能量，Φ3描述三体势能量，系统总能量E等于二体势能量与三体势能量的加和，对于微结构与水接触的异质形核问题，水分子之间的相互作用由三体势描述，水与基体原子间的相互作用由二体势描述，由于基体原子固定，则不用设定基体原子间的作用势，以减少计算量；由于基体原子不同，水分子与基体原子的势函数值也会有所不同，具体可通过赋予不同的ε值改变基体的元素种类。[0022]所述的评价方法，在步骤3)中，降温速度，降温范围与输出时间步长均可自行设说 明 书2/7页5CN109711094A5定。降温速度一般在0.1k/ns-1k/ns之间，降温范围从室温到200k附近任意值。[0023]所述的评价方法，在步骤4)中，设定的平衡温度可以设定为冰点以下的任意温度，优选选择230k，保温时间需要随保温温度改变，温度较高时，保温时间需要延长，温度较低时保温时间缩短。[0024]所述的评价方法，在步骤5)中，结冰温度为结冰分子数目结束震荡，开始明显上升的温度。[0025]所述的评价方法，在步骤6)中，结冰时间为结冰分子数目结束震荡，开始明显上升的时间。[0026]所述的评价方法，在步骤7)中，结冰温度越低，结冰时间越长，微结构对于结冰的抑制能力越优。[0027]采用本发明的材料微观结构表面防结冰性能的评价方法具有以下优势：[0028]排除了实际实验时外界不稳定因素的影响，可直接判断微观结构几何参数这一自变量对于结冰难易程度的影响，避免了试错法的资