塩津 弘之; 伊藤 裕人*; 石川 淳; 杉山 智之; 丸山 結
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 6 Pages, 2018/11
The VERDON-2 experiment for FPs transport in steam environment was analyzed with the mechanistic FPs transport code incorporating thermodynamic chemical equilibrium model in order to assess its predictive capability for transport behavior of key FPs, especially for highly volatile FPs such as Cs and I. The present analysis reproduced well the Cs deposition profile obtained from the experiment, which revealed that Cs was transported as CsOH in early phase of FP release from fuel, and then formed CsMoO after increasing Mo release. On the other hand, the deposition peak of I was predicted to appear at 720 K, which was significantly higher than the experimental result at 600 K. This discrepancy was potentially caused by the following two points: lack of the other stable species in thermodynamics database for thermodynamic chemical equilibrium model, or failure of chemical equilibrium assumption for iodide species.
Trianti, N.; 佐藤 允俊*; 杉山 智之; 丸山 結
Proceedings of 11th Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-11) (Internet), 7 Pages, 2018/11
Simulation techniques have been developed to analyze the deflagration behavior of hydrogen generated during a hypothetical severe accident in nuclear power plants. The CFD analysis was carried out on the hydrogen deflagration experiment performed at the ENACCEF2 facility composed mainly of a vertical cylindrical tube filled with hydrogen-air mixture and nine annular obstacles were placed in the lower part of the tube. The simulation was carried out by the reactingFoam solver of OpenFOAM 3.0, an open source software for the CFD analysis. The RNG (Renormalization group) k- model was applied for turbulent flow. The interaction of the chemical reaction with the turbulent flow was considered using PaSR (Partial Stirred Reactor) model with 19 elementary reactions for the hydrogen combustion. The analysis result showed the characteristic of flame acceleration by the obstacle region was qualitatively reproduced even though has discrepancy with the experiment.
伊藤 裕人*; 塩津 弘之; 田中 洋一*; 西原 慧径*; 杉山 智之; 丸山 結
JAEA-Data/Code 2018-012, 42 Pages, 2018/10
原子力施設事故時において施設内を移行する核分裂生成物(FP)の化学組成は、比較的遅い反応の影響を受けることにより化学平衡を仮定して評価した組成とは異なる場合が想定される。そのため、反応速度を考慮した化学組成評価が求められる。一方で、原子力施設事故時の複雑な反応に関する反応速度の知見は現状では限られており、実機解析に適用できるデータベースの構築に至っていない。そこで、FP化学組成評価における反応速度による不確かさの低減のため、化学平衡論及び反応速度論の部分混合モデルに基づく化学組成評価コードCHEMKEqを開発した。このモデルは、系全体の質量保存則の下、前駆平衡と見なせる化学種を化学平衡論モデルにより評価し、その後の比較的遅い反応を反応速度論モデルにより解くものである。さらにCHEMKEqは、本混合モデルに加え一般的な化学平衡論モデル及び反応速度論モデルが使用可能であり、かつ、それらモデル計算に必要なデータベースを外部ファイル形式とすることで汎用性の高い化学組成評価コードとなっている。本報は、CHEMKEqコードの使用手引書であり、モデル, 解法, コードの構成とその計算例を記す。また付録には、CHEMKEqコードを使用する上で必要な情報をまとめる。
Bentaib, A.*; Chaumeix, N.*; Grosseuvres, R.*; Bleyer, A.*; Gastaldo, L.*; Maas, L.*; Jallais, S.*; Vyazmina, E.*; Kudriakov, S.*; Studer, E.*; et al.
Proceedings of 12th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics, Operation and Safety (NUTHOS-12) (USB Flash Drive), 11 Pages, 2018/10
In the framework of the French MITHYGENE project, the new highly instrumented ENACCEF2 facility was built at the Institut de Combustion Aerothermique Reactivite et Environnement (ICARE) of the Centre National de la Recherche Scientifique (CNRS) in Orleans (France) to address the flame propagation in hydrogen combustion during a severe accident. The ENACCEF2 facility is a vertical tube of 7.65 m height and 0.23 m inner diameter. In the lower part of the tube, annular obstacles are installed to promote turbulent flame propagation. At the initiative of the MITHYGENE project consortium and the European Technical Safety Organisation Network (ETSON), a benchmark on hydrogen combustion was organised with the goal to identify the current level of the computational tools in the area of hydrogen combustion simulation under conditions typical for safety considerations for NPP. In the proposed paper, the simulation results obtained by participating organizations, using both Computational Fluid Dynamics (CFD) and lumped-parameter computer codes, are compared to experimental results and analysed.
玉置 等史; 石川 淳; 杉山 智之; 丸山 結
Proceedings of Asian Symposium on Risk Assessment and Management 2018 (ASRAM 2018) (USB Flash Drive), 6 Pages, 2018/10
石川 淳; Zheng, X.; 塩津 弘之; 杉山 智之; 丸山 結
Proceedings of Asian Symposium on Risk Assessment and Management 2018 (ASRAM 2018) (USB Flash Drive), 6 Pages, 2018/10
Japan Atomic Energy Agency is pursuing the development and application of the methodologies on fission product (FP) chemistry for source term analysis by using integrated severe accident analysis code THALES2/KICHE. Generally, specific chemical forms of iodine and cesium such as cesium iodide (CsI) and cesium hydroxide (CsOH) were assumed in the source term analysis for light water reactors using an integrated severe accident analysis code. The accident at the Fukushima Dai-ichi Nuclear Power Station leads possible chemical effects of BC control materials and atmosphere on chemical speciation of iodine and cesium such as cesium metaborate (CsBO) and hydrogen iodide (HI). The difference of chemical speciation affects not only the FP behavior in the reactor coolant system (RCS) and transport to containment but also pH value of the suppression pool water in the containment. The pH value is one of the influential factors on the release of gaseous iodine (I and organic iodine) from containment liquid phase. In the present study, the improvement of the THALES2/KICHE code in terms of FP chemistry in RCS was performed and applied to source term analysis for severe accidents at a boil water reactor with Mark-I containment vessel. This paper discusses the chemical speciation of iodine and cesium, and FP behavior and transport to containment.
Zheng, X.; 玉置 等史; 石川 淳; 杉山 智之; 丸山 結
Proceedings of 14th International Conference on Probabilistic Safety Assessment and Management (PSAM-14) (USB Flash Drive), 10 Pages, 2018/09
Several types of uncertainties exist during the simulation of a severe accident. These may result from incomplete knowledge about the plant systems, accident progression and oversimplified numerical models. Among them, parameter uncertainty can be treated via Monte-Carlo-sampling-based methods. To tackle the severe accident scenario uncertainty, we must resort to advanced dynamic probabilistic risk assessment (PRA) methods. In this paper, authors reviewed the previous dynamic PRA methods and tools, and then performed a preliminary scenario uncertainty analysis, by using an integrated SA code (THALES2) and a scenario generator (RAPID, risk assessment with plant interactive dynamics), both being developed at JAEA. THALES2 is a fast-running severe accident code for the simulation of severe accident progression and source term in light water reactors. Typical scenarios of station-blackout (SBO)-initiated accidents in boiling water reactors are generated and simulated. The dynamic event tree (DET) method is applied to consider the stochastic uncertainties during the scenario progression. Major groups of SBO sequences with the similar accident characteristics can be found. To provide a reference value for risk, a conditional core damage frequency is calculated accordingly. This is a preliminary analysis for severe accident scenario uncertainty quantification at JAEA, and further DPRA researches are in progress.
塩津 弘之; 石川 淳; 杉山 智之; 丸山 結
Journal of Nuclear Science and Technology, 55(4), p.363 - 373, 2018/04
The influences of chemical speciation for Cs-I-Te-Mo-Sn-B-C-O-H system, simulating a state in the reactor cooling system (RCS) of BWR, on pH of the suppression chamber (S/C) water pool were analytically investigated with PHREEQC code. Major conditions were chosen on the basis of the outputs from a BWR severe accident analysis by THALES2 code and chemical thermodynamic analysis with VICTORIA2.0 code. The chemical thermodynamic analysis showed that the chemical speciation of important volatile FPs, Cs and I, was strongly influenced by Mo and BC control material. As a consequence, pH of the S/C water pool was predicted to range from approximately 6 to 10, depending on the fraction of volatile FPs transported from the RCS to the S/C water pool and the H/HO ratio associated with the oxygen potential. It was implied that the formation of volatile I species such as I in the S/C water pool was larger by 3 orders at the lowest pH than that at the highest pH.
Zheng, X.; 玉置 等史; 塩津 弘之; 杉山 智之; 丸山 結
Proceedings of Asian Symposium on Risk Assessment and Management 2017 (ASRAM 2017) (USB Flash Drive), 11 Pages, 2017/11
Nuclear reactor severe accident simulation involves uncertainties, which may result from incompleteness of modeling of accident scenarios, selection of alternative models and unrealistic setting of parameters during the numerical simulation, etc. Both deterministic and probabilistic methods are required to reach reasonable estimation of risk for severe accidents. Computational codes are widely used for the deterministic accident simulations. Bayesian approaches, including both parametric and nonparametric, are applied to the simulation-based severe accident researches at Japan Atomic Energy Agency (JAEA). In the paper, an overview of these research activities is introduced: (1) Dirichlet process models, a nonparametric Bayesian approach, are applied to source term uncertainty and sensitivity analyses; (2) Gaussian process models are applied to the optimization for operations of severe accident countermeasures; (3) Nonparametric models, include models based on Dirichlet process and K-nearest neighbors algorithm, are built to predict the chemical forms of fission products. Simplified models are integrated into the integral severe accident code, THALES2/KICHE; (4) We have also launched the research of dynamic probabilistic risk assessment (DPRA), and because a great number of accident scenarios will be generated during DPRA, Bayesian approaches would be useful for the boosting of computational efficiency.
玉置 等史; 吉田 一雄; 阿部 仁; 杉山 智之; 丸山 結
Proceedings of Asian Symposium on Risk Assessment and Management 2017 (ASRAM 2017) (USB Flash Drive), 9 Pages, 2017/11
堀田 亮年*; 森田 彰伸*; 梶本 光廣*; 丸山 結
日本原子力学会和文論文誌, 16(3), p.139 - 152, 2017/09
Among twelve FCI cases conducted in the OECD/NEA/CSNI/SERENA2 test series using two facilities, six steam explosion cases, five from TROI and one from KROTOS, were analyzed by JASMINE V.3. Major model parameters were categorized into "focused zone", a core part of interest, and "peripheral zone", the initial and boundary conditions given intentionally for each test case. For the former, base values established through past validation studies of JASMINE V.3 were applied. The code was modified to implement the measured distribution of entrained droplet size acquired in TROI-VISU. For the latter, melt release histories were given as a combination of time tables of jet diameter and release velocity that were estimated based on image data and transit timing data of the melt leading edge. The base values were shown to predict impulse responses of SERENA2 systematically with a reasonable error band. A statistical analysis based on the LHS method was performed. Uncertainty ranges were given based on measurement errors and past validation studies in the JASMINE development. Underlying mechanisms causing apparent differences in the mechanical energy conversion ratio between two facilities were studied from the view point of breakup length and trigger timing.
松本 俊慶; 佐藤 允俊; 杉山 智之; 丸山 結
Proceedings of 25th International Conference on Nuclear Engineering (ICONE-25) (CD-ROM), 6 Pages, 2017/07
Hydrogen combustion including deflagration and detonation could become a significant threat to the integrity of containment vessel or reactor building in a severe accident of nuclear power stations. In the present study, numerical analyses were carried out for the ENACCEF No.153 test to develop computational techniques to evaluate the flame acceleration phenomenon during the hydrogen deflagration. This experiment investigated flame propagation in the hydrogen-air premixed gas in a vertical channel with flow obstacles. The reactingFoam solver of the open source CFD code, OpenFOAM, was used for the present analysis. Nineteen elementary chemical reactions were considered for the overall process of the hydrogen combustion. For a turbulent flow, renormalization group (RNG) k-e two-equation model was used in combination with wall functions. Three manners of nodalization were applied and its influences on the flame propagation acceleration were discussed.
石川 淳; 塩津 弘之; 杉山 智之; 丸山 結
Proceedings of 25th International Conference on Nuclear Engineering (ICONE-25) (CD-ROM), 7 Pages, 2017/07
The Japan Atomic Energy Agency (JAEA) is pursuing the development and application of the methodologies on fission product (FP) chemistry for source term analysis by using the integrated severe accident analysis code THALES2. In the present study, models for the eutectic interaction of boron carbide (BC) with steel and the BC oxidation were incorporated into THALES2 code and applied to the source term analyses for a boiling water reactor (BWR) with Mark-I containment vessel (CV). Two severe accident sequences with drywell (D/W) failure by overpressure initiated by loss of core coolant injection (TQUV sequence) and long-term station blackout (TB sequence) were selected as representative sequences. The analyses indicated that a much larger amount of species from the BC oxidation was produced in TB sequence than TQUV sequence. More than a half of carbon dioxide (CO) produced by the BC oxidation was predicted to dissolve into the water pool of the suppression chamber (S/C), which could largely influence pH of the water pool and consequent formation and release of volatile iodine species.
佐藤 允俊; 松本 俊慶; 杉山 智之; 丸山 結
Proceedings of 8th European Review Meeting on Severe Accident Research (ERMSAR 2017) (Internet), 12 Pages, 2017/05
In this study, thermofluid dynamic analyses were carried out for 3 tests, HR3, HR5 and HR12 in the OECD/NEA THAI project, with the Passive autocatalytic recombiner (PAR) manufactured by AREVA. The major parameters in these 3 tests were the initial pressure and steam concentration in the test vessel. The analyses were performed with an open source computational fluid dynamics code, OpenFOAM. The solver was modified by embedding the correlation equations of hydrogen recombination rate for the PAR. The results from the present analyses indicated that the modified solver well reproduced the measured characteristics for PAR behaviour such as hydrogen recombination rate, flow velocity and temperature distribution, hydrogen and oxygen concentration, and so on.
Zheng, X.; 石川 淳; 杉山 智之; 丸山 結
Nuclear Engineering and Technology, 49(2), p.434 - 441, 2017/03
Containment venting is one of essential measures to protect the integrity of the final barrier of a nuclear reactor during severe accidents, by which the uncontrollable release of fission products can be avoided. The authors seek to develop an optimization approach, from a simulation-based perspective, to the venting operations by using an integrated severe accident code, THALES2/KICHE. The effectiveness of containment venting strategies needs to be verified via numerical simulations based on various settings of venting conditions. The number of iterations, however, needs to be controlled for cumbersome computational burden of integrated codes. Bayesian optimization is an efficient global optimization approach. By using Gaussian process regression, a surrogate model of the "black-box" code is constructed. It can be updated simultaneously whenever new simulation results are acquired. With predictions via the surrogate model, upcoming locations of the most probable optimum can be revealed. The sampling procedure is adaptive. The number of code queries is largely reduced for the optimum finding, compared with pure random searches. One typical severe accident scenario of a boiling water reactor is chosen as an example. The research demonstrates the applicability of the Bayesian optimization approach to the design and establishment of containment-venting strategies during severe accidents.
松本 俊慶; 川部 隆平; 杉山 智之; 丸山 結
Proceedings of 10th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-10) (USB Flash Drive), 9 Pages, 2016/11
佐藤 允俊; 松本 俊慶; 杉山 智之; 丸山 結
Proceedings of 10th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-10) (USB Flash Drive), 10 Pages, 2016/11
A numerical analysis was carried out on the thermal-hydraulic behavior during the operation of the PAR for the HR-5 test conducted in the OECD/NEA THAI project. In the HR-5 test, measurements were performed in the test vessel on the volume fractions of oxygen and hydrogen, gas temperature, pressure, flow velocity at the PAR inlet and so on. The open source code OpenFOAM was used for the present study with the reactingFoam solver which is appropriate to treat thermal-hydraulic phenomena including chemical reactions. The code was implemented with the correlation equations for the PAR used in the HR-5 and was modified to be capable of calculating the gas composition change during the recombination of hydrogen and oxygen. Comparison was made between the analysis and experimental results in the gas volume fraction and so on. It was shown that the analyses well reproduced the recombination behavior at the PAR and influences of the recombination heat on the thermal-hydraulic behavior.
Zheng, X.; 石川 淳; 杉山 智之; 丸山 結
Proceedings of 13th Probabilistic Safety Assessment and Management Conference (PSAM-13) (USB Flash Drive), 10 Pages, 2016/10
Containment venting is one of essential measures to protect the integrity of the final barrier of a nuclear reactor, by which the uncontrollable release of fission products can be avoided. The authors seek to develop an optimization approach to the planning of containment-venting operations by using THALES2/KICHE. Factors that control the activation of the venting system, for example, containment pressure, amount of fission products within the containment and pH value in the suppression chamber water pool, will affect radiological consequences. The effectiveness of containment venting strategies needs to be confirmed through numerical simulations. The number of iterations, however, needs to be controlled for cumbersome computational burden of severe accident codes. Bayesian optimization is a computationally efficient global optimization approach to find desired solutions. With the use of Gaussian process regression, a surrogate model of the "black-box" code is constructed. According to the predictions through the surrogate model, the upcoming location of the most probable optimum can be revealed. The number of code queries is largely reduced for the optimum finding, compared with simpler methods such as pure random search. The research demonstrates the applicability of the Bayesian optimization approach to the design and establishment of containment-venting strategies under BWR severe accident conditions.
川口 賢司*; 丸山 結; Zheng, X.
Journal of Artificial Intelligence Research, 56, p.153 - 195, 2016/06
This paper considers global optimization with a black-box unknown objective function that can be non-convex and partly non-smooth. Such a difficult optimization problem arises in many real-world applications, such as parameter tuning in machine learning, engineering design problem, and planning with a complex physics simulator. This paper proposes a new global optimization algorithm, called Locally Oriented Global Optimization (LOGO), to achieve both fast convergence in practice and finite-time error bound in theory. The advantage and usage of the new algorithm are illustrated via theoretical analysis and an experiment conducted with 10 bench-mark test functions. Further, we modify the LOGO algorithm to specifically solve a planning problem with continuous state/action space and long time horizon while maintaining its finite-time error bound. We apply the proposed planning method to severe accident management of a nuclear power plant. The result of the application study demonstrates the practical utility of our method.
端 邦樹; 城戸 健太朗; 西山 裕孝; 丸山 結
NEA/CSNI/R(2016)5 (Internet), p.196 - 203, 2016/05
Model calculations of radiolysis of seawater with iodide were carried out to predict effects of seawater constituents on iodine chemistry at the time of a severe accident. Through the calculations based on -radiolysis of solutions at ambient temperature, it was found that the production of molecular iodine (I) was promoted by the addition of seawater constituents. Especially, Br and HCO had a large impact on its production. The production yields of radiolytic species were affected by the addition of I as well. It was also shown that the pH of seawater is an important parameter to determine radiolytic conditions of iodine. These results imply the potential that the injection of seawater has significant impact on iodine source term.