首页 / 专利库 / 物理 / 能量状态 / 一种水平规则双粗糙度表面重液滴接触角滞后的计算方法

一种平规则双粗糙度表面重液滴接触滞后的计算方法

阅读:478发布:2023-03-12

专利汇可以提供一种平规则双粗糙度表面重液滴接触滞后的计算方法专利检索,专利查询,专利分析的服务。并且一种 水 平规则双粗糙度表面重液滴 接触 角 滞后的计算方法,包括:液滴轮廓离散化;四种不同润湿接触状态的系统无量纲自由能函数的推导;计算出接触圆固定时的重液滴的润湿状态,接触圆半径和在一级(二级)机构中的渗透深度;设置初始体积增量,在接触圆半径固定的条件下计算对应的最小无量纲自由能;接触圆将变化时,计算当前的局部最小无量纲自由能;通过使用二进制 迭代 法,可以搜索到一个体积增量和表面接触角(前进角)使得当前的最小无量纲自由能和局部最小无量纲自由能只差等于 能量 势垒;设置体积初始减量。能够帮助研究人员进一步了解液滴的动态性能,对于 微 流体 系统 的设计具有重要指导意义。,下面是一种平规则双粗糙度表面重液滴接触滞后的计算方法专利的具体信息内容。

1.一种平规则双粗糙度表面重液滴接触滞后的计算方法,包括以下步骤:
(1)假设一个由液滴和微纳二级结构表面组成的模型,其中第一级纳米结构和第二级微米结构都是方柱,几何尺度周期性分布,二级结构和基底由不同的材料组成。模型的总表观面积为边长L0的正方形,边长远大于液滴直径。a1,b1,h1和x1分别是第一级结构方柱的边长、间距、高度和浸润深度,a2,b2,h2和x2分别是第二级结构方柱的边长、间距、高度和浸润深度,第一级结构的尺寸远小于第二级结构。由于水平规则的双粗糙度表面上的重液滴形状是轴对称的,所以在二维平面上考虑了这个问题。选择重液滴的对称轴为v轴,重液滴的表面接触面为u轴,并利用有限差分法将重液滴的轮廓曲线沿v方向分成N等份,得到离散点:
其中,ha是液滴的最大高度,ui是离散点的横坐标,vi是离散点的纵坐标。
(2)四种不同润湿接触状态的系统无量纲自由能函数的推导,推导过程如下:
21.通过重液滴在一级结构和二级结构中的渗透深度,将润湿接触状态分为四种,分别为W-W,C-W,W-C和C-C状态,划分依据如下:
(a)一级结构和二级结构均处于Wenzel状态,即x1=h1,x2=h2,这种接触状态命名为W-W状态。
(b)一级结构处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,二级结构处于Wenzel状态,即0≤x1(c)一级结构处于Wenzel状态,二级结构处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,即x1=h1,0≤x2(d)一级结构和二级结构均处于Cassie状态或者接触状态将从Cassie状态转变为
Wenzel状态,即0≤x1其中,x1是液滴在一级结构中的渗透深度,x2是液滴在二级结构中的渗透深度,h1是一级结构的高度,h2是二级结构的高度。
22.四种不同润湿接触状态的系统自由能EI等于EI=Ea+Eb,I=1,2,3,4。下标I=1,2,3,
4分别表示W-W,C-W,W-C和C-C状态。Ea和Eb分别代表重势能和表面自由能。
Ea可以表示为:
Eb可以表示为:
Eb=γLVSLV+γSL0SSL0+γSL1SSL1+γSL2SSL2+γSV0SSV0+γSV1SSV1+γSV2SSV2  (3)其中,ρ是液体密度,g是重力加速度,S是界面的接触面积,γ是界面张力系数,下标LV,SL和SV分别代表液-气,固-液和固-气界面,下标1,2和0分别代表一级结构,二级结构和基底。
液滴的外表面积Sext可以表示为:
式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度。
液滴与基底的表面接触面积可以表示为:
Sbase=πrb2  (5)
式中,rb为液滴与二级微纳结构表面接触圆半径。
根据每种浸润状态的不同,可以容易求得四种状态下的SLV,SSV0,SSV1,SSV2,SSL0,SSL1,SSL2,所对应的表达式:
(a)W-W状态:
SLV=Sext=2πR2(1-cosθ)  (6)
2
SSL0=Sbase(1-f1)(1-f2)=πrb(1-f1)(1-f2)  (7)
SSL1=Sbaser2[r1-(1-f1)]=πrb2r2(r1-1+f1)  (8)
SSL2=Sbase(1-f1)[r2-(1-f2)]=πrb2(1-f1)(r2-1+f2)  (9)
式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
(b)C-W状态
SSL0=0  (14)
SSL2=0  (16)
式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
(c)W-C状态:
SSL0=0  (21)
式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
(d)C-C状态:
SSL0=0  (28)
SSL2=0  (30)
式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
将参数带入公式可得四种不同润湿接触状态的系统自由能为:
(a)W-W状态
(b)C-W状态
(c)W-C状态
(d)C-C状态
式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,f1是一级结构的固体部分,f2是二级结构的固体部分,是接触表面中实面积与表面积之比。r1是一级结构的粗糙度系数,r2是二级结构的粗糙度系数,即实际总面积与表面积之比;ɑ1和ɑ2分别是一级、二级结构中每种柱阵列截面边长,rb为液滴与二级微纳结构表面接触圆半径;x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
23.四种不同润湿接触状态的系统无量纲自由能EI’的确定,I=1,2,3,4分别代表W-W,C-W,W-C和C-C状态:当给定重液滴和双粗糙表面的物理性质时,
是恒定的,并且由杨氏方程
可得等效接触角为 四种不同润湿接触状态的系统自由能可以简化为:
(a)W-W状态
(b)C-W状态
(c)W-C状态
(d)C-C状态
(3)计算出当前重液滴的润湿状态If,接触圆半径rb-f,液滴在一级结构中的渗透深度x1-f和液滴在级级结构中的渗透深度x2-f。
31四种接触状态下局部最小无量纲自由能的优化:
311.变量:ui,vi,i=1,2,...,N+1,液滴在一级结构x1和二级结构x2中的渗透深度;
312.约束:
i.
ii.
iii.0≤x1<h1,0≤x2<h2  (45)
iv.因为轮廓曲线函数u=f(v)是一个凸函数,所以
v.
其中,hɑmax是液滴的最大高度,V是液滴的体积,ui是离散点的横坐标,vi是离散点的纵坐标,rbmax是接触线半径的最大值,x1是液滴在一级结构中的渗透深度,x2是液滴在二级结构中的渗透深度,h1是一级结构的高度,h2是二级结构的高度;
313.搜索:
在约束条件下搜索变量x1,x2,h1,h2,利用Matlab函数“fmincon”得到最佳对象并计算出四种润湿接触状态下的部分最小无量纲自由能EImin’,I=1,2,3,4;
32.四种接触状态下整体最小无量纲自由能的优化:利用Matlab函数“fmincon”来搜索部分最小无量纲自由能EImin’,I=1,2,3,4,中的整体最小无量纲自由能Emin’,不断计算,直到整体最小无量纲自由能差收敛到10-4的可接受范围内;这时的无量纲自由能就认为是整体最小无量纲自由能,计下当前重液滴的润湿状态If,接触圆半径rb-f,液滴在一级结构中的渗透深度x1-f和
液滴在级级结构中的渗透深度x2-f;否则,加倍等分数N,再次进行步骤31和步骤32。
(4)设置初始体积增量为ΔV=0.5V,接触圆半径固定且u1=rb=rb-f,x1=x1-f,计算对应的最小无量纲自由能EV+ΔV-fix’,即EImin’。设u1=rb,v1=0,因为I,rb,x1和x2为定值时,所以最小无量纲自由能的优化为:
(a)变量:离散点坐标ui,vi,i=2,....,N+1;
(b)设置优化对象EI'min和约束条件i,ii,iii,iv以及v,同步骤312。
(c)在约束条件下搜索变量ui,vi,i=2,.....,N+1,在接触圆固定的条件下,利用Matlab中的fmincon函数得到最优对象并计算出整体最小无量纲自由能EImin’。不断进行计算,直到整体最小无量纲自由能差收敛到10-4的可接受范围内,这时的无量纲自由能就认为是整体最小无量纲自由能EImin’,即体积为ΔV+V的重液滴对应的最小无量纲自由能EV+ΔV-fix’。同时记下离散点pi-f[u0,v0]和pi-f[ui,vi]并得到相应的轮廓曲线以及相应的表面接触角通过公式 否则,就加倍等分数N,重复步骤(b)、(c)。
(5)当接触圆半径改变时,即接触圆变化时,局部最小无量纲自由能从变EV+ΔV-fix’为EV+ΔV’,重复步骤(3),计算出EV+ΔV’=Emin’。
(6)CAH是由接触圆上的能量势垒造成的,可以通过计算得出Ebarr=U0l0+U1l1+U2l2,式中,Ebarr是能量势垒,U0,U1和U2分别是液滴和基底、一级结构和二级结构间的能量势垒;l0,l1,l2分别是液滴和基底、一级结构和二级结构间的实际长度。当接触圆变化时,局部最小无量纲自由能从EV+ΔV-fix’变为EV+ΔV’。ΔE=EV+ΔV-fix’-EV+ΔV’>0可知,当接触圆变化时,接触圆获得了潜在的能量ΔE。当ΔE=Ebarr时,接触圆会向前,这是的ACA就是θɑ。
(7)通过使用二进制迭代法,可以搜索到ΔV和相应的θɑ=θV+ΔV-fix使得EV+ΔV-fix’-EV+ΔV’=Ebarr。
(8)设置初始体积增量为ΔV=-0.5V并使用相同的方法计算出θr,并计算出接触角滞后CAH=θɑ-θr。

