Mass flow rate sensitivity and uncertainty analysis in natural circulation boiling water reactor core from Monte Carlo simulations

Our aim was to evaluate the sensitivity and uncertainty of mass flow rate in the core on the performance of natural circulation boiling water reactor (NCBWR). This analysis was carried out through Monte Carlo simulations of sizes up to 40,000, and the size, i.e., repetition of 25,000 was considered as valid for routine applications. A simplified boiling water reactor (SBWR) was used as an application example of Monte Carlo method. The numerical code to simulate the SBWR performance considers a one-dimensional thermo-hydraulics model along with non-equilibrium thermodynamics and non-homogeneous flow approximation, one-dimensional fuel rod heat transfer. The neutron processes were simulated with a point reactor kinetics model with six groups of delayed neutrons. The sensitivity was evaluated in terms of 99% confidence intervals of the mean to understand the range of mean values that may represent the entire statistical population of performance variables. The regression analysis with mass flow rate as the predictor variable showed statistically valid linear correlations for both neutron flux and fuel temperature and quadratic relationship for the void fraction. No statistically valid correlation was observed for the total heat flux as a function of the mass flow rate although heat flux at individual nodes was positively correlated with this variable. These correlations are useful for the study, analysis and design of any NCBWR. The uncertainties were propagated as follows: for 10% change in the mass flow rate in the core, the responses for neutron power, total heat flux, average fuel temperature and average void fraction changed by 8.74%, 7.77%, 2.74% and 0.58%, respectively.     

Safety analysis of reactivity insertion accidents in a heavy water nuclear research reactor core using coupled 3D neutron kinetics thermal-hydraulic system code technique

Nuclear power plant Safety analysis using coupled 3D neutron kinetics/thermal-hydraulic codes technique is increasingly used nowadays. Actually, the use of this technique allows getting less conservatism and more realistic simulations of the physical phenomena. The challenge today is oriented toward the application of this technique to the operating conditions of nuclear research reactors. In the current study, a three-Dimensional Neutron Kinetics and best estimate Thermal-Hydraulic model based upon the coupled PARCS/RELAP5 codes has been developed and applied for a heavy water research reactor. The objective is to perform safety analysis related to design accidents of this reactor types. In the current study two positive reactivity insertion transients are considered, SCRAM protected and self-limiting power excursion cases. The results of the steady state calculations were compared with results obtained from conventional diffusion codes, while transient calculations were assessed using the point kinetic model of the RELAP5 code. Through this study, the applicability and the suitability of using the coupled code technique with respect to the classical models are emphasized and discussed.     

A coupled analysis of system thermal-hydraulics and three-dimensional reactor kinetics for a 12-finger control element assembly drop event in a PWR plant

A coupled system thermal-hydraulics (T-H) and three-dimensional reactor kinetics code, MARS/MASTER, was developed to attain more accurate predictions for nuclear system transients that involve strong interactions between neutronic and T-H phenomena. In this paper, a 12-finger control element assembly (CEA) drop event in a two-loop pressurized water reactor (PWR) plant under a full power condition was analyzed, where the 12-finger CEA that is nearest to the hot leg of Loop 2 is assumed to incidentally drop. This instantaneously results in an asymmetric radial power distribution and, in turn, asymmetric loop behavior, which may lead to a reactor trip due to a low departure from nucleate boiling (DNB) ratio at the intact side of the core or an excessive difference between the cold leg coolant temperatures. This event clearly requires a coupled calculation of system T-H and three-dimensional reactor kinetics to realistically investigate the thermal-hydraulic behavior of the reactor core. A simple theoretical modeling is also devised to evaluate the cold leg temperature difference under a quasi-steady state.     

Identification of reactor internals'' vibration modes of a Korean standard PWR using structural modeling and neutron noise analysis

The vibration characteristics of a Korean standard PWR reactor internals have been estimated through a three-dimensional finite element analyses and verified by using the mode separated power spectral density functions obtained from the ex-core neutron noise signals. Also the natural vibration modes of the fuel assembly have been identified measuring both the ex-core and the in-core neutron noise signals which are close to each other. As a result, the fundamental bending mode frequency of the reactor internal structure is found to be around 8 Hz and the fundamental shell mode frequency 14.5 Hz, respectively. It is also shown that the fundamental bending mode frequency of the fuel assembly is 2.3 Hz and the 2^nd bending mode frequency 5.8 Hz, respectively. These results can be used for the supplements of the Korean standard PWR's CVAP (Comprehensive Vibration Assessment Program) data.     

Comparison of hydrogen generation for TVSM and TVSA fuel assemblies for water water energy reactor (VVER)-1000

This paper presents the results received during investigation of hydrogen generation for both types fuel assemblies—the old modernistic type of fuel assemblies (TVSM) and recently installed new one alternative type of fuel assemblies (TVSA) in case of severe accident. There are some differences between both types FAs. They have different geometry as well as different burnable poisons. To investigate behavior of new fuel assemblies during the severe conditions it have been performed comparison of fuel behavior of old type TVSM fuel assembly to new one TVSA.To perform this investigation it has been used MELCOR “input model” for Kozloduy Nuclear Power Plant (KNPP) VVER 1000. The model was developed by Institute for Nuclear Research and Nuclear Energy-Bulgarian Academy of Sciences (INRNE-BAS) for investigation of severe accident scenarios and Probabilistic Safety Analyses (PSA) level 2. The model provides a significant analytical capability for the Bulgarian technical specialists, working in the field of the NPP safety, for analysis of core and containment damaged states and the estimation of radionuclides release outside fuel cladding.It was accepted criteria for vessel integrity about hydrogen concentration to be 8%. This criterion was based on the decision of RSK (Germany commission for reactor safety).Generally based on the received results it was made conclusion that using both types of fuel assemblies it was not disturbance safety conditions of NPP.     

Analysis of different containment models for IRIS small break LOCA, using GOTHIC and RELAP5 codes

Advanced nuclear water reactors rely on containment behaviour in realization of some of their passive safety functions. Steam condensation on containment walls, where non-condensable gas effects are significant, is an important feature of the new passive containment concepts, like the AP600/1000 ones.In this work the international reactor innovative and secure (IRIS) was taken as reference, and the relevant condensation phenomena involved within its containment were investigated with different computational tools. In particular, IRIS containment response to a small break LOCA (SBLOCA) was calculated with GOTHIC and RELAP5 codes. A simplified model of IRIS containment drywell was implemented with RELAP5 according to a sliced approach, based on the two-pipe-with-junction concept, while it was addressed with GOTHIC using several modelling options, regarding both heat transfer correlations and volume and thermal structure nodalization. The influence on containment behaviour prediction was investigated in terms of drywell temperature and pressure response, heat transfer coefficient (HTC) and steam volume fraction distribution, and internal recirculating mass flow rate. The objective of the paper is to preliminarily compare the capability of the two codes in modelling of the same postulated accident, thus to check the results obtained with RELAP5, when applied in a situation not covered by its validation matrix (comprising SBLOCA and to some extent LBLOCA transients, but not explicitly the modelling of large dry containment volumes).The option to include or not droplets in fluid mass flow discharged to the containment was the most influencing parameter for GOTHIC simulations. Despite some drawbacks, due, e.g. to a marked overestimation of internal natural recirculation, RELAP5 confirmed its capability to satisfactorily model the basic processes in IRIS containment following SBLOCA.