流体振荡器气泡碎化特性数值模拟_别海燕.pdf
《流体振荡器气泡碎化特性数值模拟_别海燕.pdf》由会员分享,可在线阅读,更多相关《流体振荡器气泡碎化特性数值模拟_别海燕.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 37 卷第 1 期 高 校 化 学 工 程 学 报 No.1 Vol.37 2023 年 2 月 Journal of Chemical Engineering of Chinese Universities Feb.2023 文章编号:1003-9015(2023)01-0053-08 流体振荡器气泡碎化特性数值模拟 别海燕1,李云霞1,安维中1,李玉龙1,薛立成1,林子昕1,郝宗睿2,刘 刚2(1.中国海洋大学 化学化工学院,山东 青岛 266100;2.齐鲁工业大学(山东省科学院)海洋仪器仪表研究所,山东 青岛 266100)摘 要:为了探究流体振荡器中气泡碎化特性,采用重整化群(R
2、NG)k-湍流模型和流体体积函数(VOF)多相流模型研究流体振荡器中的气液两相流。通过数值模拟分析了气泡在流场中的变形和破碎行为,考察不同流动参数对气泡破碎的影响。同时进行可视化实验,用实验结果验证了模拟的准确性。结果表明,流体振荡器可作为一种新型气泡发生装置。振荡频率越高,流体振荡器中的气泡碎化效果越好;随着液速增大,气泡变形程度加剧,破碎概率增加,生成子气泡也越多;剪切应力和涡流碰撞是导致气泡破碎的主要机制;射流振荡促使气泡向振荡腔壁面偏移,气泡初始直径越大,水平偏移量越大,越容易发生破碎。关键词:流体振荡器;气泡破碎;气液两相流;数值模拟 中图分类号:TQ05 文献标志码:A DOI:1
3、0.3969/j.issn.1003-9015.2023.01.007 Numerical simulation of bubble breakup characteristics in a fluidic oscillator BIE Hai-yan1,LI Yun-xia1,AN Wei-zhong1,LI Yu-long1,XUE Li-cheng1,LIN Zi-xin1,HAO Zong-rui2,LIU Gang2(1.College of Chemistry and Chemical Engineering,Ocean University of China,Qingdao 26
4、6100,China;2.Institute of Oceanographic Instrumentation,Qilu University of Technology(Shandong Academy of Sciences),Qingdao 266100,China)Abstract:In order to explore the characteristics of bubble breakup in a fluidic oscillator,the RNG k-turbulence model and the VOF multiphase flow model were used t
5、o simulate the gas-liquid two-phase flow.The deformation and breakup behaviors of bubbles in the flow field were analyzed by numerical simulation,and the effects of different flow parameters on bubble breakup were investigated.Meanwhile,a visualized experiment was carried out to verify the accuracy
6、of simulation results.The results show that the fluidic oscillator could be used as a new type of bubble generator.High oscillation frequency results in better effect of bubble breakup.With the increase of liquid velocity,the degree of bubble deformation intensifies and the possibility of bubble bre
7、akup increases with more new bubbles generated.Shear stress and eddy current are the main mechanisms of bubble breakup.In addition,jet oscillation causes bubble to deviate towards the wall of oscillating cavity.The horizontal offset of the bubble increases at relatively larger initial diameter,which
8、 indicates the bubble has a higher probability to break up.Key words:fluidic oscillator;bubble breakup;gas-liquid two-phase flow;numerical simulation 1 前 言 流体振荡器是一种能在出口产生连续振荡射流的流体器件,具有无运动部件、结构简单、运行稳定、操作方便、安全可靠等优点,可以应用于航空航天1-2、强化传热3、石油化工4-5、生物能源6-7、煤矿 收稿日期:2021-10-26;修订日期:2022-01-20。基金项目:国家自然科学基金(221
9、78329);山东省自然科学基金(ZR2020ME175,ZR2020QE192);中央高校基本科研业务费专项(202165002)。作者简介:别海燕(1982-),女,山东日照人,中国海洋大学副教授,博士。通信联系人:安维中,Email: 引用本文:别海燕,李云霞,安维中,李玉龙,薛立成,林子昕,郝宗睿,刘刚.流体振荡器气泡碎化特性数值模拟 J.高校化学工程学报,2023,37(1):53-60.Citation:BIE Hai-yan,LI Yun-xia,AN Wei-zhong,LI Yu-long,XUE Li-cheng,LIN Zi-xin,HAO Zong-rui,LIU Ga
10、ng.Numerical simulation of bubble breakup characteristic in a fluidic oscillator J.Journal of Chemical Engineering of Chinese Universities,2023,37(1):53-60.54 高 校 化 学 工 程 学 报 2023年2月 浮选8、膜除污9-10等领域。在流体振荡器中,流体流入喷嘴后形成主射流,由于射流扩散和卷吸作用,这股射流会发生偏转并被低压吸引到邻近壁面上附壁流动,形成科恩达效应11。主射流在振荡腔内来回切换,周而复始形成周期性振荡。完成一次周期性振
11、荡的时间为振荡周期,而主射流在单位时间内的振荡次数即为振荡频率 f。反馈回路的存在则有助于流体振荡器在没有任何运动部件的情况下自动实现振荡。Hao 等12对流体振荡器进行了三维大涡模拟和可视化实验,结果表明分布在振荡腔和主射流之间的涡流以及左右两侧压差会驱动主射流偏转。别海燕等13建立了正反馈式流体振荡器的三维数值模型,分析了振荡器内部流动特性。Ming 等14采用非定常雷诺平均 N-S 湍流模型研究了自激流体振荡器阵列的流体动力学和传热性能。Pandey 等15采用三维非定常雷诺平均 N-S 方程对流体振荡器进行数值研究,结果表明射流流动结构和振荡频率与已有实验结果有较好一致性,入口宽度显著
12、影响流体振荡器性能。Zheng 等16采用雷诺应力湍流模型对横向气流和静止气流中的流体振荡器进行数值研究,阐明了流体振荡器的射流振荡机制。Tesa等17利用流体振荡器搅拌气流,再将气流注入曝气器中产生微气泡,这类微气泡可用于传热传质和废水处理。李辉等18介绍了流体振荡器耦合微孔膜产生气泡的应用进展。气液两相流是现今流体力学的研究热点,含有微气泡的流体流动在环境治理、工业应用、生物医药等19-22众多领域中起着重要作用。国内外学者已经对流体振荡器的理论和数值模拟进行了广泛研究,而对流体振荡器在气液两相流中的研究较少。因此本研究采用数值模拟方法对流体振荡器中的气液两相流进行研究,借助 ANSYS
13、Fluent 模拟气泡在流场中的碎化行为,探究振荡频率、液速和气泡初始直径对流体中气泡运动和破碎的影响,揭示流体振荡器气泡碎化特性,阐明流体振荡器可作为一种新型气泡发生装置,并为后续流体振荡器的实验研究和性能提升提供理论基础。同时,本研究也进行了可视化实验,以此来验证数值模拟结果的准确性和可靠性。2 实 验 2.1 流体振荡器结构 流体振荡器结构如图1所示,主要包括气相入口、液相入口、喷嘴、反馈回路、振荡腔和出口支路。液相入口长度、宽度分别为 14、8 mm,喷嘴宽度为 4 mm,振荡腔长度为 34 mm,两侧反馈回路长度为 60 mm,两个出口支路长为 38 mm,宽为 5 mm,并形成一个
14、 34夹角。为了保证流体充分发展,采用 3D 打印技术制作一个长度为 300 mm 的射流管道。流体振荡器则用聚甲基丙烯酸甲酯(polymethyl methacrylate,PMMA)制作而成。2.2 实验设置 实验在室温下进行,以水和空气为工作介质,具体流程如图 2 所示。水由离心泵输送到流体振荡器图 2 实验装置示意图 Fig.2 Schematic diagram of the experimental setup Control and data acquisition system Air storage Water tank Gas volume flowmeter Freque
15、ncy converter Electromagnetic flowmeter Centrifugal pump Light source Fluidic oscillator High-speed camera Gas inlet Fluid inlet Nozzle Feedback loop Oscillating cavity Outlet 图 1 流体振荡器结构示意图 Fig.1 Schematic diagram of the fluidic oscillator 第 37 卷第 1 期 别海燕等:流体振荡器气泡碎化特性数值模拟 55 图 6 不同网格数量气泡界面 Fig.6 Bu
16、bble interfaces with different grid numbers(a)54 977(b)80 975(c)114 125(d)130 617(e)150 634 Grid numbers 图 5 网格无关性检验 Fig.5 Grid independence verification 图 4 流体振荡器网格划分示意图 Fig.4 Grid mesh of the fluidic oscillator Oscillation frequency/Hz 5497780975114125 1306171506342.752.802.852.902.953.003.05中,使用电
17、磁流量计调节水速,液相质量流量为 0.1 0.5 kgs1。空气由注射器注入气相入口,使用气体流量计调节气速,气相体积流量为 0.010.05 Lmin1。高速摄像系统主要由计算机、高速相机(Photron FASTCAM Mini WX)和 LED 光源组成,通过高速拍摄可以获得清晰完整的气泡运动过程。3 数值模拟 3.1 物理模型 流体振荡器的几何结构对内部流动结构和性能有重要影响,通过改变流体振荡器振荡腔的长径比(L/D),可以得到不同尺寸参数的流体振荡器。如图 3 所示,L 为振荡腔长度,D 为反馈管之间的距离。保持振荡腔长度 L=34 mm,设计 4 种不同尺寸的流体振荡器,长径比分
18、别为 1.1、1.3、1.5 和 1.85。3.2 网格划分和无关性检验 网格质量和网格数量会影响模拟结果的准确性,所以在模拟之前需要进行网格无关性检验,以此来确保计算精度和结果准确性。如图 4所示,本研究采用六面体网格对流体振荡器进行网格划分,并对主射流区域进行局部加密,在提高网格质量的同时保证了瞬态模拟的精度和准确性。此外,分别构建网格数量为 54 977、80 975、114 125、130 617 和 150 634 的 5 组流体振荡器模型进行网格无关性检验,通过模拟直径为 5 mm 的气泡在流体振荡器中的运动变化,比较网格数量对振荡频率和气泡界面的影响。在模拟时监控出口流量随时间的
19、变化情况,得到周期性变化的出口流量时均值,再通过快速傅里叶变换获得出口流量的脉动频域,脉动频域中出现的最大峰值为流体振荡器的振荡频率。如图 5 所示,随着网格数量增加,流体振荡器振荡频率逐渐趋于稳定。当网格数量大于 114 125 时,网格数对振荡频率几乎无影响。图 6 显示了气泡初始直径为 5 mm 时不同网格数量对气泡界面的影响。相较于其他网格数,网格数量为 130 617 时气泡界面最清晰。当网格数量为 150 634 时,气泡界面稍微变薄但不明显。因此,综合考虑计算精度和成本,最终选取网格数量为 130 617 的模型进行后续模拟计算。图 3 流体振荡器振荡腔结构 Fig.3 Osci
20、llating cavity structure of the fluidic oscillator L y x o D y 56 高 校 化 学 工 程 学 报 2023年2月 3.3 模型验证 由于流体振荡器内部流场复杂,流动具有强烈的各向异性,因此需要选择适当的湍流模型对其进行数值模拟。目前常用的湍流模型有 Standard k-模型、Realizable k-模型和 RNG k-模型。而 RNG k-模型引入了时均应变率,加入了旋流的影响,对计算精度有较大改善。为了验证湍流模型的准确性和可靠性,搭建实验装置,对比不同湍流模型的气泡上升速度数值模拟结果和实验结果。其中,入口液速均为0.2
21、5 ms1,气泡初始直径为 4 mm。如图 7(a)所示,图中 y 为从入口到出口的距离,相比于 Standard k-模型和 Realizable k-模型,RNG k-模型得到的气泡上升速度结果与实验结果具有较好的一致性。图 7(b)显示了 3 种模型的模拟值与实验值的相对误差,从图中可以发现 RNG k-模型相对误差最小。因此,选用 RNG k-湍流模型对流体振荡器进行模拟计算。3.4 边界条件和求解设置 采用 Fluent 15.0 的 RNG k-湍流模型对流体振荡器进行非稳态模拟,待流场稳定后采用 VOF 两相流模型对气-液相界面进行追踪。选用水和空气为流体介质,入口采用速度入口边
22、界,出口采用压力出口。湍流指定方法采用湍流强度和水力直径,其中湍流强度 I 由式(1)计算得出。压力-速度耦合采用压力耦合方程组的半隐式方法,体积分数函数采用几何重构法,压力差值采用 PRESTO,压力离散采用二阶迎风格式,瞬态方程采用一阶隐式。计算时间步长设为 1.0 105 s,每步迭代 20 次,以保证较高稳定性和收敛精度。1/80.16-IRe=(1)式中:Re 为雷诺数。4 结果与讨论 4.1 振荡频率对气泡破碎的影响 振荡频率 f 是衡量流体振荡器的重要性能参数,而流体振荡器的结构形状对振荡频率有很大影响。设计 4 种流体振荡器,长径比(L/D)分别为 1.1、1.3、1.5 和
23、1.85。如图 8 所示为不同振荡腔长径比对振荡频率的影响,随着长径比增大,振荡频率呈先下降后上升的趋势。为了研究振荡频率对气泡碎化效果的影响,对不同结构的流体振荡器进行气泡破碎数值模拟。图 9(a)中,长径比为1.1 时,f=1.67 Hz,此时初始气泡由喷嘴进入振荡腔,从气泡主体直接脱离出一个较小子气泡,随后初始气泡主体受射流偏转的影响靠近壁面移动,并再次发生较为明显的碎化。在图 9(b)、9(c)和 9(d)中,长径比为 1.3、1.5 和 1.85 时,对应的振荡频率分别为 1.45、1.88 和 3.7 Hz。L/D 图 8 振荡腔长径比对振荡频率的影响 Fig.8 Effects
24、of length to diameter ratio of oscillating cavity on oscillation frequency Oscillation frequency/Hz 1.11.31.51.851.01.52.02.53.03.54.0图 7 湍流模型验证与相对误差分析 Fig.7 Turbulence model validation and relative error analysis y/mm(a)velocity of bubbles Relative error/%020406080048121620 RNG k-Standard k-Realiza
25、ble k-y/mm(b)error analysis Velocity/(ms1)0204060800.20.40.60.81.0 experiment RNG k-Standard k-Realizable k-experiment RNG k-Standard k-Realizable k-RNG k-Standard k-Realizable k-第 37 卷第 1 期 别海燕等:流体振荡器气泡碎化特性数值模拟 57 图 9 不同振荡频率时气泡破碎情况 Fig.9 Bubble breakup at varying oscillation frequencies(a)L/D=1.1 f
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 流体 振荡器 气泡 特性 数值 模拟 海燕
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。