
摘要:采用 COMSOL Multiphysics软件分别就有孔抛光垫和实体研抛盘,对研抛区流体进行了有限元建模与仿真,以达到保持环超声振动对抛光试件和抛光垫间流场性能影响的研究目的,同时分别对试件保持环超声振动的频率、振幅对研抛区流体的速度、压力及气含率分布的影响规律进行了分析。研究结果表明:流体的速度与压力在超声振动下所产生的变化较为剧烈,保持环振子下方产生极强的微射流,试件下方的流速变化较为平缓,流体内的最大速度、最大压力以及最大气含率随着超声频率或振幅的增加呈线性增加,使用打孔抛光垫后, 抛光试件下方流体最大速度与采用实体研抛盘时最大速度相比减小数十倍,变化更加平稳,但最大压力却有百余倍的增加,使得试件底面处的最大气含率增大,空化效应更加显著, 在抛光过程中有利于形成化学反应自由基,从而使材料去除率得到提高。
基超声-化学机械研抛是将超声振动与化学机械抛光两者结合对加工过程中衬底材料进行平坦化的抛光技术,由于超声振动的存在,在振动边界区域会形成局部的暂时负压区,在该区域则会生成大量的空化气泡,发生空化现象。当空化气泡溃灭时会在流体内部产生冲击波和微射流,可不断清洗固体表面,剥除表面粘附的反应物和杂质,强烈时可侵蚀固体表面,对材料产生破坏,还可导致分子间强烈的相互碰撞或聚集,使流体内催化剂的结构、成分及反应活性发生显著性的变化[1-2]。
对于超声空化的研究主要通过仿真模拟[3]和实验探究[4]2种方法。王心迪等[5]基于Noltingk-Neppiras方程探究了超声频率对空化过程的影响,研究发现,最佳空化频率与超声作用时间存在负相关,而且超声场能的衰减可导致空化气泡半径减小,削弱空化效应。Jiao等[6]建立了一种求解超声增强液固反应传质系数的计算模型,并在不同的超声功率下对硅表面传质系数进行了数值模拟,仿真结果表明,传质系数随着超声功率的增大而增大,温度、反应器和超声源之间的距离等外在因素也会影响硅表面的传质系数。Xu等[7]在490kHz下对超声诱导流动的液体速度分布进行了数值模拟,结果表明,液体速度随超声功率的增加而增加,反应器内液体高度变化也会影响液体速度分布。
超声空化可在溶液中形成一种极端的物理化学环境,显著影响着电极动力学反应和电化学溶解速率[8-10]。王蕾等[11]深入分析了微孔电解加工过程中超声扰动的作用机理,指出在加工间隙区域内施加低频超声扰动更容易形成超声空化,有利于提高微孔电解加工精度。Bhat等[12]对电化学腐蚀过程中超声振动的影响进行了探究,发现超声振动明显有助于改善加工条件,提高加工精度和效率。杨治伟等[13]对微电解技术处理废水中的硝基苯进行了实验研究,结果表明,超声空化条件下的微电解方法降解废水中的硝基苯效果最佳。Ninomiya等[14]采用超声辅助芬顿反应处理废水中的有机物,发现超声在有机物分解过程中起明显促进作用。Sebastian等[15]从仿真和实验两方面分析了研抛间隙中的电解质流动状态,结果显示,超声振动可使加工区域内的压力、流速和空蚀强度都有明显增强。
本文利用COMSOL Multiphysics有限元分析软件,分别针对实体研抛盘和有孔抛光垫,对研抛区流体进行了有限元建模与仿真,研究保持环超声振动对抛光试件和抛光垫间流场性能的影响,分析了不同的超声频率、振幅对流体的速度、压力及气含率分布的影响规律。
1、超声振动对流场的影响
1.1 建立二维流场模型
图1(a)所示为本次仿真用到的超声-化学机械抛光的局部结构图,流场仿真的二维模型的具体尺寸如图(b)所示,求解域总长82mm,高115μm,试件与抛光盘之间的液膜厚为40μm。
如图2,模型两侧边界设置为压力出口,试件边界、夹持器边界与抛光盘边界固定,超声振子边界作纵向往复振动,速度v(t)为:
v(t)=H0·w0·cos(w0t)式中:H0为超声振幅;w0为超声频率。
图3所示为流场的网格剖分图,边界层为四边形网格,内部为自由剖分三角形网格,对角落处网格进行了细化,便于求解的准确性。整个求解域共划分单元51577个,其中边界层单元10701个。由于流体边界存在超声振动,流体的速度及压力变化极为剧烈,为了便于计算流体的压力和速度,现假设流体在界面上无滑移,即贴于表面的流速与表面滑移速度相同。
图1几何模型示意;图2流场仿真的边界条件;图3流场网格剖分
1.2 仿真结果分析
设定超声振动频率f=20kHz,振幅H=50μm,观察一个周期内t=T/4,T/2,3T/4、T这4个时刻的速度及压力分布分别如图4、图5所示。由图可知,超声边界下方速度变化极为剧烈,以横向剪切流为主,而在试件下方速度较小且变化平稳。流体中的压力随时间呈现周期性变化,变化范围主要集中在试件底部区域。选择试件底面中心点作为参考点,试件底部的直径线作为参考线,如图6所示。
图4模型内流体速度分布云图;图5模型内流体压力分布云图;图6模型参考点及参考线
1.2.1 超声频率对流场的影响
保持振幅H=50μm不变,对超声振动频率f进行参数化扫描,扫描范围:10kHz~30kHz,间隔5kHz。
在不同的振动频率下,参考点的速度及压力变化如图7所示。由图可知,流体在参考点的速度近似正铉周期性变化,随着超声频率的增加而增大。参考点处的压力也呈现周期性变化,随着超声频率的增加而增大,但是由于空化的存在,流体中无负压出现。
图7不同频率下参考点处的速度及压力变化;图8参考线上流体速度和压力以及气含率随频率变化
以f=20kHz为例,在一周期内流体速度、压力及气含率沿参考线的分布情况如图8(a)、8(c)、8(e)所示,最大速度、压力及气含率随频率的变化如图8(b)、8(d)、8(f)所示。可以看出在超声振动边界处速度变化极为剧烈,幅值可达到千米每秒,振动频率越大,速度变化越大,而在试件下方流体速度变化相对平稳。在非空化区,试件下方压力均匀分布,且随着振动频率的增加压力也逐渐增大,在空化区流体压力等于空化压力,在试件下方气含率也均匀分布,且随着振动频率的增加而增大。
1.2.2 超声振幅对流场的影响
保持振动频率f=20kHz不变,对超声振幅H进行参数化扫描,扫描范围:30~70μm,间隔10μm。
在不同的超声振幅下,参考点的速度及压力变化如图9所示。由图可知,流体在参考点的速度波动性上升,振幅越大,速度升高越多。参考点处的压力出现周期性变化,压力随着振幅的升高而增大。在t=T/4时刻,流体速度、压强及气含率沿参考线的分布情况如图10(a)、10(c)、10(e)所示,最大速度、压强及气含率随振幅的变化如图10(b)、10(d)、10(f)所示。从图10可知,在超声振子下方的流体速度变化剧烈,变化幅度随着振幅的增加而增大,而试件下方流体速度变化相对平稳。在非空化区流体的压力随着振幅的增加而增大(在空化区流体压力等于空化压力),气含率随着振幅的增加而增大。
图9不同振幅下参考点处的速度及压力变化;图10参考线上流体速度和压力以及气含率随振幅变化
2、多孔抛光垫对流场的影响
为模拟抛光过程中抛光垫的多孔效应,采用如图11所示的打孔抛光垫仿真模型,抛光垫厚度为1mm,其上均匀分布着直径2mm,中心距为4mm的圆柱孔。抛光垫底部与抛光盘粘接在一起。边界条件设置与不带抛光垫时基本相同,只是将抛光盘的固定壁边界换为抛光垫的固定壁边界。设定超声振动频率f=20kHz,振幅H=50μm,选择超声振子下方的一个液孔,观察一个周期内t=T/4,T/2,3T/4,T这4个时刻的速度及压力分布如图12、图13所示。
图11使用多孔抛光垫的流场模型;图12使用多孔抛光垫时流体速度分布矢量图;图13使用多孔抛光垫时流体压力分布图
从图12、13可知,当超声振子向上运动时,抛光垫边界处的剪切流从左侧流入液孔,形成了顺时针方向的涡流。当振子向下运动时,剪切流从右侧流入液孔,形成了逆时针方向的涡流。使用多孔抛光垫后,流体内的最大压力可达到7.68×1010Pa,相比于无抛光垫时增大了160倍左右。
在不同的超声振幅下,模型参考点的速度及压力变化分别如图14(a)、14(b)所示。对比图9可知,带抛光垫时流体在参考点的速度与不带抛光垫时相比明显减小,速度变化也更平稳,压力虽然变化趋势相同,但幅值却增大了百余倍。
在t=T/4时刻,当时用多孔抛光垫时流体速度、压力及气含率沿参考线的分布情况如图15(a)、15(c)、15(e)所示,最大速度、压力及气含率随振幅的变化如图15(b)、15(d)、15(f)所示。对比图9可以看出,使用多孔抛光垫后,无论是超声振子下方还是试件下方流体的速度相比于不带抛光垫时更加平稳,变化幅值更小。试件下方流体的压力相比于不带抛光垫时增大了百余倍,振幅越大,增大越多,而气含率也相应增加。
图14使用多孔抛光垫时参考点处的压力;图15使用多孔抛光垫时参考线上流体速度和压力以及气含率分布
3、结论
1)超声振动会使流体内的速度与压力剧烈变化,在振子下方产生以横向剪切流为主的极强微射流,而试件表面处的流速变化相对平缓。2)超声振动频率对流场的影响与振幅的影响类似,即随着超声频率或振幅的增加,流体内的最大速度和压力,以及最大气含率都会呈线性增加。3)使用打孔抛光垫后,试件下方流体速度减小数十倍,变化更加平稳,而最大压力却增大了百余倍;试件底面处的最大气含率也有所增大,空化效应更加强烈,更有利于在抛光过程中生成化学反应自由基与提高材料去除率。
参考文献:
[1]王蕾,赵志强.超声扰动在微孔电解加工中作用机理的研究[J].现代制造工程,2016(10):83-88.
[2]杨治伟,刘翠娟,慎爱民,等.超声空化微电解技术处理废水中硝基苯的研究[J].佳木斯大学学报(自然科学版),2007,25(3):408-410.
[3]罗经,李健.基于CFD的水轮机气蚀机制探讨[J].润滑与密封,2006(7):102-104.
[4]王国刚,孙冬柏,张秀丽,等.空泡溃灭过程中力学损伤行为[J].北京科技大学学报,2007,29(5):483-485.
[5]王心迪,韩智超,孟楠,等.浅谈液体中气泡在超声作用下的空化性质[J].物理实验,2013,33(12):37-40.
翟文杰,翟权.超声-化学机械研抛的流体动力学仿真[J].哈尔滨工程大学学报,2019,40(12):2050-2056.
基金项目:国家自然科学基金项目(51475119).
分享:
专业导论课程是各个工科专业都要开设的一门必修课程,旨在介绍本专业的研究对象、发展方向、关键技术及课程体系,通过该门课程的教学让学生对专业有一个总体上的了解,并为后续的专业课学习提供系统性指导,同时也需要培养学生树立投身于祖国工业和国防建设的伟大理想,在专业导论课的教学方面,多位教师学者都提出了很多教学改革思路[1,2,3,4,5]。
2024-07-17在网络信息发达、繁杂的时代背景下,青年大学生极易接触到不良的思潮和价值观,存在急功近利的价值取向,缺乏家国情怀,集体责任感以及坚定的社会主义人生理想信念。大学生在校期间的绝大多数时间是在非思政类课程学习中度过。习近平总书记在全国高校思想政治工作会议上强调指出“要用好课堂教学这个主渠道。
2020-12-15腐蚀是普遍存在于国民经济和生产生活中的重要破坏因素,学者们针对材料腐蚀现象进行了众多研究,其中,静态环境下的腐蚀研究已较为丰富。在流动体系中,流体对材料表面及腐蚀产物膜的力学冲刷作用、对腐蚀性离子的传质促进作用以及对界面反应过程的影响使得流动环境下的腐蚀问题更加严重和复杂。
2020-12-03在核反应堆中,控制棒驱动机构是反应堆控制和保护系统中重要的伺服机构,控制棒驱动机构是驱动控制棒组件运动、实现准确调节反应堆内反应性变化的重要部件,控制棒驱动机构具有实现反应堆启动、提升功率、保持功率、负荷跟踪、正常停堆和紧急事故停堆等重要的安全功能。
2020-11-27腐蚀是普遍存在于国民经济和生产生活中的重要破坏因素,学者们针对材料腐蚀现象进行了众多研究,其中,静态环境下的腐蚀研究已较为丰富。在流动体系中,流体对材料表面及腐蚀产物膜的力学冲刷作用、对腐蚀性离子的传质促进作用以及对界面反应过程的影响使得流动环境下的腐蚀问题更加严重和复杂。
2020-11-17房间内污染物来源包括:煤气泄漏、一氧化碳泄漏、室内装修后残留的甲醛等污染物。其中煤气和天然气、CO等污染物通过管道或者开关未关闭泄露至空间中,而室内装修残留物则主要为家具或者地板等使用的材料超出规定指标而不断残留在房间中。本文主要采用CFD的方法进行结构合理性的验证。
2020-11-17为研究一种高减压比组合式空气减压阀的流场特征,基于计算流体动力学(CFD)方法,建立减压阀数值计算模型,针对两级阀体特征参数组合的典型工况进行仿真计算,研究减压阀流场参数分布,分析柱阀间隙和引气位置对流动状态和减压比影响。计算结果表明,组合式减压阀流场状态复杂,第一级柱阀间隙对于流动损失和局部参数分布有重要影响,两级阀分别通过总压损失和流速变化实现减压过程,通过合理的特征参数组合可获得较大减压比及其调节范围。
2020-11-16气体在液体中的运动广泛存在于自然界和工程领域中,如发动机水下排气、水下爆炸引起的气泡运动、石油开采.其中的动力学过程在一个多世纪以来一直受到人们的关注[1,2].数值模拟是研究气体在液体中运动的一种有效方法[3,4].由于液体与气体之间的密度比一般远大于重液体与轻液体之间的密度比,所以气液相界面的捕捉一直是数值模拟的重点与难点.
2020-08-10理论和实践证明,环量控制[1]方法可以通过产生射流在后缘形成Coanda效应[2]来改变环量,从而产生高升力[3]。国内外一些学者针对环量控制参数的影响都进行了研究:Englar[3]、Liu等[4]经研究得出在小攻角、小动量系数下环量控制可产生较大升力的结论;张艳华等[5]通过计算流体力学数值仿真方法(CFD)研究了动量系数和攻角对NASA超临界翼型的升阻特性和附面层分离特性的影响规律。
2020-08-10流体力学是研究流体在静止和运动状态下的运动规律及流体与固壁相互作用力的一门学科,在机械、能源、汽车、化工、航空、航天等诸多领域有着极广泛的应用[1]。该课程是力学、机械等专业本科生的专业基础课程。由于其涉及数学、物理、力学等内容,概念多,公式推导繁琐,理论性强,实际应用广,对学生综合利用知识的要求较高,是普遍公认的教师难教、学生难学的课程之一[2]。
2020-08-10人气:4730
人气:3238
人气:2803
人气:2791
人气:2560
我要评论
期刊名称:应用数学和力学
期刊人气:1743
主管单位:重庆交通大学
主办单位:重庆交通大学
出版地方:重庆
专业分类:科学
国际刊号:1000-0887
国内刊号:50-1060/O3
邮发代号:78-21
创刊时间:1980年
发行周期:月刊
期刊开本:16开
见刊时间:一年半以上
影响因子:0.963
影响因子:1.730
影响因子:0.542
影响因子:0.744
影响因子:1.028
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!