说明书全文

一种平规则双粗糙度表面重液滴接触滞后的计算方法

技术领域

[0001] 本发明属于微纳复合结构技术领域,涉及一种水平规则双粗糙度表面重液滴接触角滞后的计算方法,具体涉及一种利用最小自由能理论和非线性优化算法计算得到水平规则双粗糙度表面重液滴的接触角滞后(CAH)。

背景技术

[0002] 对于粗糙表面上平衡的重液滴,当我们将液滴体积V减小或增加到间隔Vr流体系统的设计,具有重要的现实意义。
[0003] 目前,还未有一种明确的计算方法来计算接触角滞后(CAH)。因此,为了能够在粗糙表面上快速准确的获取液滴的接触角滞后(CAH),研究和衡量粗糙表面液滴的动态特性,迫切需要一种简单、完善的理论分析方法,计算获得液滴的接触角滞后(CAH)。

发明内容

[0004] 本发明要克服现有技术的上述缺点,提供一种水平规则双粗糙度表面重液滴接触角滞后(CAH)的计算方法,能够计算粗糙表面上重液滴的接触角滞后(CAH)。
[0005] 本发明设计了一个计算方法,利用这一方法可以得出水平规则双粗糙度表面上的重液滴接触角滞后(CAH)。对于水平规则双粗糙度表面,只要知道其液气界面张系数,液滴在一级结构、二级结构和基底上的等效接触角,液滴体积,液滴密度、在基底、一级结构和二级结构上的能量势垒以及结构尺寸,通过计算,就能够从理论上分析计算得到液滴接触角滞后(CAH),对于研究液滴在粗糙表面的动态性能和实际应用中指导微流体系统的设计具有重要指导意义。
[0006] 一种水平规则双粗糙度表面重液滴接触角滞后(CAH)的计算方法,包括以下步骤:(1)假设一个由液滴和微纳二级结构表面组成的模型,其中第一级纳米结构和第二级微米结构都是方柱,几何尺度周期性分布,二级结构和基底由不同的材料组成。模型的总表观面积为边长L0的正方形,边长远大于液滴直径。a1,b1,h1和x1分别是第一级结构方柱的边长、间距、高度和浸润深度,a2,b2,h2和x2分别是第二级结构方柱的边长、间距、高度和浸润深度,第一级结构的尺寸远小于第二级结构。由于水平规则的双粗糙度表面上的重液滴形状是轴对称的,所以在二维平面上考虑了这个问题。选择重液滴的对称轴为v轴,重液滴的表面接触面为u轴,并利用有限差分法将重液滴的轮廓曲线沿v方向分成N等份,得到离散点:
[0007]
[0008] 其中,ha是液滴的最大高度,ui是离散点的横坐标,vi是离散点的纵坐标。
[0009] (2)四种不同润湿接触状态的系统无量纲自由能函数的推导,推导过程如下:
[0010] 21.通过重液滴在一级结构和二级结构中的渗透深度,将润湿接触状态分为四种,分别为W-W,C-W,W-C和C-C状态,划分依据如下:
[0011] (a)一级结构和二级结构均处于Wenzel状态,即x1=h1,x2=h2,这种接触状态命名为W-W状态。
[0012] (b)一级结构处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,二级结构处于Wenzel状态,即0≤x1
[0013] (c)一级结构处于Wenzel状态,二级结构处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,即x1=h1,0≤x2
[0014] (d)一级结构和二级结构均处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,即0≤x1
[0015] 其中,x1是液滴在一级结构中的渗透深度,x2是液滴在二级结构中的渗透深度,h1是一级结构的高度,h2是二级结构的高度。
[0016] 22.四种不同润湿接触状态的系统自由能EI等于EI=Ea+Eb,I=1,2,3,4。下标I=1,2,3,4分别表示W-W,C-W,W-C和C-C状态。Ea和Eb分别代表重力势能和表面自由能。
[0017] Ea可以表示为:
[0018]
[0019] Eb可以表示为:
[0020] Eb=γLVSLV+γSL0SSL0+γSL1SSL1+γSL2SSL2+γSV0SSV0+γSV1SSV1+γSV2SSV2  (3)[0021] 其中,ρ是液体密度,g是重力加速度,S是界面的接触面积,γ是界面张力系数,下标LV,SL和SV分别代表液-气,固-液和固-气界面,下标1,2和0分别代表一级结构,二级结构和基底。
[0022] 液滴的外表面积Sext可以表示为:
[0023]
[0024] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度。
[0025] 液滴与基底的表面接触面积可以表示为:
[0026] Sbase=πrb2  (5)
[0027] 式中,rb为液滴与二级微纳结构表面接触圆半径。
[0028] 根据每种浸润状态的不同,可以容易求得四种状态下的SLV,SSV0,SSV1,SSV2,SSL0,SSL1,SSL2,所对应的表达式:
[0029] (a)W-W状态:
[0030] SLV=Sext=2πR2(1-cosθ)  (6)
[0031] SSL0=Sbase(1-f1)(1-f2)=πrb2(1-f1)(1-f2)  (7)
[0032] SSL1=Sbaser2[r1-(1-f1)]=πrb2r2(r1-1+f1)  (8)
[0033] SSL2=Sbase(1-f1)[r2-(1-f2)]=πrb2(1-f1)(r2-1+f2)  (9)
[0034]
[0035]
[0036]
[0037] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0038] (b)C-W状态
[0039]
[0040] SSL0=0  (14)
[0041]
[0042] SSL2=0  (16)
[0043]
[0044]
[0045]
[0046] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0047] (c)W-C状态:
[0048]
[0049] SSL0=0  (21)
[0050]
[0051]
[0052]
[0053]
[0054]
[0055] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0056] (d)C-C状态:
[0057]
[0058] SSL0=0  (28)
[0059]
[0060] SSL2=0  (30)
[0061]
[0062]
[0063]
[0064] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0065] 将参数带入公式可得四种不同润湿接触状态的系统自由能为:
[0066] (a)W-W状态
[0067]
[0068] (b)C-W状态
[0069]
[0070] (c)W-C状态
[0071]
[0072] (d)C-C状态
[0073]
[0074] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,f1是一级结构的固体部分,f2是二级结构的固体部分,是接触表面中实面积与表面积之比,r1是一级结构的粗糙度系数,r2是二级结构的粗糙度系数,即实际总面积与表面积之比,ɑ1和ɑ2分别是一级、二级结构中每种柱阵列截面边长,rb为液滴与二级微纳结构表面接触圆半径;x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0075] 23.四种不同润湿接触状态的系统无量纲自由能EI’的确定,I=1,2,3,4分别代表W-W,C-W,W-C和C-C状态:当给定重液滴和双粗糙表面的物理性质时,是恒定的,并且由杨氏方程
可得等效接触角为 四种不同润湿接触状态的系统自由能可以简化为:
[0076] (a)W-W状态
[0077]
[0078] (b)C-W状态
[0079]
[0080] (c)W-C状态
[0081]
[0082] (d)C-C状态
[0083]
[0084] 式中,θe1,θe2和θe0分别表示平坦表面上液滴在一级结构,二级结构和基底上的本征接触角;f1是一级结构的固体部分,f2是二级结构的固体部分,是接触表面中实面积与表面积之比。r1是一级结构的粗糙度系数,r2是二级结构的粗糙度系数即实际总面积与表面积之比;ɑ1和ɑ2分别是一级、二级结构中每种柱阵列截面边长,rb为液滴与二级微纳结构表面接触圆半径;x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0085] (3)计算出当前重液滴的润湿状态If,接触圆半径rb-f,液滴在一级结构中的渗透深度x1-f和液滴在级级结构中的渗透深度x2-f
[0086] 31.四种接触状态下局部最小无量纲自由能的优化:
[0087] 311.变量:离散点坐标ui,vi,i=1,2,...,N+1,液滴在一级结构x1和二级结构x2中的渗透深度;
[0088] 312.约束:
[0089] i.
[0090] ii.
[0091]
[0092] iii.0≤x1<h1,0≤x2<h2  (45)
[0093] iv.因为轮廓曲线函数u=f(v)是一个凸函数,所以
[0094]
[0095] v.
[0096] 其中,hɑmax是液滴的最大高度,V是液滴的体积,ui是离散点的横坐标,vi是离散点的纵坐标,rbmax是接触线半径的最大值,x1是液滴在一级结构中的渗透深度,x2是液滴在二级结构中的渗透深度,h1是一级结构的高度,h2是二级结构的高度;
[0097] 313.搜索:
[0098] 在约束条件下搜索变量x1,x2,h1,h2,利用Matlab函数“fmincon”得到最佳对象并计算出四种润湿接触状态下的部分最小无量纲自由能EImin’,I=1,2,3,4;
[0099] 32.四种接触状态下整体最小无量纲自由能的优化:利用Matlab函数“fmincon”来搜索部分最小无量纲自由能EImin’,I=1,2,3,4,中的整体最小无量纲自由能Emin’,不断计算,直到整体最小无量纲自由能差收敛到10-4的可接受范围内;这时的无量纲自由能就认为是整体最小无量纲自由能,计下当前重液滴的润湿状态If,接触圆半径rb-f,液滴在一级结构中的渗透深度x1-f和液滴在级级结构中的渗透深度x2-f;否则,加倍等分数N,再次进行步骤31和步骤32。
[0100] (4)设置初始体积增量为ΔV=0.5V,接触圆半径固定且u1=rb=rb-f,x1=x1-f,计算对应的最小无量纲自由能EV+ΔV-fix’,即EImin’。设u1=rb,v1=0,因为I,rb,x1和x2为定值时,所以最小无量纲自由能的优化为:
[0101] (a)变量:离散点坐标ui,vi,
[0102] (b)设置优化对象EI'min和约束条件i,ii,iii,iv以及v,同步骤312。
[0103] (c)在约束条件下搜索变量ui,vi,i=2,.....,N+1,在接触圆固定的条件下,利用Matlab中的fmincon函数得到最优对象并计算出整体最小无量纲自由能EImin’。不断进行计算,直到整体最小无量纲自由能差收敛到10-4的可接受范围内,这时的无量纲自由能就认为是整体最小无量纲自由能EImin’,即体积为ΔV+V的重液滴对应的最小无量纲自由能EV+ΔV-fix’。同时记下离散点pi-f[u0,v0]和pi-f[ui,vi]并得到相应的轮廓曲线以及相应的表面接触角通过公式 否则,就加倍等分数N,重复步骤(b),(c)。
[0104] (5)当接触圆半径改变时,即接触圆变化时,局部最小无量纲自由能从变EV+ΔV-fix’为EV+ΔV’,重复步骤(3),计算出EV+ΔV’=Emin’。
[0105] (6)CAH是由接触圆上的能量势垒造成的,可以通过计算得出Ebarr=U0l0+U1l1+U2l2,式中,Ebarr是能量势垒,U0,U1和U2分别是液滴和基底、一级结构和二级结构间的能量势垒;l0,l1,l2分别是液滴和基底、一级结构和二级结构间的实际长度。当接触圆变化时,局部最小无量纲自由能从EV+ΔV-fix’变为EV+ΔV’。ΔE=EV+ΔV-fix’-EV+ΔV’>0可知,当接触圆变化时,接触圆获得了潜在的能量ΔE。当ΔE=Ebarr时,接触圆会向前,这是的ACA就是θɑ。
[0106] (7)通过使用二进制迭代法,可以搜索到ΔV和相应的θɑ=θV+ΔV-fix使得EV+ΔV-fix’-EV+ΔV’=Ebarr。
[0107] (8)设置初始体积增量为ΔV=-0.5V并使用相同的方法计算出θr,并计算出接触角滞后CAH=θɑ-θr。
[0108] 本发明的发明点是:提出了一种计算粗糙表面上重液滴接触角滞后(CAH)的计算方法,只要知道其液气界面张力系数,液滴在一级结构、二级结构和基底上的等效接触角,液滴体积,液滴密度、在基底、一级结构和二级结构上的能量势垒以及结构尺寸,通过计算,就能够从理论上分析计算得到液滴接触角滞后(CAH),帮助研究人员进一步理解粗糙表面的动态性和粗糙表面的润湿性,对于实际应用中指导微流体系统的设计具有重要指导意义。附图说明
[0109] 图1a是本发明实施例的系统模型的侧面示意图;
[0110] 图1b是本发明实施例的三维系统模型示意图;
[0111] 图1c是本发明实施例的一级和二级结构的草图;
[0112] 图2a~图2d是本发明实施例重液滴在二级结构表面的润湿示意图,其中,图2a代表W-W状态,图2b代表C-W状态,图2c代表W-C状态,图2d代表C-C状态。
[0113] 图3是本发明实施例重液滴的半剖面图;
[0114] 图4a是重液滴体积变大后无量纲自由能和ACA的变化的示意图;
[0115] 图4b是重液滴体积变小后无量纲自由能和ACA的变化的示意图;
[0116] 图5是CAH随Bond数(Bo)的变化图。

