[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] 以上所述仅为本发明的其中一个实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何
修改、等同替换、改进等,均包含在本发明的保护范围内。