CN106096018B - 一种在多群数据库中考虑共振弹性散射的方法 - Google Patents

一种在多群数据库中考虑共振弹性散射的方法 Download PDF

Info

Publication number
CN106096018B
CN106096018B CN201610471982.1A CN201610471982A CN106096018B CN 106096018 B CN106096018 B CN 106096018B CN 201610471982 A CN201610471982 A CN 201610471982A CN 106096018 B CN106096018 B CN 106096018B
Authority
CN
China
Prior art keywords
resonance
section
openmc
rank
square
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.)
Active
Application number
CN201610471982.1A
Other languages
English (en)
Other versions
CN106096018A (zh
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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610471982.1A priority Critical patent/CN106096018B/zh
Publication of CN106096018A publication Critical patent/CN106096018A/zh
Application granted granted Critical
Publication of CN106096018B publication Critical patent/CN106096018B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90—Details of database functions independent of the retrieved data types
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2216/00—Indexing scheme relating to additional aspects of information retrieval not explicitly covered by G06F16/00 and subgroups
    • G06F2216/01—Automatic library building

Landscapes

  • Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种在多群数据库中考虑共振弹性散射效应的方法，1、对于共振核素，用蒙特卡罗程序OpenMC统计各个背景截面点下各个温度点下各个共振能群的吸收截面，散射截面和中子产生截面，整理成共振积分表；2、对于该共振核素，用OpenMC统计各个共振能群到其他能群的0阶散射矩；3、用OpenMC统计的共振积分表和0阶散射矩替换组件程序多群数据库中该核素的共振积分表和0阶散射矩；本发明方法中，用OpenMC统计产生共振积分表和0阶散射矩，使多群截面能够考虑中子的共振弹性散射效应；用OpenMC统计产生的共振积分表和0阶散射矩替换组件程序多群数据库中的共振积分表和0阶散射矩，使组件程序能够直接使用考虑共振弹性散射效应的多群截面。

Description

一种在多群数据库中考虑共振弹性散射的方法
技术领域
本发明涉及核反应堆堆芯设计和安全技术领域，具体涉及一种在多群数据库中考虑共振弹性散射的方法。
背景技术
反应堆物理设计计算一般采用确定论方法，该方法对能量变量采用多群的近似，多群数据库通常利用NJOY程序产生。而NJOY程序在两个方面忽略了共振弹性散射效应。第一个方面，NJOY制作的共振积分表忽略了共振弹性散射效应。该程序通过求解不同温度和背景截面下的中子慢化方程制作共振积分表。为了简化中子慢化方程中散射源项的计算，采用靶核静止模型，忽略了靶核的热运动和弹性散射的共振。第二个方面，NJOY制作的0阶散射矩忽略了共振弹性散射效应。0阶散射矩通过归并点截面形式的散射概率得到，计算散射概率时需要用到评价数据库4号文件中的勒让德展开系数。但4号文件提供的展开系数没有考虑靶核的热运动，即采用靶核静止模型，同样忽略了弹性散射的共振。因此这两个方面的实质都是采用靶核静止模型，不能同时考虑靶核的热运动和弹性散射的共振，这种效应被称为共振弹性散射效应。这会导致低估中子向上散射，进而导致低估吸收反应率，从而对后续的计算带来影响。
发明内容
为了克服上述现有技术存在的问题，本发明的目的在于提供一种在多群数据库中考虑共振弹性散射效应的方法，该方法采用蒙特卡罗程序OpenMC统计计算得到共振积分表和0阶散射矩，并替换掉多群数据库中相应的部分，以在多群数据库中考虑共振弹性散射效应。
为了实现上述目的，本发明采取了以下技术方案予以实施：
一种在多群数据库中考虑共振弹性散射效应的方法，该方法包括以下步骤：
步骤1：对于一个共振核素，选取一系列温度点和背景截面点，构建对应的共振核素与1H核均匀混合的问题；采用蒙特卡罗程序OpenMC统计这些问题各个共振能群的吸收截面、散射截面和中子产生截面；运行程序时，在OpenMC的输入卡片中打开共振弹性散射的开关，使OpenMC能够考虑共振弹性散射效应；统计完成时将这些截面整理成该核素的共振积分表的形式：
其中σb是背景截面；Ix,g(σb,T)是背景截面为σb温度为T时的x反应类型第g能群的共振积分；σx,g(σb,T)是背景截面为σb温度为T时x反应类型第g能群的截面；x反应类型包括吸收反应，散射反应和中子产生反应；σa,g(σb,T)是背景截面为σb温度为T时第g能群的吸收截面；
步骤2：对于该共振核素，选取一系列温度点和一个背景截面点，采用OpenMC程序统计这些温度和这个背景截面下各个共振能群到其他能群的0阶散射矩；在统计截面时，在OpenMC的输入卡片中打开共振弹性散射的开关；
步骤3：把多群数据库中的共振积分表和共振能群到其他能群的0阶散射矩替换成OpenMC统计的共振积分表和0阶散射矩。
与现有技术相比，本发明有如下突出优点：
本发明利用现有的能够考虑共振弹性散射效应的蒙特卡罗程序OpenMC，经统计计算得到共振积分表和0阶散射矩，不需要对原先用于加工多群数据库的NJOY程序做任何修改。将多群数据库中的共振积分表和0阶散射矩替换成OpenMC统计计算得到的，不需要对多群数据库做格式上的更改，使确定论程序能够直接使用。采用新的多群数据库进行计算，能够考虑共振弹性散射效应并提高计算的精度。
附图说明
图1为OpenMC程序输入卡片共振弹性散射的开关。
图2为与OpenMC计算的无限介值增值因数相比，采用原多群数据库和新的多群数据库进行多群输运计算的无限介值增值因数误差。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细说明：
本发明利用OpenMC制作多群数据库中的共振积分表和0阶散射矩。该方法具体计算流程包括以下方面：
1)对于一个共振核素A，选取4个温度点Ti(293K，600K，900K和1100K)和10个背景截面点σb,j(10barn，28barn，52barn，64barn，140barn，260barn，1000barn，3600barn， )，构建对应的40个该核素与1H核素(简单起见令为B核素)均匀混合的问题。共振核素的核子密度取为1.0，B核素的核子密度采用公式(2)计算得到
利用OpenMC统计这些温度和背景截面下各个共振能群的吸收反应率，散射反应率，中子产生反应率和通量。在统计反应率和通量时，在OpenMC的输入卡片中打开共振弹性散射的开关，即在输入卡片中增加如图1所示的内容。利用公式(3)计算得到对应的吸收截面，散射截面和中子产生截面。
其中Rx,g为反应率；φg为通量；σx,g为截面；g是能群的编号；x是反应类型，包括吸收反应，散射反应和中子产生反应。利用公式(4)计算得到各个共振能群对应的真实背景截面。
σ'b,g,j＝λAσp,A,g+NB,jσp,B,g 公式(4)
由于多群数据库中的背景截面点对于每一个共振能群都是相同的，利用平方根线性插值公式得到选取的背景截面点σb,j对应的吸收截面，散射截面和中子产生截面。然后利用公式(1)将截面整理成共振积分表的形式；
2)对于该共振核素A，选取4个温度点Ti(293K，600K，900K和1100K)和一个背景截面点28barn，采用OpenMC程序统计这些温度和该背景截面下各个共振能群到其他能群的0阶散射反应率和通量。在统计反应率和通量时，在OpenMC的输入卡片中打开共振弹性散射的开关。采用公式(5)计算得到0阶散射矩。
其中Rs0,g'→g是g’群到g群的0阶散射反应率；φg'是g’群的通量；
σs0,g'→g是g’群到g群的0阶散射矩；
3)把原多群数据库中的共振积分表和共振能群到其他能群的0阶散射矩替换成OpenMC统计计算得到的共振积分表和0阶散射矩，制成新的多群数据库。
采用原多群数据库和新的多群数据库对一系列单棒问题进行多群输运计算，得到无限介值增值因数；对于同样的问题，采用OpenMC程序进行计算，得到基准的无限增值因数。利用公式(6)计算得到多群输运计算的误差。
error＝(k∞,MG-k∞,OpenMC)×105 公式(6)
其中k∞,OpenMC是OpenMC计算的结果；k∞,MG是多群输运计算的结果。图2给出了采用不同的多群数据库得到的无限增值因数误差。可以发现采用新的多群数据库的计算误差显著小于采用原多群数据库的计算误差。

