高级检索

留言板

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

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

铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发

潘麒文 欧文澜 辜峙钘 龚政宇 戴嘉宁

潘麒文, 欧文澜, 辜峙钘, 龚政宇, 戴嘉宁. 铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发[J]. 原子核物理评论, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
引用本文: 潘麒文, 欧文澜, 辜峙钘, 龚政宇, 戴嘉宁. 铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发[J]. 原子核物理评论, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
Qiwen PAN, Wenlan OU, Zhixing GU, Zhengyu GONG, Jianing DAI. Development of a Heat Exchanger Module for a Transient Safety Analysis MPC_LBE Program for Lead-bismuth Reactors[J]. Nuclear Physics Review, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
Citation: Qiwen PAN, Wenlan OU, Zhixing GU, Zhengyu GONG, Jianing DAI. Development of a Heat Exchanger Module for a Transient Safety Analysis MPC_LBE Program for Lead-bismuth Reactors[J]. Nuclear Physics Review, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125

铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发

doi: 10.11804/NuclPhysRev.40.2022125
基金项目: 四川省自然科学基金资助项目(2022NSFSC0253)
详细信息
    作者简介:

    潘麒文(1998−),男,辽宁庄河人,硕士研究生,从事能源动力研究;E-mail: 39899382@qq.com

    通讯作者: 辜峙钘,E-mail: guzhixing17@163.com
  • 中图分类号: TL99

Development of a Heat Exchanger Module for a Transient Safety Analysis MPC_LBE Program for Lead-bismuth Reactors

