CN106203722A - 一种预测核反应堆并联通道管间脉动的方法 - Google Patents
一种预测核反应堆并联通道管间脉动的方法 Download PDFInfo
- Publication number
- CN106203722A CN106203722A CN201610566356.0A CN201610566356A CN106203722A CN 106203722 A CN106203722 A CN 106203722A CN 201610566356 A CN201610566356 A CN 201610566356A CN 106203722 A CN106203722 A CN 106203722A
- Authority
- CN
- China
- Prior art keywords
- parallel port
- conservation equation
- pulsation
- tubes
- calculate
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/02—Devices or arrangements for monitoring coolant or moderator
- G21C17/032—Reactor-coolant flow measuring or monitoring
-
- 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)
- Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Plasma & Fusion (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种预测核反应堆并联通道管间脉动的方法,步骤如下:给定各并联通道的几何及热工水力边界参数;初始化各并联通道的各相计算初值;计算守恒方程源项;对守恒方程进行离散求解;对各通道进行流量分配并判断收敛;判断计算是否结束计算;利用频谱分析方法获得各通道流量并判定管间脉动是否发生;本发明基于两流体模型,能够精确描述两相流动行为;利用频谱分析方法可以准确获得各通道流量的频域信息,从而精确预测核反应堆并联通道管间脉动现象,对核反应堆安全分析具有重要意义。
Description
技术领域
本发明属于核反应堆安全分析技术领域,具体涉及一种预测核反应堆并联通道管间脉动的方法。
技术背景
随着核能的广泛利用,科学有效地评估核反应堆的安全性越来越重要,对反应堆系统中流动不稳定性现象的研究越来越多。反应堆系统中存在大量的并联冷却剂通道,通道间相互作用会引起各通道流量的震荡,即管间脉动现象。管间脉动发生时,尽管并联通道总流量可能保持不变,某些通道的流量可能发生剧烈变化,影响反应堆的安全运行。核反应堆运行过程中,通道间的局部流量震荡难以准确测量,因此,预测核反应堆并联通道管间脉动对于核安全有着重要意义。
管间脉动是一种两相流动不稳定现象。目前国内外的管间脉动分析中存在两方面问题。第一,计算管间脉动所采用的两相流动模型多基于均相流模型。由于均相流模型将两相流动等效为单相流动,无法对复杂的两相流动不稳定性现象进行准确的描述。第二,现有的两相热工水力分析程序无法专门针对管间脉动特性进行计算分析,难以直接应用于反应堆并联通道管间脉动的预测。
基于以上背景技术,需要发明一种准确描述两相特性,针对管间脉动现象的分析计算方法,使得其能够准确预测核反应堆并联通道管间脉动现象,对反应堆实际运行和安全分析评价提供参考。
发明说明
为解决上述现有的管间脉动分析方法存在的问题,本发明提出了一种基于两流体模型的并联通道管间脉动计算分析方法,既能够对并联通道管间脉动的两相流动行为进行准确计算,又能对管间脉动的特性进行分析。
本发明的采用技术方案为:
一种预测核反应堆并联通道管间脉动的方法,包括如下步骤:
步骤一:给定各并联通道的几何及热工水力边界参数:其中几何参数包括各并联通道的流通面积、高度、水力直径和热力直径,热工水力边界参数包括各并联通道的出口压力、进口流速、进口比焓和热功率:给定计算步长和计算总时间;
步骤二:初始化各并联通道的各相计算初值,包括各并联通道的出口压力、进口流速和进口比焓;
步骤三:根据当前时层的出口压力、进口流速和进口比焓,计算各并联通道的密度、温度以及质量守恒方程、动量守恒方程及能量守恒方程的源项;
步骤四:对质量守恒方程、动量守恒方程及能量守恒方程进行离散,并对各并联通道的压力、流速、密度和比焓进行求解;
步骤五:通过判断各通道压降是否相等的方法判定各并联通道的流量分配是否收敛;是则进行下一步,否则重新进行各并联通道的流量分配后,重复步骤三;
步骤六:根据给定的计算步长和计算总时间判断是否完成计算;是则进行下一步,否则计算下一计算时间,返回步骤三进行下一时层的计算;
步骤七:使用频谱分析方法将计算时间内的各通道的流量时域信息变化为频域信息,即可获得各并联通道流量变化的频率、周期及相位;当各通道流量出现反相时,即可判定发生了管间脉动。
与现有方法相比,本发明有如下优点:
1.使用了两流体模型,通过求解各相的守恒方程,对两相行为进行精确的描述。
2.使用频谱分析方法对两流体模型求得的流量信息进行分析,可以准确获得管间脉动的细节信息。
附图说明
图1为本发明的计算流程图。
图2为一种并联通道结构示意图。
图3为一种管间脉动发生时的各通道流量变化的时域信息。
图4为一种管间脉动发生时的各通道流量变化的频域信息。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细的描述。
如图1所示,本发明一种预测核反应堆并联通道管间脉动的方法,包括如下步骤:
步骤一:取一个典型的压水堆并联双通道,如图2所示。各通道直径10mm,高度1m。给定压力为15Mpa,进口流速为1m/s,进口比焓为1040kJ/kg,热功率为258.3W。给定计算步长为0.01s,计算总时间为20s。
步骤二:本例中,各并联通道的各相计算初值化数值采用步骤一中给定的边界值。
步骤三:本例中,密度采用IAPWS-IF97水物性计算。采用的质量守恒方程为
动量守恒方程为
能量守恒方程为
式中,上标k表示为气相或液相,气相用g表示,液相用f表示。以下为本例中计算源项的方法。为壁面摩擦引起的动量源项,本例中采用Martinelli-Lockhart方法计算两相摩擦阻力。为相间摩擦阻力,本例中采用EPRI方法计算两相相间摩擦阻力。为相间换热,采用Plesset-Zwick方法计算两相之间的换热量。
步骤四:本例对守恒方程采用有限体积法进行离散,采用半隐式压力速度耦合方法进行求解。
步骤五:采用判断各通道压降是否相等的方法判定流量分配是否收敛。若不相等,流量分配的方法为
式中为各通道各相经过分配的入口速度,为上一时层计算的速度,为各通道进口压力,Pa为平均压力。若流量分配不收敛,返回步骤三计算,直至收敛,即得到当前时间的流量。
步骤六:根据给定的计算步长,计算总时间判断是否完成计算。若未完成计算,将当前时间加上给定的时间步长,返回步骤三进行下一时层的计算。直至计算完成。
步骤七:本例中计算得到的流量时域信号如图3所示。使用快速傅里叶变换方法进行频域分析,结果如图4所示。两个通道发生了相位相反,即判定发生了管间脉动。
以上实施例预测了上述形式的反应堆并联通道在上述给定的参数下,会发生管间脉动。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (3)
1.一种预测核反应堆并联通道管间脉动的方法,其特征在于:包括如下步骤:
步骤一:给定各并联通道的几何及热工水力边界参数;其中几何参数包括各并联通道的流通面积、高度、水力直径和热力直径,热工水力边界参数包括各并联通道的出口压力、进口流速、进口比焓和热功率;给定计算步长和计算总时间;
步骤二:初始化各并联通道的各相计算初值,包括各并联通道的出口压力、进口流速和进口比焓;
步骤三:根据当前时层的出口压力、进口流速和进口比焓,计算各并联通道的密度、温度以及质量守恒方程、动量守恒方程及能量守恒方程的源项;
步骤四:对质量守恒方程、动量守恒方程及能量守恒方程进行离散,并对各并联通道的压力、流速、密度和比焓进行求解;
步骤五:通过判断各通道压降是否相等的方法判定各并联通道的流量分配是否收敛;是则进行下一步,否则重新进行各并联通道的流量分配后,重复步骤三;
步骤六:根据给定的计算步长和计算总时间判断是否完成计算;是则进行下一步,否则计算下一计算时间,返回步骤三进行下一时层的计算;
步骤七:使用频谱分析方法将计算时间内的各通道的流量时域信息变化为频域信息,即可获得各并联通道流量变化的频率、周期及相位;当各通道流量出现反相时,即可判定发生了管间脉动。
2.根据权利要求1所述的一种预测核反应堆并联通道管间脉动的方法,其特征在于:步骤三所述质量守恒方程、动量守恒方程及能量守恒方程的源项采用两流体模型。
3.根据权利要求1所述的一种预测核反应堆并联通道管间脉动的方法,其特征在于:步骤四所述质量守恒方程、动量守恒方程及能量守恒方程的求解过程采用两流体模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610566356.0A CN106203722B (zh) | 2016-07-18 | 2016-07-18 | 一种预测核反应堆并联通道管间脉动的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610566356.0A CN106203722B (zh) | 2016-07-18 | 2016-07-18 | 一种预测核反应堆并联通道管间脉动的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106203722A true CN106203722A (zh) | 2016-12-07 |
CN106203722B CN106203722B (zh) | 2019-07-02 |
Family
ID=57493974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610566356.0A Active CN106203722B (zh) | 2016-07-18 | 2016-07-18 | 一种预测核反应堆并联通道管间脉动的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106203722B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107421850A (zh) * | 2017-07-31 | 2017-12-01 | 中国核动力研究设计院 | 一种并联通道密度波脉动起始点识别方法 |
CN107436274A (zh) * | 2017-07-31 | 2017-12-05 | 中国核动力研究设计院 | 多重流量脉动叠加条件下并联通道密度波脉动识别方法 |
CN111081400A (zh) * | 2019-11-25 | 2020-04-28 | 西安交通大学 | 一种基于机器学习的棒束子通道热工水力特性预测方法 |
CN113551870A (zh) * | 2021-06-23 | 2021-10-26 | 中国核动力研究设计院 | 一种对并联通道流动失稳行为影响机制的表征方法及系统 |
CN113792449A (zh) * | 2021-08-10 | 2021-12-14 | 中国船舶重工集团公司第七一九研究所 | 核反应堆监测方法及系统 |
CN114792074A (zh) * | 2022-04-25 | 2022-07-26 | 南京航空航天大学 | 一种流动不稳定性影响区域划分的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103219055A (zh) * | 2013-04-09 | 2013-07-24 | 上海交通大学 | 用于检测核反应堆中间隙流动不稳定现象的装置 |
CN104133958A (zh) * | 2014-07-28 | 2014-11-05 | 浙江中控软件技术有限公司 | 一种复杂管网模拟仿真计算方法及装置 |
CN105247622A (zh) * | 2013-04-10 | 2016-01-13 | 阿海珐核能公司 | 用于模拟核反应堆中容器内流体流动的方法以及用于计算核反应堆堆芯组件的机械形变的方法,以及相关的计算机程序产品 |
-
2016
- 2016-07-18 CN CN201610566356.0A patent/CN106203722B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103219055A (zh) * | 2013-04-09 | 2013-07-24 | 上海交通大学 | 用于检测核反应堆中间隙流动不稳定现象的装置 |
CN105247622A (zh) * | 2013-04-10 | 2016-01-13 | 阿海珐核能公司 | 用于模拟核反应堆中容器内流体流动的方法以及用于计算核反应堆堆芯组件的机械形变的方法,以及相关的计算机程序产品 |
CN104133958A (zh) * | 2014-07-28 | 2014-11-05 | 浙江中控软件技术有限公司 | 一种复杂管网模拟仿真计算方法及装置 |
Non-Patent Citations (2)
Title |
---|
左娟莉等: "液态金属内单个气泡上升行为的MPS法数值模拟", 《原子能科学技术》 * |
邓大祥: "微尺度热质输运强化槽道多孔结构制造及性能研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107421850A (zh) * | 2017-07-31 | 2017-12-01 | 中国核动力研究设计院 | 一种并联通道密度波脉动起始点识别方法 |
CN107436274A (zh) * | 2017-07-31 | 2017-12-05 | 中国核动力研究设计院 | 多重流量脉动叠加条件下并联通道密度波脉动识别方法 |
CN107436274B (zh) * | 2017-07-31 | 2018-05-08 | 中国核动力研究设计院 | 多重流量脉动叠加条件下并联通道密度波脉动识别方法 |
CN107421850B (zh) * | 2017-07-31 | 2018-06-12 | 中国核动力研究设计院 | 一种并联通道密度波脉动起始点识别方法 |
CN111081400A (zh) * | 2019-11-25 | 2020-04-28 | 西安交通大学 | 一种基于机器学习的棒束子通道热工水力特性预测方法 |
CN113551870A (zh) * | 2021-06-23 | 2021-10-26 | 中国核动力研究设计院 | 一种对并联通道流动失稳行为影响机制的表征方法及系统 |
CN113792449A (zh) * | 2021-08-10 | 2021-12-14 | 中国船舶重工集团公司第七一九研究所 | 核反应堆监测方法及系统 |
CN114792074A (zh) * | 2022-04-25 | 2022-07-26 | 南京航空航天大学 | 一种流动不稳定性影响区域划分的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106203722B (zh) | 2019-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106203722B (zh) | 一种预测核反应堆并联通道管间脉动的方法 | |
Chu et al. | Study on hydraulic and thermal performance of printed circuit heat transfer surface with distributed airfoil fins | |
Mahmoudi et al. | Analytical investigation of heat transfer enhancement in a channel partially filled with a porous material under local thermal non-equilibrium condition: effects of different thermal boundary conditions at the porous-fluid interface | |
Yaïci et al. | Three-dimensional unsteady CFD simulations of a thermal storage tank performance for optimum design | |
Dang et al. | In-tube cooling heat transfer of supercritical carbon dioxide. Part 1. Experimental measurement | |
Koohi-Fayegh et al. | An analytical approach to evaluating the effect of thermal interaction of geothermal heat exchangers on ground heat pump efficiency | |
Wang et al. | CFD simulation of flow and heat transfer characteristics in a 5× 5 fuel rod bundles with spacer grids of advanced PWR | |
CN108119761A (zh) | 一种气液两相流下天然气管道腐蚀程度预测方法 | |
CN102419827B (zh) | 基于rbf神经网络的沸腾换热预测方法 | |
Dehghan et al. | A new 1D analytical model for investigating the long term heat transfer rate of a borehole ground heat exchanger by Green's function method | |
Farhadi | Analysis of flow coastdown for an MTR-pool type research reactor | |
Li et al. | Dryout prediction with CFD model of annular two-phase flow | |
Gómez et al. | Numerical study of the thermal behaviour of a water heater tank with a corrugated coil | |
CN109783972B (zh) | 基于流固耦合分析计算的止回阀内泄漏流量的监测方法 | |
Adamsson et al. | Transient dryout prediction using a computationally efficient method for sub-channel film-flow analysis | |
Avery et al. | Critical heat flux in TRIGA-fueled reactors cooled by natural convection | |
Lu et al. | Large-eddy simulation of thermal stratification in a straight branch of a tee junction with or without leakage | |
Hashemi-Tilehnoee et al. | Benchmarking a sub-channel program based on a drift-flux model with 8× 8 NUPEC BWR rod bundle | |
Chen et al. | Modeling and experimental research on ground-source heat pump in operation by neural network | |
Talley et al. | Horizontal bubbly flow with elbow restrictions: Interfacial area transport modeling | |
Yang et al. | Les of a Space-Developing Natural Convection Boundary Layer in Water Along a Vertical Wall with High Heat Flux | |
Yoon et al. | Validation of a CFD analysis model for predicting CANDU-6 moderator temperature against SPEL experiments | |
Gong et al. | Analysis of the heat-carrying performance of natural circulation loop using modified RELAP5/MOD3. 4 | |
Laurien et al. | Two-layer model for the heat transfer to supercritical CO2 | |
Hauptmann et al. | Axial development of unusual velocity profiles due to heat transfer in variable density fluids |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |