流体振荡器气泡碎化特性数值模拟_别海燕.pdf

第 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)摘 要：为了探究流体振荡器中气泡碎化特性，采用重整化群(RNG)k-湍流模型和流体体积函数(VOF)多相流模型研究流体振荡器中的气液两相流。通过数值模拟分析了气泡在流场中的变形和破碎行为，考察不同流动参数对气泡破碎的影响。同时进行可视化实验，用实验结果验证了模拟的准确性。结果表明,流体振荡器可作为一种新型气泡发生装置。振荡频率越高，流体振荡器中的气泡碎化效果越好；随着液速增大，气泡变形程度加剧，破碎概率增加，生成子气泡也越多；剪切应力和涡流碰撞是导致气泡破碎的主要机制；射流振荡促使气泡向振荡腔壁面偏移，气泡初始直径越大，水平偏移量越大，越容易发生破碎。关键词：流体振荡器；气泡破碎；气液两相流；数值模拟 中图分类号：TQ05 文献标志码：A DOI：10.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 266100,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 to 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 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 breakup 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 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。基金项目：国家自然科学基金(22178329)；山东省自然科学基金(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 Gang.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。主射流在振荡腔内来回切换，周而复始形成周期性振荡。完成一次周期性振荡的时间为振荡周期，而主射流在单位时间内的振荡次数即为振荡频率 f。反馈回路的存在则有助于流体振荡器在没有任何运动部件的情况下自动实现振荡。Hao 等12对流体振荡器进行了三维大涡模拟和可视化实验，结果表明分布在振荡腔和主射流之间的涡流以及左右两侧压差会驱动主射流偏转。别海燕等13建立了正反馈式流体振荡器的三维数值模型，分析了振荡器内部流动特性。Ming 等14采用非定常雷诺平均 N-S 湍流模型研究了自激流体振荡器阵列的流体动力学和传热性能。Pandey 等15采用三维非定常雷诺平均 N-S 方程对流体振荡器进行数值研究，结果表明射流流动结构和振荡频率与已有实验结果有较好一致性，入口宽度显著影响流体振荡器性能。Zheng 等16采用雷诺应力湍流模型对横向气流和静止气流中的流体振荡器进行数值研究，阐明了流体振荡器的射流振荡机制。Tesa等17利用流体振荡器搅拌气流，再将气流注入曝气器中产生微气泡，这类微气泡可用于传热传质和废水处理。李辉等18介绍了流体振荡器耦合微孔膜产生气泡的应用进展。气液两相流是现今流体力学的研究热点，含有微气泡的流体流动在环境治理、工业应用、生物医药等19-22众多领域中起着重要作用。国内外学者已经对流体振荡器的理论和数值模拟进行了广泛研究，而对流体振荡器在气液两相流中的研究较少。因此本研究采用数值模拟方法对流体振荡器中的气液两相流进行研究，借助 ANSYS Fluent 模拟气泡在流场中的碎化行为，探究振荡频率、液速和气泡初始直径对流体中气泡运动和破碎的影响，揭示流体振荡器气泡碎化特性，阐明流体振荡器可作为一种新型气泡发生装置，并为后续流体振荡器的实验研究和性能提升提供理论基础。同时，本研究也进行了可视化实验，以此来验证数值模拟结果的准确性和可靠性。2 实 验 2.1 流体振荡器结构 流体振荡器结构如图1所示，主要包括气相入口、液相入口、喷嘴、反馈回路、振荡腔和出口支路。液相入口长度、宽度分别为 14、8 mm，喷嘴宽度为 4 mm，振荡腔长度为 34 mm，两侧反馈回路长度为 60 mm，两个出口支路长为 38 mm，宽为 5 mm，并形成一个 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 Frequency 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 Bubble 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中，使用电磁流量计调节水速，液相质量流量为 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 种不同尺寸的流体振荡器，长径比分别为 1.1、1.3、1.5 和 1.85。3.2 网格划分和无关性检验 网格质量和网格数量会影响模拟结果的准确性，所以在模拟之前需要进行网格无关性检验，以此来确保计算精度和结果准确性。如图 4所示，本研究采用六面体网格对流体振荡器进行网格划分，并对主射流区域进行局部加密，在提高网格质量的同时保证了瞬态模拟的精度和准确性。此外，分别构建网格数量为 54 977、80 975、114 125、130 617 和 150 634 的 5 组流体振荡器模型进行网格无关性检验，通过模拟直径为 5 mm 的气泡在流体振荡器中的运动变化，比较网格数量对振荡频率和气泡界面的影响。在模拟时监控出口流量随时间的变化情况，得到周期性变化的出口流量时均值，再通过快速傅里叶变换获得出口流量的脉动频域，脉动频域中出现的最大峰值为流体振荡器的振荡频率。如图 5 所示，随着网格数量增加，流体振荡器振荡频率逐渐趋于稳定。当网格数量大于 114 125 时，网格数对振荡频率几乎无影响。图 6 显示了气泡初始直径为 5 mm 时不同网格数量对气泡界面的影响。相较于其他网格数，网格数量为 130 617 时气泡界面最清晰。当网格数量为 150 634 时，气泡界面稍微变薄但不明显。因此，综合考虑计算精度和成本，最终选取网格数量为 130 617 的模型进行后续模拟计算。图 3 流体振荡器振荡腔结构 Fig.3 Oscillating 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.25 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 两相流模型对气-液相界面进行追踪。选用水和空气为流体介质，入口采用速度入口边界，出口采用压力出口。湍流指定方法采用湍流强度和水力直径，其中湍流强度 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 和 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 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-Realizable 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=1.67 Hz(b)L/D=1.3 f=1.45 Hz(c)L/D=1.5 f=1.88 Hz(d)L/D=1.85 f=3.7 Hz 位于振荡腔入口的初始气泡发生变形，均形成“颈部”结构，并在流场作用下发生破碎。破碎后形成的子气泡处于主射流区域，在偏转射流影响下沿着振荡腔边壁运动。由于主射流速度梯度较大，边壁处的子气泡表面产生黏性剪切应力，当剪切应力大于子气泡表面张力时，子气泡发生二次破碎。随着振荡频率增大，子气泡二次破碎后形成的小气泡数量增多，流体振荡器气泡碎化效果也随之增强。因此，选用振荡频率为 3.7 Hz，即长径比为 1.85 的流体振荡器进行后续模拟。4.2 液速对气泡破碎的影响 液速是影响气泡破碎的重要因素，Uesewa 等23研究发现液速越大，气泡碎化程度越高。在振荡器入口位置，注入初始直径为 4 mm 的气泡，设置不同的入口液速，相对应的喷嘴液速 v 分别为 1.05、2.14、3.20、4.29、5.41 ms1。发现不同喷嘴液速下气泡具有不同的破碎形式和破碎机制。如图 10 所示，当v=1.05 ms1时，因为液速较低，气泡发生变形但并未破碎。当 v 分别为 2.14 和 3.20 ms1时，气泡从喷嘴进入振荡腔后拉伸变形并发生凹陷，湍流旋涡与气泡间相互作用加强，气泡界面不稳定性增强，最后破碎形成子气泡。随着喷嘴液速增加到 4.29 和 5.41 ms1，气泡进入喷嘴后迅速变形并破碎成子图 10 不同喷嘴液速时气泡变形及破碎图 Fig.10 Deformation and breakup of bubbles at different liquid velocities in the nozzle(e)v=5.41 ms1(a)v=1.05 ms1(b)v=2.14 ms1(c)v=3.20 ms1(d)v=4.29 ms1 58 高 校 化 学 工 程 学 报 2023年2月 图 12 气泡内外区域速度云图 Fig.12 Velocity contours of bubble inside and outside regions(a)v=2.14 ms1(b)v=5.41 ms1 5.00 4.54 4.08 3.62 3.16 2.70 2.24 1.78 1.32 0.86 0.40 Velocity magnitude 13.0 11.9 10.8 9.7 8.6 7.5 6.4 5.3 4.2 3.1 2.0 Velocity magnitude 图 13 振荡腔内涡流区域 Fig.13 Eddy current regions in the oscillating cavity(a)v=1.05 ms1(b)v=2.14 ms1(c)v=3.20 ms1(d)v=4.29 ms1(e)v=5.41 m s1 气泡，其碎化程度明显。这是因为液速增大导致湍流强度和速度梯度增大，湍流脉动和旋涡碰撞使气泡更容易发生破碎。图11对比了v=2.14 ms1时流体振荡器中模拟和实验的气泡破碎结果。其中左边第一个黑色的图为模拟图，其余4 个为实验图，从图中可以看出，气泡进入流体振荡器后在偏转射流的作用下向左侧壁面斜向上运动，并在运动过程中不断发生破碎，形成粒径较小的子气泡。通过两者结果对比，发现实验与模拟得到的气泡破碎形式吻合较好，进一步验证了模拟的准确性。如图 12 所示，在流体振荡器局部放大速度云图中，圈出区域表示气相区域。对比喷嘴液速 v 为 2.14和 5.41 ms1的速度分布云图，气泡顶部与底部凹陷部位明显存在较大速度梯度。当 v=2.14 ms1时，气泡顶部与底部速度差达到 1.5 ms1；当 v=5.41 ms1时，气泡顶部与底部速度差接近 4.8 ms1。当速度梯度产生的剪切应力大于气泡表面张力时，气泡就会发生破碎。此外，由于速度梯度增大，气泡边缘会产生大小不一的涡流，涡流碰撞也会增强气泡碎化。如图 13 所示，振荡腔壁面存在大小不一的旋涡，箭头指示了涡流旋转方向。在不同液速下，气泡变形破碎位置均位于旋涡边缘与主射流之间。因为气泡与液体之间存在速度梯度差，加之旋涡碰撞，在二者共同作用下气泡最终碎化成子气泡，并且随着液速增大，气泡变形碎化程度更大。4.3 气泡初始直径对气泡运动的影响 气相流量增加，进入流体振荡器的气泡初始直径随之增大。为研究气泡初始直径对气泡运动的影响，保持 v=2.14 ms1，待流场稳定后在流体振荡器入口位置注入不同粒径的气泡，气泡初始直径分别为 2、图 11 模拟和实验气泡破碎对比 Fig.11 Comparison on bubble breakup in simulation and experiments 第 37 卷第 1 期 别海燕等：流体振荡器气泡碎化特性数值模拟 59 y/mm 图 15 不同气泡初始直径时气泡水平偏移量 Fig.15 Horizontal offsets of bubbles at different initial bubble diameters 0102030405060-3-2-10123Oscillating cavityInlet 2 mm 3 mm 4 mm 5 mmHorizontal offset/mm y/mm 图 14 不同气泡初始直径时气泡运动速度 Fig.14 Velocities of bubbles at different initial bubble diameters 0204060800.40.81.21.6OutletOscillating cavityInlet 2 mm 3 mm 4 mm 5 mmVelocity/(ms1)3、4、5 mm。如图 14 所示，在入口段气泡运动速度均呈先缓慢上升后快速上升趋势，并在振荡腔入口处达到最大速度。气泡从入口到达喷嘴处，流道横截面积缩小，液速增大，因此气泡运动速度随之增大。当气泡从喷嘴进入振荡腔后，液速逐渐降低，气泡也进入减速阶段。随着气泡初始直径增加，气泡运动速度和最大速度均减小。这是因为液速远高于气泡运动速度，初始直径较小的气泡受到液体阻力较小，更容易被液体裹挟前进。此外，当气泡在振荡腔内运动时，由于流体振荡器中主射流产生周期性振荡，导致射流在振荡腔内来回切换，流场发生偏转并左右振荡，气泡受到振荡流场的影响会向振荡器一侧发生偏转。如图 15 所示，当气泡初始直径为 2 mm 时，水平偏移量为 2.32。随着气泡初始直径增加，偏移量逐渐增大。当气泡初始直径达到 5 mm，气泡水平偏移量增大到 2.91 mm，表明气泡初始直径越大，在振荡腔内的偏移量越大，此时气泡越靠近振荡腔壁面，越容易受到振荡流场的影响而发生破碎。5 结 论 采用 RNG k-湍流模型和 VOF 两相流模型对流体振荡器进行数值模拟，研究了气液两相流中气泡碎化特性，分析得到以下结论：(1)设计了长径比为 1.1、1.3、1.5 和 1.85 的流体振荡器，对应振荡频率分别为 1.67、1.45、1.88、3.7 Hz，表明随着长径比增大，流体振荡器内振荡频率先减后增。通过 4 种振荡频率下的气泡破碎形式对比，发现振荡频率越高，初始气泡破碎越剧烈，流体振荡器气泡碎化效果越好。(2)液速显著影响气泡破碎，当液速较低时，气泡发生变形而未破碎；随着液速增大，气泡变形和碎化程度显著增强，破碎形成的子气泡增多。剪切应力和涡流碰撞是气泡破碎的主导机制。(3)流体振荡器产生的射流振荡导致气泡偏移，进而影响气泡运动过程。气泡初始直径越大，水平偏移量越大，气泡越靠近振荡腔壁面，发生破碎概率越高。参考文献：1 王忠,李甘牛,陈滨.一种单出口流体振荡器件特性实验研究 J.南京航空航天大学学报,2015,47(6):877-883.WANG Z,LI G N,CHEN B.Experimental investigation on characteristics of single-outlet fluidic oscillators J.Journal of Nanjing University of Aeronautics&Astronautics,2015,47(6):877-883.2 刘影,李春鹏,张铁军,等.后缘连续偏转机翼振荡射流控制的数值模拟研究 J.航空科学技术,2020,31(5):36-43.LIU Y,LI C P,ZHANG T J,et al.Numerical simulation of oscillating jet control for trailing edge continuous deflection wing J.Aeronautical Science&Technology,2020,31(5):36-43.3 WU Y J,YU S F,ZUO L.Large eddy simulation analysis of the heat transfer enhancement using self-oscillating fluidic oscillators J.International Journal of Heat and Mass Transfer,2019,131:463-471.4 赵钰,龚盼,刘志刚,等.新型射流振荡器工作原理与试验研究 J.流体机械,2018,46(9):13-17,83.60 高 校 化 学 工 程 学 报 2023年2月 ZHAO Y,GONG P,LIU Z G,et al.Working principle and experimental research of a new type fluidic oscillator J.Fluid Machinery,2018,46(9):13-17,83.5 HE J F,YIN K,PENG J M,et al.Design and feasibility analysis of a fluidic jet oscillator with application to horizontal directional well drilling J.Journal of Natural Gas Science and Engineering,2015,27(3):1723-1731.6 张超.生物膜反应器内微生物运动及附着特性研究 D.重庆:重庆大学,2015.ZHANG C.Study on bacterial locomotion and attachment in biofilm reactor D.Chongqing:Chongqing University,2015.7 赵世坤.集成微流体振荡器的大范围梯度稀释芯片的构建及其在生物学中的应用 D.上海:上海交通大学,2017.ZHAO S K.A mircofluidic oscillator based long-range gradient generator and its application in biology D.Shanghai:Shanghai Jiaotong University,2017.8 WANG J Y,WANG L G,HANOTU J,et al.Improving the performance of coal flotation using oscillatory air supply J.Fuel Processing Technology,2017,165:131-137.9 HARUN M H C,ZIMMERMAN W B.Membrane defouling using microbubbles generated by fluidic oscillation J.Water Science and Technology:Water Supply,2019,19(1):97-106.10 WANG J M,JIN Q Q,ZHANG Y Y,et al.Reducing the membrane fouling in cross-flow filtration using a facile fluidic oscillator J.Separation and Purification Technology,2021,272:118854.11 GREGORY J W,TOMAC M N.A review of fluidic oscillator development and application for flow control:43rd AIAA Fluid Dynamics Conference C.San Diego:American Institute of Aeronautics and Astronautics,2013.12 HAO Z R,LIU G,WANG Y,et al.Studies on the drive mechanism of the main jet deflection inside a fluidic oscillator J.Industrial&Engineering Chemistry Research,2020,59(20):9629-9641.13 别海燕,黄晨,安维中,等.正反馈式流体振荡器内部流动特性的数值模拟 J.化工学报,2021,72(3):1504-1511.BIE H Y,HUANG C,AN W Z,et al.Numerical simulation of the internal flow characteristics in a feedback fluidic oscillator J.CIESC Journal,2021,72(3):1504-1511.14 MING T Z,WANG Z Y,LIAO X W,et al.Unsteady RANS simulation of fluid dynamic and heat transfer in an oblique self-oscillating fluidic oscillator array J.International Journal of Heat and Mass Transfer,2021,177:121515.15 PANDEY R J,KIM K Y.Numerical modeling of internal flow in a fluidic oscillator J.Journal of Mechanical Science and Technology,2018,32(3):1041-1048.16 ZHENG J G,XIA L,HU J,et al.Numerical investigation of characteristics of fluidic oscillator operating in quiescent air and a crossflow J.Aerospace Science and Technology,2021,113(3):106731.17 TESA V.Fluidic oscillators mediating generation of microbubbles(Survey)J.Fluids,2021,6(2):77.18 李辉,李干禄,何峰,等.微气泡耦合生物反应器的研究进展 J.化工进展,2020,39(12):4758-4765.LI H,LI G L,HE F,et al.Research progress of microbubble coupled bioreactor J.Chemical Industry and Engineering Progress,2020,39(12):4758-4765.19 冯俊杰,孙冰,姜杰,等.射流场中单个气泡破裂的 CFD 模拟与 PIV 实验研究 J.高校化学工程学报,2018,32(2):286-294.FENG J J,SUN B,JIANG J,et al.CFD simulation and PIV investigation on single bubble breakup in turbulent jet flow J.Journal of Chemical Engineering of Chinese Universities,2018,32(2):286-294.20 MOVAHED S,SARMAH A K.Global trends and characteristics of nano-and micro-bubbles research in environmental engineering over the past two decades:A scientometric analysis J.Science of the Total Environment,2021,785:147362.21 HUANG Z L,CHEN Q Y,YAO Y,et al.Micro-bubbles enhanced removal of diesel oil from the contaminated soil in washing/flushing with surfactant and additives J.Journal of Environmental Management,2021,290:112570.22 LACERDA Q,TANTAWI M,LEEPER D B,et al.Emerging applications of ultrasound-contrast agents in radiation therapy J.Ultrasound in Medicine&Biology,2021,47(6):1465-1474.23 UESAWA S,KANEKO A,NOM