Funds: Natural Science Foundation of Sichuan Province(2022NSFSC0253)
More Information
  • 摘要: 铅铋堆作为第四代先进核能系统之一,具有优良的中子经济性、固有安全性。为提高其紧凑性、安全性,铅铋堆主冷却剂系统倾向于采用一体化池式结构设计理念,但该设计理念同时也造成了复杂的热工水力问题。铅铋堆多物理耦合瞬态安全分析程序MPC_LBE在此背景下开发,但该程序的换热器模块采用了定温简化模型,无法模拟换热器一、二次侧之间的换热过程,事故瞬态模拟过于保守,偏离实际工况。为此,开展了铅铋堆换热器模块数值模拟方法研究,换热器一次侧、管壁以及二次侧均采用了一维通道等效处理手段,构建了其数值传热模型,编制了程序代码,并与MPC_LBE实现了外部显式耦合。对于该数值传热模型,单独进行了稳态验证及时间步长敏感性分析,结果显示,显式耦合策略的时间步长敏感性较大,而隐式耦合策略时间步长设置对模拟结果几乎无影响。对耦合了该数值传热模型的新MPC_LBE程序,开展了自然循环铅铋反应堆稳态模拟应用,同时添加了以一次侧出口段的设计温度为基准的二次侧冷却剂流量调节系统。
  • 图  1  套管式换热器(在线彩图)

    图  2  换热器网格图(在线彩图)

    图  3  显式耦合算法 (在线彩图)

    图  4  隐式耦合算法 (在线彩图)

    图  5  显式耦合方案下换热器温度的演化过程

    图  6  显式耦合方案下换热器稳态温度分布

    图  7  隐式耦合方案下换热器温度的演化过程

    图  8  隐式耦合方案下换热器稳态温度分布

    图  9  MPC_LBE主程序模块的主体框架(在线彩图)

    图  10  MPC_LBE主程序4号通道温度演化过程

    图  11  换热器一次侧温度演化过程

    图  12  MPC_LBE与换热器耦合稳态时一回路铅铋温度分布(在线彩图)

    图  13  MPC_LBE未与换热器耦合稳态时一回路铅铋温度分布(在线彩图)

    图  14  稳态时堆芯4号通道燃料棒温度分布(在线彩图)

    图  15  稳态时一回路冷却剂矢量图

    表  1  铅铋堆套管式换热器的设计参数

    名称 单位壳程 管程
    入口温度 K673 488.15
    出口温度 K573 503.15
    流量 kg·s−1158.844 40.21
    压降Pa863.33 17 503
    传热面积 m 25.86 /
    换热管数量142/
    换热管规格mm外管:∅30 × 2.5/
    mm内管:∅25 × 2.5/
    换热管有效换热长度mm1932/
    下载: 导出CSV

    表  2  二次侧水的物性参数

    参数数值
    ρ2 /(g·m−3)1 000
    Cp,2 /J·(kg·K)−14.2×103
    λ2 /W·(m·K)−1 0.645
    μ2 /(Pa·s)124.6×10−6
    下载: 导出CSV

    表  3  管壁材料的物性参数

    参数数值
    ρw /(g·m−3)7.93×103
    Cp,w /[J·(kg·K)−1]4.6×103
    λw /[W·(m·K)−1] 21.5
    下载: 导出CSV
  • [1] WU Yican. Engineering, 2016, 2(1): 262. doi:  10.1016/J.ENG.2016.01.023
    [2] 吴宜灿, 王明煌,黄群英,等. 核科学与工程, 2015, 35(2): 213. doi:  10.3969/j.issn.0258-0918.2015.02.004

    WU Yican, WANG Minghuang, HUANG Qunying, et al. Nuclear Science and Engineering, 2015, 35(2): 213. (in Chinese) doi:  10.3969/j.issn.0258-0918.2015.02.004
    [3] 吴宜灿, 柏云清, 宋勇, 等. 核科学与工程, 2014, 34(2): 201. doi:  10.3969/j.issn.0258-0918.2014.02.009

    WU Yican, BAI Yunqing, SONG Yong, et al. Nuclear Science and Engineering, 2014, 34(2): 201. (in Chinese) doi:  10.3969/j.issn.0258-0918.2014.02.009
    [4] GU Zhixing, ZHANG Qingxian, GU Yi, et al. Nuclear Science and Techniques, 2021, 32(5): 86. doi:  10.1007/s41365-021-00887-x
    [5] 张玲, 辜峙钘, 戴嘉宁, 等. 核技术, 2022, 45(10): 87. doi:  10.11889/j.0253-3219.2022.hjs.45.100602

    ZHANG Ling, GU Zhixing, DAI Jianing, et al. Nuclear Techniques, 2022, 45(10): 87. (in Chinese) doi:  10.11889/j.0253-3219.2022.hjs.45.100602
    [6] ZHANG L, SONG T, GU Z, et al. Kerntechnik, 2023, 88(2): 240.
    [7] 张巍, 李净松, 施慧烈. 核动力工程, 2022, 43(3): 38.

    ZHANG Wei, LI Jingsong, SHI Huilie. Nuclear Power Engineering, 2022, 43(3): 38. (in Chinese)
    [8] 丁雪友, 陈志强, 文青龙. 核动力工程, 2021, 42(4): 21.

    DING Xueyou, CHEN Zhiqiang, WEN Qinglong. Nuclear Power Engineering, 2021, 42(4): 21. (in Chinese)
    [9] 陈森. 加速器驱动铅铋冷却自然循环次临界堆无保护瞬态分析研究[D]. 合肥: 中国科学技术大学, 2014.

    CHEN Sen. Unprotected Transient Analysis of Accelerator Driven LBE-Cooled Natural Circulation Sub-critical Reactor [D]. Hefei: University of Science and Technology of China, 2014. (in Chinese)
    [10] 席文宣. 基于特殊介质的高效换热器流动传热特性研究[D]. 合肥: 中国科学院大学(中国科学院工程热物理研究所), 2018.

    XI Wenxua. Investigation on Flow and Heat Transfer Characteristics of High Efficiency Heat Exchangers Based on Special Medium [D]. University of Chinese Academy of Sciences(Institute of Engineering Thermophysics, Chinese Academy of Sciences), 2018. (in Chinese)
    [11] 丁雪友, 文青龙, 阮神辉. 铅铋快堆螺旋管直流蒸汽发生器热工水力模型研究[C]/第十六届全国反应堆热工流体学术会议暨中核核反应堆热工水力技术重点实验室 2019 年学术年会论文集. 2019: 196.

    DING Xueyou, WEN Qinglong, RUAN Shenhui. A thermal Hydraulic Model for a Helical Coiled Tube Once Through Steam Generator of Lead Cooled Fast Reactor[C]//Proceedings of the 16th National Conference on Reactor Thermal Fluids and the 2019 Annual Conference of the Key Laboratory of Thermal and Hydraulic Technology of Nuclear Nuclear Reactors. 2019: 196. (in Chinese)
    [12] 郭政, 阴继翔, 易文杰. 中国科技论文, 2021, 16(2): 217. doi:  10.3969/j.issn.2095-2783.2021.02.015

    GUO Zheng, YIN Jixiang, YI Wenjie. China Sciencepaper, 2021, 16(2): 217. (in Chinese) doi:  10.3969/j.issn.2095-2783.2021.02.015
    [13] 邢海坤. U型管蒸汽发生器仿真模型及动态特性研究[D]. 北京: 华北电力大学, 2014.

    XING Haikun. The Simulation Model and Research on the Dynamic Characteristics of U-tube Steam Generator[D]. Beijing: North China Electric Power University, 2014. (in Chinese)
    [14] 陈佳跃. 自然循环蒸汽发生器热工水力过程的数值仿真 [D]. 广州: 华南理工大学, 2013.

    CHEN Jiayue. Numerical Simulation of Natural Circulation Stam Generator’s Thermal-hydraulic Process[D]. Guangzhou: South China University of Technology, 2013. (in Chinese)
    [15] 张羽, 孙宝芝, 张国磊. 原子能科学技术, 2012, 46(1): 57.

    ZHANG Yu, SUN Baozhi, ZHANG Guolei. Atomic Energy Science and Technology, 2012, 46(1): 57. (in Chinese)
    [16] 刘乐, 陈文振, 黄文. 核科学与工程, 2020, 40(4): 596. doi:  10.3969/j.issn.0258-0918.2020.04.012

    LIU Le, CHEN Wenzhen, HUANG Wen. Nuclear Science and Engineering, 2020, 40(4): 596. (in Chinese) doi:  10.3969/j.issn.0258-0918.2020.04.012
    [17] 缠阿芳. 压水堆核电站U型管蒸汽发生器建模与仿真 [D]. 北京: 华北电力大学, 2012.

    CHAN Afang. Modeling and Simulation of U-Tube Steam Generator in PWR Nuclear Power Plant[D]. Beijing: North China Electric Power University, 2012. (in Chinese)
    [18] 程蔡阳. 山西化工, 1982(2): 53.

    CHENG Caiyang. Shanxi Chemical Industry, 1982(2): 53. (in Chinese)
    [19] 何啸峰, 阳小华, 何金桥. 南华大学学报(自然科学版), 2006, 20(2): 47.

    HE Xiaofeng, YANG Xiaohua, HE Jinqiao. Journal of the University of South China (Natural Science Edition), 2006, 20(2): 47. (in Chinese)
    [20] 杨瑞昌, 覃世伟, 刘若雷. 工程热物理学报, 2006(1): 130. doi:  10.3321/j.issn:0253-231X.2006.01.041

    YANG Ruichang, QIN Shiwei, LIU Ruolei. Journal of Engineering Thermophysics, 2006(1): 130. (in Chinese) doi:  10.3321/j.issn:0253-231X.2006.01.041
    [21] 苏子威, 周涛, 刘梦影, 等. 核技术, 2013, 36(9): 090205. doi:  10.11889/j.0253-3219.2013.hjs.36.090205

    SU Ziwei, ZHOU Tao, LIU Mengying, et al. Nuclear Techniques, 2013, 36(9): 090205. (in Chinese) doi:  10.11889/j.0253-3219.2013.hjs.36.090205
    [22] GU Zhixing, LI Fei, GE Liangquan, et al. Annals of Nuclear Energy, 2019, 133: 491. doi:  10.1016/j.anucene.2019.05.053
    [23] XIAO Yulong, GU Zhixing, ZHANG Qingxian, et al. Annals of Nuclaer Energy, 2020, 141: 107308. doi:  10.1016/j.anucene.2020.107308
  • 加载中
