CN113434988B - 一种核动力系统中热管的传热特性分析方法 - Google Patents

一种核动力系统中热管的传热特性分析方法 Download PDF

Info

Publication number
CN113434988B
CN113434988B CN202110655604.XA CN202110655604A CN113434988B CN 113434988 B CN113434988 B CN 113434988B CN 202110655604 A CN202110655604 A CN 202110655604A CN 113434988 B CN113434988 B CN 113434988B
Authority
CN
China
Prior art keywords
heat
pipe
temperature
heat transfer
calculation
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
Application number
CN202110655604.XA
Other languages
English (en)
Other versions
CN113434988A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202110655604.XA priority Critical patent/CN113434988B/zh
Publication of CN113434988A publication Critical patent/CN113434988A/zh
Application granted granted Critical
Publication of CN113434988B publication Critical patent/CN113434988B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computational Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

一种核动力系统中热管的传热特性分析方法,主要步骤如下:1、输入热管热源、热阱工况,工质、结构、几何参数,设定计算时间;2、热管划分控制体并初始化计算;3、计算管壁控制体温度随时间变化率,考虑热管轴向传热及边界条件;4、计算吸液芯内控制体温度随时间变化率,忽略液体流动引起的对流换热,视为纯导热;5、计算蒸汽腔内控制体温度随时间变化率,连续流动的蒸汽流动简化为热阻导热;6、用Gear算法求解方程组,完成所有时间点的计算;7、热管传热极限校核及等效热阻计算;8、对比分析热管关键参数对传热性能的影响。本发明对不同的热管进行瞬态计算及传热特性分析,为核动力系统中高性能的热管换热器设计提供建议与指导。

Description

