• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.218-226
• /
• 2008

A component-scale two-phase analysis code has been developed for a realistic simulation of two-phase flow transients in a light water nuclear reactor component. In the code, a two-fluid three-field model is adopted and the governing equations are solved on an unstructured mesh. For the numerical solution scheme, the semi-implicit method used in the RELAP5 code was selected, which has been proved to be very stable and accurate for most of practical applications. However, some modifications were needed for its application to an unstructured non-staggered grid. This paper presents the modified semi-implicit numerical method for unstructured grid and the preliminary results of the calculations.

• Nuclear Engineering and Technology
• /
• v.45 no.5
• /
• pp.625-636
• /
• 2013

Pressing demands of economic competitiveness, the need for large-scale deployment, minimizing the need of human intervention, and experience from the past events and incidents at operating reactors have guided the evolution and innovations in reactor technologies. Indian innovative reactor 'AHWR' is a pressure-tube type natural circulation based boiling water reactor that is designed to meet such requirements, which essentially reflect the needs of next generation reactors. The reactor employs various passive features to prevent and mitigate accidental conditions, like a slightly negative void reactivity coefficient, passive poison injection to scram the reactor in event of failure of the wired shutdown systems, a large elevated pool of water as a heat sink inside the containment, passive decay heat removal based on natural circulation and passive valves, passive ECC injection, etc. It is designed to meet the fundamental safety requirements of safe shutdown, safe decay heat removal and confinement of activity with no impact in public domain, and hence, no need for emergency planning under all conceivable scenarios. This paper examines the role of the various passive safety systems in prevention and mitigation of severe plant conditions that may arise in event of multiple failures. For the purpose of demonstration of the effectiveness of its passive features, postulated scenarios on the lines of three major severe accidents in the history of nuclear power reactors are considered, namely; the Three Mile Island (TMI), Chernobyl and Fukushima accidents. Severe plant conditions along the lines of these scenarios are postulated to the extent conceivable in the reactor under consideration and analyzed using best estimate system thermal-hydraulics code RELAP5/Mod3.2. It is found that the various passive systems incorporated enable the reactor to tolerate the postulated accident conditions without causing severe plant conditions and core degradation.

• Nuclear Engineering and Technology
• /
• v.48 no.4
• /
• pp.847-858
• /
• 2016

An overview is given on the work of the Laboratory of Nuclear Energy Systems at ETH, Zurich (ETHZ) and of the Laboratory of Thermal Hydraulics at Paul Scherrer Institute (PSI), Switzerland on tight-lattice bundles. Two-phase flow in subchannels of a tight triangular lattice was studied experimentally and by computational fluid dynamics simulations. Two adiabatic facilities were used: (1) a vertical channel modeling a pair of neighboring sub-channels; and (2) an arrangement of four subchannels with one subchannel in the center. The first geometry was equipped with two electrical film sensors placed on opposing rod surfaces forming the subchannel gap. They recorded 2D liquid film thickness distributions on a domain of $16{\times}64$ measuring points each, with a time resolution of 10 kHz. In the bubbly and slug flow regime, information on the bubble size, shape, and velocity and the residual liquid film thickness underneath the bubbles were obtained. The second channel was investigated using cold neutron tomography, which allowed the measurement of average liquid film profiles showing the effect of spacer grids with vanes. The results were reproduced by large eddy simulation + volume of fluid. In the outlook, a novel nonadiabatic subchannel experiment is introduced that can be driven to steady-state dryout. A refrigerant is heated by a heavy water circuit, which allows the application of cold neutron tomography.

• Nuclear Engineering and Technology
• /
• v.42 no.4
• /
• pp.339-364
• /
• 2010

Following a joint OECD/NEA-IAEA-sponsored meeting to define the current role and future perspectives of the application of Computational Fluid Dynamics (CFD) to nuclear reactor safety problems, three Writing Groups were created, under the auspices of the NEA working group WGAMA, to produce state-of-the-art reports on different aspects of the subject. The work of the second group, WG2, was to document the existing assessment databases for CFD simulation in the context of Nuclear Reactor Safety (NRS) analysis, to gain a measure of the degree of quality and trust in CFD as a numerical analysis tool, and to take initiatives to extend the existing databases. The group worked over the period of 2003-2007 and produced a final state-of-the-art report. The present paper summarises the material gathered during the study, illustrating the points with a few highlights. A total of 22 safety issues were identified for which the application of CFD was considered to potentially bring real benefits in terms of better understanding and increased safety. A list of the existing databases was drawn up and synthesised, both from the nuclear area and from other parallel, non-nuclear, industrial activities. The gaps in the technology base were also identified and discussed. In order to initiate new ways of bringing experimentalists and numerical analysts together, an international workshop -- CFD4NRS (the first in a series) -- was organised, a new blind benchmark activity was set up based on turbulent mixing in T-junctions, and a Wiki-type web portal was created to offer online access to the material put together by the group giving the reader the opportunity to update and extend the contents to keep the information source topical and dynamic.

