核磁共振测速技术在燃料组件冷却剂通道流场实验研究的应用

王子瑞, 张科, 段敬添, 武俊梅, 兰治科

王子瑞, 张科, 段敬添, 等. 核磁共振测速技术在燃料组件冷却剂通道流场实验研究的应用[J]. 实验流体力学, doi: 10.11729/syltlx20230167.
引用本文: 王子瑞, 张科, 段敬添, 等. 核磁共振测速技术在燃料组件冷却剂通道流场实验研究的应用[J]. 实验流体力学, doi: 10.11729/syltlx20230167.
WANG Z R, ZHANG K, DUAN J T, et al. Application of MRV in experimental investigation of flow field characteristics in a coolant channel within a fuel assembly[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20230167.
Citation: WANG Z R, ZHANG K, DUAN J T, et al. Application of MRV in experimental investigation of flow field characteristics in a coolant channel within a fuel assembly[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20230167.

核磁共振测速技术在燃料组件冷却剂通道流场实验研究的应用

基金项目: 中核集团领创科研项目(J202210014–09)
详细信息
    作者简介:

    王子瑞: (1998—),男,陕西安康人,硕士研究生。研究方向:MRV流动测试技术。E-mail:15309153589@163.com

    通讯作者:

    张科: E-mail:kezhang@mail.xjtu.edu.cn

  • 中图分类号: TL339

Application of MRV in experimental investigation of flow field characteristics in a coolant channel within a fuel assembly

  • 摘要:

    核磁共振成像测速技术(Magnetic Resonance Velocimetry, MRV)作为复杂结构内三维全流场测量的前沿技术,是流场可视化研究、流动换热机理分析和数值模型验证的重要手段。近年来,MRV技术也被应用于核反应堆堆芯燃料组件冷却剂通道流场实验研究中,为堆芯热工水力性能研究和安全分析提供基础数据。本文整理了MRV技术在堆芯燃料组件冷却剂通道流场实验中的研究进展,在雷诺数12500和25000工况下,利用MRV技术对2 × 2简单棒束燃料组件冷却剂通道内的全空间三维流场进行了测量,并与文献中相同模型的PIV测量结果进行了定量对比,研究影响MRV测量精度的关键参数。研究结果表明:MRV技术可以精确解析棒束通道全空间三维流场,本文棒束通道内无量纲轴向速度的MRV测量结果与文献中PIV测量结果的最大误差为5.05%。

    Abstract:

    Magnetic Resonance Velocimetry (MRV) is becoming a crucial experimental technology for flow field visualization, flow heat transfer mechanism investigating, and numerical model validation. In recent years, the MRV technology has been employed in experimental investigations of coolant channel flow fields in nuclear reactor core fuel assemblies, providing fundamental data for core thermal hydraulic characteristics research and safety analysis. A literature overview of the research progress using MRV technology in experimental studies on the flow field of reactor core fuel assembly channels is presented in this paper. The three-dimensional coolant flow field within a 2 × 2 simplified rod bundle channel was measured using MRV, and quantitative comparisons were conducted with PIV measurement results from previous literature to investigate the key parameters influencing the measurement accuracy of MRV. In this study, MRV tests were performed under two conditions of Re = 12500 and 25000, respectively. The results indicate that the axial velocity from MRV within the investigated channel has a maximum discrepancy of 5.05% compared to the PIV data reported in the literature.

  • 可移动小堆具有小型紧凑、功率密度大、机动灵活、可靠性高等优点,是目前核推进技术领域的研究前沿。堆芯的高效能量转换是实现可移动小堆预定功能目标的关键,而堆芯燃料组件构型及冷却剂通道内的流动、换热性能则是堆芯热工性能的决定因素。因此,针对燃料组件内冷却剂流动特性开展相关研究至关重要。

    将易裂变材料以小颗粒形式弥散于辐照稳定性强的非裂变基体中构成燃料芯体,与外围金属包壳冶金结合制成弥散型燃料元件,将元件按照一定方式排列即可构成核燃料组件。目前,针对大型核电站压水堆采用的棒状燃料组件的研究和应用较为成熟,但在小堆紧凑性要求下,板状和环形燃料组件在小堆堆芯方案中更具优势。板状燃料元件发热面积大、冷却剂通道宽度小(窄缝通道宽度通常为1~3 mm),容易导致冷却剂流量分布不均、板内温度梯度过大,使板结构起泡、变形甚至破裂,进而增大通道阻力,使冷却剂流动产生局部不稳定性,引发燃料板流致振动。研究者还提出了一种带中心孔的双面冷却环形燃料芯块设计方案,目的是为了显著降低芯块最高温度、提高冷却均匀性,进而提高堆芯安全性和热工性能。这种结构涉及了内外双通道流动,流动特性更加复杂。为优化燃料组件设计方案,提高燃料组件的热工水力性能与安全性,需要进一步揭示其三维流场分布、边界层效应、流量分配和几何参数影响。

    燃料组件内冷却剂的三维流动特性是决定堆芯传热效率和安全性的重要因素。堆芯内冷却剂通道结构复杂、尺寸微小,给相关实验技术研究带来了巨大的挑战。近年来,有研究者尝试将核磁共振成像技术应用于微小复杂通道内的三维流场测量,即核磁共振成像测速(Magnetic Resonance Velocimetry, MRV)[1]。与激光多普勒测速(LDV)[2]、粒子成像测速(PIV)[3]等传统流场测量技术相比,核磁共振成像测速技术无粒子浸入、不受光学条件限制,能够快速解析微小复杂结构内部全空间三维三分量(3D3C)速度分布,可以精确捕捉涡旋的三维空间结构,是精确可靠的微小复杂通道内全场流动特性研究前沿实验技术[4]。此外,研究者也开始将MRV技术应用于气、固、液等多相流动的复杂流场测量研究,MRV技术应用潜力巨大。

    基于核磁共振成像原理进行流动测量并不是一个新概念,它已被广泛应用于医学领域(如生物活体中的层流研究等),但在工程湍流方面的测量研究则相对较少[5]。2003年,Elkins等[6]首先将MRV应用于高速湍流测量,认为MRV在湍流测量方面具有较高可靠性和准确度。之后,Ryan等[7]利用MRV测量了垂直旋转风涡轮的三维流场,Benson等[8]借助MRV研究了涡轮叶片气膜冷却形成的三维混合流,Freudenhammer等[9]利用MRV研究了气缸盖几何形状变化对进气体积流量的影响,Borup等[10]借助MRV对粒子分离器内的三维流动特性进行了研究,Hoffman等[11]利用MRV得到了离心叶轮间的流场分布。

    目前,MRV技术在燃料组件通道空间流场测量方面的应用还不多。本文整理了MRV技术在堆芯燃料组件通道流场实验研究中的研究进展,利用MRV技术对2 × 2简单棒束燃料组件冷却剂通道(后文或简称“棒束通道”)内全空间三维流场进行测量,并与公开文献中相同模型的PIV测量结果进行定量对比,分析影响MRV测量精度的关键参数。

    部分研究者已将MRV技术应用于燃料组件冷却剂通道内流场测量实验研究。Piro等[12]以激光烧结3D打印技术制造了以CANDU燃料棒束通道为原型并按2/3比例缩放简化后的8单元棒束通道模型,采用MRV技术对其内部流场特性开展了实验研究,得到的不同轴向位置的轴向速度与CFD模拟结果吻合良好。Piro等[13]还制作了等比例37单元棒束通道模型,在雷诺数12500工况下得到的MRV测量结果与CFD模拟结果基本一致;将燃料棒束截面细分为4个子通道环,通过分析不同轴向位置截面的径向分速度和角速度,研究了4个子通道间的流量掺混能力。此外,该研究还意外地验证了MRV技术能够检出棒束通道中卡住的材料碎片,并能够量化这些碎片对局部流动的影响。综上所述,MRV技术适用于微小复杂结构的测量,使得多单元棒束通道内的全空间三维流场测量成为可能。

    以上述研究作为基础,Bruschewski等[14]将燃料组件通道压力管内径扩大了6%,得到了发生蠕变老化的棒束通道模型,利用MRV技术测量了雷诺数Re = 11550工况下的三维流场数据,得到的轴向速度与采用LDV技术测得的局部轴向速度十分吻合。此外,Bruschewski等还指出,在径向蠕变的棒束通道中,通过旁路区域(该旁路区域的相对几何面积为25%)的流量在下游端板附近稳定在总流量的35%左右,且旁路区域相对速度稳定在主流速度的1.9倍左右。Oliveira等[15]Re = 1936工况下,利用MRV技术分别测量了7 × 7棒束的法国压水堆常规模型、模拟失水事故(Loss Of Coolant Accident, LOCA)下堵塞率为90%的通道模型的流场特性,分析了“上游过渡区–堵塞区–下游过渡区”的流量重分配现象,指出从堵塞通道偏离出去的流量主要由几何堵塞率决定,同时证明了带搅混翼的定位格架能够很好地弱化堵塞区对流动的影响。基于上述研究,Oliveira等[16]Re = 1936、4671、7223和9599工况下测试了6个棒束通道模型(分为堵塞区共面和非共面情况),共面堵塞的4个模型最终堵塞率为61%、80%、90%、90%,对应堵塞长度为100、100、100、240 mm(前3个模型为1组,对比讨论堵塞率对流动的影响;后2个模型为1组,对比讨论堵塞长度对流动的影响),非共面堵塞的2个模型最终堵塞率都为80%。文献指出:堵塞率对下游流场影响更大(具体而言,堵塞率越大,堵塞区内轴向速度越低,其下游的流量重分配越强烈);流量掺混程度对堵塞长度不敏感;随着流量增大,堵塞区下游过渡区流动介质的轴向速度下降更为剧烈。在相同工况下,非共面堵塞情况下发生的流量重分配强度较小,而下游堵塞区可能会影响上游堵塞区引起的流量重分配效应。此外,在不同堵塞区几何形状、不同实验流率工况下,上下游2个带搅混翼的定位格架总能有效地将堵塞效应限制在两者之间的区域(即能使流动通过格架后均匀化)。由此可见,MRV技术能够用于研究复杂结构内几何形状的微小变化对局部与全空间三维流场的影响。

    采用MRV技术测量高雷诺数湍流流场的研究也有很大进展。在新近研究中,Bruschewski等[17]采用MRV技术对Re = 20300工况下5 × 5方形排布压水堆燃料组件模型(通道以透明有机玻璃制作;燃料棒束模型和定位格架按照实际压水堆燃料组件放大1.26倍,以聚甲醛树脂制作;格架采用了一种法马通设计构型)内部的三维速度场进行了测量分析,测量结果与CFD模拟结果十分吻合。此外,在丘陵状通道模型中进行了二维MRV湍流测量,获得了通道中第12个和第13个“山丘”之间的雷诺应力分布,除流动略有不对称及存在一些系统误差外,MRV可得到出色的数据质量。

    本文将MRV技术用于棒束通道全空间流场测量。目前,采用非MRV技术对棒束通道流场分布开展的实验研究相对较多,有一些公认的实验结果可作为MRV技术实验方法和可靠性验证基准。本文研究模型为2 × 2简单棒束燃料组件冷却剂通道模型,参考了Hosokawa等[18]的文献。将本文MRV测量结果与文献[18]PIV测量结果进行对比分析,研究影响MRV测量精度的系统参数。

    关于MRV的基本原理详见文献[4],本文仅作简要介绍。将原子核系统置于均匀强外磁场B0中(下文均默认B0方向沿z轴正向),原子核受磁力矩作用绕B0方向发生拉莫尔进动(进动频率ω = γB0,其中γ为原子核的旋磁比,仅由原子核种类决定),并在宏观上产生纵向磁化。此时若再向系统施加频率恰为ω的射频(RF)磁场,系统将对射频脉冲能量产生强烈的吸收,发生核磁共振。脉冲结束后,系统趋于恢复至最低能态,因此纵向磁化趋于恢复、横向磁化矢量在Oxy平面内绕B0方向旋转、变化。由法拉第电磁感应定律可知,固定位置的接收线圈内的磁通量不断发生变化,线圈内将会产生感应电流,此感应电流即为核磁共振信号MR。

    在磁共振扫描中,线圈接收的是全场氢质子的MR信号,其位置信息是利用各方向的编码来区分的。例如,向质子系统施加一个平行于z轴的梯度磁场GSS,则z方向上不同位置的核的拉莫尔进动频率呈线性变化,则经过z方向上某位置的Oxy截面内的核仅能由一种频率的脉冲所激励,如此便实现了z方向上不同位置截面的选择。xy方向的位置选择也同样通过施加梯度磁场来实现。如图1所示,在1个重复时间(图中的TR)内施加RF激励、层面选择梯度磁场GSS、相位编码梯度磁场GPE、频率编码梯度场GFE,可接收到一组MR信号,重复此过程,即可由二维傅里叶变换重建信号灰度图[19]

    图  1  单方向速度编码扫描序列
    Fig.  1  MRV sequences for velocity encoding

    下面简述MRV中的速度编码方法。以x方向为例,将沿x方向施加的梯度磁场记为$ \mathop G\limits^ \to \left(t\right) $,则在该磁场中以速度$ \mathop u\limits^ \to $运动的当地氢核拉莫尔进动频率$ {\omega }_{L} $可表示为:

    $$ \begin{array}{c}{\omega }_{L}\left( \mathop {x}\limits^ \to,t\right)=\gamma \left({B}_{0} + \mathop G\limits^ \to \left(t\right) \cdot \mathop x\limits^ \to \left(t\right)\right)\end{array} $$ (1)

    式中,$ \mathop x\limits^ \to \left(t\right)=\mathop {{x}_{0}}\limits^ \to + \mathop u\limits^ \to(t-{t}_{0}) $,$ \mathop {{x}_{0}}\limits^ \to $表示质子在初始时刻$ {t}_{0} $的位置,代入式(1)可得:

    $$ {\text{Φ}} \left(\mathop {x}\limits^ \to,\tau \right)={{\text{Φ}} }_{0} + \gamma \mathop {{x}_{0}}\limits^ \to{\int }_{0}^{\tau }\mathop {G}\limits^ \to \left(t\right)\mathrm{d}t + \gamma \mathop {u}\limits^ \to{\int }_{0}^{\tau }\mathop {G}\limits^ \to \left(t\right)t\mathrm{d}t + O\left({t}^{n}\right) $$ (2)

    式中等号右侧前3项分别反映初始相位、梯度磁场对静止氢核进动相位的影响、氢核运动速度$ \mathop {u}\limits^ \to $对进动相位的影响。如图1所示,在$ {t}_{2}~{t}_{3} $期间,对运动氢核施加正负双向梯度磁场$ {G}_{\mathrm{k}} $(即图1中绿色矩形表示的流动编码中所用到的梯度磁场),可消除仅施加单向梯度磁场带来的氢核进动相位变化,则式(2)等号右侧第3项可表示为:

    $$ \begin{array}{c}{{\text{Φ}}}_{1}=-\gamma \overline{u}{N}_{1}\end{array} $$ (3)

    式中,$ {N}_{1}={G}_{\mathrm{k}}{\text{Δ}}{t}^{2} $。因此,采用MRV技术进行测量时,在$ {t}_{1}~{t}_{2} $期间施加层面选择梯度$ {G}_{\mathrm{s}\mathrm{s}} $,激发一定层面后,先进行一次流动编码,即可使式(2)等号右侧第2项为0。初始相位$ {{\text{Φ}}}_{0} $则通过重复编码消除,如图1所示,在$ {t}_{7}~{t}_{12} $期间施加第2次重复编码,两次编码的初始相位$ {{\text{Φ}} }_{0} $相同,两次编码后的相位差为$ \mathrm{{\text{Δ}} }{\text{Φ}} ={{\text{Φ}} }_{1}^{\left(2\right)}-{{\text{Φ}} }_{1}^{\left(1\right)} $,整理得:

    $$ \begin{array}{c}{\text{Δ}} {\text{Φ}} =\gamma \overline{u}\left({N}_{1}^{\left(1\right)}-{N}_{1}^{\left(2\right)}\right)=\gamma \overline{u}{\text{Δ}} {N}_{1}\end{array} $$ (4)

    根据式(4),当地x方向氢核宏观运动的速度$\overline{u} $可表示为:

    $$ \begin{array}{c}\overline{u}={E}_{\mathrm{v}}{\text{Δ}} {\text{Φ}} /\pi \end{array} $$ (5)

    式中,速度编码(velocity encoding)$ {E}_{\mathrm{v}}=\mathrm{\pi }/\gamma \mathrm{{\text{Δ}} }{N}_{1} $,反映了流动编码磁场强度引起的相位变化程度:由于$ \mathrm{{\text{Δ}} }{\text{Φ}} $介于$ -\mathrm{\pi }~\mathrm{\pi } $之间,因此式(5)中的速度编码表征了该次测量中出现的最大流速值。在MRV实验时,系统根据事先设定的速度编码(即所测最大流速值)计算所需的流动编码磁场强度。

    通过上述系列编码,接收线圈中获得MR信号,再经三维傅里叶变换,获得不同位置MR信号中的幅值与相位信息,并存储于核磁共振成像(Magnetic Resonance Imaging, MRI)设备。提取MR信号中的相位信息,结合式(5),可计算得到扫描区域内各体素(MRV扫描的最小空间单元)中氢核沿x方向的速度分布$ \overline{u} $(x, y, z)。MRV信号为体素内所有氢核共同产生,因此所测速度即为体素内的平均速度。将相位信息转换为速度信息,对空间流场进行分析的过程由自主开发的Matlab程序完成。

    本文实验在西安交通大学生物医学与健康工程研究院的核磁共振成像仪(西门子MAGNETOM Spectra 3.0T)上进行,其先进的Tim 4G技术可提供优异的信噪比和均匀度,减少误差和伪影。考虑到该设备的实际安装情况、舱内参数和使用要求,设计了如图2所示的棒束燃料组件冷却剂通道MRV全场流动测量实验平台。为防止核磁共振成像仪中的强磁体吸引器材或零件,影响信号接收甚至危及人身安全,须严格避免核磁共振仪周围5 m范围内出现铁磁性物质。因此,将核磁共振成像仪和棒束燃料组件冷却剂通道实验段(非金属材料制成)置于专用放射室内,其余设备均置于放射室外,研究人员在放射室外的计算机上进行实验操作。

    图  2  棒束燃料组件冷却剂通道MRV全场流动测量实验平台
    Fig.  2  MRV full-field flow measurement platform for coolant channel of rod bundle

    实验中,冷却剂从冷却剂箱出发,沿塑料软管流经离心泵、调节阀、涡轮流量计等,流入实验段中的沉降室,再从沉降室进入燃料组件冷却剂通道,从冷却剂通道中流出后,经软管流回冷却剂箱,构成冷却剂循环系统。冷却剂通道内的流量大小通过调节阀和流量计进行调节和控制,可满足在Re = 0~4 × 104之间工作的需求。

    置于核磁共振成像仪中的实验段以有机玻璃制成,部件之间通过法兰和尼龙螺栓连接,连接管道为塑料软管。离心泵中的金属部件以304不锈钢制成,可防止生锈引入杂质,影响测试结果。

    核磁共振成像仪舱室洞口截面直径仅约80 cm,实验段尺寸不能过大,同时需兼顾测量空间分辨率要求,实验模型尺寸不能过小,以保证能够在冷却剂流场内布置足够的数据采集点。简单棒束燃料组件冷却剂通道流场测量实验的目的是验证MRV测量的可靠性和精度,因此,将棒束通道设计为简单结构(主要参数来自Hosokawa等[18]发表的PIV测量文献,本文MRV测量结果将与该文献PIV测量结果进行对比),如图3所示。棒束通道横截面(垂直来流方向的截面)尺寸为50 mm × 50 mm,通道内对称布置2 × 2个单元棒,棒直径为20 mm。棒之间的冷却剂通道最窄处的宽度为5.0 mm,棒和通道壁面之间最窄处的宽度为2.5 mm。

    图  3  棒束燃料组件冷却剂通道模型结构
    Fig.  3  Diagram of coolant channel in rod bundle

    棒束燃料组件冷却剂通道实验段(图4)由沉降室、发展段、测试段和出口段等组成,置于核磁共振成像仪舱室内。其中,沉降室为正方体腔室(100 mm × 100 mm × 100 mm),能够改善来流稳定性并沉降可能出现的杂质。沉降室后连接棒束燃料组件冷却剂通道,通道内的燃料棒束由3个定位格架支撑。定位格架通过12个螺栓连接孔和通道以法兰连接。单个棒束通道尺寸为50 mm × 50 mm,通道4个壁面布置宽2.0 mm、高2.5 mm或1.5 mm的突起,用于支撑固定圆柱形棒束,如图5所示(图中以黑色实心圆形表示棒束)。测试区域位于第2个定位格架下游,以保证来流在通道内充分发展。实验模型(指图4中的所有部件,包括插于实验通道内部、图中未显示出的棒束)均采用有机玻璃加工制成。

    图  4  棒束燃料组件冷却剂通道实验段
    Fig.  4  Coolant channel model of rod bundle
    图  5  定位格架右视图及主要几何参数
    Fig.  5  Right elevation of spacer grid and geometric parameters

    在雷诺数Re = 12500、25000工况下进行了MRV测试(Re =$ \rho {u}_{\mathrm{b}}{D}_{\mathrm{h}}/\mu $,其中,通道水力直径$ {D}_{\mathrm{h}} $= 11 mm;ρ为冷却剂密度,单位kg/m3ub为冷却剂主流流速,单位m/s;μ为冷却剂动力黏度,单位Pa·s)。实验中,对每个工况设定x、${\textit{y}} $、${\textit{z}} $等3个方向的速度编码值相同,即Ev, x = Ev, y = Ev, z。对于Re = 12500工况,Ev, x = Ev, y = Ev, z = 4.0 m/s,对于Re = 25000工况,Ev, x = Ev, y = Ev, z = 4.5 m/s。经过多次实验对比,综合考虑所测空间大小及采集效率,在空间分辨率为1.0 mm × 0.4 mm × 0.4 mm、0.7 mm × 0.7 mm × 0.7 mm和 1.0 mm × 0.7 mm × 0.7 mm条件下进行了测试。在第2个定位格架下游600 mm处、约200.0 mm × 230.2 mm × 224.0 mm空间区域中,单次流动扫描在8 min左右测量了超过8 × 106个空间点3个方向的速度场。为提高信噪比,在水中加入CuSO4·5H2O粉末,制成摩尔浓度0.06 mol/L的CuSO4溶液作为冷却剂。对每个工况至少进行了6次流动扫描,以减小随机误差的影响。

    Re = 12500时,在6次流动扫描中测得的平均主流流速标准差为0.005 m/s。主流总平均流速$ \overline{{u}_{\mathrm{b}}} $= 0.9276 m/s,与基于流量计示数计算所得平均流速仅相差4.1%。MRV测得的速度不确定度$ {\stackrel{~}{u}}_{i} $可通过既定方向上的速度编码值Ev和所采集信号的信噪比Rsn计算[6]

    $$ \begin{array}{c}\widetilde{{u}_{i}}=\left(\sqrt{2}/\mathrm{\pi }\right){E}_{\mathrm{v}}/{R}_{\mathrm{s}\mathrm{n}}\end{array} $$ (6)

    式中,i表示x、${\textit{y}} $或${\textit{z}} $方向。多次实验结果对比表明,信噪比和空间分辨率、扫描时间及扫描区流速大小有关:空间分辨率越高,信噪比越低;扫描时间延长,信噪比相对升高;流速低的区域,MR信号较弱,则该区域信噪比较低,测量不确定度较大。

    在各次MRV实验中,Re = 12500时信噪比最高,空间分辨率为1.0 mm × 0.7 mm × 0.7 mm,设定x方向速度编码Ev, x = 4.0 m/s,由MRV结果计算得出Rsn, x = 49.69,将Ev, xRsn, x代入式(6),得到本次实验x方向速度不确定度$ \widetilde{{u}_{x}} $= 0.0362 m/s,仅为通道中心最大速度的2.4%,和文献[14]得到的精度处于同一水平(该文献中得到的速度不确定度约为子通道中心处通道最大速度的3.0%),验证了MRV实验结果的可靠性。

    Re = 12500和25000条件下,对棒束通道进行多次不同扫描序列和参数设置下的MRV速度测量实验。图6为通过实验测得的棒束通道扫描区域内的全空间轴向速度分布云图,空间分辨率分别为1.0 mm × 0.7 mm × 0.7 mm和1.0 mm × 0.4 mm × 0.4 mm,图中H为冷却剂通道横截面边长。对比可见,空间分辨率减小会导致信噪比降低和速度云图平滑性下降。从图中可以看到:云图中棒束所在位置未显示速度数据;单元棒之间区域冷却剂速度较高,而单元棒和壁面之间冷却剂速度较低,最高流速发生于通道中心。可以看出,MRV全空间速度测量使复杂结构流动全场可视化解析成为可能,是复杂结构流场全空间分布特征研究的有力工具。

    图  6  棒束通道扫描区域内轴向速度三维分布云图
    Fig.  6  3D contour of axial velocity distribution in the scanning area

    在全空间三维速度分布云图中,任意截取二维截面、一维线和单点数据,可以提取二维截面、一维线及单点三维速度分布数据,研究空间细节流场分布特征。图7图6(a)扫描区域最上游(x/H = 0处)横截面轴向速度分布云图,可以更清晰地看到:通道截面中心轴向速度最高;速度分布沿z方向对称性较好,沿y方向对称性较差,y/H = 0.5处壁面附近冷却剂速度相对y/H = −0.5处壁面附近较小。分析发现,沿y方向速度分布的非对称性可能是由棒束变形产生的几何结构非对称所引起,更换棒束后,该非对称性发生改变。MRI图像表明,沿流向不同位置横截面上的棒束圆心位置有所偏移,证实棒束沿通道发生了变形。

    图  7  Re = 12500时x/H = 0处截面轴向速度分布二维云图
    Fig.  7  2D contour of axial velocity distribution in cross section at x/H = 0 for Re = 12500

    图8图6(a)通道高度方向中心位置(z/H = 0处)的Oxy截面速度分布云图和速度矢量图。同样可以看出:速度分布沿z方向对称性较好,沿y方向对称性较差;yz方向无速度分量,或相对主流可忽略不计。

    图  8  Re = 12500时z/H = 0处Oxy截面速度云图和速度矢量图
    Fig.  8  2D contour of Oxy velocity distribution and velocity vector diagram in coronal plane at z/H = 0 for Re = 12500

    图9为扫描区域最上游(x/H = 0处)横截面轴向速度分布云图(Re = 25000)。从图中可见,通道内速度沿y方向存在非对称性,y轴正向一侧轴向速度比y轴负向一侧更高,和Re = 12500时正好相反。在各次实验中,将棒束更换位置重新安装,棒束变形引起的几何结构非对称性发生了改变,因此,不同雷诺数下流场非对称性的变化可能是由棒束变形所引起的。由此可见,棒束通道中细微的几何变形,就可能导致通道内流场分布的非对称性,而MRV技术能够清晰解析几何结构改变引起的流场变化特征。

    图  9  Re = 25000时x/H = 0处截面轴向速度分布二维云图
    Fig.  9  2D contour of axial velocity distribution in cross section at x/H = 0 for Re = 25000

    图10Re = 25000时z/H = 0处的Oxy截面速度分布云图和速度矢量图。同样地,当Re = 25000时,在下游远离定位格架的棒束通道中,流场中未出现显著的二次流,而沿y方向与Re = 12500时相反的不对称性进一步说明流场的非对称性可能由棒束变形所致。

    图  10  Re = 25000时z/H = 0处的Oxy截面速度云图和速度矢量图
    Fig.  10  2D contour of Oxy velocity distribution and velocity vector in coronal plane at z/H = 0 for Re = 25000

    本文测量的棒束通道在几何参数上与文献[18]一致。将MRV测量结果和文献PIV测量结果进行对比分析,以验证MRV技术的精确性,并讨论MRV系统参数设置等对测量结果的影响。

    图11(a)为Re = 25000时棒束通道垂直来流的横截面左上角单元棒周围区域速度分布云图。图11(b)为文献中PIV测量结果。从图中可以看到,MRV测得的横截面速度分布特征和速度大小与PIV测得的基本一致,但该次MRV测量结果的信噪比和空间分辨率较低。图11(c)为Re = 12500时的MRV测量结果,其信噪比图11(a)明显提升,除因雷诺数差异造成的绝对速度大小差异外,速度分布特征和PIV测量结果更加接近。由于文献[18]未展示Re = 12500时同一位置的速度分布云图,本文未定性对比该工况下的MRV与PIV测量结果。下文将对不同雷诺数下特定区域的速度开展MRV测量结果和PIV测量结果的定量对比。

    图  11  Re = 25000和12500时横截面左上角1/4区域内MRV轴向速度分布云图与PIV测量结果对比
    Fig.  11  Contrast for contour of axial velocity distribution in 1/4 cross section at x/H = 0 for Re = 25000 and 12500

    x/H = 0处的横截面上,选取3条路径:Line A(红色)、Line B(蓝色)、Line C(绿色),如图12(a)所示。沿3条路径提取轴向速度数据,得到如图12(b)所示的速度分布线图,并与文献中给出的PIV测量结果对比(图中散点数据为文献中的PIV测量结果,线数据为MRV测量结果;为与PIV数据进行区分,将MRV点数据拟合为线数据)。除信噪比较低外,MRV测量结果的整体趋势与PIV测量结果吻合良好。

    图  12  Re = 25000时通道不同位置速度分布线图
    Fig.  12  Line map of axial velocity distribution at different positions for Re = 25000

    图7x/H = 0处横截面内轴向速度分布云图中,提取z/H = 0和z/H = −0.25处的Line B(蓝色)、Line D(黑色)路径的轴向速度值并进行无量纲化,与文献中不同雷诺数下的PIV测量结果进行对比,如图13所示。无量纲轴向速度$ \overline{{U}_{0}} $定义为$ \overline{{U}_{0}}=\overline{u}/\overline{{u}_{0}} $,其中$ \overline{{u}_{0}} $为测量区域入口截面平均轴向速度。从图中可以看到,针对相同几何模型,通过MRV测量得到的无量纲速度分布与文献PIV测量结果完全吻合。图13给出的MRV测量结果与文献中PIV测量结果最大误差约为5.05%。

    图  13  不同位置无量纲轴向速度分布线图
    Fig.  13  Non-dimensional axial velocity linemap at different positions

    以上测量和对比结果说明,目前掌握的MRV技术在简单棒束通道内冷却剂全场速度测量方面结果精确可靠,为后续更复杂的燃料组件冷却剂通道全空间流场的测量分析打下了基础。

    在空间分辨率影响方面,分别研究了3种不同空间分辨率(高分辨率:1.0 mm × 0.4 mm × 0.4 mm;中分辨率:0.7 mm × 0.7 mm × 0.7 mm;低分辨率:1.0 mm × 0.7 mm × 0.7 mm)对棒束通道速度分布测量的影响。图14为不同空间分辨率下的棒束通道横截面速度分布对比,从左至右空间分辨率依次降低。从图中可以看到:随着空间分辨降低,MRV速度测量结果的信噪比显著增大。因为随着空间分辨率降低,磁扫描单位体素内所含的流体体积变大,MR信号的信噪比增大,导致速度测量的不确定度降低。在最低分辨率(1.0 mm × 0.7 mm × 0.7 mm)下,速度测量不确定度为0.0362 m/s,相对主流最大速度为2.4%,对应信噪比为49.69;在最高分辨率(1.0 mm × 0.4 mm × 0.4 mm)下,速度测量不确定度为0.1479 m/s,相对主流最大速度为4.9%,对应信噪比为12.18。从最低分辨率提高至最高分辨率,信噪比降低了75.5%。

    图  14  不同空间分辨率下棒束通道截面速度分布云图对比
    Fig.  14  Comparison of velocity distribution for different spatial resolutions

    x/H = 0处横截面内、Line B路径(z/H = 0)上的速度数据除以该横截面上平均流速,得到如图15所示的无量纲轴向速度沿y方向的分布曲线。从图中可以明显看到:随着空间分辨率降低,分布曲线更加光滑,说明速度测量信噪比提高、不确定度降低。影响速度测量不确定度的因素还有速度编码和信噪比。当空间分辨率降低、单位体素内所含的流体体积变大时,信号增强,信噪比增大,速度测量不确定度降低。速度编码值越小,速度测量不确定度越低。图14(a)的速度编码值为4.5 m/s,图14(b)和(c)的速度编码值为4.0 m/s。因此,图14中速度测量不确定度的降低同时受到了速度编码值降低的影响。图14(a)对应实验Re = 25000,图14(b)和(c)对应Re = 12500,实验雷诺数降低,引起流场中速度脉动降低,也是速度测量信噪比增大的一个不可忽视的原因。

    图  15  不同空间分辨率下z/H = 0处无量纲轴向速度沿y方向的分布曲线
    Fig.  15  Distribution of dimensionless velocity along y direction at z/h = 0 for different spatial resolutions

    在梯度磁场方向影响方面,图14(b)和(c)分别对应2种不同的梯度磁场方向。图14(b)测量截面选择沿通道主流方向,图14(c)测量截面选择垂直通道主流方向。从图中可以看出,不同截面选择即不同测量梯度磁场方向对测量结果影响不大。从图14(b)到(c),测量信噪比增大主要由空间分辨率的不确定度引起。但是,当测量截面垂直主流方向时,由于截面内流体所占区域较小,固体(棒束截面)区域较大,在固体边界会出现由于截断导致的边缘信号振荡(已在数据处理过程中将振荡部分剖去),在分辨率较高的图14(a)中该现象不明显。

    本文利用MRV技术实现了对2 × 2棒束通道内冷却剂全场三维流场信息的解析,并与公开文献PIV测量结果进行了定量对比。

    1)2 × 2简单棒束通道内冷却剂流场分布特征方面,在远离定位格架的通道下游区域中,流场中未出现显著的二次流。

    2)在相同几何结构和相同雷诺数下,由MRV测量得到的无量纲轴向速度分布与文献中相应的PIV测量结果最大误差为5.05%,验证了MRV在棒束通道内进行全空间三维三分量速度测量的准确度和可靠性。

    3)在影响MRV速度测量精度的关键系统参数方面:空间分辨率越小,信噪比越低;速度编码值越小,速度不确定度越低;MRV测量梯度磁场方向对测量结果影响不大。

    致谢:感谢西安交通大学生物医学与健康工程研究院健康与康复科学研究所王珏教授团队为本文研究工作提供MRI设备和实验指导。

  • 图  10   Re = 25000时z/H = 0处的Oxy截面速度云图和速度矢量图

    Fig.  10   2D contour of Oxy velocity distribution and velocity vector in coronal plane at z/H = 0 for Re = 25000

    图  1   单方向速度编码扫描序列

    Fig.  1   MRV sequences for velocity encoding

    图  2   棒束燃料组件冷却剂通道MRV全场流动测量实验平台

    Fig.  2   MRV full-field flow measurement platform for coolant channel of rod bundle

    图  3   棒束燃料组件冷却剂通道模型结构

    Fig.  3   Diagram of coolant channel in rod bundle

    图  4   棒束燃料组件冷却剂通道实验段

    Fig.  4   Coolant channel model of rod bundle

    图  5   定位格架右视图及主要几何参数

    Fig.  5   Right elevation of spacer grid and geometric parameters

    图  6   棒束通道扫描区域内轴向速度三维分布云图

    Fig.  6   3D contour of axial velocity distribution in the scanning area

    图  7   Re = 12500时x/H = 0处截面轴向速度分布二维云图

    Fig.  7   2D contour of axial velocity distribution in cross section at x/H = 0 for Re = 12500

    图  8   Re = 12500时z/H = 0处Oxy截面速度云图和速度矢量图

    Fig.  8   2D contour of Oxy velocity distribution and velocity vector diagram in coronal plane at z/H = 0 for Re = 12500

    图  9   Re = 25000时x/H = 0处截面轴向速度分布二维云图

    Fig.  9   2D contour of axial velocity distribution in cross section at x/H = 0 for Re = 25000

    图  11   Re = 25000和12500时横截面左上角1/4区域内MRV轴向速度分布云图与PIV测量结果对比

    Fig.  11   Contrast for contour of axial velocity distribution in 1/4 cross section at x/H = 0 for Re = 25000 and 12500

    图  12   Re = 25000时通道不同位置速度分布线图

    Fig.  12   Line map of axial velocity distribution at different positions for Re = 25000

    图  13   不同位置无量纲轴向速度分布线图

    Fig.  13   Non-dimensional axial velocity linemap at different positions

    图  14   不同空间分辨率下棒束通道截面速度分布云图对比

    Fig.  14   Comparison of velocity distribution for different spatial resolutions

    图  15   不同空间分辨率下z/H = 0处无量纲轴向速度沿y方向的分布曲线

    Fig.  15   Distribution of dimensionless velocity along y direction at z/h = 0 for different spatial resolutions

  • [1]

    MORAN P R. A flow velocity zeugmatographic interlace for NMR imaging in humans[J]. Magnetic Resonance Imaging, 1982, 1(4): 197–203. doi: 10.1016/0730-725x(82)90170-9

    [2]

    GEORGE W K Jr. An analysis of the laser-doppler velocimetry and its application to the measurement of turbulence[D]. Baltimore, Washington D. C. : Johns Hopkins University, 1971.

    [3]

    ADRIAN R J. Particle-imaging techniques for experimental fluid-mechanics[J]. Annual Review of Fluid Mechanics, 1991, 23: 261–304. doi: 10.1146/annurev.fl.23.010191.001401

    [4] 张科, 段敬添, 雷蒋, 等. 基于MRV的菱形肋柱冷却通道三维全流场分析[J]. 航空动力学报, 2023, 38(12): 2837–2847. DOI: 10.13224/j.cnki.jasp.20220051

    ZHANG K, DUAN J T, LEI J, et al. Analysis on 3D full flow field in a cooling channel with diamond pin fins based on MRV[J]. Journal of Aerospace Power, 2023, 38(12): 2837–2847. doi: 10.13224/j.cnki.jasp.20220051

    [5]

    ELKINS C J, ALLEY M T. Magnetic resonance velocimetry: applications of magnetic resonance imaging in the measurement of fluid motion[J]. Experiments in Fluids, 2007, 43(6): 823–858. doi: 10.1007/s00348-007-0383-2

    [6]

    ELKINS C J, MARKL M, PELC N, et al. 4D Magnetic resonance velocimetry for mean velocity measurements in complex turbulent flows[J]. Experiments in Fluids, 2003, 34(4): 494–503. doi: 10.1007/s00348-003-0587-z

    [7]

    RYAN K J, COLETTI F, ELKINS C J, et al. Three-dimensional flow field around and downstream of a subscale model rotating vertical axis wind turbine[J]. Experiments in Fluids, 2016, 57(3): No. 38. doi: 10.1007/s00348-016-2122-z

    [8]

    BENSON M J, ELKINS C J, EATON J K. Measurements of 3D velocity and scalar field for a film-cooled airfoil trailing edge[J]. Experiments in Fluids, 2011, 51(2): 443–455. doi: 10.1007/s00348-011-1062-x

    [9]

    FREUDENHAMMER D, PETERSON B, DING C P, et al. The influence of cylinder head geometry variations on the volumetric intake flow captured by magnetic resonance velocimetry[J]. SAE International Journal of Engines, 2015, 8(4): 1826–1836. doi: 10.4271/2015-01-1697

    [10]

    BORUP D D, ELKINS C J, EATON J K. Experimental analysis of a particle separator design with full-field 3D measurements[R] ASME Paper GT2019-90400, 2019. doi: 10.1115/GT2019-90400

    [11]

    HOFFMAN D W, VILLAFAÑE L, ELKINS C J, et al. Experimental study of flow inside a centrifugal fan using magnetic resonance velocimetry[J]. Journal of Engineering for Gas Turbines and Power, 2020, 142(4): 041019. doi: 10.1115/1.4045064

    [12]

    PIRO M H A, WASSERMANN F, GRUNDMANN S, et al. Progress in on-going experimental and computational fluid dynamic investigations within a CANDU fuel channel[J]. Nuclear Engineering and Design, 2016, 299: 184–200. doi: 10.1016/j.nucengdes.2015.07.009

    [13]

    PIRO M H A, WASSERMANN F, GRUNDMANN S, et al. Fluid flow investigations within a 37 element CANDU fuel bundle supported by magnetic resonance velocimetry and computational fluid dynamics[J]. International Journal of Heat and Fluid Flow, 2017, 66: 27–42. doi: 10.1016/j.ijheatfluidflow.2017.04.010

    [14]

    BRUSCHEWSKI M, PIRO M H A, TROPEA C, et al. Fluid flow in a diametrally expanded CANDU fuel channel – Part 1: experimental study[J]. Nuclear Engineering and Design, 2020, 357: 110371. doi: 10.1016/j.nucengdes.2019.110371

    [15]

    OLIVEIRA A V S, STEMMELEN D, LECLERC S, et al. Velocity field and flow redistribution in a ballooned 7 × 7 fuel bundle measured by magnetic resonance velocimetry[J]. Nuclear Engineering and Design, 2020, 369: 110828. doi: 10.1016/j.nucengdes.2020.110828

    [16]

    OLIVEIRA A V S, STEMMELEN D, LECLERC S, et al. Parametric effects on the flow redistribution in ballooned bundles evaluated by magnetic resonance velocimetry[J]. Experimental Thermal and Fluid Science, 2021, 125: 110383. doi: 10.1016/j.expthermflusci.2021.110383

    [17]

    BRUSCHEWSKI M, JOHN K, WÜSTENHAGEN C, et al. Commissioning of an MRI test facility for CFD-grade flow experiments in replicas of nuclear fuel assemblies and other reactor components[J]. Nuclear Engineering and Design, 2021, 375: 111080. doi: 10.1016/j.nucengdes.2021.111080

    [18]

    HOSOKAWA S, YAMAMOTO T, OKAJIMA J, et al. Measurements of turbulent flows in a 2 × 2 rod bundle[J]. Nuclear Engineering and Design, 2012, 249: 2–13. doi: 10.1016/j.nucengdes.2011.11.035

    [19] 熊国欣, 李立本. 核磁共振成像原理[M]. 北京: 科学出版社, 2007.
图(15)
计量
  • 文章访问数:  73
  • HTML全文浏览量:  23
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-12-27
  • 修回日期:  2024-03-24
  • 录用日期:  2024-04-02
  • 网络出版日期:  2024-06-12

目录

/

返回文章
返回
x 关闭 永久关闭