发明内容
有鉴于此,本发明所要解决的技术问题是提出了一种非能动核电厂钢制安全壳热移出过程热移出过程的比例分析方法,从机理上分析AP600,AP1000和CAP1400等非能动钢制安全壳的传热传质过程,并以比例的形式(π群)量化传热传质过程对安全壳系统热移出能力的贡献。
为实现上述目的,本发明提供了一种非能动核电厂钢制安全壳热移出过程的比例分析方法,包括如下步骤:
a)将钢制安全壳热阱分为蒸发区子热阱、干壁区子热阱和过冷区子热阱;
b)分别确定非能动安全壳冷却系统投入前和非能动安全壳冷却系统投入后的蒸发区子热阱、干壁区子热阱和过冷区子热阱的面积;
c)对蒸发区子热阱、干壁区子热阱和过冷区子热阱的传热传质过程进行分析,并分别计算蒸发区子热阱、干壁区子热阱和过冷区子热阱的传输质量比例群和传输能量比例群。
进一步地,步骤a)中根据钢制安全壳的外壁面覆盖水膜与否和水膜换热特性,将钢制安全壳热阱分为三部分:干壁区子热阱(该部分没有覆盖水膜)、过冷区子热阱(该部分覆盖水膜,水膜储热温升为主要换热方式)和蒸发区子热阱(该部分覆盖水膜,水膜蒸发换热为主要换热方式。
进一步地,还包括建立钢制安全壳热移出过程比例分析的模型;模型模拟钢制安全壳内壁面冷凝、对流、辐射传热传质过程,模拟钢制安全壳外壁面蒸发、对流、辐射传热传质过程,以及模拟钢制安全壳的储热/释热和温升/温降过程。
进一步地,步骤b)中确定非能动安全壳冷却系统投入前和非能动安全壳冷却系统投入后蒸发区子热阱、干壁区子热阱和过冷区子热阱的面积的计算过程包括:
步骤(401)获得钢制安全壳总面积Ashell;
步骤(402)在非能动安全壳冷却系统投入前,获得干壁区子热阱的面积Adry、蒸发区子热阱的面积Aevap和过冷区子热阱的面积Asubc;
步骤(403)假定非能动安全壳冷却系统投入后钢制安全壳外壁面水膜润湿比为Fwet’;
步骤(404)根据钢制安全壳总面积Ashell和步骤(403)中钢制安全壳外壁面水膜润湿比Fwet’,计算非能动安全壳冷却系统投入后干壁区子热阱的面积Adry=(1–Fwet’)Ashell;
步骤(405)根据步骤(404),计算非能动安全壳冷却系统投入后钢制安全壳的湿区面积Awet=Fwet’Ashell;
步骤(406)根据过冷水膜能量方程,确定非能动安全壳冷却系统投入后过冷区子热阱的面积Asubc;
步骤(407)根据步骤(406),计算非能动安全壳冷却系统投入后蒸发区子热阱的面积Aevap=Awet-Asubc;
步骤(408)根据非能动安全壳冷却系统投入后蒸发区子热阱的面积Aevap和换热系数,计算水膜的蒸发质量流mevap;
步骤(409)判断非能动安全壳冷却系统投入后蒸发质量流量mevap是否小于或等于非能动安全壳冷却系统的冷却水量mPCS;如果水膜的蒸发质量流量mevap小于或等于非能动安全壳冷却系统的冷却水量mPCS,即:mevap≤mPCS,则Aevap即为蒸发区子热阱的面积;如果水膜的蒸发质量流量mevap大于非能动安全壳冷却系统的冷却水量mPCS,即:mevap>mPCS,则需减小外壁面水膜润湿比Fwet’,重新进行计算,直至水膜的蒸发质量流量mevap等于非能动安全壳冷却系统的冷却水量mPCS;
步骤(410)根据步骤(409),确定非能动安全壳冷却系统的投入后钢制安全壳的外壁面水膜润湿比Fwet。
进一步地,步骤(402)中干壁区子热阱的面积Adry等于钢制安全壳总面积Ashell,蒸发区子热阱的面积和过冷区子热阱的面积为0m2。
进一步地,根据步骤(410)的外壁面水膜润湿比Fwet,确定非能动安全壳冷却系统投入后干壁区子热阱的面积Adry=(1–Fwet)Ashell;确定非能动安全壳冷却系统投入后钢制安全壳的湿区面积Awet=Fwet Ashell;根据过冷水膜能量方程,确定非能动安全壳冷却系统投入后过冷区子热阱的面积Asubc;确定非能动安全壳冷却系统投入后蒸发区子热阱的面积Aevap=Awet-Asubc。
进一步地,步骤c)中确定钢制安全壳热阱的传热传质过程包括:
将钢制安全壳热阱分为蒸发区子热阱、干壁区子热阱和过冷区子热阱;
对蒸发区子热阱、干壁区子热阱和过冷区子热阱进行一维网格离散;
采用一维非稳态导热的数值方法确定蒸发区子热阱、干壁区子热阱和过冷区子热阱的温度场;
步骤c)中对钢制安全壳热阱的质量比例群和能量比例群的计算过程包括:根据确定的蒸发区子热阱、干壁区子热阱和过冷区子热阱的温度场,分别确定蒸发区子热阱、干壁区子热阱和过冷区子热阱的传输质量比例群和传输能量比例群。
进一步地,采用一维非稳态导热的数值方法确定设计基准事故后蒸发区子热阱、干壁区子热阱和过冷区子热阱的温度场包括:
步骤(301)假设设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的壁面液膜的温度为Tf’;
步骤(302)获得设计基准事故后安全壳内的气体的温度T和安全壳内的气体的总压P;
步骤(303)根据设计基准事故后安全壳气体的温度T和蒸汽分压Pstm,结合饱和水蒸气表,确定蒸发区子热阱、干壁区子热阱或过冷区子热阱的辐射传热系数hr、对流传热系数hc和冷凝传热系数hm;
步骤(304)根据步骤(303)的结果,计算设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的辐射传热量Qr、对流传热量Qc和冷凝传热量Qm;
步骤(305)确定设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的壁面液膜传热系数hf;
步骤(306)根据设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的辐射传热系数hr、对流传热系数hc和冷凝传热系数hm以及蒸发区子热阱、干壁区子热阱或过冷区子热阱的壁面液膜传热系数hf,确定液膜传热的综合换热系数he=1/[1/(hc+hm+hr)+1/hf];
步骤(307)获得设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的壁面的初始温度
;
步骤(308)根据一维非稳态导热的数值方法,确定设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的壁面温度Tw;
步骤(309)根据傅里叶导热定律,确定设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的液膜表面温度Tf;
步骤(310)判断温度Tf和假设的温度Tf’是否相同;如果热阱液膜表面温度Tf和假设的温度Tf’相差很小,则温度Tf即为蒸发区子热阱、干壁区子热阱或过冷区子热阱的液膜表面温度;如果温度Tf和假设的温度Tf’相差较大,则令Tf’=Tf,重新进行迭代计算;
步骤(311)根据步骤(310),确定设计基准事故后蒸发区子热阱、干壁区子热阱或过冷区子热阱的液膜表面温度Tf。
在本发明中,具体所涉及到的钢制安全壳热移出过程的比例分析还可以包括破口源比例分析、安全壳内气相比例分析、安全壳内固体热阱比例分析、钢制安全壳热阱比例分析、导流板热阱比例分析和屏蔽厂房热阱比例分析六部分。钢制安全壳热阱比例分析包括钢制安全壳的干壁区子热阱比例分析、钢制安全壳的过冷区子热阱热阱比例分析、钢制安全壳的蒸发区子热阱比例分析;屏蔽厂房热阱包括烟囱子热阱比例分析、屋顶子热阱比例分析、筒身子热阱比例分析。其中,钢制安全壳的干壁区子热阱比例分析、钢制安全壳的蒸发区子热阱比例分析、导流板热阱比例分析、烟囱子热阱比例分析、屋顶子热阱比例分析、筒身子热阱比例分析还包括环腔空气流动和传热分析。
本发明还针对上述安全壳系统热移出过程,给出了非能动核电厂安全壳系统发生设计基准事故(DBA)事故后,模拟安全壳系统热移出过程比例分析的简化模型,该模型的特征如下:
安全壳内空间为一个控制容积
模拟蒸汽释放源
模拟气相空间的储热/释热过程
模拟蒸汽在内部热阱壁面的传热传质过程
模拟内部热阱的储热/释热和温升/温降过程
模拟钢制安全壳内壁面冷凝、对流、辐射传热传质过程
模拟钢制安全壳外壁面蒸发、对流、辐射传热传质过程
模拟钢制安全壳的储热/释热和温升/温降过程
模拟导流板热阱的冷凝、对流、辐射传热传质过程
模拟导流板热阱的储热/释热和温升/温降过程
模拟屏蔽厂房热阱的冷凝、对流、辐射传热传质过程
模拟屏蔽厂房热阱的储热/释热和温升/温降过程
模拟环腔流体的流动特性
进一步,对于上述非能动核电厂安全壳系统热移出过程的比例分析方法,建立度量传热传质过程的质量比例群和能量比例群表达式,并以比例的形式(π群)量化传热传质过程对安全壳系统热移出能力的贡献。包括:
(1)质量比例群
定义通过蒸汽从破口进入安全壳的方式传输质量的能力为1,即:
πm,brk=1
通过蒸汽储存在气相空间的方式传输质量的能力为:
其中ρ为安全壳内气相密度,ρg,brk为破口处蒸汽的密度。
通过蒸汽在热阱壁面冷凝的方式传输质量的能力为:
其中
为破口处蒸汽的质量流量,
为蒸汽在热阱壁面的冷凝质量流量。
(2)能量比例群
定义通过蒸汽从破口进入安全壳的方式传输能量的能力为1,即:
πe,brk=l
通过蒸汽储存在气相空间的方式传输能量的能力为:
其中ug为安全壳气相内能,hg,brk为破口处蒸汽焓值,uf,T0和hf,T0为初始温度T0对应的液态水内能和焓值。
通过蒸汽冷凝后传给热阱热量的方式传输能量的能力为:
其中hstm,j为气相空间蒸汽的焓值,hif,j为冷凝液膜的焓值。
通过蒸汽冷凝后形成液膜的方式传输能量的能力为:
通过对流和辐射换热方式传输能量的能力为:
其中hq,j为对流和辐射换热系数,Aj为热阱壁面面积,ΔT为热阱壁面和气相空间的温差。
在质量比例群和能量比例群的表达式中,下角标“j”表示某类热阱,后续用到的热阱字母缩写和热阱的对应关系详见表1。
表1 安全壳系统热移出过程的热阱字母缩写和热阱的对应关系
缩写 |
热阱名称 |
hs |
内部热阱 |
ish,evap |
钢制安全壳蒸发区子热阱内壁面 |
xsh,evap |
钢制安全壳蒸发区子热阱外壁面 |
ish,dry |
钢制安全壳干壁区子热阱内壁面 |
xsh,dry |
钢制安全壳干壁区子热阱外壁面 |
ish,subc |
钢制安全壳过冷区子热阱内壁面 |
xsh,subc |
钢制安全壳过冷区子热阱外壁面 |
bf |
导流板热阱 |
br |
屏蔽厂房筒身子热阱 |
rf |
屏蔽厂房屋顶子热阱 |
ch |
屏蔽厂房烟囱子热阱 |
由此可见,本发明提出了一种非能动核电厂安全壳系统热移出过程的比例分析方法,将安全壳热移出过程中可以评价的热阱分为安全壳内气相空间、安全壳内固体热阱、钢制安全壳热阱、导流板热阱和屏蔽厂房热阱等,可以评价的PCS现象包括:对流传热、辐射传热、蒸发传热传质、冷凝传热传质、热阱储热和PCS冷却水储热等。本发明还分别建立了安全壳系统热移出过程中的热阱比例分析,包括破口源比例分析、气相比例分析、内部固体热阱比例分析、钢制安全壳热阱比例分析、导流板热阱比例分析和屏蔽厂房比例分析,并建立了相应的安全壳系统热移出过程比例分析的模型,系统地分析了上述六类安全壳系统热移出过程中的传热传质过程,给出了度量传热传质过程对安全壳系统热移出能力的比例群,包括质量比例群和能量比例群两类;并以比例的形式(π群)量化传热传质过程对安全壳系统热移出能力的贡献。
进一步,本发明适用于AP600、AP1000和CAP1400等非能动核电厂安全壳系统热移出过程的比例分析,可以评价任一时刻下钢制安全壳传热传质过程对安全壳系统热移出能力的贡献。另外,本发明还适用于热移出过程和AP1000、CAP1400相似的系统、试验或台架,指导试验或台架参数的设计,以及评价试验或台架模拟原型核电厂的逼真度。
以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明,以充分地了解本发明的目的、特征和效果。
具体实施方式
如图1所示,给出了本发明的非能动核电厂安全壳系统发生假想的DBA事故后,PCS热移出过程的示意图。当非能动核电厂安全壳系统发生假想的大破口失水事故(LOCA)或主蒸汽管道断裂事故(MSLB)后,高温、高压的流体从破口101释放进入安全壳的气相空间102。气相空间102储存质量和能量,安全壳压力迅速升高。当安全壳压力达到设定值并延迟适当时间后PCS投入,水从PCCWST水箱(图中未示)流出至钢制安全壳103的外壁面,钢制安全壳103外壁面大部分面积将覆盖一层水膜。在安全壳的气相空间102中,蒸汽一方面在安全壳内固体热阱104的壁面冷凝形成液膜,并将热量传输至固体热阱104的内部;另一方面,蒸汽在钢制安全壳103的内壁面冷凝形成液膜,并将热量传输至钢制安全壳103的内部和钢制安全壳103的外壁面。钢制安全壳103的外壁面主要通过液膜蒸发的方式将热量传输至钢制安全壳103和导流板105所围成的环腔上升段106,环腔上升段的气相密度变小。在密度差的驱动下,来自环境的空气从屏蔽厂房107的空气入口108进入由屏蔽厂房107的筒身和导流板105围成的环腔下降段109,然后经导流板105折返180°进入由导流板105和钢制安全壳103围成环腔上升段106,最终从一个高位排气口110排至外界大气环境。当PCS热移出能力超过热源释热能力时,安全壳压力开始下降,直至PCS热移出能力和热源释热能力逐渐匹配后,安全壳压力和温度趋于稳定。
结合图1,本发明具体所涉及到的安全壳系统热移出过程的比例分析包括破口源比例分析、安全壳内气相比例分析、安全壳内固体热阱比例分析、钢制安全壳热阱比例分析、导流板热阱比例分析和屏蔽厂房热阱比例分析六部分。钢制安全壳热阱比例分析包括钢制安全壳干壁区子热阱比例分析、钢制安全壳过冷区子热阱比例分析、钢制安全壳蒸发区子热阱比例分析;屏蔽厂房热阱包括烟囱子热阱比例分析、屋顶子热阱比例分析、筒身子热阱比例分析。其中,钢制安全壳干壁区子热阱比例分析、钢制安全壳蒸发区子热阱比例分析、导流板热阱比例分析、烟囱子热阱比例分析、屋顶子热阱比例分析、筒身子热阱比例分析还包括环腔空气流动和传热分析。
本发明针对上述安全壳系统热移出过程,给出了一种非能动核电厂安全壳系统热移出过程的比例分析方法,分析了不同热阱(热源)的传热传质过程以及计算各个热阱的质量比例群和能量比例群。
1、蒸汽释放源(破口源)
蒸汽释放源(尤其是蒸汽注入)是安全壳压力升高的根本原因。在比例分析中,需要选择一个重要的传热传质途径作为基准。一般地,选取蒸汽注入的质量和能量作为比例分析的基准,即:
πm,brk=1和πe,brk=1
2、安全壳内气相空间
安全壳具有很大的气空间容积,DBA事故后,气相空间将存储或释放质量和能量,同时伴随安全壳压力的升高或降低。特别是在喷放阶段,主要通过气相空间的储热能力来缓解安全壳压力的升高。
2.1 分析过程
发生DBA事故前,安全壳内全部为空气,初始的温度、压力均已知。在整个DBA事故进程中,认为安全壳内气相(空气和蒸汽)处于准热力学平衡态,气相空间的蒸汽处于充分搅混状态和饱和状态。
图2给出了DBA事故后安全壳的温度T计算过程,具体说明如下:
步骤(201)获得安全壳气相空间热阱中的气体的初始温度T0和初始压力P0;
步骤(202)假设DBA事故后的安全壳气相空间热阱中的气体的温度为T’;
步骤(203)根据气体状态方程和步骤(201)的安全壳气相空间热阱中的气体的初始温度T0和初始压力P0,确定步骤(202)的DBA事故后的安全壳气相空间热阱中的空气的压力Pair;
步骤(204)获得DBA事故后的安全壳气相空间热阱的气体的总压P;
步骤(205)根据步骤(203)的空气的压力Pair和步骤(204)的气体的总压P,计算DBA事故后的安全壳气相空间热阱中的蒸汽的分压Pstm=P-Pair;
步骤(206)根据步骤(205)的蒸汽的分压Pstm,查饱和水蒸气表,确定DBA事故后的安全壳气相空间热阱中的气体的温度T;
步骤(207)判断温度T和假设的温度T’是否相同;如果温度T和温度T’相差很小,则温度T为安全壳气相空间热阱中的气体的温度;如果温度T和温度T’相差较大,则令T’=T,重新进行迭代计算;
步骤(208)根据步骤(207),确定DBA事故后的安全壳气相空间的气体的温度T;
步骤(209)在DBA事故后,根据步骤(208)确定的安全壳气相空间热阱的气体的温度T和气体的总压P,结合饱和水蒸汽表确定安全壳气相空间热阱的气相密度ρ和气相内能ug等参数。
2.2 质量、能量比例群
本发明可以量化的气相空间传输质量和能量的比例群包括:
(1)通过储存在气相空间的方式传输质量的能力,πm,τ
(2)通过储存在气相空间的方式传输能量的能力,πe,τ
3、安全壳内部热阱
在DBA事故下,高能流体通过破口释放进入安全壳气相空间,使气相空间升温、升压。由于气相空间和与之接触的热阱壁面存在蒸汽浓度差,在热阱壁面将出现冷凝过程;由于气相空间和与之接触的热阱壁面存在温度差,在热阱壁面将出现自然对流和辐射传热过程。
3.1 分析过程
本发明计算内部热阱壁面传热传质过程和热阱温升的过程如下:
(1)查找安全壳内所有的固体热阱,诸如楼板、设备、楼梯和金属平台等,并进行编号。
(2)将每一个已经编号的内部热阱简化为下述热阱中的一类:
钢板子热阱——由纯钢组成
混凝土子热阱——由纯混凝土组成
钢模块墙子热阱——单面或两面衬有钢板,中间为混凝土层。
(3)对简化后的内部热阱进行一维网格离散。
(4)采用一维非稳态导热的数值方法求解内部热阱的温度场,计算内部热阱壁面的蒸汽冷凝质量、对流和辐射换热量以及热阱的储热量。
图3给出了DBA事故后安全壳内部热阱的一个壁面的传热传质计算过程:
步骤(301)假设事故后某一时刻内部热阱壁面液膜的温度为Tf’;
步骤(302)获得事故后安全壳内该时刻的气相温度T;
步骤(303)根据安全壳气相温度T和蒸汽分压Pstm,结合饱和水蒸气表,计算辐射、对流和冷凝传热系数hr,hc和hm;
步骤(304)根据步骤(303)的结果,计算辐射、对流和冷凝的传热量Qr,Qc和Qm;
步骤(305)获得液膜传热系数hf;
步骤(306)根据计算辐射、对流、冷凝和液膜传热的综合换热系数he,计算公式如下:he=1/[1/(hc+hm+hr)+1/hf];
步骤(308)根据一维非稳态导热的数值方法(第三类边界条件)计算该时刻热阱壁面温度Tw;
步骤(309)根据傅里叶导热定律,计算出液膜表面温度Tf;
步骤(310)判断温度Tf和假设的温度Tf’是否相同;如果液膜表面温度Tf和假设的温度Tf’相差很小,则温度Tf即为该时刻液膜表面温度;如果温度Tf和假设的温度Tf’相差较大,则令Tf’=Tf,重新进行迭代计算;
步骤(311)根据步骤(310),确定DBA事故后内部热阱的液膜表面温度Tf。
3.2 质量、能量比例群
本发明可以量化的内部热阱传输质量和能量的比例群包括:
(1)通过蒸汽在热阱壁面冷凝的方式传输质量的能力,πm,hs
(2)通过蒸汽在热阱壁面冷凝传给热阱热量的方式传输能量的能力,πe,fg,hs
(3)通过蒸汽在热阱壁面对流和辐射换热的方式传输能量的能力,πe,q,hs
(4)通过蒸汽在热阱壁面冷凝后形成液膜的方式传输能量的能力,πe,f,hs
4、钢制安全壳热阱
4.1 钢制安全壳分区方法
根据钢制安全壳外壁面是否覆盖水膜和水膜换热特性,将钢制安全壳分为三部分:干壁区子热阱(该部分没有覆盖水膜)、过冷区子热阱(该部分覆盖水膜,水膜储热温升为主要换热方式)和蒸发区子热阱(该部分覆盖水膜,水膜蒸发换热为主要换热方式
如图4所示,给出了钢制安全壳的蒸发区子热阱、干壁区子热阱和过冷区子热阱面积的计算过程:
步骤(401)获得钢制安全壳总面积Ashell;
步骤(402)在PCS投入前,得到干壁区子热阱面积Adry、蒸发区子热阱面积Aevap和过冷区子热阱面积Asubc;具体来说,安全壳外壁面未被水膜覆盖,干区面积Adry等于钢制安全壳总面积Ashell,蒸发区子热阱和过冷区子热阱面积为0m2,即:
Adry=Ashell
Asubc=0m2
Aevap=0m2
步骤(403)假定PCS投入后,钢制安全壳外壁面水膜润湿比为Fwet;
步骤(404)根据钢制安全壳总面积Ashell和步骤(403)中钢制安全壳外壁面水膜润湿比Fwet,计算干区面积Adry=(1–Fwet)Ashell;
步骤(405)根据步骤(404),计算湿区面积Awet=Fwet Ashell;
步骤(406)根据过冷水膜能量方程,计算出过冷区子热阱面积Asubc;
步骤(407)根据步骤(406),计算出蒸发区子热阱面积Aevap=Awet-Asubc;
步骤(408)根据蒸发区子热阱面积Aevap和换热系数,计算水膜的蒸发质量流mevap;
步骤(409)判断蒸发质量流量mevap是否小于或等于PCS冷却水量mPCS;如果水膜的蒸发质量流量mevap小于或等于PCS冷却水量mPCS,即:mevap≤mPCS,则Aevap即为蒸发区子热阱面积;如果水膜的蒸发质量流量mevap大于PCS冷却水量mPCS,即:mevap>mPCS,则需减小外壁面水膜润湿比Fwet,重新进行计算,直至水膜的蒸发质量流量mevap等于PCS冷却水量mPCS;
步骤(410)根据步骤(409),获得PCS投入后钢制安全壳外壁面水膜润湿比为Fwet。
进一步地,根据(410)确定的PCS投入后钢制安全壳外壁面水膜润湿比为Fwet,获得钢制安全壳的蒸发区子热阱、干壁区子热阱和过冷区子热阱的面积。
4.2 传热传质过程
对钢制安全壳蒸发区子热阱、干壁区子热阱和过冷区子热阱的传热传质过程描述如下:
(1)蒸发区子热阱。
气相空间通过自然对流、辐射和冷凝传热方式和液膜表面进行热量传输。
热量以导热方式从液膜表面传至钢制安全壳内壁涂层和钢制安全壳。
进入钢制安全壳的热量一部分被钢制安全壳吸收,另一部分通过导热方式传输至钢制安全壳外壁涂层和蒸发液膜外表面。
热量通过强制对流、辐射和蒸发传热方式将热量传输至屏蔽厂房环腔的上升段。
(2)干壁区子热阱。
气相空间通过自然对流、辐射和冷凝传热方式和液膜表面进行热量传输。
热量以导热方式从液膜表面传至钢制安全壳内壁涂层和钢制安全壳。
进入钢制安全壳的热量一部分被钢制安全壳吸收,另一部分通过导热方式传输至钢制安全壳外壁涂层。
热量通过强制对流和辐射传热方式将热量从钢制安全壳外壁涂层面传输至屏蔽厂房环腔的上升段。
(3)过冷区子热阱。
气相空间通过自然对流、辐射和冷凝传热方式和液膜表面进行热量传输。
热量以导热方式从液膜表面传至钢制安全壳内壁涂层和钢制安全壳。
进入钢制安全壳的热量一部分被钢制安全壳吸收,另一部分通过导热方式传输至钢制安全壳外壁涂层。
传输至钢制安全壳外壁涂层的热量全部被水膜储存,水膜温度上升至可引起液膜大量蒸发的温度。
如图3所示,还可以用于钢制安全壳的蒸发区子热阱、干壁区子热阱和过冷区子热阱的传热传质计算过程。
下面以钢制安全壳的蒸发区子热阱的壁面为例,说明其传热传质计算过程。
步骤(301)假设事故后某一时刻蒸发区子热阱壁面液膜的温度为Tf’;
步骤(302)获得事故后安全壳内该时刻的气相温度T;
步骤(303)根据安全壳气相温度T和蒸汽分压Pstm,结合饱和水蒸气表,计算辐射、对流和冷凝传热系数hr,hc和hm;
步骤(304)根据步骤(303)的结果,计算辐射、对流和冷凝的传热量Qr,Qc和Qm;
步骤(305)获得液膜传热系数hf;
步骤(306)根据计算辐射、对流、冷凝和液膜传热的综合换热系数he,计算公式如下:he=1/[1/(hc+hm+hr)+1/hf];
步骤(308)根据一维非稳态导热的数值方法(第三类边界条件)计算该时刻蒸发区子热阱壁面温度Tw;
步骤(309)根据傅里叶导热定律,计算出蒸发区子热阱液膜表面温度Tf;
步骤(310)判断温度Tf和假设的温度Tf’是否相同;如果蒸发区子热阱液膜表面温度Tf和假设的温度Tf’相差很小,则温度Tf即为该时刻蒸发区子热阱的液膜表面温度;如果温度Tf和假设的温度Tf’相差较大,则令Tf’=Tf,重新进行迭代计算;
步骤(311)根据步骤(310),确定DBA事故后蒸发区子热阱的液膜表面温度Tf。
进一步,下面以钢制安全壳的蒸发区子热阱的壁面为例,说明钢制安全壳蒸发区子热阱的储热、释热过程,包括:
(1)对蒸发区子热阱进行一维网格离散。
(2)采用一维非稳态导热的数值方法求解蒸发区子热阱的温度场,计算蒸发区子热阱壁面的蒸汽冷凝质量、对流和辐射换热量以及热阱的储热量。
进一步,钢制安全壳的干壁区子热阱和过冷区子热阱的传热传质过程以及钢制安全壳储热、释热过程与上述蒸发区子热阱相似,在此不做赘述。
4.3 质量、能量比例群
本发明可以量化的钢制安全壳热阱传输质量和能量的比例群包括:
(1)蒸发区子热阱
通过在钢制安全壳内/外壁面冷凝/蒸发的方式传输质量的能力
通过在钢制安全壳内/外壁面冷凝/蒸发的方式传输能量的能力
通过在钢制安全壳内/外壁面对流和辐射换热的方式传输能量的能力
通过钢制安全壳内/外壁面液膜储能的方式传输能量的能力
(2)干壁区子热阱
通过在钢制安全壳内/外壁面冷凝/蒸发的方式传输质量的能力
πm,xsh,dry=0
通过在钢制安全壳内/外壁面冷凝/蒸发的方式传输能量的能力
πe,fg,xsh,dry=0
通过在钢制安全壳内/外壁面对流和辐射换热的方式传输能量的能力
通过钢制安全壳内/外壁面液膜储能的方式传输能量的能力
πe,f,xsh,dry=0
(3)过冷区子热阱
通过在钢制安全壳内/外壁面冷凝/蒸发的方式传输质量的能力
πm,xsh,subc=0
通过在钢制安全壳内/外壁面冷凝/蒸发的方式传输能量的能力
πe,fg,xsh,subc=0
通过在钢制安全壳内/外壁面对流和辐射换热的方式
πe,q,xsh,subc=0
通过钢制安全壳内/外壁面液膜储能的方式传输能量的能力
5、导流板热阱
5.1 传热传质过程
如图1所示,导流板热阱和钢制安全壳组成了环腔上升段,导流板和屏蔽厂房筒身子热阱组成了环腔下降段。对于导流板热阱而言,其传热传质过程为:
(1)钢制安全壳干壁区子热阱以辐射传热方式将热量传输至导流板;
(2)钢制安全壳蒸发区子热阱水膜以辐射传热方式将热量传输至导流板;
(3)导流板以辐射传热方式将热量传输至屏蔽厂房筒身子热阱;
(4)导流板以对流传热方式将热量传输至环腔上升段;
(5)导流板以对流传热方式将热量传输至环腔和下降段;
(6)环腔上升段蒸汽以冷凝方式将热量传输至导流板。
采用集总参数法,根据能量守恒方程,可计算任意时刻导流板热阱的温度和对流、蒸发、辐射换热量。
5.2 质量、能量比例群
本发明可以量化的导流板热阱传输质量和能量的比例群包括:
(1)通过上升段环腔蒸汽在导流板壁面冷凝的方式传输质量的能力,πm,bf
(2)通过在导流板壁面冷凝的方式传输能量的能力,πe,fg,bf
(3)通过壁面对流和辐射换热的方式传输能量的能力
上升段对流,
上升段辐射,
下降段对流,
下降段辐射,
6、屏蔽厂房热阱
6.1 传热传质过程
屏蔽厂房是环腔流道的重要组成部分,如图1所示,沿气体流动方向,屏蔽厂房热阱可划分为筒身子热阱、屋顶子热阱和烟囱子热阱三部分:
(1)筒身子热阱。位于空气入口下方,和导流板围成环腔下降段。其传热过程为:
环腔下降段气体通过对流方式将热量传输至筒身子热阱壁面
导流板通过辐射方式将热量传输至筒身子热阱壁面
筒身子热阱的储热温升
(2)屋顶子热阱。安全壳穹顶上方、非能动安全壳冷却系统储水箱下方的锥形屋顶。其传热过程为:
环腔气体通过对流方式将热量传输至屋顶液膜表面
安全壳穹顶干壁区子热阱通过辐射方式将热量传输至屋顶液膜表面
环腔气体通过冷凝方式将热量传输至屋顶液膜表面
热量以导热方式从液膜表面传至屋顶子热阱壁面
屋顶子热阱的储热温升
(3)烟囱子热阱。位于屋顶的中心位置,非能动安全壳冷却系统储水箱的组成部分,该部分热阱围成环腔排气段。其传热传质过程为:
环腔气体通过对流方式将热量传输至烟囱液膜表面
环腔排气段气体通过冷凝方式将热量传输至烟囱液膜表面
热量以导热方式从液膜表面传至烟囱子热阱壁面
烟囱子热阱的储热温升
屏蔽厂房热阱壁面的传热传质过程以及热阱储热过程和内部热阱相似,在此不做赘述。
6.2 质量、能量比例群
本发明可以量化的屏蔽厂房热阱传输质量和能量的比例群包括:
(1)筒身子热阱
通过环腔蒸汽在热阱壁面冷凝的方式传输质量的能力,πm,br
在屏蔽厂房壁面冷凝的方式传输能量的能力,πe,fg,br
通过屏蔽厂房壁面液膜储能的方式传输能量的能力,πe,f,br
通过屏蔽厂房壁面对流和辐射换热的方式传输能量的能力,πe,q,br
(2)屋顶子热阱
通过环腔蒸汽在热阱壁面冷凝的方式传输质量的能力,πm,rf
在屏蔽厂房壁面冷凝的方式传输能量的能力,πe,fg,rf
通过屏蔽厂房壁面液膜储能的方式传输能量的能力,πe,f,rf
通过屏蔽厂房壁面对流和辐射换热的方式传输能量的能力,πe,q,rf
(3)烟囱子热阱
通过环腔蒸汽在热阱壁面冷凝的方式传输质量的能力,πm,ch
在屏蔽厂房壁面冷凝的方式传输能量的能力,πe,fg,ch
通过屏蔽厂房壁面液膜储能的方式传输能量的能力,πe,f,ch
通过屏蔽厂房壁面对流和辐射换热的方式传输能量的能力,πe,q,ch
为了进一步说明安全壳系统热移出过程的比例分析方法,本部分以钢制安全壳蒸发区子热阱为例,给出了能量比例分析的示例。
表2给出了利用本发明非能动核电厂钢制安全壳热阱热移出过程的比例分析方法,获得LOCA事故后1600s的时刻,钢制安全壳热阱传输能量的比例值。
表2 钢制安全壳热阱传输能量的比例值(t=1600s)
图5给出了利用本发明非能动核电厂钢制安全壳热阱热移出过程的比例分析方法,获得LOCA事故后,图中实线代表钢制安全壳蒸发区子热阱外壁面水膜蒸发传输热量(πe,fg,xsh,evap)的比例曲线,图中虚线代表过冷水储热传输热量(πe,f,xsh,subc)的比例曲线。
可以看出,本发明提出了一种非能动核电厂安全壳系统热移出过程的比例分析方法,将安全壳热移出过程中可以评价的热阱分为安全壳内气相空间、安全壳内固体热阱、钢制安全壳热阱、导流板热阱和屏蔽厂房热阱等,可以评价的PCS现象包括:对流传热、辐射传热、蒸发传热传质、冷凝传热传质、热阱储热和PCS冷却水储热等。本发明还分别建立了安全壳系统热移出过程中的热阱比例分析,包括破口源比例分析、气相比例分析、内部固体热阱比例分析、钢制安全壳热阱比例分析、导流板热阱比例分析和屏蔽厂房比例分析,并建立了相应的安全壳系统热移出过程比例分析的模型,系统地分析了上述六类安全壳系统热移出过程中的传热传质过程,给出了度量传热传质过程对安全壳系统热移出能力的比例群,包括质量比例群和能量比例群两类;并以比例的形式(π群)量化传热传质过程对安全壳系统热移出能力的贡献。
进一步,本发明适用于AP600、AP1000和CAP1400等非能动核电厂安全壳系统热移出过程的比例分析,可以评价任一时刻下多种传热传质过程对安全壳系统热移出能力的贡献。另外,本发明还适用于热移出过程和AP1000、CAP1400相似的系统、试验或台架,指导试验或台架参数的设计,以及评价试验或台架模拟原型核电厂的逼真度。
以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。