Processing math: 100%

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

浸没射流高韦伯数条件下气溶胶水洗去除

吕强 田超 佟立丽 曹学武

吴峥嵘, 施龙波, 江国栋, 等. 谐振环和矩形软波导性能研究[J]. 强激光与粒子束, 2025, 37: 023004. doi: 10.11884/HPLPB202537.240310
引用本文: 吕强, 田超, 佟立丽, 等. 浸没射流高韦伯数条件下气溶胶水洗去除[J]. 强激光与粒子束, 2025, 37: 026003. doi: 10.11884/HPLPB202537.240150
Wu Zhengrong, Shi Longbo, Jiang Guodong, et al. Performance study of resonant ring and rectangular flexible waveguide[J]. High Power Laser and Particle Beams, 2025, 37: 023004. doi: 10.11884/HPLPB202537.240310
Citation: Lü Qiang, Tian Chao, Tong Lili, et al. Aerosol scrubbing removal under high Weber number immersed jet condition[J]. High Power Laser and Particle Beams, 2025, 37: 026003. doi: 10.11884/HPLPB202537.240150

浸没射流高韦伯数条件下气溶胶水洗去除

doi: 10.11884/HPLPB202537.240150
详细信息
    作者简介:

    吕 强,lvqiang1227@sjtu.edu.cn

    通讯作者:

    佟立丽,lltong@sjtu.edu.cn

  • 中图分类号: TN124

