CN112052579A - 一种基于浮动网格的核-热-力多物理耦合计算方法 - Google Patents

一种基于浮动网格的核-热-力多物理耦合计算方法 Download PDF

Info

Publication number
CN112052579A
CN112052579A CN202010888363.9A CN202010888363A CN112052579A CN 112052579 A CN112052579 A CN 112052579A CN 202010888363 A CN202010888363 A CN 202010888363A CN 112052579 A CN112052579 A CN 112052579A
Authority
CN
China
Prior art keywords
calculation
axial
plane
temperature
neutron
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.)
Pending
Application number
CN202010888363.9A
Other languages
English (en)
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.)
Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Original Assignee
Shanghai Nuclear Engineering Research and Design Institute Co Ltd
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 Shanghai Nuclear Engineering Research and Design Institute Co Ltd filed Critical Shanghai Nuclear Engineering Research and Design Institute Co Ltd
Priority to CN202010888363.9A priority Critical patent/CN112052579A/zh
Publication of CN112052579A publication Critical patent/CN112052579A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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

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)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明公开了一种基于浮动网格的核‑热‑力多物理耦合计算方法,包括建立燃料棒几何模型,进行核‑热‑力多物理迭代计算流程,开展平面传热计算、平面力学计算和平面中子输运计算和进行轴向一维中子扩散计算。本发明基于可以变形的浮动网格耦合求解中子学、热工水力、机械力学问题,模拟燃料棒在反应堆内极端环境下的热传导、机械变形、中子输运、核素燃耗等行为,同时获取燃料棒的温度、应力、变形、棒内压力、中子群通量、功率、燃耗、同位素密度等多物理状态。

Description

