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.
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.
Sugihara, Kenta; Onodera, Naoyuki; Idomura, Yasuhiro; Yamashita, Susumu
Dai-36-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 5 Pages, 2022/12
The conventional Allen-Cahn type multi-phase field method was modified to conserve not only the sum of the masses of all phases but also the mass of each phase. The interface advection calculations within a two-dimensional rotational velocity field were performed as a verification problem, and the conservation was successfully achieved. The proposed method was used to calculate the horizontally aligned pair of bubbles rising, and it was found that the bouncing phenomenon between bubbles can be calculated at 1/50 resolution of the high-resolution calculation by Zhang et al. using the volume of fluid method.
Uesawa, Shinichiro; Yoshida, Hiroyuki
Konsoryu Shimpojiumu 2022 Koen Rombunshu (Internet), 2 Pages, 2022/08
no abstracts in English
Sugihara, Kenta; Onodera, Naoyuki; Idomura, Yasuhiro; Yamashita, Susumu
Keisan Kogaku Koenkai Rombunshu (CD-ROM), 27, 5 Pages, 2022/06
The phase-field method has been successfully applied to various multi-phase flow problems as an interface tracking method for gas-liquid interfaces. However, the accuracy of the phase-field method depends on hyper-parameters, which are empirically adjusted for each problem. The phase-field method sustains sharp interfaces by the balance between the numerical viscosity of the advection term and the interface modification by the diffusion and anti-diffusion terms. Based on this fact, we propose a method for deriving the optimal hyper-parameters in a non-empirical manner by performing a basic error analysis of the interface advection.
Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Yamashita, Susumu; Shimokawabe, Takashi*; Aoki, Takayuki*
Proceedings of International Conference on High Performance Computing in Asia-Pacific Region (HPC Asia 2021) (Internet), p.120 - 128, 2021/01
Times Cited Count:0 Percentile:0.01(Computer Science, Hardware & Architecture)We develop a multigrid preconditioned conjugate gradient (MG-CG) solver for the pressure Poisson equation in a two-phase flow CFD code JUPITER. The MG preconditioner is constructed based on the geometric MG method with a three-stage V-cycle, and a RB-SOR smoother and its variant with cache-reuse optimization (CR-SOR) are applied at each stage. The numerical experiments are conducted for two-phase flows in a fuel bundle of a nuclear reactor. The MG-CG solvers with the RB-SOR and CR-SOR smoothers reduce the number of iterations to less than 15% and 9% of the original preconditioned CG method, leading to 3.1- and 5.9-times speedups, respectively. The obtained performance indicates that the MG-CG solver designed for the block-structured grid is highly efficient and enables large-scale simulations of two-phase flows on GPU based supercomputers.
Akiba, Miyuki*; Hotta, Akitoshi*; Abe, Yutaka*; Sun, Haomin
Nihon Genshiryoku Gakkai Wabun Rombunshi, 19(1), p.1 - 15, 2020/02
Tests at three different scales were conducted in order to understand the mechanisms of pool scrubbing. In the small-scale separate effect test, high resolution two-phase flow measurement techniques such as a high-speed camera, wire mesh sensor and PIV were applied to capture the behaviors of a single bubble and two-phase flow structures. In the large-scale integral effect test, the dependence of the aerosol removal efficiency on submergence and pool temperature was measured under constant pressure and depressurized conditions. To clarify relationships between individual phenomena and combined phenomena observed in two tests, the mid-scale integral effect test was undertaken.
Nagayoshi, Takuji*; Yoshida, Hiroyuki; Onuki, Akira; Akimoto, Hajime
Nihon Genshiryoku Gakkai Wabun Rombunshi, 4(1), p.16 - 24, 2005/03
A detailed gas-liquid two-phase flow analysis code based on an advanced interface-tracking method has been developed. It is expected that the developed code would be able to simulate two-phase cross flow behavior within tight-lattice fuel bundles without relying on any empirical correlations. In order to verify the applicability of the code to simulate two-phase cross flow behavior in such situations, numerical analyses of 2-channel model tests were conducted to compare the air slug deformation and separation behavior caused by cross flow through a narrow interconnection between channels. Although the code underestimated the ascending velocity of the slug, the calculated slug deformation and separation behavior were shown to be quite similar to those observed by a high-speed video camera. Moreover the minimum differential pressure between the subchannels through the interconnection, causing channel-to-channel air transfer to occur could be predicted to within 20Pa. However, further studies of modeling and implementation related to the interface-channel wall interaction, such as a contact angle of a gas-liquid interface at the channel wall, are required for prediction improvements. Nevertheless, the qualitative capability of the developed code to simulate two-phase cross flow phenomena was demonstrated.
Takase, Kazuyuki; Yoshida, Hiroyuki; Ose, Yasuo*
Hai Pafomansu Komputingu To Keisan Kagaku Shimpojium (HPCS 2005) Rombunshu, P. 16, 2005/01
no abstracts in English
Takase, Kazuyuki; Yoshida, Hiroyuki; Ose, Yasuo*; Tamai, Hidesada; Akimoto, Hajime
Dai-22-Kai Nihon Shimyureshon Gakkai Taikai Happyo Rombunshu, p.85 - 88, 2003/06
no abstracts in English
Mishima, Kaichiro*; Nakamura, Hideo
JAERI-Review 2001-012, 51 Pages, 2001/03
no abstracts in English
Suzuki, Toru; Tobita, Yoshiharu
JNC TN9400 2000-019, 35 Pages, 2000/03
The transition phase analysis code SIMMER-III has been developed to appropriately evaluate the core disruptive accident in a fast breeder reactor. The momentum exchange model used in the fluid dynamics portion of the code uses the conventional correlation based on ordinary flows such as air-water flows. lt has already been confirmed that this code can represent the experimental results of ordinary flows. However, more detailed research is needed to confirm that this code is applicable to two-phase flow with large liquid density, which would be formed in an actual molten core pool. In addition, since the shapes of bubbles affect their drag in the bubbly flow where the liquid and the gas form continuous and dispersed phases, respectively, it is necessary to take this effect of bubble shape into account to improve SIMMER-III's analytical precision. ln this study, using experimental results obtained through a joint research program with Kyoto university, the momentum exchange model of SIMMER-III is assessed with regard to the bubbly flow regime of two-phase flow with large liquid density, on which experimental data and information on bubble shapes had been lacking. This study suggests that the original SIMMER-III can appropriately represent the characteristics of bubbly flows containing ellipsoidal bubbles with relatively small gas flux. Moreover, this study shows that the precision of SIMMER-III for bubbly flows containing cap bubbles with relatively large gas flux is much improved by using Kataoka-lshii's correlation to determine the drag coefficient of bubbles in the momentum exchange model.
Takase, Kazuyuki; Ose, Yasuo*; Suzuki, Takayuki*; Akimoto, Hajime
Nihon Kikai Gakkai Yamanashi Koenkai Koen Rombunshu (000-4), p.173 - 174, 2000/00
no abstracts in English
Mishima, Kaichiro*; Hibiki, Takashi*; *; Tobita, Yoshiharu
PNC TY9604 96-003, 10 Pages, 1996/05
no abstracts in English
*; *; *
PNC TJ9614 94-001, 59 Pages, 1994/03
Crossflow of a two-phase mixture between vertical subchannels is subdivided into three components in the literature; turbulent mixing, void drift and diversion crossflow. Of these, turbulent mixing alone occurs in an equiliblium flow, in which flow rates of both phases in each subchannel do not change in the axial direction. In a general non-equilibrium flow, however, all three components occur simultaneously. In this report, effect of pressure differential between subchannels on flow redistribution process along the channel axis has been studied experimentally. In the experiment, a multiple channel, consisting of two identical circular subchannels of 16 mm I.D., were used as a test channel. And, air and water were introduced unevenly into the two subchannels at the inlet to get several non-equilibrium flows with and without the pressure differential between subchannels. For each flow, we have obtained the axial distributions data of pressure differential between the subchannels, the air and water flow rates, the void fractions, and the tracer concentrations for both phases when gas and liquid tracers were injected into one of the two subchannels. From these experimental data, we have estimated lateral velocities of the air and water corresponding to each crossflow component, and analyzed the effect of the pressure differential on the lateral velocities.
The Japan Society of Multiphase Flow*; Special Committee for Examination of Thermohydraulic Analysis Code based on Three-Fluid Model*
PNC TJ9565 92-001, 173 Pages, 1992/05
The purpose of the present study is to improve a numerical analysis method for multiphase flows based on the three-fluid model. The improvement will be conducted for the following three items; (1)modeling of the multiphase flow, (2)constitutive equations (3)numerical method. A systematic survey was carried out for the aforementioned three items. As a result, a reliable plan for the improvement was proposed. Experimental database for air-steam and air-water two-phase annular flow was also presented in this report, which will be useful for the future improvement of the constitutive equations. As for the numerical method, a method for improving the computational efficiency was proposed. It was confirmed that the CPU time decreases about 95% with the proposed method.
Sobajima, Makoto
Konsoryu Rekucha Shirizu Dai-3-Kai; Kiso Kara Saizensen Made, p.55 - 67, 1989/00
no abstracts in English