具体实施方式

[0117] 下面结合附图并通过具体实施例对本发明作进一步说明,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
[0118] 一种水平规则双粗糙度表面重液滴接触角滞后的计算方法,包括以下步骤:
[0119] (1)我们考虑了一个规则的双粗糙度表面,一级结构和二级结构都是正方形柱,并且它们的几何尺寸都是周期性的。如图1所示,该系统由一个液滴和一个方形表面区域组成,并且方形表面区域的边长L0恒定且大于液滴直径。a1,b1,h1和x1分别是第一级结构方柱的边长、间距、高度和浸润深度,a2,b2,h2和x2分别是第二级结构方柱的边长、间距、高度和浸润深度,第一级结构的尺寸远小于第二级结构,外部界面面积和表面接触区域。f1是一级结构的固体部分,f2是二级结构的固体部分,即接触表面中实面积与表面积之比,r1是一级结构的粗糙度系数,r2是二级结构的粗糙度系数,
即实际总面积与表面积之比, 由于水平规则的双粗糙度
表面上的重液滴形状是轴对称的,所以在二维平面上考虑了这个问题。如图3所示,选择重液滴的对称轴为v轴,重液滴的表面接触面为u轴,并利用有限差分法将重液滴的轮廓曲线沿v方向分成N等份,得到离散点:
[0120]
[0121] 其中,ha是液滴的最大高度,ui是离散点的横坐标,vi是离散点的纵坐标。
[0122] (2)四种不同润湿接触状态的系统无量纲自由能函数的推导,推导过程如下:
[0123] 21.通过重液滴在一级结构和二级结构中的渗透深度,将润湿接触状态分为四种,分别为W-W,C-W,W-C和C-C状态,划分依据如下:
[0124] (a)一级结构和二级结构均处于Wenzel状态,即x1=h1,x2=h2,这种接触状态命名为W-W状态。
[0125] (b)一级结构处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,二级结构处于Wenzel状态,即0≤x1
[0126] (c)一级结构处于Wenzel状态,二级结构处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,即x1=h1,0≤x2
[0127] (d)一级结构和二级结构均处于Cassie状态或者接触状态将从Cassie状态转变为Wenzel状态,即0≤x1
[0128] 其中,x1是液滴在一级结构中的渗透深度,x2是液滴在二级结构中的渗透深度,h1是一级结构的高度,h2是二级结构的高度。
[0129] 22.四种不同润湿接触状态的系统自由能EI等于EI=Ea+Eb,I=1,2,3,4。下标I=1,2,3,4分别表示W-W,C-W,W-C和C-C状态。Ea和Eb分别代表重力势能和表面自由能。
[0130] Ea可以表示为:
[0131]
[0132] Eb可以表示为:
[0133] Eb=γLVSLV+γSL0SSL0+γSL1SSL1+γSL2SSL2+γSV0SSV0+γSV1SSV1+γSV2SSV2  (3)[0134] 其中,ρ是液体密度,g是重力加速度,S是界面的接触面积,γ是界面张力系数,下标LV,SL和SV分别代表液-气,固-液和固-气界面,下标1,2和0分别代表一级结构,二级结构和基底。
[0135] 液滴的外表面积Sext可以表示为:
[0136]
[0137] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度。
[0138] 液滴与基底的表面接触面积可以表示为:
[0139] Sbase=πrb2(5)
[0140] 式中,rb为液滴与二级微纳结构表面接触圆半径。
[0141] 根据每种浸润状态的不同,可以容易求得四种状态下的SLV,SSV0,SSV1,SSV2,SSL0,SSL1,SSL2,所对应的表达式:
[0142] (a)W-W状态:
[0143] SLV=Sext=2πR2(1-cosθ)  (6)
[0144] SSL0=Sbase(1-f1)(1-f2)=πrb2(1-f1)(1-f2)  (7)
[0145] SSL1=Sbaser2[r1-(1-f1)]=πrb2r2(r1-1+f1)  (8)
[0146] SSL2=Sbase(1-f1)[r2-(1-f2)]=πrb2(1-f1)(r2-1+f2)  (9)
[0147]
[0148]
[0149]
[0150] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0151] (b)C-W状态
[0152]
[0153] SSL0=0  (14)
[0154]
[0155] SSL2=0  (16)
[0156]
[0157]
[0158]
[0159] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0160] (c)W-C状态:
[0161]
[0162] SSL0=0  (21)
[0163]
[0164]
[0165]
[0166]
[0167]
[0168] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0169] (d)C-C状态:
[0170]
[0171] SSL0=0  (28)
[0172]
[0173] SSL2=0  (30)
[0174]
[0175]
[0176]
[0177] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0178] 将参数带入公式可得四种不同润湿接触状态的系统自由能为:
[0179] (a)W-W状态
[0180]
[0181] (b)C-W状态
[0182]
[0183] (c)W-C状态
[0184]
[0185] (d)C-C状态
[0186]
[0187] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,Sext为球冠液滴的外表面积,Sbase为液滴与基底的接触面积,r1为纳米结构的粗糙度因子,r2为微米结构的粗糙度因子,f1为纳米结构的面积分数,f2为微米结构的面积分数,L0为方形表面区域的边长,rb为液滴与二级微纳结构表面接触圆半径,x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0188] 23.四种不同润湿接触状态的系统无量纲自由能EI’的确定,I=1,2,3,4分别代表W-W,C-W,W-C和C-C状态:当给定重液滴和双粗糙表面的物理性质时,
[0189] 是恒定的,并且由杨氏方程可得等效接触角为 四种不同润湿接触状态的系统自由能可以简
化为:
[0190] (a)W-W状态
[0191]
[0192] (b)C-W状态
[0193]
[0194] (c)W-C状态
[0195]
[0196] (d)C-C状态
[0197]
[0198] 式中,u是离散点的横坐标,v是离散点的纵坐标,ha是液滴的最大高度,f1是一级结构的固体部分,f2是二级结构的固体部分,是接触表面中实面积与表面积之比。r1是一级结构的粗糙度系数,r2是二级结构的粗糙度系数即实际总面积与表面积之比,ɑ1和ɑ2分别是一级、二级结构中每种柱阵列截面边长,rb为液滴与二级微纳结构表面接触圆半径;x1为纳米结构下的浸润深度,x2为微米结构下的浸润深度。
[0199] (3)计算出当前重液滴的润湿状态If,接触圆半径rb-f,液滴在一级结构中的渗透深度x1-f和液滴在级级结构中的渗透深度x2-f。
[0200] 31.四种接触状态下局部最小无量纲自由能的优化:
[0201] 311.变量:离散点坐标ui,vi,i=1,2,...,N+1,液滴在一级结构x1和二级结构x2中的渗透深度。
[0202] 312.约束:
[0203] i.
[0204] ii.
[0205]
[0206] iii.0≤x1<h1,0≤x2<h2  (45)
[0207] iv.因为轮廓曲线函数u=f(v)是一个凸函数,所以
[0208]
[0209] v.
[0210] 其中,hɑmax是液滴的最大高度,V是液滴的体积,ui是离散点的横坐标,vi是离散点的纵坐标,rbmax是接触线半径的最大值,x1是液滴在一级结构中的渗透深度,x2是液滴在二级结构中的渗透深度,h1是一级结构的高度,h2是二级结构的高度。
[0211] 313.搜索:
[0212] 在约束条件下搜索变量x1,x2,h1,h2,利用Matlab函数“fmincon”得到最佳对象并计算出四种润湿接触状态下的部分最小无量纲自由能EI’,I=1,2,3,4。
[0213] 32.四种接触状态下整体最小无量纲自由能的优化:利用Matlab函数“fmincon”来搜索部分最小无量纲自由能EImin’,I=1,2,3,4,中的整体最小无量纲自由能Emin’,不断计-4算,直到整体最小无量纲自由能差收敛到10 的可接受范围内。这时的无量纲自由能就认为是整体最小无量纲自由能,并记下对应整体最小无量纲自由Emin’能时的离散点,并通过它们可以获得轮廓曲线和表面接触角 以及相应的润湿状态If和接触圆
半径rb-f,液滴在一级结构中的渗透深度x1-f和液滴在级级结构中的渗透深度x2-f。否则,加倍等分数N,再次进行步骤311,312和313。
[0214] (4)设置初始体积增量为ΔV=0.5V,接触圆半径固定且u1=rb=rb-f,x1=x1-f,计算对应的最小无量纲自由能EV+ΔV-fix’,即EImin’。设u1=rb,v1=0,因为I,rb,x1和x2为定值时,所以最小无量纲自由能的优化为:
[0215] (a)变量:离散点坐标ui,vi,i=2,....,N+1;
[0216] (b)设置优化对象EI'min和约束条件i,ii,iii,iv以及v,同步骤312。
[0217] (c)在约束条件下搜索变量ui,vi,i=2,.....,N+1,,在接触圆固定的条件下,利用Matlab中的fmincon函数得到最优对象并计算出整体最小无量纲自由能EImin’。不断进行计算,直到整体最小无量纲自由能差收敛到10-4的可接受范围内,这时的无量纲自由能就认为是整体最小无量纲自由能EImin’,即体积为ΔV+V的重液滴对应的最小无量纲自由能EV+ΔV-fix’。同时记下离散点pi-f[u0,v0]和pi-f[ui,vi],,得到相应的轮廓曲线以及相应的表面接触角通过公式 否则,就加倍等分数N,重复步骤(b),(c)。
[0218] (5)当接触圆半径改变时,即接触圆变化时,局部最小无量纲自由能从变EV+ΔV-fix’为EV+ΔV’,重复步骤(3),计算出EV+ΔV’=Emin’。
[0219] (6)CAH是由接触圆上的能量势垒造成的,可以通过计算得出Ebarr=U0l0+U1l1+U2l2,式中,Ebarr是能量势垒,U0,U1和U2分别是液滴和基底、一级结构和二级结构间的能量势垒;l0,l1,l2分别是液滴和基底、一级结构和二级结构间的实际长度。当接触圆变化时,局部最小无量纲自由能从EV+ΔV-fix’变为EV+ΔV’。ΔE=EV+ΔV-fix’-EV+ΔV’>0可知,当接触圆变化时,接触圆获得了潜在的能量ΔE。当ΔE=Ebarr时,接触圆会向前,这是的ACA就是θɑ。
[0220] (7)通过使用二进制迭代法,可以搜索到ΔV和相应的θɑ=θV+ΔV-fix使得EV+ΔV-fix’-EV+ΔV’=Ebarr。
[0221] (8)设置初始体积增量为ΔV=-0.5V并使用相同的方法计算出θr,并计算出接触角滞后CAH=θɑ-θr。
[0222] 实施例:
[0223] 二级结构的几何结构尺寸参数为:纳米结构方柱边长是ɑ1=100nm、规则方阵的周期距离b1=88nm、方柱的高度h1=400nm,微米结构方柱边长是ɑ2=10μm、规则方阵的周期距离b2=8.26μm、方柱的高度h2=13μm。液滴在一级结构、二级结构和基底上的本征接触角为θe0=θe1=θe2=106°,液滴和基底、一级结构和二级结构间的能量势垒为U0=U1=U2=4.2×10-5。
[0224] 如图5所示为是CAH随Bond数(Bo)的变化图,即不同的体积对应的CAH值。说明本方法对于二级结构固体粗糙表面接触角和液滴轮廓计算的可行性,同时对于不同的液体和演变的二阶结构都是适用的。
[0225] 以上所述仅为本发明的其中一个实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