基于格子BoltzMann方法的自驱动Janus颗粒扩散泳力∗

周光雨 陈力 张鸿雁 崔海航

(西安建筑科技大学环境与市政工程学院,西安710055)

基于格子BoltzMann方法的自驱动Janus颗粒扩散泳力∗

周光雨 陈力†张鸿雁 崔海航

(西安建筑科技大学环境与市政工程学院,西安710055)

(2016年10月28日收到;2017年1月23日收到修改稿)

Janus颗粒的自驱动力研究对于纳微米尺度驱动力课题具有重要意义,本文针对Pt-SiO2型Janus颗粒,基于格子Boltzmann模型及动量交换法提出了计算其扩散泳力的方法,通过与实验数据对比修正验证了模型准确性,并通过分析证明了此类Janus颗粒的扩散泳力与速度场无关,进一步模拟比较了不同形状颗粒的自驱运动.分析发现,对于体积相等形状不同的Janus颗粒,扩散泳力主要由轴线投影面积决定,此外反应面积也会对扩散泳力产生影响.

格子Boltzmann方法,Janus颗粒,数值模拟,扩散泳力

1 引言

发动机等动力系统对于现代生活是必不可少的,在微观尺度上,纳米技术的发展使得微尺度马达的实现成为可能.近些年来,科学家已开始尝试仿照生物大分子来制造纳米机器或分子机器[1].2000年,Soong等[2]利用ATP合成酶与螺旋桨组装,使其能够利用ATP驱动旋转.然而,设计合成纳微米尺度的动力系统具有很大难度,其难点不仅在于这种尺度的复杂,更是因为宏观的热工转换的驱动原理在微观的尺度下很难实现[3].在此背景下,以自生梯度为驱动力的马达开始引起人们的重视,本文的研究对象Janus颗粒就是此类一个典型的例子.

由具有不同物理化学性质的两部分所组成的微米级别的粒子被称之为Janus颗粒[4].当Janus颗粒在特定条件下的流域中发生反应时,由于两侧性质不同,颗粒周围会形成不对称的浓度场、温度场、电磁场等,其中利用不对称浓度场形成的驱动力称为扩散泳力,它可以推动Janus颗粒运动,将化学能转化为机械能[5].利用Janus颗粒自驱性制造的纳微米尺度马达在药物的靶向输运[6]及微流控装置中的自主载体[7]等方面具有重要前景.在实验方面,国内外学者如Howse等[8],Ke等[9],Zheng等[10]以及宫春亮[11]在运动机理方面进行了大量研究,但仍没有统一的认识.相比之下,针对Janus颗粒的数值模拟研究相对较少.Córdova-Figueroa和Brady[12]引入渗透压的概念,认为周围溶质粒子浓度的改变产生了渗透力,从而导致了Janus颗粒的运动.de Buyl和Kapral[13]则提出了介观粒子模型,认为Janus颗粒是一种球状体组合颗粒并结合分子动力学与多体碰撞力学模拟了Janus颗粒在溶液发生放热反应和分解反应时颗粒的运动机理.但Córdova-Figueroa、Brady及de Buyl等提出描述理论进行的数值模拟并未与具体实验结果对比;胡静等[14]进行了Janus微球分数布朗运动的颗粒动力学模拟,崔海航等[15]进行了不同形状Janus颗粒的自驱动特性模拟,但其采用基于实验数据的半经验的计算方法无法从机理方面很好地分析扩散泳力.

格子Boltzmann方法(lattice boltzmann method,LBM)是自20世纪90年代兴起的流体计算和建模的一种方法,是基于流体微观粒子特性和微观动理论方程的介观统计方法,在一定条件下可还原为Navier-Stokes方程[16].LBM的微观粒子背景使得它可以比较直观地处理流体内部以及流体与周围环境的相互作用,因而LBM在多相和多组分流及微尺度流等领域发展迅速.1993年,Shan和Chen[17]提出了一种能够直接刻画粒子间相互作用的多组分格子模型,随后对该模型做出进一步改进[18],并得到了广泛应用.张任良等[19]分析了基于Shan-Chen模型的格子Boltzmann方法在微米尺度范围内流动模拟问题的有效性,拓展了格子方法在微流动模拟方面的应用.在边界处理方面,史冬岩等[20]研究了任意复杂流固边界的格子处理方法,提出了具有较高精度的方法.综上所述,在Janus颗粒的驱动这一微观界面动力问题的描述方面,格子Boltzmann方法比传统的数值方法更有优势.