图(15) / 表 (3)
计量
  • 文章访问数:  275
  • HTML全文浏览量:  121
  • PDF下载量:  10
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-12-07
  • 修回日期:  2023-03-08
  • 网络出版日期:  2024-02-04
  • 刊出日期:  2023-12-20

铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发

doi: 10.11804/NuclPhysRev.40.2022125
    基金项目:  四川省自然科学基金资助项目(2022NSFSC0253)
    作者简介:

    潘麒文(1998−),男,辽宁庄河人,硕士研究生,从事能源动力研究;E-mail: 39899382@qq.com

    通讯作者: 辜峙钘,E-mail: guzhixing17@163.com
  • 中图分类号: TL99

摘要: 铅铋堆作为第四代先进核能系统之一,具有优良的中子经济性、固有安全性。为提高其紧凑性、安全性,铅铋堆主冷却剂系统倾向于采用一体化池式结构设计理念,但该设计理念同时也造成了复杂的热工水力问题。铅铋堆多物理耦合瞬态安全分析程序MPC_LBE在此背景下开发,但该程序的换热器模块采用了定温简化模型,无法模拟换热器一、二次侧之间的换热过程,事故瞬态模拟过于保守,偏离实际工况。为此,开展了铅铋堆换热器模块数值模拟方法研究,换热器一次侧、管壁以及二次侧均采用了一维通道等效处理手段,构建了其数值传热模型,编制了程序代码,并与MPC_LBE实现了外部显式耦合。对于该数值传热模型,单独进行了稳态验证及时间步长敏感性分析,结果显示,显式耦合策略的时间步长敏感性较大,而隐式耦合策略时间步长设置对模拟结果几乎无影响。对耦合了该数值传热模型的新MPC_LBE程序,开展了自然循环铅铋反应堆稳态模拟应用,同时添加了以一次侧出口段的设计温度为基准的二次侧冷却剂流量调节系统。