• Nuclear Engineering and Technology
• /
• v.51 no.7
• /
• pp.1738-1748
• /
• 2019

Previous research has shown that the inclined condition has an impact on the natural circulation (natural circulation) mode operation of Floating Nuclear Power Plant (FNPP) mounted on the movable marine platform. Due to its compact structure, small volume, strong maneuverability, the Integral Pressurized Water Reactor (IPWR) is adopted as marine reactor in general. The OTSGs of IPWR are symmetrically arranged in the annular region between the reactor vessel and core support barrel in this paper. Therefore, many parallel natural circulation loops are built between the core and the OTSGs primary side when the main pump is stopped. and the inclined condition would lead to discrepancies of the natural circulation drive head among the OTSGs in different locations. In addition, the flow rate and temperature nonuniform distribution of the core caused by inclined condition are coupled with the thermal hydraulics parameters maldistribution caused by OTSG group operating mode on low power operation. By means of the RELAP5 codes were modified by adding module calculating the effect of inclined, heaving and rolling condition, the simulation model of IPWR in inclined condition was built. Using the models developed, the influences on natural circulation operation by inclined angle and OTSG position, the transitions between forced circulation (forced circulation) and natural circulation and the effect on natural circulation operation by different OTSG grouping situations in inclined condition were analyzed. It was observed that a larger inclined angle results the temperature of the core outlet is too high and the OTSG superheat steam is insufficient in natural circulation mode operation. In general, the inclined angle is smaller unless the hull is destroyed seriously or the platform overturn in the ocean. In consequence, the results indicated that the IPWR in the movable marine platform in natural circulation mode operation is safety. Selecting an appropriate average temperature setting value or operating the uplifted OTSG group individually is able to reduce the influence on natural circulation flow of IPWR by inclined condition.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2311-2320
• /
• 2022

The main outcomes of the experiments H2P6 performed in the thermal-hydraulics large-scale PANDA facility at PSI in the frame of the OECD/NEA HYMERES-2 project are presented in this article. The experiments of the H2P6 series consists of two PANDA tests characterized by the activation of three (H2P6_1) or one (H2P6_2) cooler(s) in an initially stratified and pressurized containment atmosphere. The initial stratification is defined by a helium-rich region located in the upper part of the vessel and a steam/air atmosphere in the lower part. The activation of the cooler(s) results i) in the condensation of the steam in the vicinity of the cooler(s), ii) the corresponding activation of large scale natural circulation currents in the vessel atmosphere, with the result of iii) the re-distribution and mixing of the Helium stratification initially located in the upper half of the vessel and iv) the continuous pressure decay. The initial helium layer represents hydrogen generated in a postulated severe accident. The main question to be answered by the experiments is whether or not the interaction of the different, localized cooler units would be important for the application of numerical methods. The paper describes the initial and boundary conditions and the experimental results of the H2P6 series with the suggestion of simple scaling laws for both experiments in terms of i) the temperature difference(s) across the cooler(s), ii) the transient steam and helium content and iii) the pressure decay in the vessel. The outcomes of this scaling indicate that the interaction between separate, closely localized units does not play a prominent role for the present experiments. It is therefore reasonable to model several units as one large component with equivalent heat transfer area and total water flow rate.

• Journal of the Korean Society of Systems Engineering
• /
• v.18 no.2
• /
• pp.75-93
• /
• 2022

