Prediction of fission product release during the LOFC experiments at the HTTR

HTTRのLOFC試験における核分裂生成物の放出予測

Shi, D.*; Xhonneux, A.*; 植田 祥平   ; Verfondern, K.*; Allelein, H.-J.*

Shi, D.*; Xhonneux, A.*; Ueta, Shohei; Verfondern, K.*; Allelein, H.-J.*

高温ガス炉技術の安全性と当該事象下で想定される物理現象の確証を目的として、高温工学試験研究炉(HTTR)を用いた実証試験が行われた。OECD主導によるLOFC(冷却材喪失事故)プロジェクトの一部として、HTTRを用いた3つのシリーズ試験運転が計画された。本試験運転では、ガス循環機の全台停止によって1次冷却材流量を異常低下させ、原子炉をスクラム不能かつ反応度制御しないATWS(原子炉スクラム失敗を伴う予期された過渡事象)とする計画であり、本事象下において原子炉の再臨界を伴う。本論では、ユーリッヒ研究所が新たに開発したソースターム解析コードSTACYについて、またLOFC試験条件下でのHTTR炉心における核分裂生成物挙動の計算結果について述べる。STACYは元来、核分裂生成物の移行放出挙動シミュレーションの検証・妥当性確認モデルであり、本研究において中空円筒形状の燃料コンパクトや六角柱状ブロック型炉心を取扱えるよう拡張改良した。本論では、3次元連続エネルギーモンテカルロ燃焼計算コードSERPENTならびに事故・過渡事象時の核分裂生成物移行放出挙動の多群解析モデルMGTによる時間依存の核熱流動解析値にもとづくSTACYの計算結果について述べる。

Demonstration tests were conducted using the High Temperature Engineering Test Reactor (HTTR) in Oarai, Japan, to confirm the safety of HTGR technologies and assure the expected physical phenomena to occur under given conditions. As part of the OECD directed LOFC (loss of forced cooling) project, a series of three tests at the HTTR has been planned with tripping of all gas circulators while deactivating all reactor reactivity control to disallow reactor scram due to abnormal reduction of primary coolant flow rate. The tests fall into anticipated transient without scram (ATWS) with occurrence of reactor recriticality. The paper will describe the Source Term Analysis Code System (STACY) newly developed at the Research Center Jlich and present the results of fission product behavior in the HTTR core under the LOFC test conditions. STACY encompasses the original verified and validated computer models for simulating fission product transport and release. For verification of the modernized and extended version, it was assured that results obtained with the original tools could be reproduced. One of the new features of STACY is its ability to also treat fuel compacts of (full) cylindrical or annular shape and a complete prismatic block reactor core, respectively, supposed sufficient input data be available. In the paper, calculations are based on time-dependent neutronics and fluid dynamics results obtained with the Serpent and MGT models.