由于设计尺寸不同或加工精度的偏差，组件内部冷却剂通道缝隙宽度可能是不一致的；或者由于燃料辐照肿胀起泡、燃料表面结垢、堆内材料碎片进入等原因导致某个流道堵塞。使得燃料板两侧热工水力条件不同，两侧冷却剂冷却燃料板能力不同，这样燃料最高温度就不在燃料几何中心线上，此时即为不对称冷却。而现在的仿真程序没有考虑燃料元件两侧的不对称冷却情况，这在实际的仿真中不能完全模拟出实际的堆芯情况。不对称冷却问题的难点在于燃料最高温度的位置难以确定，不能通过直接积分得到燃料芯块平均温度与燃料内部各界面温度的关系式。本文利用Hermite插值积分，给出边界条件和积分平均值之间的关系，提出了一种改进的集总参数模型来求解板状燃料元件不对称冷却问题。

对[a, b]区间上的连续函数f(x)，Hermite两点插值积分表达式为：

H_0，0近似：

H_1，1近似：

燃料芯块内裂变释放出来的热量经导热传给燃料包壳，再由对流换热通过包壳传给冷却剂。图1是板状燃料元件导热示意图。燃料芯块厚度为δ_u，燃料包壳厚度为δ_c。由于燃料元件一般为长而薄的平板，轴向导热热阻要远大于横向导热热阻，因此可以忽略其轴向导热，仅考虑横向导热，因此可以用一维瞬态热传导方程求解燃料温度。

图  1  板状燃料元件温度分布示意图

Figure 1.  Distribution of Temperature in Plate Type Fuel Pellet

对于燃料芯块有：

对于燃料包壳有：

式中：ρ为密度，kg/m³；c_p为比热容，J/kg·K；λ为热导率，W/(m·K)；为燃料芯块功率密度，W/m³。

板状燃料芯块和包壳之间不存在气隙，在这里假设芯块和包壳压实紧密接触，两者之间没有接触热阻。则包壳内表面和燃料芯块之间为连续性边界条件：

左侧：

右侧：

式中：T_cil为燃料包壳左侧内表面温度，K；T_cir为燃料包壳右侧内表面温度，K。

包壳外表面和冷却剂之间为第三类边界条件：

左侧：

右侧：

式中：T_col为燃料包壳左侧外表面温度，K；T_cor为燃料包壳右侧外表面温度，K；T_fl为燃料包壳左侧冷却剂温度，K；T_fr为燃料包壳右侧冷却剂温度，K。

利用Hermite两点插值积分表达式即可给出温度平均值、边界热流密度和温度之间的关系。

计算时做如下假设：

1)忽略燃料轴向导热，仅考虑横向导热。

2)燃料芯块的热物性取决于芯块的平均温度，燃料包壳的热物性取决于包壳的平均温度。

3)芯块温度分布用二次曲线近似，包壳内温度分布用直线近似。

4)包壳和芯块紧密结合在一起，中间没有气隙。

采用H_1，1近似计算燃料芯块的积分平均温度：

采用H_0，0近似计算燃料包壳积分平均温度和积分平均热流密度：

对导热方程燃料芯块和燃料包壳在空间坐标进行积分，并应用边界条件(5)～(7)同时代入式(2) 、式(10)～式(13)，可以得到燃料芯块、左侧包壳和右侧包壳平均温度的微分方程：

联立方程可得：

式中：。

求解方程(17)～式(21)即可求得燃料内部各界面温度。若假设燃料芯块内温度分布为二次曲线形式T_u＝ax²＋bx＋c，则可根据δ_c和δ_c＋δ_u处温度和热流密度值，求出燃料最高温度和所在位置。

为了研究中国先进研究堆燃料组件仿真，首先需要建立起反应堆板状燃料组件的仿真模型，根据仿真结果的需求选取了美国GSE公司开发的核动力装置热工水力仿真软件THEATRe(Thermal Hydraulic Engineering Analysis Tools in Real Time)程序是由美国GSE公司开发的反应堆热工水力实时仿真分析工具。将所编写的FORTRAN算法程序移植到THEATRe仿真软件中，准备填写CARR堆的输入卡。

中国先进研究堆(CARR)是我国自主设计和建造的一座多用途、高性能研究反应堆。其堆芯由17个几何结构相同的标准燃料组件(编号1～17)和4个几何结构相同的控制棒-跟随体组件(编号18～21)构成。其中标准燃料组件平行排列了21块燃料板，从燃料组件左右外侧向内侧开始冷却剂流道间隙不同，从左右外侧开始间隙分别为1.56 mm、2.59 mm、2.45 mm、2.32 mm，其余冷却剂通道间隙都为2 mm。跟随体燃料组件平行排列了17块燃料板，除了最外侧冷却剂通道间隙为3.94 mm，其余间隙通道都为2.25 mm。根据CARR的结构与功率份额参数结合以上数据填写THEATRe程序的输入卡。

利用所编制程序计算CARR堆稳态工况下，各燃料组件间流量分配、燃料组件内流量分配问题和组件内部温度场分布。计算时每个燃料组件划分为一个通道，每个冷却剂通道轴向上均匀划分为17个控制体。

在计算标准燃料组件选择组件9，跟随体燃料组件选择组件18，在计算时假设组件内燃料径向功率均匀分布。

图2和图3分别给出了CARR中标准燃料组件和跟随体燃料组件内部流量分配计算结果，计算中标准燃料组件选择组件9，跟随体燃料组件选择组件18，在计算时假设组件内燃料径向功率均匀分布。从计算结果可以看出，流量分配主要受流道的尺寸影响，间隙宽的流道，冷却剂流量分配份额较大。从图中还可以看出本文计算结果和现有参考结果符合较好。

