开式布雷顿移动堆电源热电转换运行特性稳态分析.pdf

SCIENCE&TECHNOLOGY INFORMATION科技资讯 2023 NO.19 能 源 与 环 境科技资讯SCIENCE&TECHNOLOGY INFORMATION开式布雷顿移动堆电源热电转换运行特性稳态分析吕知讯(中国原子能科学研究院 北京 102400)摘要：该文针对开式布雷顿移动堆电源发电系统热力循环过程及其运行特性，开展系统稳态运行特性建模研究，建立系统运行工况稳态计算模型。在此基础上，研究开式布雷顿循环系统输出功率、效率以及系统内空气流量随转子转速、堆芯功率的变化规律，为系统运行控制研究提供基础。研究表明：在典型工况中，低转速下空气流量随转速增加而明显增加，高转速下空气流量随转速变化不明显；低转速下输出功率随转速增加而增加，高转速下输出功率随转速增加而减小。换热器出口温度对反应堆功率的敏感性随转速增加而降低；输出功率和输出功效率都随反应堆功率增加而增加。关键词：移动堆电源 开式布雷顿循环系统 压气机 涡轮机 稳态分析中图分类号：TL351.1文献标识码：A 文章编号：1672-3791(2023)19-0181-07Steady-State Analysis of Thermoelectric Conversion Operation Characteristics of Open Brayton Cycle Mobile Nuclear ReactorLYU Zhixun(China Institute Of Atomic Energy,Beijing,102400 China)Abstract:Aiming at the thermal cycle process and operation characteristics of the small reactor power generation system of the open Brayton system,the steady state operation characteristics modeling research of the system is carried out,and the steady state calculation model of the system operation conditions is established.On this basis,the output power and efficiency of the open Brayton cycle power generation system as well as the variation of the air flow in the system with the rotor speed and reactor power are studied.These studies can provide a basis for the system operation control research.The research shows that in typical working conditions,the air flow rate increases significantly with the increase of speed at low speed,while the air flow rate does not change significantly with speed at high speed;The output power increases with the increase of speed at low speed,and decreases with the increase of speed at high speed.The sensitivity of heat exchanger outlet temperature to reactor power decreases with the increase of rotating speed;The output power and efficiency increase with the increase of reactor power.Key Words:Mobile nuclear reactor power supply;Open Brayton system;Compressor;Turbine;Steady-state analysis气冷移动堆间接开式布雷顿循环是满足小型可移动核电源的一种优选方案，这种电源具备结构紧凑、启动快、重量轻等优点。系统由两个回路组成：一是高压气体闭式循环系统，由闭式风机驱动氦氙冷却剂，将反应堆热量通过换热器导出；二是开式布雷顿循环系统，以环境空气直接作为做功工质，通过布雷顿发动机发电。相比于闭式布雷顿循环，开式循环系统不需要配置热排放系统与二回路冷却剂储存系统，系统整体重量轻、体积小，可以更好地与移动载具进行匹配。王朝蓬1建立了单轴微型燃气轮机的性能计算模型，编制相应计算程序，对某型微燃机进行了仿真计算，并对该微燃机的气候、高度特性进行了计算分析。DOI：10.16661/ki.1672-3791.2303-5042-6575作者简介：吕知讯（1998），男，硕士在读，研究方向为反应堆安全。181SCIENCE&TECHNOLOGY INFORMATION科技资讯能 源 与 环 境 2023 NO.19 SCIENCE&TECHNOLOGY INFORMATION科技资讯钟亚飞2利用MATLAB软件，建立了包括C30微型燃气轮机部件和性能参数换算模块的整体模型，计算了微型燃气轮机在不同环境下的发电量和发电效率。都静3建立回热循环微燃机动态模型，并利用 Bowman TG80微燃机的实际运行数据开展数据验证工作验证了其动态模型的可信性。马同玲等人4建立闭式布雷顿循环发电系统系统计算模型，研究闭式布雷顿循环发电系统比功率和效率随涡轮入口总温及效率、压气机入口总温、压气机压比与效率、累积总压恢复系数等的变化规律，对系统的比功率和效率进行灵敏度分析。郭凯伦5对兆瓦级核电推进系统的动态布雷顿热电转换方式进行特性分析，比较兆瓦级核电推进系统气体透平循环在采用不同比例混合物作为工质时的循环效率，并对参数变化对循环效率的影响进行了研究。相比于闭式布雷顿循环移动堆电源，开式布雷顿移动堆电源的质量轻体积小，更能与车载进行匹配，但目前尚未对基于高温换热器热源的开式布雷顿循环移动堆电源开展研究，对于不同工况、参数变化时开式布雷顿循环系统特性缺少定量分析，其特性只能参考微型燃气轮机开式布雷顿循环和核反应堆闭式布雷顿循环研究成果。鉴于此，本文将以开式布雷顿循环原理为基础，综合压气机与涡轮机特性，建立基于换热器热源的开式布雷顿循环系统稳态计算模型，对开式布雷顿移动堆电源热电转换运行特性进行研究。1 方法1.1 建模方法及假设建模对象为开式布雷顿循环发电系统（二回路系统和反应堆一回路冷却系统），其系统简化如图1所示，开式布雷顿系统温熵图如图2所示。以带回热的开式布雷顿循环为基础进行建模计算，系统流程如下：压气机对空气进行加压，高压空气进入回热器冷管段进行一次加热，一次加热空气进入换热器进行二次加热，高温高压气体进入涡轮机内膨胀做功，带动涡轮机和发电机转子转动，涡轮机喷出的高温空气进入回热器热管段作为一次加热热源，经过压气机压气机涡轮机涡轮机G回热器回热器Tt0,pt0Tt1,pt1Tt1,pt15Tt2,pt2Tt2,pt25Tt3,pt3Tt3,pt35Tt4,pt4Tt4,pt45Tt5,pt5 p1 p2 p3 p4换热器换热器核反应堆核反应堆TinTout风机风机 图1 开式布雷顿循环发电系统原理图图2 开式布雷顿系统温熵图182SCIENCE&TECHNOLOGY INFORMATION科技资讯 2023 NO.19 能 源 与 环 境科技资讯SCIENCE&TECHNOLOGY INFORMATION冷却的空气从回热器出口排入环境中。以二回路计算结果作为设计输入，对一回路系统流量、冷却剂温度等参数进行计算6-8。开式布雷顿移动堆电源系统稳态计算模型流程图如图3所示。依据循环系统内各部件特性，建立数学物理模型，以空气质量守恒、能量守恒以及转子转速相等为条件进行计算求解，获得主要设备及关键部位在不同稳态工况下的状态参数，模拟全工况下的开式布雷顿循环系统实际工作状态。建立的系统稳态分析程序包括：压气机模型、涡轮机模型、换热器模型、回热器模型、二回路压降模型、一回路模型、空气物性模型和氦氙气体物性模型。为简化模型计算，做出以下基本假设条件：工质为干空气，在系统内做一维定常流动；忽略工质泄漏，系统内各截面空气质量流量相等；系统除空气出口与入口外，与外界无能量和质量交换；压气机、涡轮机转子机械能损失不随温度、压强、转速变化，机械能损失为2%。1.2 关键设备物理模型1.2.1 压气机模型压气机在变工况和设计工况下的特性由一些基本参数流量、压比、转速和效率来反映，对特性曲线进行拟合，可根据空气流量和转子转速插值得出压气机增压比和绝热效率9。c=f(qmnn0)（1）c=f(qmnn0)（2）根据压气机特性参数计算出当前压比和质量流量下的压气机出口截面空气压强p1和出口截面空气温度T1：p1=p0c（3）T2=T1(1+(-1)/c-1c)（4）根据出口空气温度计算出空气的比焓h2：h2=1.01(T2-273.15)（5）为模拟由摩擦造成的能量损失，引入参数压气机机械效率，计算出压气机比耗功lc：lc=()h2-h1mc（6）1.2.2 回热器模型为了减少计算量并准确描述回热器的换热过程，假定：不考虑回热器与外界的相互传热；气体工质做一维流动；不考虑流体与金属壁的热传导和热辐射。初始化回热器热管段进口温度T4，设定回热器效率x，计算回热器冷管段出口温度T2。T2=x(T4-T1)+T1（7）T5=T4-x(T4-T1)（8）开始输入边界条件压气机特性计算压气机出口空气参数计算回热器冷管段计算换热器计算涡轮机特性计算涡轮机出口空气参数计算回热器热管段计算回热器热管段与冷管段参数是否相符合对回热器热管段出口温度迭代计算否计算系统出口空气质量流量系统内空气质量流量是否守恒对空气质量流量进行迭代计算否是二回路运算结束一回路冷却剂参数计算一回路闭式风机耗功计算堆芯换热计算堆芯出口冷却剂温度与换热器入口冷却剂温度是否相等是对一回路冷却剂流量迭代计算一回路运算结束否是 图3 开式布雷顿移动堆电源系统稳态计算模型流程图183SCIENCE&TECHNOLOGY INFORMATION科技资讯能 源 与 环 境 2023 NO.19 SCIENCE&TECHNOLOGY INFORMATION科技资讯式（7）、式（8）中：T1为回热器冷管段进口截面空气温度；T2为回热器冷管段出口截面空气温度；T4为回热器热管段进口截面空气温度；T5为回热器热管段出口截面空气温度；x为回热器度，可根据回热器设计和所处环境不同而选取不同数值，在稳定运行条件下，回热度的取值一般在0.80.9之间。在此设计中选取为0.910-11。1.2.3 涡轮机模型涡轮机特性包括两个自变量空气质量流量qm和转速n，膨胀比t和绝热效率t。根据空气质量流量守恒和三机一体机转子转速相等可以得出涡轮机膨胀比和绝热效率，计算得出涡轮机出口截面空气温度和压强。T4=T3 1-t(1-t-k-1k)（9）pt4=pt3t（10）计算得涡轮机进出口比焓分别为h3、h4，引入机械效率mt，计算出涡轮机输出比功lt：lt=h3-h4mt（11）涡轮机输出功减去压气机耗功等于轴输出比功ln，输出功等于空气质量流量乘以输出比功：PG=qmln（12）1.2.4 管道压降模型二回路系统通过管道将各设备进行连接，为简化计算只考虑空气在管道内流动产生的摩擦压降，采用Darcy公式选取合理的摩擦系数，选定各段管道管长和水力直径，摩擦压降12计算为pf=fLiDeiiu22（13）式（13）中：f为管道摩擦阻力系数；Li为第i段管道长度；i为第i段管道入口空气密度；u为空气流速；Dei为第i段管道的当量直径。1.2.5 一回路模型一回路模型包括堆芯模型和冷却剂流动模型，堆芯模型包括堆芯压降计算和堆芯传热计算。本研究不需要计算出燃料芯块温度和芯块冷却剂通道内冷却剂温度分布，因此可将堆芯作为轴向功率密度和径向功率密度均匀分布的简单热源进行简化计算。Pw=qmcoolcpcool(Tout-Tin)（14）式（14）中：Pw为反应堆功率；cp,cool为一回路冷却剂定压比热容；Tout为堆芯出口冷却剂温度；Tin为堆芯入口冷却剂温度；堆芯压降模型计算包括入口突缩压降pin、出口突扩压降pout、堆内冷却剂通道压降pc。pcore=pin+pc+pout（15）堆芯冷却剂通道内压降计算只包括摩擦压降pf和加速压降pa：pc=pf+pa（16）当设定堆芯出口冷却剂温度为固定值时，流量需要随功率变化。通过设计工况确定一回路闭式风机性能系数，在稳态工况下，一回路闭式风机压头与一回路系统压降相等，闭式风机计算计算为Pfeng=qmcool(pcore+pex1+ppip1)（17）式（17）中：pex1为一回路换热器压降；ppaip1为一回路管道压降。1.2.6 物性模型假设二回路系统冷却工质为干空气，计算时需要计算空气的密度、比焓以及定压比热容：air=1.293pair0.101 325273.15Tair（18）hair=1.005(Tair-273.15)（19）cpair=0.995 6+0.000 093(Tair-273.15)（20）式（18）式（20）中：pair为空气压强，单位为MPa；Tair为空气绝对温度，单位为K。根据空气热物性散点图13对空气比热比进行曲面拟合，拟合出数学模型：=9.01710-5p2+2.65910-7T2-1.52210-4Tp+0.0693p-2.10410-4T+1.43（21）相比于纯氦气，分子量为15 g/mol的氦氙气体混合物的传热系数高7%，使用15 g/mol分子量的氦氙可以减小热交换器尺寸，因此选用15 g/mol的氦氙冷却剂作为一回路传热工质进行计算14-15。一回路系统为封闭恒压系统，压强可以根据实现功能变化和一回路压力边界制造水平选定不同参数，此设计设为定压3.3 MPa。一回路氦氙冷却剂压降计算物性参数主要包括密度He,Xe、动力粘度He,Xe，可由拟合公式16计算得到：HeXe=47.268 86-0.046 01THeXe+1.471 0810-5T2HeXe（22）HeXe=1.924 2210-5+4.623 6910-8THeXe（23）1.2.7 换热器模型换热器模型输入值为换热功率、二回路进出口温度，以及一回路入温度，通过平均对数温差法，进行迭代计算可以得到一回路换热器出口温度。平均对数温差tm可以通过热功率与换热器参数求得。tm=kA（24）tm=Tout-T3-Tin+T2lnTout-T3Tin-T2（25）式（24）、式（25）中：为换热器热流量，等于堆芯184SCIENCE&TECHNOLOGY INFORMATION科技资讯 2023 NO.19 能 源 与 环 境科技资讯SCIENCE&TECHNOLOGY INFORMATION功率；k为换热系数，A为换热面积。通过换热器热流量、二回路换热器进出口空气温度、堆芯出口空气温度，可以计算出堆芯入口冷却剂温度与一回路所需冷却剂流量qm,cool。2 结果与讨论2.1 转速对运行特性影响分析通过电力管理和分配子系统（Power Management And Distribution subsystem，PMAD）对转速进行调节，通过调节转速对输出功率进行控制17。根据系统稳态工况计算程序计算输出结果确定典型工况下的输出轴功率和输出功效率随三机一体机转子转速变化规律，以及通过维持换热器温度所需要的堆芯功率（如图4图6所示）。通过计算结果曲线得出如下结论，当转速增加时，输出功率变化会根据转速不同呈现不同的变化趋势，且存在最值轴功率转速临界点（以下统称为转速临界点）。计算结果得到，在同一转速下，空气质量流量随换 10 00015 00020 00025 00030 00035 00040 0000.00.51.01.52.0输出功率/MW转速/rmp二回路换热器出口空气温度 800K 900K 1 000K 1 100K 1 200K 1 300K 1 400K图4 不同换热器出口温度轴功率与转速关系 10 00015 00020 00025 00030 00035 00040 000012345反应堆功率/MW转速/rmp二回路换热器出口温度 800K 900K 1 000K 1 100K 1 200K 1 300K 1 400K图6 不同换热器出口温度堆芯功率与转速关系图 10 00015 00020 00025 00030 00035 00040 0000.00.20.40.6输出功效率转速/rmp二回路换热器出口温度 800K 900K 1 000K 1 100K 1 200K 1 300K 1 400K图5 不同换热器出口温度输出功效率与转速关系185SCIENCE&TECHNOLOGY INFORMATION科技资讯能 源 与 环 境 2023 NO.19 SCIENCE&TECHNOLOGY INFORMATION科技资讯热器出口温度变化很微小，同一转速下，流量的极差与平均值的比值不超过0.3%，因此假设空气质量流量只受转速影响，空气质量流量随转速变化规律情况具体如图7所示。当发电系统在受到流量微小扰动时，输出功率受到如下影响：转速低于转速临界点，流量瞬间有微小提升，导致输出功率提升，要想稳定输出功率则应降低转速；转速低于转速临界点，流量瞬间有微小降低，转速减小，轴输出功率减小；转速高于转速临界点，流量瞬间有微小提升，转速增加，轴输出功率降低；转速高于转速临界点，流量瞬间有微小降低，转速减小，轴输出功率增加。因此，在设置PMAD系统根据输出功受扰动变化对转速进行调节时需要预先判断实际转速与当前工况转速临界值的关系。2.2 反应堆功率对运行特性影响分析通过控制堆芯功率的方式对一回路冷却剂温度进行调节，可能导致由于不断改变堆芯温度而引起的反应堆燃料可靠性问题的增加，且调节堆芯功率需要控制元件的频繁和连续运动，从而引入更多的寿命问题与可靠性问题，并且通过控制反应性调节输出功率的控制方法是缓慢的，难以及时克服对外界微小扰动对输出功率的影响，因此改变反应堆功率的控制方法通常被考虑用于启动、停机、长期功率水平变化、事故以及系统的非正常响应，或用于燃料消耗引起的补偿。根据计算结果可以得出，转子转速分别为16 000、20 000、23 000、28 000、33 000、38 000 rmp时，堆芯功率每提升10 kW，输出换热器出口空气温度升高38.92 K、17.46 K、9.95 K、5.03 K、3.36 K、3.01 K。转速越大，换热器二回路侧温度对于反应堆功率变化敏感性越低。由于受压气机喘振特性影响，转速从 16 00038 000 rmp区间内的系统输出功率随反应堆功率变化范围如图8所示。输出功效率和换热器出口温度随堆芯功率 10 00015 00020 00025 00030 00035 00040 0000.51.01.52.02.53.03.5空气流量/(kg/s)转速/rmp图7 空气流量与转速关系图 012340.00.20.40.60.81.01.2输出功率/MW反应堆功率/MW图8 输出功率范围与反应堆功率关系186SCIENCE&TECHNOLOGY INFORMATION科技资讯 2023 NO.19 能 源 与 环 境科技资讯SCIENCE&TECHNOLOGY INFORMATION变化情况具体如图9、图10所示。3 结论本文通过研究，得到如下结论。（1）开式布雷顿循环发电系统的空气质量流量随转速增加而增大，典型工况下输出功率和输出效率随转速变化存在最大值，低转速下输出功率和输出功效率随转速增加而升高，高转速下输出功率和输出功效率随转速增加而降低。（2）相同转速下，输出功率、输出功效率以及二回路换热器出口空气温度随堆芯功率提升而增大。转子转速越大，二回路换热器出口空气温度对于反应堆功率变化敏感性越低，输出功对反应堆功率变化的敏感性越高。本文研究基于稳态工况下对开式布雷顿系统运行特性进行分析，因此在瞬态分析研究方面具有局限性，可以为基于高温换热器热源的开式布雷顿热电转换模块的瞬态系统分析与系统设计优化提供参考。PMAD系统对转速进行调节需要预先判断实际转速与当前工况转速临界值的关系，因此后续计算不同工况的准确转速临界值对系统控制研究具有重要意义。参考文献1 王朝蓬.微型燃气轮机稳态控制规律研究与特性计算J.科学技术与工程,2010(9):2255-2259.2 钟亚飞.基于MATLAB的微型燃气轮机建模与仿真研究J.发电设备,2019,33(5):325-330.3 都静.回热循环微燃机动态模拟与验证J.中国电机工程学报,2021,41(9):3207-3215.4 马同玲,张扬军,王正等.闭式布雷顿循环发电系统热力过程建模及其参数影响研究J.推进技术,2022,43(7):332-340.5 郭凯伦.兆瓦级核电推进系统布雷顿循环热电转换特性分析J.原子能科学技术,2019,53(1):16-23.6 袁春,陈彬兵.微型燃气轮机发电技术M.北京:机械工业出版社,2012:57-72.7 骆广琦,桑增产.航空燃气涡轮发动机数值仿真M.北京:国防工业出版社,2007:132-189.8 沈维道,童钧耕.工程热力学M.北京:高等教育出版社,2016:304-314.9 唐丽丽.燃气轮机系统的建模与仿真研究D.重庆:重庆大学,2008.10 杨帆.回热循环微型燃气轮机动态特性研究D.镇江:江苏大学,2022.11 都静.微型燃气轮机循环的动态模拟与性能分析D.北京:中国科学院,2020.12 孙中宁.反应堆热工水力学M.哈尔滨:哈尔滨工程大出版社,2017:74-84.13 张家荣,赵廷元.工程常用物质的热物理性质手册M.北京:新时代出版社,1987:203-204.14 MALIK A,ZHENG Q,QURESHI S R,et al.Effect of Helium Xenon as Workingfluid on Centrifugal Compressor of Power Conversion Unit of Closed Brayton Cycle Power PlantJ.International Journal of Hydrogen Energy,2020,46:7546-7557.15 MALIK A,ZHENG Q,QURESHI S R,et al.Effect of Helium Xenon as Workingfluid on the Compressor of Power Conversion Unit of Closed Brayton Cycle HTGR Power PlantJ.International Journal of Hydrogen Energy,2020,45:10119-10129.16 李杨柳.氦氙混合气体冷却反应堆单通道程序开发J.原子能科学技术,2017,51(1):41-45.17 WRIGHT S A,LIPINSKI R J,VERNON,M E,et al.Closed Brayton Cycle Power Conversion Systems for Nuclear Reactors:United StatesR/OL.(2006-04-01)2023-01-23.https:/www.osti.gov/servlets/purl/1177051.0.00.51.01.52.02.53.03.54.00.00.10.20.30.40.5输出功效率反应堆功率/MW转子转速 16 000rmp 20 000rmp 23 000rmp 28 000rmp 33 000rmp 38 000rmp图9 输出功效率与反应堆功率关系图10 换热器出口温度与反应堆功率关系187