CN112178621B - 一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法 - Google Patents
一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法 Download PDFInfo
- Publication number
- CN112178621B CN112178621B CN202011108877.4A CN202011108877A CN112178621B CN 112178621 B CN112178621 B CN 112178621B CN 202011108877 A CN202011108877 A CN 202011108877A CN 112178621 B CN112178621 B CN 112178621B
- Authority
- CN
- China
- Prior art keywords
- evaporation
- steam
- deviation
- heat exchange
- degree
- 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
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F22—STEAM GENERATION
- F22G—SUPERHEATING OF STEAM
- F22G3/00—Steam superheaters characterised by constructional features; Details of component parts thereof
- F22G3/008—Protection of superheater elements, e.g. cooling superheater tubes during starting-up periods, water tube screens
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F22—STEAM GENERATION
- F22G—SUPERHEATING OF STEAM
- F22G5/00—Controlling superheat temperature
- F22G5/20—Controlling superheat temperature by combined controlling procedures
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法,根据换热设备的结构尺寸、流量、进出口温度和运行压力建立设备的热质平衡计算模型,根据经典蒸干质量含汽率计算关联式计算蒸干传热恶化现象发生时的质量含汽率;根据换热设备实际运行过程,通过实验或数值模拟确定其实际质量含汽率分布规律;解析偏离热力平衡传热机理,定义衡量蒸干后偏离热力平衡程度的量化指标‑偏离度,并提出偏离度与缺液区蒸汽过热温度的计算方法;根据换热设备实际运行状态,得到蒸干传热恶化现象发生后的偏离热力平衡程度与蒸汽过热温度。本发明为缺液区蒸汽过热温度的预测以及涉及蒸干传热恶化现象的汽液两相换热设备热力计算合理分区提供指导。
Description
技术领域
本发明涉及一种计算方法,尤其涉及一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法。
背景技术
流动沸腾是工业用换热器、核反应堆棒束、核动力蒸汽发生器中的一种常见现象。为提高工质作功能力,通常增大换热设备出口质量含汽率或使蒸汽在出口处达到过热状态。但是当流道内质量含汽率达到相应运行状态下的临界值时将发生沸腾危机-蒸干传热恶化现象。蒸干发生在流动沸腾环状流向雾状流转变的区域。在环状流域,随着流动沸腾的进行,质量含汽率逐渐增大,壁面处连续液膜在液滴夹带和沉积、自身蒸发及表面扰动等多重作用下逐渐变薄,最终当质量含汽率达到临界值时液膜被撕裂成液滴、液膜消失,壁面处直接被蒸汽覆盖。由于相同状态下蒸汽的传热性能远低于液相,因此液膜的消失意味着壁面处传热性能的急剧下降、壁温的急剧升高。蒸干发生后,传热分区转变为缺液区。该区域内连续蒸汽夹带着饱和离散液滴进行流动与换热。尽管该区域仍然处于流动沸腾区,但是由于壁面温度较高,并且换热方式为壁面与蒸汽间的对流换热,因此出现偏离热力平衡现象,如图1所示。可以发现蒸干发生后壁温发生显著的非线性变化,蒸干位置附近出现高应力区域(即应力集中)。传热管束长期在这种环境下工作时,该区域内首先会出现微裂纹,裂纹尖端位置的高度应力集中促使裂纹逐渐扩展,传热管变薄,最终在偶然的超载瞬间,传热管出现疲劳破坏(断裂),从而导致一、二次侧流体的混合,进而破坏反应堆的工作状态。在核反应堆、蒸汽发生器和其它换热设备设计和运行过程中确定壁温的轴向变化率十分重要。如果不对其进行有效的控制,较大的壁温轴向变化率将引发疲劳破坏、应力腐蚀和老化失效等问题,直接影响换热设备的安全和稳定运行。因此必须对蒸干后偏离热力平衡现象进行准确预测以实现有效的控制,避免壁温突然升高引起的传热管应力腐蚀、老化失效以及可能引发的安全事故,这也是对蒸干后偏离热力平衡程度量化的意义所在。
目前研究人员针对蒸干传热恶化现象主要关注的是蒸干发生的位置及蒸干发生时壁温、表面传热系数等的变化规律,对蒸干及蒸干后区域的偏离热力平衡传热关注很少。但偏离热力平衡现象对于蒸干及蒸干后壁温的控制又非常重要。因此,研究蒸干后偏离热力平衡程度及蒸汽过热温度计算方法,提出有价值的建议对直流蒸汽发生器等换热设备的安全、稳定运行具有非常重要的意义。
发明内容
本发明的目的在于提供一种蒸干后偏离热力平衡程度及缺液区蒸汽过热温度的计算方法,以实现管内流动沸腾过程中蒸干及蒸干后偏离热力平衡程度的量化及汽液两相换热设备热力计算的合理分区。
本发明的目的是这样实现的:
一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法,包括如下步骤:
步骤1:根据换热设备的结构尺寸、流量、进出口温度和运行压力建立设备的热质平衡计算模型,根据经典蒸干质量含汽率计算关联式计算蒸干传热恶化现象发生时的质量含汽率;根据换热设备实际运行过程,通过实验或数值模拟确定其实际质量含汽率分布规律;
步骤2:解析偏离热力平衡传热机理,定义衡量蒸干后偏离热力平衡程度的量化指标-偏离度,并提出偏离度与缺液区蒸汽过热温度的计算方法;
步骤3:根据换热设备实际运行状态,得到蒸干传热恶化现象发生后的偏离热力平衡程度与蒸汽过热温度;
步骤4:利用已有公开数据对比验证本发明提出的计算方法的准确性。
本发明还包括这样一些特征:
步骤1所述的热质平衡计算模型如下:
qAz=qm[(xeh"+(1-xe)h')-hin]
式中:q为壁面热流密度,W/m2;Az为从入口到任一轴向位置的换热面积,m2;qm为质量流量,kg/s;xe为热平衡含汽率;h"为饱和蒸汽焓,J/kg;h'为饱和液焓,J/kg;hin为入口焓,J/kg。
步骤1所述的蒸干质量含汽率计算如下:
式中:xcrit为蒸干质量含汽率;G为质量流速,kg/(m2·s);de为流域当量直径,m;p为运行压力,Pa。
步骤2所述的偏离度计算如下:
式中:xa为实际质量含汽率;zpost-dryout为蒸干位置到蒸干后任一位置的距离,m。
步骤2所述的缺液区蒸汽过热温度(焓值)计算如下:
Q=Ql+Qv
Ql=(xa-xcrit)(M'+M")hlv
Qv=(xa-xcrit)(M'+M")(hv-h")
式中:Q为蒸干位置到蒸干后任一位置的总换热量,W;Ql为蒸干位置到蒸干后任一位置液相汽化吸收的热量,W;Qv为蒸干位置到蒸干后任一位置蒸汽过热吸收的热量,W;hv为过热蒸汽焓值,J/kg。
与现有技术相比,本发明的有益效果是:
1、本发明定义了衡量缺液区偏离热力平衡程度的量化指标-偏离度,为流动沸腾过程中蒸干后偏离热力平衡现象的量化研究提供了很大的便捷;
2、本发明提出了一种基于上述定义的偏离度确定从缺液区向单相汽对流区过渡时的过热蒸汽焓值(温度)的方法,进而计算出更符合实际换热的缺液区和单相蒸汽对流区所需的传热面积(传热管长度);
3、本发明提出的蒸干后偏离热力平衡程度及缺液区蒸汽过热温度的计算方法可以更加准确地研究两相流动沸腾蒸干及蒸干后区域的应力集中现象,这样可以发现、进而避免由于蒸干传热恶化现象带来的传热管壁面烧毁问题。
附图说明
图1是蒸干后偏离热力平衡示意图;
图2是蒸干后偏离热力平衡程度及蒸汽过热温度计算流程图;
图3是直流蒸汽发生器实际运行状态图;
图4是直流蒸汽发生器轴向含汽率数值模拟结果图;
图5是缺液区偏离度沿传热管轴向高度的分布图;
图6是一、二次侧耦合换热下缺液区过热蒸汽温度预测图;
图7a-b是针对公开实验数据的缺液区过热蒸汽温度预测图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
为量化地评估蒸干后偏离热力平衡程度,定义新参数-偏离度post-dl(post-dryout deviation level):
上述提出的新参数-偏离度可用于预测缺液区的过热蒸汽温度。以蒸干位置为预测起始点,各参数的计算式如下:
Q=Ql+Qv (3)
Ql=(xa-xcrit)(M'+M")hlv (4)
Qv=(xa-xcrit)(M'+M")(hv-h") (5)
在公式(6)中,Q、M'、M"、h'均已知,xcrit可以根据适用于给定条件下的经典临界质量含汽率计算关联式计算得到,xe可根据热质平衡计算得到。
qAz=qm[(xeh"+(1-xe)h')-hin] (9)
这样,当给定任一位置处的实际质量含汽率或偏离热力平衡程度值时,即可根据偏离热力平衡程度的定义式(1)计算出对应的偏离度或实际质量含汽率。然后根据式(6)就可以得到该位置处过热蒸汽的焓值,进而利用水和水蒸汽参数表就可以确定相应压力和焓值下的过热蒸汽温度。
在蒸干传热恶化现象发生的位置上游,热平衡含汽率基本等于实际质量含汽率,因此不存在偏离热力平衡现象。在蒸干传热恶化现象发生后(蒸干后充分发展区域),壁面直接与蒸汽接触,此时壁面处传递的热量先全部被蒸汽吸收,连续蒸汽中夹带的饱和液滴进而从蒸汽中吸收部分热量,用于液滴本身的汽化,由于该区域内传热方式的复杂性导致偏离度发生非线性变化。当液滴全部汽化,传热进入单相蒸汽对流区时,偏离度降低到0。也就是说,当偏离度为零时,即为缺液区向单相汽对流区的过渡,此时的过热蒸汽温度(过热蒸汽焓值) 可作为划分缺液区与单相汽对流区的判据,代替直流蒸汽发生器等汽液两相换热设备热力计算中采用干饱和蒸汽焓作为判据的传统做法,为传热区域的合理划分提供依据。
热力计算是直流蒸汽发生器等汽液两相换热设备设计的关键步骤,计算时通常根据流体与壁面间换热的特点将其分为单相液对流区、核态沸腾区、液膜强制对流蒸发区、缺液区和单相汽对流区。缺液区与单相汽对流区的过渡以干饱和蒸汽的焓值为基准进行划分。然而通过前述研究发现,在实际流动沸腾过程中缺液区发生了偏离热力平衡现象,该区域内蒸汽吸收热量进入过热状态。目前暂未提出缺液区内过热蒸汽温度的确定方法,所以工程中多以干饱和蒸汽的焓值为依据进行传热区域的划分,很明显采用干饱和蒸汽焓值作为基准的分区方法带有一定的近似性。在定义衡量缺液区偏离热力平衡程度的量化指标-偏离度这一参数后,就可以根据式(3)~(9)所述的方法确定从缺液区向单相汽对流区过渡时的过热蒸汽焓值(温度),进而计算出更符合实际换热的缺液区和单相蒸汽对流区所需的传热面积(传热管长度),为直流蒸汽发生器等汽液两相换热设备传热区域的合理确定提供一定的工程参考。
图7a中公开数据是针对均匀电加热竖直管内的蒸干进行的数值模拟研究(工况:压力 7.01MPa,入口过冷度11.7K,入口质量流速1000.9kg/(m2·s),壁面热流密度765kW/m2),给出了轴向四个高度(名称分别为z2、z5、z8、z10)处蒸汽温度的径向分布,其中后三个位置均位于缺液区内;图7b中公开数据是针对非均匀热流密度下管束间的蒸干进行的数值模拟研究(工况:压力6.38MPa,入口过冷度41.82K,入口质量流速190.19kg/(m2·s),缺液区壁面热流密度约为205kW/m2),提供了缺液区内过热蒸汽温度的轴向分布数据。
本发明以换热器等设备中典型的蒸干传热恶化现象为研究对象,提出一种蒸干后偏离热力平衡程度及缺液区蒸汽过热温度的计算方法。
下面结合附图和具体实施方式进一步说明本发明技术方案。本发明以换热设备实际结构与运行状态为基础,借助流道内质量含汽率的分布规律进行蒸干后偏离热力平衡程度和蒸汽过热温度的计算,具体流程图见图2。步骤如下:
(1)确定换热设备结构及运行状态(例如附图3),通过热质平衡(公式(9))计算热平衡质量含汽率的分布(图4);
(2)确定换热设备内实际质量含汽率的轴向分布(实验测量或数值模拟,图4);
(3)通过蒸干质量含汽率计算关联式(公式(7)~(8))确定蒸干传热恶化现象发生时的质量含汽率;
(4)计算本发明定义的蒸干后偏离热力平衡程度的量化指标-偏离度(公式(1)~(2),图5);
(5)根据本发明提出的式(3)~(6)计算缺液区蒸汽过热温度(图6~7)。从图中可以看出,采用本方法能较为准确地预测所研究工况下缺液区内的蒸汽过热温度。
上述实施例并非是对本发明的限制,上述实施例和说明书中描述的只是本发明的基本原理、主要特征和本发明的优点,在不脱离本发明精神和范围的前提下,本领域的技术人员能从本发明公开的内容直接导出或联想到的所有变形,均被应认为是处于本发明的保护范围之内。
蒸干传热恶化及蒸干后缺液区偏离热力现象的发生导致壁温发生非线性变化,并且缺液区尽管仍存在饱和离散液滴,但蒸汽已经进入过热状态。通过本方法计算蒸干后缺液区偏离热力平衡程度和蒸汽过热温度分布规律,据此对换热设备传统热力计算方法进行改进,得到更符合实际换热的缺液区和单相蒸汽对流区所需的传热面积(传热管长度)。
综上所述:本发明的目的在于提出一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法。通过运用热质平衡基础理论,综合考虑偏离热力平衡传热机理,得到流动沸腾蒸干后偏离热力平衡程度及缺液区蒸汽过热温度的计算式。针对工程中变热流密度和定热流密度加热等情形,运用本发明提出的方法进行了计算,精度较高。本发明可以为缺液区蒸汽过热温度的预测以及涉及蒸干传热恶化现象的汽液两相换热设备热力计算合理分区提供指导。
Claims (1)
1.一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法,其特征在于,包括如下步骤:
步骤1:根据换热设备的结构尺寸、流量、进出口温度和运行压力建立设备的热质平衡计算模型,根据经典蒸干质量含汽率计算关联式计算蒸干传热恶化现象发生时的质量含汽率;根据换热设备实际运行过程,通过实验或数值模拟确定其实际质量含汽率分布规律;
步骤2:解析偏离热力平衡传热机理,定义衡量蒸干后偏离热力平衡程度的量化指标-偏离度,并提出偏离度与缺液区蒸汽过热温度的计算方法;
步骤3:根据换热设备实际运行状态,得到蒸干传热恶化现象发生后的偏离热力平衡程度与蒸汽过热温度;
步骤1所述的热质平衡计算模型如下:
qAz=qm[(xeh"+(1-xe)h')-hin]
式中:q为壁面热流密度,W/m2;Az为从入口到任一轴向位置的换热面积,m2;qm为质量流量,kg/s;xe为热平衡含汽率;h"为饱和蒸汽焓,J/kg;h'为饱和液焓,J/kg;hin为入口焓,J/kg;
步骤1所述的蒸干质量含汽率计算如下:
式中:xcrit为蒸干质量含汽率;G为质量流速,kg/(m2·s);de为流域当量直径,m;p为运行压力,Pa;
步骤2所述的偏离度计算如下:
式中:xa为实际质量含汽率;zpost-dryout为蒸干位置到蒸干后任一位置的距离,m;
步骤2所述的缺液区蒸汽过热温度(焓值)计算如下:
Q=Ql+Qv
Ql=(xa-xcrit)(M'+M")hlv
Qv=(xa-xcrit)(M'+M")(hv-h")
式中:Q为蒸干位置到蒸干后任一位置的总换热量,W;Ql为蒸干位置到蒸干后任一位置液相汽化吸收的热量,W;Qv为蒸干位置到蒸干后任一位置蒸汽过热吸收的热量,W;hv为过热蒸汽焓值,J/kg。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011108877.4A CN112178621B (zh) | 2020-10-16 | 2020-10-16 | 一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011108877.4A CN112178621B (zh) | 2020-10-16 | 2020-10-16 | 一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112178621A CN112178621A (zh) | 2021-01-05 |
CN112178621B true CN112178621B (zh) | 2022-08-02 |
Family
ID=73950620
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011108877.4A Active CN112178621B (zh) | 2020-10-16 | 2020-10-16 | 一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112178621B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113490391B (zh) * | 2021-05-27 | 2024-06-21 | 合肥通用机械研究院有限公司 | 电子冷却用矩形微通道单元一维温度分布计算方法和*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101216164A (zh) * | 2007-12-29 | 2008-07-09 | 西安交通大学 | 一种水冷壁在线安全评价方法 |
WO2015032614A1 (fr) * | 2013-09-03 | 2015-03-12 | Cockerill Maintenance & Ingenierie S.A. | Procede et dispositif pour prevenir l'assechement dans une chaudiere de centrale solaire a concentration de type tour |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9069532B2 (en) * | 2011-07-25 | 2015-06-30 | International Business Machines Corporation | Valve controlled, node-level vapor condensation for two-phase heat sink(s) |
-
2020
- 2020-10-16 CN CN202011108877.4A patent/CN112178621B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101216164A (zh) * | 2007-12-29 | 2008-07-09 | 西安交通大学 | 一种水冷壁在线安全评价方法 |
WO2015032614A1 (fr) * | 2013-09-03 | 2015-03-12 | Cockerill Maintenance & Ingenierie S.A. | Procede et dispositif pour prevenir l'assechement dans une chaudiere de centrale solaire a concentration de type tour |
Non-Patent Citations (3)
Title |
---|
Numerical simulation of dryout and post-dryout heat transfer in a straight-pipe once-through steam generator;Jianxin Shi;《Applied Thermal Engineering》;20160524;第135页公式(20) * |
直流蒸汽发生器全范围流动与传热数值模拟;史建新;《原子能科学技术》;20170331;第51卷(第3期);第427页 * |
直管式直流蒸汽发生器蒸干及蒸干后传热数值;史建新;《工程科技Ⅱ辑》;20210316(第04期);第2.4.3节、第3.4.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN112178621A (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104834773B (zh) | 一种直管式直流蒸汽发生器换热性能的仿真方法 | |
Davidson et al. | Studies of heat transmission through boiler tubing at pressures from 500 to 3300 pounds | |
González-Gómez et al. | Thermo-economic optimization of molten salt steam generators | |
Tao et al. | Parametric analysis of heat transfer rate of passive residual heat removal heat exchanger submerged in water tank | |
CN112178621B (zh) | 一种蒸干后偏离热力平衡程度及蒸汽过热温度的计算方法 | |
O'Donovan et al. | Pressure drop analysis of steam condensation in air-cooled circular tube bundles | |
Walter et al. | How can the heat transfer correlations for finned-tubes influence the numerical simulation of the dynamic behavior of a heat recovery steam generator? | |
Shen et al. | Simulation of dryout phenomenon and transient heat transfer performance of the once-through steam generator based on heat transfer partition | |
Köhler et al. | Influence of the wetting state of a heated surface on heat transfer and pressure loss in an evaporator tube | |
González-Gómez et al. | Steam generator design for solar towers using solar salt as heat transfer fluid | |
Panchal et al. | Experimental investigation of single-phase, condensation, and flow boiling heat transfer for a spirally fluted tube | |
Zhang et al. | Heat transfer characteristics in super-low finned-tube bundles of moisture separator reheaters | |
Shao et al. | Numerical investigation on the thermal-hydraulic performance of separated structure steam generator with different tube arrangements | |
Zhang | Local and overall condensation heat transfer behavior in horizontal tube bundles | |
Seubert et al. | Analysis of a helical coil once-through molten salt steam generator: Experimental results and heat transfer evaluation | |
Sun et al. | Transient analysis and dynamic modeling of the steam generator water level for nuclear power plants | |
González-Gómez et al. | Assessment of evaporators using solar salt as heat transfer fluid | |
Su et al. | Frictional pressure drop correlation of steam-water two-phase flow in helically coiled tubes | |
Gómez et al. | Design and dynamic analysis of steam generators for concentrating solar power plants | |
Jeon et al. | Development of heat transfer model package for horizontal U-shaped heat exchanger submerged in pool of passive safety system | |
CN112417681B (zh) | 一种蒸汽发生器一二次侧对流换热系数分布的估计方法 | |
Charles et al. | Loop thermosyphon design for solar thermal desalination | |
Unal et al. | Pressure drop at rod-bundle spacers in the post-CHF dispersed flow regime | |
Hanson et al. | ECC performance in the semiscale geometry | |
Manabe et al. | Moisture separator reheater for NPP turbines |
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 |