Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
相馬 秀; 石垣 将宏*; 柴本 泰照
Annals of Nuclear Energy, 219, p.111455_1 - 111455_12, 2025/09
Containment venting is one of the accident mitigation measures during severe accidents in nuclear power plants for preventing overpressure failure of the containment vessels. Because of the capability of releasing hydrogen generated in the containment vessel, the hydrogen risk can be also reduced. In this study, we conducted experiments with the large-scale test facility CIGMA to investigate the light gas transport during the venting action, mainly focusing on the effect of sump water boiling caused by the vent. The CIGMA test vessel initially pressurized by steam, air, and helium (hydrogen simulant) that formed a helium-rich density stratification was depressurized with and without sump water, with different venting flow rates, and at different venting positions. As the sump water became a steam source due to flash boiling, the helium stratification was diluted and the venting time increased twofold compared to the case without sump water, which significantly affected the amount of helium discharged to the atmosphere. Especially for the high venting flow rate condition, the amount of helium remaining in the vessel at the end of depressurization was half that of the case without sump water. Lowering the venting position from within the initial stratification to 3 m below its interface led to a threefold increase in the amount of helium remaining at the same low pressure, because of the longer time until the helium-rich stratification reached the venting position.
Rizaal, M.; 中島 邦久; 鈴木 恵理子; 三輪 周平
Annals of Nuclear Energy, 218, p.111433_1 - 111433_10, 2025/08
The release of iodine in a case of severe nuclear accident is directly linked to short-term radiological consequences. This concern raises issues in understanding the chemical forms of the transported iodine to devise proper accident management measures/strategies. In contributing to such efforts, this work presents experimental and theoretical approaches to defining the impact of molybdenum as a semi-volatile fission product toward iodine speciation in the gas phase. Given humid atmospheric conditions with different oxygen potentials, the interactions were revealed through the reaction products consisting of both gas and aerosols upon their transport and condensation in the temperature range of 1150-450 K. Through thermodynamic equilibrium calculations where new thermodynamic data of cesium molybdates have been incorporated, the experimental observation was reproduced, hence general interaction mechanism was proposed in this work.
荒木 祥平; 會澤 栄寿; 村上 貴彦; 新垣 優; 多田 裕太; 神川 豊; 長谷川 健太; 吉川 智輝; 住谷 正人; 関 真和; et al.
Annals of Nuclear Energy, 217, p.111323_1 - 111323_8, 2025/07
被引用回数：0原子力機構では、臨界集合体STACYを均質溶液体系から非均質軽水減速体系へと更新した。STACY更新炉においても最大熱出力は200Wと定められており、熱出力校正は運転を行う上で重要である。熱出力測定においては、溶液系STACYで用いていたFPの分析による熱出力の評価が適応できなかったため、放射化法をベースとする実験データと数値計算を組み合わせて出力を評価する手法をSTACY更新炉の体系に適応し、測定を実施した。測定データを基に出力校正を実施した結果、校正後の指示値は放射化法による測定結果と3%以内で一致した。
小川 達彦
Annals of Nuclear Energy, 216, p.111256_1 - 111256_12, 2025/06
被引用回数：0AmBeに代表されるような、アルファ線を放出するアクチノイドと軽元素から構成される複合中性子源の性能を、多面的かつ絶対値でシミュレートできる新しい方法を開発した。この手法はモンテカルロ放射線輸送コードPHITSを用いて、
粒子用JENDL-5断面データライブラリと、ATIMAによる阻止能計算、ラティスジオメトリ計算機能を組み合わせることで、複合線源の様々な観測量を再現できる。従来の複合中性子源シミュレーション法は、測定値や近似的な断面積、多数イベントを平均した積分量などを使用するため、線源の特性パラメータ(アクチノイド粒径や線源のマクロなサイズ)を変更できなかったり、イベント内に放出される粒子間の相関を考慮できないなどの問題があったが、本手法ではそれらの問題が克服されている。この手法を用いて複合中性子源の様々な量、アクチノイド粒径依存性、中性子放出絶対強度、中性子および光子のエネルギースペクトル、中性子多重度、中性子放出の時間構造などを予測できることを示した。特に中性子の放出絶対強度とエネルギースペクトルについては、
AmBe、AmF、AmB、
PuBe、
CmBe、
Pu
Cの6種類の線源で実験値と比較し、実験値の絶対値に不確定性があるPuBeを除いて矛盾ない結果が得られた。この方法を用いれば、複合中性子源から生じる中性子のシミュレーションにおいて、カウントレート、コインシデンスイベントの量、ガンマ線によるノイズなど、実用的な定量指標の計算が可能になる。
青山 高士; 上野 文義; 佐藤 智徳; 加藤 千明; 佐野 成人; 山下 直輝; 大谷 恭平; 五十嵐 誉廣
Annals of Nuclear Energy, 214, p.111229_1 - 111229_6, 2025/05
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)To elucidate the effect of dissolved radionuclides on corrosion of carbon steels and on formation of corrosion products of carbon steel, corrosion tests and imaging plate analysis were conducted. Carbon steel samples immersed in 10 mM NaCl containing 
Sr and 
Cs were analyzed using an imaging plate. As a result, the distribution of Sr or Cs in the corrosion products formed on carbon steel was successfully visualized. Furthermore, the radioactivity of the corroded specimens was calculated from calibration curves prepared using a Sr standard.
Li, X.; 山路 哲史*; 佐藤 一憲*; 山下 拓哉
Annals of Nuclear Energy, 214, p.111217_1 - 111217_13, 2025/05
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)The decommissioning of Fukushima Daiichi NPP Unit-2 requires understanding of reactor damage and fuel debris distribution for effective debris retrieval. This study numerically analyzes potential Reactor Pressure Vessel (RPV) boundary failure due to eutectic melting of Control Rod Drive (CRD) housings during reheating after debris bed dryout. The Moving Particle Semi-implicit (MPS) method, with an enthalpy-based temperature algorithm and Boussinesq approximation, is applied to simulate melt/solid interactions in a 2-D model of the lower plenum. The CRD housing melting temperature is set at 1523 K based on a quasi-binary phase diagram of 304 Stainless Steel (SS) and Zirconium (Zr) and ELSA experiments. Results suggest local RPV failure at CRD housings, leading to melt release and refreezing. The estimated failure occurs 8-12 hours post-dryout (ca. 12:00-16:00 on 3/15/2011), providing insights into melt progression and boundary breach scenarios in Unit-2.
大塚 直彦*; 多田 健一; Cabellos, O.*; 岩本 修
Annals of Nuclear Energy, 212, p.110977_1 - 110977_9, 2025/03
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)ロスアラモス国立研究所のLANCEにて測定された新しい測定値を考慮して、3keVから1MeVの範囲のU-233の中性子捕獲断面積を評価した。得られた捕獲断面積は、JENDL-5の捕獲断面積よりも系統的に小さく、20keV付近では50%近く減少することとなった。新しく評価された断面積の妥当性を確認するため、ISCBEPハンドブックから選択された166の臨界実験を対象に、U-233の中性子捕獲断面積について、JENDL-5のデータを新たに評価した値に置き換えた上で、モンテカルロ中性子輸送計算を実行した。その結果、新しく評価された捕獲断面積の採用により、JENDL-5のU-233を用いた場合に比べてカイ2乗値がわずかに改善することが分かった。
elik, Y.*; Stankovskiy, A.*; 岩元 大樹; 岩元 洋介; Van den Eynde, G.*
Annals of Nuclear Energy, 212, p.111048_1 - 111048_12, 2025/03
被引用回数：1 パーセンタイル：68.64(Nuclear Science & Technology)The MCNP, PHITS, and FLUKA are general-purpose Monte Carlo radiation transport codes that are widely used for many real-world shielding problems at accelerator facilities around the world. For high beam energy and high beam current accelerator applications, neutron emission through the vacuum pipe along the reverse direction of incident proton beam is an important factor for a shielding design in order to correctly assess the dose rates for workers and the structural materials of the accelerator and handle with the waste activated by the backscattered neutron fluxes. In this work, neutron-production cross sections and thick target yield predictions from MC codes relying on physics models and nuclear data libraries are benchmarked against the experimental data, in order to assess their accuracy in predicting neutron emission and furthermore to assess the corresponding impact on shielding design. The results of this study demonstrate that the nuclear data libraries and physics models, which are not expected to give good results at lower energies (
MeV) but are used anyhow when there is no nuclear data available or above the energy range where the data tables end in the so-called "mix-and-match" strategy, need further improvements. Among the investigated proton induced nuclear data libraries, JENDL-4.0/HE produces the most satisfactory agreement to experimental data for all target materials, but may still benefit from refinement. Concerning the physics models of the codes, FLUKA V4-4.0 has the best performance in terms of reproducibility of the experimental values. It is also shown that all discrepancies between the calculations and the experiments for the energy range 
MeV are up to factor of two. This might be considered as an acceptable figure as it is equivalent to a normal safety margin (
) considered in shielding calculations of accelerator facilities around the world.
Lee, J.; Rossi, F.; 児玉 有; 弘中 浩太; 小泉 光生; 佐野 忠史*; 松尾 泰典*; 堀 順一*
Annals of Nuclear Energy, 211, p.111017_1 - 111017_7, 2025/02
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)Silica glass has been used as a base and host material in vitrified radioactive waste and lithium glass scintillator for neutron detection because of its superb transparency, high heat resistance, and excellent chemical inertness. Therefore, an accurate total cross section of the silica glass is important to evaluate the criticality safety for the vitrified wastes and to understand the neutron response for lithium glass scintillators accurately. In the present study, to provide the accurate total cross section in the thermal and epithermal energy range, the neutron transmission measurements were carried out by a pulsed neutron beam with the time-of-flight method at the Kyoto University Institute for Integrated Radiation and Nuclear Science - Linear Accelerator. We obtained the neutron total cross section of the silica glass in the energy region from 0.002 eV to 25 eV. The obtained results were compared and discussed with the previous results and the evaluated data.
Brumm, S.*; Gabrielli, F.*; Sanchez Espinoza, V.*; Stakhanova, A.*; Groudev, P.*; Petrova, P.*; Vryashkova, P.*; Ou, P.*; Zhang, W.*; Malkhasyan, A.*; et al.
Annals of Nuclear Energy, 211, p.110962_1 - 110962_16, 2025/02
被引用回数：3 パーセンタイル：86.32(Nuclear Science & Technology)The completed Horizon-2020 project on "Management and Uncertainties of Severe Accidents (MUSA)" has reviewed uncertainty sources and Uncertainty Quantification methodology for the purpose of assessing Severe Accidents (SA). The key motivation of the project has been to bring the advantages of the Best Estimate Plus Uncertainty approach to the field of Severe Accident. The applications brought together a large group of participants that set out to apply uncertainty analysis (UA) within their field of SA modelling expertise, in particular reactor types, but also SA code used (ASTEC, MELCOR, etc.), uncertainty quantification tools used (DAKOTA, RAVEN, etc.), detailed accident scenarios, and in some cases SAM actions. This paper synthesizes the reactor-application work at the end of the project. Analyses of 23 partners are sorted into different categories, depending on whether their main goal is/are (i) uncertainty bands of simulation results; (ii) the understanding of dominating uncertainties in specific sub-models of the SA code; (iii) improving the understanding of specific accident scenarios, with or without the application of SAM actions; or, (iv) a demonstration of the tools used and developed, and of the capability to carry out an uncertainty analysis in the presence of the challenges faced. The partners' experiences made during the project have been evaluated and are presented as good practice recommendations. The paper ends with conclusions on the level of readiness of UA in SA modelling, on the determination of governing uncertainties, and on the analysis of SAM actions.
郡司 智; 荒木 祥平; 井澤 一彦; 須山 賢也
Annals of Nuclear Energy, 209, p.110783_1 - 110783_7, 2024/12
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)燃料デブリの組成や性状は不確実であるため、臨界安全性評価に使用される計算コードや核データを検証するには臨界実験が必要である。このため、原子力機構は、臨界集合体STACYの改造を行っている。STACY更新炉の初臨界は2024年春に予定されている。本稿では、STACY更新炉の初臨界時の基本炉心構成仕様の特性について事前解析の結果を報告する。初臨界時には中性子減速条件の異なる2種類の格子板(格子間隔は1.50cmと1.27cm)を用意される。一方で、利用可能なUO
燃料棒の数には制限がある。これらの実験的制約を満たす最初の臨界のための炉心構成は、計算解析によって設計された。最適減速条件に近い1.50cmピッチの格子板を備えた円柱形の炉心構成では、臨界に達するには253本の燃料棒が必要となる。1.27cmピッチの格子板については、ピッチを2倍にして2.54cmピッチの炉心構成を検討した。この場合、臨界に達するには213本の燃料棒が必要となる。さらに、燃料デブリの様態をシミュレートするために、鉄またはコンクリート模擬棒を使用した実験炉心構成についても検討した。本稿では、これらの炉心構成と炉心特性の解析結果を示す。
福田 航大
Annals of Nuclear Energy, 208(1), p.110748_1 - 110748_10, 2024/12
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)The Windscale Works criticality accident in 1970 resulted from mixing an aqueous solution with an organic solvent with different plutonium concentrations and densities. Although this accident has been studied using improved computer capabilities in recent years, a precise criticality scenario has not yet been identified. This study aims to clarify a possible criticality scenario of the accident-the time variation of reactivity and its mechanism. The accident was simulated by combining the multiphase computational fluid dynamics solver of OpenFOAM and the delta-tracking-based Monte Carlo neutron transport code Serpent2. Consequently, the periodic uneven arrangement of fluids might have caused oscillations in neutron leakage and absorption, resulting in periodic wavy reactivity changes. Furthermore, the emulsion, which was thought to be the primary cause, might not be the dominant mechanism for reactivity change, although it contributed to the criticality of the accident.
今泉 悠也; 神山 健司; 松場 賢一
Annals of Nuclear Energy, 206, p.110658_1 - 110658_10, 2024/10
被引用回数：1 パーセンタイル：68.64(Nuclear Science & Technology)In severe accidents of SFRs, molten core materials can discharge from the core, and the jet can impinge on the lower structure plate. After the jet impingement, fragmented discharged materials can form ring-shape solidification. A fundamental experiment was conducted to simulate the behavior. In order to simulate the behavior of solid body creation and motion, a new solid body formation model by inter-particle attraction force in particle method was developed. The advantage of the new model is that it can simulate creation, formation, and motion of solid bodies without any artificial treatment as solid bodies. The movable solid bodies by the new model have any size, shape, and number, and they are created and grown by solidification, and diminish and disappear by melting. The mechanism based on the inter-particle attraction force is common with that in real world where interatomic attraction force is the cause of solid body formation.
丸山 修平; 山本 章夫*; 遠藤 知弘*
Annals of Nuclear Energy, 205, p.110591_1 - 110591_13, 2024/09
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)This study developed a new method for evaluating the uncertainty in reactor core/shielding characteristics attributable to the scattering angle distribution, employing a random sampling (RS) technique integrated with continuous energy Monte Carlo (CEMC) calculations. The impact of neutron scattering angle is not negligible in the analysis of fast reactor cores and shielding. Recent advancements have enabled the high-accuracy assessment of nuclear data-induced uncertainty by merging CEMC calculations and the RS technique. Nonetheless, a method to quantify uncertainty due to scattering angle distribution remains unestablished. This study introduces a new approach for uncertainty quantification related to scattering angle distribution in CEMC-RS, utilizing the maximum entropy method. The effectiveness of this method was verified through comparison with results from the classical deterministic uncertainty quantification approach based on generalized perturbation theory. Overall, this method offers a more accurate tool for nuclear engineers and researchers in evaluating and managing uncertainties in reactor design and safety analysis.
安部 諭; 柴本 泰照
Annals of Nuclear Energy, 202, p.110461_1 - 110461_16, 2024/07
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)During a severe accident in a nuclear containment vessel, jets released from the primary system exhibit complex thermohydraulic behavior due to buoyancy effects and impingement on internal obstacles such as inner walls and floors. Thus, the obstacle-influenced jets are of interest in recent research activities. This paper describes an experimental investigation of the behavior of jets passing through a grid-type obstacle. The flow field was acquired by a particle image velocimetry system. The experiment captured the jet fragmentation by the grid-type obstacle and their recoupling. The mean velocity field obtained by postprocessing indicates a "Rectifying effect," with the axial velocity increasing at the center and the magnitude of the radial velocity decreasing. The meandering flow was suppressed due to this effect. In the near grid-obstacle region, the axial turbulence intensity was relatively large at the edge of each fragmented region due to shear stress. Moreover, the spatial distribution of the radial turbulence fluctuation became more complex. Further investigation is required to clarify the budget of the transport equation for turbulence fluctuation. The experimental data shown in this paper is useful for computational fluid dynamics validation.
嶋田 和真; 櫻原 達也*; Farshadmanesh, P.*; Reihani, S.*; Mohagehgh, Z.*
Annals of Nuclear Energy, 197, p.110243_1 - 110243_12, 2024/03
被引用回数：1 パーセンタイル：30.19(Nuclear Science & Technology)本研究は、原子力発電所に対するレベル3確率論的リスク評価(PRA)において住民の避難行動を設定する際の主観的な専門家判断を回避したレベル3PRA手法を開発する。そのために、交通シミュレーションコードMATSimで出力した避難速度をレベル3PRAコードMACCSに入力した。さらに、道路封鎖を検討する箇所の優先順位を設定するために、自然災害リスク評価コードHAZUSを用いて地震による道路封鎖リスクを評価する手法を開発した。そして、米国原子力規制委員会が実施した最先端の原子炉事故影響研究において採用されたSequoyah原子力発電所のケーススタディに対して、住民の避難経路と放射線被ばく線量の関係を評価した。その結果、地震封鎖リスクは小さいが、封鎖されると住民の被ばく線量が増加する避難経路を見出した。この結果は、提案するレベル3PRA手法が避難経路を強化する意思決定を支援することを示した。
谷口 良徳; 三原 武; 垣内 一雄; 宇田川 豊
Annals of Nuclear Energy, 195, p.110144_1 - 110144_11, 2024/01
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)A reactivity-initiated accident (RIA)-simulated test CN-1 on a high-burnup 64 GWd/t mixed-oxide fuel rod sheathed with M5
cladding was conducted at the Nuclear Safety Research Reactor, resulting in fuel failure. A small opening with slight ballooning deformation characterized the post-test visual appearance of the test fuel rod. Simulation using fuel performance codes FEMAXI-8/RANNS predicted rod survival under early phase loading induced by pellet-cladding mechanical interaction and subsequent boiling transition, and the cladding surface temperature measured online confirmed the occurrence of boiling transition. The experimental observation and simulation indicate that the failure was caused by a high-temperature rupture following increased rod-internal pressure. The RANNS sensitivity analysis revealed that a mechanical state parameter dedicated to predicting plastic instability might be an effective index for evaluating the risk of rupture failure during RIAs.
浜瀬 枝里菜; 大釜 和也; 河村 拓己*; 堂田 哲広; 田中 正暁; 山野 秀将
Annals of Nuclear Energy, 195, p.110157_1 - 110157_14, 2024/01
被引用回数：2 パーセンタイル：30.19(Nuclear Science & Technology)高速炉プラント動特性解析コードSuper-COPDのスクラム不作動流量喪失事象に対する妥当性確認のため、FFTFの受動的安全性試験LOFWOS No.13試験を対象としたIAEAベンチマークに参加した。ブラインドフェーズで課題として抽出された燃料集合体出口温度及び全反応度の評価精度向上のため、集合体間熱移行及び集合体間ギャップ部流れを考慮した全炉心モデル及び炉心湾曲反応度簡易評価モデルを導入した。最終フェーズ解析の結果、2次ピーク時の集合体出口温度を良好に再現するとともに、全反応度の実測値の挙動を概ね評価できたことから、LOFWOSに対するSuper-COPDの妥当性を確認した。
Sahboun, N. F.; 松本 俊慶; 岩澤 譲; Wang, Z.; 杉山 智之
Annals of Nuclear Energy, 195, p.110145_1 - 110145_12, 2024/01
被引用回数：2 パーセンタイル：30.19(Nuclear Science & Technology)Relocated corium into the Primary Containment Vessel needs to be properly cooled to avoid or mitigate molten core concrete interactions in the PCV in order to maintain its supporting capability for the reactor pressure vessel and to suppress combustible or non-condensable gas releases. To know how effective the cooling is, it became important to know the geometry of the relocated corium. The present study chooses to focus on the "Wet Cavity" strategy and to build a reliable tool to evaluate the corium coolability in such a case. To achieve this goal, a previously developed formulation built to predict the corium geometry under the "Dry Cavity" strategy was extended to the conditions used in the "Wet Cavity" strategy. This extension includes the effect of solidification and cooling from the water by using a newly developed expression for the dimensionless thickness s, the water subcooling, and the melts super heat. After the validation of the extended formulation was confirmed, potential restrictions and limitations were investigated.
今泉 悠也; 青柳 光裕; 神山 健司; 松場 賢一; Akaev, A.*; Mikisha, A.*; Baklanov, V.*; Vurim, A.*
Annals of Nuclear Energy, 194, p.110107_1 - 110107_11, 2023/12
被引用回数：2 パーセンタイル：51.90(Nuclear Science & Technology)In severe accidents of SFRs, the cooling of the residual core materials, which is called "in-place cooling", is one of the important factors for In-Vessel Retention (IVR). For the evaluation, behavior of the in-place cooling was analyzed by the SIMMER-III code. In order to understand the in-place cooling, method of Phenomena Identification and Ranking Table (PIRT) was applied. Based on the result, an out-of-pile experiment which focused on the extracted factors was conducted. In the experiment, continuous oscillation of sodium level was observed by sodium vaporization and condensation. Analysis for the out-of-pile experiment was conducted by SIMMER-III, but the results were different between the experiment and the analysis. By investigation of the analysis result, it was revealed that the difference was due to occupation of non-condensable gas. Therefore, an analysis model of inter-cell gas mixing was newly developed, and the agreement was significantly improved by the new model.
