【技术实现步骤摘要】
本专利技术涉及核反应堆多群核数据库和反应堆物理计算领域,具体涉及一种获取多群核数据库中的多群Pn散射矩阵的方法。
技术介绍
反应堆数值计算中,多群Pn散射矩阵(其中n=0,1,2,…)起着至关重要的作用。目前针对多群Pn散射矩阵(其中n=0,1,2,…),传统采用的方法是由评价核数据库(以下简称ENDF)通过一系列处理得到,此方法最大缺陷在于针对弹性散射,在中子入射能量在200eV以下的能量范围的弹性散射反应,没有考虑到靶核的热运动情况。随着反应堆计算要求的逐渐提高,传统方法所给出多群Pn散射矩阵(其中n=0,1,2,…)计算方法已不能满足要求。传统方法计算某温度下的多群Pn散射矩阵(其中n=0,1,2,…)分为以下几步:1.对于某核素,从ENDF出发,经评价核数据处理程序得到该温度下关于能量相关的连续能量中子反应截面数据库(以下简称PENDF)。从PENDF中获得温度为TK下的连续能量弹性散射截面σs,T(E)和温度为TK下的连续能量总截面σt,T(E),其中,s代表弹性散射反应,t代表总反应,T代表温度,E代表中子入射能量;再从ENDF中获得温度为0K下的弹性散射角余弦分布由f(E,μ)表示,其中,E代表入射能量,μ代表弹性散射角余弦。2.计算多群Pn散射矩阵(其中n=0,1,2,…)产生时所需要的温度为TK下的连续能量中子通量密度φn,T(E)(其中n=0,1,2,…),其中σs0,T(E′→E)(即令n=0)是基于渐进散射模型,求解0维中子慢化方程所得出,即: ...
【技术保护点】
一种获取反应堆多群数据库中的多群Pn散射矩阵的方法,其特征在于:包括如下步骤:步骤1:读取ENDF以及PENDF,获得温度为0K下的连续能量弹性散射截面σs,0(E)、温度为0K下弹性散射角余弦分布f(E,μ)和温度为TK下连续能量总反应截面σt,T(E);其中,s代表弹性散射反应,t代表总反应,0代表温度为0K,T代表温度为TK,μ代表弹性散射角余弦;步骤2:在中子入射能量在200eV以下的能量范围引入自由气体模型,以步骤1所获取的温度为0K下的连续能量弹性散射截面σs,0(E)和温度为0K下的弹性散射角余弦分布f(E,μ)为基础,结合0K下的弹性散射角余弦分布f(E,μ),随后计算出各阶基于自由气体模型的温度为TK下的连续能量Pn散射矩阵表达式σsn,T(E→E′),其中n=0,1,2,…;基于自由气体模型的温度为TK下的连续能量多群Pn散射矩阵表达式为σsn,T(E→E′)=β5/24Eexp(E/kT)∫0∞tσs,0(kTAt2)×exp(-t2/A)ψn(t ...
【技术特征摘要】
1.一种获取反应堆多群数据库中的多群Pn散射矩阵的方法,其特征在于:包括如下步骤:步骤1:读取ENDF以及PENDF,获得温度为0K下的连续能量弹性散射截面σs,0(E)、温度为0K下弹性散射角余弦分布f(E,μ)和温度为TK下连续能量总反应截面σt,T(E);其中,s代表弹性散射反应,t代表总反应,0代表温度为0K,T代表温度为TK,μ代表弹性散射角余弦;步骤2:在中子入射能量在200eV以下的能量范围引入自由气体模型,以步骤1所获取的温度为0K下的连续能量弹性散射截面σs,0(E)和温度为0K下的弹性散射角余弦分布f(E,μ)为基础,结合0K下的弹性散射角余弦分布f(E,μ),随后计算出各阶基于自由气体模型的温度为TK下的连续能量Pn散射矩阵表达式σsn,T(E→E′),其中n=0,1,2,…;基于自由气体模型的温度为TK下的连续能量多群Pn散射矩阵表达式为 σ s n , T ( E → E ′ ) = β 5 / 2 4 E exp ( E / k T ) ∫ 0 ∞ tσ s , 0 ( k T A t 2 ) × exp ( - t 2 / A ) ψ n ( t ) d t , β = ( A + 1 ) / A - - - ( 6 ) ]]>其中,E——中子入射能量E′——中子出射能量A——靶核与中子的质量比T——开氏温度k——玻尔兹曼常数 ψ n ( t ) = H ( t + - t ) H ( t - t - ) × ∫ ϵ max - t t + ϵ min exp ( - x 2 ) Q n ( x , t ) d x + H ( t - t + ) × ∫ t - ϵ min t + ϵ min exp ( - x 2 ) Q n ( x , t ) d x - - - ( 7 ) ]]> ϵ max = ( A + 1 ) max ( E , E ′ ) / k T ϵ min = ( A + 1 ) min ( E , E ′ ) / k T - - - ( 8 ) ]]> t ± = ϵ max ± ϵ min 2 - - - ( 9 ) ]]> Q n ( x , t ) = 4 π ∫ 0 2 π P n ( μ l a b ) P ( μ C M ) d φ μ C M = 1 4 x 2 ...
【专利技术属性】
技术研发人员:刘宙宇,徐嘉隆,吴宏春,祖铁军,
申请(专利权)人:西安交通大学,
类型:发明
国别省市:陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。