Janus颗粒运动的扩散泳力计算模拟是这一课题的核心问题之一.如何准确地模拟微尺度下巨大的梯度量产生的扩散泳力是本课题的主要研究内容.本文将采取LBM中特有的反弹边界格式中动量交换的方式,根据动量定理计算浓度分布函数对壁面的作用力.此外,由于Janus颗粒处于微米尺度,其运动不可避免地受到布朗力的影响,这将与扩散泳力叠加而形成更加复杂的运动.本文的工作主要关注扩散泳力的计算,这将为后续Janus颗粒运动的进一步研究奠定基础.

本文选择Pt-SiO2型Janus颗粒作为研究对象.首先简要描述微球在H2O2溶液中的自驱动现象、模拟采用的各类模型及扩散泳力计算方法;之后进行数值模拟并结合实验数据验证模型的合理性;最后分析Janus颗粒自驱运动,计算比较不同形状的Janus颗粒自驱动.

2 Pt-SiO2型Janus颗粒自驱运动分析

Pt-SiO2型Janus颗粒通过在SiO2微粒一侧镀一层Pt膜获得,将Janus颗粒置于H2O2溶液中,Pt侧将作为催化剂参与化学反应侧不发生反应.当H2O2浓度较低,Janus颗粒较小时,生成的O2以溶质分子状态存在于溶液中[21].使得Janus颗粒两侧溶质浓度产生显著不对称变化,高浓度端分子数多,分子对颗粒的碰撞强度大于低浓度端,提供给颗粒的动量多于低浓度端,导致Janus颗粒朝低浓度端,即SiO2一侧移动.原理如图1所示.

图1 (网刊彩色)Pt-SiO2型Janus颗粒驱动原理示意Fig.1.(color on line)ScheMatic d iagraMof selfp ropu lsion of Janus particle.

通过实验观察可以看出Janus颗粒的运动轨迹具有明显的随机性,但在局部区域内则具有定向运动的特征.经统计得出微球的运动速度VJanus与观察时间间隔∆t的关系[22],如图2(a)所示.根据VJanus的变化趋势可以将微球的运动划分为三个阶段,分别由布朗运动、自驱动及类布朗运动所主导.可以看出,当观察时间间隔∆t为0.1—1.0 s量级时处于自驱动阶段,此时微球的时均速度VJanus近似保持恒定,约为3—5µm/s.从所对应的运动轨迹则可以看出,此时微球近似做匀速直线运动.而当观察时间间隔很短或很长时,微球的运动都呈现出明显的随机性.

由于Janus颗粒的密度高于外部环境溶液,颗粒会在靠近下壁面处的平面上运动,由此产生的偏转角φ、距底面高度δ如图2(b),即壁面效应也会对Janus颗粒运动产生影响[23].

图2 (a)Janus颗粒在不同下的平均速度V Janus[22];(b)Janus颗粒近壁姿态Fig.2.(a)The average speed of Janus particles under d iff erent tiMe intervals[22];(b)the d iagraMof Janus particles near the wall.

3 计算模型

针对扩散泳型Janus颗粒自驱运动这类复杂问题,需要涉及流场模型、浓度场模型、表面化学反应模型及扩散泳力模型,本部分将依序介绍各类所需模型.

3.1 双分布函数格子BoltzMann模型

对于Janus颗粒表面反应,单位H2O2分解放热98.2 kJ[24],溶液比热较大,通过对反应温度场的预先数值模拟验证,发现温度梯度过小,远不能产生实验中观察到的粒子运动,因此认为温度不是引起Janus颗粒运动的原因,可以忽略反应放热对粒子运动的影响.主要关注浓度场的模拟,模拟采用双分布函数模型[18],其中流场模拟演化方程为

式中fi(x,t)是流场粒子速度分布函数,c=δx/δt是粒子速度,δx和δt是单位网格长度和时间步长,ei是离散速度,τ是无量纲松弛时间,是平衡态分布函数,本文采用D2Q9模型,此时

平衡态分布函数为

其中ωi为权系数,在本模型中ω0=4/9,当i=1,2,3,4时,ωi=1/9,当i=5,6,7,8时,ωi=1/36;ρ为流体密度,u为流体速度,该模型对应的宏观方程为