一种基于浮动网格的核-热-力多物理耦合计算方法
技术领域
本发明属于核电数值模拟技术领域,特别是涉及模拟燃料棒在反应堆内多物理耦合状态的技术。
背景技术
核燃料棒在反应堆内极端条件(高温、高压、强辐照)下产生的物理现象和物理过程十分复杂。尽管如此,在现象学研究的基础上,根据对燃料棒行为、现象的研究和理解,可以把燃料棒的多物理问题归结为传热问题、力学问题、中子学问题以及化学腐蚀问题。其中,由于化学腐蚀的变化非常缓慢,数值模拟比较困难,暂不考虑。传热问题研究冷却剂和燃料棒的温度分布,力学问题研究芯块和包壳的变形和应力分布,中子学问题研究中子通量密度分布、功率分布和燃耗分布。这三个问题虽然分属不同的学科领域,但相互之间存在不可忽略的耦合效应。其耦合关系主要表现为:燃料棒的轴向功率分布影响流道内冷却剂的纵向温度分布,进而影响燃料棒内部的平面温度分布;温度场的变化影响热膨胀等变形效应,并改变应力场;温度、应力、几何位置的变化又使材料物性发生改变,反馈到物理方程中,导致中子通量密度分布、功率分布随之改变;最后,由于燃料棒的功率分布发生了变化,再度影响温度分布。这三个问题的解,也就是温度、位移、应力应变和中子通量密度,共同主导了燃料棒的基本物理状态和变化规律。它们互为因果,相互影响,共同构成了燃料棒多物理耦合问题。
在传统反应堆设计中,不会同时考虑完整的多物理耦合,而是通过一些假设,使热传导问题、力学问题和中子学问题部分脱耦然后分别计算求解。计算工具包括热工水力设计程序、燃料性能分析程序和中子物理计算程序等,它们之间不产生即时、双向的信息交换。于是,复杂的多物理耦合问题就被分解为几个较简单的独立问题,求解所需的计算量大大减少。在计算资源较贫乏的年代,这种学科之间相互独立的反应堆设计方法的确具有重要意义。但同时也注意到,单一设计程序的结果并不能完全反映堆内复杂的多物理状态,计算量的减少来源于反应堆物理模型的简化、计算精度的降低。为了保证程序的计算结果保守,最终要以牺牲核电厂的经济性为代价。随着计算机技术的发展,计算能力的提高,利用多物理耦合技术来研究反应堆中复杂现象的条件也越来越成熟。
本发明涉及一种基于浮动网格的核-热-力多物理耦合计算方法,该方法为燃料棒模型建立可以变形的浮动网格,基于该网格耦合求解中子学、热工水力、机械力学问题,模拟燃料棒在反应堆内极端环境下的热传导、机械变形、中子输运、核素燃耗等行为,同时获取燃料棒的温度、应力、变形、棒内压力、中子群通量、功率、燃耗、同位素密度等多物理状态。
发明内容
本发明的目的是为燃料棒性能分析提供一种基于浮动网格的核-热-力多物理耦合计算方法,该方法为燃料棒模型建立可以变形的浮动网格,基于该网格耦合求解中子学、热工水力、机械力学问题,模拟燃料棒在反应堆内极端环境下的热传导、机械变形、中子输运、核素燃耗等行为,同时获取燃料棒的温度、应力、变形、棒内压力、中子群通量、功率、燃耗、同位素密度等多物理状态。其方法如下:
建立燃料棒几何模型:先在燃料棒轴向上划分为若干层离散结构,其中每一层再在水平方向上分别划分为独立的二维有限元网格;
进行核-热-力多物理迭代计算流程:从轴向冷却剂温度计算开始,根据能量守恒原则,由进口冷却剂温度和燃料线功率密度计算得到所述轴向冷却剂温度分布,再通过壁面流体对流传热计算得到包壳外表面轴向温度分布,作为平面传热计算中将作为边界条件使用;
开展平面传热计算、平面力学计算和平面中子输运计算:依次对每一层轴向离散结构得到温度、位移、功率等物理量的平面分布;所述平面传热计算后更新温度分布,受温度影响的材料物性参数、反应截面、热膨胀变形会随之变化,进而反馈到后续计算中;所述平面力学计算后更新单元轴向变形,然后对整个平面作均匀化处理,得到轴向离散结构的平均轴向变形,在轴向中子扩散计算中将作为几何条件使用;所述轴向中子扩散计算后更新轴向功率分布,在下一次迭代的轴向冷却剂温度计算中将作为条件使用;所述平面中子输运计算后更新功率分布,据此求出更新后的芯块发热率,在所述平面传热计算中将作为内热源条件使用;
进行轴向一维中子扩散计算:求得功率轴向偏移,并判断轴向偏移是否收敛;若轴向偏移不收敛则返回迭代开始,重复执行以上步骤;若收敛则计算完成,即可同时获取温度、应力、变形、棒内压力、中子群通量、功率、燃耗、同位素密度等多物理状态参数计算结果。
优选的,所述平面传热计算和所述平面力学计算采用有限元方法,所述平面中子输运计算采用特征线方法。
本发明为燃料棒性能分析提供一种基于浮动网格的核-热-力多物理耦合计算方法。采用该方法,可以精确模拟燃料棒在堆内的核-热-力多物理过程。相比传统反应堆设计中不完全耦合的计算方法,该方法可以提高计算精度,并能进一步揭示反应堆内复杂的物理现象。
附图说明
图1为“2+1”维几何模型示意图;
图2为多物理耦合迭代流程图;
图3为多物理耦合计算结果示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
下面利用本发明所述方法针对某1000MWe级核电站压水堆的典型燃料棒进行核-热-力多物理耦合计算。
该反应堆的核功率为1000MWe,系统压力为15.5MPa,;冷却剂入口流量为15170kg/s,冷却剂入口温度为290℃。反应堆堆芯共有157盒燃料组件,每盒组件中共有264根燃料棒。该典型燃料棒基于正方形栅元,采用UO2芯块、Zirlo锆合金包壳,间隙填充氦气,芯块直径、包壳内外径、栅距分别为0.8192,0.8357,0.95,1.26cm,有效高度为426.72cm。
针对上述典型燃料棒,实施基于浮动网格的核-热-力多物理耦合计算的具体步骤如下:
1)首先对其进行几何建模,在轴向上将其均匀划分为10层,每层高度为42.672cm,再依次将每一层划分为独立的二维网格,如图1所示。
2)进入核-热-力多物理迭代计算流程,依次开展轴向冷却剂温度计算、平面传热计算、平面力学计算、平面中子输运计算、轴向一维中子扩散计算,当轴向偏移收敛则迭代结束,如图2所示。
3)将计算结果进行可视化处理,如图3所示。相比传统计算方法,核-热-力多物理耦合计算方法所得结果更加精细、全面,适合开展多物理过程分析。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。上述实施例或实施方式只是对本发明的举例说明,本发明也可以以其它的特定方式或其它的特定形式实施,而不偏离本发明的要旨或本质特征。因此,描述的实施方式从任何方面来看均应视为说明性而非限定性的。本发明的范围应由附加的权利要求说明,任何与权利要求的意图和范围等效的变化也应包含在本发明的范围内。

Claims (2)

