Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
福田 航大
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.
寺田 敦彦; 永石 隆二
Journal of Nuclear Science and Technology, 61(8), p.1135 - 1154, 2024/08
被引用回数：0 パーセンタイル：0.00(Nuclear Science & Technology)天井と側面に換気口を設けた部屋の床面から漏洩した水素の挙動について、CFD解析を用いて実験結果との照合を行うことで、部屋周囲の風況による影響を明らかにした。室内の水素濃度の経時変化は、天井と床面付近の間で成層化がみられ、計算結果については局所的な流れ場の特性変化を反映できる乱流モデルと比較的大きな乱流シュミット数(1.0)を組み合わせた場合に実験結果の再現性が向上した。風況によって、室内に水素が滞留したケースについて、天井開口部にシンプルな板を設置することで換気の改善が図れることを示した。
Hamdani, A.; 安部 諭; 石垣 将宏; 柴本 泰照; 与能本 泰介
Progress in Nuclear Energy, 153, p.104415_1 - 104415_16, 2022/11
被引用回数：3 パーセンタイル：52.93(Nuclear Science & Technology)This paper describes the computational fluid dynamics (CFD) analysis and validation works from the previous experimental study on the natural convection driven by outer surface cooling in the presence of density stratification consisting of air and helium (as a mimic gas of hydrogen). The experiment was conducted in the Containment InteGral effects Measurement Apparatus (CIGMA) facility at Japan Atomic Energy Agency (JAEA). The numerical simulation was carried out to analyze the detailed effect of the cooling region on the erosion of the helium stratification layer. The temporal and spatial evolution of the helium concentration and the gas temperature inside the containment vessel was predicted and validated against the experimental data. In addition, two stratification behaviors that depend on the cooling location were presented and discussed. The CFD simulation confirmed that an upper head cooling caused two counter-rotating vortexes in the helium-rich zone. Meanwhile, the upper half body cooling caused two counter-rotating vortexes in the helium-poor zone. These findings are important to understand the mechanism of the density stratification process driven by natural convection in the containment vessel.
中村 秀夫; Bentaib, A.*; Herranz, L. E.*; Ruyer, P.*; Mascari, F.*; Jacquemain, D.*; Adorni, M.*
Proceedings of International Conference on Topical Issues in Nuclear Installation Safety; Strengthening Safety of Evolutionary and Innovative Reactor Designs (TIC 2022) (Internet), 10 Pages, 2022/10
The WGAMA activity achievements have been published as technical reports, becoming reference materials to discuss innovative methods, materials and technologies in the fields of thermal-hydraulics, computational fluid dynamics (CFD) and severe accidents (SAs). The International Standard Problems (ISPs) and Benchmarks of computer codes have been supported by a huge amount of the databases for the code validation necessary for the reactor safety assessment with accuracy. The paper aims to review and summarize the recent WGAMA outcomes with focus on new advanced reactor applications including small modular reactors (SMRs). Particularly, discussed are applicability of major outcomes in the relevant subjects of passive system, modelling innovation in CFD, severe accident management (SAM) countermeasures, advanced measurement methods and instrumentation, and modelling robustness of safety analysis codes. Although large portions of the outcomes are considered applicable, design-specific subjects may need careful considerations when applied. The WGAMA efforts, experiences and achievements for the safety assessment of operating nuclear power plants including SA will be of great help for the continuous safety improvements required for the advanced reactors including SMRs.
岡垣 百合亜; 与能本 泰介; 石垣 将宏; 廣瀬 意育
Fluids (Internet), 6(2), p.80_1 - 80_17, 2021/02
Many thermohydraulic issues about the safety of light water reactors are related to complicated two-phase flow phenomena. In these phenomena, computational fluid dynamics (CFD) analysis using the volume of fluid (VOF) method causes numerical diffusion generated by the first-order upwind scheme used in the convection term of the volume fraction equation. Thus, in this study, we focused on an interface compression (IC) method for such a VOF approach; this technique prevents numerical diffusion issues and maintains boundedness and conservation with negative diffusion. First, on a sufficiently high mesh resolution and without the IC method, the validation process was considered by comparing the amplitude growth of the interfacial wave between a two-dimensional gas sheet and a quiescent liquid using the linear theory. The disturbance growth rates were consistent with the linear theory, and the validation process was considered appropriate. Then, this validation process confirmed the effects of the IC method on numerical diffusion, and we derived the optimum value of the IC coefficient, which is the parameter that controls the numerical diffusion.
齋藤 慎平*; De Rosis, A.*; Fei, L.*; Luo, K. H.*; 海老原 健一; 金子 暁子*; 阿部 豊*
Physics of Fluids, 33(2), p.023307_1 - 023307_21, 2021/02
被引用回数：40 パーセンタイル：98.33(Mechanics)流れ場で沸騰が発生する現象は強制対流沸騰として知られている。今回、飽和条件の流れ場中の円柱上での沸騰システムを数値的に調査した。複雑な気液相変化現象に対処するために、擬ポテンシャル格子ボルツマン法(LBM)に基づく数値スキームを開発した。高いレイノルズ数の数値安定性を高めるため、衝突項を中心モーメント(CM)の空間で解いた。CMベースのLBMに適した力場スキームを採用することで簡潔でありながら堅牢なアルゴリズムとなっている。さらに、熱力学的一貫性を確保するために必要な追加項をCMの枠組みにおいて導出した。現在のスキームの有効性は、核形成,成長、および30-30000の間で変化するレイノルズ数の蒸気泡の離脱を含む一連の沸騰プロセスに対してテストされた。開発したCMベースのLBMは、初期気相などの人工的な入力なしに、核沸騰,遷移沸騰、および膜沸騰の沸騰様式を再現できる。結果からプール沸騰ではなく強制対流システムでも抜山曲線として知られる典型的な沸騰曲線が現れることが分かった。また、今回のシミュレーションは膜沸騰領域でも断続的な直接固液接触の実験的観察を支持することが分かった。
Herranz, L. E.*; Jacquemain, D.*; Nitheanandan, T.*; Sandberg, N.*; Barr
, F.*; Bechta, S.*; Choi, K.-Y.*; D'Auria, F.*; Lee, R.*; 中村 秀夫
Progress in Nuclear Energy, 127, p.103432_1 - 103432_14, 2020/09
被引用回数：4 パーセンタイル：19.21(Nuclear Science & Technology)WGAMA started on Dec. 31st 1999 to assess and strengthen the technical basis needed for the prevention, mitigation and management of potential accidents in NPP and to facilitate international convergence on safety issues and AM analyses and strategies. WGAMA addresses reactor thermal-hydraulics (Thys), in-vessel behavior of degraded cores, containment behavior and protection, and FP release, transport, deposition and retention, for both current and advanced reactors. This paper summarizes such WGAMA contributions in Thys, CFD and severe accidents, which include the Fukushima-Daiichi accident impacts on the WGAMA activities and their substantial outcomes. Around 50 technical reports have become reference in the related fields, which appear in References. Recommendations in these reports include further research, some of which have given rise to the joint projects conducted or underway within the OECD framework. Ongoing WGAMA activities are numerous and a number of them are to be launched in the near future, which are shortly mentioned too.
Trianti, N.; 茂木 孝介; 杉山 智之; 丸山 結
Proceedings of 2020 International Conference on Nuclear Engineering (ICONE 2020) (Internet), 9 Pages, 2020/08
The computational fluid dynamics (CFD) have been developed to analyze the correlation equation for laminar flame speed of hydrogen-air mixtures. This analysis was carried out on the combustion of hydrogen-air mixtures performed at the spherical bomb experiment facility consists of a spherical vessel equipped (563 mm internal diameter). The facility has been designed and built at CNRS-ICARE laboratory. The simulation was carried out using the reactingFoam solver, one of a transient chemical reaction solver in OpenFOAM 5.0. The LaunderSharmaKE model was applied for turbulent flow. The interaction of the chemical reaction with the turbulent flow was taken into account using PaSR (Partial Stirred Reactor) model with 19 elementary reactions for the hydrogen combustion. The initial condition of spherical flame analysis was set so as to be consistent with those of the experiment. The position of the flame front was detected by the steep drop of hydrogen mass fraction in the spherical radii, and the flame propagation velocity was estimated from the time-position relationship. The analysis result showed the characteristic of spherical flame acceleration was qualitatively reproduced even though it has a discrepancy with the experiment. After validating the calculation of spherical experiments, a laminar burning velocity correlation is presented using the same boundary conditions with the variation of hydrogen concentration, temperature, and pressure. The calculation of laminar flame speed of hydrogen-air mixtures by reactingFoam use reference temperature T
= 293 K and reference pressure P = 1 atm with validated in the range of hydrogen concentration 6-20%; range of temperature 293-493 K; and range of pressure 1-3 atm.
内堀 昭寛; 渡部 晃*; 高田 孝; 大島 宏之
Nuclear Technology, 205(1-2), p.119 - 127, 2019/01
被引用回数：3 パーセンタイル：28.63(Nuclear Science & Technology)ナトリウム冷却高速炉の蒸気発生器におけるナトリウム-水反応現象を評価するため、圧縮性多成分多相流及びNa-水化学反応を対象とした数値解析コードSERAPHIMを開発している。従来のSERAPHIMコードは差分法を用いているが、本研究では、伝熱管の存在する複雑形状領域に対して解析精度を向上することを目的に非構造格子に対応した解析手法を開発した。解析手法妥当性確認の一環として不足膨張噴流実験の解析を実施した結果、解析結果における圧力分布やマッハディスクの形成位置が実験結果と一致する結果を得た。また、Na中へ水蒸気が噴出する現象を対象とした解析を実施し、実現象に対する適用性を確認した。
小野 綾子; 鈴木 貴行*; 吉田 啓之
Proceedings of 26th International Conference on Nuclear Engineering (ICONE-26) (Internet), 6 Pages, 2018/07
液膜蒸発モデルに基づいた限界熱流束予測評価手法の開発において必要な伝熱面ごく近傍の気泡の挙動や気泡同士の接合現象について、詳細二相流解析コードTPFITによる発泡,気泡生長および接合をシミュレートすることで知見を得る。本発表では、高圧域での沸騰現象の理解と、気泡形状の影響の理解のために、圧力をパラメタとした解析結果に基づいて、気泡の挙動に圧力が与える影響について報告する。
内堀 昭寛; 渡部 晃*; 高田 孝; 大島 宏之
Journal of Nuclear Science and Technology, 54(10), p.1036 - 1045, 2017/10
被引用回数：5 パーセンタイル：41.83(Nuclear Science & Technology)ナトリウム冷却高速炉の蒸気発生器におけるナトリウム-水反応現象を評価するため、圧縮性多成分多相流及びNa-水化学反応を対象とした数値解析コードSERAPHIMを開発している。従来のSERAPHIMコードは差分法を用いているが、本研究では、伝熱管の存在する複雑形状領域に対して解析精度を向上することを目的に非構造格子に対応した解析手法を開発した。解析手法妥当性確認の一環として不足膨張噴流実験の解析を実施した結果、解析結果における圧力分布が実験結果と一致する結果を得た。また、Na中へ水蒸気が噴出する現象を対象とした解析を実施し、実現象に対する適用性を確認した。構造格子及び非構造格子を用いた解析から、非構造格子を適用することの効果についても確認した。
内堀 昭寛; 高田 孝; 大島 宏之; 渡部 晃*
Proceedings of 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-17) (USB Flash Drive), 12 Pages, 2017/09
ナトリウム冷却高速炉の蒸気発生器におけるナトリウム-水反応現象を評価するため、圧縮性多成分多相流及びNa-水化学反応を対象とした数値解析コードSERAPHIMを開発している。従来のSERAPHIMコードは差分法を用いているが、本研究では、伝熱管の存在する複雑形状領域に対して解析精度を向上することを目的に非構造格子に対応した解析手法を開発した。解析手法妥当性確認の一環として不足膨張噴流実験の解析を実施した結果、解析結果における圧力分布やマッハディスクの形成位置が実験結果と一致する結果を得た。また、Na中へ水蒸気が噴出する現象を対象とした解析を実施し、実現象に対する適用性を確認した。
Guo, L.*; 守田 幸路*; 飛田 吉春
Journal of Nuclear Science and Technology, 53(2), p.271 - 280, 2016/02
被引用回数：9 パーセンタイル：62.19(Nuclear Science & Technology)For the analysis of debris behavior in core disruptive accidents of liquid metal fast reactors, a hybrid computational tool was developed using the discrete element method (DEM) for calculation of solid particle dynamics and a multi-fluid model of a reactor safety analysis code, SIMMER-III, to reasonably simulate transient behavior of three-phase flows of gas-liquid-particle mixtures. A coupling numerical algorithm was developed to combine the DEM and fluid-dynamic calculations, which are based on an explicit and a semi-implicit method, respectively. The developed method was validated based on experiments of water-particle dam break and fluidized bed in systems of gas-liquid-particle flows. Reasonable agreements between the simulation results and experimental data demonstrate the validity of the present method for complicated three-phase flows with large amounts of solid particles.
加藤 由幹; 吉田 啓之; 横山 諒太郎*; 金川 哲也*; 金子 暁子*; 文字 秀明*; 阿部 豊*
Proceedings of 23rd International Conference on Nuclear Engineering (ICONE-23) (DVD-ROM), 8 Pages, 2015/05
In this study, visualization experiment for the bubbly flow in a horizontal pipe excited by oscillation acceleration was performed to understand bubbly flow behavior under earthquake acceleration. Liquid pressure was also measured at upstream and downstream of the test section. In addition, to consider detailed effects of pressure gradient on bubble motion, numerical simulation of two-phase flow in horizontal pipe with vibration was performed by a detailed two-phase flow simulation code with an advanced interface tracking method: TPFIT. Based on observed images and calculated results, bubble velocity was evaluated. It was confirmed that the pressure gradient amplitude increased with the increase of the frequency of the table. In addition, it was observed that the bubble velocity amplitude also increases with the increase of the frequency of the table. It was concluded that the bubble motion was strongly affected by the pressure gradient in the test section.
吉田 啓之; 永武 拓; 高瀬 和之; 金子 暁子*; 文字 秀明*; 阿部 豊*
Mechanical Engineering Journal (Internet), 1(4), p.TEP0025_1 - TEP0025_11, 2014/08
In this study, to develop the predictive technology of two-phase flow dynamics under earthquake acceleration, a detailed two-phase flow simulation code with an advanced interface tracking method TPFIT was expanded. In addition, the bubbly flow in a horizontal pipe excited by oscillation acceleration and under the fluctuation of the liquid flow was simulated by using the modified TPFIT. In the results, it was confirmed that the modified TPFIT can predict time dependent velocity distribution around the bubbles and shapes of bubbles qualitatively. The main cause of bubble deformation observed is large shear stress at the lower part of the bubble, and this large shear stress is induced by the velocity difference between the liquid phase and bubble. Moreover, we discussed about the difference between both effects of flow rate fluctuation and structure vibration on two-phase flow. In the results, bubble acceleration of the structure vibration case was larger than that of the flow rate fluctuation case. Finally, it was concluded that unsteady shear stress induced by vibration of the pipe wall was one of the main driving forces of bubble motion in structure vibration case.
海老原 健一
JAERI-Research 2005-004, 121 Pages, 2005/03
JAERI-Research-2005-004.pdf:19.79MB
本報告書は、まず、格子ガス法及び格子ボルツマン法の二相流体モデルの二相流シミュレーションへの適用妥当性及び有用性を検討した。その結果に基づき、格子ボルツマン法の二相流体モデルであるHCZモデルを、基本的かつ重要な二相流流動様式の1つである水平層状二相流のシミュレーションへ適用した。その結果、HCZモデルの界面は、Kelvin-Helmholtz不安定性を満足することが確認された。またTaitelとDuklerの理論的流動様式線図(T-D線図)の界面成長に関する曲線が再現された。さらに、幅が狭い流路では、流れの三次元性の影響が顕著となり、T-D線図が与える以上の低密度相のみかけ流速が界面成長に必要であることが明らかとなった。次に、より複雑な界面現象を伴う液滴発生のシミュレーションを行い、Ishii-Grolmesの実験相関式を、実験データの分布範囲において、再現できることがわかった。
海老原 健一; 渡辺 正
Proceedings of 2004 ASME International Mechanical Engineering Congress and Exposition (IMECE '04) (CD-ROM), 9 Pages, 2004/11
本論文は、水平矩形断面流路中における成層波状二相流の界面成長の数値シミュレーションについて記述する。界面波の成長への流路幅の影響を、異なる流路幅に対するシミュレーションを行うことによって評価した。本論文で採用した数値解析モデルは、自発的な界面現象をシミュレーションすることができる、格子ボルツマン法の一成分二相流体モデルである。シミュレーションでは、波の成長を観察し、その波の成長と流動状態を特徴づける無次元数との関係を、TaitelとDuklerによって提案された流動様式線図と比較した。流路幅が広い場合には、計算結果の関係がTaitel-Duklerの流動様式線図とほぼ一致することが確認された。流路幅が狭くなると、界面が成長のために、より多くの低密度相の質量流量が必要になり、流動様式線図の関係からのずれが生ずることが示された。
海老原 健一
筑波大学大学院システム情報工学研究科博士学位論文, 134 Pages, 2004/09
本論文では、格子法の気液モデルの有用性を検討し、さらに格子ボルツマン法気液モデル(HCZモデル)を、格子ボルツマン法の適用例が見られていない水平層状二相流の三次元シミュレーションへ適用した。その結果、以下のことが明らかになった。(1)HCZモデルによってシミュレーションされる二相界面が、Kelvin-Helmholtz不安定性理論を、二次元及び三次元の場合について満足することがわかった。(2)矩形断面流路中の界面成長シミュレーションにおいて、界面成長と流動状態の関係が、理論的考察によって得られるTaitelとDuklerの流動様式線図を再現できることがわかった。(3)流路幅が異なる場合の界面成長シミュレーションによって、流路幅が狭くなった時、流れの三次元性の影響が顕著となり、界面成長に必要な低密度相の流量が、理論的流動様式線図から得られる値より大きくなることが明らかとなった。(4)界面成長より複雑な界面現象を示す液滴発生シミュレーションにおいて、液滴発生とその時の流動状態の関係が、IshiiとGrolmesによって提案された実験相関式を、ほぼ再現できることがわかった。
服部 博文*; 佐藤 博; 長野 靖尚*
日本機械学会論文集,B, 70(696), p.1919 - 1926, 2004/08
衝突噴流はよどみ点近傍で高い熱及び物質伝達率が得られるため工業上の広い分野で多用されている。また、衝突噴流は物質・生命科学実験施設に設置される減速材容器の冷却に使用される。そのため、その流動及び伝熱特性を明らかにすることは各種機器の性能を向上させるために重要である。本研究では、平面衝突噴流熱伝達場に着目し、高精度差分法による直接数値シミュレーションを行った。熱伝達場に最も影響する衝突距離をパラメータとして計算を行い、特に局所熱伝達場の第2極大値の出現に対する熱輸送機構を精査した。
海老原 健一; 渡辺 正
日本機械学会論文集,B, 70(694), p.1393 - 1399, 2004/06
一成分二相格子ボルツマン法によって、その幅が高さより狭い矩形断面流路内における水平層状二相流のシミュレーションを行った。3種類の異なる流路幅の場合について、二相間の界面成長を測定した。さらに二相流動を特徴付ける無次元数を測定し、それらと界面成長との関係とTaitelとDuklerによって提案された流動様式線図とを比較した。流路幅が狭くなるに従い、シミュレーションにより得られた界面の成長・非成長の境界が大きくなることがわかった。
