Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Kamiya, Tomohiro; Yoshida, Hiroyuki
Proceedings of the Symposium on Shock Waves in Japan (Internet), 7 Pages, 2024/03
We developed a ghost fluid method satisfying conservation laws to simulate steam explosions that can occur at the accident of a nuclear power plant. In the developed method, a first-order approximation is applied to interface effect regions, and a high-order approximation is applied to bulk regions. In other words, the algorithm of the developed method is not consistent. Therefore, we modify the way of getting ghost fluids and propose a comprehensive algorithm that applies a high-order approximation to interface effect regions. In the presentation, we will report the outlines and results of the numerical tests of it.
Kamiya, Tomohiro; Yoshida, Hiroyuki
Dai-37-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 8 Pages, 2023/12
We developed a sharp-interface method satisfying a conservation law for a compressible two-phase flow. In this presentation, the outline and numerical test results of the developed method in multi-dimension were reported. The ghost fluid method does not cause numerical diffusion at a gas-liquid interface because difference between gas and liquid phases is avoided. It cannot satisfy the conservation law because cells in which liquid and gas coexist are not prepared although in fact an interface crosses a cell. Hence, we developed the ghost fluid method satisfying a conservation law by preparing cells in which liquid and gas coexist by VOF method. Multi-dimensional basic equations are solved by a split method which is one of the geometric VOF methods. We solved an underwater explosion problem and confirmed that gas bubble expansion and compressible wave propagation which are observed in the steam explosion can be represented and developed method satisfies the conservation law.
unheated rod bundleNagatake, Taku; Yoshida, Hiroyuki
Journal of Nuclear Science and Technology, 60(11), p.1417 - 1430, 2023/11
Times Cited Count：0 Percentile：0.01(Nuclear Science & Technology)In recent years, computational fluid dynamics (CFD) codes have been used to evaluate two-phase flow behavior inside a fuel bundle for nuclear core design and accident management. A space-time distribution of void fraction and interfacial velocity in bundle systems at high temperatures and pressures, are important for validation of two-phase flow CFD codes. However, it is difficult to obtain a space-time distribution of void fraction and interfacial velocity in a bundle system at high temperature and pressure conditions. We have so far developed an experimental apparatus with a unheated rod bundle by adapting a through-rod WMS to measure distributions of a void fraction and an interfacial velocity in high pressure conditions. We newly measured distributions of the void fraction and interfacial velocity in the water-vapor system under high pressure up to 2.6 MPa by the developed apparatus. It has been confirmed that reasonable results were obtained by the experimental apparatus.
Horiguchi, Naoki; Yoshida, Hiroyuki; Kaneko, Akiko*; Abe, Yutaka*
Nihon Kikai Gakkai Kanto Shibu Dai-29-Ki Sokai, Koenkai Koen Rombunshu (Internet), 5 Pages, 2023/10
To elucidate the behavior of molten fuels as a liquid jet in a shallow pool, which is assumed in a core meltdown accident of an LWR, and develop the evaluation method, we investigated the behavior of the vortical liquid film of the simulated wall-impinging liquid jet using 3-dimensional interface shape data obtained by the experiment in a liquid-liquid system.
4 simulated fuel bundle for validation of thermal-hydraulics simulation codesNagatake, Taku; Shibata, Mitsuhiko; Uesawa, Shinichiro; Ono, Ayako; Yoshida, Hiroyuki
Dai-27-Kai Doryoku, Enerugi Gijutsu Shimpojiumu Koen Rombunshu (Internet), 5 Pages, 2023/09
JAEA is developing a neutronics/thermal-hydraulics coupling simulation code for light-water reactors. Thermal-hydraulic simulation codes applied to the platform are expected to evaluate void fraction distributions in fuel assemblies under operational conditions, which is necessary for neutron transport simulation, and need to be validated using void fraction distribution data in a rod bundle under high-temperature and high-pressure conditions. To obtain the data for code validation, we have been measuring the instantaneous void fraction distribution in a 44 simulated fuel assembly by a wire mesh sensor. In this paper, we report the results of the experiments with pressure and flow rate as parameters at a maximum pressure of 2.6 MPa.
Zhang, H.*; Umehara, Yutaro*; Yoshida, Hiroyuki; Mori, Shoji*
International Journal of Heat and Mass Transfer, 211, p.124253_1 - 124253_13, 2023/09
Times Cited Count：2 Percentile：59.37(Thermodynamics)Yamashita, Susumu; Sato, Takumi; Nagae, Yuji; Kurata, Masaki; Yoshida, Hiroyuki
Journal of Nuclear Science and Technology, 60(9), p.1029 - 1045, 2023/09
Times Cited Count：0 Percentile：0.01(Nuclear Science & Technology)Yamashita, Susumu; Uesawa, Shinichiro; Ono, Ayako; Yoshida, Hiroyuki
Mechanical Engineering Journal (Internet), 10(4), p.22-00485_1 - 22-00485_25, 2023/08
A detailed evaluation for air cooling of fuel debris in actual reactors will be essential in fuel debris retrieval under dry conditions. To understand the heat transfer in and around fuel debris, which is assumed as a porous medium in the primary containment vessel (PCV) mechanistically, we newly applied the porous medium model to the multiphase and multicomponent computational fluid dynamics code named JUPITER (JAEA Utility Program for Interdisciplinary Thermal-hydraulics Engineering and Research). We applied the Darcy-Brinkman model as for the porous medium model. This model has high compatibility with JUPITER because it can treat both a pure fluid and a porous medium phase simultaneously in the same manner as the one-fluid model in multiphase flow simulation. We addressed the case of natural convection with a high-velocity flow standing out nonlinear effects by implementing the Forchheimer model, including the term of the square of the velocity as a nonlinear effect to the momentum transport equation of JUPITER. We performed some simple verification and validation simulations, such as the natural convection simulation in a square cavity and the natural convective heat transfer experiment with the porous medium, to confirm the validity of the implemented model. We confirmed that the result of JUPITER agreed well with these simulations and experiments. In addition, as an application of the updated JUPITER, we performed the preliminary simulation of air cooling of fuel debris in the condition of the Fukushima Daiichi Nuclear Power Station unit 2 including the actual core materials. As a result, JUPITER calculated the temperature and velocity field stably in and around the fuel debris inside the PCV. Therefore, JUPITER has the potential to estimate the detailed and accurate thermal-hydraulics behaviors of fuel debris.
Horiguchi, Naoki; Yoshida, Hiroyuki; Kaneko, Akiko*; Abe, Yutaka*
Physics of Fluids, 35(7), p.073309_1 - 073309_17, 2023/07
Times Cited Count：0 Percentile：0.01(Mechanics)The atomization of a liquid jet in an immiscible liquid-liquid system is significant for the safety in the nuclear industry field. The Japan Atomic Energy Agency has developed an evaluation method of a melt fuel behavior as a liquid jet in an immiscible liquid-liquid system for subsequence using mechanistic numerical simulation and has investigated liquid jet behavior in a shallow pool through numerical simulations and experiments. The paper clarifies the atomization mechanism in the wall-impinging liquid jet. Herein, the atomization behavior in the wall-impinging liquid jet in a shallow pool in an immiscible liquid-liquid system was studied in terms of droplet formation and flow field using numerical simulation and the dispersed-phase tracking method. The results show that the droplet formation in the liquid film flow of the wall-impinging liquid jet had the three patterns, and we obtained the droplet properties immediately after droplet formation and developed the theoretical criterion regions using the dimensionless numbers for droplet formation. We characterized the patterns by comparing them with the regions and elucidated the droplet formation mechanisms depending on their sources. Moreover, we elucidated that the relationship between droplet formation as the local behaviors of the jet and atomization as the whole behavior.
Horiguchi, Naoki; Yoshida, Hiroyuki; Kitatsuji, Yoshihiro; Hasegawa, Makoto*; Kishimoto, Tadafumi*
Proceedings of 30th International Conference on Nuclear Engineering (ICONE30) (Internet), 7 Pages, 2023/05
From the viewpoint of energy security in Japan and reduction of the environmental load, continuous operation of light water reactors is essential. Since a pH adjuster with enriched Li-7 ions is required for water quality control on PWR, the development of Li-7 enrichment technology is one of the key issues. The multi-channel counter-current electrophoresis (MCCCE) method has been developed as the technology with a low environmental load. To put this method into practical use, it is necessary to understand Li-7 ion behavior in the channel flow and optimize the experimental condition to separate Li-7 and its isotope. In this paper, to understand Li-7 ion behavior in a single channel of the experimental apparatus, a numerical simulation method based on a computational fluid dynamics (CFD) code with a particle tracking method, TPFIT-LPT, was developed. In the method, the motion of multiple ions under the electric field was simulated as a particle with an added velocity by the electric field. The difference in the isotopes was represented by changing of the magnitude of the added velocity. We also considered that although it is impossible to measure the behavior of each ion, it is important to measure the flow velocity of the bulk fluid for the validation of the numerical simulation. We developed a lab-scale experimental apparatus in which the single channel of the actual apparatus was simplified to measure the flow velocity by Particle Image Velocimetry (PIV). We set a pulsation flow condition on the lab-scale experiment, which is one of difficult conditions for the numerical simulation, and measured the velocity. As the result, we confirmed that the pulsation flow was reproduced. We set the measured data as the inlet boundary condition of the numerical simulation and conducted it. As the numerical result, we confirmed the ions affected by the electric field moved upstream with pulsation. We also confirmed the effect of the electric field on the motion of the isotope.
Uesawa, Shinichiro; Yamashita, Susumu; Shibata, Mitsuhiko; Yoshida, Hiroyuki
Proceedings of 30th International Conference on Nuclear Engineering (ICONE30) (Internet), 9 Pages, 2023/05
For contaminated water management in decommissioning Fukushima Daiichi Nuclear Power Stations, reduction in water injection, intermittent injection water and air cooling are considered. However, since there are uncertainties of fuel debris in the PCV, it is necessary to examine and evaluate optimal cooling methods according to the distribution state of the fuel debris and the progress of the fuel debris retrieval work in advance. We have developed a method for estimating the thermal behavior in the air cooling, including the influence of the position, heat generation and the porosity of fuel debris. Since a large-scale thermal-hydraulics analysis of natural convection is necessary for the method, JUPITER developed independently by JAEA is used. It is however difficult to perform the large-scale thermal-hydraulics analysis with JUPITER by modeling the internal structure of the debris which may consist of a porous medium. Therefore, it is possible to analyze the heat transfer of the porous medium by adding porous models to JUPITER. In this study, we report the validation of JUPITER applied the porous model and discuss which heat transfer models are most effective in porous models such as series, parallel and geometric mean models. To obtain validation data of JUPITER for the natural convective heat transfer analysis around the porous medium, we performed the heat transfer and the flow visualization experiments of the natural convection in the experimental system including the porous medium. In the comparison between the experiment and the numerical analysis with each model, the numerical result with the geometric mean model was the closest of the models to the experimental results. However, the numerical results of the temperature and the velocity were overestimated for those experimental results. In particular, the temperature near the interface between the porous medium and air was more overestimated.
4 simulated fuel bundleNagatake, Taku; Shibata, Mitsuhiko; Yoshida, Hiroyuki; Nemoto, Yoshiyuki; Kaji, Yoshiyuki
Journal of Nuclear Science and Technology, 60(3), p.320 - 333, 2023/03
Times Cited Count：0 Percentile：0.01(Nuclear Science & Technology)In Fukushima Daiichi Nuclear Power Plant accident, failure of cooling system for spent fuel pool occurred and there was a concern that the spent fuels were damaged. Then a safety measures for SFP cooling in severe accident condition is required. As a countermeasure for SFP severe accident, it is considered that a portable spray is used for SFP cooling in such condition. In this research project, the numerical simulation methods have been developed in order to evaluate the applicability of portable spray system for cooling SFPs. And experiments were also performed in order to get a knowledge of spray cooling phenomena and validation data for the numerical simulation methods. As one of the experiments, a cooling experiment by using 44 simulated fuel assembly were performed and temperature distribution during spray cooling process were measured. In this paper, the results of a cooling experiment are reported.
Uesawa, Shinichiro; Shibata, Mitsuhiko; Yoshida, Hiroyuki
Konsoryu, 37(1), p.55 - 64, 2023/03
In decommissioning Fukushima Daiichi Nuclear Power Station, the issue is confinement of radioactive aerosols in the primary containment vessel. Although a High Efficiency Particulate Air (HEPA) filter is used to collect the aerosol particles, pretreatment equipment such as a scrubber may be applied to reduce the load of HEPA filters. In the scrubber, the aerosol particles are removed by moving from gas to liquid through gas-liquid interface. Since the collection efficiency (CE) depends on gas-liquid interfacial area, fine bubbles are necessary to obtain high collection efficiency. JAEA developed a new particle removal technique by using bubble breakup phenomenon in a Venturi tube. To confirm usefulness of the technique, we performed the CE measurements and observed gas-liquid two-phase flow in the Venturi tube. In comparison with a straight pipe type, the Venturi type can have removed particles more 1,000 than it. The CE is almost the same as a HEPA filter. In addition, the Venturi type has the enough CE as the pretreatment equipment for various materials of particles such as Kanto loam, SUS and oil. Besides, we clarified that the CE of the Venturi type depended on the gas and liquid flow rates. The CE increases with the increase of the liquid flow rate but decreases with the increase of the gas flow rate. This is because the CE is affected by the bubble breakup phenomenon in the Venturi tube. In the experiment, we confirmed that cavitation number which is the parameter of the bubble breakup was related to the CE of the Venturi type.
Zhang, H.*; Mori, Shoji*; Hisano, Tsutomu*; Yoshida, Hiroyuki
International Journal of Multiphase Flow, 159, p.104342_1 - 104342_15, 2023/02
Times Cited Count：5 Percentile：65.31(Mechanics)Yamashita, Susumu; Kondo, Nao; Sugawara, Takanori; Monji, Hideaki*; Yoshida, Hiroyuki
Journal of Nuclear Science and Technology, 22 Pages, 2023/00
Times Cited Count：0 Percentile：0.01(Nuclear Science & Technology)To confirm the validity of the thermal-hydraulics design tool based on the Ansys Fluent, we used a detailed computational fluid dynamics code named JAEA Utility Program for Interdisciplinary Thermal-hydraulics Engineering and Research (JUPITER) for the thermal-hydraulics around the beam window (BW) of the Accelerator-Driven System (ADS). The Fluent uses the Reynolds-Averaged Navier-Stokes (RANS) model and can quickly calculate the turbulent flow around the BW as a BW design tool. At first, we compared the results of JUPITER with the experimental results using a mock-up BW system in water to confirm the validity of JUPITER. As a result, we confirmed that numerical results are in good agreement with the experimental results. Thus, we showed that JUPITER could be used as a benchmark code. We also performed a benchmark simulation for the Fluent calculation using validated JUPITER to show the applicability of JUPITER as an alternative of experiments. As a result, the mean values around the BW agreed with each other, e.g., the mean velocity profile for stream and horizontal directions. Therefore, we confirmed that JUPITER showed a good performance in validating the thermal-hydraulics design tool as a fluid dynamics solver. Moreover, Fluent has enough accuracy as a thermal-hydraulics design tool for the ADS.
Kamiya, Tomohiro; Ono, Ayako; Tada, Kenichi; Akie, Hiroshi; Nagaya, Yasunobu; Yoshida, Hiroyuki; Kawanishi, Tomohiro
Proceedings of 29th International Conference on Nuclear Engineering (ICONE 29) (Internet), 8 Pages, 2022/11
JAEA started to develop the advanced reactor analysis code JAMPAN (JAEA advanced multi-physics analysis platform for nuclear systems). The current version of JAMPAN handles the continuous energy Monte Carlo code MVP and the detailed thermal-hydraulics analysis code for multiphase and multicomponent JUPITER. JAMPAN is designed to consider the extensibility and it does not depend on the analysis codes. All calculations in JAMAPAN are not directly connected. JAMPAN has data containers, and all input and output data of each analysis code are set in these data containers. JAMPAN will easily exchange the calculation codes and add the other calculations, e.g., structure calculation and irradiation calculation since the input and the output format of each code has no impact on the other calculation codes. The 4 by 4 pin-cell geometry was used as the demonstration calculation of JAMPAN and the physically reasonable calculation results were obtained.
Yamashita, Susumu; Uesawa, Shinichiro; Ono, Ayako; Yoshida, Hiroyuki
Proceedings of 12th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS12) (Internet), 8 Pages, 2022/10
no abstracts in English
Uesawa, Shinichiro; Yamashita, Susumu; Shibata, Mitsuhiko; Yoshida, Hiroyuki
Proceedings of 12th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS12) (Internet), 8 Pages, 2022/10
Ono, Ayako; Sakashita, Hiroto*; Yamashita, Susumu; Suzuki, Takayuki*; Yoshida, Hiroyuki
Proceedings of 12th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS12) (Internet), 7 Pages, 2022/10
The new prediction method of critical heat flux (CHF) of the fuel assemblies based on the mechanism is proposed in this study. The prediction method of CHF based on the mechanism has been needed for a long time to enhance the safety analysis and reduce the design cost. From several experimental findings of the liquid-vapor behavior near the heating surface from the nucleate boiling to the CHF, the authors consider that the macrolayer dryout model will be appropriate to predict the CHF under the reactor condition. The prediction method of the macrolayer thickness and the passage period of vapor mass in the fuel assemblies are needed to predict CHF from the macrolayer dryout model. In this study, the CHF under the forced convection is evaluated by combining the prediction methods for the macrolayer thickness and passage period of vapor mass, which are proposed by authors. The prediction of the CHF under the forced convection is examined and compared with the experimental data.
Horiguchi, Naoki; Yoshida, Hiroyuki; Kaneko, Akiko*; Abe, Yutaka*
Proceedings of 12th Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS12) (Internet), 6 Pages, 2022/10
For safety evaluation of nuclear reactors in severe accidents, it is important to estimate physical quantities of fragments generated from the molten fuel jet, which falls in a pool and breaks up. The evaluation method has been developed for the behavior as liquid jet with hydrodynamic interaction including fuel coolant interaction (FCI). In case of a shallow pool assumed in ex-vessel, the molten fuel jet is assumed to behave as wall-impinging liquid jet and to form liquid film flow spreading on the floor with/without fragmentation. In our research, focusing on hydrodynamic interaction and the transient 3-dimensional spreading on the floor, we have developed the evaluation method by numerical simulation using the two-phase flow simulation code with interface tracking method (TPFIT) developed by JAEA and, the experimental method using the 3D-LIF method in liquid-liquid system for the validation data. In our previous studies, we investigated the wall-impinging liquid jet behavior with fragmentation and observed that the liquid film flow had some characteristic parts transiently. Since it indicates that the quantities change depending on the parts and affect the safety evaluation, it is important to measure the quantities of the fragments generated from each part. This paper explains the measurement of the physical quantities of the fragments generated from each part of the wall-impinging liquid jet in a shallow pool for the validation of the numerical simulation. We conducted an experiment with the 3D-LIF method and segmented the experimental data based on the fragmentation point over the liquid film flow using the dispersed phase tracking method, developed by JAEA. Then, we measured the diameter and amount of the fragments from the segmented experimental data and investigated their changing trend.