图  2  CARR堆芯燃料布置图

Figure 2.  Core Fuel Position Arrangement of CARR

图  3  CARR标准燃料组件流量分配

Figure 3.  The Standard Fuei Assembly Flow Distribution of CARR

从以上分析和计算结果可以看出，由于CARR燃料组件中边界流道尺寸不一致，各流道内冷却剂流量是不一样的，使得燃料元件两侧冷却条件不一样，这时燃料元件为非对称冷却。图4和图5分别显示了标准燃料组件9内冷却剂温度和燃料芯体最高温度沿流动方向分布。

图  4  CARR跟随体燃料组件流量分配

Figure 4.  The Following Fuel Assembly Flow Distribution of CARR

图  5  CARR标准燃料组件冷却剂温度

Figure 5.  Coolant Temperature Standard Fuel Components of CARR

从图中可见，组件最边缘通道1冷却剂温度最低，但与其它流道温度相差不大，这是由于边缘通道只有单侧燃料板加热，但同时流道间隙最窄，流量分配份额最小。由外而内流道冷却剂温度逐渐升高，这是由于冷却剂通道逐渐变窄，通道冷却剂流量减小，在燃料功率相同的情况下，温度升高较多。燃料组件内部的通道5和通道6的温度分布基本相同。燃料板5和燃料板6的燃料最高温度也分布基本相同，这说明组件边缘流道由于流量分配不均而导致为非对称换热，内部间隙相同的流道流量分配均匀为对称换热。表1和表2为本文计算结果和参考结果比较，从比较可以看出，冷却剂温度计算结果与参考结果符合较好，但燃料温度整体偏高，这可能是获得的燃料芯体物性不够准确造成的。

在板状燃料堆芯中，冷却剂通道之间彼此相隔，通道之间没有流动交混。因此当某个燃料组件内流道传热恶化时，可能由于热量不能及时导出而造成恶劣通道热工水力条件更加恶劣，直至燃料烧毁。板状燃料堆芯中冷却剂通道大都比较狭窄，在运行中由于燃料辐照起泡，燃料板结垢或异物堵塞通道等工况，使得冷却剂通道局部堵塞而造成传热恶化，因此有必要对这一现象进行分析。

事故描述：如图6所示，假设在额定工况下，10s时CARR堆芯标准组件11、12的中间流道(编号12)发生流道堵塞事故，在1s内流道流通面积被堵塞90%。假设在堵流事故中不能停堆。由于点堆模型不能反映局部热工水力参数变化对局部功率的影响，因此在瞬态计算时采用保守估计，认为组件功率在堵流后保持不变，同时各燃料板间功率均匀分布。

图  6  CARR标准燃料组件燃料芯块最高温度

Figure 6.  The Highest Temperature Standard Fuel of CARR

图7为堵流事故后，被堵塞流道所在的燃料组件11冷却剂流量变化曲线。从图中可以看出，组件内一个小流道堵塞对整个组件流量影响很小，通过组件的冷却剂流量仅由11.88kg/s下降至11.86kg/s。从下图7中可以看出被堵的通道12流量迅速降低，其它通道的流量均有所增加，与被堵通道相邻的通道11和通道13流量增加的更多一些。这是由于通道12堵塞后两侧燃料板的热量传给另一侧冷却剂通道11和13，使得这两个通道受热功率增加，而且从以上计算结果看出，对于CARR功率份额高的通道，其流量份额也高，因此与被堵流道相邻两个流道，流量增加的更多些。

图  7  单一通道堵塞示意图

Figure 7.  Single Channel Check Diagrammatic Sketch

图8为堵流事故中组件内各通道出口冷却剂温度计算结果。堵流后通道12流量减小，冷却剂温度升高。相邻通道11和13尽管流量有所增加，但是堵流通道两侧燃料板对通道11和13加热功率也增加，因此这两个通道冷却剂出口温度增加，其余通道冷却剂出口温度由于冷却剂流量增加而都有所降低。图9为被堵通道右侧燃料板12温度场计算结果。堵流前燃料板12两侧为对称传热，其温度场为对称分布有：T_cil＝T_cir，T_col＝T_cor，T_fl＝T_fr，且T_ci>T_co>T_f。堵流后，燃料板12左侧对流传热减弱，两侧为不对称传热，温度场为不对称分布，其中T_cil>T_col>T_cir>T_cor>T_fl>T_fr。其中：T_fl——燃料左侧通道冷却剂温度，T_col——燃料板左侧包壳外表面温度，T_cil——燃料板左侧包壳内表面温度，T_cir——燃料板右侧包壳内表面温度，T_cor——燃料板右侧包壳外表面温度，T_fr——燃料右侧通道冷却剂温度。

图  8  组件内流量分配

Figure 8.  Flow Distribution Within the Fuel Subassembly

图  9  通道出口、燃料板12和堵流后燃料板12温度变化曲线

Figure 9.  Temperature Variation Curve of the Channel Outlet, the Fuel Slab 12 Before and After Blockage

本文通过Hermite插值积分理论列出求解方程，利用Fortran语言将方程转换为算法程序，并将算法程序移植加入反应堆热工水力实时仿真程序THEATRe中。通过计算中国先进研究堆(CARR)的标准燃料组件和跟随体燃料组件的稳态温度分布，并与参考数值对比验证了此种方法的可行性。应用该方法模拟板型燃料组件流道堵塞事故，可以更准确的分析出不对称冷却条件下燃料元件的温度分布。