一种核动力系统中热管的传热特性分析方法
技术领域
本发明涉及核动力技术领域换热技术,具体涉及一种海上核动力系统中热管的传热特性分析方法。
背景技术
核反应堆非能动的排出堆芯余热的安全装置被学者广泛研究,海上核动力与陆地核动力相比余热排出系统较为落后,且工作环境恶劣,安全性要求更高,但它工作的海洋环境是余热排出最理想、最方便的热阱,是天然优势。典型的核反应堆AP1000中非能动余热排出系统中的C型换热器,其原理为反应堆冷却液流过管内,冷却水流过壳侧换热,此类换热器应用于海上恶劣环境下的反应堆则存在安全隐患:(1)流体循环需要动力,其非能动性不足;(2)主回路中的冷却剂受核燃料辐射,如果换热管道或连接管道和挡板有损坏,放射性冷却剂会泄漏而产生危险;(3)核反应堆停堆或事故工况时,非能动余热排出系统投入运行,冷却剂突然流入换热器管内与管外冷水换热,易产生水锤而振动损坏换热器;(4)若有一根管失效,则会导致整个换热器无法正常运行。
发明内容
为了克服上述现有技术存在的问题,本发明的目的在于提供一种核动力系统中热管的传热特性分析方法,该方法针对不同结构特点及应用环境的热管依次进行各部分传热计算、传热极限校核、等效热阻计算、关键参数敏感性分析。首先计算初步设计参数时热管运行的瞬态温度分布,并校核其设计是否合理,进行传热特性分析。随后,改变热管边界条件或改变几何条件等关键参数进行计算,对比分析热管瞬态温度分布、传热极限、等效热阻,得出关键参数对热管传热特性的影响。本发明为核动力系统中高性能的换热技术提供理论建议与指导。
为了实现上述目的,本发明采取了以下技术方案予以实施:
一种核动力系统中热管的传热特性分析方法,步骤如下:
步骤1:输入热管热源及热阱条件,热源是给定热流密度的电加热丝加热或给定表面换热系数和热流体温度的流动加热;热阱是给定表面换热系数和冷流体温度的横掠冷却或套管冷却;输入热管的工质、结构、几何参数,设定总计算时间;
步骤2:将热管管壁和吸液芯沿厚度方向分别划分m层和n层,蒸汽腔划作单独一层,每层中蒸发段e个控制体,绝热段a个控制体,冷凝段c个控制体;对每个控制体进行初始化计算,获得初始时刻管壁、吸液芯和蒸汽腔的温度值;
步骤3:管壁温度变化率
Figure GDA0003782575170000021
计算:依据傅里叶定律在管壁建立二维导热方程,考虑轴向传热及周围介质的导热热阻;
Figure GDA0003782575170000022
式中:
T——温度,K
ρw——管壁材料密度,kg·m-3
Cp——管壁材料比热容,J·kg-1·K-1
t——时间,S
r——沿计算区域径向方向
z——沿计算区域径向方向
k——管壁材料导热率,W·m-1·K-1
导热方程式(1)的边界条件如下:
蒸发段:
Figure GDA0003782575170000031
电加丝加热:Qe=UI (3)
流动加热:Qe=hAe(Te,w-Te) (4)
绝热段:
Figure GDA0003782575170000032
冷凝段:
Figure GDA0003782575170000033
流体冷却:Qc=hAc(Tc,w-Tc) (7)
式中:Ae——蒸发段换热面积,m2
Ac——冷凝段换热面积,m2
U——电压,m2
I——电流,A
h——管外流体换热系数,W·m-2·K-1
Te,w——蒸发段壁面平均壁温,K
Tc,w——冷凝段壁面平均壁温,K
Te——蒸发段热流体温度,K
Tc——冷凝段壁温,K
不同边界条件时,式(4),式(7)中管外流体换热系数h计算:
Figure GDA0003782575170000041
直接横掠加热或冷却:
Figure GDA0003782575170000042
套管加热或冷却:
当Re<2300时:
Figure GDA0003782575170000043
当Re>2300时:Nu=0.023Re4/5Pr0.4 (11)
式中:
Re——雷诺数,
Figure GDA0003782575170000044
μ——流体动力粘度,pa·s
u——流体速度,m·s-1
Pr——普朗特数,
Figure GDA0003782575170000045
λ——管外流体导热系数,W·m-1·K-1
dp——套管外径,m
d0——管壁外径,m
cp——管外液体比热容,J·kg-1·K-1
dh——特征长度,m
ρ——管外液体密度,kg·m-3
公式(10)中套管内流动Nu数:
d<sub>p</sub>/d<sub>o</sub> Nu
0.00 -
0.05 17.46
0.10 11.56
0.25 7.73
0.50 5.74
1.00 4.86
步骤4:吸液芯内温度变化率
Figure GDA0003782575170000051
计算:此处忽略吸液芯结构中液体流动引起的对流换热,视为纯导热模型;因此吸液芯和管壁导热方程相同,差异仅在于有效导热系数的计算:
Figure GDA0003782575170000052
式中:
ρeff——吸液芯有效密度,kg·m-3
Cpeff——吸液芯有效比容,J·kg-1·K-1
keff——吸液芯的有效热导率,W·m-1·K-1
吸液芯的有效导热率keff计算:
Figure GDA0003782575170000053
式中:
kl——吸液芯工质导热系数,W·m-2·K-1
kw——吸液芯材料导热系数,W·m-2·K-1
ε——多孔介质结构孔隙率
步骤5:蒸汽腔内蒸汽温度变化率
Figure GDA0003782575170000061
计算,对于连续流动的蒸汽流蒸汽导热,简化为热阻网络:
Figure GDA0003782575170000062
式中:
ρv——蒸汽密度,kg·m-3
Cv——蒸汽比热容,KJ·kg-1·K-1
T+——气液界面温度,K
Rr——气液表面相变热阻
Rz——蒸汽腔轴向热阻,K·W-1
蒸汽轴向热阻Rz非常小以致忽略不计,相变热阻Rr计算公式如下:
Figure GDA0003782575170000063
式中:
Rn——气体常数,J·mol-1·K-1
Tv——蒸汽温度,K
hfg——汽化潜热,kJ·kg-1
p——蒸汽压力,pa
Av——蒸汽腔截面积,m2
步骤6:由于将热管各部分划分为控制体,则存在一系列的传热控制体,因此有一系列一阶线性方程,离散后,方程转化为带有初值的非线性微分方程组,统一体形式为:
Figure GDA0003782575170000071
此方程刚性大,调用Gear算法求解,若计算结果不收敛,输出不收敛,重新输入数据;若计算结果收敛,将其一个时间节点解出的值作为初值,重复步骤2-6过程,进行下一时间节点计算,直到设定的总计算时间,计算停止;
步骤7:以蒸汽腔温度稳定后的工作温度为基础,校核热管传热极限,包括沸腾极限、毛细极限和夹带极限;若传热极限大于输入功率,输出热管各部分瞬态温度计算结果、传热极限和有效热阻;若任意一个传热极限小于设计功率,则重复步骤1-7,重新输入数据进行计算;
沸腾极限:
Figure GDA0003782575170000072
夹带极限:
Figure GDA0003782575170000073
毛细极限:吸液芯设计应满足产生的毛细压力Δpc足以克服冷凝液体从冷凝段回流到蒸发段的压力降Δpl和重力对流体引起的压降Δpg,即Δpc≥Δpl+Δpg,从而推出毛细极限;此处考虑的结构为丝网结构、凹槽、丝网覆盖凹槽、丝网覆盖环形空隙、不同目数丝网复合;对于不同丝网目数复合结构,毛细压差为大目数丝网的毛细压力,流动压降为小目数丝网内的压降;丝网覆盖凹槽毛细压力为丝网的毛细压力,流动压降为凹槽内压降;丝网覆盖环形空隙毛细压力为丝网毛细压力,流动压降为环形空隙流动压降:
Figure GDA0003782575170000081
Figure GDA0003782575170000082
不同目数丝网结构:
Figure GDA0003782575170000083
凹槽覆盖丝网结构:
Figure GDA0003782575170000084
环形空隙覆盖丝网结构:
Figure GDA0003782575170000085
式中:
K——吸液芯丝网结构渗透率,
Figure GDA0003782575170000086
ε——多孔介质结构孔隙率
Aw——吸液芯结构截面积,m2
N——丝网目数,表示丝与丝之间的疏密程度,目数越大,丝与丝越密
σ——工质表面张力,N·m-1
re——丝径,mm
dv——蒸汽腔直径,mm
ρl——液态工质密度,kg·m-3
Figure GDA0003782575170000087
——热管倾角
l——热管长度,m
μl——工质动力粘度,pa·s
Figure GDA0003782575170000092
——毛细极限,W
x——凹槽个数
w——环形空隙宽度,mm
de——凹槽宽度,mm
le——热管蒸发段长度,m
ri/rv——热管内径和蒸汽腔直径之比
rb——汽泡的汽化核心半径,m
Δpc——最大毛细压力,Mpa
rhs——毛细吸液芯的水力半径,m
等效热阻计算:等效热阻是反映热管性能的重要参数之一,等效热阻越小,热管蒸发段和冷凝段温差越小,等温性能越好:
Figure GDA0003782575170000091
Te,w——蒸发段壁面平均壁温,K
Tc,w——冷凝段壁面平均壁温,K
R——等效热阻,K·W-1
Qe——传热功率,W
步骤8:对管径、倾角、各段长度重要参数进行敏感性分析,改变输入参数,进行步骤2-7计算,对比不同工况的传热极限及等效热阻大小,最终分析热管不同参数对传热特性的影响。
和现有技术相比较,本发明具备如下优点:
(1)对于热管壁面吸液芯传热分析,采用二维分析,考虑轴向传热。(2)考虑热管蒸发段及冷凝段外不同类型的边界条件。(3)热管传热能力采用传热极限进行校核,且给出一套不同类型复合型吸液芯的毛细极限计算理论。(4)满足传热极限后的热管,分析不同参数对其热阻的影响,进行敏感性分析,得出此应用环境下最优热管参数。
针对已存在的问题,本发明提出将热管式换热器应用于开式余热排出系统中实现非能动换热,热管是一种有效导热系数极高的非能动换热装置,能以较小的单位尺寸,在较低温降下传递大量的热量。热管换热器可以被认为是由两个耦合的热交换器组成,将其与开式余热排出系统耦合独特的优势在于:(1)热管换热是依靠管内的工质汽化潜热带走热量,热管内部工质的蒸发、凝结而非能动换热。事故工况时,即使所有电源都被切断的最坏情况下,热量也可以通过热管换热器有效的排出;(2)热管换热器将热源及热阱在物理上分离,热源与热阱均流过热管换热器的蒸发段外及冷凝段外,避免放射性冷却剂泄漏而交叉污染,高效地转移热能;(3)冷却剂冲刷热管换热器蒸发段外实现热量传递,而不是流入管内,避免了冷却剂突然流入管内与冷却水换热而产生水锤振动,引起管道破裂;(4)每根热管都是相互独立的,即使有热管发生故障,系统仍将安全运行,克服了换热器单点失效问题。
附图说明
图1为本发明的计算流程框图。
图2为热管壁面、吸液芯、蒸汽腔控制体划分及边界条件。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明。
如图1所示,本发明一种核动力系统中热管的传热特性分析方法,步骤如下:
步骤1:输入热管热源及热阱条件,热源是给定热流密度的电加热丝加热或给定表面换热系数和热流体温度的流动加热;热阱是给定表面换热系数和冷流体温度的横掠冷却或套管冷却;输入热管的工质、结构、几何参数,设定总计算时间;
步骤2:如图2所示,将热管管壁和吸液芯沿厚度方向分别划分m层和n层,蒸汽腔划作单独一层,每层中蒸发段e个控制体,绝热段a个控制体,冷凝段c个控制体;对每个控制体进行初始化计算,获得初始时刻管壁、吸液芯和蒸汽腔的温度值;
步骤3:管壁温度变化率
Figure GDA0003782575170000111
计算:依据傅里叶定律在管壁建立二维导热方程,考虑轴向传热及周围介质的导热热阻;
Figure GDA0003782575170000112
式中:
T——温度,K
ρw——管壁材料密度,kg·m-3
Cp——管壁材料比热容,J·kg-1·K-1
t——时间,S
r——沿计算区域径向方向
z——径向区域轴向方向
k——管壁材料导热率,W·m-1·K-1
导热方程式(1)的边界条件如下:
蒸发段:
Figure GDA0003782575170000121
电加丝加热:Qe=UI (3)
流动加热:Qe=hAe(Te,w-Te) (4)
绝热段:
Figure GDA0003782575170000122
冷凝段:
Figure GDA0003782575170000123
流体冷却:Qc=hAc(Tc,w-Tc) (7)
式中:Ae——蒸发段换热面积,m2
Ac——冷凝段换热面积,m2
U——电压,m2
I——电流,A
h——管外流体换热系数,W·m-2·K-1
Te,w——蒸发段壁面平均壁温,K
Tc,w——冷凝段壁面平均壁温,K
Te——蒸发段热流体温度,K
Tc——冷凝段壁温,K
不同边界条件时,式(4),式(7)中管外流体换热系数h计算:
Figure GDA0003782575170000124
直接横掠加热或冷却:
Figure GDA0003782575170000131
套管加热或冷却:
当Re<2300时:
Figure GDA0003782575170000132
当Re>2300时:Nu=0.023Re4/5Pr0.4 (11)
式中:
Re——雷诺数,
Figure GDA0003782575170000133
μ——流体动力粘度,pa·s
u——流体速度,m·s-1
Pr——普朗特数,
Figure GDA0003782575170000134
λ——管外流体导热系数,W·m-1·K-1
dp——套管外径,m
d0——管壁外径,m
cp——管外液体比热容,J·kg-1·K-1
dh——特征长度,m
ρ——管外液体密度,kg·m-3
公式(10)中套管内流动Nu数:
Figure GDA0003782575170000135
Figure GDA0003782575170000141
步骤4:吸液芯内温度变化率
Figure GDA0003782575170000142
计算:此处忽略吸液芯结构中液体流动引起的对流换热,视为纯导热模型;因此吸液芯和管壁导热方程相同,差异仅在于有效导热系数的计算:
Figure GDA0003782575170000143
式中:
ρeff——吸液芯有效密度,kg·m-3
Cpeff——吸液芯有效比容,J·kg-1·K-1
keff——吸液芯的有效热导率,W·m-1·K-1
吸液芯的有效导热率keff计算:
Figure GDA0003782575170000144
式中:
kl——吸液芯工质导热系数,W·m-2·K-1
kw——吸液芯材料导热系数,W·m-2·K-1
ε——多孔介质结构孔隙率
步骤5:蒸汽腔内蒸汽温度变化率
Figure GDA0003782575170000145
计算,对于连续流动的蒸汽流蒸汽导热,简化为热阻网络:
Figure GDA0003782575170000146
式中:
ρv——蒸汽密度,kg·m-3
Cv——蒸汽比热容,KJ·kg-1·K-1
T+——气液界面温度,K
Rr——气液表面相变热阻
Rz——蒸汽腔轴向热阻,K·W-1
蒸汽轴向热阻Rz非常小以致忽略不计,相变热阻Rr计算公式如下:
Figure GDA0003782575170000151
式中:
Rn——气体常数,J·mol-1·K-1
Tv——蒸汽温度,K
hfg——汽化潜热,kJ·kg-1
p——蒸汽压力,pa
Av——蒸汽腔截面积,m2
步骤6:由于将热管各部分划分为控制体,则存在一系列的传热控制体,因此有一系列一阶线性方程,离散后,方程转化为带有初值的非线性微分方程组,统一体形式为:
Figure GDA0003782575170000152
此方程刚性大,调用Gear算法求解,若计算结果不收敛,输出不收敛,重新输入数据;若计算结果收敛,将其一个时间节点解出的值作为初值,重复步骤2-6过程,进行下一时间节点计算,直到设定的总计算时间,计算停止;
步骤7:以蒸汽腔温度稳定后的工作温度为基础,校核热管传热极限,包括沸腾极限、毛细极限和夹带极限;若传热极限大于输入功率,输出热管各部分瞬态温度计算结果、传热极限和有效热阻;若任意一个传热极限小于设计功率,则重复步骤1-7,重新输入数据进行计算;
沸腾极限:
Figure GDA0003782575170000161
夹带极限:
Figure GDA0003782575170000162
毛细极限:吸液芯设计应满足产生的毛细压力Δpc足以克服冷凝液体从冷凝段回流到蒸发段的压力降Δpl和重力对流体引起的压降Δpg,即Δpc≥Δpl+Δpg,从而推出毛细极限;此处考虑的结构为丝网结构、凹槽、丝网覆盖凹槽、丝网覆盖环形空隙、不同目数丝网复合;对于不同丝网目数复合结构,毛细压差为大目数丝网的毛细压力,流动压降为小目数丝网内的压降;丝网覆盖凹槽毛细压力为丝网的毛细压力,流动压降为凹槽内压降;丝网覆盖环形空隙毛细压力为丝网毛细压力,流动压降为环形空隙流动压降:
Figure GDA0003782575170000163
Figure GDA0003782575170000164
不同目数丝网结构:
Figure GDA0003782575170000171
凹槽覆盖丝网结构:
Figure GDA0003782575170000172
环形空隙覆盖丝网结构:
Figure GDA0003782575170000173
式中:
K——吸液芯丝网结构渗透率,
Figure GDA0003782575170000174
ε——多孔介质结构孔隙率
Aw——吸液芯结构截面积,m2
N——丝网目数,表示丝与丝之间的疏密程度,目数越大,丝与丝越密
σ——工质表面张力,N·m-1
re——丝径,mm
dv——蒸汽腔直径,mm
ρl——液态工质密度,kg·m-3
Figure GDA0003782575170000175
——热管倾角
l——热管长度,m
μl——工质动力粘度,pa·s
Figure GDA0003782575170000176
——毛细极限,W
x——凹槽个数
w——环形空隙宽度,mm
de——凹槽宽度,mm
le——热管蒸发段长度,m
ri/rv——热管内径和蒸汽腔直径之比
rb——汽泡的汽化核心半径,m
Δpc——最大毛细压力,Mpa
rhs——毛细吸液芯的水力半径,m
等效热阻计算:等效热阻是反映热管性能的重要参数之一,等效热阻越小,热管蒸发段和冷凝段温差越小,等温性能越好:
Figure GDA0003782575170000181
Te,w——蒸发段壁面平均壁温,K
Tc,w——冷凝段壁面平均壁温,K
R——等效热阻,K·W-1
Qe——传热功率,W
步骤8:对管径、倾角、各段长度重要参数进行敏感性分析,改变输入参数,进行步骤2-7计算,对比不同工况的传热极限及等效热阻大小,最终分析热管不同参数对传热特性的影响。

