基于分子动力学模拟的纳米液滴接触角获取方法技术

本发明专利技术适用于测量技术领域，提供了一种基于分子动力学模拟的纳米液滴接触角获取方法，包括：利用LAMMPS软件输出液滴分子平衡时的位置和数密度，导入TECPLOT软件；TECPLOT软件将导入的数密度分成不同的数密度区域，将数密度中值所在的数密度等值线作为液气表面和固气表面，提取位于所述数密度等值线上的数据点，数据点由湿润长度及液滴高度组成，导入MATLAB；利用MATLAB多项式拟合法对导入的数据点进行液气拟合曲线及固气拟合曲线拟合，并确定两拟合曲线的交点；确定交点在两拟合曲线上的倾斜角，两倾斜角差值的绝对值即为液滴的静态接触角。通过捕捉合适的液气和固气界面使得计算液滴接触角更加可信，解决了人为测量静态接触角存在的测量误差大，数据可信度低的问题。

A method for obtaining contact angle of nano droplets based on molecular dynamics simulation

【技术实现步骤摘要】
基于分子动力学模拟的纳米液滴接触角获取方法
本专利技术属于测量
，提供了一种基于分子动力学模拟的纳米液滴接触角获取方法。
技术介绍
液滴润湿涵盖多个尺度(宏观、介观以及微观)，是基础研究和工程应用的重点研究方向，例如农药沉积、飞机机身防冰、内燃机内燃油黏壁以及微纳米机器的黏附失效等问题。在农药喷雾中，农药液滴越亲水，则农药与植物叶面的接触面积越大，除虫效果越明显。而对于飞机机身防冰、内燃机内燃油黏壁以及微纳米机器的黏附失效，液滴若能在壁面静止时呈现疏水或者超疏水形貌，则液滴在撞击壁面过程中容易发生反弹并离开壁面，这将有利于工程需求。近年来，采用分子动力学模拟方法,利用Large-scaleAtomic/MolecularMassivelyParallelSimulator(LAMMPS，一个分子动力学的模拟软件)软件研究纳米液滴在固体壁面上的润湿特性、自组织行为以及固体表面流体的凝结换热问题受到大批学者的关注。液滴的静态润湿性采用其在固体表面所形成的接触角(该角是指气-固-液三相交界处，固液界面到液气表面所经过液滴内部的夹角)来区分。对于宏观液滴常采用接触角测量仪截取液滴平衡时的图片，确定三相接触线的液气界面，然后利用接触角测量仪内置软件选取切线，读出相应的角度即可。另一方面也有部分学者采用相场动力学数值方法捕捉相界面，然后利用界面拟合法或者界面切线法得到接触角。然而，随着液滴尺度的缩小，一方面受制于人眼生理的观察限制，另一方面得益于利用分子动力学研究润湿性的日益发展，对于纳米液滴数值仿真时纳米液滴静态接触角的获取还没很好的处理方法。大多数研究者常采用截取纳米液滴静态图片，利用软件量取角度的方法获取。此外，也有部分学者考虑液气界面存在一定的厚度，并且在接近固壁处液滴分子会形成“固溶区”，这两个因素将极大影响纳米液滴接触角的测量。综上所述，对纳米液滴静态接触角的测量既有人为因素影响又有其本身固有的难度，这导致现有方法测量误差较大，数据可信度较小。
技术实现思路
本专利技术实施例提供一种基于分子动力学模拟的纳米液滴接触角获取方法，旨在解决人为测量静态接触角存在的测量误差大，数据可信度低的问题。本专利技术是这样实现的，一种基于分子动力学模拟的纳米液滴接触角获取方法，所述方法具体包括如下步骤：S1、利用LAMMPS软件输出液滴分子平衡时的位置和数密度，并导入TECPLOT软件；S2、TECPLOT软件将导入的数密度分成不同的数密度区域，将数密度中值所在的数密度等值线作为液气表面和固气表面，提取位于所述数密度等值线上的数据点，数据点由湿润长度及液滴高度组成，将提取的数据点导入MATLAB；S3、利用MATLAB多项式拟合法对导入的数据点进行液气拟合曲线及固气拟合曲线拟合，并确定两拟合曲线的交点；S4、确定上述交点在液气拟合曲线及固气拟合曲线上的倾斜角，两倾斜角差值的绝对值即为液滴的静态接触角。进一步的，液滴分子在平衡时的位置和数密度获取方法具体如下：S11、基于LAPPMS软件模拟液滴在基底上润湿时，利用郎之万热浴法使液滴驰豫50ps；S12、之后每隔1ps提取一组液滴分子的位置和数密度数据，连续提取L组；S13、针对每组位置和数密度数据分别执行步骤S1至步骤S4，获取L个静态接触角，将L个静态接触角的平均值作为静态接触角进行输出。进一步的，液气拟合曲线及固气拟合曲线的获取方法具体如下：S31、确定多项式阶数s，4≤s≤6，固气拟合曲线及液气拟合曲线分别为P(x)和Q(x)具体表达式如下：S32、利用最小二乘法拟合多项式系数向量a＝[a0,a1,...,am]和b＝[b0,b1,...,bn]，其中，n及m取值为s。本专利技术提供了一种的可信的液滴静态接触角的测量方法，同时通过捕捉合适的液气和固气界面使得计算液滴接触角更加可信，解决了人为测量静态接触角存在的测量误差大，数据可信度低的问题。附图说明图1为本专利技术实施例提供的分子动力学模拟的纳米液滴接触角获取方法流程图；图2为本专利技术实施例提供的TECPLOT处理后是数密度云图；图3为本专利技术实施例提供的MATLAB处理后的固气和液气拟合曲线图相交图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。图1为本专利技术实施例提供的分子动力学模拟的纳米液滴接触角获取方法流程图，该方法具体包括如下步骤：S1、利用LAMMPS软件中输出液滴分子在平衡时的位置和数密度数据，存入txt文本中，然后添加相应的文件头，导入进TECPLOT中；S2、TECPLOT将数密度数据划分不同的数密度区域，将数密度中值所在的数密度等值线作为为液气表面和固气界面，提取位于所述数密度等值线上的数据点(xi,yi)，xi为湿润长度，yi为液滴高度组成，如图2所示，将提取的数据点存为txt文本后导入到MATLAB中；S3、基于MATLAB多项式拟合法对导入的数据点进行液气拟合曲线和固气拟合曲线的拟合，确定两拟合曲线的交点；S4、定上述交点在液气拟合曲线及固气拟合曲线上的倾斜角，两倾斜角差值的绝对值即为液滴的静态接触角。在本专利技术实施例中，步骤1中输出相应的液滴分子在平衡时的位置和数密度数据，具体操作包含以下操作：S11、使用LAPPMS软件模拟液滴在基底上润湿时，利用郎之万热浴法使液滴驰豫50ps。然后，在in文件后添加以下命令：computecc1watchunk/atombin/cylinderzlower2000.070.035discardyesunitsbox；fix6watave/chunk100011000cc1density/numberaveonefilevel.profile；computecc1表示对计算盒子进行分割取片，计算数密度；fix6表示输出坐标信息和数密度信息；S12、之后每隔1ps提取一组液滴分子的位置和数密度数据，连续提取L组；L组位置和数密度数据作为可信数据，导入TECPLOT中进行数据的处理，有利于降低随机性，取平均后增加数据可信度；S13、针对每组位置和数密度数据分别执行步骤S1至步骤S4，获取L个静态接触角，将L个静态接触角的平均值作为静态接触角进行输出。在本专利技术实施例中，步骤2中txt文本中添加相应的文件头具体操作包含以下操作：S22、TECPLOT具有官方的数据文件格式，导入外部数据时数据源文件必须按照相应的格式填写，否则TECPLOT提示出错。具体文件头表达如下：Title＝"datafornumberdensity",表示题目，一般用作提示作用；VARIABLES＝"X","Y","P"表示变量名称，X为横坐标(液滴湿润长度)，Y为纵坐标(液滴高度)，P为数密度本文档来自技高网...

