CN116502555A - 一种碱金属热管吸液芯内部流速压力计算方法 - Google Patents
一种碱金属热管吸液芯内部流速压力计算方法 Download PDFInfo
- Publication number
- CN116502555A CN116502555A CN202310476490.1A CN202310476490A CN116502555A CN 116502555 A CN116502555 A CN 116502555A CN 202310476490 A CN202310476490 A CN 202310476490A CN 116502555 A CN116502555 A CN 116502555A
- Authority
- CN
- China
- Prior art keywords
- liquid
- gas
- formula
- suction core
- working medium
- 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
Links
- 239000007788 liquid Substances 0.000 title claims abstract description 175
- 229910052783 alkali metal Inorganic materials 0.000 title claims abstract description 20
- 150000001340 alkali metals Chemical class 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 230000004907 flux Effects 0.000 claims abstract description 23
- 239000000463 material Substances 0.000 claims abstract description 16
- 238000010521 absorption reaction Methods 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims abstract description 5
- 238000009833 condensation Methods 0.000 claims description 13
- 230000005494 condensation Effects 0.000 claims description 13
- 238000001704 evaporation Methods 0.000 claims description 13
- 230000000704 physical effect Effects 0.000 claims description 12
- 230000008020 evaporation Effects 0.000 claims description 10
- 230000007704 transition Effects 0.000 claims description 6
- 238000009413 insulation Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000009834 vaporization Methods 0.000 claims description 3
- 230000008016 vaporization Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 2
- 230000001052 transient effect Effects 0.000 abstract description 2
- 230000007613 environmental effect Effects 0.000 abstract 1
- 239000007791 liquid phase Substances 0.000 abstract 1
- 238000002844 melting Methods 0.000 abstract 1
- 230000008018 melting Effects 0.000 abstract 1
- 239000012071 phase Substances 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005555 metalworking Methods 0.000 description 1
- 239000002689 soil Substances 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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/12—Timing analysis or timing 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
一种碱金属热管吸液芯内部流速压力计算方法,1、确定碱金属热管几何结构,各部件材料及外部环境条件,建立热管模型;2、确定吸液芯中工质熔化界面位置,确定蒸气区连续锋面位置;3、计算气液界面温度分布;4、计算气液界面质量通量以及相变导致的能量通量,更新气液界面边界条件;5、构建吸液芯内液相工质流动方程;6、计算管壁区域和吸液芯区域温度分布;7、迭代求解吸液芯压力和流速分布、管壁和吸液芯温度分布;8、计算完成,输出计算结果。本发明考虑高温热管运行过程中吸液芯内工质流动特性,计算吸液芯内液态工质压力和流速分布,进而获得碱金属热管内部工质的瞬态传热流动特性。为碱金属高温热管设计和工程应用提供建议和指导。
Description
技术领域
本发明涉及相变换热设备技术领域,具体涉及一种碱金属热管吸液芯内部流速压力计算方法。
背景技术
热管是一种高效非能动热量传输装置,具有很高的热导率和良好的等温性,热管在稳态工况下整体温差很小。高温热管的主要构造包括装有碱金属工质的真空管以及管内壁的复合丝网吸液芯。在轴向上,热管一般被划分为蒸发段、绝热段和冷凝段三个部分,热量通过蒸发段输入到热管中,通过冷凝段输出,绝热段仅起连接和传质作用,各个区域长度可灵活布置。热管壁为整个热管系统外边界,吸液芯产生毛细驱动力,驱动液态工质从冷凝段回流到蒸发段,蒸气区内为气态工质,其中热量以蒸气潜热的形式从蒸发段输送到冷凝段。热管已广泛应用于航空航天、高原冻土等领域,在先进核能系统中具有独特的优势和广阔应用前景。碱金属高温热管内部吸液芯和蒸气区存在传热传质耦合特性,吸液芯内液体流动情况较为复杂,导致针对热管的数值模拟计算难以同时保证计算速度和精度。
发明内容
为了克服上述现有技术存在的问题,本发明的目的在于提供一种碱金属热管吸液芯内部流速压力计算方法,该方法考虑高温热管运行过程中吸液芯内工质流动特性,针对不同结构及工作状况的碱金属高温热管进行数值模拟,计算吸液芯内液态工质压力和流速分布,计算吸液芯和管壁温度分布,进而获得碱金属热管内部工质的瞬态传热流动特性。为碱金属高温热管设计和工程应用提供建议和指导。
为了实现上述目的,本发明采取以下技术方案:
一种碱金属热管吸液芯内部流速压力计算方法,步骤如下:
步骤1:确定碱金属热管几何结构,各部件材料及外部环境条件,建立热管模型:划分蒸发段、绝热段和冷凝段,根据热管内部状态划分管壁、吸液芯、蒸气区三部分计算域;蒸发段给定外壁温度或给定热流密度、绝热段给定零温度梯度、冷凝段给定换热系数和参考温度;设置热管整体温度初始条件,设置计算时间步长及计算总步数;
步骤2:根据吸液芯内部温度分布,确定吸液芯内工质液固两相分界面位置;通过气液界面转变温度确定蒸气区中连续流动气体与分子流态气体分界面位置;通过式(1)所示转变温度计算关系式,划分蒸气区流态;
式(1)中:
Ttrans——转变温度/K;
M——相对分子质量/g·mol-1;
Ru——通用气体常数/J·mol-1·K-1;
μ——动力粘度系数/Pa·s;
ρ——气相工质密度/kg·m-3;
D——蒸气区直径/m;
步骤3:计算气液界面温度分布:分子流态部分气液界面采用绝热边界,连续流动气体部分采用能量平衡边界;
通过式(2)所示的绝热边界计算气液界面温度:
式(2)中
T——吸液芯温度分布/K;
f——气液界面位置;
通过式(3)所示边界计算气液界面温度:
式(3)中:
kwi——吸液芯等效导热系数/W·m-2·K-1;
mi——气液界面质量通量/kg·m-2·s-1;
hfg——汽化潜热/J·kg-1;
步骤4:计算气液界面质量通量以及相变导致的能量通量,更新气液界面边界条件:
通过式(4)所示的拟合函数关系计算气液界面饱和压力:
式(4)中:
Pf——气液界面饱和压力/Pa;
A、B、C——拟合关系式系数;
通过式(5)、式(6)所示气体动力学关系式计算气液界面质量通量和能量通量:
式(5)中:
ε——吸液芯孔隙率;
——蒸发冷凝调节系数;
M——蒸气区蒸气的摩尔质量/(kg·mol-1);
Tli——气液界面上的节点液态工质温度/K;
Pli——气液界面上的节点液态工质饱和压力/Pa;
Tv——蒸气区蒸气的温度/K;
Pv——蒸气区蒸气的饱和压力/Pa;
qi——气液界面能量通量/J·m-2·s-1;
步骤5:构建吸液芯内液态工质流动物理现象及气液界面边界条件构建吸液芯质量、动量守恒方程,得到吸液芯内液态工质压力分布和流速分布描述方程组:
吸液芯内液态工质流动连续方程如式(7)所示:
式(7)中:
ρl——吸液芯内液态工质密度/kg·m-3;
t——时间/s;
——液态工质表观流速矢量/m·s-1;
A——气液界面上面积/m2;
吸液芯动量方程如式(8)所示:
式(8)中:
μl——液态工质动力粘度系数/Pa·s;
K——吸液芯结构渗透率/m2;
Pl——吸液芯内部相对压力分布/Pa;
步骤6:计算管壁区域和吸液芯区域温度分布,对于管壁使用纯导热控制方程,对于吸液芯采用流动换热控制方程;
吸液芯部分流动换热控制方程如式(9)所示:
式(13)中:
(ρCP)eff——等效质量热容/kg·m-3
T——热管温度分布/K
keff——等效导热系数/W·m-2·K-1;
对于管壁区域,流速设为零,各项物性均采用材料自身关系式计算;而对于吸液芯区域,采用Chi模型计算各项物性,物性计算关系式如式(10)、式(11)所示;
(ρc)eff=ερlcl+(1-ε)ρwscws (10)
keff——吸液芯的等效热导率/(W·m-1·K-1)
(ρc)eff——吸液芯内工质与多孔介质混合等效热容(J·kg-1·K-1);
cl——吸液芯内液态工质的热容/(J·m-3·K-1);
ρws——吸液芯材料的密度/(kg·m-3);
cws——吸液芯材料的热容/(J·m-3·K-1);
kl——吸液芯内液态工质的热导率/(W·m-1·K-1);
kws——吸液芯材料的热导率/(W·m-1·K-1);
步骤7:计算吸液芯内液态工质压力分布和流速分布,计算管壁及吸液芯温度分布:通过迭代求解吸液芯内流动方程组和管壁吸液芯温度方程组,得到吸液芯内压力和流速分布,以及管壁吸液芯温度分布;采用超松弛迭代法进行线性方程组的求解;
对于式(12)所示的方程组:
改写为式(13)的形式:
anm——控制方程左侧系数矩阵第n行第m列元素;
bn——控制方程右侧向量第n行元素;
x(k)——解向量第k次迭代计算结果;
对得出结果进行修改,如式(14)所示:
式(14)中:
——迭代残向量;
对式(13)的结果进行变换得;
式(15)中:
ω——松弛因子;
迭代时,预设迭代精度ε,判断是否满足如果满足,则停止迭代,;如果不满足,则重复此步骤,将x(k)代入式(15)计算得到x(k+1),直到n次迭代后满足停止迭代,得到吸液芯压力分布及温度分布的近似值;
步骤8:当前时间步计算完成,输出计算结果;通过计算结果更新物性,更新气液界面边界条件,进行下一时间步计算。
和现有技术相比较,本发明具备如下优点:
(1)考虑热管吸液芯流动特性,可以准确计算吸液芯内液态工质压力流速分布。
(2)构造的迭代方程形式简洁且收敛性较好。
(3)由理论公式计算得到气液界面质量通量和能量通量,并以边界条件形式作用于吸液芯流动传热计算。
(4)可适用于多种几何结构和工质的碱金属热管,具有广泛的适用性。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明。
如图1所示,本发明一种碱金属热管吸液芯内部流速压力计算方法,步骤如下:
步骤1:确定碱金属热管几何结构,各部件材料及外部环境条件,建立热管模型:划分蒸发段、绝热段和冷凝段,根据热管内部状态划分管壁、吸液芯、蒸气区三部分计算域;蒸发段给定外壁温度或给定热流密度、绝热段给定零温度梯度、冷凝段给定换热系数和参考温度;设置热管整体温度初始条件,设置计算时间步长及计算总步数;
步骤2:根据吸液芯内部温度分布,确定吸液芯内工质液固两相分界面位置;通过气液界面转变温度确定蒸气区中连续流动气体与分子流态气体分界面位置;通过式(1)所示转变温度计算关系式,划分蒸气区流态;
按照无量纲数克努森数Kn=0.01划分自由分子流态和连续流态
式(1)中:
Ttrans——转变温度/K;
M——相对分子质量/g·mol-1;
Ru——通用气体常数/J·mol-1·K-1;
μ——动力粘度系数/Pa·s;
ρ——气相工质密度/kg·m-3;
D——蒸气区直径/m;
步骤3:计算气液界面温度分布:分子流态部分气液界面采用绝热边界,连续流动气体部分采用能量平衡边界;
通过式(2)所示的绝热边界计算气液界面温度:
式(2)中
T——吸液芯温度分布/K;
f——气液界面位置;
通过式(3)所示边界计算气液界面温度:
式(3)中:
kwi——吸液芯等效导热系数/W·m-2·K-1;
mi——气液界面质量通量/kg·m-2·s-1;
hfg——汽化潜热/J·kg-1;
通过吸液芯内部温度分布和气液界面能量通量计算得到气液界面温度分布
步骤4:计算气液界面质量通量以及相变导致的能量通量,更新气液界面边界条件:
通过式(4)所示的拟合函数关系计算气液界面饱和压力:
式(4)中:
Pf——气液界面饱和压力/Pa;
A、B、C——拟合关系式系数;
通过式(5)、式(6)所示气体动力学关系式计算气液界面质量通量和能量通量:
式(5)中:
ε——吸液芯孔隙率;
——蒸发冷凝调节系数;
M——蒸气区蒸气的摩尔质量/(kg·mol-1);
Tli——气液界面上的节点液态工质温度/K;
Pli——气液界面上的节点液态工质饱和压力/Pa;
Tv——蒸气区蒸气的温度/K;
Pv——蒸气区蒸气的饱和压力/Pa;
qi——气液界面能量通量/J·m-2·s-1;
步骤5:构建吸液芯内液态工质流动物理现象及气液界面边界条件构建吸液芯质量、动量守恒方程,得到吸液芯内液态工质压力分布和流速分布描述方程组:
吸液芯内液态工质流动连续方程如式(7)所示:
式(7)中:
ρl——吸液芯内液态工质密度/kg·m-3;
t——时间/s;
——液态工质表观流速矢量/m·s-1;
A——气液界面上面积/m2;
吸液芯动量方程如式(8)所示:
式(8)中:
μl——液态工质动力粘度系数/Pa·s;
K——吸液芯结构渗透率/m2;
Pl——吸液芯内部相对压力分布/Pa
步骤6:计算管壁区域和吸液芯区域温度分布,对于管壁使用纯导热控制方程,对于吸液芯采用流动换热控制方程;
吸液芯部分流动换热控制方程如式(9)所示:
式(13)中:
(ρCP)eff——等效质量热容/kg·m-3
T——热管温度分布/K
keff——等效导热系数/W·m-2·K-1;
对于管壁区域,流速设为零,各项物性均采用材料自身关系式计算;而对于吸液芯区域,采用Chi模型计算各项物性,物性计算关系式如式(11)、式(11)所示;
(ρc)eff=ερlcl+(1-ε)ρwscws (10)
keff——吸液芯的等效热导率/(W·m-1·K-1)
(ρc)eff——吸液芯内工质与多孔介质混合等效热容(J·kg-1·K-1);
cl——吸液芯内液态工质的热容/(J·m-3·K-1);
ρws——吸液芯材料的密度/(kg·m-3);
cws——吸液芯材料的热容/(J·m-3·K-1);
kl——吸液芯内液态工质的热导率/(W·m-1·K-1);
kws——吸液芯材料的热导率/(W·m-1·K-1);
步骤7:计算吸液芯内液态工质压力分布和流速分布,计算管壁及吸液芯温度分布:通过迭代求解吸液芯内流动方程组和管壁吸液芯温度方程组,得到吸液芯内压力和流速分布,以及管壁吸液芯温度分布;采用超松弛迭代法进行线性方程组的求解;
对于式(12)所示的方程组:
改写为式(13)的形式:
anm——控制方程左侧系数矩阵第n行第m列元素;bn——控制方程右侧向量第n行元素;
x(k)——解向量第k次迭代计算结果;
对得出结果进行修改,如式(14)所示:
式(14)中:
——迭代残向量;
对式(13)的结果进行变换可得;
式(15)中:
ω——松弛因子;
迭代时,预设迭代精度ε,判断是否满足如果满足,则停止迭代,;如果不满足,则重复此步骤,将x(k)代入式(15)计算得到x(k+1),直到n次迭代后满足停止迭代,得到吸液芯压力分布及温度分布的近似值;
步骤8:当前时间步计算完成,输出计算结果;通过计算结果更新物性,更新气液界面边界条件,进行下一时间步计算。
优选的,步骤7中,对于压力和流速求解构建的流动方程,迭代精度ε取10-4;对于温度求解构建的传热方程,迭代精度ε取10-7。
Claims (2)
1.一种碱金属热管吸液芯内部流速压力计算方法,其特征在于:步骤如下:
步骤1:确定碱金属热管几何结构,各部件材料及外部环境条件,建立热管模型:划分蒸发段、绝热段和冷凝段,根据热管内部状态划分管壁、吸液芯、蒸气区三部分计算域;蒸发段给定外壁温度或给定热流密度、绝热段给定零温度梯度、冷凝段给定换热系数和参考温度;设置热管整体温度初始条件,设置计算时间步长及计算总步数;
步骤2:根据吸液芯内部温度分布,确定吸液芯内工质液固两相分界面位置;通过气液界面转变温度确定蒸气区中连续流动气体与分子流态气体分界面位置;通过式(1)所示转变温度计算关系式,划分蒸气区流态;
式(1)中:
Ttrans——转变温度/K;
M——相对分子质量/g·mol-1;
Ru——通用气体常数/J·mol-1·K-1;
μ——动力粘度系数/Pa·s;
ρ——气相工质密度/kg·m-3;
D——蒸气区直径/m;
步骤3:计算气液界面温度分布:分子流态部分气液界面采用绝热边界,连续流动气体部分采用能量平衡边界;
通过式(2)所示的绝热边界计算气液界面温度:
式(2)中
T——吸液芯温度分布/K;
f——气液界面位置;
通过式(3)所示边界计算气液界面温度:
式(3)中:
kwi——吸液芯等效导热系数/W·m-2·K-1;
mi——气液界面质量通量/kg·m-2·s-1;
hfg——汽化潜热/J·kg-1;
步骤4:计算气液界面质量通量以及相变导致的能量通量,更新气液界面边界条件:
通过式(4)所示的拟合函数关系计算气液界面饱和压力:
式(4)中:
Pf——气液界面饱和压力/Pa;
A、B、C——拟合关系式系数;
通过式(5)、式(6)所示气体动力学关系式计算气液界面质量通量和能量通量:
式(5)中:
ε——吸液芯孔隙率;
——蒸发冷凝调节系数;
M——蒸气区蒸气的摩尔质量/(kg·mol-1);
Tli——气液界面上的节点液态工质温度/K;
Pli——气液界面上的节点液态工质饱和压力/Pa;
Tv——蒸气区蒸气的温度/K;
Pv——蒸气区蒸气的饱和压力/Pa;
qi——气液界面能量通量/J·m-2·s-1;
步骤5:构建吸液芯内液态工质流动物理现象及气液界面边界条件构建吸液芯质量、动量守恒方程,得到吸液芯内液态工质压力分布和流速分布描述方程组:
吸液芯内液态工质流动连续方程如式(7)所示:
式(7)中:
ρl——吸液芯内液态工质密度/kg·m-3;
t——时间/s;
——液态工质表观流速矢量/m·s-1;
A——气液界面上面积/m2;
吸液芯动量方程如式(8)所示:
式(8)中:
μl——液态工质动力粘度系数/Pa·s;
K——吸液芯结构渗透率/m2;
Pl——吸液芯内部相对压力分布/Pa;
步骤6:计算管壁区域和吸液芯区域温度分布,对于管壁使用纯导热控制方程,对于吸液芯采用流动换热控制方程;
吸液芯部分流动换热控制方程如式(9)所示:
式(13)中:
(ρCP)eff——等效质量热容/kg·m-3
T——热管温度分布/K
keff——等效导热系数/W·m-2·K-1;
对于管壁区域,流速设为零,各项物性均采用材料自身关系式计算;而对于吸液芯区域,采用Chi模型计算各项物性,物性计算关系式如式(10)、式(11)所示;
(ρc)eff=ερlcl+(1-ε)ρwscws (10)
keff——吸液芯的等效热导率/(W·m-1·K-1)
(ρc)eff——吸液芯内工质与多孔介质混合等效热容(J·kg-1·K-1);
cl——吸液芯内液态工质的热容/(J·m-3·K-1);
ρws——吸液芯材料的密度/(kg·m-3);
cws——吸液芯材料的热容/(J·m-3·K-1);
kl——吸液芯内液态工质的热导率/(W·m-1·K-1);
kws——吸液芯材料的热导率/(W·m-1·K-1);
步骤7:计算吸液芯内液态工质压力分布和流速分布,计算管壁及吸液芯温度分布:通过迭代求解吸液芯内流动方程组和管壁吸液芯温度方程组,得到吸液芯内压力和流速分布,以及管壁吸液芯温度分布;采用超松弛迭代法进行线性方程组的求解;
对于式(12)所示的方程组:
改写为式(13)的形式:
anm——控制方程左侧系数矩阵第n行第m列元素;
bn——控制方程右侧向量第n行元素;
x(k)——解向量第k次迭代计算结果;
对得出结果进行修改,如式(14)所示:
式(14)中:
ri (k+1)——迭代残向量;
对式(13)的结果进行变换得;
式(15)中:
ω——松弛因子;
迭代时,预设迭代精度ε,判断是否满足如果满足,则停止迭代,;如果不满足,则重复此步骤,将x(k)代入式(15)计算得到x(k+1),直到n次迭代后满足停止迭代,得到吸液芯压力分布及温度分布的近似值;
步骤8:当前时间步计算完成,输出计算结果;通过计算结果更新物性,更新气液界面边界条件,进行下一时间步计算。
2.根据权利要求1所述的一种碱金属热管吸液芯内部流速压力计算方法,其特征在于:步骤7中,对于压力和流速求解构建的流动方程,迭代精度ε取10-4;对于温度求解构建的传热方程,迭代精度ε取10-7。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310476490.1A CN116502555B (zh) | 2023-04-28 | 2023-04-28 | 一种碱金属热管吸液芯内部流速压力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310476490.1A CN116502555B (zh) | 2023-04-28 | 2023-04-28 | 一种碱金属热管吸液芯内部流速压力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116502555A true CN116502555A (zh) | 2023-07-28 |
CN116502555B CN116502555B (zh) | 2024-04-09 |
Family
ID=87326125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310476490.1A Active CN116502555B (zh) | 2023-04-28 | 2023-04-28 | 一种碱金属热管吸液芯内部流速压力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116502555B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015110599A1 (en) * | 2014-01-24 | 2015-07-30 | Ledaflow Technologies Da | Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines |
CN112528435A (zh) * | 2020-12-07 | 2021-03-19 | 西安交通大学 | 一种高温热管设计优化方法 |
CN113434988A (zh) * | 2021-06-11 | 2021-09-24 | 西安交通大学 | 一种核动力系统中热管的传热特性分析方法 |
CN114154438A (zh) * | 2021-12-07 | 2022-03-08 | 西安交通大学 | 一种碱金属热管冷态启动三阶段计算方法 |
WO2022067498A1 (zh) * | 2020-09-29 | 2022-04-07 | 中南大学 | 一种气液相变的介观模拟方法 |
-
2023
- 2023-04-28 CN CN202310476490.1A patent/CN116502555B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015110599A1 (en) * | 2014-01-24 | 2015-07-30 | Ledaflow Technologies Da | Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines |
WO2022067498A1 (zh) * | 2020-09-29 | 2022-04-07 | 中南大学 | 一种气液相变的介观模拟方法 |
CN112528435A (zh) * | 2020-12-07 | 2021-03-19 | 西安交通大学 | 一种高温热管设计优化方法 |
CN113434988A (zh) * | 2021-06-11 | 2021-09-24 | 西安交通大学 | 一种核动力系统中热管的传热特性分析方法 |
CN114154438A (zh) * | 2021-12-07 | 2022-03-08 | 西安交通大学 | 一种碱金属热管冷态启动三阶段计算方法 |
Non-Patent Citations (2)
Title |
---|
ZEQIN ZHANG 等: "Numerical investigation on startup characteristics of high temperature heat pipe for nuclear reactor", NUCLEAR ENGINEERING AND DESIGN, pages 1 - 11 * |
李文涛 等: "基于自适应松弛Picard法的高效非饱和渗流有限元分析", 岩土力学, pages 262 - 268 * |
Also Published As
Publication number | Publication date |
---|---|
CN116502555B (zh) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Czajkowski et al. | Experimental study on a large scale pulsating heat pipe operating at high heat loads, different adiabatic lengths and various filling ratios of acetone, ethanol, and water | |
Jafari et al. | Phase change heat transfer characteristics of an additively manufactured wick for heat pipe applications | |
CN112528435B (zh) | 一种高温热管设计优化方法 | |
Kemme | Ultimate heat-pipe performance | |
Jiang et al. | Flow and heat transfer characteristics of r22 and ethanol at supercritical pressures | |
Tang et al. | Experimental determination of the discharge coefficients for critical flow through an axisymmetric nozzle | |
CN106650120A (zh) | 一种液烃管道泄漏量测算方法 | |
Chen et al. | Verification at Mach 4 heat conditions of an annular microtube-typed precooler for hypersonic precooled engines | |
Cao et al. | Transient multidimensional analysis of nonconventional heat pipes with uniform and nonuniform heat distributions | |
Huang et al. | Thermal and hydraulic analysis of ultra-thin vapor chamber with copper columns considering Marangoni effect | |
Zimmermann et al. | An advanced conduction based heat pipe model accounting for vapor pressure drop | |
Li et al. | Theoretical and experimental investigations on the supercritical startup of a cryogenic axially Ω-shaped grooved heat pipe | |
CN116502555B (zh) | 一种碱金属热管吸液芯内部流速压力计算方法 | |
Sun et al. | Numerical investigation of nitrogen spontaneous condensation flow in cryogenic nozzles using varying nucleation theories | |
Fu et al. | Airside heat transfer analysis using Wilson plot method of three analogous serpentine tube heat exchangers for aero-engine cooling | |
Singh et al. | Numerical investigation of a nitrogen based cryogenic pulsating heat pipe | |
CN114154438A (zh) | 一种碱金属热管冷态启动三阶段计算方法 | |
Ma et al. | Capillary limit of a sodium screen-wick heat pipe | |
Chen et al. | Investigation of interface profiles in meshed wicks and related evaporation characteristics | |
CN113007482B (zh) | 一种长距离低温输送管线的自动降温方法 | |
Adams et al. | Drawn-polymer microchannel heat exchangers with application to recuperative cryocoolers | |
CN116541640B (zh) | 一种碱金属热管的归一化蒸气传热计算方法 | |
Naraghi et al. | A simple approach for thermal analysis of regenerative cooling of Rocket engines | |
Zimmermann et al. | A volume-of-fluid heat pipe model resolving pressure discontinuities at the wick-vapor interface | |
Mameli et al. | Towards quantitative validation of a closed loop pulsating heat pipe numerical model |
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 |