Claims (1)

1.一种核动力系统中热管的传热特性分析方法,其特征在于:步骤如下:
步骤1:输入热管热源及热阱条件,热源是给定热流密度的电加热丝加热或给定表面换热系数和热流体温度的流动加热;热阱是给定表面换热系数和冷流体温度的横掠冷却或套管冷却;输入热管的工质、结构、几何参数,设定总计算时间;
步骤2:将热管管壁和吸液芯沿厚度方向分别划分m层和n层,蒸汽腔划作单独一层,每层中蒸发段e个控制体,绝热段a个控制体,冷凝段c个控制体;对每个控制体进行初始化计算,获得初始时刻管壁、吸液芯和蒸汽腔的温度值;
步骤3:管壁温度变化率
Figure FDA0003782575160000011
计算:依据傅里叶定律在管壁建立二维导热方程,考虑轴向传热及周围介质的导热热阻;
Figure FDA0003782575160000012
式中:
T——温度,K
ρw——管壁材料密度,kg·m-3
Cp——管壁材料比热容,J·kg-1·K-1
t——时间,S
r——沿计算区域径向方向
z——沿径向区域轴向方向
k——管壁材料导热率,W·m-1·K-1
导热方程式(1)的边界条件如下:
蒸发段:
Figure FDA0003782575160000021
电加丝加热:Qe=UI (3)
流动加热:Qe=hAe(Te,w-Te) (4)
绝热段:
Figure FDA0003782575160000022
冷凝段:
Figure FDA0003782575160000023
流体冷却:Qc=hAc(Tc,w-Tc) (7)
式中:Ae——蒸发段换热面积,m2
Ac——冷凝段换热面积,m2
U——电压,m2
I——电流,A
h——管外流体换热系数,W·m-2·K-1
Te,w——蒸发段壁面平均壁温,K
Tc,w——冷凝段壁面平均壁温,K
Te——蒸发段热流体温度,K
Tc——冷凝段壁温,K
不同边界条件时,式(4),式(7)中管外流体换热系数h计算:
Figure FDA0003782575160000024
直接横掠加热或冷却:
Figure FDA0003782575160000025
套管加热或冷却:
当Re<2300时:
Figure FDA0003782575160000031
当Re>2300时:Nu=0.023Re4/5Pr0.4 (11)
式中:
Re——雷诺数,
Figure FDA0003782575160000032
μ——流体动力粘度,pa·s
u——流体速度,m·s-1
Pr——普朗特数,
Figure FDA0003782575160000033
λ——管外流体导热系数,W·m-1·K-1
dp——套管外径,m
d0——管壁外径,m
cp——管外液体比热容,J·kg-1·K-1
dh——特征长度,m
ρ——管外液体密度,kg·m-3
公式(10)中套管内流动Nu数:
d<sub>p</sub>/d<sub>o</sub> Nu 0.00 - 0.05 17.46 0.10 11.56 0.25 7.73 0.50 5.74 1.00 4.86
步骤4:吸液芯内温度变化率
Figure FDA0003782575160000041
计算:此处忽略吸液芯结构中液体流动引起的对流换热,视为纯导热模型;因此吸液芯和管壁导热方程相同,差异仅在于有效导热系数的计算:
Figure FDA0003782575160000042
式中:
ρeff——吸液芯有效密度,kg·m-3
Cpeff——吸液芯有效比容,J·kg-1·K-1
keff——吸液芯的有效热导率,W·m-1·K-1
吸液芯的有效导热率keff计算:
Figure FDA0003782575160000043
式中:
kl——吸液芯工质导热系数,W·m-2·K-1
kw——吸液芯材料导热系数,W·m-2·K-1
ε——多孔介质结构孔隙率
步骤5:蒸汽腔内蒸汽温度变化率
Figure FDA0003782575160000044
计算,对于连续流动的蒸汽流蒸汽导热,简化为热阻网络:
Figure FDA0003782575160000045
式中:
ρv——蒸汽密度,kg·m-3
Cv——蒸汽比热容,KJ·kg-1·K-1
T+——气液界面温度,K
Rr——气液表面相变热阻
Rz——蒸汽腔轴向热阻,K·W-1
蒸汽轴向热阻Rz非常小以致忽略不计,相变热阻Rr计算公式如下:
Figure FDA0003782575160000051
式中:
Rn——气体常数,J·mol-1·K-1
Tv——蒸汽温度,K
hfg——汽化潜热,kJ·kg-1
p——蒸汽压力,pa
Av——蒸汽腔截面积,m2
步骤6:由于将热管各部分划分为控制体,则存在一系列的传热控制体,因此有一系列一阶线性方程,离散后,方程转化为带有初值的非线性微分方程组,统一体形式为:
Figure FDA0003782575160000052
此方程刚性大,调用Gear算法求解,若计算结果不收敛,输出不收敛,重新输入数据;若计算结果收敛,将其一个时间节点解出的值作为初值,重复步骤2-6过程,进行下一时间节点计算,直到设定的总计算时间,计算停止;
步骤7:以蒸汽腔温度稳定后的工作温度为基础,校核热管传热极限,包括沸腾极限、毛细极限和夹带极限;若传热极限大于输入功率,输出热管各部分瞬态温度计算结果、传热极限和有效热阻;若任意一个传热极限小于设计功率,则重复步骤1-7,重新输入数据进行计算;
沸腾极限:
Figure FDA0003782575160000061
夹带极限:
Figure FDA0003782575160000062
毛细极限:吸液芯设计应满足产生的毛细压力Δpc足以克服冷凝液体从冷凝段回流到蒸发段的压力降Δpl和重力对流体引起的压降Δpg,即Δpc≥Δpl+Δpg,从而推出毛细极限;此处考虑的结构为丝网结构、凹槽、丝网覆盖凹槽、丝网覆盖环形空隙、不同目数丝网复合;对于不同丝网目数复合结构,毛细压差为大目数丝网的毛细压力,流动压降为小目数丝网内的压降;丝网覆盖凹槽毛细压力为丝网的毛细压力,流动压降为凹槽内压降;丝网覆盖环形空隙毛细压力为丝网毛细压力,流动压降为环形空隙流动压降:
Figure FDA0003782575160000063
Figure FDA0003782575160000064
不同目数丝网结构:
Figure FDA0003782575160000065
凹槽覆盖丝网结构:
Figure FDA0003782575160000066
环形空隙覆盖丝网结构:
Figure FDA0003782575160000071
式中:
K——吸液芯丝网结构渗透率,
Figure FDA0003782575160000072
ε——多孔介质结构孔隙率
Aw——吸液芯结构截面积,m2
N——丝网目数,表示丝与丝之间的疏密程度,目数越大,丝与丝越密
σ——工质表面张力,N·m-1
re——丝径,mm
dv——蒸汽腔直径,mm
ρl——液态工质密度,kg·m-3
Figure FDA0003782575160000073
——热管倾角
l——热管长度,m
μl——工质动力粘度,pa·s
Figure FDA0003782575160000074
——毛细极限,W
x——凹槽个数
w——环形空隙宽度,mm
de——凹槽宽度,mm
le——热管蒸发段长度,m
ri/rv——热管内径和蒸汽腔直径之比
rb——汽泡的汽化核心半径,m
Δpc——最大毛细压力,Mpa
rhs——毛细吸液芯的水力半径,m
等效热阻计算:等效热阻是反映热管性能的重要参数之一,等效热阻越小,热管蒸发段和冷凝段温差越小,等温性能越好:
Figure FDA0003782575160000081
Te,w——蒸发段壁面平均壁温,K
Tc,w——冷凝段壁面平均壁温,K
R——等效热阻,K·W-1
Qe——传热功率,W
步骤8:对管径、倾角、各段长度重要参数进行敏感性分析,改变输入参数,进行步骤2-7计算,对比不同工况的传热极限及等效热阻大小,最终分析热管不同参数对传热特性的影响。
CN202110655604.XA 2021-06-11 2021-06-11 一种核动力系统中热管的传热特性分析方法 Active CN113434988B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110655604.XA CN113434988B (zh) 2021-06-11 2021-06-11 一种核动力系统中热管的传热特性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110655604.XA CN113434988B (zh) 2021-06-11 2021-06-11 一种核动力系统中热管的传热特性分析方法