【技术保护点】
1.一种基于分子动力学模拟的纳米液滴接触角获取方法，其特征在于，所述方法具体包括如下步骤：/nS1、利用LAMMPS软件输出液滴分子平衡时的位置和数密度，并导入TECPLOT软件；/nS2、TECPLOT软件将导入的数密度分成不同的数密度区域，将数密度中值所在的数密度等值线作为液气表面和固气表面，提取位于所述数密度等值线上的数据点，数据点由湿润长度及液滴高度组成，将提取的数据点导入MATLAB；/nS3、利用MATLAB多项式拟合法对导入的数据点进行液气拟合曲线及固气拟合曲线拟合，并确定两拟合曲线的交点；/nS4、确定上述交点在液气拟合曲线及固气拟合曲线上的倾斜角，两倾斜角差值的绝对值即为液滴的静态接触角。/n

【技术特征摘要】
1.一种基于分子动力学模拟的纳米液滴接触角获取方法，其特征在于，所述方法具体包括如下步骤：
S1、利用LAMMPS软件输出液滴分子平衡时的位置和数密度，并导入TECPLOT软件；
S2、TECPLOT软件将导入的数密度分成不同的数密度区域，将数密度中值所在的数密度等值线作为液气表面和固气表面，提取位于所述数密度等值线上的数据点，数据点由湿润长度及液滴高度组成，将提取的数据点导入MATLAB；
S3、利用MATLAB多项式拟合法对导入的数据点进行液气拟合曲线及固气拟合曲线拟合，并确定两拟合曲线的交点；
S4、确定上述交点在液气拟合曲线及固气拟合曲线上的倾斜角，两倾斜角差值的绝对值即为液滴的静态接触角。


2.如权利要求1所述基于分子动力学模拟的纳米液滴接触角获取方法，其特征在于，液滴分子在平衡时的位置和...

【专利技术属性】
技术研发人员：尹宗军，王玲芝，许星月，洪雨，孙鹏，
申请(专利权)人：安徽信息工程学院，
类型：发明
国别省市：安徽;34