本发明涉及核反应堆多群核数据库和反应堆物理计算领域，具体涉及一种获取多群核数据库中的多群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维中子慢化方程所得出，即：Σt,T(E)φ0,T(E)=∫EE/αΣs0,T(E′→E)φ0,T(E′)dE′,α=(A-1A+1)2Σt,T(E)=σt,T(E)·N;Σs0(E′→E)=σs0,T(E′→E)·Nσs0,T(E′→E)=σs,T(E)(1-α)E′---(1)]]>其中，E——中子入射能量E′——中子出射能量N——核子密度A——靶核与中子质量比φ0,T(E)——温度TK下的0阶连续能量中子通量密度σt,T(E)——温度为TK下的连续能量总截面σs,T(E)——温度为TK下的连续能量弹性散射截面高阶通量按如下公式计算：φn,T(E)=φn-1,T(E)σt,T(E)---(2)]]>其中，σt,T(E)为温度为TK下的连续能量总反应截面，φn-1,T(E)为温度为TK下的n-1阶连续能量通量，φn,T(E)为温度为TK下的n阶连续能量通量。3.由如下公式得到多群Pn散射矩阵(其中n＝0,1,2,…)：σsn,T(g→g′)=∫ΔEgdE∫ΔEg′φn,T(E)σsn,T(E→E′)dE′∫ΔEgφn,T(E)dE=∫ΔEgdE∫ΔEg′φn,T(E)σs,T(E)fn(E→E′)dE′∫ΔEgφn,T(E)dE---(3)]]>∫ΔEg′fn(E→E′)dE′=∫μ1μ2f(E,μ)Pn(μ)dμ---(4)]]>E′E=A2+2Aμ+1(A+1)2---(5)]]>其中，E——中子入射能量E′——中子出射能量g——中子入射能群g′——中子出射能群ΔEg——中子入射能群的能量间隔ΔEg′——中子出射能群的能量间隔φn,T(E)——温度TK下的n阶中子通量密度σs,T(E)——温度TK下的连续能量弹性散射截面σsn,T(E→E′)——温度为TK下的连续能量Pn散射矩阵(其中n＝0,1,2,…)σsn,T(g→g′)——温度为TK下的多群Pn阶散射矩阵μ——弹性散射角余弦μ1——对应出射能群下边界的散射角余弦μ2——对应出射能群上边界的散射角余弦f(E,μ)——中子弹性散射角余弦分布Pn(μ)——n阶勒让德多项式A——靶核与中子的质量比该方法被广泛使用，该方法仅对于连续能量弹性散射截面σs,T(E)作了关于靶核自由气体模型的处理，而在计算温度为TK下的连续能量Pn散射矩阵(其中n＝0,1,2,…)σsn,T(E→E′)时，认为靶核是处于静止的，这与实际情况不符；同时，在推导自由气体模型的过程中，会使用到0K质心系下的弹性散射角分布，但是现有的推导方法对其做出了散射角严格各向同性假设(即认为f(E,μ)＝0.5)，这与ENDF中所给出的数据是不一致的。因此，针对以上所存在的问题，需要发明一种准确、可行的多群Pn阶散射矩阵计算方法。技术实现要素：为了克服上述现有技术存在的问题，本发明的目的在于提供一种获取反应堆多群数据库中的多群Pn散射矩阵的方法，为了得到准确的多群Pn阶散射矩阵(其中n＝0,1,2,…)，本发明方法基于自由气体散射模型，通过公式推导并且在推导过程中利用了ENDF中所提供的温度为0K下的弹性散射角分布数据，得出了关于温度TK下的连续能量Pn散射矩阵表达式σsn,T(E→E′)，并基于温度TK下的连续能量Pn散射矩阵求解中子慢化方程得到温度TK下的n阶中子通量密度，进而通过并群计算计算出温度TK下的多群Pn散射矩阵。为了实现上述目的，本发明采取了以下技术方案：一种获取反应堆多群数据库中的多群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)dt,β=(A+1)/A---(6)]]>其中，E——中子入射能量E′——中子出射能量A——靶核与中子的质量比T——开氏温度k——玻尔兹曼常数ψn(t)=H(t+-t)H(t-t-)×∫ϵmax-tt+ϵminexp(-x2)Qn(x,t)dx+H(t-t+)×∫t-ϵmint+ϵminexp(-x2)Qn(x,t)dx---(7)]]>ϵmax=(A+1)max(E,E′)/kTϵmin=(A+1)min(E,E′)/kT---(8)]]>t±=ϵmax±ϵmin2---(9)]]>Qn(x,t)=4π∫02πPn(μlab)P(μCM)dφμCM=14x2t2(A+Bcosφ)μlab=14x2ϵmaxϵmin(C+Bcosφ)---(10)]]>A=(ϵmax2-x2-t2)(ϵmin2-x2-t2)C=(ϵmax2+x2-t2)(ϵmin2+x2-t2)B=[(t+x)2-ϵmax2][(t+x)2-ϵmin2]×[ϵmax2-(t-x)2][ϵmin2-(t-x)2]---(11)]]>其中，H为Heaviside阶跃函数，Pn(μlab)为n阶勒让德多项式，μlab为实验系下散射角，P(μCM)为在质心系下散射方位角为μCM的散射概率分布；严格依照ENDF中所提供的温度为0K下弹性散射角余弦分布f(E,μ)，即将并带入式中，通过数值积分求解温度为TK下的连续能量Pn散射矩阵；步骤3：使用步骤2计算出的温度为TK下的连续能量Pn散射矩阵，令n＝0得到温度为TK下的连续能量P0散射矩阵σs0,T(E→E′)，并基于此温度为TK下的续能量P0散射矩阵σs0,T(E→E′)，在中子入射能量在200eV以下的能量范围，获得温度为TK下的各阶中子通量密度φn,T(E)，其中n＝0,1,2,…。为了得到温度为TK下的各阶中子通量密度φn,T(E)，其中n＝0,1,2,…，首先求解温度TK下的0阶连续能量中子通量密度φ0,T(E)，即令φn,T(E)中的n＝0，建立温度TK下的0维中子慢化方程：σt,T(E)·N·φ0,T(E)=∫0∞σs0,T(E′→E)·N·φ0,T(E′)dE′---(12)]]>其中，E——中子入射能量E′——中子出射能量N——靶核核子密度φ0,T(E)——温度TK下的0阶连续能量中子通量密度σs0,T(E′→E)——温度为TK下的连续能量P0散射矩阵σt,T(E)——温度为TK下的连续能量总截面使用超细群方法求解该中子慢化方程，在超细群方法中，把共振能区分割成为非常精细的能量间隔，每一个这样的能量间隔称为一个超细群，认为每一个超细群宽度远小于中子与核素碰撞所获得的最大对数能降，即认为不可能发生超细群的自散射，这样只要给定最高能群的散射源后便能够依次由高能向低能逐群求解获得精细通量；在200eV以下，考虑自由气体模型，该模型会导致中子的上散射效应，在计算能谱时不能一次性由高能向低能逐群求解得到通量，需要通过迭代计算直至φ0,T(E)收敛；求得了温度TK下的0阶连续能量中子通量密度φ0,T(E)后，温度TK下的n阶连续能量中子通量密度，其中n＝1,2,…，按公式(2)计算，即：φn,T(E)=φn-1,T(E)σt,T(E)---(2)]]>其中，σt,T(E)为温度为TK下的连续能量总反应截面，φn-1,T(E)为温度为TK下的n-1阶连续能量通量，φn,T(E)为温度为TK下的n阶连续能量通量；步骤4：获得了由步骤2所计算出的各阶基于自由气体模型的连续能量Pn散射矩阵表达式σsn,T(E→E′)和步骤3所计算出的各阶中子通量密度φn,T(E)，进行并群计算；按下列并群公式计算出多群Pn散射矩阵：σsn,T(g-g′)=∫ΔEgdE∫ΔEg′φn,T(E)σsn,T(E→E′)dE′∫ΔEgφn,T(E)dE---(13)]]>其中，E——中子入射能量E′——中子出射能量g——中子入射能群g′——中子出射能群ΔEg——中子入射能群的能量间隔ΔEg′——中子出射能群的能量间隔φn,T(E)——温度TK下的n阶中子通量密度σs,T(E)——温度TK下的连续能量弹性散射截面σsn，T(E→E′)——温度为TK下的连续能量Pn散射矩阵σsn,T(g→g′)——温度为TK下的多群Pn散射矩阵。与现有技术相比，本发明有如下突出优点：1、考虑了靶核的热运动效应，能够基于自由气体模型计算出各阶连续能量Pn散射矩阵，有效的考虑了中子的共振上散射效应和共振吸收增加效应。2、基于所计算出的连续能量P0散射矩阵(即令连续能量Pn散射矩阵中的n＝0)，更加准确的描述中子与靶核碰撞后的运动并且能够得出更为真实的各阶中子通量密度，使并群计算后的多群Pn散射矩阵更加准确，并最终提高反应堆物理计算中共振计算的精度。3、结合了ENDF中真实的质心系下的弹性散射角分布，提高了自由气体模型的精度，减少了对于质心系下的弹性散射角余弦分布假设。具体实施方式下面结合具体实施方式对本发明作进一步详细说明：本发明一种获取反应堆多群数据库中的多群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散射矩阵(其中n＝0,1,2,…)表达式σsn,T(E→E′)。经相关推导，基于自由气体模型的连续能量多群Pn散射矩阵(其中n＝0,1,2,…)表达式为σsn,T(E→E′)=β5/24Eexp(E/kT)∫0∞tσs,0(kTAt2)×exp(-t2/A)ψn(t)dt,β=(A+1)/A---(6)]]>其中，E——中子入射能量E′——中子出射能量A——靶核与中子的质量比T——开氏温度k——玻尔兹曼常数ψn(t)=H(t+-t)H(t-t-)×∫ϵmax-tt+ϵminexp(-x2)Qn(x,t)dx+H(t-t+)×∫t-ϵmint+ϵminexp(-x2)Qn(x,t)dx---(7)]]>ϵmax=(A+1)max(E,E′)/kTϵmin=(A+1)min(E,E′)/kT---(8)]]>t±=ϵmax±ϵmin2---(9)]]>Qn(x,t)=4π∫02πPn(μlab)P(μCM)dφμCM=14x2t2(A+Bcosφ)μlab=14x2ϵmaxϵmin(C+Bcosφ)---(10)]]>A=(ϵmax2-x2-t2)(ϵmin2-x2-t2)C=(ϵmax2+x2-t2)(ϵmin2+x2-t2)B=[(t+x)2-ϵmax2][(t+x)2-ϵmin2]×[ϵmax2-(t-x)2][ϵmin2-(t-x)2]---(11)]]>其中，H为Heaviside阶跃函数，Pn(μlab)为n阶勒让德多项式，μlab为实验系下散射角，P(μCM)为在质心系下散射方位角为μCM的散射概率分布。本发明中，严格依照ENDF中所提供的温度为0K下弹性散射角余弦分布f(E,μ)，即将并带入式中，通过数值积分求解温度为TK下的连续能量Pn散射矩阵(其中n＝0,1,2,…)。3.使用步骤2计算出的温度为TK下的连续能量Pn散射矩阵(其中n＝0,1,2,…)，令n＝0得到温度为TK下的连续能量P0散射矩阵σs0,T(E→E′)，并基于此连温度为TK下的续能量P0散射矩阵σs0,T(E→E′)，在中子入射能量在200eV以下的能量范围，获得温度为TK下的各阶中子通量密度φn,T(E)(其中n＝0,1,2,…)。为了得到温度为TK下的各阶中子通量密度φn,T(E)(其中n＝0,1,2,…)，首先求解温度TK下的0阶连续能量中子通量密度φ0,T(E)(即令φn,T(E)中的n＝0)，建立温度TK下的0维中子慢化方程：σt,T(E)·N·φ0,T(E)=∫0∞σs0,T(E′→E)·N·φ0,T(E′)dE′---(12)]]>其中，E——中子入射能量E′——中子出射能量N——靶核核子密度φ0,T(E)——温度TK下的0阶连续能量中子通量密度σs0,T(E′→E)——温度为TK下的连续能量P0散射矩阵σt,T(E)——温度为TK下的连续能量总截面使用超细群方法求解该中子慢化方程，在超细群方法中，把共振能区分割成为非常精细的能量间隔，每一个这样的能量间隔称为一个超细群，认为每一个超细群宽度远小于中子与核素碰撞所获得的最大对数能降，即认为不可能发生超细群的自散射，这样只要给定最高能群的散射源后便可以依次由高能向低能逐群求解获得精细通量。在200eV以下，考虑自由气体模型，该模型会导致中子的上散射效应，在计算能谱时不能一次性由高能向低能逐群求解得到通量，需要通过迭代计算直至φ0,T(E)收敛。求得了温度TK下的0阶连续能量中子通量密度φ0,T(E)后，温度TK下的n阶连续能量中子通量密度(其中n＝1,2,…)按公式(2)计算，即：φn,T(E)=φn-1,T(E)σt,T(E)---(2)]]>其中，σt,T(E)为温度为TK下的连续能量总反应截面，φn-1,T(E)为温度为TK下的n-1阶连续能量通量，φn,T(E)为温度为TK下的n阶连续能量通量。4.获得了由步骤2所计算出的各阶基于自由气体模型的连续能量Pn散射矩阵(其中n＝0,1,2,…)表达式σsn,T(E→E′)和步骤3所计算出的各阶中子通量密度φn,T(E)(其中n＝0,1,2,…)，进行并群计算。按下列并群公式计算出多群Pn散射矩阵(其中n＝0,1,2,…)：σsn,T(g-g′)=∫ΔEgdE∫ΔEg′φn,T(E)σsn,T(E→E′)dE′∫ΔEgφn,T(E)dE---(13)]]>其中，E——中子入射能量E′——中子出射能量g——中子入射能群g′——中子出射能群ΔEg——中子入射能群的能量间隔ΔEg′——中子出射能群的能量间隔φn,T(E)——温度TK下的n阶中子通量密度σs,T(E)——温度TK下的连续能量弹性散射截面σsn,T(E→E′)——温度为TK下的连续能量Pn散射矩阵(其中n＝0,1,2,…)σsn,T(g→g′)——温度为TK下的多群Pn散射矩阵(其中n＝0,1,2,…)。当前第1页1 2 3