1.一种基于浮动网格的核-热-力多物理耦合计算方法,其特征在于,其方法如下:
建立燃料棒几何模型:先在燃料棒轴向上划分为若干层离散结构,其中每一层再在水平方向上分别划分为独立的二维有限元网格;
进行核-热-力多物理迭代计算流程:从轴向冷却剂温度计算开始,根据能量守恒原则,由进口冷却剂温度和燃料线功率密度计算得到所述轴向冷却剂温度分布,再通过壁面流体对流传热计算得到包壳外表面轴向温度分布,作为平面传热计算中将作为边界条件使用;
开展平面传热计算、平面力学计算和平面中子输运计算:依次对每一层轴向离散结构得到温度、位移、功率等物理量的平面分布;所述平面传热计算后更新温度分布,受温度影响的材料物性参数、反应截面、热膨胀变形会随之变化,进而反馈到后续计算中;所述平面力学计算后更新单元轴向变形,然后对整个平面作均匀化处理,得到轴向离散结构的平均轴向变形,在轴向中子扩散计算中将作为几何条件使用;所述轴向中子扩散计算后更新轴向功率分布,在下一次迭代的轴向冷却剂温度计算中将作为条件使用;所述平面中子输运计算后更新功率分布,据此求出更新后的芯块发热率,在所述平面传热计算中将作为内热源条件使用;
进行轴向一维中子扩散计算:求得功率轴向偏移,并判断轴向偏移是否收敛;若轴向偏移不收敛则返回迭代开始,重复执行以上步骤;若收敛则计算完成,即可同时获取温度、应力、变形、棒内压力、中子群通量、功率、燃耗、同位素密度等多物理状态参数计算结果。
2.如权利要求1所述的一种基于浮动网格的核-热-力多物理耦合计算方法,其特征在于,所述平面传热计算和所述平面力学计算采用有限元方法,所述平面中子输运计算采用特征线方法。
CN202010888363.9A 2020-08-28 2020-08-28 一种基于浮动网格的核-热-力多物理耦合计算方法 Pending CN112052579A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010888363.9A CN112052579A (zh) 2020-08-28 2020-08-28 一种基于浮动网格的核-热-力多物理耦合计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010888363.9A CN112052579A (zh) 2020-08-28 2020-08-28 一种基于浮动网格的核-热-力多物理耦合计算方法

Publications (1)

Publication Number Publication Date
CN112052579A true CN112052579A (zh) 2020-12-08

Family

ID=73608351

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010888363.9A Pending CN112052579A (zh) 2020-08-28 2020-08-28 一种基于浮动网格的核-热-力多物理耦合计算方法

Country Status (1)

Country Link
CN (1) CN112052579A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114913936A (zh) * 2022-07-18 2022-08-16 西安交通大学 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法
CN115982956A (zh) * 2022-12-07 2023-04-18 上海交通大学 氦氙冷却移动式核反应堆确定论多物理场耦合仿真模拟方法
CN116956770A (zh) * 2023-07-03 2023-10-27 西安交通大学 一种热管反应堆堆芯多物理场耦合方法
CN117454627A (zh) * 2023-10-30 2024-01-26 上海交通大学 基于统一网格的数值反应堆核-热-材耦合模拟方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003161796A (ja) * 2001-11-27 2003-06-06 Global Nuclear Fuel-Japan Co Ltd 燃料棒熱的特性評価方法及び装置
US20110224966A1 (en) * 2008-11-10 2011-09-15 Global Nuclear Fuel - Japan Co., Ltd. System and method for evaluating nuclear reactor fueling plan
US20110264426A1 (en) * 2008-02-11 2011-10-27 Westinghouse Electric Company Llc Methodology for modeling the fuel rod power distribution within a nuclear reactor core
CN105468832A (zh) * 2015-11-19 2016-04-06 中国核动力研究设计院 一种超临界水堆堆芯核热耦合迭代方法
CN107122546A (zh) * 2017-04-27 2017-09-01 西安交通大学 一种压水堆稳态计算的多物理耦合方法
CN110598324A (zh) * 2019-09-12 2019-12-20 西安交通大学 一种核反应堆弥散型板型燃料元件堆芯流固耦合计算方法
CN111414722A (zh) * 2020-03-19 2020-07-14 西安交通大学 一种核反应堆堆芯物理与热工耦合的模拟方法
US20200265963A1 (en) * 2017-10-23 2020-08-20 Alpha Ring International, Ltd. Retrofit for fission reactor

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003161796A (ja) * 2001-11-27 2003-06-06 Global Nuclear Fuel-Japan Co Ltd 燃料棒熱的特性評価方法及び装置
US20110264426A1 (en) * 2008-02-11 2011-10-27 Westinghouse Electric Company Llc Methodology for modeling the fuel rod power distribution within a nuclear reactor core
US20110224966A1 (en) * 2008-11-10 2011-09-15 Global Nuclear Fuel - Japan Co., Ltd. System and method for evaluating nuclear reactor fueling plan
CN105468832A (zh) * 2015-11-19 2016-04-06 中国核动力研究设计院 一种超临界水堆堆芯核热耦合迭代方法
CN107122546A (zh) * 2017-04-27 2017-09-01 西安交通大学 一种压水堆稳态计算的多物理耦合方法
US20200265963A1 (en) * 2017-10-23 2020-08-20 Alpha Ring International, Ltd. Retrofit for fission reactor
CN110598324A (zh) * 2019-09-12 2019-12-20 西安交通大学 一种核反应堆弥散型板型燃料元件堆芯流固耦合计算方法
CN111414722A (zh) * 2020-03-19 2020-07-14 西安交通大学 一种核反应堆堆芯物理与热工耦合的模拟方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
费敬然 等: ""金属燃料变形多物理效应分析"", 《核动力工程》, vol. 38, no. 1, pages 67 - 71 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114913936A (zh) * 2022-07-18 2022-08-16 西安交通大学 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法
CN114913936B (zh) * 2022-07-18 2022-10-25 西安交通大学 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法
CN115982956A (zh) * 2022-12-07 2023-04-18 上海交通大学 氦氙冷却移动式核反应堆确定论多物理场耦合仿真模拟方法
CN115982956B (zh) * 2022-12-07 2023-08-22 上海交通大学 氦氙冷却移动式核反应堆确定论多物理场耦合仿真模拟方法
CN116956770A (zh) * 2023-07-03 2023-10-27 西安交通大学 一种热管反应堆堆芯多物理场耦合方法
CN116956770B (zh) * 2023-07-03 2024-02-20 西安交通大学 一种热管反应堆堆芯多物理场耦合方法
CN117454627A (zh) * 2023-10-30 2024-01-26 上海交通大学 基于统一网格的数值反应堆核-热-材耦合模拟方法
CN117454627B (zh) * 2023-10-30 2024-05-28 上海交通大学 基于统一网格的数值反应堆核-热-材耦合模拟方法