υ为运动黏性系数,定义为

流体密度及动量可通过下式计算得到

双分布函数模型的浓度场模拟演化方程为

其中gi(x,t)是浓度场速度分布函数,τs是无量纲松弛时间,是平衡态分布函数,本文采用D2Q9模型,平衡态分布函数为

其中u为流场速度,C为摩尔浓度.该模型对应的宏观方程为

M为溶质对应摩尔质量,ur为溶质流速.

3.2 边界条件格式

流场边界采用Ladd提出的Half-Way反弹格式[25,26].该格式可表示为

其中ub是壁面速度.

对于浓度场边界条件,一般可表示为

其中∂C/∂n表示的是界面法线方向上的浓度梯度.类似Half-Way反弹格式,本文采用张婷[27]提出的一类半反弹格式.首先对于狄利克雷边界条件,即b1=0,b2̸=0.可通过插值得到边界流体格点xf处的未知分布函数

其中Cf是与边界结点相邻的流体结点上的宏观浓度,n表示界面的法向.求解(15)式,则边界上的浓度为

3.3 化学反应模型

发生在流固界面上的表面反应可以通过以下边界条件来描述:

其中D为反应物/生成物的扩散系数,k为反应速率即边界上反应物/生成物的通量值.对于Pt催化H2O2分解反应

其反应过程可以分解为两步[8].第一步,H2O2分子被吸附在Pt的表面,吸附速率为k1;第二步,被吸附在Pt表面的H2O2分子发生分解反应生成H2O和O2,分解速率为k2.整个反应的H2O2反应速率k为

式中,[H2O2]vol为H2O2的体积百分数,k1=4.4×1011/(µm2·s),k2=4.8×1010/(µm2·s)换算成摩尔数即除以阿伏伽德罗常数6.02×1023,得k1=7.3×10−13mol/(µm2·s),k2=7.9×10−14mol/(µm2·s).O2反应速率为H2O2的0.5倍.

3.4 扩散泳力的计算

1994年由Ladd提出动量交换法用以计算流体-固体颗粒间相互作用[25,26],是通过Half-Way反弹格式实现的.

如图3所示,格点xf流体粒子碰撞后的分布函数为经过δt/2时间后,流体粒子运动到边界点xb处并与之发生反弹碰撞

在流体粒子和固体壁面的碰撞过程中,流体粒子的动量发生变化,碰撞前后的动量变化是

按照动量定理,流体粒子对固体粒子施加的作用力是

图3 Half-W ay反弹格式示意图Fig.3.ScheMatic d iagraMof Half-W ay bounce-back scheMe.

Janus颗粒的运动是由于溶质粒子对Janus颗粒的碰撞而产生的作用.溶质粒子与固体颗粒的相互作用可类比于流固体颗粒的相互作用.格点xf溶质粒子碰撞后的分布函数为经过δt/2时间后,溶质粒子运动到边界点xb处并与之发生反弹碰撞

溶质分子的动量可通过(10)式计算得出,溶质粒子和固体壁面的碰撞过程中,碰撞前后的动量变化为

同理,溶质粒子对固体粒子作用力计算公式如下:

其中α为动量修正系数,大小与溶液浓度、微尺度效应等有关,具体值可根据实验结果确定.

积分计算出最终三维球体的扩散泳力.

3.5 计算流程

结合上述各计算模型,整个模拟过程可分为以下几步:

1)设定各初始物理参数,并将其转为格子系统参数;

2)计算流场,得出流场速度信息;

3)执行浓度场碰撞过程,其中平衡态分布函数中流速为流场速度;

4)执行浓度场流动过程,并利用(18)式分别计算H2O2,O2反应通量、利用(13)式计算浓度场边界条件、利用(22)式和(23)式计算扩散泳力;

5)计算宏观摩尔浓度;

6)进入浓度场的下一个时间步,回到步骤3);

7)进入流场的下一个时间步,回到步骤2).

4 数值模拟

4.1 Janus颗粒自驱动扩散泳力模拟

