二次元拡散有限要素法コード"DIFEM"

A Tow dimensional finite element diffusion code "DIFEM"

西山 征夫*; 大竹 俊英*; 松野 義明*

not registered; not registered; Matsuno, Yoshiaki*

有限要素法にょる中性子拡散コード"DIFEM"を作成し,従来の差分コードでは扱い難い複雑な体系の核計算を行えるようにした。▲有限要素法は従来の差分法に較べて次の利点がある。▲1.計算体系のメッシュの粗密化が可能なため,計算精度の向上が図れる。▲2.任意に計算体系のメッシュをとれるので,複雑な計算体系が扱える。▲コードは,巨視断面積,微視断面積を与えて,2次元R―Z,X―Y,三角,R―体系の固有値計算が行える。▲プログラムはバリアブルディメンションを採用し,FORTRAN―4で書かれている。▲入力はFIDOフォーマットで行う。有限要素法の欠点である複雑な計算体系の入力の手間を軽減するため,オートメッシュ機能による入力も行える。又,入力チェック用に要素をプロッタで出力する機能を持っている。▲計算時間短縮のため,フラックスゲスが使え,再スタート機能を持っている。▲本文では主に計算方法の説明を記し,付録では入力の方法とコードの構成を記した。▲

A computer code ; "DIFEM" has been developed for solving neutron diffusion problem with a Finite Element Method, which has the following advantage points as conpared with the Finite Difference Method : (1)to perform a calculation with any mesh size in arbitrary regions (2)to treat any geometry. DIFEM code is programmed to treat two dimensional (RZ, XY, RO, Triangle) eigenvalue problems and to use the variable dimension to minimize the allocated memoty. A auto-mesh routine, plotting routine option and restart option are involved in the computor code to make mesh data and to save computing time. Input data should be written with "FIDO" format. Programs are written in "FORTRAN IV" language. The code description, and the input/output instruction will be stated in this report.