CN112464472A - 提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法 - Google Patents
提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法 Download PDFInfo
- Publication number
- CN112464472A CN112464472A CN202011340317.1A CN202011340317A CN112464472A CN 112464472 A CN112464472 A CN 112464472A CN 202011340317 A CN202011340317 A CN 202011340317A CN 112464472 A CN112464472 A CN 112464472A
- Authority
- CN
- China
- Prior art keywords
- heat exchange
- chf
- sodium
- grid
- steam generator
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000009835 boiling Methods 0.000 claims abstract description 43
- 239000007788 liquid Substances 0.000 claims abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 42
- 229910052708 sodium Inorganic materials 0.000 claims description 24
- 239000011734 sodium Substances 0.000 claims description 24
- 238000012546 transfer Methods 0.000 claims description 24
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 claims description 23
- 239000007791 liquid phase Substances 0.000 claims description 17
- 239000012530 fluid Substances 0.000 claims description 14
- 239000012808 vapor phase Substances 0.000 claims description 13
- 239000012071 phase Substances 0.000 claims description 10
- 239000003112 inhibitor Substances 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 238000009834 vaporization Methods 0.000 claims description 3
- 230000008016 vaporization Effects 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000036760 body temperature Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 230000008859 change Effects 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- KEAYESYHFKHZAL-UHFFFAOYSA-N Sodium Chemical compound [Na] KEAYESYHFKHZAL-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
本发明公开了一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,采用固定网格法划分钠冷快堆直流蒸汽发生器,采用固定网格法确定蒸汽发生器热构件网格划分,计算CHF发生的位置。判定各网格换热模式,针对单相液换热、核态沸腾、膜态沸腾和过热蒸汽工况,分别计算各网格换热系数。基于固定网格,使用CHF可移动边界模型对发生CHF现象的控制体再次网格划分,重新计算上、下热构件的换热系数。设置边界条件,重复以上步骤,分别计算每个网格下一时间步换热系数,最终求得每一时刻蒸汽发生器的换热量。本发明既保证钠冷快堆系统分析程序中直流蒸汽发生器换热计算精度,又提升反应堆数值计算程序的求解速度。
Description
技术领域
本发明属于核反应堆系统安全分析计算技术领域,具体涉及一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法。
背景技术
钠冷快堆由于冷却剂载热性能好、热惰性大、温度水平高,因此一般采用直流蒸汽发生器来减小反应堆尺寸,提高经济性。三回路直流蒸汽发生器管侧为水,经二回路液态钠加热,与壁面的换热模式包括过冷水对流、过冷沸腾、核态沸腾、膜态沸腾和过热蒸汽对流。各模式换热系数相差很大,尤其是CHF点前后的核态沸腾和膜态沸腾。固定网格内可能出现两种不同的换热模式,若只选取其中一种模式进行换热计算,会引起换热量的计算误差。同时,如果网格的换热模式发生变化,那么由于差分方程的变化会人为造成解的不连续。这种固定网格内换热模式转变引起的热量剧烈变换,将导致直流蒸汽发生器模拟的震荡问题。
目前直流蒸汽发生器换热计算方法有两种,一种是固定网格法,另一种是移动边界法。固定网格法是指在计算前将系统划分多个网格,在瞬态计算过程中网格不变,利用守恒方程计算流体的热力学状态。固定网格法适用于系统分析程序,但会造成解的不连续。可以通过增加网格的数量和减少时间不长的方法降低不连续性,但会降低计算效率。移动边界法中,网格边界随换热模式变化,因此网格界面是时间的连续函数,避免计算结果的突跳。同时,移动边界法可以减少网格数目,提高计算效率,但移动边界法不适用于系统分析程序。
钠冷快堆系统分析程序中直流蒸汽发生器仍采用传统固定网格法计算换热,计算效率和计算结果稳定性难以两全。亟待为钠冷快堆系统分析程序开发一种快速准确求解直流蒸汽发生器换热的方法,提升钠冷快堆系统程序瞬态计算的精度和速度。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,既保留系统程序的所有特点,又提升钠冷快堆系统程序瞬态计算的精度和速度。
本发明采用以下技术方案:
提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,包括以下步骤:
S1、根据蒸汽发生器水力部件控制体划分方式,采用固定网格法确定蒸汽发生器热构件网格划分,计算CHF位置;
S2、判定各网格换热模式,针对单相液换热、核态沸腾、膜态沸腾和过热蒸汽工况,分别计算各网格换热系数h;
S3、基于固定网格,使用CHF可移动边界模型对发生CHF现象的控制体再次网格划分,重新计算上、下热构件的换热系数;
S4、设置边界条件,重复以上步骤,分别计算每个网格下一时间步换热系数,最终求得每一时刻蒸汽发生器的换热量。
具体的,步骤S1中,CHF点在网格中的具体位置为:
其中,p表示钠侧,ΔL为CHF位置距前一个网格的距离,Tp,CHF为发生CHF时钠侧控制体温度,Ts为水侧饱和温度,Tp,i为钠侧控制体i的温度,U为综合传热系数,ri为传热管内半径,Wp为钠侧流体质量流量,Cp为钠的比热容。
具体的,步骤S2中,固定网格水侧换热模式判定采用平衡含汽率Xe和CHF位置,逻辑如下:
(1)Xe≤0,网格为单相对流换热;
(2)0<Xe<1,且网格在CHF发生前,网格为核态沸腾换热;
(3)0<Xe<1,且网格在CHF发生后,网格为膜态沸腾换热;
(4)Xe≥1,网格为过热蒸汽换热。
进一步的,水侧单相对流换热关系式根据流型选择,具体为:
对于层流,即Re<2300,选择Sellars关系式进行计算:
对于湍流,即Re≥2300,选取Dittus-Boelter公式进行计算:
Nu=0.023Re0.8Pr0.4
其中,kf为液相热导率;Dh为管道的水力直径,h为内表面传热系数,Nu为努塞尔数,Re为雷诺数,Pr为普朗特数。
进一步的,水侧核态沸腾换热采用Chen关系式,具体为:
hTP=Fhmac+Shmic
其中,hmac为宏观对流项换热系数;F为雷诺数因子;hmic为微观沸腾项换热系数;S为抑制因子。
更进一步的,对于饱和核态沸腾,雷诺数因子F计算如下:
对于过冷核态沸腾,雷诺数因子F计算如下:
其中,G为流体质量流速,ρf为液相密度,ρg为汽相密度,Gg为液相流量,Gf为汽相流量,μg为液相粘性系数,μf为汽相粘性系数,Ts为水侧饱和温度,Tf为液相温度。
进一步的,水侧膜态沸腾换热关系式根据流量进行选择,具体为:
G≤271.24kg/(m2·s),采用Bromely关系式计算如下:
G>271.24kg/(m2·s),采用Groneveld关系式计算如下:
其中,h为内表面传热系数,λg为汽相导热系数,g为重力加速度,ρg为汽相密度,μg为液相粘性系数,ρf为液相密度,Tw为壁面温度,Ts为水侧饱和温度,h'fg为汽化比焓,D为管道直径,G为流体质量流速,cpg为汽相比热容,Y为修正系数,x为平衡含汽率,De为当量直径。
进一步的,水侧过热蒸汽换热采用Sider-Tate关系式,具体为:
其中,Nu为努塞尔数,Re为雷诺数,Pr为普朗特数,μ为粘性系数,μw为壁面粘性系数。
具体的,步骤S2中,钠侧换热系数采用Westinghouse棒束关系式:
具体的,步骤S3中,采用移动边界法对发生CHF的网格再划分,将随时间变化的CHF位置作为移动边界,网格边界下为核态沸腾,边界上为膜态沸腾,并重新计算边界上、下区域各自的换热。
与现有技术相比,本发明至少具有以下有益效果:
本发明一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,使用固定网格加移动边界法计算直流蒸汽发生器的换热,保留了固定网格在系统分析程序中的优势,可以通过差分形式计算热工水力特性。以CHF位置作为移动边界,并只对发生CHF网格重新计算换热,极大地简化了单一移动边界法在系统程序中的编程困难问题,同时有效减弱传统固定网格法造成的计算结果不稳定性,对移动边界上下区域重新计算换热,提高粗网格的计算精度,同时提高了计算效率。
进一步的,确定CHF点在网格中的具体位置,热构件的网格划分依据水力部件控制体划分方式,满足钠冷快堆系统分析程序划分方式。
进一步的,水侧换热逻辑依据平衡含汽率Xe和CHF位置,简化了固定网格水侧换热模式判定过程。
进一步的,通过水侧单相对流换热关系式考虑了层流和湍流的影响。
进一步的,水侧核态沸腾采用更精确的Chen关系式。
进一步的,水侧膜态沸腾关系式分别考虑了高低流量的影响。
进一步的,水侧过热蒸汽关系式可通过粘性系数比、Re和Pr直接求得。
进一步的,钠侧换热系数关系式考虑了棒束,且棒束可以任意排布。
进一步的,在固定网格下再次划分移动边界,可以增加网格的连续性。
综上所述,本发明采用固定网格划分蒸汽发生器热构件,并在发生CHF的网格划分可移动边界,针对不同的工况,壁面换热关系式采用了不同的计算方式,最后可以更快地得到更精确的钠冷快堆蒸汽发生器的壁面换热量。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明示意图,其中(a)为固定网格划分,(b)为固定网格加移动边界划分;
图2为本发明流程图;
图3为CHF轴向位置随时间变化的曲线图。
具体实施方式
本发明提供了一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,采用固定网格法划分钠冷快堆直流蒸汽发生器,采用固定网格法确定蒸汽发生器热构件网格划分,计算CHF发生的位置。判定各网格换热模式,针对单相液换热、核态沸腾、膜态沸腾和过热蒸汽工况,分别计算各网格换热系数。基于固定网格,使用CHF可移动边界模型对发生CHF现象的控制体再次网格划分,重新计算上、下热构件的换热系数。设置边界条件,重复以上步骤,分别计算每个网格下一时间步换热系数,最终求得每一时刻蒸汽发生器的换热量。既保证钠冷快堆系统分析程序中直流蒸汽发生器换热计算精度,又提升反应堆数值计算程序的求解速度。
请参阅图2,本发明一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,包括以下步骤:
S1、根据蒸汽发生器水力部件控制体划分方式,采用固定网格法确定蒸汽发生器热构件网格划分,计算CHF位置;
首先建立蒸汽发生器正常运行时一、二次侧工质换热量计算公式:
其中,下标i为管道内工质参数;下标o为管道外工质参数;T为流体温度/K;h为内表面传热系数/J·m-2·K-1;l为控制体长度/m;λ为传热管热导率/J·m-1·K-1。
对上式进行进一步简化:
其中,
根据上式写出发生CHF时,蒸汽发生器一、二次侧工质换热量的计算关系式:
其中,U为综合传热系数/J·m-2·K-1;Tp.CHF为发生CHF时钠侧控制体温度/K;Ts为水侧饱和温度/K。
其中,综合传热系数U为:
当蒸汽发生器达到稳定状态时,单位时间内钠侧工质能量的减少等于水侧工质能量增加,如下式所示:
CpWp(Tp.CHF-Tpi)=Ws(Hs,CHF-Hsi)
其中,Cp为钠的比热容/J·kg-1·K-1;Tp.CHF为发生CHF时钠侧温度/K;Wp为钠侧流体质量流量/kg·s-1;Tpi为钠侧控制体i的温度/K;Ws为钠侧工质质量流量/kg·s-1;Hs,CHF为发生CHF时水侧工质焓值/J·kg-1;Hsi为水侧控制体i的焓值/J·kg-1。
水在两相区的焓值变化还可以用含汽率来表示:
Ws(HCHF-Hi)=WsHfg(xCHF-xi)
水侧含汽率与钠侧温度的关系:
WpCp(Tp,CHF-Tp,i)=WsHfg(xCHF-xi)
其中,Hfg为水的汽化潜热/J·kg-1;xCHF为发生CHF时水侧的含汽率;xi为水侧控制体i的含汽率。
对上式进行化简:
可以得到发生CHF时,水侧含汽率关于钠侧温度的关系式:
根据CHF查询表,临界热流密度qCHF与含汽率xCHF,流体压力和流体质量流速G有关,如下式所示:
qchf=f(G,xchf,P)
迭代求解得到临界热流密度qCHF,计算出CHF点处的钠侧温度Tp.CHF。
将换热计算关系式转成微分形式:
dQ=WpCpdT=-2πU(Tp-Ts)ridL
并沿蒸汽发生器轴向进行积分,得:
其中,T1、T2为蒸汽发生器钠侧两端温度/K;Ts为水侧饱和温度/K;Wp为钠侧质量流量/kg·s-1;ri为传热管内半径/m;U为蒸汽发生器传热管一、二次侧之间的综合传热系数/W·(m2·K)-1;Cp为钠侧比热容/J·(kg·K)-1;ΔL为水侧两端距离/m。
如果T1=Tp.chf,T2=Tp.i,那么ΔL就是CHF发生位置与前一个控制体i之间的距离,这样就可以计算出CHF实际位置。接下来将发生CHF的控制体进行进一步网格划分,细化该控制体的换热过程,从而提高计算精度。
请参阅图1,(a)图为采用固定网格划分的节点,CHF边界只能在控制体边界上,(b)图采用固定网格和可移动边界划分节点,CHF边界可以出现在控制体内。
S2、针对不同工况,换热关系式采用不同的计算方式;
对钠侧采用Westinghouse棒束关系式:
对水侧,则进行换热模式判断,再根据不同工况选取换热系数计算公式。
(1)Xe≤0,网格为单相对流换热;
对于层流,即Re<2300,选择Sellars关系式进行计算:
其中,kf为液相热导率/J·(m·K)-1;Dh为管道的水力直径/m。
对于湍流,即Re≥2300,选取Dittus-Boelter公式进行计算:
Nu=0.023Re0.8Pr0.4
(2)0<Xe<1,且网格在CHF发生前,网格为核态沸腾换热;
核态沸腾换热采用Chen关系式:
hTP=Fhmac+Shmic
其中,hmac为宏观对流项换热系数/J·m-2·K-1;F为雷诺数因子;hmic为微观沸腾项换热系数/J·m-2·K-1;S为抑制因子。
其中,宏观对流换热系数hmac计算公式为:
微观微观沸腾项换热系数为hmic计算公式为:
其中,Cpf为液相比热容/J·(kg·K)-1;ρf为液相密度/kg·m-3;σ为流体表面张力/N·m-1;ρg为汽相密度/kg·m-3。
对于饱和核态沸腾,雷诺数因子F计算关系式如下:
其中,G为流体质量流速/kg·m-2·s-1。
对于过冷核态沸腾,雷诺数因子F计算关系式如下:
抑制因子计算关系式如下:
其中,
Retp=min(10-4RefF1.25,70)
(3)0<Xe<1,且网格在CHF发生后,网格为膜态沸腾换热;
膜态沸腾换热关系式根据流量进行选择。
对小流量(G≤271.24kg/(m2·s))采用Bromely关系式:
对大流量(G>271.24kg/(m2·s))采用Groneveld关系式:
(4)Xe≥1,网格为过热蒸汽换热。
水侧过热蒸汽换热采用Sider-Tate关系式如下:
S3、基于固定网格,使用CHF可移动边界模型对发生CHF现象的控制体再次网格划分,将随时间变化的CHF位置作为移动边界,网格边界下为核态沸腾,边界上为膜态沸腾,并重新计算边界上、下区域各自的换热,关系式选择方式与步骤S2中相同。
S4、设置边界条件,重复以上步骤,分别计算每个网格下一时间步换热系数,最终求得每一时刻蒸汽发生器的换热量。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中的描述和所示的本发明实施例的组件可以通过各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
根据ETEC实验的相关参数,选取了A2、A5、A7和A11实验工况采用NUSOL-LMR程序进行了模拟。表1给出了水侧出口温度在是否考虑CHF可移动边界模型前后的对比。表2为考虑CHF可移动边界模型前、后程序计算ETEC稳态工况的水侧进出口焓差与实验测量值的对比表。表3为本程序分别使用固定网格和CHF可移动边界模型计算ETEC蒸汽发生器时在不同网格数目下计算精度和计算时间的对比表。
表1考虑CHF可移动边界模型前、后程序计算的水侧出口温度对比
表2考虑CHF可移动边界模型前、后程序计算的水侧进出口焓差与实验值的对比
表3在不同网格数目下使用两种网格划分方法计算ETEC蒸汽发生器的结果对比表
通过将本发明方法的结果与实验结果进行对比,使用了CHF可移动边界法之后,程序计算结果与实验测量值之间的误差得到减小,尤其在控制体中CHF发生的位置靠后,使用CHF可移动边界模型后,控制体的换热量得到较大提升,蒸汽发生器的计算结果有较大改进。增加蒸汽发生器轴向节点数目和使用CHF可移动边界模型都可以提高程序的计算结果的准确性,但CHF可移动边界模型可以提高程序的计算速度。
图3为固定网格和CHF可移动边界模型下程序计算ETEC蒸汽发生器一次侧工质入口温度由450K升高到770K时CHF轴向位置随时间变化曲线。该瞬态工况下固定网格中CHF位置在控制体边界上发生移动时会在较大范围内发生反复振荡。而使用了CHF可移动边界模型之后,因为把CHF边界定义在了控制体内部,所以当蒸汽发生器工质参数发生变化时,程序可以在每一个时间步长上都计算出CHF发生时其边界在控制体中的具体位置,CHF边界在控制体中的位置变化是连续的,不会发生突跳。因此CHF可移动边界模型会增加程序的稳定性,不易发生振荡。
综上所述,本发明一种提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,结合了系统程序固定网格划分方式,采用CHF可移动边界模型计算壁面换热,既保留了系统程序建模方式,同时又能提高蒸汽发生器换热计算的精度、稳定性和速度。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (10)
1.提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,其特征在于,包括以下步骤:
S1、根据蒸汽发生器水力部件控制体划分方式,采用固定网格法确定蒸汽发生器热构件网格划分,计算CHF位置;
S2、判定各网格换热模式,针对单相液换热、核态沸腾、膜态沸腾和过热蒸汽工况,分别计算各网格换热系数h;
S3、基于固定网格,使用CHF可移动边界模型对发生CHF现象的控制体再次网格划分,重新计算上、下热构件的换热系数;
S4、设置边界条件,重复以上步骤,分别计算每个网格下一时间步换热系数,最终求得每一时刻蒸汽发生器的换热量。
3.根据权利要求1所述的提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,其特征在于,步骤S2中,固定网格水侧换热模式判定采用平衡含汽率Xe和CHF位置,逻辑如下:
(1)Xe≤0,网格为单相对流换热;
(2)0<Xe<1,且网格在CHF发生前,网格为核态沸腾换热;
(3)0<Xe<1,且网格在CHF发生后,网格为膜态沸腾换热;
(4)Xe≥1,网格为过热蒸汽换热。
5.根据权利要求3所述的提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,其特征在于,水侧核态沸腾换热采用Chen关系式,具体为:
hTP=Fhmac+Shmic
其中,hmac为宏观对流项换热系数;F为雷诺数因子;hmic为微观沸腾项换热系数;S为抑制因子。
10.根据权利要求1所述的提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法,其特征在于,步骤S3中,采用移动边界法对发生CHF的网格再划分,将随时间变化的CHF位置作为移动边界,网格边界下为核态沸腾,边界上为膜态沸腾,并重新计算边界上、下区域各自的换热。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011340317.1A CN112464472B (zh) | 2020-11-25 | 2020-11-25 | 提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011340317.1A CN112464472B (zh) | 2020-11-25 | 2020-11-25 | 提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464472A true CN112464472A (zh) | 2021-03-09 |
CN112464472B CN112464472B (zh) | 2024-04-02 |
Family
ID=74808181
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011340317.1A Active CN112464472B (zh) | 2020-11-25 | 2020-11-25 | 提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464472B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113686917A (zh) * | 2021-08-17 | 2021-11-23 | 西安交通大学 | 一种液态铅铋合金7棒束组件流动换热实验装置及实验方法 |
CN115169262A (zh) * | 2022-07-18 | 2022-10-11 | 重庆大学 | 一种考虑汽液相间温差变化的汽泡分析方法 |
CN116030906A (zh) * | 2023-03-28 | 2023-04-28 | 西安交通大学 | 一种铅铋快堆多组分流体的物性参数计算方法 |
CN116362155A (zh) * | 2023-03-22 | 2023-06-30 | 西安交通大学 | 一种液态金属直流蒸汽发生器腔室换热系数计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020476A (zh) * | 2019-04-08 | 2019-07-16 | 西安交通大学 | 一种反应堆u型管式蒸汽发生器全三维耦合模型建立方法 |
CN111680417A (zh) * | 2020-06-03 | 2020-09-18 | 西安交通大学 | 一种适用于钠水直流蒸汽发生器的热工设计方法 |
CN111680458A (zh) * | 2020-06-03 | 2020-09-18 | 西安交通大学 | 一种适用于钠水直流蒸汽发生器的热工水力瞬态计算方法 |
CA3080721A1 (en) * | 2019-05-14 | 2020-11-14 | Cenovus Energy Inc. | Methods and systems for two-stage steam generation |
-
2020
- 2020-11-25 CN CN202011340317.1A patent/CN112464472B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110020476A (zh) * | 2019-04-08 | 2019-07-16 | 西安交通大学 | 一种反应堆u型管式蒸汽发生器全三维耦合模型建立方法 |
CA3080721A1 (en) * | 2019-05-14 | 2020-11-14 | Cenovus Energy Inc. | Methods and systems for two-stage steam generation |
CN111680417A (zh) * | 2020-06-03 | 2020-09-18 | 西安交通大学 | 一种适用于钠水直流蒸汽发生器的热工设计方法 |
CN111680458A (zh) * | 2020-06-03 | 2020-09-18 | 西安交通大学 | 一种适用于钠水直流蒸汽发生器的热工水力瞬态计算方法 |
Non-Patent Citations (2)
Title |
---|
叶尚尚;杨红义;刘一哲;杨晓燕;王晓坤;齐少璞;王利霞;: "钠冷快堆直流式蒸汽发生器固定网格与滑移网格模型对比分析", 原子能科学技术, no. 07 * |
叶桂珍, 王洲, 杨献勇, 罗锐: "快堆螺旋管式直流蒸汽发生器的热工水力分析模型", 清华大学学报(自然科学版), no. 02 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113686917A (zh) * | 2021-08-17 | 2021-11-23 | 西安交通大学 | 一种液态铅铋合金7棒束组件流动换热实验装置及实验方法 |
CN115169262A (zh) * | 2022-07-18 | 2022-10-11 | 重庆大学 | 一种考虑汽液相间温差变化的汽泡分析方法 |
CN115169262B (zh) * | 2022-07-18 | 2023-06-06 | 重庆大学 | 一种考虑汽液相间温差变化的汽泡分析方法 |
CN116362155A (zh) * | 2023-03-22 | 2023-06-30 | 西安交通大学 | 一种液态金属直流蒸汽发生器腔室换热系数计算方法 |
CN116362155B (zh) * | 2023-03-22 | 2024-01-30 | 西安交通大学 | 一种液态金属直流蒸汽发生器腔室换热系数计算方法 |
CN116030906A (zh) * | 2023-03-28 | 2023-04-28 | 西安交通大学 | 一种铅铋快堆多组分流体的物性参数计算方法 |
CN116030906B (zh) * | 2023-03-28 | 2023-06-27 | 西安交通大学 | 一种铅铋快堆多组分流体的物性参数计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112464472B (zh) | 2024-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112464472B (zh) | 提高钠冷快堆系统程序中蒸汽发生器换热计算性能的方法 | |
CN109657372B (zh) | 一种新型的管壳式换热器多尺度耦合协同换热模拟方法 | |
Li et al. | Modeling and dynamic simulation of a steam generation system for a parabolic trough solar power plant | |
CN104834773A (zh) | 一种直管式直流蒸汽发生器换热性能的仿真方法 | |
CN111680417B (zh) | 一种适用于钠水直流蒸汽发生器的热工设计方法 | |
CN114398840B (zh) | 一种蒸发器通用水动力计算方法、系统、设备及存储介质 | |
CN118246201A (zh) | 涡轮叶片多孔介质液-气相变发汗冷却的高鲁棒性数值计算方法 | |
CN113935259A (zh) | 一种反应堆一回路系统瞬态水锤波的计算方法及系统 | |
Wu et al. | A three-region movable-boundary helical coil once-through steam generator model for dynamic simulation and controller design | |
CN116757108B (zh) | 换热器的模拟方法、装置及设备 | |
CN117371343A (zh) | 一种用于高温气冷堆换热装置温度分布的分析方法及系统 | |
Galindo-García et al. | CFD simulations of heat recovery steam generators including tube banks | |
Wang et al. | Effect of shell on performance of heat exchanger with trefoil-hole baffle | |
Souyri et al. | Pressurized water reactor modelling with Modelica | |
CN112446097A (zh) | 一种蒸汽发生器体积和负荷的多目标优化方法 | |
Louw | Discrete and porous computational fluid dynamics modelling of an air-rock bed thermal energy storage system | |
Wang et al. | Optimal design of supercritical He–H2 PCHE in SABER system by multi-objective genetic algorithm | |
Yuan et al. | Numerical Simulation of Flow Boiling Heat Transfer in Helical Tubes Under Marine Conditions | |
Wang et al. | Numerical investigation on thermal hydraulic and transit time characteristics of density wave oscillations | |
Ding et al. | Experimental and Numerical Study on Convective Heat Transfer Characteristic in the Turbulent Region of Molten Salt in Shell-Side of Shell and Tube Heat Exchanger | |
CN114896910B (zh) | 流动沸腾传热精细化模拟方法 | |
Dong et al. | Relap5 MOD3. 2 Lead and LBE Model Addition and Program Validation | |
Ghadge | Analytical and numerical study of instability of supercritical water flowing upward in two heated parallel channels with wall thermal energy storage effects | |
CN118313306A (zh) | Lng船航行中储罐内lng传热和日蒸发率预测方法 | |
Chenyang et al. | Numerical Simulation Study of Steam Condensation in Vertical Tube |
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 |