A method of generating online feature line, which comprises the following steps; characteristic line angle and radial width of 1, entered by the user, determining the characteristic line scanning starting point position; 2, to determine the current characteristics of the line scan position; 3, to obtain the feature line information; 4, MOC computing storage network standard of coarse flux; 5 and the nearest surface distance calculation; 6, determine the lattice characteristic line will arrive; the method of the invention does not need to repeat the unit feature line at the same time generating process of scanning in the feature line, namely, side side scanning. There is no limit to the geometry of the core, which greatly broadens the scope of the traditional method of feature lines.
【技术实现步骤摘要】
一种在线生成特征线的方法
本专利技术涉及核反应堆堆芯数值计算领域,具体涉及一种在线生成特征线并将特征线信息应用于堆芯特征线计算的方法。
技术介绍
中子输运方程又称玻尔兹曼输运方程,描述的是处于运动状态的中子的守恒关系,如式(1)所示:式中:ψ——中子角通量密度;r——中子的空间位置向量;Ω——中子运动的角度向量;E——中子的能量;t——时间;v——中子运动的速度;——数学算符,其中i,j,k分别为x,y,z方向的方向向量;Σt——反应总截面;Q——中子源项;当讨论稳态问题时,认为中子通量密度不随时间变化,即:且ψ(r,Ω,E,t)=ψ(r,Ω,E),此时方程变成:在特征线上,式(2)可写成:式中:s——特征线坐标;对能量进行离散可得特征线上的多群中子输运方程:将求解区域划分为许多源区i,用离散后的角度m代替Ω。在源区内式(4)可以写成:式中:k——方向为m穿过源区的特征线编号;对源项采取平源近似,有:Qi,m(s)=Qi,m(6)联立式(5)和式(6),可以解得特征线上的中子角通量密度分布:在实际特征线计算中,令s=sk即可得到第k条特征线的出射中子角通量:由式( ...
【技术保护点】
一种在线生成特征线的方法,其特征在于:包括如下步骤;步骤1:由用户输入的特征线辐角
【技术特征摘要】
1.一种在线生成特征线的方法,其特征在于:包括如下步骤;步骤1:由用户输入的特征线辐角和径向宽度由式(13)确定特征线扫描起点的位置;式中:n——特征线扫描的起点编号dn——特征线扫描的第n个起点距离堆芯底面的距离——特征线宽度步骤2:确定目前特征线扫描进行的位置,步骤如下:1)确定组件:存储最靠近每个组件边界的所有栅元编号,当特征线扫描进入到这些栅元,则判定当前处于相应的组件;2)确定栅元:栅元的边界由若干个面组成,假设其方程分别为:f1(x,y)=0(14)f2(x,y)=0(15)……fm(x,y)=0(16)式中:m——组成栅元边界面的数目f1(x,y)——第1个面的方程f2(x,y)——第2个面的方程fm(x,y)——第m个面的方程由于一个面能够将整个空间划分为为两个半平面,而每个半平面由一个不等式表示,这样一来,栅元内部区域就由若干个不等式的交运算确定;3)确定源区:与栅元的构成类似,每个源区也是有若干半平面的交集构成,每个半平面有一个不等式确定,故通过若干个不等式的交运算确定源区;步骤3:获取特征线信息:通过定位到的源区ID,获取相应的材料区编号;步骤4:MOC计算存储粗网标通量:在每个源区内特征线方程为:式中:s——特征线上的位置坐标i——源区编号m——角度编号k——角度为m的特征线编号ψi,m,k(s)——源区i内,角度为m的第k条特征线上的角通量密度函数Σt,i——源区i内的材料总截面Qi,m(s)——源区i内,角度为m的中子源强关于特征线的函数对源项应用平源近似,公式为:Qi,m(s)=Qi,m(18)式中:Qi,m——不是s的函数,表示其值不随特征线变化联立式(17)和式(18)解得:将s用特征线长度si,m,k代入,解得的ψi,m,k(si,m,k)即为出射角通量:
【专利技术属性】
技术研发人员:刘宙宇,刘阳,曹良志,吴宏春,
申请(专利权)人:西安交通大学,
类型:发明
国别省市:陕西,61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。