选择Janus颗粒在初始浓度为2.5%的H2O2溶液中的自驱动性能进行研究.计算域为边长为20µm的正方形,网格数目为800×800,即网格尺寸∆x=2.5×10−8m.取直径为2µm的Janus颗粒位于正方形中心,粒子左侧为反应侧,右侧为不反应侧,如图4(a)所示.曲边界采用格线反弹格式,如图4(b)所示,设距圆心距离小于半径的格点为固体格点,大于半径的格点为流体格点,边界点为流固格点中点,Janus颗粒直径为80个格子,可以提供足够的精度.

图4 (a)几何模型;(b)曲面边界格点分布图Fig.4.(a)GeoMetric Model;(b)scheMatic d iagraMof surface boundary grid points.

流场采用2.2节的格子BoltzMann模型,对于H2O2溶液,其浓度较小,假定运动黏度与水相等,即υ=1.01×10−6m2/s,令∆t=4.53×10−10s,根据松弛时间与运动黏度关系式υ=(τ−1/2)c2∆t/3,松弛时间为1.24.

实验中溶液静止,Janus颗粒做匀速直线运动,在数值模拟中可采用相对坐标进行研究,流场右侧流体流入,流速为3µm/s,左侧流体流出,上下两侧为滑动边界条件,Janus颗粒边界为无滑移边界条件.模拟得出速度场.

浓度场采用2.2节的格子Boltzmann模型,O2与H2O2的扩散系数分别为2.3×10−9,1.4×10−9m2/s[12],设O2浓度场的松弛时间为1,根据松弛时间与扩散系数的关系D=(τ−1/2)c2∆t/3,可得时间步长∆t=4.53×10−8s,并由上述关系计算得出H2O2浓度场的松弛时间为0.804.

H2O2初始浓度为2.5%,即735.29mol/m3,O2初始浓度为0mol/m3.边界采用3.2节的边界处理格式,Janus颗粒左侧边界条件为−D(∂C/∂n)=k,反应速率参考3.3节计算方法.Janus颗粒右侧及浓度场其余边界为0浓度梯度边界.

由此可以模拟得出浓度分布,并得出最终稳定时的浓度场,如图5所示.按3.4节方法可计算Janus颗粒的扩散泳力.

Janus颗粒运动速率较小,Pe数(Pe=ul/D)为10−2量级[22],表明溶质输运过程中,扩散输运要远大于对流输运,假定流场速度均为0,重新计算浓度场及扩散泳力,并比较两种方法下的扩散泳力,计算结果如图6所示.可以看到两者结果一致,故可以忽略流场对浓度场的影响,在后面的模拟中假定流场速度均为0,仅使用浓度场模型以简化模拟.

图5 (网刊彩色)稳定时O 2浓度场(Mol/M3)Fig.5.(color on line)The graph of steady O 2 concentration field.

图6 扩散泳力随时间步的变化Fig.6.The relationship curves of the diff usiophores and tiMe.

4.2 Janus颗粒自驱动扩散泳力模型的验证

由于Janus颗粒在局部区域可看作定向匀速运动,此时认为颗粒的随机布朗运动可以忽略,从而简化分析过程,底壁面对Janus颗粒的影响以壁面影响力表示.受力分析如图7所示,可以列出力平衡方程:水平方向

垂直方向

图7 Janus颗粒受力分析图Fig.7.ScheMatic d raw ing of forces analysis to Janus particle.

表1 Janus 颗粒不同初始浓度模拟结果Table 1. Simulation results about Janus particles with different initial concentrations.

因此,可用水平方向上黏滞阻力验证模拟得出的扩散泳力结果.对此阶段Janus颗粒的自驱动实验分析[22,23],在2.5%,5.0%和10.0%的溶液中统计速度分别为3.0,4.4,4.8µm/s,偏转角分别为15.5°,12.0°,9.5°.

此时颗粒是在流态为低雷诺数的流动,黏滞阻力可用斯托克斯公式计算,即

式中µ为流体动力黏度,r为颗粒半径.

分别计算粒径2µMJanus颗粒在初始浓度2.5%,5.0%,10.0%的H2O2溶液中的黏滞阻力及扩散泳力,最终计算结果整理在表1中.

分析表1结果可发现,数值模拟得出的扩散泳力与黏滞阻力在数量级上一致,模拟结果可用修正系数α来修正,且修正系数α大小随溶液浓度增大而减小.据模拟结果可以得出修正系数α与溶液浓度C间的函数关系:

4.3 Janus颗粒自驱运动分析

