CN108734181B - 一种加速在线生成核反应堆堆芯特征线的方法 - Google Patents
一种加速在线生成核反应堆堆芯特征线的方法 Download PDFInfo
- Publication number
- CN108734181B CN108734181B CN201810499010.2A CN201810499010A CN108734181B CN 108734181 B CN108734181 B CN 108734181B CN 201810499010 A CN201810499010 A CN 201810499010A CN 108734181 B CN108734181 B CN 108734181B
- Authority
- CN
- China
- Prior art keywords
- boundaries
- characteristic curve
- boundary
- point
- line
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种加速在线生成核反应堆堆芯特征线的方法,通过选定典型几何单元,将其每个几何边界划分为若干区间,使得每个区间内对应一个目标边界,存储区间和目标边界的对应关系。使得在生成特征线段时不必全局检索所有的几何边界以确定目标边界,只需根据当前特征线的位置,判断其属于的边界号码以及区间号码,即可从内存中检索当前区间所对应的目标边界。对于实际的反应堆几何布置,可以以较少的存储代价显著加速特征线扫描过程;本发明方法使得生成特征线这一过程的计算代价节省了成千上万倍。
Description
技术领域
本发明涉及核反应堆堆芯数值计算领域,具体涉及一种加速在线 生成核反应堆堆芯特征线的方法。
背景技术
中子输运方程又称玻尔兹曼输运方程,描述的是处于运动状态的 中子的守恒关系,如式(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条特征线的出射中 子角通量:
由式(8)可知,欲求解特征线方程首先应该知道特征线信息, 包括特征线段长度sk和特征线段对应的截面信息。为获得这些信息, 通常有两种方法可以选择,一种是生成模块化特征线并存储各模块的 特征线参数;另一种是在线生成特征线并计算特征线参数。当采用第 一种方法时,对于模块几何结构比较复杂、几何尺寸比较大、或者模 块种类比较多的情况,都会大大增加计算机的存储负担。但是对于几 何较为复杂的反应堆堆芯,我们采用在线生成特征线的方法,其核心 思想是在特征线扫描的同时针对局部几何获取特征线参数。但是传统 的在线生成特征线方法计算效率非常低,原因是在特征线扫描之前不存储任何关于特征线的参数,所有参数都必须在线计算,使用过后随 机释放内存。这大大地限制了在线生成特征线方法在实际工程中的应 用。
由前文可知,传统的模块化特征线方法需要在进行特征线计算之 前存储模块内所有特征线的信息,这导致了它对于几何的适应性较 差。在线生成特征线方法可以有效地拓展特征线方法的几何适应性, 但是由于传统的在线生成特征线方法不存储特征线信息,每条长特征 线上的计算参数都在扫描时在线计算生成,并在扫描结束后抛弃以释 放计算机内存。在真实的核反应堆堆芯中,往往有成千上万个基本几 何单元(如轻水堆中的栅元结构、球床堆中的燃料球结构等)排列在 其中,每条长特征线往往会穿过几十至上百个这种基本几何单元。直 接将在线生成特征线方法应用到这类几何会带来大量的重复计算,会大大降低程序的计算效率。
发明内容
为了克服上述现有技术存在的问题,本发明的目的在于提供一种 加速在线生成核反应堆堆芯特征线的方法,该方法适用于存在重复基 本单元的核反应堆堆芯几何结构,考虑到现有反应堆堆型以及当前所 规划的下一代反应堆堆型,或者包含燃料栅元,或者包含燃料球,而 燃料栅元和燃料球在堆芯中都是大量重复存在的基本几何单元;因此 本发明方法对于待模拟的核反应堆堆型没有要求;为了方便,以包含 3个平源区的二维压水堆燃料栅元为例说明本发明的技术方案,但是 本发明同样适合其他包含多个平源区的二维或三维几何。
为了达到上述目的,本发明采用如下技术方案:
一种加速在线生成核反应堆堆芯特征线的方法,包括如下步骤:
步骤1:典型压水堆燃料栅元由内到外分别为燃料区、气隙区和 慢化剂区;描述一个典型压水堆燃料栅元的几何共需要8个面,分别 为:上边界面、下边界面、左边界面、右边界面、前边界面、后边界 面、气隙内表面和气隙外表面;其中气隙内表面和气隙外表面为圆柱 面,其余面均为平面;若考虑二维压水堆燃料栅元,认为栅元几何在 z方向无限,则上边界面和下边界面消失,同时,三维燃料栅元中的 “面”在二维几何中变成了“线”;所以二维压水堆燃料栅元几何由6条 线描述,分别为前边界线、后边界线、左边界线、右边界线以及气隙 的内边界线、外边界线;其中气隙内边界线和外边界线为同心圆,其 余边界线均为直线段。
为了叙述方便,给6条线分别编号如下:
其中左边界线为1号边界,前边界线为2号边界,右边界线为3 号边界,后边界线为4号边界,气隙内边界线为5号边界,气隙外边 界线为6号边界;
若以气隙内边界线和外边界线的圆心为坐标原点建立二维直角 坐标系,根据实际压水堆燃料栅元的尺寸,写出各条边界线的方程如 下:
1号边界:x=-0.63
2号边界:y=-0.63
3号边界:x=0.63
4号边界:y=0.63
5号边界:r=0.41
6号边界:r=0.47;
步骤2:在几何描述完成之后,对目标几何应用特征线方法进行 计算之前,应该产生特征线;特征线是穿过目标几何的一簇平行的等 距的直线,每条直线都被目标几何截成若干段,记作n段(n ≥1), 这些特征线段是进行特征线计算的基本单位;特征线段的空间位置由 角度和起始点坐标唯一确定,在二维几何中,角度被定义为特征线 与x正方向的夹角,特征线的角度由用户指定,此处以说明; 特征线段的起始点坐标为(x,y)的形式,对于同一条特征线,第n+1 段特征线的起始点即为第n段特征线的终点;由于存在大量的特征线 段,而每一条特征线的起始点各不相同,所以并不提前存储所有特征 线段的起始点坐标,仅存储每条特征线第一段的起始点坐标。由 经过数学换算,将穿过压水堆燃料栅元的平行特征线段统一 表示为:
y=0.5774x+b (9)
其中b特征线在边界上的截距,由特征线段的起始点确定;
步骤3:对于角度为30°的特征线,根据起始点的不同,从1号边 界出发的特征线的终点有两种可能:到达4号边界或者到达6号边界, 最终到达哪个边界取决于特征线在1号边界上的起始点坐标与1号边 界上分界点之间的关系;从2号边界出发的特征线的终点有两种可 能:到达6号边界或者到达3号边界,最终到达哪个边界取决于特征 线在2号边界上的起始点坐标与2号边界上分界点之间的关系;从3 号边界和4号边界出发的特征线将离开此基本几何单元,故不需考 虑;从5号边界出发的特征线的终点有两种可能:到达5号边界或者 到达6号边界,最终到达哪个边界取决于特征线在5号边界上的起始 点坐标与5号边界上分界点之间的关系;从6号边界出发的特征线的 终点有四种可能:到达6号边界、到达5号边界、到达4号边界以及 到达3号边界,最终到达哪个边界取决于特征线在6号边界上的起始 点坐标与6号边界上分界点之间的关系;具体对应关系见步骤4说明;
步骤4:求解各个边界上的分界点,所谓分界点,是指在特征线 上从该点的两侧出发,将到达不同的边界;具体如下:
1号边界:联立式(10)和式(11)
(y-y0)=0.5774(x+0.63) (10)
x2+y2=0.472 (11)
考虑直线和圆相切的情况,解得:
y0=0.1689或y0=-0.9166(舍去)
当-0.63≤y<0.1689时,特征线将和6号边界相交;当 0.1689≤y<0.63时,特征线将和4号边界相交;
2号边界:联立式(12)和式(13)
(y+0.63)=0.5774(x-x0) (12)
x2+y2=0.472 (13)
考虑直线和圆相切的情况,解得:
x0=-0.1512或x0=-2.0309(舍去)
当-0.63≤x<-0.1512时,特征线将和6号边界相交;当 -0.1512≤y<0.63时,特征线将和3号边界相交;
5号边界:过原点作与特征线方向垂直的辅助线a,由于以特征 线角度说明,所以辅助线a的角度为120°,辅助线a与5号边界 有2个交点,分别在第二象限与第四象限,记第二象限的交点为A, 第四象限的交点为B;将从A点逆时针到达B点所表示的圆弧记为将从B点逆时针到达A点所表示的圆弧记为从5号边界 出发的特征线,根据其出发点的位置不同,将到达不同的目标边界; 具体来说,从出发的特征线,将到达5号边界;从出发的特征 线,将到达6号面边界;已知5号边界是半径为0.41的圆,所以点A 和点B的坐标用极坐标表示为:A:B:
点A和点B确定之后,将各区间表示为:
对于从5号边界出发的特征线,设其起始点坐标为:(x,y),则: 当时,特征线将和5号边界相交;当 或时,特征线将和6号边界相 交;
6号边界:过原点作与特征线方向垂直的辅助线a,由于以特征 线角度说明,所以辅助线a的角度为120°,辅助线a与6号边界 有2个交点,分别在第二象限与第四象限,记第二象限的交点为A, 第四象限的交点为D;将从A点逆时针到达D点所表示的圆弧记为将从D点逆时针到达A点所表示的圆弧记为直线a与5 号边界也有2个交点,作经过这2个交点的5号边界的切线b和c, 切线b和c与的交点记为B和C。过3号边界和4号边界的交点 作平行于特征线方向的辅助线d,直线d与的交点记为E;这样, 6号边界被分为等5段;已知6号边界是半 径为0.47的圆,利用各条辅助线的方程,通过几何运算获得点A~点 E的坐标,用极坐标表示为:A:B:(0.47,4.1762), C:D:E:(0.47,1.0364)。
点A至点E确定之后,将各区间表示为:
对于从6号边界出发的特征线,设其起始点坐标为:(x,y),则: 当或时,特征线将和6号 边界相交;当时,特征线将和5号边界相交; 当或时,特征线将和3号边 界相交;当时,特征线将和5号边界相交;
步骤5:存储上述各个边界的坐标区间与特征线目标边界的对应 关系;
步骤6:在特征线扫描的过程中,根据所存储的对应关系,由特 征线在某个边界上的起点坐标,迅速检索到其目标边界;即能快速地 计算下一段特征线段的长度,将此参数提供给式(8)即能完成此段 特征线计算;
下面针对典型压水堆栅元,假设特征线辐角为给定特征线的 起始点坐标,描述此特征线的扫描以说明此过程:
假设特征线从1号边界的点A出发,给定A的坐标为: (-0.63,-0.31),根据所存储的1号边界上的对应关系:-0.63≤y<0.1689,可知特征线的目标边界是6号边界;通过特征线与6号边界求交点得到点B的坐标(-0.4282,-0.1935)以及AB的长 度0.233,与AB段的材料参数一起,即能够进行后续的特征线方程 计算;
由点B的坐标,根据所存储的6号边界上的对应关系: 可知特征线的目标边界是5号边界;通过特 征线与5号边界求交点得到点C的坐标(-0.376,-0.1634)以及BC的长 度0.0601,与BC段的材料参数一起,即能够进行后续的特征线方程 计算;
由点C的坐标,根据所存储的5号边界上的对应关系: 可知特征线的目标边界是5号边界;通过特征 线与5号边界求交点得到点D的坐标(0.3296,0.244)以及CD的长度 0.8148,与CD段的材料参数一起,即能够进行后续的特征线方程计算;
由点D的坐标,根据所存储的5号边界上的对应关系: 或可知特征线的目标边界是6 号边界;通过特征线与6号边界求交点能够得到点E的坐标 (0.3818,0.2741)以及DE的长度0.0601,与DE段的材料参数一起,即能够进行后续的特征线方程计算;
由点E的坐标,根据所存储的6号边界上的对应关系: 或可知特征线的目标边界 是3号边界;通过特征线与3号边界求交点得到点F的坐标 (1.26,0.4174)以及EF的长度0.8898,与EF段的材料参数一起,即能 够进行后续的特征线方程计算;
通过以上过程,即完成此燃料栅元几何中的特征线扫描计算。
和现有技术相比较,本发明具备如下优点:
采取传统的在线生成特征线方法对全堆芯生成特征线时,一条长 特征线被复杂的堆芯几何分割成了大量的特征线段,获取每条特征线 段的长度信息的方式为:已知特征线起始点,全局检索特征线的目标 边界(特征线的起始点坐标和方向求解与每个边界的距离,距离最短 者为目标边界)。由于压水堆堆芯中存在成千上万个边界面,所以此 过程将耗费海量的计算代价。
本发明考虑到在真实的压水堆堆芯中,充满了大量的燃料栅元结 构,这些燃料栅元结构大量重复出现;将这些大量重复出现的几何结 构视为典型几何单元,存储典型几何单元的特征信息,在获取目标边 界的时候,不需要全局搜索目标边界,而是根据所存储的对应关系直 接获取。这样使得生成特征线这一过程的计算代价节省了成千上万 倍。
附图说明
图1为5号边界上的分界点。
图2为6号边界上的分界点。
图3为特征线扫描过程。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明:
本实施例一种加速在线生成核反应堆堆芯特征线的方法,包括如 下步骤:
步骤1:典型压水堆燃料栅元由内到外分别为燃料区、气隙区和 慢化剂区。描述一个典型压水堆燃料栅元的几何共需要8个面,分别 为:上边界面、下边界面、左边界面、右边界面、前边界面、后边界 面、气隙内表面和气隙外表面。其中气隙内、外表面为圆柱面,其余面均为平面。若考虑二维压水堆燃料栅元,认为栅元几何在z方向无 限,则上、下边界面消失。同时,三维燃料栅元中的“面”在二维几何 中变成了“线”。所以可以说:二维压水堆燃料栅元几何可由6条线描 述。分别为前、后、左、右边界线以及气隙的内、外边界线。其中气 隙内、外边界线为同心圆,其余边界线均为直线段。
为了叙述方便,不妨给6条线分别编号。
其中左边界线为1号边界,前边界线为2号边界,右边界线为3 号边界,后边界线为4号边界,气隙内边界线为5号边界,气隙外边 界线为6号边界。
若以气隙内、外边界线的圆心为坐标原点建立二维直角坐标系, 根据实际压水堆燃料栅元的尺寸,可以写出各条边界线的方程如下。
1号边界:x=-0.63
2号边界:y=-0.63
3号边界:x=0.63
4号边界:y=0.63
5号边界:r=0.41
6号边界:r=0.47
步骤2:在几何描述完成之后,对目标几何应用特征线方法进行 计算之前,应该产生特征线。特征线是穿过目标几何的一簇平行的等 距的直线,每条直线都被目标几何截成n段(n≥1)。这些特征线段 是进行特征线计算的基本单位。特征线段的空间位置由角度和起始 点坐标唯一确定。在二维几何中,角度被定义为特征线与x正方向的 夹角。特征线的角度由用户指定,此处以为例。特征线段的起 始点坐标为(x,y)的形式,对于同一条特征线,第n+1段特征线的起 始点即为第n段特征线的终点。由于存在大量的特征线段,而每一条 特征线的起始点各不相同,所以并不提前存储所有特征线段的起始点 坐标。仅存储每条特征线第一段的起始点坐标。由经过简单 的数学换算,可以将穿过压水堆燃料栅元的平行特征线段统一表示 为:
y=0.5774x+b (9)
其中b由特征线段的起始点确定。
步骤3:对于角度为30°的特征线,根据起始点的不同,从1号边 界出发的特征线的终点有两种可能:到达4号边界或者到达6号边界, 最终到达哪个边界取决于特征线在1号边界上的起始点坐标与1号边 界上分界点之间的关系;从2号边界出发的特征线的终点有两种可 能:到达6号边界或者到达3号边界,最终到达哪个边界取决于特征 线在2号边界上的起始点坐标与2号边界上分界点之间的关系;从3 号边界和4号边界出发的特征线将离开此基本几何单元,故不需考 虑;从5号边界出发的特征线的终点有两种可能:到达5号边界或者 到达6号边界,最终到达哪个边界取决于特征线在5号边界上的起始 点坐标与5号边界上分界点之间的关系;从6号边界出发的特征线的 终点有四种可能:到达6号边界、到达5号边界、到达4号边界以及 到达3号边界,最终到达哪个边界取决于特征线在6号边界上的起始 点坐标与6号边界上分界点之间的关系。具体对应关系见步骤4说明;
步骤4:求解各个边界上的分界点,所谓分界点,是指在特征线 从该点的两侧出发,将到达不同的边界;具体如下:
1号边界:联立式(10)和式(11)
(y-y0)=0.5774(x+0.63) (10)
x2+y2=0.472 (11)
考虑直线和圆相切的情况,解得:
y0=0.1689或y0=-0.9166(舍去)
当-0.63≤y<0.1689时,特征线将和6号边界相交;当 0.1689≤y<0.63时,特征线将和4号边界相交。
2号边界:联立式(12)和式(13)
(y+0.63)=0.5774(x-x0) (12)
x2+y2=0.472 (13)
考虑直线和圆相切的情况,可以解得:
x0=-0.1512或x0=-2.0309(舍去)
当-0.63≤x<-0.1512时,特征线将和6号边界相交;当 -0.1512≤y<0.63时,特征线将和3号边界相交。
5号边界:如图1所示,过原点作与特征线方向垂直的辅助线a, 由于以特征线角度为例,所以辅助线a的角度为120°,辅助线a与 5号边界有2个交点,分别在第二象限与第四象限,记第二象限的交 点为A,第四象限的交点为B。将从A点逆时针到达B点所表示的圆弧记为将从B点逆时针到达A点所表示的圆弧记为从5 号边界出发的特征线,根据其出发点的位置不同,将到达不同的目标 边界。具体来说,从出发的特征线,将到达5号边界;从出发 的特征线,将到达6号面边界。已知5号边界是半径为0.41的圆, 所以点A和点B的坐标可以用极坐标表示为:A: B:
点A和点B确定之后,可以将各区间表示为:
对于从5号边界出发的特征线,设其起始点坐标为:(x,y),则: 当时,特征线将和5号边界相交;当 或时,特征线将和6号边界相交。
6号边界:如图2所示,过原点作与特征线方向垂直的辅助线a, 由于以特征线角度为例,所以辅助线a的角度为120°,辅助线a与 6号边界有2个交点,分别在第二象限与第四象限,记第二象限的交 点为A,第四象限的交点为D。将从A点逆时针到达D点所表示的圆弧记为将从D点逆时针到达A点所表示的圆弧记为直 线a与5号边界也有2个交点,作经过这2个交点的5号边界的切线 b和c,切线b和c与的交点记为B和C。过3号边界和4号边界 的交点作平行于特征线方向的辅助线d,直线d与的交点记为E。 这样,6号边界被分为等5段。已知6号边 界是半径为0.47的圆,利用各条辅助线的方程,通过简单的几何运 算可以获得点A~点E的坐标,用极坐标可以表示为:A: B:(0.47,4.1762),C:D:E:(0.47,1.0364)。
点A至点E确定之后,可以将各区间表示为:
对于从6号边界出发的特征线,设其起始点坐标为:(x,y),则: 当或时,特征线将和6号 边界相交;当时,特征线将和5号边界相交; 当或时,特征线将和3号边 界相交;当时,特征线将和5号边界相交。
步骤5:存储上述各个边界的坐标区间与特征线目标边界的对应 关系;
步骤6:在特征线扫描的过程中,根据所存储的对应关系,由特征线在某个边界上的起点坐标,迅速检索到其目标边界。即可快速地 计算下一段特征线段的长度,将此参数将提供给式(8)即可完成此段特征线计算。下面针对典型压水堆栅元,假设特征线辐角为给定特征线的起始点坐标,描述此特征线的扫描以说明此过程。
如图3所示,假设特征线从1号边界的点A出发,给定A的坐标为:(-0.63,-0.31),根据所存储的1号边界上的对应关系:-0.63≤y<0.1689,可知特征线的目标边界是6号边界。通过特征线与6号边界求交点可以得到点B的坐标(-0.4282,-0.1935)以及AB 的长度0.233,与AB段的材料参数一起,即可进行AB段的特征线 方程计算。
由点B的坐标,根据所存储的6号边界上的对应关系: 可知特征线的目标边界是5号边界。通过特 征线与5号边界求交点可以得到点C的坐标(-0.376,-0.1634)以及BC 的长度0.0601,与BC段的材料参数一起,即可进行BC段的特征线方程计算。
由点C的坐标,根据所存储的5号边界上的对应关系: 可知特征线的目标边界是5号边界。通过特征 线与5号边界求交点可以得到点D的坐标(0.3296,0.244)以及CD的长 度0.8148,与CD段的材料参数一起,即可进行CD段的特征线方程 计算。
由点D的坐标,根据所存储的5号边界上的对应关系: 或可知特征线的目标边界是6 号边界。通过特征线与6号边界求交点可以得到点E的坐标 (0.3818,0.2741)以及DE的长度0.0601,与DE段的材料参数一起,即可进行DE段的特征线方程计算。
由点E的坐标,根据所存储的6号边界上的对应关系: 或可知特征线的目标边界 是3号边界。通过特征线与3号边界求交点可以得到点F的坐标 (1.26,0.4174)以及EF的长度0.8898,与EF段的材料参数一起,即可进行EF段的特征线方程计算。
通过以上过程,即可完成此燃料栅元几何中的特征线扫描计算。
Claims (1)
1.一种加速在线生成核反应堆堆芯特征线的方法,其特征在于:包括如下步骤:
步骤1:典型压水堆燃料栅元由内到外分别为燃料区、气隙区和慢化剂区;描述一个典型压水堆燃料栅元的几何共需要8个面,分别为:上边界面、下边界面、左边界面、右边界面、前边界面、后边界面、气隙内表面和气隙外表面;其中气隙内表面和气隙外表面为圆柱面,其余面均为平面;若考虑二维压水堆燃料栅元,认为栅元几何在z方向无限,则上边界面和下边界面消失,同时,三维燃料栅元中的“面”在二维几何中变成了“线”;所以二维压水堆燃料栅元几何由6条线描述,分别为前边界线、后边界线、左边界线、右边界线以及气隙的内边界线、外边界线;其中气隙内边界线和外边界线为同心圆,其余边界线均为直线段;
为了叙述方便,给6条线分别编号如下:
其中左边界线为1号边界,前边界线为2号边界,右边界线为3号边界,后边界线为4号边界,气隙内边界线为5号边界,气隙外边界线为6号边界;
若以气隙内边界线和外边界线的圆心为坐标原点建立二维直角坐标系,根据实际压水堆燃料栅元的尺寸,写出各条边界线的方程如下:
1号边界:x=-0.63
2号边界:y=-0.63
3号边界:x=0.63
4号边界:y=0.63
5号边界:r=0.41
6号边界:r=0.47;
步骤2:在几何描述完成之后,对目标几何应用特征线方法进行计算之前,应该产生特征线;特征线是穿过目标几何的一簇平行的等距的直线,每条直线都被目标几何截成n段,n≥1,这些特征线段是进行特征线计算的基本单位;特征线段的空间位置由角度和起始点坐标唯一确定,在二维几何中,角度被定义为特征线与x正方向的夹角,当特征线的角度特征线段的起始点坐标为(x,y)的形式,对于同一条特征线,第n+1段特征线的起始点即为第n段特征线的终点;由于存在大量的特征线段,而每一条特征线的起始点各不相同,所以并不提前存储所有特征线段的起始点坐标,仅存储每条特征线第一段的起始点坐标;由经过数学换算,将穿过压水堆燃料栅元的平行特征线段统一表示为:
y=0.5774x+b (9)
其中b为特征线在边界上的截距,由特征线段的起始点确定;
步骤3:对于角度为30°的特征线,根据起始点的不同,从1号边界出发的特征线的终点有两种可能:到达4号边界或者到达6号边界,最终到达哪个边界取决于特征线在1号边界上的起始点坐标与1号边界上分界点之间的关系;从2号边界出发的特征线的终点有两种可能:到达6号边界或者到达3号边界,最终到达哪个边界取决于特征线在2号边界上的起始点坐标与2号边界上分界点之间的关系;从3号边界和4号边界出发的特征线将离开此基本几何单元,故不需考虑;从5号边界出发的特征线的终点有两种可能:到达5号边界或者到达6号边界,最终到达哪个边界取决于特征线在5号边界上的起始点坐标与5号边界上分界点之间的关系;从6号边界出发的特征线的终点有四种可能:到达6号边界、到达5号边界、到达4号边界以及到达3号边界,最终到达哪个边界取决于特征线在6号边界上的起始点坐标与6号边界上分界点之间的关系;具体对应关系见步骤4说明;
步骤4:求解各个边界上的分界点,所谓分界点,是指在特征线上从该点的两侧出发,将到达不同的边界;具体如下:
1号边界:联立式(10)和式(11)
(y-y0)=0.5774(x+0.63) (10)
x2+y2=0.472 (11)
考虑直线和圆相切的情况,解得:
y0=0.1689或y0=-0.9166,舍去
当-0.63≤y<0.1689时,特征线将和6号边界相交;当0.1689≤y<0.63时,特征线将和4号边界相交;
2号边界:联立式(12)和式(13)
(y+0.63)=0.5774(x-x0) (12)
x2+y2=0.472 (13)
考虑直线和圆相切的情况,解得:
x0=-0.1512或x0=-2.0309,舍去
当-0.63≤x<-0.1512时,特征线将和6号边界相交;当-0.1512≤y<0.63时,特征线将和3号边界相交;
5号边界:过原点作与特征线方向垂直的辅助线a,由于特征线角度所以辅助线a的角度为120°,辅助线a与5号边界有2个交点,分别在第二象限与第四象限,记第二象限的交点为A,第四象限的交点为B;将从A点逆时针到达B点所表示的圆弧记为将从B点逆时针到达A点所表示的圆弧记为从5号边界出发的特征线,根据其出发点的位置不同,将到达不同的目标边界;具体来说,从出发的特征线,将到达5号边界;从出发的特征线,将到达6号面边界;已知5号边界是半径为0.41的圆,所以点A和点B的坐标用极坐标表示为:
点A和点B确定之后,将各区间表示为:
对于从5号边界出发的特征线,设其起始点坐标为:(x,y),则:当时,特征线将和5号边界相交;当或时,特征线将和6号边界相交;
6号边界:过原点作与特征线方向垂直的辅助线a,由于特征线角度所以辅助线a的角度为120°,辅助线a与6号边界有2个交点,分别在第二象限与第四象限,记第二象限的交点为A,第四象限的交点为D;将从A点逆时针到达D点所表示的圆弧记为将从D点逆时针到达A点所表示的圆弧记为直线a与5号边界也有2个交点,作经过这2个交点的5号边界的切线b和c,切线b和c与的交点记为B和C;过3号边界和4号边界的交点作平行于特征线方向的辅助线d,直线d与的交点记为E;这样,6号边界被分为等5段;已知6号边界是半径为0.47的圆,利用各条辅助线的方程,通过几何运算获得点A~点E的坐标,用极坐标表示为:
点A至点E确定之后,将各区间表示为:
对于从6号边界出发的特征线,设其起始点坐标为:(x,y),则:当或时,特征线将和6号边界相交;当时,特征线将和5号边界相交;当或时,特征线将和3号边界相交;当时,特征线将和5号边界相交;
步骤5:存储上述各个边界的坐标区间与特征线目标边界的对应关系;
步骤6:在特征线扫描的过程中,根据所存储的对应关系,由特征线在某个边界上的起点坐标,迅速检索到其目标边界;即能快速地计算下一段特征线段的长度,将此参数提供给式(8)即能完成此段特征线计算;
下面针对典型压水堆栅元,当特征线辐角为给定特征线的起始点坐标,描述此特征线的扫描以说明此过程:
假设特征线从1号边界的点A出发,给定A的坐标为:(-0.63,-0.31),根据所存储的1号边界上的对应关系:-0.63≤y<0.1689,可知特征线的目标边界是6号边界;通过特征线与6号边界求交点得到点B的坐标(-0.4282,-0.1935)以及AB的长度0.233,与AB段的材料参数一起,即能够进行后续的特征线方程计算;
由点B的坐标,根据所存储的6号边界上的对应关系:可知特征线的目标边界是5号边界;通过特征线与5号边界求交点得到点C的坐标(-0.376,-0.1634)以及BC的长度0.0601,与BC段的材料参数一起,即能够进行后续的特征线方程计算;
由点C的坐标,根据所存储的5号边界上的对应关系:可知特征线的目标边界是5号边界;通过特征线与5号边界求交点得到点D的坐标(0.3296,0.244)以及CD的长度0.8148,与CD段的材料参数一起,即能够进行后续的特征线方程计算;
由点D的坐标,根据所存储的5号边界上的对应关系:或可知特征线的目标边界是6号边界;通过特征线与6号边界求交点能够得到点E的坐标(0.3818,0.2741)以及DE的长度0.0601,与DE段的材料参数一起,即能够进行后续的特征线方程计算;
由点E的坐标,根据所存储的6号边界上的对应关系:或可知特征线的目标边界是3号边界;通过特征线与3号边界求交点得到点F的坐标(1.26,0.4174)以及EF的长度0.8898,与EF段的材料参数一起,即能够进行后续的特征线方程计算;
通过以上过程,即完成此燃料栅元几何中的特征线扫描计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810499010.2A CN108734181B (zh) | 2018-05-23 | 2018-05-23 | 一种加速在线生成核反应堆堆芯特征线的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810499010.2A CN108734181B (zh) | 2018-05-23 | 2018-05-23 | 一种加速在线生成核反应堆堆芯特征线的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108734181A CN108734181A (zh) | 2018-11-02 |
CN108734181B true CN108734181B (zh) | 2019-04-16 |
Family
ID=63935188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810499010.2A Active CN108734181B (zh) | 2018-05-23 | 2018-05-23 | 一种加速在线生成核反应堆堆芯特征线的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108734181B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111696197B (zh) * | 2020-05-27 | 2023-10-27 | 北京科技大学 | 一种快堆六角形堆芯中子输运计算的射线生成与追踪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096183A (zh) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | 一种基于特征线法的多重并行方法 |
CN106202613A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种适用于耦合输运计算的大规模并行方法 |
CN107085563A (zh) * | 2017-04-27 | 2017-08-22 | 西安交通大学 | 一种在线生成特征线的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4664399B2 (ja) * | 2008-07-31 | 2011-04-06 | 財団法人発電設備技術検査協会 | 高Cr系鋼構造物のじん性評価方法 |
EP3062312A4 (en) * | 2013-10-21 | 2017-06-14 | Hitachi, Ltd. | Nuclear reactor core |
-
2018
- 2018-05-23 CN CN201810499010.2A patent/CN108734181B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096183A (zh) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | 一种基于特征线法的多重并行方法 |
CN106202613A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种适用于耦合输运计算的大规模并行方法 |
CN107085563A (zh) * | 2017-04-27 | 2017-08-22 | 西安交通大学 | 一种在线生成特征线的方法 |
Non-Patent Citations (5)
Title |
---|
NECP-X中角度与特征线并行的研究与实现;马党伟等;《原子能科学技术》;20180331;第52卷(第3期) |
三维模块化特征线方法;刘宙宇等;《核动力工程》;20101231;第31卷(第S2期);21-24 |
以栅元为模块进行特征线跟踪的中子输运方程解法;汤春桃等;《核动力工程》;20090831;第30卷(第4期);32-36 |
先进中子学栅格程序KYLIN-2中特征线方法模块的开发与验证;柴晓明等;《核动力工程》;20160830;第37卷(第4期);154-159 |
基于模块化特征线方法的二维/一维祸合输运程序开发;梁亮等;《核动力工程》;20141231;第35卷(第S2期);131-134 |
Also Published As
Publication number | Publication date |
---|---|
CN108734181A (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103578133A (zh) | 一种将二维图像信息进行三维重建的方法和设备 | |
CN105844067B (zh) | 一种船体曲板焊接变坡口自动加工数据的获取方法 | |
CN105718643A (zh) | 一种基于最优化视角的船舶生产设计审图器实现方法 | |
CN109344459A (zh) | 一种蒙特卡罗程序的可视建模与转换方法及系统 | |
Liu et al. | Constant scallop-height tool path generation for three-axis discrete data points machining | |
CN113673186A (zh) | 一种基于stl文件的笛卡尔网格快速生成方法 | |
CN110362913A (zh) | 基于Phi函数进行干涉计算的卫星组件布局优化设计方法 | |
CN109859287A (zh) | 一种自动生成矢量扇形的方法 | |
CN108734181B (zh) | 一种加速在线生成核反应堆堆芯特征线的方法 | |
CN105785919B (zh) | 点云五轴无全局干涉刀轴矢量可行域计算方法 | |
CN104392030B (zh) | 一种基于stl三维模型的曲面加工方法 | |
An et al. | Self-adaptive polygon mesh reconstruction based on ball-pivoting algorithm | |
CN109598785A (zh) | 一种三维网格模型视图转换方法 | |
US8180607B2 (en) | Acoustic modeling method | |
CN106846483B (zh) | 基于共享直线段的平面片束排序方法和系统 | |
CN105929792A (zh) | 点云模型五轴无局部干涉刀轴偏角范围计算方法 | |
Liu et al. | Accelerating Large-Scale CFD Simulations with Lattice Boltzmann Method on a 40-Million-Core Sunway Supercomputer | |
Wang et al. | User controllable anisotropic shape distribution on 3D meshes | |
CN109979007A (zh) | 一种建筑体的几何造型方法和装置 | |
CN103985150B (zh) | 一种图像处理方法及电子设备 | |
CN109658489A (zh) | 一种基于神经网络的立体网格数据处理方法和系统 | |
CN116542183B (zh) | 非结构二阶曲面多面体网格重构方法、系统、设备及介质 | |
Jiang et al. | Three-dimensional garment surface reconstruction based on ball-pivoting algorithm | |
CN116051782B (zh) | 基于正交网格曲线插值的数据处理及重构建模方法、设备和存储介质 | |
US20060020456A1 (en) | Object oriented acoustic modeling tools and methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |