压水堆核电站汽轮机循环的热效率

压水堆核电站热效率计算

压水堆核电站汽轮机循环的热效率

作者:Valery Ochkov (http://twt.mpei.ac.ru/ochkov/v_ochkov.htm)

引言

本应用程序计算压水堆核电站中汽轮机循环的热效率。

加压水从核反应堆(核电站一回路)流至蒸汽发生器,并将热量传递给二回路的沸水。饱和蒸汽从蒸汽发生器供应至汽轮机的高压缸。来自蒸汽发生器的新蒸汽的一部分进入过热器。在过热之前,来自高压缸的排汽在分离器中被部分干燥,分离出部分水分,然后绕过汽轮机的低压缸返回回路。

核电站的热力回路引入了一个化石燃料发电厂所没有的风险因素。在核电站中,如果将饱和蒸汽供给出口压力为5-4 kPa的汽轮机,则汽轮机末级的湿度将高得不可接受;这可能导致汽轮机故障。

水性质

开始

restart:

with(ThermophysicalData):

with(Units[Standard]):

水/蒸汽性质

作为压力和温度函数的水或蒸汽的比焓

wspTSP:=(P)→ Property(temperature,pressure=P,Q=1,water):

作为温度函数的饱和线上蒸汽的比焓

wspHSST:=(T)→Property(enthalpy,temperature=T,Q=1,water):

作为温度函数的饱和线上蒸汽的比熵

wspSSST:=(T)→Property(entropy,temperature=T,Q=1,water):

作为温度函数的饱和线上水的比焓

wspHSWT:=(T)→Property(enthalpy,temperature=T,Q=0,water):

作为压力函数的饱和线上水的比熵

wspSSWP:=(P)→Property(entropy,pressure=P,Q=0,water):

作为压力函数的饱和线上水的比焓

wspHSWP:=(P)→Property(enthalpy,pressure=P,Q=0,water):

作为压力和比焓函数的蒸汽干度

wspXSPH:=(P, H)→ Property(Q, pressure=P, enthalpy=H, water):

作为压力和温度函数的水或蒸汽的比焓

wspHPT:=(P,T)→Property(enthalpy,pressure=P,temperature=T,water):

作为压力和温度函数的水或蒸汽的比熵

wspSPT:=(P,T)→Property(entropy,pressure=P,temperature=T,water):

作为压力和温度函数的水或蒸汽的比容

wspXPS:=(P,S)→Property(Q,pressure=P,entropy=S,water):

作为温度和干度函数的水或蒸汽的比焓

wspHSTX:=(T,X)→Property(enthalpy,temperature=T,Q=X,water):

作为压力和比熵函数的水或蒸汽的温度

wspTPS:=(P,S)→Property(temperature, pressure=P,entropy=S,water):

作为压力和比熵函数的水或蒸汽的比焓

wspHPS:=(P,S)→Property(enthalpy,pressure=P,entropy=S,water):

作为压力和比焓函数的水或蒸汽的温度

wspTPH:=(P,H)→Property(temperature, pressure=P,enthalpy=H,water):

用于蒸汽膨胀和水压缩过程的函数

wspHEXPANSIONPPEFF:= proc( p0, p1, eff)
#作为初始点压力p0(饱和线)、终点压力p1和涡轮内效率eff的函数,蒸汽膨胀过程结束时蒸汽的比焓L

local t0, h0, s0, s1, t1, x1, h1;

t0:= wspTSP(p0):

h0:= wspHSST(t0):

s0:= wspSSST(t0):

s1:= s0:

t1:= wspTSP(p1):

x1:= wspXPS(p1,s1):

h1:= wspHSTX(t1,x1):

return h0-(h0-h1)·eff

end proc:

wspHEXPANSIONPTPEFF:= proc( p0, t0, p1, eff)

#作为初始点压力p0(饱和线)、终点压力p1和泵内效率eff的函数,压缩过程结束时水的比焓

local h0, s0, s1, t1, x1, h1;

h0:= wspHPT(p0, t0):

s0:= wspSPT(p0, t0):

s1:= s0:

t1:= wspTPS(p1, s1):

if Property(PhaseString, pressure=p1, entropy=s1, water)="twophase"

then

x1:= wspXPS(p1, s1): h1:= wspHSTX(t1, x1)

else

h1:= wspHPT(p1, t1)

end if:

return h0-(h0-h1)·eff

end proc:

wspHCOMPRESSIONPPEFF:= proc(p0, p1, eff)

#作为初始点压力p0(饱和线)、终点压力p1和泵内效率eff的函数,压缩过程结束时水的比焓

local h0, s0, s1, h1;

h0:= wspHSWP(p0):

s0:= wspSSWP(p0):

s1:= s0:

h1:= wspHPS(p1, s1):

return h0 + (h1 - h0)/eff

end proc:

wspHCOMPRESSIONPTPEFF:= proc( p0, t0, p1, eff)

#作为初始点压力p0、初始点温度t0、终点压力p1和泵内效率eff的函数,压缩过程结束时水的比焓

local h0, s0, s1, h1;

h0:= wspHPT(p0, t0):

s0:= wspSPT(p0, t0):

s1:= s0:

h1:= wspHPS(p1, s1):

return h0 + (h1 - h0)/eff

end proc:

输入数据

p_1 := 6.2 MPa: # 高压缸入口饱和蒸汽压力

p_2 := 0.5 MPa: # 分离器入口湿蒸汽压力

p_4 := 4 kPa: # 冷凝器压力

t_3 := 250 deg C: # 低压缸入口蒸汽温度

p_10 := (1+0.05)·p_1: # 给水压力

η_i_HPC := 0.89: # 高压缸内相对效率

η_i_LPC := 0.86: # 低压缸内效率

η_i_p := 0.75: # 泵内效率

计算

高压缸入口蒸汽参数

t_1 := wspTSP(p_1): t_1 - 273.15 K;

277.7 °C (4.1)

h_1 := wspHSST(t_1);

2782.37 kJ/kg (4.2)

来自过热器水的比焓

hw_1 := wspHSWT(t_1);

1225.06 kJ/kg (4.3)

分离器入口蒸汽参数

t_2 := wspTSP(p_2): t_2 - 273.15 K;

151.8 °C (4.4)

hw_2 := wspHSWT(t_2);

640.1 kJ/kg (4.5)

s_2 := s_1: # 假设等熵?

h_2 := wspHEXPANSIONPPEFF(p_1, p_2, η_i_HPC);

2394.01 kJ/kg (4.6)

x_2 := wspXSPH(p_2, h_2);

83.20 % (4.7)

低压缸入口参数

p_3 := p_2:

h_3 := wspHPT(p_3, t_3);

2961.04 kJ/kg (4.8)

蒸汽在低压缸中膨胀过程终点的参数

t_4 := wspTSP(p_4): t_4 - 273.15 K;

28.96 °C (4.9)

h_4 := wspHEXPANSIONPTPEFF(p_3, t_3, p_4, η_i_LPC);

2298.68 kJ/kg (4.10)

x_4 := wspXSPH(p_4, h_4);

89.52 % (4.11)

冷凝器出口参数

h_5 := wspHSWP(p_4);

121.39 kJ/kg (4.12)

分离器中的抽汽份额

alpha_SH := solve(h_1*alpha_SH + h_2 = (1 - x_2)*hw_2 + alpha_SH*hw_1 + x_2*h_3, alpha_SH);

p_5 := p_4:

冷凝泵1出口水的比焓

h_6 := wspHCOMPRESSIONPPEFF(p_5, p_2, η_i_p);

122.01 kJ/kg (4.14)

冷凝泵2入口水的比焓

h_7 := solve(x_2·h_5 + (1-x_2)·hw_2 = 1·h_7, h_7);

208.52 kJ/kg (4.15)

冷凝泵2入口水温度

t_7 := wspTPH(p_2, h_7): t_7 - 273.15 K;

49.70 °C (4.16)

冷凝泵2出口水的比焓

h_8 := wspHCOMPRESSIONPTPEFF(p_2, t_7, p_1, η_i_p);

215.72 kJ/kg (4.17)

冷凝泵2出口水温度

t_8 := wspTPH(p_1, h_8): t_8 - 273.15 K;

50.25 °C (4.18)

给水泵入口水的比焓

h_9 := solve(α_SH · hw_1 + 1 · h_8 = (1 + α_SH) · h_9, h_9);

318.81 kJ/kg (4.19)

给水泵入口水温度

t_9 := wspTPH(p_1, h_9): t_9 - 273.15 K;

74.95 °C (4.20)

给水的比焓

h_10 := wspHCOMPRESSIONPTPEFF(p_1, t_9, p_10, η_i_p);

319.21 kJ/kg (4.21)

泵的比功

w_p := x_2 · (h_6 - h_5) + 1 · (h_8 - h_7) + (1 + α_SH) · (h_10 - h_9);

8.160 kJ/kg (4.22)

高压缸的比功

w_HPC := 1 · (h_1 - h_2);

388.4 kJ/kg (4.23)

低压缸的比功

w_LPC := x_2 · (h_3 - h_4);

551.1 kJ/kg (4.24)

蒸汽发生器的比热

q_SG := (1 + α_SH) · (h_1 - h_7);

2866.7 kJ/kg (4.25)

结论

因此,热效率为

η_T := (w_HPC + w_LPC - w_p) / q_SG;

32.49 %

Maple文件下载:https://www.maplesoft.com/applications/Detail.aspx?id=154033

【事件触发一致性】研究多智能体网络如何通过分布式事件驱动控制实现有限时间内的共识(Matlab代码实现)内容概要:本文围绕多智能体网络中的事件触发一致性问题,研究如何通过分布式事件驱动控制实现有限时间内的共识,并提供了相应的Matlab代码实现方案。文中探讨了事件触发机制在降低通信负担、提升系统效率方面的优势,重点分析了多智能体系统在有限时间收敛的一致性控制策略,涉及系统模型构建、触发条件设计、稳定性与收敛性分析等核心技术环节。此外,文档还展示了该技术在航空航天、电力系统、机器人协同、无人机编队等多个前沿领域的潜在应用,体现了其跨学科的研究价值和工程实用性。; 适合人群:具备一定控制理论基础和Matlab编程能力的研究生、科研人员及从事自动化、智能系统、多智能体协同控制等相关领域的工程技术人员。; 使用场景及目标:①用于理解和实现多智能体系统在有限时间内达成一致的分布式控制方法;②为事件触发控制、分布式优化、协同控制等课题提供算法设计与仿真验证的技术参考;③支撑科研项目开发、学术论文复现及工程原型系统搭建; 阅读建议:建议结合文中提供的Matlab代码进行实践操作,重点关注事件触发条件的设计逻辑与系统收敛性证明之间的关系,同时可延伸至其他应用场景进行二次开发与性能优化。
【四旋翼无人机】具备螺旋桨倾斜机构的全驱动四旋翼无人机:建模与控制研究(Matlab代码、Simulink仿真实现)内容概要:本文围绕具备螺旋桨倾斜机构的全驱动四旋翼无人机展开,重点研究其动力学建模与控制系统设计。通过Matlab代码与Simulink仿真实现,详细阐述了该类无人机的运动学与动力学模型构建过程,分析了螺旋桨倾斜机构如何提升无人机的全向机动能力与姿态控制性能,并设计相应的控制策略以实现稳定飞行与精确轨迹跟踪。文中涵盖了从系统建模、控制器设计到仿真验证的完整流程,突出了全驱动结构相较于传统四旋翼在欠驱动问题上的优势。; 适合人群:具备一定控制理论基础和Matlab/Simulink使用经验的自动化、航空航天及相关专业的研究生、科研人员或无人机开发工程师。; 使用场景及目标:①学习全驱动四旋翼无人机的动力学建模方法;②掌握基于Matlab/Simulink的无人机控制系统设计与仿真技术;③深入理解螺旋桨倾斜机构对飞行性能的影响及其控制实现;④为相关课题研究或工程开发提供可复现的技术参考与代码支持。; 阅读建议:建议读者结合提供的Matlab代码与Simulink模型,逐步跟进文档中的建模与控制设计步骤,动手实践仿真过程,以加深对全驱动无人机控制原理的理解,并可根据实际需求对模型与控制器进行修改与优化。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值