基于三角形网格的穿透概率方法研讨.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 基于三角形网格的穿透概率方法研究· 刘国明刘萍萍吴宏春 (西安交通大学，710049) 摘要：在现代反应堆的几何设计中很多采用了非结构几何的设计，对于非结构几何环境下中子输运方程的 求解，已有的蒙特卡罗方法和有限元方法存在速度不快的缺点。本文研究了基于三角形网格的穿透概率方法， 其三角形网格子区内部中子源采用平源分布，子区界面中子角通量密度空间分布设为均匀、角度分布采用简 化4只近似，采用一种按照非结构网格编码顺序进行网格扫描的新方法，并编制了程序1PTRI．对一些二维 能够很好的计算二维非结构几何的问题，结果准确，而且耗时更少． 关键词：中子输运穿透几率法三角形网格 1引言 穿透概率法(TPM：TransmissionMethod)是求解中子输运方程的基本方法之 Probability 一，已经普遍应用于核反应堆燃料组件参数的计算，它的基本思想是将系统分为若干子区， 且认为每个子区的材料性质是均匀的，并在每个子区建立中子输运方程；各子区之间通过界 面上中子流必须连续的条件来耦合，每个子区都只和其相邻的子区耦合，所以也称为界面流 方法。迄今为止，国内外已有许多用穿透几率法进行组件计算的程序，如TPFAP[¨、TPHEXl2J、 子区划分也多采用矩形p～sl和正六边形I盈。 随着核技术应用的不断发展和人们对计算精度要求的不断提高，越来越多的核技术问题 需要在非结构几何下进行求解中子输运方程，现有的矩形和六角形网格均不再适用，必须研 究基于非结构网格的中子输运方程的计算模型。R．T．AcKgd曾将有限元方法与PⅣ方法19J相结 合，研究了非结构网格下中子输运方程的求解，但由于方法计算复杂，计算量大，所以有必 要研究计算速度较快、模型相对简单的穿透几率方法在非结构网格中的应用。 由于任意形状的三角形网格可以对曲边或者多边形区域很好地逼近，因此，本文研究了 基于三角形网格的穿透概率法。本文建立了数值计算模型，研制了适合于二维任意不规则几 何区域中子输运计算程序TPTRI，并通过一系列的基准题的数值计算结果验证了该方法的正 确性和适用性。本文研究成果对复杂环境下的中子输运计算方法研究具有较高的学术价值， 对复杂反应堆堆芯设计具有良好的应用前景。 2理论模型 ’国家自然科学基金资助项目 第十一届反厘堆致值计算和粒子输运学术会议暨2006年反应堆物理会议 2．1基本方程 设将所研究的区域或系统划分为IEE个子区，且认为子区内的材料性质是均匀的。对于 任一子区可写出分群近似下的中子积分输运方程【1川： og(隅=他(只函也肛口而d／+①；(暑即州～ (1) 式中：。譬(只点)为尸处墨方向的中子角通量密度，Qg(P，点)为户处益方向的中子源密度， U ①暑一(￡，墨)为界面s上E处的入射中子角通量，，、to R、R在x、y平面上的投影 R=旷一Pl，R，=I尸一￡1)．如图1所示。所选用的坐标系和式中其他各量为习惯表示。 图l三角形子区的中子方向矢量 对每个三角形子区各界面和各项角按图2所示的方法进行编码。约定编码按照逆时针顺 序，顶角的编码就是其对边的相应编码。以三角形子区界面外法线萨为基轴，中子飞行方向 Ⅲ w 儒 Q在X-Y平面上的投影Q朋与扩夹角为矽，并规定逆时针为正，顺时针为负，依据缈值将 三角形界面上中子飞行方向划分为四个象限， 现将各象限描述如下(见图2)。 第l象限(q=1)：0≤缈≤州2 0 第2象限(q=2)：一州2≤伊s 第3象限(q=3)：州2≤伊≤石 第4象限(q=4)：--R≤伊≤一刀／2 w I,Mo t,U ． 分，可得子区界面k象限q的出射中子通量密度和出射中子流密度公式： 81 第十一届反应堆