Publications (2)

Publication Number Publication Date
CN113434988A CN113434988A (zh) 2021-09-24
CN113434988B true CN113434988B (zh) 2022-10-28

Family

ID=77755810

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110655604.XA Active CN113434988B (zh) 2021-06-11 2021-06-11 一种核动力系统中热管的传热特性分析方法

Country Status (1)

Country Link
CN (1) CN113434988B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113948229B (zh) * 2021-09-30 2023-08-01 西安交通大学 一种海洋核动力热管式余热排出系统稳定运行条件分析方法
CN114154438A (zh) * 2021-12-07 2022-03-08 西安交通大学 一种碱金属热管冷态启动三阶段计算方法
CN114239435B (zh) * 2021-12-20 2023-05-16 西安交通大学 一种三维平板水热管数值计算方法
CN114893715B (zh) * 2022-04-02 2023-11-21 安徽宇航派蒙健康科技股份有限公司 加热控制方法及其装置、系统、计算机设备和存储介质
CN115659908B (zh) * 2022-12-22 2023-05-09 西安交通大学 一种印刷电路板换热器的三单元非平衡多孔介质方法
CN116045717B (zh) * 2023-02-13 2023-07-04 中国核动力研究设计院 热质传递装置、换热系数的计算方法、装置、设备及介质
CN116502555B (zh) * 2023-04-28 2024-04-09 西安交通大学 一种碱金属热管吸液芯内部流速压力计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991002096A1 (de) * 1989-07-28 1991-02-21 Paul Braisch Verfahren zur werkstoffabhängigen steuerung von wärmebehandlungsprozessen von metallen und vorrichtung zur durchführung des verfahrens
JP2006010636A (ja) * 2004-06-29 2006-01-12 Canon Inc 熱流体解析方法およびその方法を実行するためのプログラムおよびそのプログラムを記憶した記憶媒体および熱流体解析装置
CN107180657A (zh) * 2017-06-12 2017-09-19 中广核工程有限公司 一种核电站多样化热阱系统传热性能试验系统和方法
CN108599137A (zh) * 2017-12-28 2018-09-28 东南大学 一种考虑区域热网暂态传热特性的多能流系统优化运行方法
CN110345787A (zh) * 2019-07-24 2019-10-18 西安交通大学 一种用于一体化高温碱金属热管的设计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080099569A1 (en) * 2006-10-31 2008-05-01 Husky Injection Molding Systems Ltd. Thermal Analysis of Apparatus having Multiple Thermal Control Zones
CN102494325B (zh) * 2011-12-19 2014-07-09 上海望特能源科技有限公司 电站锅炉高温管系炉内动态壁温监测的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1991002096A1 (de) * 1989-07-28 1991-02-21 Paul Braisch Verfahren zur werkstoffabhängigen steuerung von wärmebehandlungsprozessen von metallen und vorrichtung zur durchführung des verfahrens
JP2006010636A (ja) * 2004-06-29 2006-01-12 Canon Inc 熱流体解析方法およびその方法を実行するためのプログラムおよびそのプログラムを記憶した記憶媒体および熱流体解析装置
CN107180657A (zh) * 2017-06-12 2017-09-19 中广核工程有限公司 一种核电站多样化热阱系统传热性能试验系统和方法
CN108599137A (zh) * 2017-12-28 2018-09-28 东南大学 一种考虑区域热网暂态传热特性的多能流系统优化运行方法
CN110345787A (zh) * 2019-07-24 2019-10-18 西安交通大学 一种用于一体化高温碱金属热管的设计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Conceptual design and analysis of heat pipe cooled silo cooling system for the transportable fluoride-salt-cooled high-temperature reactor;Chenglong Wang, Limin Liu, Minghao Liu, Dalin Zhang;《Annals of Nuclear Energy》;20170606;第458-468页 *
基于SIMPLEC算法的高温热管启动特性数值模拟;葛攀和等;《原子能科学技术》;20170907(第11期);第1974-1981页 *
热管技术及性能分析程序研究;陈其昌;《原子能科学技术》;20200528(第07期);第1176-1184页 *
高温钾热管稳态运行传热特性研究;田智星,刘逍,王成龙,苏光辉,田文喜,秋穗正;《原子能科学技术》;20201031;第54卷(第10期);第1771-1778页 *