Claims (1)

1.一种在多群数据库中考虑共振弹性散射效应的方法，其特征在于：该方法包括以下步骤：
步骤1：对于一个共振核素，选取一系列温度点和背景截面点，构建对应的共振核素与1H核均匀混合的问题；采用蒙特卡罗程序OpenMC统计这些问题各个共振能群的吸收截面、散射截面和中子产生截面；运行程序时，在OpenMC的输入卡片中打开共振弹性散射的开关，使OpenMC能够考虑共振弹性散射效应；统计完成时将这些截面整理成该核素的共振积分表的形式：
其中σb是背景截面；Ix,g(σb,T)是背景截面为σb温度为T时的x反应类型第g能群的共振积分；σx,g(σb,T)是背景截面为σb温度为T时x反应类型第g能群的截面；x反应类型包括吸收反应，散射反应和中子产生反应；σa,g(σb,T)是背景截面为σb温度为T时第g能群的吸收截面；
步骤2：对于该共振核素，选取一系列温度点和一个背景截面点，采用OpenMC程序统计这些温度和这个背景截面下各个共振能群到其他能群的0阶散射矩；在统计截面时，在OpenMC的输入卡片中打开共振弹性散射的开关；
步骤3：把多群数据库中的共振积分表和共振能群到其他能群的0阶散射矩替换成OpenMC统计的共振积分表和0阶散射矩。
CN201610471982.1A 2016-06-24 2016-06-24 一种在多群数据库中考虑共振弹性散射的方法 Active CN106096018B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610471982.1A CN106096018B (zh) 2016-06-24 2016-06-24 一种在多群数据库中考虑共振弹性散射的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610471982.1A CN106096018B (zh) 2016-06-24 2016-06-24 一种在多群数据库中考虑共振弹性散射的方法

