CN114093432A - 核反应堆事故工况下耦合传热传质的包壳氧化分析方法 - Google Patents
核反应堆事故工况下耦合传热传质的包壳氧化分析方法 Download PDFInfo
- Publication number
- CN114093432A CN114093432A CN202111400921.3A CN202111400921A CN114093432A CN 114093432 A CN114093432 A CN 114093432A CN 202111400921 A CN202111400921 A CN 202111400921A CN 114093432 A CN114093432 A CN 114093432A
- Authority
- CN
- China
- Prior art keywords
- phase
- interface
- oxygen atom
- equation
- cladding
- 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
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/90—Programming languages; Computing architectures; Database systems; Data warehousing
-
- 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
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种核反应堆事故工况下耦合传热传质的包壳氧化分析方法,步骤如下:1.建立包含β‑Zr,α‑Zr(O)和ZrO2的包壳几何结构;2.建立计算模型的控制方程,包括氧原子扩散方程、锆原子守恒方程和导热方程;3.建立相界面的迁移方程;4.设置相界面处的边界条件;5.对包壳几何划分网格,离散控制方程,通过迭代求解获得下一个时刻收敛的温度场、氧原子浓度场、锆原子浓度场和变形速度场;6.计算在此次时间步长内相界面的位移,更新包壳的几何结构;7.若计算时间大于终止时间,则结束计算,否则对更新后的包壳几何划分网格,并执行步骤5~7。本发明可分析核反应堆事故工况下包壳的复杂氧化行为,对核反应堆安全分析具有重要意义。
Description
技术领域
本发明涉及核反应堆安全分析方法领域,具体涉及一种核反应堆事故工况下耦合传热传质的包壳氧化分析方法。
背景技术
核反应堆燃料包壳作为核电站设计中的第二道屏障,采用锆合金制成,在长期运行的条件下不会使放射性裂变产物逸出,对反应堆安全具有重要作用。当反应堆发生事故例如冷却剂丧失事故(Loss-Of-Coolant Accident,LOCA)时,由于无法有效消除衰变热,包壳温度迅速上升,并发生剧烈的锆-水或锆-蒸汽氧化反应。当温度达到1500K时,反应热可与衰变热相当;当温度达到1800K时,反应热可达衰变热的10倍。因此包壳氧化现象将加速堆芯升温,过程中释放的氢气将引入爆炸风险,同时反应后包壳生成脆性的氧化物,更易产生裂纹或发生破碎,引起放射性裂变产物的逸出。因此,包壳氧化行为成为核反应堆事故分析中的重要内容,开发高效精确的包壳氧化模型对事故缓解措施的制定具有重要意义。
目前的反应堆事故分析软件广泛采用实验数据拟合的抛物线动力学关系式计算包壳氧化行为。但该方法基于包壳吸氧增重和化学计量二氧化锆假设,仅能获得等效的氧化层厚度,无法分析包壳内的氧原子浓度分布、物相演变过程及长时间、蒸汽不足等工况下的复杂行为。
发明内容
为克服上述现有技术存在的问题,本发明的目的在于提供一种核反应堆事故工况下耦合传热传质的包壳氧化分析方法,可以计算包壳内部的氧原子扩散、热传导和膨胀变形,获得核反应堆事故瞬态下的复杂氧化行为。
为达到上述目的,本发明采用如下技术方案:
一种核反应堆事故工况下耦合传热传质的包壳氧化分析方法,包括如下步骤:
步骤1:建立包壳的几何结构:环形包壳由三个相构成,从内侧到外侧分别为β-Zr,α-Zr(O)和ZrO2;包壳外侧与水或水蒸汽接触,包壳内侧与包壳和芯块间的气隙接触;氧原子由外侧进入包壳;
步骤2:建立计算模型的控制方程。具体包括如下内容:
1)氧原子扩散方程:
式中:FO为氧原子通量密度,v为变形速度,CO为氧原子浓度,D为氧原子扩散系数,r为空间坐标,t为时间;
2)锆原子守恒方程:
式中:FZr为锆原子通量密度,CZr为锆原子浓度;
3)导热方程:
式中:Q为热流密度,ρ为密度,Cp为比热容,T为温度,λ为导热系数;
步骤3:建立相界面迁移方程,具体包括如下内容:
1)在β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,分别建立如下方程:
式中:p和p-1分别为相邻两相的编号,CO,p为相界面处p相的氧原子浓度,CO,p-1为相界面处p-1相的氧原子浓度,ξ为相界面位置,vp为相界面处p相的变形速度,vp-1为相界面处p-1相的变形速度,Dp为相界面处p相的氧原子扩散系数,Dp-1为相界面处p-1相的氧原子扩散系数,CZr,p为相界面处p相的锆原子浓度,CZr,p-1为相界面处p-1相的锆原子浓度;
2)在气隙与β-Zr的界面和ZrO2与水蒸汽的界面,分别建立如下方程:
式中:ξ1为气隙与β-Zr的界面位置,ξ4为ZrO2与水蒸汽的界面位置,v1为气隙与β-Zr界面处的变形速度,v3为ZrO2与水蒸汽界面处的变形速度;
步骤4:设置相界面处的边界条件,具体包括如下内容:
1)建立氧原子扩散方程的边界条件:
对于β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,在恒温条件下建立的边界条件为:式中:为相界面处p-1相的平衡氧原子浓度,为相界面处p相的平衡氧原子浓度;对于温度瞬变条件,边界条件修正为:式中:τ为弛豫时间;
对于ZrO2与水蒸汽的界面,若蒸汽不足,建立的边界条件为:式中:FO,0为已知的氧原子通量密度,可随时间变化;若蒸汽充足,在恒温条件下建立的边界条件为:式中:为ZrO2与水蒸汽界面处ZrO2的平衡氧原子浓度;对于温度瞬变条件,边界条件修正为:
2)建立导热方程的边界条件:
对于β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,建立的边界条件为:Tp-1=Tp,式中:λp-1为相界面处p-1相的导热系数,λp为相界面处p相的导热系数,Tp-1为相界面处p-1相的温度,Tp为相界面处p相的温度,为相界面处p-1相的温度梯度,为相界面处p相的温度梯度,Qch为相界面处的化学反应热流密度;
对于气隙与β-Zr的界面和ZrO2与水蒸汽的界面,根据不同的模拟条件有三种可能的边界条件,分别为:①第一类边界条件:T=T0,式中:T0为已知的温度,可随时间变化;②第二类边界条件:式中:Q0为已知的热流密度,可随时间变化;③第三类边界条件:式中:hG为对流换热系数,TG为包壳外部流体主流温度;
步骤5:离散步骤2建立的控制方程,基于步骤4建立的边界条件求解获得氧原子浓度场、锆原子浓度场、温度场、变形速度场;具体包括以下子步骤:
1)对初始时刻的包壳几何结构划分网格;
2)采用控制容积积分法分别离散氧原子扩散方程、锆原子守恒方程和导热方程,得到:
式中:aO,P1、aO,W、aO,E、aO,P2、aZr,w、aZr,e、aZr,b、aT,P1、aT,W、aT,E、aT,P2为离散方程中的系数,分别为t+Δt时刻离散节点P、W、E处的氧原子浓度,为t时刻离散节点P处的氧原子浓度,分别为t+Δt时刻网格界面w、e处的变形速度,分别为t+Δt时刻离散节点P、W、E处的温度,为t时刻离散节点P处的温度;
3)采用三对角阵算法分别求解导热方程和氧原子扩散方程的离散方程,获得t+Δt时刻的温度场和氧原子浓度场;
4)根据温度场和氧原子浓度场更新平衡氧原子浓度、氧原子扩散系数、导热系数、密度、比热容等物性;
5)氧原子浓度和锆原子浓度的数量关系为:COMO+CZrMZr=ρ,式中MO和MZr分别为氧原子和锆原子的摩尔质量;根据氧原子浓度和锆原子浓度的数量关系,计算得到t+Δt时刻的锆原子浓度场;
6)根据锆原子浓度场,计算锆原子守恒方程的离散方程中的系数aZr,w、aZr,e、aZr,b;联立离散方程组,求解t+Δt时刻的变形速度场;
7)若本次迭代计算得到的温度场、氧原子浓度场、变形速度场与上一次迭代计算结果间的相对误差小于设定值,则执行步骤6,否则再次执行子步骤2)~7);
步骤6:基于步骤3建立的相界面迁移方程,计算各个相界面的迁移速度。采用梯形求积公式计算Δt时间内各个相界面的位移,获得t+Δt时刻各个相界面的位置,更新包壳的几何结构;
步骤7:若此时的计算时间大于设定的终止时间,则结束计算,否则对更新后的包壳几何结构划分网格,执行步骤5~7。
与现有技术相比,本发明具备如下优点:
1.建立了氧原子扩散方程、锆原子守恒方程和导热方程,可以获得包壳氧化过程中的氧原子浓度场、锆原子浓度场、温度场;
2.可以计算包壳氧化时由于温度变化、氧原子浓度变化、相变等引起的变形;
3.可以获得β-Zr,α-Zr(O)和ZrO2的厚度,进行精细的物相演化分析;
4.可以分析蒸汽不足等复杂事故工况下的包壳氧化行为。
附图说明
图1为核反应堆事故工况下耦合传热传质的包壳氧化分析方法的计算流程图。
图2为核反应堆事故工况下耦合传热传质的包壳氧化分析方法建立的包壳几何模型。
图3为核反应堆事故工况下耦合传热传质的包壳氧化分析中离散和求解控制方程的计算流程图。
具体实施方式
下面结合附图和具体实施方式对本发明作详细的说明:
如图1所示,本发明一种核反应堆事故工况下耦合传热传质的包壳氧化分析方法,包括如下步骤:
步骤1:建立如图2所示的包壳几何结构。环形包壳由三个相构成,从内侧到外侧分别为β-Zr,α-Zr(O)和ZrO2。包壳外侧与水或水蒸汽接触,包壳内侧与包壳和芯块间的气隙接触。氧原子由外侧进入包壳。包壳几何结构中不考虑α-Zr(O)和ZrO2的产生过程,因此对于未发生氧化的包壳,将α-Zr(O)和ZrO2的厚度设置为相对于β-Zr厚度可忽略的值。
步骤2:忽略轴向和周向上的原子传递和热传导,不考虑晶界扩散的影响,建立计算模型的控制方程。具体包括如下内容:
1)氧原子扩散方程:
式中:FO为氧原子通量密度,v为变形速度,CO为氧原子浓度,D为氧原子扩散系数,r为空间坐标,t为时间。
2)锆原子守恒方程:
式中:FZr为锆原子通量密度,CZr为锆原子浓度。
3)导热方程:
式中:Q为热流密度,ρ为密度,Cp为比热容,T为温度,λ为导热系数。
步骤3:根据相界面处的原子守恒规律,分别建立相界面迁移方程。具体包括如下内容:
1)在β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,建立如下方程:
式中:p和p-1分别为相邻两相的编号,CO,p为相界面处p相的氧原子浓度,CO,p-1为相界面处p-1相的氧原子浓度,ξ为相界面位置,vp为相界面处p相的变形速度,vp-1为相界面处p-1相的变形速度,Dp为相界面处p相的氧原子扩散系数,Dp-1为相界面处p-1相的氧原子扩散系数,CZr,p为相界面处p相的锆原子浓度,CZr,p-1为相界面处p-1相的锆原子浓度。
2)在气隙与β-Zr的界面和ZrO2与水蒸汽的界面,分别建立如下方程:
式中:ξ1为气隙与β-Zr的界面位置,ξ4为ZrO2与水蒸汽的界面位置,v1为气隙与β-Zr界面处的变形速度,v3为ZrO2与水蒸汽界面处的变形速度。
步骤4:设置相界面处的边界条件。具体包括如下内容:
1)建立氧原子扩散方程的边界条件:
对于β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,在恒温条件下,假设相界面处于热力学平衡状态,建立的边界条件为:式中:为相界面处p-1相的平衡氧原子浓度,为相界面处p相的平衡氧原子浓度。对于温度瞬变条件,边界条件修正为:式中:τ为弛豫时间。
对于ZrO2与水蒸汽的界面,若蒸汽不足,建立的边界条件为:式中:FO,0为已知的氧原子通量密度,可随时间变化。若蒸汽充足,在恒温条件下,假设相界面处于热力学平衡状态,建立的边界条件为:式中:为ZrO2与水蒸汽界面处ZrO2的平衡氧原子浓度;对于温度瞬变条件,边界条件修正为:
2)建立导热方程的边界条件:
对于β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,建立的边界条件为:Tp-1=Tp,式中:λp-1为相界面处p-1相的导热系数,λp为相界面处p相的导热系数,Tp-1为相界面处p-1相的温度,Tp为相界面处p相的温度,为相界面处p-1相的温度梯度,为相界面处p相的温度梯度,Qch为相界面处的化学反应热流密度。
对于气隙与β-Zr的界面和ZrO2与水蒸汽的界面,根据不同的模拟条件有三种可能的边界条件,分别为:①第一类边界条件:T=T0,式中:T0为已知的温度,可随时间变化;②第二类边界条件:式中:Q0为已知的热流密度,可随时间变化;③第三类边界条件:式中:hG为对流换热系数,TG为包壳外部流体主流温度。
步骤5:离散步骤2建立的控制方程,基于步骤4建立的边界条件求解获得氧原子浓度场、锆原子浓度场、温度场、变形速度场。如图3所示,具体包括以下子步骤:
1)对初始时刻的包壳几何结构划分网格。
2)采用控制容积积分法分别离散氧原子扩散方程、锆原子守恒方程和导热方程,得到:
式中:aO,P1、aO,W、aO,E、aO,P2、aZr,w、aZr,e、aZr,b、aT,P1、aT,W、aT,E、aT,P2为离散方程中的系数,分别为t+Δt时刻离散节点P、W、E处的氧原子浓度,为t时刻离散节点P处的氧原子浓度,分别为t+Δt时刻网格界面w、e处的变形速度,分别为t+Δt时刻离散节点P、W、E处的温度,为t时刻离散节点P处的温度。
3)采用三对角阵算法分别求解导热方程和氧原子扩散方程的离散方程,获得t+Δt时刻的温度场和氧原子浓度场。
4)根据温度场和氧原子浓度场更新平衡氧原子浓度、氧原子扩散系数、导热系数、密度、比热容等物性。
5)氧原子浓度和锆原子浓度的数量关系为:COMO+CZrMZr=ρ,式中MO和MZr分别为氧原子和锆原子的摩尔质量。根据氧原子浓度和锆原子浓度的数量关系,计算得到t+Δt时刻的锆原子浓度场。
6)根据锆原子浓度场,计算锆原子守恒方程的离散方程中的系数aZr,w、aZr,e、aZr,b。联立离散方程组,求解t+Δt时刻的变形速度场。
7)若本次迭代计算得到的温度场、氧原子浓度场、变形速度场与上一次迭代计算结果间的相对误差小于设定值,则执行步骤6,否则再次执行子步骤2)~7)。
步骤6:基于步骤3建立的相界面迁移方程,计算各个相界面的迁移速度。采用梯形求积公式计算Δt时间内各个相界面的位移,获得t+Δt时刻各个相界面的位置,更新包壳的几何结构。
步骤7:若此时的计算时间大于设定的终止时间,则结束计算,否则对更新后的包壳几何结构划分网格,执行步骤5~7。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施方式仅限于此,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单的推演或替换,都应当视为属于本发明由所提交的权利要求书确定专利保护范围。
Claims (1)
1.一种核反应堆事故工况下耦合传热传质的包壳氧化分析方法,其特征在于:包括如下步骤:
步骤1:建立包壳的几何结构:环形包壳由三个相构成,从内侧到外侧分别为β-Zr,α-Zr(O)和ZrO2;包壳外侧与水或水蒸汽接触,包壳内侧与包壳和芯块间的气隙接触;氧原子由外侧进入包壳;
步骤2:建立计算模型的控制方程,具体包括如下内容:
1)氧原子扩散方程:
式中:FO为氧原子通量密度,v为变形速度,CO为氧原子浓度,D为氧原子扩散系数,r为空间坐标,t为时间;
2)锆原子守恒方程:
式中:FZr为锆原子通量密度,CZr为锆原子浓度;
3)导热方程:
式中:Q为热流密度,ρ为密度,Cp为比热容,T为温度,λ为导热系数;
步骤3:建立相界面迁移方程,具体包括如下内容:
1)在β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,分别建立如下方程:
式中:p和p-1分别为相邻两相的编号,CO,p为相界面处p相的氧原子浓度,CO,p-1为相界面处p-1相的氧原子浓度,ξ为相界面位置,vp为相界面处p相的变形速度,vp-1为相界面处p-1相的变形速度,Dp为相界面处p相的氧原子扩散系数,Dp-1为相界面处p-1相的氧原子扩散系数,CZr,p为相界面处p相的锆原子浓度,CZr,p-1为相界面处p-1相的锆原子浓度;
2)在气隙与β-Zr的界面和ZrO2与水蒸汽的界面,分别建立如下方程:
式中:ξ1为气隙与β-Zr的界面位置,ξ4为ZrO2与水蒸汽的界面位置,v1为气隙与β-Zr界面处的变形速度,v3为ZrO2与水蒸汽界面处的变形速度;
步骤4:设置相界面处的边界条件,具体包括如下内容:
1)建立氧原子扩散方程的边界条件:
对于β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,在恒温条件下建立的边界条件为:式中:为相界面处p-1相的平衡氧原子浓度,为相界面处p相的平衡氧原子浓度;对于温度瞬变条件,边界条件修正为:式中:τ为弛豫时间;
对于ZrO2与水蒸汽的界面,若蒸汽不足,建立的边界条件为:式中:FO,0为已知的氧原子通量密度,随时间变化;若蒸汽充足,在恒温条件下建立的边界条件为:式中:为ZrO2与水蒸汽界面处ZrO2的平衡氧原子浓度;对于温度瞬变条件,边界条件修正为:
2)建立导热方程的边界条件:
对于β-Zr与α-Zr(O)的界面和α-Zr(O)与ZrO2的界面,建立的边界条件为:Tp-1=Tp,式中:λp-1为相界面处p-1相的导热系数,λp为相界面处p相的导热系数,Tp-1为相界面处p-1相的温度,Tp为相界面处p相的温度,为相界面处p-1相的温度梯度,为相界面处p相的温度梯度,Qch为相界面处的化学反应热流密度;
对于气隙与β-Zr的界面和ZrO2与水蒸汽的界面,根据不同的模拟条件有三种可能的边界条件,分别为:①第一类边界条件:T=T0,式中:T0为已知的温度,随时间变化;②第二类边界条件:式中:Q0为已知的热流密度,随时间变化;③第三类边界条件:式中:hG为对流换热系数,TG为包壳外部流体主流温度;
步骤5:离散步骤2建立的控制方程,基于步骤4建立的边界条件求解获得氧原子浓度场、锆原子浓度场、温度场、变形速度场;具体包括以下子步骤:
1)对初始时刻的包壳几何结构划分网格;
2)采用控制容积积分法分别离散氧原子扩散方程、锆原子守恒方程和导热方程,得到:
式中:aO,P1、aO,W、aO,E、aO,P2、aZr,w、aZr,e、aZr,b、aT,P1、aT,W、aT,E、aT,P2为离散方程中的系数,分别为t+Δt时刻离散节点P、W、E处的氧原子浓度,为t时刻离散节点P处的氧原子浓度,分别为t+Δt时刻网格界面w、e处的变形速度,分别为t+Δt时刻离散节点P、W、E处的温度,为t时刻离散节点P处的温度;
3)采用三对角阵算法分别求解导热方程和氧原子扩散方程的离散方程,获得t+Δt时刻的温度场和氧原子浓度场;
4)根据温度场和氧原子浓度场更新平衡氧原子浓度、氧原子扩散系数、导热系数、密度和比热容;
5)氧原子浓度和锆原子浓度的数量关系为:COMO+CZrMZr=ρ,式中MO和MZr分别为氧原子和锆原子的摩尔质量;根据氧原子浓度和锆原子浓度的数量关系,计算得到t+Δt时刻的锆原子浓度场;
6)根据锆原子浓度场,计算锆原子守恒方程的离散方程中的系数aZr,w、aZr,e、aZr,b;联立离散方程组,求解t+Δt时刻的变形速度场;
7)若本次迭代计算得到的温度场、氧原子浓度场、变形速度场与上一次迭代计算结果间的相对误差小于设定值,则执行步骤6,否则再次执行子步骤2)~7);
步骤6:基于步骤3建立的相界面迁移方程,计算各个相界面的迁移速度;采用梯形求积公式计算Δt时间内各个相界面的位移,获得t+Δt时刻各个相界面的位置,更新包壳的几何结构;
步骤7:若此时的计算时间大于设定的终止时间,则结束计算,否则对更新后的包壳几何结构划分网格,执行步骤5~7。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111400921.3A CN114093432B (zh) | 2021-11-19 | 2021-11-19 | 核反应堆事故工况下耦合传热传质的包壳氧化分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111400921.3A CN114093432B (zh) | 2021-11-19 | 2021-11-19 | 核反应堆事故工况下耦合传热传质的包壳氧化分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114093432A true CN114093432A (zh) | 2022-02-25 |
CN114093432B CN114093432B (zh) | 2023-03-21 |
Family
ID=80303738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111400921.3A Active CN114093432B (zh) | 2021-11-19 | 2021-11-19 | 核反应堆事故工况下耦合传热传质的包壳氧化分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114093432B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004020411A (ja) * | 2002-06-18 | 2004-01-22 | Toshiba Corp | 原子力発電プラントおよびその運転方法 |
US20140211903A1 (en) * | 2013-01-28 | 2014-07-31 | Areva Np Inc. | Method of operating a power generator based on noble metal induced oxidation of a heat transfer surface |
CN107451398A (zh) * | 2017-07-07 | 2017-12-08 | 西安交通大学 | 压水堆核电厂严重事故分析方法 |
CN109767853A (zh) * | 2018-12-19 | 2019-05-17 | 西安交通大学 | 模拟核反应堆失水事故后包壳骤冷破裂实验装置及方法 |
CN112989651A (zh) * | 2021-02-06 | 2021-06-18 | 西安交通大学 | 反应堆堆芯多物理场耦合方法 |
CN113191066A (zh) * | 2021-04-30 | 2021-07-30 | 西安交通大学 | 基于无网格法的核反应堆燃料元件失效分析方法 |
-
2021
- 2021-11-19 CN CN202111400921.3A patent/CN114093432B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004020411A (ja) * | 2002-06-18 | 2004-01-22 | Toshiba Corp | 原子力発電プラントおよびその運転方法 |
US20140211903A1 (en) * | 2013-01-28 | 2014-07-31 | Areva Np Inc. | Method of operating a power generator based on noble metal induced oxidation of a heat transfer surface |
CN107451398A (zh) * | 2017-07-07 | 2017-12-08 | 西安交通大学 | 压水堆核电厂严重事故分析方法 |
CN109767853A (zh) * | 2018-12-19 | 2019-05-17 | 西安交通大学 | 模拟核反应堆失水事故后包壳骤冷破裂实验装置及方法 |
CN112989651A (zh) * | 2021-02-06 | 2021-06-18 | 西安交通大学 | 反应堆堆芯多物理场耦合方法 |
CN113191066A (zh) * | 2021-04-30 | 2021-07-30 | 西安交通大学 | 基于无网格法的核反应堆燃料元件失效分析方法 |
Non-Patent Citations (3)
Title |
---|
何晓强等: "基于扩散方程的包壳氧化模型", 《核动力工程》 * |
张喜燕: "UO_2-锆合金化学反应条件下的包壳水侧氧化动力学", 《核科学与工程》 * |
张晓宇: "高温低温度梯度场下UO_2芯块气孔迁移计算研究", 《科技创新导报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN114093432B (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Beyer et al. | GAPCON-THERMAL-2: a computer program for calculating the thermal behavior of an oxide fuel rod | |
Weisman et al. | Methods for detailed thermal and hydraulic analysis of water-cooled reactors | |
Fee et al. | Fuel-cladding chemical interaction in uranium-plutonium oxide fast reactor fuel pins | |
CN114091310B (zh) | 反应堆严重事故中包壳行为多尺度多物理场耦合分析方法 | |
Saito et al. | Measurement and evaluation on pulsing characteristics and experimental capability of NSRR | |
JP4854654B2 (ja) | 炉心性能計算装置 | |
CN114093432B (zh) | 核反应堆事故工况下耦合传热传质的包壳氧化分析方法 | |
Godesar et al. | Comethe II-a computer code for predicting the mechanical and thermal behavior of a fuel pin | |
CN113192566B (zh) | 反应堆严重事故下熔池瞬态相变模拟方法 | |
Notley | A computer program to predict the performance of uo2 fuel elements irradiated at high power outputs to a burnup of 10 000 mwd/mtu | |
CN118114581A (zh) | 铅铋冷却核反应堆核热力材耦合仿真方法 | |
Warner et al. | A statistical fuel swelling and fission gas release model | |
CN116644628B (zh) | 判定弥散型板燃料元件起泡与沸腾临界先后顺序的数值模拟方法 | |
CN116070424A (zh) | 一种铅铋快堆铅池热分层降阶分析方法 | |
Notley | Calculation of Fission-Product Gas Pressures in Operating UO2 Fuel Elements | |
Aliev et al. | Software package for safety analysis of dry storage regimes of spent russian nuclear fuel | |
Cox et al. | Performance analysis of a mixed-oxide LMFBR fuel pin | |
Beyer et al. | User's guide for GAPCON-THERMAL-2: A computer program for calculating the thermal behavior of an oxide fuel rod | |
CN117408182A (zh) | 一种铅铋燃料泄漏定量评估的分析方法 | |
Jha | Plenum Gas Effect on Fuel Temperature | |
Holt et al. | Simulation of nuclear fuel sheath deformation following fuel power transients | |
JP4323901B2 (ja) | 沸騰水型原子炉のボイド率 | |
Lassmann et al. | The development of fuel performance models at the European institute for transuranium elements | |
Pautz | A time-dependent neutron transport model and its coupling to thermal-hydraulics | |
CN118114600A (zh) | 一种适用于环形燃料的双通道物理热工耦合方法 |
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 |