Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
藤村 統一郎*; 奥村 啓介
JAERI-Research 2002-024, 27 Pages, 2002/11
JAERI-Research-2002-024.pdf:1.04MB
低減速スペクトル炉等の6角形状の炉心を解析する拡散コードの原型版を開発し、その反復解法を高速化するため、さまざまな収束加速法の適用性について検討した。本3次元コードMOSRA-Prismは、6角形状の炉心を正3角柱に分割し、その中の中性子束分布を3次の多項式で近似する多項式展開ノード法に基づいている。多群拡散コードとしての反復解法は、通常の内側反復法と外側反復法を採用するが、内側反復に適応的加速法、外側反復に中性子源外挿法を適用し、その有効性を確認した。本報告書では、コードの数値解法の元となる多項式展開ノード法の定式化の概要を説明するとともに、さまざまなサンプル計算で得られた、収束加速法の局所的な効率及び全体的な効率について検討する。また、コード開発過程で新たに導出した真空境界条件の一般的な記述法についても述べる。
黒石 武; 野村 靖
Proceedings of International Conference on the New Frontiers of Nuclear Technology; Reactor Physics, Safety and High-Performance Computing (PHYSOR 2002) (CD-ROM), 10 Pages, 2002/10
核分裂源分布の収束緩慢性を加速するため、マトリクスK計算が開発され、従来のモンテカルロ計算に導入されてきた。核分裂源が未収束であるモンテカルロ計算の途中段階において核的結合係数を近似的に求めることが出来れば、核分裂源行列方程式の固有ベクトルを用いて核分裂源を補正することにより加速が実施される。本論文では弱結合相互干渉系に対するマトリクスK計算の効果的な2つの適用手法、即ち、繰返し加速法とソース生成法を提案する。前者はマトリクスK計算による加速手順を単純に繰返すものであり、照射済みピンセル体系に対する計算結果は、臨界性を統計評価する上での信頼できる核分裂源を得るという十分な加速効果を示した。しかしながら、ある種の弱結合マルチユニット体系に対しては、ソースレベルの低いユニットが多数あるために、収束に至るための2回以上のマトリクスK計算を繰返す手順が実施できないかも知れない。後者は、このような場合に適用すべく新たに検討したものである。チェッカーボード燃料貯蔵ラック体系はそのような典型例の一つであり、計算結果により本手法の有効性が示された。
藤村 統一郎
JAERI-Data/Code 95-019, 54 Pages, 1996/01
JAERI-Data-Code-95-019.pdf:1.47MB
複雑形状の炉心を解析するため、二重有限要素法による3次元中性子輸送コードDFEMを開発した。基本となる解法としては、空間と角度の変数に有限要素法を適用する二重有限要素法を用いている。炉心体系を多角柱でモデル化し、その中の中性子束の空間分布を1次の基底関数で近似するが、角度変数に対しては様々な基底関数について比較し、特徴を明らかにする。また、高精度化と高速化のため、反射境界打ち切り誤差について一般的な修正法を導出するとともに、反復解法に改良を加えた。本プログラムの構造について概説し、使用法を記すとともに、検証計算例として実規模の原子炉を含む様々な問題を採り上げ、中性子束分布等について他のコードによる計算結果と比較し、検討を行った。
藤村 統一郎; 松村 正弘*; 中原 康明
JAERI-M 83-144, 40 Pages, 1983/09
本稿は、三次元幾何形状における定常、多群中性子輸送問題に対する二重有限要素法に基づくガレルキン法のアルゴリズムとその特徴について述べたものである。定式化においては、現実の原子炉の形状をできる限り正確に模擬するため、空間要素として三角柱要素と四角柱要素の組合せを採用すると共に、中性子束の角度分布を滑らかに表現するため、角度空間において重なりをもつ六つの基底を採用している。本解法の特徴は境界条件を陽に記述すること、および平面の層に沿って次々と節点を走査する反復法にあり、その収束加速法の新しい提案もなされる。この解法に基づく計算コードが開発され、その概要についての解説も示される。また、今迄に実施した実在規模の問題を含む計算の経験に基づき、差分法のコードとの比較をもとに二重有限要素解の特徴が示される。
藤村 統一郎; 松井 泰*
Nuclear Science and Engineering, 77, p.360 - 367, 1981/00
被引用回数：2 パーセンタイル：44.84(Nuclear Science & Technology)中性子拡散コードにおける内側反復解法に対して適応的加速法を応用した場合の有効性が論じられる。 この方法はアルゴリズムが簡単で、反復行列が非負定値である場合の定常一階線形反復法の殆んどに対してその収束を加速することが知られている。 ここでは3次元中性子拡散コードとして、その内側反復解法にSOR法を用いた有限要素法によるものと、ADI法を用いた有限差分法によるものがとり上げられる。 これらに適応的加速法を応用した結果に対して数値的な検討が行なわれ、元の反復法の加速因子が不適切にとられたとき特に効果的なことが示される。 また、新たに、反復行列が対称で非負定値でない場合にもこの方法が拡大して応用できることが示される。
朝岡 卓見; 中原 康明; 堀上 邦彦; 西田 雄彦; 鈴木 忠和; 田次 邑吉; 宮坂 駿一; 弘田 実彌
Nuclear Science and Engineering, 59(4), p.326 - 336, 1976/04
被引用回数：3モンテカルロ法で原子炉体系の固有値を求める際の反復計算過程の収束加速のため、粗メッシュ再釣合法が導入された。1バッチの中性子ヒストリーについての計算終了毎に、原子炉の各粗メッシュ領域で中性子の釣合が保たれるように中性子束に掛けられる因子が計算される。この再釣合因子を、次のバッチ計算の最初に、各粗メッシュ領域の核分裂中性子源の重み(強度)に掛ける。この粗メッシュ再釣合法を使った計算例は、この方法が各粗メッシュについての中性子源とか損失を、通常のモンテカルロ計算より正しく求める、新しいサンプリング法であることを示している。