Also Published As

Publication number Publication date
CN113434988A (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
CN113434988B (zh) 一种核动力系统中热管的传热特性分析方法
Salem et al. Experimental investigation on the thermal performance of a double pipe heat exchanger with segmental perforated baffles
Elkholy et al. Experimental investigation of geyser boiling in a small diameter two-phase loop thermosyphon
Wang et al. A computationally derived heat transfer correlation for in-tube cooling turbulent supercritical CO2
Chen et al. Experimental investigation on performances of trisection helical baffled heat exchangers for oil/water–water heat transfer
Gou et al. An assessment of heat transfer models of water flow in helically coiled tubes based on selected experimental datasets
Liu et al. Experimental investigation on shell-side performance of a novel shell and tube oil cooler with twisted oval tubes
Sadhu et al. Steady-state analysis of a high-temperature natural circulation loop based on water-cooled supercritical CO2
Ma et al. Experimental study on sodium Screen-Wick heat pipe capillary limit
Zhao et al. Structural effect of internal composite wick on the anti-gravity heat transfer performance of a concentric annular high-temperature heat pipe
Khanmohammadi et al. Numerical investigation of heat transfer and fluid flow characteristics inside tube with internally star fins
CN114154438A (zh) 一种碱金属热管冷态启动三阶段计算方法
CN110195990B (zh) 一种蒸汽温度控制阀门蒸汽发生器
Smusz et al. Coil heat exchanger with the nanofluid filled buffer layer
CN113948229B (zh) 一种海洋核动力热管式余热排出系统稳定运行条件分析方法
Birajdar et al. Effect of adiabatic length on the performance of closed‐loop thermosyphon system
Zhang et al. Optimal research of annular baffle for heat‐discharging performance of single thermal storage tank with molten salt
Al-Kayiem et al. Ribbed double pipe heat exchanger: analytical analysis
Ma et al. A transient and steady-state network model forannular-wick heat pipes in continuum flow pattern
Kibushi et al. Heat transfer and pressure drop correlations for a gas flow in unidirectional porous copper tubes
Birajdar et al. Experimental investigations of pump‐driven closed‐loop thermosyphon system
Liu et al. Numerical Model Development for the Transient Heat and Mass Transfer Performance of a Water-Copper Wicked Heat Pipe
Pabón et al. Performance of a small diameter two-phase closed thermosyphon in geyser boiling condition
Kızılırmak et al. Three-dimensional turbulent flow, heat transfer and second-law analysis in a circular duct with single baffle
Xie et al. Numerical analysis of flow instability in the water wall of a supercritical CFB boiler with annular furnace

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