Publications (2)

Publication Number Publication Date
CN106096018A CN106096018A (zh) 2016-11-09
CN106096018B true CN106096018B (zh) 2019-04-16

Family

ID=57252625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610471982.1A Active CN106096018B (zh) 2016-06-24 2016-06-24 一种在多群数据库中考虑共振弹性散射的方法

Country Status (1)

Country Link
CN (1) CN106096018B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013093168A1 (es) * 2011-12-21 2013-06-27 Nanobiomatters Research & Development, S. L. Materiales nanocompuestos activos basados en sales generadoras de so2 y edta y el procedimiento para su obtención
CN104614715A (zh) * 2015-03-05 2015-05-13 北京航空航天大学 一种可用于目标双站雷达散射截面测量定标与极化校准装置及其测量校准方法
CN105426342A (zh) * 2015-10-30 2016-03-23 西安交通大学 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013093168A1 (es) * 2011-12-21 2013-06-27 Nanobiomatters Research & Development, S. L. Materiales nanocompuestos activos basados en sales generadoras de so2 y edta y el procedimiento para su obtención
CN104614715A (zh) * 2015-03-05 2015-05-13 北京航空航天大学 一种可用于目标双站雷达散射截面测量定标与极化校准装置及其测量校准方法
CN105426342A (zh) * 2015-10-30 2016-03-23 西安交通大学 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于WIMS格式多群核数据制作过程的优化分析;谢明亮;《核技术》;20151231;第38卷(第11期);第110502页
确定论方法中共振弹性散射的修正方法研究;贺清明;《核动力工程》;20141231;第35卷(第S2期);第123-126页

Also Published As

Publication number Publication date
CN106096018A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
Koning et al. TALYS: modeling of nuclear reactions
Koning Bayesian Monte Carlo method for nuclear data evaluation
Belgya et al. Handbook for calculations of nuclear reaction data, RIPL-2
Cristallo et al. Effects of nuclear cross sections on 19F nucleosynthesis at low metallicities
CN108717479A (zh) 一种连续能量的中子动力学加权蒙特卡罗计算方法
CN106355017B (zh) 一种基于多群虚截面的连续能点蒙特卡罗粒子输运方法
Sprouse et al. Following nuclei through nucleosynthesis: A novel tracing technique
Voinov et al. Test of nuclear level density inputs for Hauser-Feshbach model calculations
Avrigeanu et al. Analysis of uncertainties in α-particle optical-potential assessment below the Coulomb barrier
Ortiz et al. Recent radiation damage studies and developments of the Marlowe code
CN106096018B (zh) 一种在多群数据库中考虑共振弹性散射的方法
Kavun et al. Effects of some level density models and γ-ray strength functions on production cross-section calculations of 16, 18O and 24, 26Mg radioisotopes
CN113609099A (zh) 一种基于蒙特卡罗方法制作聚变堆多群屏蔽数据库的方法
CN107092028A (zh) 一种密闭容器核材料质量的计算方法
Hall Monte Carlo perturbation theory in neutron transport calculations
Tantillo et al. Adjoint neutron flux estimator implementation and verification in the continuous energy Monte Carlo code MONK
Irfan et al. Determination of cluster size in particle-nucleus interactions at 50 and 400 Gev
Unger et al. Effective parameters and the renormalization group in grand unified theories
Alton et al. Search for light-to-heavy quark flavor changing neutral currents in ν μ N and ν¯ μ N scattering at the Fermilab Tevatron
Saha A study on the metrics of neutron induced primary radiation damage in structural materials using the recent basic evaluated nuclear data libraries
Karaca A Network Account of Models in High Energy Physics Experiments
Margaroli et al. Latest top results from CDF
Williams Uncertainty Quantification Framework for Design Optimization
Hayato Brief introduction of the neutrino event generators
MARGAROLI University “Alma Mater Studiorum.” Bologna, Via B. Pichat 6/2, Bologna, 40.127, Italy E-mail: margaroli (0bo. infn. it

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