4.4 不同形状的Janus颗粒自驱动模拟

球形Janus颗粒由于制备简单,一直是Janus颗粒研究的重点.但随着研究的深入,不同的形状对于Janus颗粒的自驱特性的影响已越来越引起人们的关注,且因其制备上的复杂性,使得数值模拟研究具有重要意义.本部分将选取相同体积下不同形状的Janus颗粒,研究形状对扩散泳力及自驱动速度的影响.

不考虑偏转角及其他壁面效应和重力的影响,对水平方向上的受力分析如图8所示,扩散泳力与黏滞阻力相等.溶液浓度、微尺度效应等影响因素通过动量修正系数修正.建立如图9所示的不同形状Pt-SiO2型Janus颗粒模型,其表面组成均为一半催化(即Pt材料)与一半非催化材料(即SiO2材料),使5种颗粒的体积相同,为4π/3µm3,即对应于直径为2µm的圆球.椭球a各轴长度分别为m=2/3µm,n=2/3µm,l=9/4µm;椭球b各轴长度分别为m=1/2µm,n=1/2µm,l=4µm;圆柱a取其半径为r=1/2µm,长为l=16/3µm;圆柱b取其半径为1µm,长为l=4/3µm.圆柱b与圆球轴线投影面积相同,椭球b与圆柱a轴线投影面积相同.

图8 不同形状Janus颗粒受力分析图Fig.8.ScheMatic draw ing of forces analysis to d iff erent-shaped Janus particle.

设初始溶液浓度为2.5%,修正系数α为2.67.经过相同的LBM模拟过程可以得到扩散泳力,对于黏滞阻力,非球形不能使用斯托克斯公式计算,本部分使用CoMsol数值模拟平台快速建模,模拟流场计算粒子自驱动速度.

计算域为70µm×70µm×75µm的长方体,粒子位于长方体中心,数值模型采用相对坐标进行研究,粒子保持静止,流场速度为自驱动速度VJanus,边界条件设置为:长方体右侧壁面流体流入,流速u=−VJanus,长方体左侧壁面流体流出,其他四个壁面为滑动壁面.网格设置为标准网格,进行稳态计算.

图9 不同形状Janus颗粒(a)圆球;(b)椭球a;(c)椭球b;(d)圆柱a;(e)圆柱bFig.9.d iff erent shapes of Janus particles:(a)Ball;(b)spheroid a;(c)spheroid b;(d)cy linder a;(e)cylinder b.

对于球型Janus颗粒,自驱动速度为3µm/s,阻力模拟结果为5.65×10−14N,与斯托克斯公式算出的结果相同.对于其他类型Janus颗粒,通过改变自驱动速度的大小,使得阻力等于扩散泳力,得出自驱动速度.

比较不同形状Janus颗粒的反应面积、沿轴线投影面积及扩散泳力、自驱动速度,模拟结果整理在表2中.

可以看到对于体积相同、形状不同的Janus颗粒,自驱动速度、扩散泳力主要与轴线方向上的投影面积呈正相关.Janus颗粒径向方向上的扩散泳力对称,会相互抵消掉,只有相当于轴线方向上投影面积的扩散泳力会驱动Janus颗粒运动.

进一步分析相同体积相同投影面积的Janus颗粒,比较圆柱b与圆球及椭球b与圆柱a,发现圆柱形要大于球形和椭球.反应面积更大的形状,扩散泳力更大,自驱动速度更快,但相对于轴线投影面积的影响较小.

表2 不同形状Janus颗粒模拟结果Tab le 2.SiMulation results about Janus particlesw ith diff erent shapes.

5 结论

Janus颗粒所受的扩散泳力及其运动速率对于Janus颗粒自驱动课题的研究具有重要意义,本文采用格子Boltzmann方法建立了Janus颗粒的自驱动模型,并针对这一问题进行了研究,得到以下主要结论.

1)本文采用动量交换法计算扩散泳力,模拟分析了粒径2µMJanus颗粒的自驱运动,比较了不同初始溶液浓度下,壁面稳定浓度、反应速率及扩散泳力的大小.并通过与实验数据对比,验证了扩散泳力模型的合理性.

2)对于体积相同、形状不同的Janus颗粒,扩散泳力、自驱动速度主要与轴线投影面积呈正相关.对于相同体积相同投影面积的Janus颗粒,反应面积大的形状,扩散泳力、自驱动速度更大.

