Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Kurisaka, Kenichi
Mechanical Engineering Journal (Internet), 11(2), p.23-00377_1 - 23-00377_14, 2024/04
This study aims to understand the time-dependent change in the occurrence rate of leak from steam generator (SG) tubes in sodium-cooled fast reactors (SFRs). The target SFRs in the present paper are Phenix in France and BN-600 in Russia. By reviewing publicly available literature that show data from the SFRs, we have investigated the numbers of tube-to-tubeplate welds and tube-to-tube welds, heat transfer areas of tube base metal, operating hours of SGs, dates when SG tube leak occurred, locations of leak, and corrective actions taken after tube leak events, such as replacement of the module, in which a leak occurred. Based on these, we have estimated the time to leak and quantitatively analyzed the time-dependent change of the occurrence rates of SG tube leak for each of the above-mentioned parts by hazard plotting method. The results show that the rates of both Phenix and BN-600 decreased over time. For Phenix, this is probably thanks to improved welding and SG operating conditions. For BN-600, it seems that in many cases, the probable cause of the leak was initial defects that developed to failure during the early stage of reactor operation, and that no special countermeasure was taken in the later stages. Therefore, it would be natural to assume that the rate simply decreased over time. The rate of leak at tube-to-tube welds in Phenix shows significant increase in a short term after a certain period of time. This can be caused by thermal stress repeatedly exerted on the materials.
Kubo, Kotaro; Zheng, X.; Tanaka, Yoichi; Tamaki, Hitoshi; Sugiyama, Tomoyuki; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Proceedings of the Institution of Mechanical Engineers, Part O; Journal of Risk and Reliability, 237(5), p.947 - 957, 2023/10
Times Cited Count:5 Percentile:56.65(Engineering, Multidisciplinary)Probabilistic risk assessment (PRA) is a method used to assess the risks associated with large and complex systems. However, the timing at which nuclear power plant structures, systems, and components are damaged is difficult to estimate if the risk of an external event is evaluated using conventional PRA based on event trees and fault trees. A methodology coupling thermal-hydraulic analysis with external event simulations using Risk Assessment with Plant Interactive Dynamics (RAPID) is therefore proposed to overcome this limitation. A flood propagation model based on Bernoulli's theorem was applied to represent internal flooding in the turbine building of the pressurized water reactor. Uncertainties were also taken into account, including the flow rate of the floodwater source and the failure criteria for the mitigation systems. The simulated recovery actions included the operator isolating the floodwater source and using a drainage pump; these actions were modeled using several simplifications. Overall, the results indicate that combining isolation and drainage can reduce the conditional core damage probability upon the occurrence of flooding by approximately 90%.
Kubo, Kotaro; Tanaka, Yoichi; Hakuta, Yuto*; Arake, Daisuke*; Uchiyama, Tomoaki*; Muramatsu, Ken
Mechanical Engineering Journal (Internet), 10(4), p.23-00051_1 - 23-00051_17, 2023/08
The significance of probabilistic risk assessments (PRAs) of nuclear power plants against external events was re-recognized after the Fukushima Daiichi Nuclear Power Plant accident. Regarding the seismic PRA, handling correlated failures of systems, components, and structures (SSCs) is very important because this type of failure negatively affects the redundancy of accident mitigation systems. The Japan Atomic Energy Research Institute initially developed a fault tree quantification methodology named the direct quantification of fault tree using Monte Carlo simulation (DQFM) to handle SSCs' correlated failures in detail and realistically. This methodology allows quantifying the top event occurrence probability by considering correlated uncertainties related to seismic responses and capacities with Monte Carlo sampling. The usefulness of DQFM has already been demonstrated. However, improving its computational efficiency would allow risk analysts to perform several analyses. Therefore, we applied quasi-Monte Carlo and importance sampling to the DQFM calculation of simplified seismic PRA and examined their effects. Specifically, the conditional core damage probability of a hypothetical pressurized water reactor was analyzed with some assumptions. Applying the quasi-Monte Carlo sampling accelerates the convergence of results at intermediate and high ground motion levels by an order of magnitude over Monte Carlo sampling. The application of importance sampling allows us to obtain a statistically significant result at a low ground motion level, which cannot be obtained through Monte Carlo and quasi-Monte Carlo sampling. These results indicate that these applications provide a notable acceleration of computation and raise the potential for the practical use of DQFM in risk-informed decision-making.
Kubo, Kotaro
Science and Technology of Nuclear Installations, 2023, p.7402217_1 - 7402217_12, 2023/06
Times Cited Count:1 Percentile:41.04(Nuclear Science & Technology)Kurisaka, Kenichi
Proceedings of 30th International Conference on Nuclear Engineering (ICONE30) (Internet), 10 Pages, 2023/05
This study aims to understand a time trend of the occurrence rate of steam generator (SG) tube leak in the existing sodium-cooled fast reactors (SFRs) based on the observed data. The target on SFRs in the present paper is Phenix in France and BN600 in Russia. From the open literature review, we investigated the number of tube-to-tube plate weld, the number of tube-to-tube weld, heat transfer area of tube base metal, operating time of SGs, dates when SG tube leak occurred, leaked location, corrective action after tube leak such as replacement of leaked module. Based on these observed data, time to leak is estimated and then time trend of the occurrence rate of SG tube leak for each of the above-mentioned parts was quantitatively analyzed by the hazard plotting method. As a result, the rate of leak at tube-to-tube weld in Phenix shows increase with time due to probable cause of cyclic thermal stress in a short term. As for a long-term trend, the rate of tube leak in both Phenix and BN600 SGs indicated decrease with time probably thanks to improvement in welding and in SG operating condition and to removal of initial failure.
Kurisaka, Kenichi; Nishino, Hiroyuki; Yamano, Hidemasa
Proceedings of 30th International Conference on Nuclear Engineering (ICONE30) (Internet), 8 Pages, 2023/05
The objective of this study is to develop an effectiveness evaluation methodology of the measures for improving resilience of nuclear structures against excessive earthquake by applying the failure mitigation technology. This study regarded those measures for improving resilience of important structures, systems, and components for safety to enlarge their seismic safety margin. To evaluate effectiveness of those measures, seismic core damage frequency (CDF) is selected as an index. Reduction of CDF as an effectiveness index is quantified by applying seismic PRA technology. Accident sequences leading to loss of decay heat removal are significant contributor to seismic CDF of sodium-cooled fast reactors (SFRs), and those sequences result in core damage via ultra-high temperature condition. This study improved the methodology to evaluate not only the measures against shaking due to excessive earthquake but also the measures at the ultra-high temperature condition. To examine applicability of the improved methodology, a trial calculation was implemented with some assumptions for a loop-type SFR. Within the assumption, the measures for improving resilience were significantly effective for decreasing CDF in excessive earthquake up to several times of a design basis ground motion. Through the applicability examination, the methodology for the effectiveness evaluation was developed successfully.
Kubo, Kotaro; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Journal of Nuclear Science and Technology, 60(4), p.359 - 373, 2023/04
Times Cited Count:8 Percentile:83.23(Nuclear Science & Technology)Probabilistic risk assessment (PRA) is an essential approach to improving the safety of nuclear power plants. However, this method includes certain difficulties, such as modeling of combinations of multiple hazards. Seismic-induced flooding scenario includes several core damage sequences, i.e., core damage caused by earthquake, flooding, and combination of earthquake and flooding. The flooding fragility is time-dependent as the flooding water propagates from the water source such as a tank to compartments. Therefore, dynamic PRA should be used to perform a realistic risk analysis and quantification. This study analyzed the risk of seismic-induced flooding events by coupling seismic, flooding, and thermal-hydraulics simulations, considering the dependency between multiple hazards explicitly. For requirements of safety improvement, especially in light of the Fukushima Daiichi Nuclear Power Plant accident, sensitivity analysis was performed on the seismic capacity of systems, and the effectiveness of alternative steam generator injection by a portable pump was estimated. We demonstrate the use of this simulation-based dynamic PRA methodology to evaluate the risk induced by a combination of hazards.
Kubo, Kotaro; Tanaka, Yoichi*; Ishikawa, Jun
Proceedings of the Institution of Mechanical Engineers, Part O; Journal of Risk and Reliability, 11 Pages, 2023/00
Times Cited Count:1 Percentile:33.61(Engineering, Multidisciplinary)Yamaguchi, Akinori*; Yokotsuka, Muneyuki*; Furuta, Masayo*; Kubota, Kazuo*; Fujine, Sachio*; Mori, Kenji*; Yoshida, Naoki; Amano, Yuki; Abe, Hitoshi
Nihon Genshiryoku Gakkai Wabun Rombunshi (Internet), 21(4), p.173 - 182, 2022/09
Risk information obtained from probabilistic risk assessment (PRA) can be used to evaluate the effectiveness of measures against severe accidents in nuclear facilities. The PRA methods used for reprocessing facilities are considered immature compared to those for nuclear power plants, and to make the methods mature, reducing the uncertainty of accident scenarios becomes crucial. In this paper, we summarized the results of literature survey on the event progression of evaporation to dryness caused by boiling of high-level liquid waste (HLLW) which is a severe accident in reprocessing facilities and migration behavior of associated radioactive materials. Since one of the important characteristics of Ru is its tendency to form volatile compounds over the course of the event progression, the migration behavior of Ru is categorized into four stages based on temperature. Although no Ru has been released in the waste in the high temperature region, other volatile elements such as Cs could be released. Sufficient experimental data, however, have not been obtained yet. It is, therefore, necessary to further clarify the migration behavior of radioactive materials that predominantly depends on temperature in this region.
Zheng, X.; Tamaki, Hitoshi; Takahara, Shogo; Sugiyama, Tomoyuki; Maruyama, Yu
Proceedings of Probabilistic Safety Assessment and Management (PSAM16) (Internet), 10 Pages, 2022/09
Kubo, Kotaro; Fujiwara, Keita*; Tanaka, Yoichi; Hakuta, Yuto*; Arake, Daisuke*; Uchiyama, Tomoaki*; Muramatsu, Ken*
Proceedings of 29th International Conference on Nuclear Engineering (ICONE 29) (Internet), 8 Pages, 2022/08
After the Fukushima Daiichi Nuclear Power Plant accident, the importance of conducting probabilistic risk assessments (PRAs) of external events, especially seismic activities and tsunamis, was recognized. The Japan Atomic Energy Agency has been developing a computational methodology for seismic PRA, called the direct quantification of fault tree using Monte Carlo simulation (DQFM). When appropriate correlation matrices are available for seismic responses and capacities of components, the DQFM makes it possible to consider the effect of correlated failures of components connected through AND and/or OR gates in fault trees, which is practically difficult when methods using analytical solutions or multidimensional numerical integrations are used to obtain minimal cut set probabilities. The usefulness of DQFM has already been demonstrated. Nevertheless, a reduction of the computational time of DQFM would allow the large number of analyses required in PRAs conducted by regulators and/or operators. We; therefore, performed scoping calculations using three different approaches, namely quasi-Monte Carlo sampling, importance sampling, and parallel computing, to improve calculation efficiency. Quasi-Monte Carlo sampling, importance sampling, and parallel computing were applied when calculating the conditional core damage probability of a simplified PRA model of a pressurized water reactor, using the DQFM method. The results indicated that the quasi-Monte Carlo sampling works well at assumed medium and high ground motion levels, importance sampling is suitable for assumed low ground motion level, and that parallel computing enables practical uncertainty and importance analysis. The combined implementation of these improvements in a PRA code is expected to provide a significant acceleration of computation and offers the prospect of practical use of DQFM in risk-informed decision-making.
Zheng, X.; Tamaki, Hitoshi; Sugiyama, Tomoyuki; Maruyama, Yu
Reliability Engineering & System Safety, 223, p.108503_1 - 108503_12, 2022/07
Times Cited Count:19 Percentile:86.93(Engineering, Industrial)Kubo, Kotaro; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Journal of Nuclear Science and Technology, 59(3), p.357 - 367, 2022/03
Times Cited Count:6 Percentile:56.19(Nuclear Science & Technology)Dynamic probabilistic risk assessment (PRA), which handles epistemic and aleatory uncertainties by coupling the thermal-hydraulics simulation and probabilistic sampling, enables a more realistic and detailed analysis than conventional PRA. However, enormous calculation costs are incurred by these improvements. One solution is to select an appropriate sampling method. In this paper, we applied the Monte Carlo, Latin hypercube, grid-point, and quasi-Monte Carlo sampling methods to the dynamic PRA of a station blackout sequence in a boiling water reactor and compared each method. The result indicated that quasi-Monte Carlo sampling method handles the uncertainties most effectively in the assumed scenario.
Kubo, Kotaro; Tanaka, Yoichi
Proceedings of Asian Symposium on Risk Assessment and Management 2021 (ASRAM 2021) (Internet), 13 Pages, 2021/10
Probabilistic risk assessment (PRA) is extensively used, e.g., in periodical safety review and the reactor oversight process, in nuclear regulation systems to improve the safety of nuclear power plants; however, one limitation of classical PRA is the handling of temporal information such as system failure and core damage timings. To resolve this limitation, the dynamic PRA method has been developed and applied for multiple safety issues; however, its improvement is accompanied by considerable computational costs. In this study, we applied the polynomial chaos expansion (PCE) technique to dynamic PRA with the expectation of reduction in computational cost. In particular, to estimate core damage timing, a PCE-based surrogate model was developed. Then, the surrogate model was applied to dynamic PRA to calculate the conditional core damage probability and core damage timing. Consequently, applying the PCE might efficiently perform these analyses without considerable reduction in accuracy.
Kubo, Kotaro; Tanaka, Yoichi
Proceedings of 31st European Safety and Reliability Conference (ESREL 2021) (Internet), p.810 - 817, 2021/09
Probabilistic risk assessment (PRA) is a method of effectively evaluating risks in nuclear power plants and is used in various agencies. Dynamic PRA is attracting considerable attention, as it enables realistic assessment by reducing the assumptions and engineering judgments related to time-dependent failure probability and/or human action reliability. However, it is difficult to remove all assumptions and engineering judgments. Therefore, their effects on assessment results should be understood. This study focuses on the "risk dilution effect," which arises from assumptions about uncertainty. Results showed that this effect causes a difference of about 10% to 20% in the relative change of the conditional core damage probability in the station blackout scenario. This effect should be fully considered when using dynamic PRA in critical decision-making, such as that on regulations.
Maruyama, Yu; Yoshida, Kazuo
Nihon Genshiryoku Gakkai-Shi ATOMO, 63(7), p.517 - 522, 2021/07
no abstracts in English
Shimada, Kazumasa; Sakurahara, Tatsuya*; Reihani, S.*; Mohagehgh, Z.*
Proceedings of Asian Symposium on Risk Assessment and Management 2020 (ASRAM 2020) (Internet), 12 Pages, 2020/11
Level 3 Probabilistic Risk Assessment (Level 3 PRA) and Traffic simulation were integrated to evaluate the effects of evacuation more realistically on radiation exposure to residents in the offsite consequence analysis. In this study, WinMACCS was used as the Level 3 PRA code. As a test case, the Sequoyah Nuclear Power Plant(NPP) site, which was targeted by the State-of-the-Art Reactor Consequence Analyzes (SOARCA) issued by U.S. Nuclear Regulatory Commission in 2017, was adopted. The MultiAgent Transport Simulation (MATSim) was used to simulate the evacuation of a Sequoyah NPP's 10-mile Emergency Planning Zone. For the transportation route choice, the route where each vehicle chooses the shortest distance and the route where the total evacuation time is shortened by iterative calculation were chosen. In the calculation of MACCS, the source term with the shortest release start time in the SOARCA report was adopted. As an example of the results, the radiation dose of the residents when the evacuation time was optimized was reduced by about 30% from the dose when the shortest distance was selected. Furthermore, a sensitivity analysis was conducted, and it was shown that the evacuation preparation time was the largest factor that contributed to the radiation dose to residents.
Kubo, Kotaro; Zheng, X.; Tanaka, Yoichi; Tamaki, Hitoshi; Sugiyama, Tomoyuki; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.308 - 315, 2020/10
Dynamic probabilistic risk assessment (PRA) is a method for improving the realism and completeness of conventional PRA. However, enormous calculation costs are incurred by these improvements. One solution is to select an appropriate sampling method. In this paper, we applied the Monte Carlo, Latin hypercube, grid-point, and quasi-Monte Carlo sampling methods to the dynamic PRA of a simplified accident sequence and compared the results for each method. Quasi-Monte Carlo sampling was found to be the most effective method in this case.
Yamano, Hidemasa; Kurisaka, Kenichi; Nishino, Hiroyuki; Okano, Yasushi; Naruto, Kenichi*
Proceedings of 12th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics, Operation and Safety (NUTHOS-12) (USB Flash Drive), 15 Pages, 2018/10
Spent fuels are transferred from a reactor core to a spent fuel pool through an external vessel storage tank (EVST) filled with sodium in sodium-cooled fast reactors in Japan. This paper describes identification of dominant accident sequences leading to fuel failure, which was achieved through probabilistic risk assessment for the EVST designed for a next sodium-cooled fast reactor plant system in Japan to improve the EVST design. The safety strategy for the EVST involves whole core refueling (early transfer of all core fuel assemblies into the EVST) assuming a severe situation that results in sodium level reduction leading finally to the top of the reactor core fuel assemblies in a long time. This study introduces the success criteria mitigation along the decay heat decrease over time. Based on the design information, this study has carried out identification of initiating events, event and fault tree analyses, a probability analysis for human error, and quantification of accident sequences. The fuel damage frequency of the EVST was evaluated to be approx. 10/year. The dominant accident sequence resulted from the static failure and human error for the switching from the stand-by to operation mode in the three stand-by cooling circuits after loss of one circuit for refueling heat removal operation as an initiating phase.
Nishino, Hiroyuki; Yamano, Hidemasa; Kurisaka, Kenichi
Mechanical Engineering Journal (Internet), 5(4), p.18-00079_1 - 18-00079_17, 2018/08