Yamano, Hidemasa; Takai, Toshihide; Emura, Yuki; Fukuyama, Hiroyuki*; Nishi, Tsuyoshi*; Morita, Koji*; Nakamura, Kinya*; Pellegrini, M.*
Nihon Kikai Gakkai 2023-Nendo Nenji Taikai Koen Rombunshu (Internet), 5 Pages, 2023/09
A research project has been conducting thermophysical property measurement of a eutectic melt, eutectic melting reaction and relocation experiments, eutectic reaction mechanism investigation, and physical model development on the eutectic melting reaction for reactor application analysis in order to simulate the eutectic melting reaction and relocation behavior of boron carbide as a control rod material and stainless steel during a core disruptive accident in an advanced sodium-cooled fast reactor designed in Japan. This paper describes the project overview and progress until JFY2022.
Noi, Hiromi*; Watanabe, Sota*; Kubo, Koji*; Okajima, Satoshi; Ando, Masanori
Nihon Kikai Gakkai M&M 2023 Zairyo Rikigaku Kanfarensu Koen Rombunshu (Internet), p.CL0712_1 - CL0712_5, 2023/09
no abstracts in English
Doda, Norihiro; Uwaba, Tomoyuki; Ohgama, Kazuya; Yoshimura, Kazuo; Nemoto, Toshiyuki*; Tanaka, Masaaki; Yamano, Hidemasa
Nihon Kikai Gakkai Kanto Shibu Dai-29-Ki Sokai, Koenkai Koen Rombunshu (Internet), 5 Pages, 2023/03
An evaluation method for reactivity feedback due to core deformation during reactor power increase in sodium-cooled fast reactors is being developed for realistic core design evaluation. In this evaluation method, fuel assembly bowing was modeled with a beam element of the finite element method, and the assembly's pad contact between adjacent assemblies was modeled with a dedicated element which could consider the wrapper tube cross-sectional distortion and the pad stiffness depending on pad contact conditions. This fuel assembly bowing analysis model was verified for thermal bowing of a single assembly and assembly pad contact between adjacent assemblies in a core as past benchmark problems. The calculation results by this model showed good agreement with those of reference solutions of theoretical solutions or results by participating institutions in the benchmark. This study confirmed that the analysis model was able to calculate thermal assembly bowing appropriately.
Yamano, Hidemasa; Morita, Koji*
Nihon Kikai Gakkai 2022-Nendo Nenji Taikai Koen Rombunshu (Internet), 5 Pages, 2022/09
It is necessary to simulate a eutectic melting reaction and relocation behavior of boron carbide (B4C) as a control rod material and stainless steel (SS) during a core disruptive accident (CDA) in an advanced large-scale sodium-cooled fast reactor (SFR) designed in Japan. A physical model simulating the eutectic reaction and relocation of the eutectic melt was developed to incorporate into the fast reactor severe accident analysis code SIMMER-IV for the CDA numerical analysis of SFRs. This study applied the SIMMER-IV code with the newly developed model to the CDA analysis of the SFR. This analysis indicated that the SIMMER-IV code using the eutectic reaction model has successfully simulated the eutectic reaction and the upward motion of the eutectic melt in the molten core pool as well as the reactivity transient behavior caused by the molten core material relocation.
Ishida, Shinya; Fukano, Yoshitaka
Nihon Kikai Gakkai Rombunshu (Internet), 88(911), p.21-00304_1 - 21-00304_11, 2022/07
In previous studies, the reliability and validity of the SAS4A code was enhanced by applying Phenomena Identification and Ranking Table (PIRT) approach to the Unprotected Loss of Flow (ULOF). SAS4A code has been developed to analyze the early stage of Core Disruptive Accident (CDA), which is named Initiating Phase (IP). In this study, PIRT approach was applied to Unprotected Transient over Power (UTOP), which was one of the most important and typical events in CDA as well as ULOF. The phenomena were identified by the investigation of UTOP event progression and physical phenomena relating to UTOP were ranked. 8 key phenomena were identified and the differences in ranking between UTOP and ULOF were clarified. The code validation matrix was completed and an SAS4A model, which was not validated in ULOF, was identified and validated. SAS4A code became applicable to various scenarios by using PIRT approach to UTOP and the reliability and validity of SAS4A code were significantly enhanced.
Kosaka, Wataru; Uchibori, Akihiro; Yanagisawa, Hideki*; Takata, Takashi; Jang, S.*
Nihon Kikai Gakkai Rombunshu (Internet), 88(905), p.21-00310_1 - 21-00310_9, 2022/01
If a pressurized water/water-vapor leaks from a heat transfer tube in a steam generator (SG) in a sodium-cooled fast reactor (SFR), sodium-water reaction forms high-velocity, high-temperature, and corrosive jet. It would damage the other tubes and might propagate the tube failure in the SG. Thus, it is important to evaluate the effect of the tube failure propagation for safety assessment of SFR. The computational code LEAP-III can evaluate water leak rate during the tube failure propagation with short calculation time, since it consists of empirical formulae and one-dimensional equations of conservation. One of the empirical models, temperature distribution evaluation model, evaluates the temperature distribution in SG as circular arc isolines determined by experiments and preliminary analyses instead of complicated real distribution. In order to improve this model to get more realistic temperature distribution, we have developed the Lagrangian particle method based on engineering approaches. In this study, we have focused on evaluating gas flow in a tube bundle system, and constructed new models for the gas-particles behavior around a tube to evaluate void fraction distribution near the tube. Through the test analysis simulating one target tube system, we confirmed the capability of the models and next topic to improve the models.
Matsushita, Kentaro; Ezure, Toshiki; Imai, Yasutomo*; Fujisaki, Tatsuya*; Tanaka, Masaaki
Nihon Kikai Gakkai Kanto Shibu Ibaraki Koenkai 2021-Nendo Koen Rombunshu (Internet), 4 Pages, 2021/08
For evaluation of gas entrainment phenomenon at free surface in reactor vessel of sodium-cooled fast reactor, the gas entrainment evaluation tool named "Stream Viewer" has been developed. In Stream Viewer, depth of surface vortex dimple is predicted by calculating pressure decrease at the vortex center using velocity distribution around the vortex and Burgers vortex model. In this report, a method to identify continuous vortex center lines from a velocity distribution is newly developed. It becomes possible to evaluate three-dimensional distribution of pressure decrease along vortex center line. Then, the method is validated by applying Stream Viewer to an open channel experiment. As the result, it was confirmed that vortex center lines were successfully identified by the improved Stream Viewer. Moreover, it was also shown that the evaluation accuracy of gas entrainment was expected to be improved by considering distribution of pressure decrease along vortex center line.
Fukasawa, Tsuyoshi*; Miyagawa, Takayuki*; Uchita, Masato*; Yamamoto, Tomohiko; Miyazaki, Masashi; Okamura, Shigeki*; Fujita, Satoshi*
Nihon Kikai Gakkai Rombunshu (Internet), 87(898), p.21-00007_1 - 21-00007_17, 2021/06
This paper describes a fundamental study on the seismic safety margin for the isolated structure using laminated rubber bearings. The variation of the seismic response assumed in the isolated structure will occur under the superposition of "Variations in seismic response due to input ground motions" and "Error with design value accompanying manufacture of the isolation devices ". The seismic response analysis which allows to their conditions is important to assess the seismic safety margin for the isolated structure. This paper clarifies that the seismic safety margin of the isolated structure, which consists of rubber bearings, for Sodium-cooled Fast Reactor (SFR) is ensured against the basis ground motions of Japan Electric Association Guide 4601 (JEAG4601) and SFR through the seismic response analysis considering the variation factors of seismic response. In addition, a relationship between the seismic safety margin and the excess probability of linearity limits is discussed using the results of seismic response analysis.
Furuyama, Taisei*; Thwe Thwe, A.; Katsumi, Toshiyuki; Kobayashi, Hideaki*; Kadowaki, Satoshi
Nihon Kikai Gakkai Rombunshu (Internet), 87(898), p.21-00107_1 - 21-00107_12, 2021/06
The effects of steam addition on the unstable behavior of hydrogen-air lean premixed flames under adiabatic and non-adiabatic conditions were investigated by numerical calculations. Adopting a detailed chemical reaction mechanism of hydrogen-oxyfuel combustion modeled by 17 reversible reactions of 8 active species and diluents, a two-dimensional unsteady reaction flow was treated based on the compressible Navier-Stokes equation. As the steam addition and heat loss increased, the burning velocity of a planar flame decreased and the normalized burning velocity increased. The addition of water vapor promotes the unstable behavior of the hydrogen-air lean premixed flame. This is because the thermal diffusivity of the gas decreases and the diffusion-thermal instability increases. The effect of adding water vapor on the instability of hydrogen premixed flames is a new finding, and it is expected to connect with hydrogen explosion-prevention measures as in NPP.
Suzuki, Tamaki*; Okawa, Teppei*; Harjo, S.; Sasaki, Toshihiko*
Nihon Kikai Gakkai Rombunshu (Internet), 87(894), p.20-00377_1 - 20-00377_15, 2021/02
Kikuchi, Shin; Sakamoto, Kan*; Takai, Toshihide; Yamano, Hidemasa
Nihon Kikai Gakkai 2020-Nendo Nenji Taikai Koen Rombunshu (Internet), 4 Pages, 2020/09
In a postulated severe accidental condition of sodium-cooled fast reactor (SFR), eutectic melting between boron carbide (BC) as control rod element and stainless steel (SS) as control rod cladding or related structure may occur. Thus, behavior of BC-SS eutectic melting is one of the phenomena to evaluate the core disruptive accidents in SFR. In order to clarify the kinetic feature of BC-SS eutectic melting process in the interface, the thinning test for SS crucibles using the pellets of BC or SS with low BC concentration were performed to obtain the rate constant with dependence of BC concentration against SS. It was found that the rate constants of eutectic melting between SS and SS low BC concentration were smaller than that of BC-SS in the high temperature range. Besides, the rate constant of eutectic melting between SS and BC containing SS became small when decreasing the BC concentration against SS.
Yamano, Hidemasa; Takai, Toshihide; Furukawa, Tomohiro
Nihon Kikai Gakkai Rombunshu (Internet), 86(883), p.19-00360_1 - 19-00360_13, 2020/03
It is necessary to simulate a eutectic melting reaction and relocation behavior of boron carbide (BC) as a control rod material and stainless steel (SS) during a core disruptive accident in an advanced sodium-cooled fast reactor designed in Japan because the BC-SS eutectic relocation behavior has a large uncertainty in the reactivity history based on a simple calculation. A physical model simulating the eutectic melting reaction and relocation was developed and implemented into a severe accident simulation code. The developed model must be validated by using test data. To validate the physical model, therefore, the visualization tests of SS-BC eutectic melting reaction was carried out by contacting SS melts of several kg with a BC pellet heated up to about 1500 C. The tests have shown the eutectic reaction visualization as well as freezing and relocation of the BC-SS eutectic in upper part of the solidified test piece due to the density separation. Post-test material analyses by using X-ray diffraction and transmission electron microscope techniques have indicated that FeB appeared at the BC-SS contact interface and (Fe,Cr)B at the top surface of the test piece. Glow discharge optical emission spectrometry has been applied to quantitative analysis of boron concentration distributions. The boron concentration was high at the upper surface and near the original position of the BC pellet.
Uchibori, Akihiro; Yanagisawa, Hideki*; Takata, Takashi; Ohshima, Hiroyuki
Nihon Kikai Gakkai Rombunshu (Internet), 86(883), p.19-00353_1 - 19-00353_6, 2020/03
Evaluation of occurrence possibility of tube failure propagation under sodium-water reaction accident is an important issue. In this study, a numerical analysis method to predict occurrence of failure propagation by overheating rupture was constructed to expand application range of an existing computer code. Applicability of the method was constructed through the numerical analysis of the experiment on water vapor discharging in liquid sodium.
Aoyagi, Mitsuhiro; Uchibori, Akihiro; Kikuchi, Shin; Takata, Takashi; Ohno, Shuji; Ohshima, Hiroyuki
Nihon Kikai Gakkai Rombunshu (Internet), 86(883), p.19-00366_1 - 19-00366_8, 2020/03
Sodium fire is one of key issues in sodium-cooled fast reactor plant. JAEA has developed sodium fire analysis codes, such as AQUA-SF and SPHINCS, to evaluate the consequence of sodium fire events. This paper describes the PIRT (Phenomena Identification and Ranking Table) process for sodium fire events. Ranking table for important phenomena and an assessment matrix are completed. As a part of comprehensive validation based on the assessment matrix, experimental analyses using the AQUA-SF and SPHINCS codes for a sodium spray fire experiment Run-E1 show good agreement with the experimental result.
Machida, Hideo*; Koizumi, Yu*; Wakai, Takashi; Takahashi, Koji*
Nihon Kikai Gakkai M&M 2019 Zairyo Rikigaku Kanfarensu Koen Rombunshu (Internet), p.OS1307_1 - OS1307_5, 2019/11
This paper describes the fracture test and fracture analysis of a pipe under displacement control load. In order to grasp the fracture behavior of the circumferential through-wall cracked pipe, which is important in evaluating the feasibility of leak before break (LBB) in sodium cooled reactor piping, a fracture test in case of a circumferential throughwall crack in the weld line between an elbow and a straight pipe was carried out. From this test, it was found that no pipe fracture occurs in the displacement control loading condition even if a large circumferential through-wall crack (180) was assumed. The fracture analysis of the pipe was carried out using Gurson's parameters set based on the tensile test results of the tested pipe material. The analytic results agree well with the test results, and it was found that it will be possible to predict the fracture behavior of sodium cooled reactor piping.
Kadowaki, Satoshi; Nogami, Masato*; Thwe Thwe, A.; Katsumi, Toshiyuki*; Yamazaki, Wataru*; Kobayashi, Hideaki*
Nihon Kikai Gakkai Rombunshu (Internet), 85(879), p.19-00274_1 - 19-00274_13, 2019/11
We dealt with three-dimensional cellular premixed flames generated by hydrodynamic and diffusive-thermal instabilities to elucidate the effects of unburned-gas temperature and heat loss by adopting the three-dimensional compressible Navier-Stokes equation. As the unburned-gas temperature became lower and the heat loss became larger, the growth rate decreased and the unstable range narrowed. With a decrease of unburned-gas temperature, the normalized growth rate increased and the normalized unstable range widened, which was because the temperature ratio of burned and unburned gases became larger. The obtained hexagonal cellular fronts were qualitatively consistent with the experimental results. As the heat loss became larger, the burning velocity of a cellular flame normalized by that of a planar flame increased. This was because diffusive-thermal effects became stronger owing to the increase of apparent Zeldovich number caused by the decrease of flame temperature.
Muramatsu, Toshiharu; Sato, Yuji; Kamei, Naomitsu; Aoyagi, Yuji*; Shobu, Takahisa
Nihon Kikai Gakkai Dai-13-Kai Seisan Kako, Kosaku Kikai Bumon Koenkai Koen Rombunshu (No.19-307) (Internet), p.157 - 160, 2019/10
no abstracts in English
Shintaku, Yuichi*; Shinozaki, Yuto*; Fujiwara, Takaki*; Takahashi, Akiyuki*; Kikuchi, Masanori
Nihon Kikai Gakkai Rombunshu (Internet), 85(876), p.19-00141_1 - 19-00141_15, 2019/08
The contribution of this paper is to develop two kinds of numerical simulation method for fatigue crack propagation with plastic-induced crack closure under different cyclic loading conditions. One of the developed methods is Direct Numerical Simulation (DNS) using S-version FEM that allow us to simulate by combining with global mesh only representing whole structure and local mesh including crack. After stress intensity factor is calculated by S-version FEM, crack opening level due to plastic-induced crack closure is determined by elastic-plastic analysis using local mesh which is enough subdivided to realize small plastic zone around crack tip. The crack growth rate considering effect of plastic-induced crack closure is predicted by modified Paris' law in which the stress intensity factor range under cyclic loading is converted into the effective value by the crack opening level. Then, the local mesh is updated by new crack shape determined from crack growth rate. By repeating these processes, our developed method can provide us to simulate fatigue crack propagation with plastic-induced crack closure directly. Another method is simplified one that the effective stress intensity factor range is approximately determined by the relationship between the maximum stress intensity factor and crack opening level as a result of preanalysis using two-dimensional DNS. By comparison of experimental results, it can be confirmed that our developed methods predict propagation of surface crack in specimen under bending and tensile loading conditions.
Nihon Kikai Gakkai Rombunshu (Internet), 85(874), p.18-00348_1 - 18-00348_9, 2019/06
Yada, Hiroki; Takaya, Shigeru; Wakai, Takashi; Nakai, Satoru; Machida, Hideo*
Nihon Kikai Gakkai Rombunshu (Internet), 84(859), p.17-00389_1 - 17-00389_15, 2018/03
no abstracts in English