基于对称群论的输运方程响应矩阵分块对角化方法技术

基于对称群论的输运方程响应矩阵分块对角化方法，包括以下步骤：1、对于整个堆芯求解区域按照规则的几何进行划分，得到规则的节块，基于变分节块方法，对中子输运方程建立泛函，构造输运方程的弱解形式；采用球谐函数以及正交多项式对中子角通量密度进行离散，球谐函数与空间正交多项式共同构成了通量的基函数；2、根据节块的几何形状，确定节块所有的对称性变换构成的对称群；利用对称群的投影算符对基函数进行投影，从而构造出不可约对称化的基函数；3、将不可约对称化的基函数替换步骤1中的球谐函数以及空间正交多项式，采用新的基函数对变量进行离散，新的响应矩阵为对角分块的矩阵；采用本发明专利技术方法大大降低矩阵的浮点运算量，存储响应矩阵所需的计算内存也大大降低。

【技术实现步骤摘要】
基于对称群论的输运方程响应矩阵分块对角化方法
本专利技术涉及核反应堆堆芯设计和安全
，具体涉及一种基于对称群论的输运方程响应矩阵分块对角化方法。
技术介绍
全堆芯输运节块计算在中子学的计算中具有重要的意义。在中子输运方程的求解中，若能够提前构造节块内部边界处中子出入射流的响应关系，则可根据该响应关系，通过节块的扫描迭代直至收敛，即可获得节块的中子泄漏，从而求解出各个节块中的中子通量分布。这种方法称为响应矩阵方法。输运方程的角度项常常采用球谐函数进行离散。若角度项展开取前N阶的球谐函数，并将输运方程在前N阶的球谐函数上进行投影，即可获得(N+1)2个线性代数方程组，即PN近似方程组。由于节块内部未知量的总数与球谐函数的个数是正相关的，因此，未知量的总数随着角度展开阶数呈平方量级增长。响应矩阵的构建通常涉及矩阵的乘法以及求逆运算，这些浮点运算量大致与未知量个数的三次方呈正比关系。因此，当未知量的个数较多时，计算内存以及浮点运算总量将会大大增加。在实际堆芯计算中，对于泄漏较强、具有强烈各向异性散射效应的小型快中子反应堆，若想得到更为精确的数值计算结果，通常要采用较高的角度阶数进行逼近本文档来自技高网...

【技术保护点】
基于对称群论的输运方程响应矩阵分块对角化方法，其特征在于：包括以下步骤：步骤一，对于整个堆芯求解区域按照规则的几何进行划分，从而得到规则的节块，基于变分节块方法，对中子输运方程建立泛函，构造输运方程的弱解形式；采用球谐函数以及正交多项式对中子角通量密度进行离散，球谐函数与空间正交多项式共同构成了通量的基函数，如公式(8)所示

【技术特征摘要】
1.基于对称群论的输运方程响应矩阵分块对角化方法，其特征在于：包括以下步骤：步骤一，对于整个堆芯求解区域按照规则的几何进行划分，从而得到规则的节块，基于变分节块方法，对中子输运方程建立泛函，构造输运方程的弱解形式；采用球谐函数以及正交多项式对中子角通量密度进行离散，球谐函数与空间正交多项式共同构成了通量的基函数，如公式(8)所示其中ψ(r,Ω)为偶通量密度；J(r,Ω)为节块边界处的拉格朗日乘子项；g(Ω)＝[Y0,0,Y2,-2,Y2,-1,Y2,0,Y2,1,Y2,2...]T为偶阶的归一化的球谐函数；f(r)＝[f1,f2,f3...]T为定义在节块v内部的正交多项式；h(r)＝[h1,h2...]T为定义在节块v边界处的正交多项式；表示位于1号边界面处的边界基函数，表示位于2号边界面处的边界基函数；以及ξ为分别为偶中子角通量密度展开矩以及中子流密度展开矩；T为转置符号；将离散形式的通量以及拉格朗日乘子项带入泛函中整理可得响应矩阵方程：j+＝Bs+Rj-(16)其中为偶中子角通量密度展开矩；j+为出射偏中子流密度展开矩；j-为入射偏中子流密度展开矩；s为源项展开矩；R,H,C,B为变分节块方法的四个响应矩阵；步骤二，根据步骤一中的节块的几何形状，确定节块所有的对称性变换构成的对称群；利用对称群的投影算符对步骤一中的基函数进行投影，从而构造出不可约对称化的基函数；投影算符由对称算符组成式中n——群的阶数，即有限群中的群元个数；sα——第α个不可约表示矩阵的维数；R——群元素；PR——与R群元对应...

【专利技术属性】
技术研发人员：李云召，李志鹏，吴宏春，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西,61