基于非结构网格的中子输运方程的计算方法研究.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 基于非结构网格的中子输运方程的 计算方法研究水 吴宏春曹良志巨海涛卢浩亮刘国明 (西安交通大学，710049) 摘要：本文针对非结构网格下的中子输运方程求解方法进行了研究，推导了Pn有限元，Sn有限元。三角 形穿透几率方法，和三角形节块方法的计算模型，并分别研制了相应的计算程序，利用基准题进行了验证， 分析比较了各个方法的优缺点． 关键词：中子输运非结构网格计算方法 1引言 中子输运方程的数值计算是先进核反应堆设计的核心研究课题，然而，中子输运方程是 极其复杂的，它是一个与介质特性、空间位置、中子能量、和中子运动方向等白变量有关的 高阶偏微分方程。特别是在强非均匀介质和强各向异性散射的情况下，更难精确求解。 近几十年来，由于核电事业的发展需要，中子输运方程的数值求解方法有了较大的发展， 并有了精度较高的计算软件。比较有代表性的有细网差分离散坐标方法程序DOT和粗网节 块离散坐标方法程序DNTM。但这些程序都是以现有核反应堆中的中子输运模拟为目的而研 制的，只能运用于像现有核反应堆这样的介质稳定(不流动)、结构化几何区域、中子散射各 向异性不强等一类问题的求解。 为了能满足我国核电和核动力的发展需求，必须尽快研制新一代先进核反应堆，这种新 型核反应堆的几何形状无疑将更加复杂，介质的不稳定性加强，中子散射的各向异性增强等。 这些因素使得模拟中子运动的中子输运方程的求解变得非常复杂，现有的中子输运计算软件 均无法满足要求。加上国际之间的竞争非常激烈，不可能在近期内从国外引进。所以开展基 于非结构网格的中子输运方程的计算方法研究已成为迫切的基础研究课题，该研究对新一代 核反应堆的设计研究具有重要意义。 2非结构网格中子输运计算方法 2．1有限元Pn方法 SN方法数值过程相对简单。发展比较成熟，在反应堆中应用较广，如ANISN[11、DOTi2l 和DNTM[31等程序就是采用sN方法，但它有固有的缺点【4l：(1)对强各向异性介质问题精度 不够，因此必须通过增加方向离散点的个数来提高计算精度。但是这样一来，为了保证方法 的收敛，就必须同时要求把空间网格划分得非常小．这无疑将导致计算量的急剧增加：(2) ’国家自然科学基金资助项目 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 存在‘‘射线效应”，即由于方法本身的固有特性所导致的在强吸收介质中解的严重失真．这些 缺点严重限制了离散纵标方法在更广的范围内(特别是复杂情况)的应用。PN方法能够避免 ‘‘射’线效应”，对于非均匀性较强、各项异性强的介质等复杂输运问题，理论上只要提高阶数 N，可以达到任意高阶的精度。但当阶数提高时．方程将变得非常复杂，而且阶数不同，利 用常规方法导出的方程和求解过程均不相同，编制可适用于不同阶数N的通用程序变得十分 困难。 球谐函数方法最早由Wick和Marshak应用于中子输运问题中，后经Mark详加发展，在 20世纪50年代一度成为了数值求解中子输运方程的主流方法15J。但是由于早期对球谐函数 方法的研究I¨】主要限于特定的阶数如Pl或者P3，限制了PN方法的发展，加之常规反应堆中 sN方法能保证足够的精度，因此PN方法曾经一度未能得到很好地发展。近年来，随着计算 条件的不断改善，复杂环境下的输运问题对计算精度的要求越来越高，国际上对PN方法的研 究又重新重视起来。 。传统的PN方法研究主要是基于一阶输运方程或者二阶奇偶自共轭输运方程，1999年， Flux)，并指出与传统的一阶输运方程和二阶奇偶自共轭方程相比，SAAF方程具有很多独特 的数值求解优点。这为PN方法的研究开辟了一条新的道路。 SAAF方程的基本形式为， 一力·vcz,一s)一1口·re(，’E，Q)+(五一s)≯(，，E，口)2 (1) q(r，E，．q)-D·V(互-s)一g(，，E，口) 中子角通量密度用球谐函数展开为 ∞ 疗