三角柱形状拡散ノード法コードにおける収束加速法の適用性の検討

藤村 統一郎*; 奥村 啓介

JAERI-Research 2002-024, 27 Pages, 2002/11

JAERI-Research-2002-024.pdf:1.04MB

低減速スペクトル炉等の6角形状の炉心を解析する拡散コードの原型版を開発し、その反復解法を高速化するため、さまざまな収束加速法の適用性について検討した。本3次元コードMOSRA-Prismは、6角形状の炉心を正3角柱に分割し、その中の中性子束分布を3次の多項式で近似する多項式展開ノード法に基づいている。多群拡散コードとしての反復解法は、通常の内側反復法と外側反復法を採用するが、内側反復に適応的加速法、外側反復に中性子源外挿法を適用し、その有効性を確認した。本報告書では、コードの数値解法の元となる多項式展開ノード法の定式化の概要を説明するとともに、さまざまなサンプル計算で得られた、収束加速法の局所的な効率及び全体的な効率について検討する。また、コード開発過程で新たに導出した真空境界条件の一般的な記述法についても述べる。

A double finite element method with accurate reflective boundary condition treatment for three-dimensional transport

藤村 統一郎

Computer Physics Communications, 82, p.111 - 119, 1994/00

https://doi.org/10.1016/0010-4655(94)90160-0

三次元中性子輸送問題を高精度で解くための定式化を行い、検証計算を実施した。二重有限要素法に基づいた変分法による定式化では、空間変数に一次関数を適用するとともに、角度変数に階段関数を適用している。ここに、多角柱形状の原子炉体系を反射条件を用いて解く場合、新しく開発した体系打切り誤差解消アルゴリズムにより、厳密な離散式が導出される。この解法の検証のため、二つの問題について計算を行った。MOZART炉心をモデル化した問題では、不規則な要素分割に対する非等方散乱計算においても、この解法が有効に働くことが示された。また、NEA-CRPのLWRをモデル化したベンチマーク問題では、粗いメッシュでも固有値や中性子束が、他の輸送方程式の解法による解と同精度に求まることが示された。