Similar Documents

Publication Publication Date Title
CN112052579A (zh) 一种基于浮动网格的核-热-力多物理耦合计算方法
CN111291494B (zh) 用于核反应堆triso燃料颗粒的多尺度多物理场耦合模拟方法
CN111414722B (zh) 一种核反应堆堆芯物理与热工耦合的模拟方法
CN110020476B (zh) 一种反应堆u型管式蒸汽发生器全三维耦合模型建立方法
CN105653869B (zh) 一种超临界水堆堆芯三维瞬态性能分析方法
CN111523233B (zh) 一种针对三维压水堆堆芯的中子输运计算方法
CN109448799B (zh) 金属冷却快堆金属燃料多物理场模型耦合方法
CN115982956B (zh) 氦氙冷却移动式核反应堆确定论多物理场耦合仿真模拟方法
CN114913936B (zh) 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法
Aufiero et al. Development of multiphysics tools for fluoride-cooled high-temperature reactors
Weng et al. Development and application of Monte Carlo and COMSOL coupling code for neutronics/thermohydraulics coupled analysis
Laureau et al. Coupled neutronics and thermal-hydraulics transient calculations based on a fission matrix approach: application to the Molten Salt Fast Reactor
Wang et al. Application of a new OpenFOAM-based neutron diffusion kinetics solver to pebble-type VHTRs
Xu et al. Neutronics/thermal-hydraulics/fuel-performance coupling for light water reactors and its application to accident tolerant fuel
Romanenko et al. Complex modeling of VVER-1000 fuel assembly using codes MCU/ATHLET
Kim et al. Preliminary coupling of MATRA code for multi-physics analysis
Gu et al. Verification of a HC-PK-CFD coupled program based a benchmark on beam trip transients for XADS reactor
de Lambert et al. Semi-analytical fluid-structure model for the analysis of fuel assembly bow in full PWR cores
Wang et al. Parametric study of effective thermal conductivity for VHTR fuel pebbles based on a neutronic and thermal coupling method
Zare et al. Introducing and validating a new method for coupling neutronic and thermal-hydraulic calculations
Lee et al. Preliminary Multi-Physics Analysis of a 2x2 Rod Array Using CUPID/GIFT Coupled Code
Girieud et al. Science Version 2: the most recent capabilities of the Framatome 3-D nuclear code package
Dong et al. The development of nuclear reactor three-dimensional neutronic thermal–hydraulic coupling code: CorTAF-2.0
Choi et al. Preliminary Multi-Physics Analysis of OPR1000 Reactor Core using coupled CUPID and nTER
CN116956770B (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
CB02 Change of applicant information

Address after: No. 29 Hong Cao Road, Xuhui District, Shanghai

Applicant after: Shanghai Nuclear Engineering Research and Design Institute Co.,Ltd.

Address before: No. 29 Hong Cao Road, Xuhui District, Shanghai

Applicant before: SHANGHAI NUCLEAR ENGINEERING RESEARCH & DESIGN INSTITUTE Co.,Ltd.

CB02 Change of applicant information
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20201208

WD01 Invention patent application deemed withdrawn after publication