CN113821998B - 利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法 - Google Patents
利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法 Download PDFInfo
- Publication number
- CN113821998B CN113821998B CN202111011367.XA CN202111011367A CN113821998B CN 113821998 B CN113821998 B CN 113821998B CN 202111011367 A CN202111011367 A CN 202111011367A CN 113821998 B CN113821998 B CN 113821998B
- Authority
- CN
- China
- Prior art keywords
- shell side
- steam
- pressure
- condenser
- enthalpy
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005094 computer simulation Methods 0.000 title claims abstract description 10
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 18
- 239000007791 liquid phase Substances 0.000 claims abstract description 9
- 239000007788 liquid Substances 0.000 claims abstract description 8
- 239000012808 vapor phase Substances 0.000 claims abstract description 6
- 238000004134 energy conservation Methods 0.000 claims abstract description 5
- 239000012071 phase Substances 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000000498 cooling water Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000005494 condensation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明的目的在于提供利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法,包括以下步骤:设定模型边界条件、结构参数及壳侧蒸汽压力初值;计算壳侧汽空间的平均比焓Hs和平均比容Vs;根据P0确定A(p)、B(p)、和计算f(p)和P1;判定f(p)是否满足偏差要求,否则重新查表、计算得到新的壳侧蒸汽压力;调用饱和水蒸汽表求得压力P1下的饱和蒸汽参数:饱和温度、气相焓、液相焓、气相比容和液相比容。本发明基于两相集总参数法建立凝汽器壳侧蒸汽区的动力学模型,通过气液两项的质量和能量守恒方程计算得到壳侧汽空间的平均比焓和平均比容,利用牛顿迭代法理论求解凝汽器壳侧的蒸汽压力,调用饱和水蒸汽表求得压力下的饱和蒸汽参数。
Description
技术领域
本发明涉及的是一种凝汽器仿真方法。
背景技术
凝汽器壳侧真空是设备运行的关键指标,直接影响汽轮机组的正常运行。目前的凝汽器动态仿真模型主要是将壳侧蒸汽视为理想气体,理想气体是一种假想没有粘性的气体,其分子是一些有弹性的、不具体积的质点,分子间相互没有作用力,对此简化的物理模型,不但可定性地分析气体某些热力现象,而且可定量地导出状态参数间存在的简单函数关系,因此通常利用理想气体状态方程求解凝汽器内壳侧蒸汽压力。但凝汽器壳侧蒸汽处于饱和状态,离液态不远,饱和蒸汽的比体积较气体小得多,分子本身体积不容忽略,分子间的内聚力随距离减小急剧增大,实际分子运动规律极其复杂,宏观上反映为状态参数的函数关系式繁复,凝汽器壳侧蒸汽不能看作理想气体进行压力求解。
表面式、管壳凝汽器作为热力系统冷源,在热电厂、核电厂以及蒸汽动力舰船中被广泛应用。汽轮机及其它系统排汽、疏水进入凝汽器的壳体,在绕流过冷却水管的同时,将热量传给管中的冷却水而被冷凝下来。随着蒸汽的凝结过程,工质的比容急剧缩小,因此形成真空,凝结下来的冷凝水流入热井内,内部工作原理见附图1。针对凝汽器设备的结构和运行特点,根据流体流动守恒方程,建立凝汽器热工水力动态仿真模型,兼顾仿真研究所需的实时性和仿真结果高精确度,真实的反映设备内部各主要参数随时间的变化规律、真实的流动和换热过程。
目前开发的凝汽器实时动态仿真模型主要为简化的理想气体物理模型,利用理想气体状态方程对凝汽器壳侧蒸汽压力进行求解。经查阅相关文献,还有一种凝汽器动态仿真模型的技术方案是基于非均匀和非平衡态的瞬态、两相流体模型,分别建立凝汽器壳侧蒸汽区气相和液相的质量、动量和能量6个守恒方程,并整场求解压力等变量,该模型需要将凝汽器壳侧蒸汽区分为多压力节点,可较真实的反映设备内部实际的物理过程,仿真计算结果精度较高,但求解过程较复杂,无法满足实时性的要求,一般不利用该方案进行凝汽器实时动态仿真模型的开发。
发明内容
本发明的目的在于提供利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法。
本发明的目的是这样实现的:
本发明利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法,其特征是:
1)设定模型边界条件、结构参数及壳侧蒸汽压力初值P0;
2)通过气液两项的质量和能量守恒方程计算壳侧汽空间的平均比焓Hs和平均比容Vs;
3)根据P0和下表,确定A(p)、B(p)、和/>
k=1,……,20;
4)根据下列公式计算f(p)和P1:
f(p)=A(p)+B(p)Vs-Hs
5)判定f(p)是否满足偏差要求,否则重新查表、计算得到新的壳侧蒸汽压力P1;
6)调用饱和水蒸汽表求得压力P1下的饱和蒸汽参数:饱和温度ts、气相焓hs、液相焓hl、气相比容νs和液相比容νl。
本发明的优势在于:本发明基于两相集总参数法建立凝汽器壳侧蒸汽区的动力学模型,通过气液两项的质量和能量守恒方程计算得到壳侧汽空间的平均比焓和平均比容,利用牛顿迭代法理论求解凝汽器壳侧的蒸汽压力,调用饱和水蒸汽表求得压力下的饱和蒸汽参数:饱和温度、气相焓、液相焓、气相比容和液相比容。
附图说明
图1为凝汽器工作原理图;
图2为迭代公式示意图;
图3为本发明的流程图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1-3,凝汽器壳侧蒸汽的平均比焓和平均比容的计算公式:
Hs=(1-x)hl+xhg (1)
Vs=(1-x)vl+xvg (2)
式中:x——蒸汽干度
hl——饱和水比焓,kJ/kg;
hg——饱和汽比焓,kJ/kg;
vl——饱和水比容,m3/kg;
vg——饱和汽比容,m3/kg。
将式(1)、(2)联立,消去x得到(3)式,满足该式的压力,便是凝汽器壳侧压力的解。
A(p)+B(p)Vs=Hs (3)
式中:
牛顿迭代法理论使用函数f(x)的泰勒级数的前面几项来寻找方程f(x)=0的根。牛顿迭代法是求方程根的重要方法之一,其最大的优点是在方程f(x)=0的单根附近具有平方收敛。用牛顿迭代法解非线性方程,是把非线性方程线性化的一种近似方法,用计算机计算十分方便。牛顿法的本质是“以直代曲”,首先猜测一个值x1,用它近似方程的根c,用过(x1,f(x1))点的切线y=f(x1)+f'(x1)(x-x1)近似代替曲线f(x),然后用切线方程y=f(x1)+f'(x1)(x-x1)=0的根x=x2=x1-f(x1)/f'(x1)近似代替曲线方程的根c,这样得到f(x)=0的第二个近似值,依次类推可得到迭代公式xn=xn-1-f(xn-1)/f'(xn-1),见附图2。
利用牛顿迭代法求解式(3)
令:f(p)=A(p)+B(p)Vs-Hs (4)
得:
利用国际通用的Fortran语言,编写蒸汽压力p、A(p)、B(p)、和/>的参数表(如下表),压力范围0.1~3180磅。
k=1,……,20;
具体计算流程见附图3:
1)设定模型边界条件、结构参数及壳侧蒸汽压力初值P0;
2)通过气液两项的质量和能量守恒方程计算得到壳侧汽空间的平均比焓Hs和平均比容Vs;
3)根据P0和参数表,确定A(p)、B(p)、和/>
4)根据公式(4)~(6)计算得到f(p)和P1;
5)判定f(p)是否满足偏差要求,否则重新查表、计算得到新的壳侧蒸汽压力P1;
6)调用饱和水蒸汽表求得压力P1下的饱和蒸汽参数:饱和温度ts、气相焓hs、液相焓hl、气相比容νs和液相比容νl。
Claims (1)
1.利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法,其特征是:
1)设定模型边界条件、结构参数及壳侧蒸汽压力初值P0;
2)通过气液两项的质量和能量守恒方程计算壳侧汽空间的平均比焓Hs和平均比容Vs;
A(p)+B(p)Vs=Hs
式中:
3)利用牛顿迭代法理论使用函数f(x)的泰勒级数来寻找方程f(x)=0的根,首先猜测一个值x1,用它近似方程的根c,用过(x1,f(x1))点的切线y=f(x1)+f'(x1)(x-x1)近似代替曲线f(x),然后用切线方程y=f(x1)+f'(x1)(x-x1)=0的根x=x2=x1-f(x1)/f'(x1)近似代替曲线方程的根c,得到f(x)=0的第二个近似值,依次类推可得到迭代公式xn=xn-1-f(xn-1)/f'(xn-1);
利用牛顿迭代法求解式A(p)+B(p)Vs=Hs
令:f(p)=A(p)+B(p)Vs-Hs
得:
根据P0和下表,确定A(p)、B(p)、和/>
k=1,……,20;
4)根据下列公式计算f(p)和P1:f(p)=A(p)+B(p)Vs-Hs
5)判定f(p)是否满足偏差要求,否则重新查表、计算得到新的壳侧蒸汽压力P1;
6)调用饱和水蒸汽表求得压力P1下的饱和蒸汽参数:饱和温度ts、气相焓hs、液相焓hl、气相比容νs和液相比容νl。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111011367.XA CN113821998B (zh) | 2021-08-31 | 2021-08-31 | 利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111011367.XA CN113821998B (zh) | 2021-08-31 | 2021-08-31 | 利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113821998A CN113821998A (zh) | 2021-12-21 |
CN113821998B true CN113821998B (zh) | 2024-03-29 |
Family
ID=78913916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111011367.XA Active CN113821998B (zh) | 2021-08-31 | 2021-08-31 | 利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113821998B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298866A (zh) * | 2014-09-30 | 2015-01-21 | 杭州电子科技大学 | 一种克劳斯硫磺回收过程中反应炉动态建模方法 |
CN104806230A (zh) * | 2015-02-16 | 2015-07-29 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
CN105069262A (zh) * | 2015-09-02 | 2015-11-18 | 广东电网有限责任公司电力科学研究院 | 电力仿真系统中节点水和蒸汽性质的监测方法及系统 |
CN107917089A (zh) * | 2017-11-17 | 2018-04-17 | 上海交通大学 | 循环水泵泵叶开度优化调节方法 |
CN108170997A (zh) * | 2018-02-07 | 2018-06-15 | 中国人民武装警察部队学院 | 一种立式固定顶油罐的预警方法及系统 |
CN109635499A (zh) * | 2018-12-29 | 2019-04-16 | 华润电力技术研究院有限公司 | 汽包建模方法、装置和计算机设备 |
CN110080884A (zh) * | 2018-10-31 | 2019-08-02 | 南京航空航天大学 | 一种涡扇发动机热端虚拟传感器信号生成与气路故障诊断方法 |
CN111159875A (zh) * | 2019-12-25 | 2020-05-15 | 东南大学 | 一种电站凝汽器壳侧和管侧动态耦合数学模型及建模方法 |
CN111652418A (zh) * | 2020-05-25 | 2020-09-11 | 西安西热节能技术有限公司 | 一种火电机组动态精细化复合参数滑压曲线生成方法 |
CN111680458A (zh) * | 2020-06-03 | 2020-09-18 | 西安交通大学 | 一种适用于钠水直流蒸汽发生器的热工水力瞬态计算方法 |
CN112613158A (zh) * | 2020-11-26 | 2021-04-06 | 中国核电工程有限公司 | 一种严重事故下安全壳内控制体热力学响应综合分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120115709A1 (en) * | 2010-11-09 | 2012-05-10 | Sienna Technologies, Inc. | High Temperature Catalysts for Decomposition of Liquid Monopropellants and Methods for Producing the Same |
-
2021
- 2021-08-31 CN CN202111011367.XA patent/CN113821998B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104298866A (zh) * | 2014-09-30 | 2015-01-21 | 杭州电子科技大学 | 一种克劳斯硫磺回收过程中反应炉动态建模方法 |
CN104806230A (zh) * | 2015-02-16 | 2015-07-29 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
CN105069262A (zh) * | 2015-09-02 | 2015-11-18 | 广东电网有限责任公司电力科学研究院 | 电力仿真系统中节点水和蒸汽性质的监测方法及系统 |
CN107917089A (zh) * | 2017-11-17 | 2018-04-17 | 上海交通大学 | 循环水泵泵叶开度优化调节方法 |
CN108170997A (zh) * | 2018-02-07 | 2018-06-15 | 中国人民武装警察部队学院 | 一种立式固定顶油罐的预警方法及系统 |
CN110080884A (zh) * | 2018-10-31 | 2019-08-02 | 南京航空航天大学 | 一种涡扇发动机热端虚拟传感器信号生成与气路故障诊断方法 |
CN109635499A (zh) * | 2018-12-29 | 2019-04-16 | 华润电力技术研究院有限公司 | 汽包建模方法、装置和计算机设备 |
CN111159875A (zh) * | 2019-12-25 | 2020-05-15 | 东南大学 | 一种电站凝汽器壳侧和管侧动态耦合数学模型及建模方法 |
CN111652418A (zh) * | 2020-05-25 | 2020-09-11 | 西安西热节能技术有限公司 | 一种火电机组动态精细化复合参数滑压曲线生成方法 |
CN111680458A (zh) * | 2020-06-03 | 2020-09-18 | 西安交通大学 | 一种适用于钠水直流蒸汽发生器的热工水力瞬态计算方法 |
CN112613158A (zh) * | 2020-11-26 | 2021-04-06 | 中国核电工程有限公司 | 一种严重事故下安全壳内控制体热力学响应综合分析方法 |
Non-Patent Citations (1)
Title |
---|
水声信号多普勒参数估计及其应用研究;徐灵基;中国博士学位论文全文数据库 工程科技II辑(第08期);C036-2 * |
Also Published As
Publication number | Publication date |
---|---|
CN113821998A (zh) | 2021-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108691585B (zh) | 一种凝汽式汽轮机低压缸效率的计算方法 | |
CN104834773A (zh) | 一种直管式直流蒸汽发生器换热性能的仿真方法 | |
CN108613565A (zh) | 一种干湿联合冷却系统的背压计算方法 | |
Alsagri et al. | Thermodynamic analysis and multi-objective optimizations of a combined recompression sCO2 brayton cycle: tCO2 rankine cycles for waste heat recovery | |
CN106096310A (zh) | 一种抽汽供热机组供热工况热力性能计算方法 | |
CN114357380A (zh) | 一种基于实时测量数据的热偏差计算方法 | |
CN105303037A (zh) | 一种引起凝汽器真空恶化的临界漏空气量的计算方法 | |
CN113821998B (zh) | 利用牛顿迭代法进行凝汽器实时动态仿真模型壳侧压力求解的方法 | |
Alsagri et al. | Performance comparison and parametric analysis of sCO2 power cycles configurations | |
CN104483152B (zh) | 非再热回热复合循环机组的热耗率测定方法 | |
Eller et al. | Techno-economic analysis of novel working fluid pairs for the Kalina cycle | |
Laskowski et al. | Determining steam condensation pressure in a power plant condenser in off-design conditions | |
Xi et al. | Optimal design of large scale dry cooling tower with consideration of off-design operation | |
Tong et al. | Dynamic simulation of an ICE-ORC combined system under various working conditions | |
CN109783904A (zh) | 一种宽参数范围二氧化碳物性求解方法 | |
CN109187036B (zh) | 一种母管制背压式汽轮机主蒸汽流量计算方法 | |
CN112100751A (zh) | 一种抽凝机组背压变化对机组功率影响的计算方法及系统 | |
Jinrong et al. | Transient Simulation Modeling of Condenser Based on Thermal System Dynamics | |
CN113076632B (zh) | 饱和蒸汽发电机组汽轮机排汽焓软测量方法 | |
Wu et al. | Dynamic Operational Characteristics for Condenser and Correlative Equipment | |
CN113049279B (zh) | 一种汽液分离式中高温地热流体实验测试系统 | |
Zhou et al. | Simulation of Waste Heat Recovery From Ship Boiler Exhaust Gas | |
CN111396847A (zh) | 一种锅炉的汽温控制方法 | |
Abril et al. | Analysis of the energetic and exergetic performance in condensers of thermal generation plant | |
CN118211385A (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 |