Machine learning (ML) data-driven meta-model is proposed as a surrogate model to reduce the excessive computational cost of the physics-based model and facilitate the real-time prediction of a nuclear power plant's transient response. To forecast the transient response three machine learning (ML) meta-models based on recurrent neural networks (RNNs); specifically, Long Short Term Memory (LSTM), Gated Recurrent Unit (GRU), and a sequence combination of Convolutional Neural Network (CNN) and LSTM are developed. The chosen accident scenario is a control element assembly withdrawal at power concurrent with the Loss Of Offsite Power (LOOP). The transient response was obtained using the best estimate thermal hydraulics code, MARS-KS, and cross-validated against the Design and control document (DCD). DAKOTA software is loosely coupled with MARS-KS code via a python interface to perform the Best Estimate Plus Uncertainty Quantification (BEPU) analysis and generate a time series database of the system response to train, test and validate the ML meta-models. Key uncertain parameters identified as required by the CASU methodology were propagated using the non-parametric Monte-Carlo (MC) random propagation and Latin Hypercube Sampling technique until a statistically significant database (181 samples) as required by Wilk's fifth order is achieved with 95% probability and 95% confidence level. The three ML RNN models were built and optimized with the help of the Talos tool and demonstrated excellent performance in forecasting the most probable NPP transient response. This research was guided by the Systems Engineering (SE) approach for the systematic and efficient planning and execution of the research.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3423-3440
• /
• 2023

In this work, a multivariate time-series machine learning meta-model is developed to predict the transient response of a typical nuclear power plant (NPP) undergoing a steam generator tube rupture (SGTR). The model employs Recurrent Neural Networks (RNNs), including the Long Short-Term Memory (LSTM), Gated Recurrent Unit (GRU), and a hybrid CNN-LSTM model. To address the uncertainty inherent in such predictions, a Bayesian Neural Network (BNN) was implemented. The models were trained using a database generated by the Best Estimate Plus Uncertainty (BEPU) methodology; coupling the thermal hydraulics code, RELAP5/SCDAP/MOD3.4 to the statistical tool, DAKOTA, to predict the variation in system response under various operational and phenomenological uncertainties. The RNN models successfully captures the underlying characteristics of the data with reasonable accuracy, and the BNN-LSTM approach offers an additional layer of insight into the level of uncertainty associated with the predictions. The results demonstrate that LSTM outperforms GRU, while the hybrid CNN-LSTM model is computationally the most efficient. This study aims to gain a better understanding of the capabilities and limitations of machine learning models in the context of nuclear safety. By expanding the application of ML models to more severe accident scenarios, where operators are under extreme stress and prone to errors, ML models can provide valuable support and act as expert systems to assist in decision-making while minimizing the chances of human error.

• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.452-459
• /
• 2023

The core power control is an important issue for the study of dynamic characteristics in China initiative accelerator driven subcritical system (CiADS), which has direct impact on the control strategy and safety analysis process. The CiADS is an experimental facility that is only controlled by the proton beam intensity without considering the control rods in the current engineering design stage. In order to get the optimized operation scheme with the stable and reliable features, the variation of beam intensity using the continuous and periodic control approaches has been adopted, and the change of collimator and the adjusting of duty ratio have been proposed in the power control process. Considering the neutronics and the thermal-hydraulics characteristics in CiADS, the physical model for the core power control has been established by means of the point reactor kinetics method and the lumped parameter method. Moreover, the multi-inputs single-output (MISO) logical structure for the power control process has been constructed using proportional integral derivative (PID) controller, and the meta-heuristic algorithm has been employed to obtain the global optimized parameters for the stable running mode without producing large perturbations. Finally, the verification and validation of the control method have been tested based on the reference scenarios in considering the disturbances of spallation neutron source and inlet temperature respectively, where all the numerical results reveal that the optimization method has satisfactory performance in the CiADS core power control scenarios.

• Nuclear Engineering and Technology
• /
• v.18 no.3
• /
• pp.200-208
• /
• 1986

The loss of coolant accident based on a double-ended cold leg break is analyzed with the discharge coefficient (Ca) of 0.4. This analysis covers the whole transient period from the start of depressurization to the complete refilling of the core by using RELAP4/MOD6-EM and RELAP4/ MOD6-HOT CHANNEL for the system thermal-hydraulics and the fuel performance during the blowdown phase respectively, and RELAP4/MOD6-FLOOD and TOODEE2 during the reflood phase. A simple analytical method has been developed to account for the lower plenum filling by approximating steam-water countercurrent flows and superheated wall effects at the downcomer during the refill period. Based on the informations. at the time of EOB (end-of-bypass), the refill duration time and the initial reflooding temperature were estimated and compared with the results from the RELAP4/MOD6, resulting in a good agreement. In addition, some parametric studies on the EOB were performed. The form loss coefficient between upper head and upper downcomer was found to be sensitive to the occurrence of the spurious EOB. Appropriate form loss coefficients should be taken into account to avoid the flow oscillations at the downcomer. The analyses with the six and three volume core nodalizations, respectively, show much similar trends in the system thermal-hydraulic performance, but the former case is recommended to obtain good results.