CN107784156A - 核动力装置蒸汽排放***参数的计算方法 - Google Patents

核动力装置蒸汽排放***参数的计算方法 Download PDF

Info

Publication number
CN107784156A
CN107784156A CN201710800214.0A CN201710800214A CN107784156A CN 107784156 A CN107784156 A CN 107784156A CN 201710800214 A CN201710800214 A CN 201710800214A CN 107784156 A CN107784156 A CN 107784156A
Authority
CN
China
Prior art keywords
mrow
msup
msub
steam
msubsup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710800214.0A
Other languages
English (en)
Other versions
CN107784156B (zh
Inventor
刘现星
罗小雨
孙海军
程堃
张靖
刘磊
王珺
唐征明
李广围
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
719th Research Institute of CSIC
Original Assignee
719th Research Institute of CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 719th Research Institute of CSIC filed Critical 719th Research Institute of CSIC
Priority to CN201710800214.0A priority Critical patent/CN107784156B/zh
Publication of CN107784156A publication Critical patent/CN107784156A/zh
Application granted granted Critical
Publication of CN107784156B publication Critical patent/CN107784156B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明涉及核动力装置蒸汽排放***参数的计算方法,包括如下步骤:初始化;计算蒸汽流量、减压装置系数、排放阀开度和管路压降;计算当前时刻的蒸汽流量、给水流量、循环倍率、平均换热系数、反应堆功率、主冷却剂平均温度、再循环水流量和循环水流量;计算得到当前时刻的蒸汽温度;获取当前时刻的饱和蒸汽压力;获取当前时刻n的传热区出口水焓值;计算当前时刻的汽水分离区温度;计算当前时刻的下降区温度;获取当前时刻的传热区入口焓值;计算当前时刻的蒸汽发生器一次侧换热系数;计算当前时刻的反应堆出口冷却剂温度;重复执行步骤S2至步骤S17;确定反应堆出口冷却剂温度和蒸汽压力的上限值。本发明简单易用。

Description

核动力装置蒸汽排放***参数的计算方法
技术领域
本发明涉及蒸汽排放,具体而言是核动力装置蒸汽排放***参数的计算方法。
背景技术
蒸汽排放***是保障核动力装置安全运行的一个重要***。在汽轮机负荷突然大幅度下降后,反应堆功率不能像汽轮发电机负荷那样快速地改变,蒸汽排放***的作用是在此过程中有控制地将蒸汽直接排放至冷凝器,从而通过为反应堆提供一个“人为”负荷的办法来降低核蒸汽供应***反应堆主冷却剂的温度瞬态和压力瞬态幅度,同时限定蒸汽***最高工作压力,确保核电厂的安全运行。
蒸汽排放***的参数设计要充分考虑***的动态效应。蒸汽排放***运行的动态过程涉及反应堆、主冷却剂***、蒸汽发生器特性,以及相应的控制策略与控制逻辑,由于核动力装置运行的耦合效应,参数的设计难以通过静态的分析计算完成。通常核动力装置蒸汽排放***的参数计算方法为:先依据经验初步确定排放容量和压力,然后建立核动力装置仿真***,通过***仿真研究对初始设计的排放容量和压力进行优化,确定具体的排放参数整定值,最终通过实际核动力装置的调试试验进行验证。而核动力装置仿真***涵盖反应堆物理、主冷却剂***、二回路***、电力***以及控制***,开发周期长、难度大、成本高,影响蒸汽排放***与相关***的迭代设计,从而影响核动力装置的总体设计进度。
因此,提供一种简单易用的可计算蒸汽排放***不同设计排放负荷下蒸汽压力和一次侧冷却剂温度上限值的核动力装置蒸汽排放***参数的计算方法十分必要。
发明内容
本发明的目的是提供一种简单易用的可计算蒸汽排放***不同设计排放负荷下蒸汽压力和一次侧冷却剂温度上限值的核动力装置蒸汽排放***参数的计算方法。
为实现上述目的,本发明采用如下技术方案:一种核动力装置蒸汽排放***参数的计算方法,包括如下步骤:
S1.初始化:
S1.1.令n=1,得到0时刻的反应堆功率P0,蒸汽发生器平均换热系数K0、循环倍率N0,蒸汽发生器二次侧传热区的出口水流量Gof 0、出口水焓值h2f 0,传热区出口蒸汽流量Gg 0、出口蒸汽焓值h2g 0,传热区入口流量G2i 0、入口焓值h2i 0,汽水分离区的温度T3 0、再循环水焓值h3f 0,给水流量Gf 0,下降区的温度T1 0
S1.2.获取0时刻的主冷却剂平均温度
通过以下公式计算得到0时刻的主冷却剂平均温度
式中:
——反应堆初始功率P0对应的主冷却剂***平均温度,℃;
△Tc——反应堆冷却剂控制死区,℃;
△Tm——反应堆冷却剂温度测量误差,℃;
△Tc由控制***设计给出,△Tm由温度测量仪表选型厂家给出;
S1.3.获取0时刻的蒸汽温度Ts 0
通过以下公式得到0时刻的蒸汽压力ps 0
ps 0=ps0+△ps(△Tc+△Tm)
式中:
ps0——主冷却剂设计平均温度下的蒸汽发生器出口蒸汽压力,Pa;
△ps——反应堆功率为P0时主冷却剂平均温度较设计值升高1℃时对应的蒸汽压力增加值,Pa;
查水蒸气物性参数表,得到与ps 0对应的饱和蒸汽温度,即0时刻的蒸汽温度Ts 0
S2.计算蒸汽流量、减压装置系数、上一时刻排放阀开度和蒸汽出口至排放阀前管路压降:
根据蒸汽排放***减温减压装置设计容量和压力,以及主冷却剂平均温度偏差,通过以下公式计算得到汽轮机目标负荷值对应蒸汽流量减压装置系数k1、上一时刻n-1的排放阀开度k2 (n-1)和蒸汽发生器蒸汽出口至排放阀前管路压降
式中:
——汽轮机目标负荷值对应蒸汽流量,kg/s;
P——目标负荷值,W;
G100——100%负荷对应蒸汽流量,kg/s;
式中:
k1——减压装置系数;
Gset——排放阀总的设计容量,kg/s;
pset——排放阀设计入口压力,Pa;
式中:
k2——排放阀开度;
——主冷却剂实际平均温度与目标负荷对应平均温度的差值;
△T100——排放阀开度为100%时对应的主冷却剂平均温度偏差设计值;
△T0——排放阀开度为0%时对应的主冷却剂平均温度偏差设计值;
如果小于△T0,则k2=0;如果大于△T0,则k2=1;
式中:
——蒸汽发生器蒸汽出口至排放阀前管路压降,Pa;
△pf——蒸汽发生器蒸汽出口至排放阀前沿程摩擦压降,Pa;
△pζ——蒸汽发生器蒸汽出口至排放阀前局部阻力件形阻压降,Pa;
其中△pf和△pζ的计算方法是如果管路中存在多段并联管路,只对并联管路中的一条工作支路进行计算;
式中:
△Pf,i——排放阀前第i段管道摩擦压降,Pa;
△Pζ,j——排放阀前第j个局部阻力件形阻压降,Pa;
λ——管道摩擦系数,无量纲;
ζ——局部阻力件形阻系数;
L——管道长度,m;
ρ——流体密度,kg/m3
π——圆周率;
d——管道直径,m,对局部阻力件,为接口管道直径;
G——蒸汽流量,kg/s;
i——管段号;
j——局部阻力件号;
其中:
式中:
其中
η——流体的动力粘度,kg/(m·s);
S3.计算当前时刻的蒸汽流量:
根据步骤S2所得汽轮机目标负荷值对应蒸汽流量减压装置系数k1,以及上一时刻n-1的排放阀开度k2 (n-1)、蒸汽发生器蒸汽出口至排放阀前管路压降通过以下公式计算得到当前时刻n的蒸汽流量:
t——时间,s;
tg——排放阀动作延迟时间,s;
——汽轮机初始负荷值对应蒸汽流量,kg/s;
ps——蒸汽发生器出口压力,Pa;
t(n)=t(n-1)+△t
其中:
△t——时间步长,s;
S4.计算给水流量:
通过以下公式计算得到给水流量:
式中:
——汽轮机初始负荷值对应给水流量,kg/s;
——汽轮机目标负荷值对应给水流量,kg/s;
tf——给水流量变化时间,s;
其中
——汽轮机初始负荷对应的给水调节阀阀位;
——汽轮机目标负荷对应的给水调节阀阀位;
△tvp——给水调节阀全行程时间,s;
S5.计算当前时刻的循环倍率:
根据蒸汽发生器静态特性,通过以下公式计算得到当前时刻n的循环倍率:
其中Gg为蒸汽流量,f1()为对静态特性数据进行拟合得到的三阶或三阶以上表达式;
S6.计算当前时刻的平均换热系数:
根据蒸汽发生器静态特性,通过以下公式计算得到当前时刻n的蒸汽发生器平均换热系数:
其中Gg为蒸汽流量,f2()为对静态特性数据进行拟合得到的三阶或三阶以上表达式;
S7.计算当前时刻的反应堆功率:
通过以下公式计算当前时刻n的反应堆功率:
式中:
t——时间,s;
tP——控制棒动作延迟时间,s;
P0——初始负荷,W;
c1——计算系数;
其中:
P100——100%负荷,W;
P90——90%负荷,W;
L100——100%负荷对应控制棒棒位,mm;
L90——90%负荷对应控制棒棒位,mm;
v——控制棒正常动作速率,mm/s;
S8.计算当前时刻的主冷却剂平均温度:
根据当前时刻n的反应堆功率P(n)和蒸汽发生器平均传热系数K(n),通过以下公式计算当前时刻n的主冷却剂平均温度:
式中:
P——反应堆功率,W;
△t——时间步长,s;
M——一回路冷却剂或结构的质量,kg;
CP——一回路冷却剂或结构的比热容,J/(kg·℃);
——主冷却剂平均温度,℃;
Ts——饱和蒸汽温度,℃;
K——蒸汽发生器平均换热系数,W/(m2·℃);
A——蒸汽发生器传热面积,m2
S9.计算当前时刻的再循环水流量和循环水流量:
根据当前时刻n的蒸汽流量Gg (n)、循环倍率N(n)和给水流量Gf (n),通过以下公式计算当前时刻n的再循环水流量Gof (n)和传热区入口流量G2i (n)
Gof (n)=(N(n)-1)Gg (n)
G2i (n)=Gof (n)+Gf (n)
其中:
Gof——再循环水流量,kg/s;
Gf——给水流量,kg/s;
G2i——传热区入口流量,kg/s;
S10.计算得到当前时刻的蒸汽温度:
通过蒸汽发生器二次侧传热区能量平衡方程,计算得到当前时刻n的蒸汽温度:
式中:
m2——二次侧传热区流体或结构的质量,kg;
cp2——二次侧传热区流体或结构的比热容,J/(kg·℃);
h2i——传热区入口焓值,J/kg;
h2g——传热区出口蒸汽焓值,J/kg;
h2f——传热区出口水焓值,J/kg;
S11.获取当前时刻的饱和蒸汽压力:
根据当前时刻n的蒸汽温度Ts (n),查阅水和水蒸气物性参数数据库,获取当前时刻的饱和蒸汽压力;
S12.获取当前时刻n的传热区出口水焓值:
根据当前时刻n的蒸汽温度Ts (n),查阅水和水蒸气物性参数数据库,获取当前时刻n的传热区出口水焓值h2of (n)
S13.计算当前时刻的汽水分离区温度:
通过汽水分离区能量平衡方程,计算得到当前时刻n的汽水分离区温度:
式中:
m3——二次侧汽水分离区流体或结构的质量,kg;
cp3——二次侧汽水分离区流体或结构的比热容,J/(kg·℃);
T3——汽水分离区的液相温度,℃;
h3f——再循环水焓值,J/kg;
S14.计算当前时刻的下降区温度:
根据当前时刻n的汽水分离区温度T3 (n)和蒸汽压力Ps (n),查阅水蒸气物性参数数据库,获取当前时刻n的再循环水焓值h3f (n),并通过下降区能量平衡方程,计算得到当前时刻n的下降区温度:
式中:
m1——二次侧下降区流体或结构的质量,kg;
cp1——二次侧下降区流体或结构的比热容,J/(kg·℃);
T1——下降区的液相温度,℃;
hf——给水焓值,J/kg;
S15.获取当前时刻的传热区入口焓值:
根据当前时刻n的下降区温度T1 (n)和蒸汽压力Ps (n),查阅水蒸气物性参数数据库,获取当前时刻n的传热区入口焓值h2i (n)
S16.计算当前时刻的蒸汽发生器一次侧换热系数:
根据蒸汽发生器静态特性,通过以下公式计算当前时刻n的蒸汽发生器一次侧换热系数;
其中Gg为蒸汽流量,f3()为对静态特性数据进行拟合得到的三阶或三阶以上表达式;
S17.计算当前时刻的反应堆出口冷却剂温度:
根据当前时刻n的主冷却剂平均温度T(n)、蒸汽发生器平均传热系数K(n)和蒸汽发生器一次侧换热系数α1 (n),通过下列方程组联立计算得到当前时刻n的反应堆出口冷却剂温度
反应堆出口温度计算表达式:
式中:
α1——蒸汽发生器一次侧换热系数,W/(m2·℃);
——主冷却剂***热段温度,℃;
——主冷却剂***冷段温度,℃;
S18.重复执行步骤S2至步骤S17:
令n=n+1,重复执行步骤S2至步骤S17,直到反应堆功率降至Pn时为止;
其中:
S19.确定反应堆出口冷却剂温度和蒸汽压力的上限值:
根据以上计算结果,得到冷却剂温度和蒸汽压力随时间变化的曲线,上述两条曲线的最大值分别为冷却剂温度和蒸汽压力的上限值,如果蒸汽压力和反应堆冷却剂温度的上限值均不超过设计允许值,则设计参数满足要求,否则调整排放容量重新计算。
本发明基于蒸汽发生器一次侧能量平衡方程,得到一次侧主冷却剂的温度表达式,通过上一时刻的蒸汽温度和主冷却剂平均温度,以及此刻的反应堆功率,求得此刻的主冷却剂平均温度;将蒸汽发生器内二次侧分为三个区域,分别为下降区、传热区和分离区,建立蒸汽发生器一、二次侧传热方程及二次侧各分区能量守恒方程组,将主冷却剂平均温度带入,并对这些方程进行时间离散化,得到关于二次侧传热区、下降区和汽水分离区的三个线性方程组,联立求解得到汽水分离区、下降区和传热区内蒸汽和水的温度。
采用本发明,可计算蒸汽排放***不同设计排放负荷下蒸汽压力和一次侧冷却剂温度最大值,作为判别蒸汽排放***参数设计能否满足主冷却剂***和二回路***安全运行的依据,也可指导蒸汽排放***的设计。本发明简单易用。
具体实施方式
以下结合具体实施例对本发明作进一步的详细描述,但该实施例不应理解为对本发明的限制。
本实施例基于以下假设:不考虑慢化剂和燃料的温度负反馈效应;功率与控制棒棒位呈线性相关;考虑控制棒延迟效应;主冷却剂与相邻结构温度保持一致;采用集总参数计算主冷却剂平均温度;将蒸汽发生器数量集总为一台,集总参数包括流量、流体与结构质量、换热面积;变负荷后给水流量随时间线性变化,达到目标负荷对应的流量后保持不变;变负荷后考虑排放阀信号与动作延迟,延迟时间内蒸汽流量随时间线性变化;蒸汽发生器循环倍率和平均换热系数表示为蒸汽发生器静态特性下蒸汽流量的函数。
本实施例中,蒸汽压力的设计允许值4.8MPa,反应堆冷却剂温度的设计允许值279.0℃,蒸汽发生器的传热面积856.82m2,一次侧总热容1.23×105kJ/℃(包括冷却剂、压力容器、主管道及蒸汽发生器一次侧结构),蒸汽发生器二次侧传热区总热容1.01×104kJ/℃(包括水、蒸汽以及蒸汽发生器二次侧结构),分离区总热容2.67×104kJ/℃,下降区总热容5.3×103kJ/℃,控制棒动作延迟时间1s,排放阀动作延迟时间0.5s,给水流量变化时间5s,给水温度142℃,给水压力4.2MPa,辅机蒸汽流量8.28kg/s。时间步长选取0.1s,初始满负荷运行反应堆功率100MW,满负荷运行时蒸汽发生器出口蒸汽压力3.9MPa,蒸汽流量和给水流量均为45kg/s,冷却剂平均温度271℃;目标负荷值为55%。
按照本发明提供的步骤进行,每一时刻下主要参数的计算结果如表1所示:
表1
本实施例选取了较小的时间步长,具有较高的计算精度,数据量较大,从上表可见设计排放负荷为55%时,蒸汽压力上限值可达到4.62MPa,反应堆冷却剂温度最大值为274.54℃,均未超过设计允许值。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
本说明书中未作详细描述的内容,属于本专业技术人员公知的现有技术。

Claims (1)

1.一种核动力装置蒸汽排放***参数的计算方法,包括如下步骤:
S1.初始化:
S1.1.令n=1,得到0时刻的反应堆功率P0,蒸汽发生器平均换热系数K0、循环倍率N0,蒸汽发生器二次侧传热区的出口水流量Gof 0、出口水焓值h2f 0,传热区出口蒸汽流量Gg 0、出口蒸汽焓值h2g 0,传热区入口流量G2i 0、入口焓值h2i 0,汽水分离区的温度T3 0、再循环水焓值h3f 0,给水流量Gf 0,下降区的温度T1 0
S1.2.获取0时刻的主冷却剂平均温度
通过以下公式计算得到0时刻的主冷却剂平均温度
<mrow> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mn>0</mn> </msup> <mo>=</mo> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>&amp;Delta;T</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>&amp;Delta;T</mi> <mi>m</mi> </msub> </mrow>
式中:
——反应堆初始功率P0对应的主冷却剂***平均温度,℃;
ΔTc——反应堆冷却剂控制死区,℃;
ΔTm——反应堆冷却剂温度测量误差,℃;
ΔTc由控制***设计给出,ΔTm由温度测量仪表选型厂家给出;
S1.3.获取0时刻的蒸汽温度Ts 0
通过以下公式得到0时刻的蒸汽压力ps 0
ps 0=ps0+Δps(ΔTc+ΔTm)
式中:
ps0——主冷却剂设计平均温度下的蒸汽发生器出口蒸汽压力,Pa;
Δps——反应堆功率为P0时主冷却剂平均温度较设计值升高1℃时对应的蒸汽压力增加值,Pa;
查水蒸气物性参数表,得到与ps 0对应的饱和蒸汽温度,即0时刻的蒸汽温度Ts 0
S2.计算蒸汽流量、减压装置系数、上一时刻排放阀开度和蒸汽出口至排放阀前管路压降:
根据蒸汽排放***减温减压装置设计容量和压力,以及主冷却剂平均温度偏差,通过以下公式计算得到汽轮机目标负荷值对应蒸汽流量减压装置系数k1、上一时刻n-1的排放阀开度k2 (n-1)和蒸汽发生器蒸汽出口至排放阀前管路压降
<mrow> <msubsup> <mi>G</mi> <mi>g</mi> <mi>&amp;infin;</mi> </msubsup> <mo>=</mo> <mfrac> <msup> <mi>P</mi> <mi>&amp;infin;</mi> </msup> <msub> <mi>P</mi> <mn>100</mn> </msub> </mfrac> <msub> <mi>G</mi> <mn>100</mn> </msub> </mrow>
式中:
——汽轮机目标负荷值对应蒸汽流量,kg/s;
P——目标负荷值,W;
G100——100%负荷对应蒸汽流量,kg/s;
<mrow> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>G</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> <msub> <mi>p</mi> <mrow> <mi>s</mi> <mi>e</mi> <mi>t</mi> </mrow> </msub> </mfrac> </mrow>
式中:
k1——减压装置系数;
Gset——排放阀总的设计容量,kg/s;
pset——排放阀设计入口压力,Pa;
<mrow> <msubsup> <mi>k</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msub> <mi>&amp;Delta;T</mi> <mn>0</mn> </msub> </mrow> <mrow> <msub> <mi>&amp;Delta;T</mi> <mn>100</mn> </msub> <mo>-</mo> <msub> <mi>&amp;Delta;T</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
式中:
k2——排放阀开度;
——主冷却剂实际平均温度与目标负荷对应平均温度的差值;
ΔT100——排放阀开度为100%时对应的主冷却剂平均温度偏差设计值;
ΔT0——排放阀开度为0%时对应的主冷却剂平均温度偏差设计值;
如果小于ΔT0,则k2=0;如果大于ΔT0,则k2=1;
<mrow> <mi>&amp;Delta;</mi> <msup> <mover> <mi>p</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msubsup> <mi>&amp;Delta;p</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;Delta;p</mi> <mi>&amp;zeta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> </mrow>
式中:
——蒸汽发生器蒸汽出口至排放阀前管路压降,Pa;
Δpf——蒸汽发生器蒸汽出口至排放阀前沿程摩擦压降,Pa;
Δpζ——蒸汽发生器蒸汽出口至排放阀前局部阻力件形阻压降,Pa;
其中Δpf和Δpζ的计算方法是如果管路中存在多段并联管路,只对并联管路中的一条工作支路进行计算;
<mrow> <msubsup> <mi>&amp;Delta;p</mi> <mrow> <mi>f</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfrac> <mrow> <mn>8</mn> <msubsup> <mi>&amp;lambda;</mi> <mi>i</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <msub> <mi>L</mi> <mi>i</mi> </msub> <msup> <mrow> <mo>(</mo> <msubsup> <mi>G</mi> <mi>i</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <msubsup> <mi>&amp;rho;</mi> <mi>i</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msubsup> <mi>d</mi> <mi>i</mi> <mn>5</mn> </msubsup> </mrow> </mfrac> </mrow>
<mrow> <msubsup> <mi>&amp;Delta;p</mi> <mrow> <mi>&amp;zeta;</mi> <mo>,</mo> <mi>j</mi> </mrow> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfrac> <mrow> <mn>8</mn> <msub> <mi>&amp;zeta;</mi> <mi>j</mi> </msub> <msup> <mrow> <mo>(</mo> <msubsup> <mi>G</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <msubsup> <mi>&amp;rho;</mi> <mi>j</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <msup> <mi>&amp;pi;</mi> <mn>2</mn> </msup> <msubsup> <mi>d</mi> <mi>j</mi> <mn>4</mn> </msubsup> </mrow> </mfrac> </mrow>
式中:
ΔPf,i——排放阀前第i段管道摩擦压降,Pa;
ΔPζ,j——排放阀前第j个局部阻力件形阻压降,Pa;
λ——管道摩擦系数,无量纲;
ζ——局部阻力件形阻系数;
L——管道长度,m;
ρ——流体密度,kg/m3
π——圆周率;
d——管道直径,m,对局部阻力件,为接口管道直径;
G——蒸汽流量,kg/s;
i——管段号;
j——局部阻力件号;
其中:
<mrow> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mfrac> <mn>64</mn> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> </mfrac> </mtd> <mtd> <mrow> <mn>0</mn> <mo>&lt;</mo> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>&lt;</mo> <mn>2320</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>0.3164</mn> <msup> <mrow> <mo>(</mo> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mn>0.25</mn> </msup> </mfrac> </mtd> <mtd> <mrow> <mn>2320</mn> <mo>&amp;le;</mo> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>&amp;le;</mo> <msup> <mn>10</mn> <mn>5</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0.0032</mn> <mo>+</mo> <mfrac> <mn>0.221</mn> <msup> <mrow> <mo>(</mo> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mn>0.237</mn> </msup> </mfrac> </mrow> </mtd> <mtd> <mrow> <msup> <mn>10</mn> <mn>5</mn> </msup> <mo>&amp;le;</mo> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>&amp;le;</mo> <mn>3</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mn>6</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中:
<mrow> <msup> <mi>Re</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mn>4</mn> <mrow> <mi>&amp;pi;</mi> <mi>d</mi> </mrow> </mfrac> <mfrac> <msup> <mi>G</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <msup> <mi>&amp;eta;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> </mfrac> </mrow>
其中
η——流体的动力粘度,kg/(m·s);
S3.计算当前时刻的蒸汽流量:
根据步骤S2所得汽轮机目标负荷值对应蒸汽流量减压装置系数k1,以及上一时刻n-1的排放阀开度k2 (n-1)、蒸汽发生器蒸汽出口至排放阀前管路压降通过以下公式计算得到当前时刻n的蒸汽流量:
<mrow> <msubsup> <mi>G</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>G</mi> <mi>g</mi> <mn>0</mn> </msubsup> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>G</mi> <mi>g</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msubsup> <mi>G</mi> <mi>g</mi> <mi>&amp;infin;</mi> </msubsup> </mrow> <msub> <mi>t</mi> <mi>g</mi> </msub> </mfrac> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&amp;le;</mo> <msub> <mi>t</mi> <mi>g</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <msubsup> <mi>k</mi> <mn>2</mn> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msubsup> <mo>-</mo> <mi>&amp;Delta;</mi> <msup> <mover> <mi>p</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>G</mi> <mi>g</mi> <mi>&amp;infin;</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&gt;</mo> <msub> <mi>t</mi> <mi>g</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中:
t——时间,s;
tg——排放阀动作延迟时间,s;
——汽轮机初始负荷值对应蒸汽流量,kg/s;
ps——蒸汽发生器出口压力,Pa;
t(n)=t(n-1)+Δt
其中:
Δt——时间步长,s;
S4.计算给水流量:
通过以下公式计算得到当前时刻n的给水流量:
<mrow> <msubsup> <mi>G</mi> <mi>f</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>G</mi> <mi>f</mi> <mn>0</mn> </msubsup> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>G</mi> <mi>f</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msubsup> <mi>G</mi> <mi>f</mi> <mi>&amp;infin;</mi> </msubsup> </mrow> <msub> <mi>t</mi> <mi>f</mi> </msub> </mfrac> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&amp;le;</mo> <msub> <mi>t</mi> <mi>f</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>G</mi> <mi>f</mi> <mi>&amp;infin;</mi> </msubsup> </mtd> <mtd> <mrow> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&gt;</mo> <msub> <mi>t</mi> <mi>f</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中:
——汽轮机初始负荷值对应给水流量,kg/s;
——汽轮机目标负荷值对应给水流量,kg/s;
tf——给水流量变化时间,s;
<mrow> <msub> <mi>t</mi> <mi>f</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msubsup> <mi>V</mi> <mi>p</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msubsup> <mi>V</mi> <mi>p</mi> <mi>&amp;infin;</mi> </msubsup> <mo>)</mo> </mrow> <msub> <mi>&amp;Delta;t</mi> <mrow> <mi>v</mi> <mi>p</mi> </mrow> </msub> </mrow>
其中
——汽轮机初始负荷对应的给水调节阀阀位;
——汽轮机目标负荷对应的给水调节阀阀位;
Δtvp——给水调节阀全行程时间,s;
S5.计算当前时刻的循环倍率:
根据蒸汽发生器静态特性,通过以下公式计算得到当前时刻n的循环倍率:
其中Gg为蒸汽流量,f1()为对静态特性数据进行拟合得到的三阶或三阶以上表达式;
S6.计算当前时刻的平均换热系数:
根据蒸汽发生器静态特性,通过以下公式计算得到当前时刻n的蒸汽发生器平均换热系数:
其中Gg为蒸汽流量,f2()为对静态特性数据进行拟合得到的三阶或三阶以上表达式;
S7.计算当前时刻的反应堆功率:
通过以下公式计算当前时刻n的反应堆功率:
<mrow> <msup> <mi>P</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msup> <mi>P</mi> <mn>0</mn> </msup> </mtd> <mtd> <mrow> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&amp;le;</mo> <msub> <mi>t</mi> <mi>P</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mi>P</mi> <mn>0</mn> </msup> <mo>-</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> </mrow> </mtd> <mtd> <mrow> <msup> <mi>t</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&gt;</mo> <msub> <mi>t</mi> <mi>P</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中:
t——时间,s;
tP——控制棒动作延迟时间,s;
P0——初始负荷,W;
c1——计算系数;
<mrow> <msub> <mi>c</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mn>100</mn> </msub> <mo>-</mo> <msub> <mi>P</mi> <mn>90</mn> </msub> </mrow> <mrow> <msub> <mi>L</mi> <mn>100</mn> </msub> <mo>-</mo> <msub> <mi>L</mi> <mn>90</mn> </msub> </mrow> </mfrac> <mi>v</mi> </mrow>
其中:
P100——100%负荷,W;
P90——90%负荷,W;
L100——100%负荷对应控制棒棒位,mm;
L90——90%负荷对应控制棒棒位,mm;
v——控制棒正常动作速率,mm/s;
S8.计算当前时刻的主冷却剂平均温度:
根据当前时刻n的反应堆功率P(n)和蒸汽发生器平均传热系数K(n),通过以下公式计算当前时刻n的主冷却剂平均温度:
<mrow> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <msup> <mi>P</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>+</mo> <msub> <mi>&amp;Sigma;MC</mi> <mi>P</mi> </msub> <mo>&amp;CenterDot;</mo> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>+</mo> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mi>A</mi> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>&amp;CenterDot;</mo> <msup> <msub> <mi>T</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> </mrow> <mrow> <msub> <mi>&amp;Sigma;MC</mi> <mi>P</mi> </msub> <mo>+</mo> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mi>A</mi> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> </mrow>
式中:
P——反应堆功率,W;
Δt——时间步长,s;
M——一回路冷却剂或结构的质量,kg;
CP——一回路冷却剂或结构的比热容,J/(kg·℃);
——主冷却剂平均温度,℃;
Ts——饱和蒸汽温度,℃;
K——蒸汽发生器平均换热系数,W/(m2·℃);
A——蒸汽发生器传热面积,m2
S9.计算当前时刻的再循环水流量和循环水流量:
根据当前时刻n的蒸汽流量Gg (n)、循环倍率N(n)和给水流量Gf (n),通过以下公式计算当前时刻n的再循环水流量Gof (n)和传热区入口流量G2i (n)
Gof (n)=(N(n)-1)Gg (n)
G2i (n)=Gof (n)+Gf (n)
其中:
Gof——再循环水流量,kg/s;
Gf——给水流量,kg/s;
G2i——传热区入口流量,kg/s;
S10.计算得到当前时刻的蒸汽温度:
通过蒸汽发生器二次侧传热区能量平衡方程,计算得到当前时刻n的蒸汽温度:
<mrow> <msup> <msub> <mi>T</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mi>A</mi> <mo>&amp;CenterDot;</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mi>T</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msup> <msub> <mi>G</mi> <mrow> <mi>o</mi> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <msup> <msub> <mi>h</mi> <mrow> <mn>2</mn> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>+</mo> <msup> <msub> <mi>G</mi> <mi>g</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <msup> <msub> <mi>h</mi> <mrow> <mn>2</mn> <mi>g</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mi>G</mi> <mrow> <mn>2</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <msup> <msub> <mi>h</mi> <mrow> <mn>2</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&amp;Sigma;m</mi> <mn>2</mn> </msub> <msub> <mi>c</mi> <mrow> <mi>p</mi> <mn>2</mn> </mrow> </msub> </mrow> </mfrac> <mo>+</mo> <msup> <msub> <mi>T</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> </mrow>
式中:
m2——二次侧传热区流体或结构的质量,kg;
cp2——二次侧传热区流体或结构的比热容,J/(kg·℃);
h2i——传热区入口焓值,J/kg;
h2g——传热区出口蒸汽焓值,J/kg;
h2f——传热区出口水焓值,J/kg;
S11.获取当前时刻的饱和蒸汽压力:
根据当前时刻n的蒸汽温度Ts (n),查阅水和水蒸气物性参数数据库,获取当前时刻的饱和蒸汽压力;
S12.获取当前时刻n的传热区出口水焓值:
根据当前时刻n的蒸汽温度Ts (n),查阅水和水蒸气物性参数数据库,获取当前时刻n的传热区出口水焓值h2f (n)
S13.计算当前时刻的汽水分离区温度:
通过汽水分离区能量平衡方程,计算得到当前时刻n的汽水分离区温度:
<mrow> <msup> <msub> <mi>T</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>&amp;CenterDot;</mo> <msup> <msub> <mi>G</mi> <mrow> <mi>o</mi> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <msup> <msub> <mi>h</mi> <mrow> <mn>2</mn> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mi>h</mi> <mrow> <mn>3</mn> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&amp;Sigma;m</mi> <mn>3</mn> </msub> <msub> <mi>c</mi> <mrow> <mi>p</mi> <mn>3</mn> </mrow> </msub> </mrow> </mfrac> <mo>+</mo> <msup> <msub> <mi>T</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> </mrow>
式中:
m3——二次侧汽水分离区流体或结构的质量,kg;
cp3——二次侧汽水分离区流体或结构的比热容,J/(kg·℃);
T3——汽水分离区的液相温度,℃;
h3f——再循环水焓值,J/kg;
S14.计算当前时刻的下降区温度:
根据当前时刻n的汽水分离区温度T3 (n)和蒸汽压力Ps (n),查阅水蒸气物性参数数据库,获取当前时刻n的再循环水焓值h3f (n),并通过下降区能量平衡方程,计算得到当前时刻n的下降区温度:
<mrow> <msup> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msup> <msub> <mi>G</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <msub> <mi>h</mi> <mi>f</mi> </msub> <mo>+</mo> <msup> <msub> <mi>G</mi> <mrow> <mi>o</mi> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <msup> <msub> <mi>h</mi> <mrow> <mn>3</mn> <mi>f</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mi>G</mi> <mrow> <mn>2</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <msup> <msub> <mi>h</mi> <mrow> <mn>2</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&amp;Sigma;m</mi> <mn>1</mn> </msub> <msub> <mi>c</mi> <mrow> <mi>p</mi> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mo>+</mo> <msup> <msub> <mi>T</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> </mrow>
式中:
m1——二次侧下降区流体或结构的质量,kg;
cp1——二次侧下降区流体或结构的比热容,J/(kg·℃);
T1——下降区的液相温度,℃;
hf——给水焓值,J/kg;
S15.获取当前时刻的传热区入口焓值:
根据当前时刻n的下降区温度T1 (n)和蒸汽压力Ps (n),查阅水蒸气物性参数数据库,获取当前时刻n的传热区入口焓值h2i (n)
S16.计算当前时刻的蒸汽发生器一次侧换热系数:
根据蒸汽发生器静态特性,通过以下公式计算当前时刻n的蒸汽发生器一次侧换热系数;
其中Gg为蒸汽流量,f3()为对静态特性数据进行拟合得到的三阶或三阶以上表达式;
S17.计算当前时刻的反应堆出口冷却剂温度:
根据当前时刻n的主冷却剂平均温度蒸汽发生器平均传热系数K(n)和蒸汽发生器一次侧换热系数α1 (n),通过下列方程组联立计算得到当前时刻n的反应堆出口冷却剂温度
<mrow> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mi>T</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> <mo>=</mo> <msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mrow> <mo>(</mo> <msup> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>o</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> </mrow>
<mrow> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <mfrac> <mrow> <msup> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>+</mo> <msup> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>o</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> </mrow> <mn>2</mn> </mfrac> </mrow>
反应堆出口温度计算表达式:
<mrow> <msup> <msub> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mi>o</mi> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>=</mo> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <mfrac> <msup> <mi>K</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mrow> <mn>2</mn> <msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> </mrow> </mfrac> <mrow> <mo>(</mo> <msup> <mover> <mi>T</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </msup> <mo>-</mo> <msup> <msub> <mi>T</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </msup> <mo>)</mo> </mrow> </mrow>
式中:
α1——蒸汽发生器一次侧换热系数,W/(m2·℃);
——主冷却剂***热段温度,℃;
——主冷却剂***冷段温度,℃;
S18.重复执行步骤S2至步骤S17:
令n=n+1,重复执行步骤S2至步骤S17,直到反应堆功率降至Pn时为止;
其中:
<mrow> <msub> <mi>P</mi> <mi>n</mi> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mi>P</mi> <mi>&amp;infin;</mi> </msup> <mo>-</mo> <msub> <mi>P</mi> <mn>100</mn> </msub> <mo>&amp;CenterDot;</mo> <mn>10</mn> <mi>%</mi> </mrow> </mtd> <mtd> <mrow> <mfrac> <msup> <mi>P</mi> <mi>&amp;infin;</mi> </msup> <msub> <mi>P</mi> <mn>100</mn> </msub> </mfrac> <mo>&amp;GreaterEqual;</mo> <mn>15</mn> <mi>%</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>P</mi> <mn>100</mn> </msub> <mo>&amp;CenterDot;</mo> <mn>5</mn> <mi>%</mi> </mrow> </mtd> <mtd> <mrow> <mfrac> <msup> <mi>P</mi> <mi>&amp;infin;</mi> </msup> <msub> <mi>P</mi> <mn>100</mn> </msub> </mfrac> <mo>&lt;</mo> <mn>15</mn> <mi>%</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
S19.确定反应堆出口冷却剂温度和蒸汽压力的上限值:
根据以上计算结果,得到冷却剂温度和蒸汽压力随时间变化的曲线,上述两条曲线的最大值分别为冷却剂温度和蒸汽压力的上限值,如果蒸汽压力和反应堆冷却剂温度的上限值均不超过设计允许值,则设计参数满足要求,否则调整排放容量重新计算。
CN201710800214.0A 2017-09-07 2017-09-07 核动力装置蒸汽排放***参数的计算方法 Expired - Fee Related CN107784156B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710800214.0A CN107784156B (zh) 2017-09-07 2017-09-07 核动力装置蒸汽排放***参数的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710800214.0A CN107784156B (zh) 2017-09-07 2017-09-07 核动力装置蒸汽排放***参数的计算方法

Publications (2)

Publication Number Publication Date
CN107784156A true CN107784156A (zh) 2018-03-09
CN107784156B CN107784156B (zh) 2020-11-03

Family

ID=61437588

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710800214.0A Expired - Fee Related CN107784156B (zh) 2017-09-07 2017-09-07 核动力装置蒸汽排放***参数的计算方法

Country Status (1)

Country Link
CN (1) CN107784156B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271686A (zh) * 2018-08-30 2019-01-25 中国船舶重工集团公司第七〇九研究所 核动力装置淡水冷却***海水入口温度上限值的计算方法
CN109783874A (zh) * 2018-12-18 2019-05-21 中国人民解放军海军工程大学 一种船用混合调节核汽轮机组耦合变工况计算模型
CN109992826A (zh) * 2019-02-20 2019-07-09 中国船舶重工集团公司第七一九研究所 一种减温减压装置的建模方法
CN111737856A (zh) * 2020-05-25 2020-10-02 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) 一种核动力装置变频给水泵转速的计算方法
CN112380713A (zh) * 2020-11-19 2021-02-19 上海交通大学 一种蒸汽发生器倒u型管金属壁温度分布估计方法
CN112417780A (zh) * 2020-11-19 2021-02-26 上海交通大学 蒸汽发生器二回路再循环水质量流量估计方法及***
CN112417681A (zh) * 2020-11-19 2021-02-26 上海交通大学 一种蒸汽发生器一二次侧对流换热系数分布的估计方法
CN112417782A (zh) * 2020-11-19 2021-02-26 上海交通大学 一种蒸汽发生器二回路工质循环倍率估计方法
CN112699523A (zh) * 2020-11-19 2021-04-23 上海交通大学 蒸汽发生器二回路工质汽化起始高度估计方法及装置
CN117558472A (zh) * 2024-01-11 2024-02-13 深圳大学 核反应堆冷却***及其冷却控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110172980A1 (en) * 2009-11-12 2011-07-14 Westinghouse Electric Company Llc Method of Modeling Steam Generator and Processing Steam Generator Tube Data of Nuclear Power Plant
CN105550394A (zh) * 2015-12-02 2016-05-04 中国船舶重工集团公司第七一九研究所 一种给水泵的建模方法
CN106528985A (zh) * 2016-11-03 2017-03-22 哈尔滨工程大学 一种核动力装置冷凝器的分区化仿真方法
CN106653113A (zh) * 2016-10-25 2017-05-10 核动力运行研究所 一种蒸汽发生器在线疲劳寿命监测装置及方法
CN206694815U (zh) * 2017-04-20 2017-12-01 邓策 蒸汽凝液可控温排放***

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110172980A1 (en) * 2009-11-12 2011-07-14 Westinghouse Electric Company Llc Method of Modeling Steam Generator and Processing Steam Generator Tube Data of Nuclear Power Plant
CN105550394A (zh) * 2015-12-02 2016-05-04 中国船舶重工集团公司第七一九研究所 一种给水泵的建模方法
CN106653113A (zh) * 2016-10-25 2017-05-10 核动力运行研究所 一种蒸汽发生器在线疲劳寿命监测装置及方法
CN106528985A (zh) * 2016-11-03 2017-03-22 哈尔滨工程大学 一种核动力装置冷凝器的分区化仿真方法
CN206694815U (zh) * 2017-04-20 2017-12-01 邓策 蒸汽凝液可控温排放***

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孔夏明 等: "蒸汽排放***蒸汽冷凝器动态特性仿真研究", 《原子能科学技术》 *
秋穗正 等: "核动力装置蒸汽发生器数值模拟计算", 《核科学与工程》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271686A (zh) * 2018-08-30 2019-01-25 中国船舶重工集团公司第七〇九研究所 核动力装置淡水冷却***海水入口温度上限值的计算方法
CN109783874A (zh) * 2018-12-18 2019-05-21 中国人民解放军海军工程大学 一种船用混合调节核汽轮机组耦合变工况计算模型
CN109783874B (zh) * 2018-12-18 2023-04-07 中国人民解放军海军工程大学 一种船用混合调节核汽轮机组耦合变工况计算模型
CN109992826A (zh) * 2019-02-20 2019-07-09 中国船舶重工集团公司第七一九研究所 一种减温减压装置的建模方法
CN111737856A (zh) * 2020-05-25 2020-10-02 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) 一种核动力装置变频给水泵转速的计算方法
CN112699523A (zh) * 2020-11-19 2021-04-23 上海交通大学 蒸汽发生器二回路工质汽化起始高度估计方法及装置
CN112417681A (zh) * 2020-11-19 2021-02-26 上海交通大学 一种蒸汽发生器一二次侧对流换热系数分布的估计方法
CN112417782A (zh) * 2020-11-19 2021-02-26 上海交通大学 一种蒸汽发生器二回路工质循环倍率估计方法
CN112417780A (zh) * 2020-11-19 2021-02-26 上海交通大学 蒸汽发生器二回路再循环水质量流量估计方法及***
CN112380713B (zh) * 2020-11-19 2022-03-22 上海交通大学 一种蒸汽发生器倒u型管金属壁温度分布估计方法
CN112417681B (zh) * 2020-11-19 2022-03-22 上海交通大学 一种蒸汽发生器一二次侧对流换热系数分布的估计方法
CN112417782B (zh) * 2020-11-19 2022-04-26 上海交通大学 一种蒸汽发生器二回路工质循环倍率估计方法
CN112699523B (zh) * 2020-11-19 2022-04-26 上海交通大学 蒸汽发生器二回路工质汽化起始高度估计方法及装置
CN112417780B (zh) * 2020-11-19 2023-03-24 上海交通大学 蒸汽发生器二回路再循环水质量流量估计方法及***
CN112380713A (zh) * 2020-11-19 2021-02-19 上海交通大学 一种蒸汽发生器倒u型管金属壁温度分布估计方法
CN117558472A (zh) * 2024-01-11 2024-02-13 深圳大学 核反应堆冷却***及其冷却控制方法
CN117558472B (zh) * 2024-01-11 2024-03-15 深圳大学 核反应堆冷却***及其冷却控制方法

Also Published As

Publication number Publication date
CN107784156B (zh) 2020-11-03

Similar Documents

Publication Publication Date Title
CN107784156A (zh) 核动力装置蒸汽排放***参数的计算方法
Lei et al. Experimental study on convection heat transfer of supercritical CO2 in small upward channels
CN102081059B (zh) 强内热源多孔介质通道强迫对流换热实验装置
CN108469744B (zh) 一种建立核电机组蒸汽发生器机理模型的方法及其***
CN106066235A (zh) 超临界水窄通道自然循环实验装置及方法
CN104458063B (zh) 一种节能型热量表流量检定装置及方法
CN106548812B (zh) 一种换热器、反应堆模拟***及其最大最小非能动运行能力的试验方法
CN106096310B (zh) 一种抽汽供热机组供热工况热力性能计算方法
CN106024079A (zh) 一种非能动余热排出循环性能测试***以及测试方法
Hao et al. Experimental research on reverse flow critical point among parallel U-tubes in SG
CN104777008B (zh) 一种电厂烟气余热利用***性能仿真测试实验装置
CN106952669A (zh) 一种熔融物堆内滞留压力容器外部冷却试验台架
CN105719708A (zh) 基于反应堆第4级降压***t型管夹带试验***及方法
CN201917536U (zh) 强内热源多孔介质通道强迫对流换热实验装置
CN204514628U (zh) 电厂排烟余热利用***动态响应仿真测试实验装置
Li et al. Experimental and numerical investigations on effect of reverse flow on transient from forced circulation to natural circulation
Hui et al. Experimental investigation of pressurized packing saturator for humid air turbine cycle
CN102285462A (zh) 一种飞机气动附件综合试验台
CN106297918B (zh) 一种超临界电加热模拟瞬态核释热的功率控制装置
CN104409118B (zh) 研究聚变堆腐蚀产物沉积的热工水力实验***及方法
CN106769159A (zh) 一种空冷器性能检测***
CN202180935U (zh) 一种飞机气动附件综合试验台
CN109187036B (zh) 一种母管制背压式汽轮机主蒸汽流量计算方法
CN107014861A (zh) 一种翅片管单管传热性能测试装置
CN112728985A (zh) 一种塔式余热高效回收利用***

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201103

Termination date: 20210907