Janus颗粒运动因受布朗运动影响,整体会呈现一定随机性,因而局部的定向速度为统计结果,使得模拟结果无法达到一个准确的修正系数α.此外,对于不同形状的Janus粒子忽略了壁面效应及重力等影响使得定量分析不够精确.未来的工作将发展为三维、动网格、叠加布朗作用的Janus颗粒运动,使得模拟更贴近真实实验,为Janus粒子研究的进一步发展奠定基础.

[1]Zhao Y P 2012 Physica l Mechanics of Surfaces and Interfaces(Beijing:Science Press)p399(in Chinese)[赵亚溥2012表面与界物理力学(北京:科学版社)第399页]

[2]Soong R K,Bachand G D,Neves H P,O lkhovets A G,C raighead H G,MonteMagno C D 2000 Science 290 1555

[3]W ang W,Duan W,AhMed S,Mallouk T E,Sen A 2013 Nano Today 8 531

[4]Jiang S,G ranick S,Schneider H J 2012 Janus Particle Syn thesis,Self AsseMbly and Applications(USA:RSC Pub lishing Press)pp1–25

[5]Chernyak V G,Starikov SA,Beresnev SA 2001 J.Appl.Mech.Tech.Phys.42 445

[6]Patra D,Sengup ta S,Duan W,Zhang H,Pav lick R,Sen A 2013 Nanoscale 5 1273

[7]Rückner G,K ap ral R 2007 Phys.Rev.Lett.98 150603

[8]How se J R,Jones R A,Ryan A J,Gough T,Vafabakhsh R,Golestanian R 2007 Phys.Rev.Lett.99 048102

[9]K e H,Ye S,Carroll R L,Showalter K 2010 J.Phys.Chem.A 114 5462

[10]Zheng X,Hagen B T,Kaiser A,W u M,Cui H H,Silber-Li Z,Löwen H 2013 Phys.Rev.E 88 032304

[11]Gong C L 2013 M.S.Thesis(X ian:X i’an University of A rchitectu re and Technology)(in Chinese)[宫春亮2013硕士学位论文(西安:西安建筑科技大学)]

[12]Córdova-Figueroa U M,B rady J F 2008 Phys.Rev.Lett.100 158303

[13]de Buy l P,Kap ral R 2013 Nanosca le 5 1337

[14]Hu J,Zhang H Y,Zheng X,Cui H H 2014 Chinese J.HydrodynaMics 04 377(in Chinese)[胡静,张鸿雁,郑旭,崔海航2014水动力学研究与进展04 377]

[15]Cui H H,Tan X J,Zhang H Y 2014 Nanotechno logy and Precision Engineering 12 340(in Chinese)[崔海航,谭晓君,张鸿雁2014纳米技术与精密工程12 340]

[16]Guo Z L,Zheng C G 2009 Theory and Applications of Lattice BoltzMann Method(Beijing:Science Press)p10(in Chinese)[郭照立,郑楚光2009格子BoltzMann方法的原理及应用(北京:科学出版社)第10页]

[17]Shan X,Chen H 1993 Phys.Rev.E 47 1815

[18]Shan X,Doolen G 1995 J.Stat.Phys.81 379

[19]Zhang R L,D i Q W,Wang X L,D ing W P,Gong W 2012 Mechanics in Engineering 2 10(in Chinese)[张任良,狄勤丰,王新亮,丁伟朋,龚玮2012力学与实践2 10]

[20]Shi D Y,Wang Z K,Zhang A M2014 Acta Phys.Sin.63 07403(in Chinese)[史冬岩,王志凯,张阿漫2014物理学报63 074703]

[21]W ang L L,Cui H H,Zhang J,Zheng X,W ang L,Chen L 2016 Acta Phys.Sin.65 220201(in Chinese)[王雷磊,崔海航,张静,郑旭,王磊,陈力2016物理学报65 220201]

[22]W u ML 2014 M.S.Thesis(X i’an:X i’an University of A rchitectu re and Technology)(in Chinese)[武美玲2014硕士学位论文(西安:西安建筑科技大学)]

[23]Cui H H,Tan X J,Zhang H Y,Chen L 2015 Acta Phys.Sin.64 134705(in Chinese)[崔海航,谭晓君,张鸿雁,陈力2015物理学报64 134705]