English Abstract

潘麒文, 欧文澜, 辜峙钘, 龚政宇, 戴嘉宁. 铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发[J]. 原子核物理评论, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
引用本文: 潘麒文, 欧文澜, 辜峙钘, 龚政宇, 戴嘉宁. 铅铋堆瞬态安全分析程序MPC_LBE换热器模块开发[J]. 原子核物理评论, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
Qiwen PAN, Wenlan OU, Zhixing GU, Zhengyu GONG, Jianing DAI. Development of a Heat Exchanger Module for a Transient Safety Analysis MPC_LBE Program for Lead-bismuth Reactors[J]. Nuclear Physics Review, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
Citation: Qiwen PAN, Wenlan OU, Zhixing GU, Zhengyu GONG, Jianing DAI. Development of a Heat Exchanger Module for a Transient Safety Analysis MPC_LBE Program for Lead-bismuth Reactors[J]. Nuclear Physics Review, 2023, 40(4): 660-667. doi: 10.11804/NuclPhysRev.40.2022125
    • 以铅基材料(铅或铅铋共晶合金)作为冷却剂的铅铋反应堆,因其独特的优良经济性和固有安全性[1],被纳入了第四代先进核能系统候选堆型,并作为其主力堆型之一[2-3]。为提高其紧凑型、安全性,铅铋堆主冷却剂系统倾向采用一体化池式结构设计理念,但这一设计理念同时也造成了复杂的热工水力问题,如上/下腔室热分层问题和气体夹带问题等。为应对池式堆中涉及的这些特殊热工水力学问题,成都理工大学课题组[4-6]开发了基于CFD方法的多物理耦合瞬态安全分析程序MPC _LBE。但由于MPC_LBE程序换热器模型过于粗糙,采用了定温简化模型,即无论换热器一次侧入口温度如何变化,一旦其进入换热器后则立即变为人为设定值,造成模拟结果过于保守,尤其是事故瞬态工况,其模拟结果与实际情况偏差较大。为解决此问题,本文开展了铅铋堆换热器模块开发,构建了其数值传热模型,并将其与MPC_LBE程序进行了外部显式耦合。

      本文建立的数值传热模型中的一次侧、管壁及二次侧均采用一维通道等效处理。首先对换热器仿真模块单独进行稳态验证,评价该程序的可行性和准确性。并基于显式和隐式两种耦合策略对换热器仿真模块进行时间步长敏感性分析。最后,对耦合了该换热器仿真模块的MPC_LBE程序进行了自然循环铅铋堆稳态模拟应用。

    • 铅铋堆换热器主要分为套管式和螺管式换热器[7-8]两大类,本文以中国科学院核能安全技术研究所设计的双层套管式换热器进行数值计算模型构建,其简图如图1所示,其相关设计参数可见表1[9-10]。在模型构建时,对换热器的管道进行了简化处理,即将142根管道等效为单管通道,同时,将其中的一次侧、管壁和二次侧以一维通道的形式表示,二次侧则采用开式通道。一次侧铅铋和二次侧高压过冷水互为逆向流动,且均与外界绝热,如图2所示。一回路的高温铅铋冷却剂在换热器套管的壳程(一次侧)内流动,将热量传递给套管管程(二次侧)内的二回路高压过冷水,即实现一次侧与二次侧间的传热耦合。其中,二次侧采用压力为4 MPa的加压过冷水,入口温度为488.15 K,出口温度为503.15 K。一次侧铅铋的入口设计温度为673 K,出口的设计温度为573 K。

      图  1  套管式换热器(在线彩图)

      表 1  铅铋堆套管式换热器的设计参数

      名称 单位壳程 管程
      入口温度 K673 488.15
      出口温度 K573 503.15
      流量 kg·s−1158.844 40.21
      压降Pa863.33 17 503
      传热面积 m 25.86 /
      换热管数量142/
      换热管规格mm外管:∅30 × 2.5/
      mm内管:∅25 × 2.5/
      换热管有效换热长度mm1932/

      图  2  换热器网格图(在线彩图)

    • 换热器套管流动换热的控制方程仅考虑能量守恒方程,一次侧、二次侧及管壁的控制方程分别如下所示:

      $$ \frac{{\partial {\rho _1}{C_{p,\,1}}{T_1}}}{{\partial t}} + \frac{{\partial {\rho _1}{C_{p,\,1}}{T_1}{v_1}}}{{\partial z}} = {\lambda _1}\frac{{{\partial ^2}{T_1}}}{{\partial {z^2}}} - {q_{V,\,1 \to w}}\text{,} $$ (1)
      $$ \frac{{\partial {\rho _2}{C_{p,\,2}}{T_2}}}{{\partial t}} + \frac{{\partial {\rho _2}{C_{p,\,2}}{T_2}{v_2}}}{{\partial z}} = {\lambda _2}\frac{{{\partial ^2}{T_2}}}{{\partial {z^2}}} + {q_{V,\,w \to 2}}\text{,} $$ (2)
      $$ \frac{{\partial {\rho _w}{C_{p,\,w}}{T_w}}}{{\partial t}} = {\lambda _w}\frac{{{\partial ^2}{T_w}}}{{\partial {z^2}}}{\text{ + }}{q_{V,\,1 \to w}} - {q_{V,\, w \to 2}}\text{,} $$ (3)

      式中:ρ为密度,单位为kg/m3Cp为比热容,J/(kg·K);T为温度;z为单位结点长度,单位为m;v为流体流动速度,单位为kg/m;λ为热导率,单位为J/(m·K);t为时间;q为对流换热项(下标1、w和2分别对应换热器一次侧铅铋、管壁和二次侧高压过冷水)。

    • 在套管式换热器中有两个对流换热模型,分别为一次侧与管壁之间的对流换热模型式(4)及管壁与二次侧之间的对流换热模型式(5)[11]

      $$ q_{V,\, 1 \rightarrow w}=h_{1} S_{1,\,w}\left(T_{1}-T_{w}\right) , $$ (4)
      $$ q_{V,\, w \rightarrow 2}=h_{1} S_{2,\,w}\left(T_{w}-T_{2}\right) 。 $$ (5)

      其中:铅铋与管壁之间的努塞尔数关系式表示如下[12-15]

      $$ N{u_1} = 4.5 + 0.018 \times P{e_1}^{0.8} \text{。} $$ (6)

      过冷水与管壁之间的努塞尔数关系式如下[16-20]

      $$ N{u_2} = 0.023 \times R{e_2}^{0.8}P{r_2}^{0.4} \text{,} $$ (7)

      式中:h为对流换热系数;Nu为努塞尔数;Re为雷诺数;Pe为贝克莱数;Pr为普朗特数。

    • 液态铅铋的密度、动力粘度及热导率等物性参数均来自公开发布的拟合公式[21],依次表示如下:

      $$ {\rho _1}{\text{ = }}11113.6 - 1.34{T_1}\text{,} $$ (8)
      $$ {\mu _1} = 4.94 \times {10^{ - 4}} \times {\rm{exp}}\left( {\frac{{757.1}}{{{T_1}}}} \right)\text{,} $$ (9)
      $$ {\lambda _1} = 4.21 + 1.2 \times {10^{ - 2}} \times {T_1}\text{,} $$ (10)
      $$ {C_{p,1}} = 143.919。 $$ (11)

      管壁材质为不锈钢,二次侧过冷水和管壁的相关物性参数均为常数,分别如表2表3所列[7-8]

      表 2  二次侧水的物性参数

      参数数值
      ρ2 /(g·m−3)1 000
      Cp,2 /J·(kg·K)−14.2×103
      λ2 /W·(m·K)−1 0.645
      μ2 /(Pa·s)124.6×10−6

      表 3  管壁材料的物性参数

      参数数值
      ρw /(g·m−3)7.93×103
      Cp,w /[J·(kg·K)−1]4.6×103
      λw /[W·(m·K)−1] 21.5
    • 从前述方程(1)、(2)及(3)中可看出,一次侧液态铅铋、管壁和二次侧水三者温度是相互耦合关联的,对此本文采用了显式耦合、隐式耦合两种方案。显式耦合和隐式耦合在算法上的不同主要表现在三个能量守恒方程中对流换热项离散化的形式不同,且在同一时间步长下显式耦合算法是一次侧、管壁和二次侧依次单独迭代计算求解,如图3所示。而隐式耦合则是将三个能量守恒方程同时进行迭代计算求解,因此隐式耦合算法比显式耦合算法更复杂,但计算结果更准确,算法步骤如图4所示。

      图  3  显式耦合算法 (在线彩图)

      图  4  隐式耦合算法 (在线彩图)

    • 在验证方法上采用铅铋堆换热器的稳态运行参数进行程序验证,换热器一、二次侧入口设置为边界条件,因此验证的关注点便是换热器的一、二次侧出口温度。为评价耦合方案的性能,采用了0.01,0.10 s两种时间步长开展时间步长敏感性分析。

    • 图5给出了在两种时间步长下换热器出口面一次侧、管壁、二次侧的温度随时间的变化情况,当时间达到50 s时,三者的温度基本达到稳定。稳态下换热器一次侧、管壁、二次侧的轴向温度分布如图6所示,均与设计值符合较好。同时,由图6可以看出显式耦合方案的时间步长敏感性较为显著,尤其是一次侧及管壁温度。

      图  5  显式耦合方案下换热器温度的演化过程

      图  6  显式耦合方案下换热器稳态温度分布

    • 图7给出了在两种时间步长下换热器出口面一次侧、管壁、二次侧的温度随时间的变化情况,当时间达到50 s时,三者的温度基本达到稳定。稳态下换热器一次侧、管壁、二次侧的轴向温度分布如图8所示,同样,均与设计值符合较好,且可以看出隐式耦合方案几乎不存在时间步长敏感性问题。

      图  7  隐式耦合方案下换热器温度的演化过程

      图  8  隐式耦合方案下换热器稳态温度分布

    • MPC_LBE是对铅铋池式反应堆进行整体热工水力模拟的多物理耦合程序[22],该程序耦合了中子动力学、燃料棒导热以及冷却剂流动传热三大物理过程,如图9所示,中子动力学模块采用点堆动力学理论,通过基于欧拉-麦克劳林近似的半解析方法求解[23]。燃料棒导热模块采用基于二维轴对称圆柱坐标系的非稳态、非线性热传导微分方程理论,使用了变物性全隐式迭代算法。冷却剂流动传热模块采用不可压缩单相流体动力学理论,使用了典型的SIMPLE求解算法。在前期工作中,完成了该程序的验证,具体包括中子动力学模块、CFD模块以及整个耦合系统的验证[4]

      图  9  MPC_LBE主程序模块的主体框架(在线彩图)

      将开发完成的铅铋堆换热器程序模块与MPC_LBE主程序进行外部显式耦合,两个程序均独立运行,分别在上述两程序内添加数据交换模块及控制模块,在外部开展数据的交换。添加控制模块的目的是保障两个程序物理运行时间保持一致,两个程序的时间步长均控制在0.01 s。就数据交换而言,一回路程序MPC_LBE向换热器程序一次侧模块传递铅铋的温度数据,换热器模块的程序对其开展运算,得到新的一次侧、管壁、二次侧的温度数据,运算完成后换热器将新的一次测温度再传递给MPC_LBE程序,再进行下一时间步的运算。

    • 由于MPC_LBE主程序在换热器一次侧部分一共设立了10条通道,因此换热器二次侧通道数量也应设置为10条,同时为了使自动控制技术得到应用,换热器的二次侧程序添加了冷却剂流速调节系统,调节系统在反应堆运行到130.0 s时开始运行,每条换热管道的流速调节系统均独立运行,其调节标准为换热器一次侧出口温度与设计温度573 K的差值。如果差值为正,调节系统就会增大当前换热管道二次侧的流速;如差值为负,则减小当前管道二次侧的流速。根据差值的绝对值大小将调节机制分为三个档位,粗调、细调和微调。当绝对值大于10 K时进行粗调,粗调是以当前管道冷却水流速的10 %对其进行调节,每10.0 s进行一次;当绝对值大于1 K时进行细调,细调则是以当前二次侧管道冷却水流速的0.3 %对其进行调节,每1.0 s进行一次;当绝对值大于0.01 K时进行微调,微调则是以当前二次侧管道冷却水流速的0.02 %对其进行调节,每0.1 s进行一次。

      MPC_LBE程序的时间步长与换热器模块的时间步长保持一致,均为0.01 s。图10给出了稳态建立过程中堆芯4号通道中平面位置燃料芯块内表面、燃料芯块外表面温度,以及4号通道出口平面冷却剂温度随时间的变化情况。如图所示,在耦合框架下各参数在600 s时已经近似趋于稳定了。

      图  10  MPC_LBE主程序4号通道温度演化过程

      图11给出了稳态建立过程中该铅铋堆换热器区域内一次侧液态铅铋温度随时间的变化情况,具体包括一次侧入口、出口以及中平面温度。如图所示,在耦合框架下各参数在600 s时已经近似趋于稳定了,出口温度与设计值符合得较好。

      图  11  换热器一次侧温度演化过程

      图10图11中可知在物理时间达到600 s时稳态基本建立,图12给出了稳态条件下铅铋堆主回路系统冷却剂液态铅铋温度的分布情况,图13为MPC_LBE原始程序主回路稳态温度分布情况,值得注意的是,铅铋堆主容器为池式结构,两个图展示的并非当前模拟的反应堆主容器的真实尺寸比例,而是其二维轴对称情况。图13中的靛蓝色区域是当前模拟的铅铋堆的冷池部分,发现在没有耦合二回路程序的情况下,可以很清楚地看出冷池区域内的温度基本无差异。与图13相比,图12中的换热器处的10条管道呈现出温度梯度,表明换热器模块与MPC_LBE主程序耦合后实现了从一次侧到二次侧的热量传递,在冷池的右下角的靛蓝色区域,由于此区域的铅铋流量很低,导致换热过程很慢。

      图  12  MPC_LBE与换热器耦合稳态时一回路铅铋温度分布(在线彩图)

      图  13  MPC_LBE未与换热器耦合稳态时一回路铅铋温度分布(在线彩图)

      图14给出了稳态时堆芯4号通道燃料棒温度的分布情况,其同样也是以二维轴对称形式展示,如图所示,右侧温度偏低的部分是燃料棒的包壳部分,包壳部分的温度沿着轴向逐渐上升,且燃料芯块的最高温度并没有出现在中心位置而是中心偏上的位置,这是由于在轴向方向燃料棒与冷却剂的温差不同所导致的,在燃料棒中心偏下区域比中心偏上区域的温差大,因此燃料棒中心偏下区域的传热能力更强,从而使燃料棒的最高温区域上移。

      图  14  稳态时堆芯4号通道燃料棒温度分布(在线彩图)

      图15给出了稳态工况下反应堆主回路系统内冷却剂液态铅铋速度矢量的分布情况,换热器下部由于泵的存在冷却剂速度加快,在堆芯中通道截面减小,冷却剂速度加快。一次侧冷却剂在冷池与热池之间循环,将堆芯的热量通过换热器传递给了二次侧高压过冷水,实现了堆芯的冷却。

      图  15  稳态时一回路冷却剂矢量图

    • 鉴于铅铋堆多物理耦合瞬态安全分析程序MPC_LBE换热器模块采用了定温简化模型,其物理模型过于粗糙,模拟结果过于保守,无法模拟换热器内换热过程,为此本文开展了换热器模块开发,并与MPC_LBE实现了外部显式耦合。

      通过稳态模拟的结果可得出,使用显式耦合方案时间步长取0.01 s时,一次侧及二次侧出口温度的相对误差分别为0.467 57%, 0.151 96%;取0.10 s时,一次侧、二次侧出口温度的相对误差分别为0.467 59%, 0.151 96%。使用隐式耦合方案时,时间步长取0.01 s时,一次侧及二次侧出口温度的相对误差均分别为0.467 57%, 0.151 96%;时间步长取0.10 s时,一次侧及二次侧出口温度的相对误差均分别为0.467 57%, 0.151 96%,两者基本无差异。因此显式耦合方案的时间步长敏感性偏大,隐式耦合方案的时间步长敏感性几乎不存在,且隐式耦合算法更加精确。在换热器与MPC_LBE主程序进行外部显式耦合方面,一次侧出口温度的相对误差为0.058 28%,与设计值符合得较好。

参考文献 (23)

目录

    /

    返回文章
    返回