Aerosol scrubbing removal under high Weber number immersed jet condition

  • 摘要: 基于射流穿透长度与夹带液滴分数等水力模型,建立了混合气体射流条件下气溶胶惯性碰撞去除模型,通过对注入区空间离散构建了高韦伯数射流流型含气芯夹带特征的气溶胶水洗去除分析方法。采用气溶胶水洗机理实验装置开展的两组蒸汽份额为64%,淹没深度为0.7 m,质量通量分别为217 kg/(m2·s)和120 kg/(m2·s)实验,以及引用淹没深度为0.5 m、质量通量为95 kg/(m2·s)的纯不凝性气体载带气溶胶的RCA2实验结果,对分析方法进行了验证。结果表明:考虑射流水力学特征的气溶胶水洗去除模型的预测结果与实验值符合较好,通过参数分析发现随着射流韦伯数的增加,射流穿透长度和夹带液滴分数均增加,增强了气溶胶与液滴的惯性碰撞作用。
  • 在加速器驱动嬗变研究装置(CiADS)[1-2]的超导直线加速器中,单根传输线的功率需求达到了上百千瓦,此时已经不适合继续选择同轴馈管进行功率传输。在这种情况下,考虑到实际应用频率为650 MHz,本文采用WR1500矩形波导作为主馈线[3-4]。为了消除矩形波导引入带来的误差,实现微波设备和负载端的指向微调,必须要考虑软波导[5]。为提高射频系统的运行稳定性,本文对WR1500软波导的电性能和热力学性能进行了综合研究,旨在为CiADS传输线设计提供参考。

    为了测试软波导在实际应用场景中的性能,需要对软波导进行几十甚至上百千瓦功率的老炼测试。这对于一般的连续波设备及负载提出了很高的要求。而在高功率输入耦合器和介质窗的老炼中,采用了一种基于行波谐振环或者驻波谐振环的老炼方法,可以有效降低离线老炼的功率需求[6-8]。单端输入的行波谐振环等效功率增益低,增益一般仅为20倍左右[9],而驻波谐振老炼方法虽然增益高,但是该方法中使用的可移动短路板结构复杂且高功率承受能力差,会导致驻波谐振老炼的稳定性较差。基于此,我们设计了一种基于谐振环的双端口输入的高功率老炼测试平台,该方案结构简单可靠,且有望能大幅度提高等效功率增益。测试平台结构如图1所示,谐振环主体为一个10 dB波导定向耦合器。

    图  1  基于谐振环的双端口测试平台结构示意图
    Figure  1.  Schematic structure of the resonant ring-based dual-port test platform

    波导定向耦合器可以等效为一个四端口器件,其S参量是4×4的矩阵。我们所使用的定向耦合器耦合度为9.95 dB,插损小于0.1 dB,方向隔离度大于40 dB,将其视为理想定向耦合器,则有

    [S]=[0TC0T00CC00T0CT0] (1)

    式中C为定向耦合器的耦合系数;T为定向耦合器的传输系数。

    射频信号从定向耦合器的端口3经过谐振环到达端口4,相当于经过了一个双端口器件,则有

    S43=eα+jφ (2)

    式中:α为谐振环的全环损耗;φ为整个谐振环的相移。

    当谐振环达到稳态后,存在如下关系[10]

    V3S43T+V1C=V3 (3)

    则该谐振环的电压增益为

    GV=V3V1=C1S43T (4)

    当谐振环的周长为2πn时,该谐振环的最大增益为

    GVm=C1Teα (5)
    GPm=(C1Teα)2 (6)

    式中GVm为最大电压增益;GPm为最大功率增益。

    当端口1和端口2同时输入功率,则在环内形成驻波,波腹位置的电压增益为

    GV=V3+V4V1=2C1S43T (7)

    此时谐振环内部形成驻波场,在驻波场波腹点处,谐振环的增益为

    GVm=2GVm (8)
    GPm=4GPm (9)

    通过改变两个输入信号的相位关系,可以移动波腹位置,重复老炼过程可以使待测设备的所有位置都被驻波场的波腹老炼到。

    为了确保谐振环在650 MHz时实现最大增益,在设计测试平台时需要特别关注谐振环的相移。考虑到WR1500波导自重较大,不易安装和拆卸,我们采取了固定谐振环主体的设计方案,在测试中只需要在固定节点处更换待测产品。谐振环主要部件的结构如图2所示。在靠近10 dB定向耦合器的位置我们设置了参考软波导,不同的测试产品除了记录自身的温度变化外,都需要与参考软波导进行对比。同时参考软波导为谐振环提供了±12 mm的长度调节余量,以简化整个环路相位调节过程。

    图  2  650 MHz谐振环的主要部件结构图
    Figure  2.  Diagram of the main components of the 650 MHz resonant ring

    测试平台实物如图3所示。我们使用两个氮化镓(GaN)功放模块分别为谐振环的两个输入端口提供功率,每个功放模块的最大输出功率为2800 W,这为测试平台提供了足够的功率裕量,能够支持大功率条件下的各种实验[11]。负载与环形器集成在功放模块内部,以防止测试过程中功放模块被烧毁。低电平系统可以输出两路同步信号,通过调节两路信号的相位差达到移动驻波波腹的目的。

    图  3  基于谐振环的双端口测试平台实物图
    Figure  3.  Physical diagram of the resonant ring-based dual-port test platform

    测试平台搭建完成后,使用矢量网络分析仪对整个系统进行小信号测试。网分的端口1接10 dB定向耦合器的输入,端口2接谐振环上功率监测定向耦合器的端口,该定向耦合器的耦合度为−65.12 dB。测试结果见图4,计算得环内行波功率增益为13.42 dB。

    图  4  谐振环增益测试
    Figure  4.  Test plot of the resonant ring gain

    为了确保高功率传输线的长期稳定性和可靠性,我们选择了来自不同制造商的软波导样品进行高功率测试与对比。除了国内制造商,我们还特别引入了来自国外的同类产品作为参照,以国际产品的性能表现来标定和对比国内产品。待测产品分别由铜和铝材料组成。铜因其优良的电导率和导热性,在传输线中被广泛应用,而铝则因其质轻和良好的机械强度成为国外制造商的首选材料。

    高功率测试中,每个待测产品都会单独进行测试,待一个产品测试完成后,再更换另一个产品重新测试。测试为两个主要部分,首先是单端口输入功率的行波测试。在这项测试中,只使用一个功放输出功率,另一个功放处于关闭状态,通过上位机设置低电平系统单路输出,并同时记录功放与谐振环上功率监测定耦的读数。当功放输出达到最大时,记录谐振环内的功率。经过30 min的稳定运行后,整个系统的发热趋于稳定,此时记录每个样品和参考软波导的热分析图。该测试旨在评估每个样品在行波条件下的电热性能。

    第二部分测试为双端口输入功率的驻波测试。在这一测试中,首先将一个功放的输出功率设置为1.5 kW。待系统稳定后,缓慢增加另一个功放的输出至相同功率,并调整其相位,以使环内驻波波腹移动。同时通过低电平系统监测前向功率和反向功率的相位,通过两路监测信号的相位差来计算波腹的具体位置。当波腹位置位于待测产品中心后,等待30 min记录测试产品和参考软波导的热分析图。这种测试方式旨在评估样品在实际功率应用场景下的性能表现。

    高功率测试结果如表1 所示,环境温度均为25 ℃,其中产品4为来自国外制造商的参考软波导。在行波测试中,功放输出功率为2 kW,环内功率为46 kW。在驻波测试中,每个功放的输出均为1.5 kW,环内驻波波腹处等效功率为138 kW。图5为4组产品的驻波测试热分析图。通过测试我们发现,虽然国产的WR1500软波导已经达到了基本可用的水平,能够满足一般应用的需求,但与国外制造商的样品相比,仍存在一定的性能差距。具体来说,在高功率条件下,国产软波导的电热损耗明显高于国外产品。这种较高的电热损耗意味着,在相同的功率传输条件下,国产软波导会消耗更多的能量,转化为不必要的热量,这不仅会降低系统的整体效率,还可能引发热管理问题。造成这种差距的原因可能是多方面的,包括材料选择、制造工艺的精细程度、结构设计的合理性等。国外制造商在软波导设计和制造方面拥有更为成熟的技术和经验,其产品在高功率应用中表现出更好的电热性能和稳定性。为此,我们认为有必要对国产软波导的结构设计进行进一步的优化和改进。

    表  1  待测软波导
    Table  1.  Flexible waveguide under test
    material manufacturing location max. temperature rise (46 kW)/℃ max. temperature rise (138 kW)/℃
    copper domestic 9.7 18.5
    copper domestic 8.2 17.5
    aluminum domestic 6.5 16.8
    aluminum foreign 4.8 11.3
    下载: 导出CSV 
    | 显示表格
    图  5  驻波测试中4组产品的热分析图
    Figure  5.  Thermograms of 4 products in standing wave testing

    软波导结构如图6所示,图中a1b1a2b2分别为软波导内外截面尺寸,rR分别为软波导内外截面的圆角半径。软波导的表面具有周期性波纹结构,H为波纹深度,T为波纹宽度,S为波纹间距,电磁波在这种结构中会产生一定的反射[12]。软波导的柔韧性是由波纹结构的波纹参数决定的:波纹深度H越大,波峰宽度T和波谷间距S越小,软波导的柔软性越好。为了满足软波导电性能要求,HTS必须均小于λg/10λg为650 MHz的导波波长0.58 m)。

    图  6  矩形软波导结构尺寸示意图
    Figure  6.  Schematic dimensions of rectangular flexible waveguide structure

    波纹结构的尺寸选取决定了软波导是否单模传输、传输带宽内的驻波系数、衰减以及可承受的功率。根据测试经验,本文选择H=9 mm、T=4 mm、S=4 mm。波纹结构虽然使波导柔软,但也造成了略大于对应硬波导的驻波系数。但是由于其HTS远小于工作波长,所以其不均匀性很小,由此引起的反射也很小。

    为了减小软波导与硬波导连接时的反射以及阻抗匹配的问题,软波导的法兰盘设计要确保与硬波导连接时具备相同的口径,以保证传输线的整体性能[13]。标准WR1500波导的内口径尺寸的长边和短边分别为381 mm和190.5 mm。考虑到生产设备的实际情况,我们优化时选取a1=379mmb1=185mmr=13.5mm,软波导的外口径尺寸a2b2分别为a2=a1+2H=397mmb2=b1+2H=203mmR=31mm。矩形软波导内表面的周期性波纹结构会导致纵向电流的传导路径增长,进而增加损耗。优化后的结构与之前的产品相比,纵向电流的传导路径缩短了50%,降低了损耗。

    根据所选择的软波导尺寸结构,软波导的驻波比VSWR理论值可表示为

    RVSW=a×(b1+H)b×(a1+H)×1(λg2a1)21(λg2a11.717×r2b1)2<1.01 (10)

    所设计的软波导RVSW<1.01,说明此尺寸设计比较合理。但是由于实际生产过程中的种种因素,实际生产出的软波导的RVSW会高于设计值。

    此外矩形软波导的工作频率与其内口径的尺寸有关,波导尺寸的选择应保证其单模传输。矩形软波导中TEmn模式的截止波长为

    λc,mn=2(ma1)2+(nb1)2 (11)

    式中:m代表电磁场沿波导宽边a的驻波最大数;n代表电磁场沿波导窄边b的驻波最大数。矩形波导的基模为TE10模式,第一高阶模为TE20模式,则单模工作条件为

    λc,20<λ0<λc,10 (12)

    通过计算,该软波导的单模工作频率为450 MHz< f <800 MHz,满足WR1500波导的工作频带要求。使用CST软件 [14]对优化后的结构进行了仿真,对应的传输系数如图7所示。在目标频带内,优化后的结构显示出优异的性能。

    图  7  矩形软波导工作频带内传输系数仿真曲线
    Figure  7.  Simulation curves of transmission coefficients in the operating band of rectangular flexible waveguide

    根据以上的结构和尺寸设计,成型了400 mm长的WR1500软波导,制造完成后的软波导实物照片如图8所示。为了对比测试,本文选择了两个不同的制造商,分别使用铜和铝制造样件。在现有的测试平台上进行高功率测试,测试结果如表2所示。由于更换了软波导组件,环内损耗降低,增益增大。单端口输入时功放模块输出2 kW,环内功率50 kW。双端口输入功率时每个功放模块功放输出为1.4 kW,环内驻波波腹等效功率为140 kW。产品2和产品3的热分析图如图9所示。

    图  8  结构优化后的软波导实物图
    Figure  8.  Image of the optimized flexible waveguide
    表  2  待测器件
    Table  2.  Flexible waveguide under test
    material manufacturing location max. temperature rise (50 kW)/℃ max. temperature rise (140 kW)/℃
    copper domestic 10.7 20.2
    aluminum domestic 4.3 9.4
    aluminum foreign 5.3 10.7
    下载: 导出CSV 
    | 显示表格
    图  9  140 kW驻波测试时两款软波导热分布对比
    Figure  9.  Comparison of heat distribution of two flexible waveguides at 140 kW standing wave test

    对大尺寸的波纹周期结构来说,由于铜质地较软,很难在加工过程中保持其自身的形状[15]。接缝焊接时,铜需要更高的温度,而由于结构限制,无法通过增加壁厚的方式来提高结构强度。在壁厚0.5 mm的情况下很容易造成焊接后的产品变形,不利于大批量生产。在测试中,测试产品1的内壁具有明显的不均匀变形,这也是结构和功率相同时,产品1反而热损耗更大的原因。而所制造的铝材质WR1500软波导电热性能优于国外制造商。

    本文搭建了一套基于谐振环的高功率测试平台,单端口输入时的功率增益为13.4 dB;双端口输入时,波腹位置的等效功率增益提升至19.4 dB。通过这一测试平台,我们能够在实验室条件下评估软波导在高功率输入下的性能表现。

    在高功率测试中,我们发现了现有WR1500软波导设计的一些缺陷。在高功率条件下,现有国产软波导的电热损耗远高于预期。这种损耗主要源于软波导结构设计,这些问题显著影响了传输线的整体性能和可靠性。为了解决这些问题,我们对软波导的结构进行了深入的优化设计。通过调整波导的结构参数以及选择更适合的材料,显著提升了WR1500软波导的电热性能,使其达到了国际标准。

    在材料选择方面,我们特别关注了铜和铝的性能表现。测试结果表明,虽然铜的电导率较高,但由于其质地较软和自重较大,在加工过程中容易发生变形。这种变形增大了内壁的不均匀性,从而增大了电磁损耗。相比之下,铝材料尽管电导率略低于铜,但由于其较高的硬度和较轻的重量,能够在加工及长期使用过程中保持形状稳定,从而显著减少了由于变形带来的电磁损耗,在高功率测试中表现出了更优异的稳定性和可靠性。

    在后续研究中,我们计划进一步研究应力形变对铜制WR1500软波导的影响,并寻找可能的改进措施。此外,我们还将改进现有的加工工艺,力图在进一步减小电热损耗的基础上,增强软波导的结构强度和耐用性,从而为高功率射频系统的稳定运行提供更为可靠的组件选择。

  • 图  1  射流区域示意图

    Figure  1.  Jet area diagram

    图  2  射流空间节点划分图

    Figure  2.  Jet space node division diagram

    图  3  SCAPS-1实验参数图

    Figure  3.  Experimental parameter diagram of SCAPS-1

    图  4  SCAPS-2实验参数图

    Figure  4.  Experimental parameter diagram of SCAPS-2

    图  5  气溶胶净化因子计算结果比对

    Figure  5.  Comparison with DF experimental data and model calculation results

    图  6  水力学模型计算结果

    Figure  6.  Hydraulic model calculation results

    表  1  实验条件及结果

    Table  1.   Test conditions and results

    experiment mass flux/
    (kg·m−2·s−1)
    Weber
    number
    steam
    fraction/%
    submerged
    depth/m
    material mass median
    diameter/μm
    decontamination
    factor
    SCAPS-1 217 3.5×106 64 0.7 TiO2 0.5 215±62.4
    SCAPS-2 120 1.8×106 64 0.7 TiO2 0.5 58±16.8
    RCA2[6] 95 2.9×105 0 0.5 Ni 1.3 28±12.3
    下载: 导出CSV
  • [1] Dong Shihang, Yang Jun. Overview of the experimental studies and numerical simulations on the filtered containment venting systems with wet scrubbers[J]. Annals of Nuclear Energy, 2019, 132: 461-485. doi: 10.1016/j.anucene.2019.04.042
    [2] Berzal M E, Marcos C, Swiderska-kowalczyk M, et al. State-of-the-art review on fission products’ aerosol pool scrubbing under severe accident conditions[R]. Luxembourg: European Commission, 1995.
    [3] Owczarski P C, Burk K W. SPARC-90: a code for calculating fission product capture in suppression pools[R]. Washington: Nuclear Regulatory Commission, 1991.
    [4] Ramsdale S A, Guentay S, Friederichs H G. BUSCA-JUN91 reference manual[R]. Villigen: Paul Scherrer Institute (PSI), 1995.
    [5] 周彦, 佟立丽, 曹学武. 低流速排放下气溶胶水洗模型[J]. 原子能科学技术, 2020, 54(5):923-929 doi: 10.7538/yzk.2019.youxian.0405

    Zhou Yan, Tong Lili, Cao Xuewu. Aerosol pool scrubbing model under low velocity discharge[J]. Atomic Energy Science and Technology, 2020, 54(5): 923-929 doi: 10.7538/yzk.2019.youxian.0405
    [6] Herranz L E, Lopez C, Penalva J. Investigation on jet scrubbing in nuclear reactor accidents: from experimental data to an empirical correlation[J]. Progress in Nuclear Energy, 2018, 107: 72-82. doi: 10.1016/j.pnucene.2018.04.019
    [7] Lind T, Dehbi A, Güntay S. Aerosol retention in the flooded steam generator bundle during SGTR[J]. Nuclear Engineering and Design, 2011, 241(1): 357-365. doi: 10.1016/j.nucengdes.2010.10.025
    [8] 李玉祥, 吴艳, 曹学武. 不同流型下气溶胶水洗效果影响因素实验研究[J]. 原子能科学技术, 2023, 57(2):276-283 doi: 10.7538/yzk.2022.youxian.0149

    Li Yuxiang, Wu Yan, Cao Xuewu. Experimental study on influencing factor of aerosol pool scrubbing under different flow regimes[J]. Atomic Energy Science and Technology, 2023, 57(2): 276-283 doi: 10.7538/yzk.2022.youxian.0149
    [9] Berna C, Escrivá A, Muñoz-Cobo J L, et al. Enhancement of the SPARC90 code to pool scrubbing events under jet injection regime[J]. Nuclear Engineering and Design, 2016, 300: 563-577. doi: 10.1016/j.nucengdes.2016.02.027
    [10] Allelein J H, Neu K, Skorek T, et al. European validation of the integral code ASTEC (EVITA)[R]. Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) mbH, 2004.
    [11] Hoefele E O, Brimacombe J K. Flow regimes in submerged gas injection[J]. Metallurgical Transactions B, 1979, 10(4): 631-648. doi: 10.1007/BF02662566
    [12] Berna C, Escrivá A, Muñoz-Cobo J L, et al. Review of droplet entrainment in annular flow: characterization of the entrained droplets[J]. Progress in Nuclear Energy, 2015, 79: 64-86. doi: 10.1016/j.pnucene.2014.11.011
    [13] Ghiaasiaan S M, Yao G F. A theoretical model for deposition of aerosols in rising spherical bubbles due to diffusion, convection, and inertia[J]. Aerosol Science and Technology, 1997, 26(2): 141-153. doi: 10.1080/02786829708965420
    [14] Slinn W G N. Precipitation scavenging: some problems, approximate solutions and suggestions for future research[R]. Richland: Battelle Pacific Northwest Laboratory, 1974.
    [15] Flagan R C, Seinfeld J H. Fundamentals of air pollution engineering[M]. Englewood Cliffs: Prentice-Hall, 1988.
    [16] Kaneko I, Fukasawa M, Naito M, et al. Experiment study on aerosol removal effect by pool scrubbing[C]//22nd DOE/NRC Nuclear Air Cleaning and Treatment Conference. 1992.
    [17] Li Yuxiang, Tong Lili, Cao Xuewu. Experimental study on influencing factors of aerosol retention by pool scrubbing[J]. Frontiers in Energy Research, 2021, 9: 675841. doi: 10.3389/fenrg.2021.675841
    [18] Moffat R J. Describing the uncertainties in experimental results[J]. Experimental Thermal and Fluid Science, 1988, 1(1): 3-17. doi: 10.1016/0894-1777(88)90043-X
  • 加载中
图(6) / 表(1)
计量
  • 文章访问数:  183
  • HTML全文浏览量:  66
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-05-08
  • 修回日期:  2024-07-23
  • 录用日期:  2024-09-24
  • 网络出版日期:  2024-12-20
  • 刊出日期:  2025-02-15

目录

/

返回文章
返回