[24]Casson V,Maschio G 2011 Ind.Eng.Chem.Res.51 7526

[25]Ladd A J C 1994 J.F luid Mech.271 285

[26]Ladd A J C 1994 J.Fluid Mech.271 311

[27]Zhang T 2001 Ph.D.D issertation(W uhan:Huazhong University of Science and Technology)(in Chinese)[张婷2012博士学位论文(武汉:华中科技大学)]

(Received 28 October 2016;revised Manuscrip t received 23 January 2017)

PACS:47.63.Mf,07.10.Cm,02.60.Cb,47.70.FwDOI:10.7498/aps.66.084703

*Pro ject supported by the National Natural Science Foundation of China for EMergency ManageMent Pro jects(G rant No.11447133),the Young Scientists Fund of the National Natural Science Foundation of China(G rant No.11602187),the Pro ject of the Natu ral Science Foundation of Shaanxi Province for You th Talent,China(G rant No.2016JQ 1008),the Scientifi c Research PrograMFunded by Shanxi Provincial Education DepartMent,China(G rant No.15JK 1385),and the Pro ject froMState Key Laboratory of Building Science and Technology in W estern China.

†Corresponding author.E-Mail:jasonchencl@163.com

Research on d iff usiophoresis of self-p ropu lsion Janus particles based on lattice BoltzMann Method∗

Zhou Guang-Yu Chen Li†Zhang Hong-Yan Cui Hai-Hang

(School of EnvironMent and Municipal Engineering,X i’an University of Architecture and Technology,X i’an 710055,China)

Studies of the driving force of the self-p ropu lsion Janus particles are very iMportant in the fields ofMicro-power and nano-Motor.In this paper,we choose the Micron Pt-SiO2-type Janus particle as a research ob ject,which is propelled by self-generated concentration gradient in the dilute solution of H2O2,focusing on the self-p ropu lsion of the single particle.According to the force analysis of the Janus particle,the surface force can be decoMposed into the viscous resistance of the fluid,the Brownian force derived froMthe Molecular therMal fl uctuation,and the diff usiophoresis caused by the diff usion of the solute coMponent.TheMain aiMof this paper is to find theway to accurately simu late the diff usiophoresis generated by the huge concentration gradient on aMicroscale.The lattice Boltzmann method(LBM)is aModern MesoscopicMethod based on theMicroscopic particle characteristics of the fl uid,which Makes itMore intuitive to dealw ith the interaction between the fluid and solid.It ismore advantageous than the traditional nuMericalMethod in the description of thisMicro-interface dynaMic p roblem,i.e.,the self-propu lsion of Janus particle.On a certain time scale,when the Janus particle show s the directionalMotion,the influence of the Brownian force can be ignored.Thus,the analytical process can be siMp lified.Based on themomentuMtheorem,themethod of calculating the diff usiophoresis produced by concentration diff usion is proposed.We introduce the MoMentuMexchange in the half-way bounce-back scheMe of LBMinto theModel of themulticoMponent diff usion and reaction.Through counting the surface forcewe can obtain the diff usiophoresis acting on the Janus particle.Moreover,this diff usiophoresismodel ismodified by coMparing the experiMental fluid resistance w ith simu lated one.This coMparision verifies the validity of the diff usiophoresisModel.Then,the analysis of the variation of diff usiophoresis p roves that the value of diff usiophoresis is independent of the fluid velocity.Through the further app lication of thismodel,the diff erent shapes of Janus particles w ith the same volume are coMpared in simu lations.The results show that the self-diff usiophoresis isMain ly deterMined by the axial projection area.In addition,the reaction area of the particle also aff ects the value of the diff usiophoresis.

lattice Boltzmann method,Janus particles,numerical simulation,diffusiophoresis

10.7498/aps.66.084703

∗国家自然科学基金应急管理项目(批准号:11447133)、国家自然科学基金青年科学基金(批准号:11602187)、陕西省自然科学基础研究计划青年人才项目(批准号:2016JQ 1008)、陕西省教育厅专项科研计划(批准号:15JK 1385)和西部绿色建筑国家重点实验室培育基地自主科研项目资助的课题.

†通信作者.E-Mail:jasonchencl@163.com

©2017中国物理学会C h inese P hysica l Society

http://w u lixb.iphy.ac.cn