Parametric Analysis of a Turbine Trip Event in a BWR Using a 3D Nodal Code

Two essential thermal hydraulics safety criteria concerning the reactor core are that even during operational transients there is no fuel melting and impermissible cladding temperatures are avoided. A common concept for boiling water reactors is to establish a minimum critical power ratio (MCPR) for steady state operation. For this MCPR it is shown that only a very small number of fuel rods suffers a short-term dryout during the transient. It is known from experience that the limiting transient for the determination of the MCPR is the turbine trip with blocked bypass system. This fast transient was simulated for a German BWR by use of the three-dimensional reactor analysis transient code SIMULATE-3K. The transient behaviour of the hot channels was used as input for the dryout calculation with the transient thermal hydraulics code FRANCESCA. By this way the maximum reduction of the CPR during the transient could be calculated. The fast increase in reactor power due to the pressure increase and to an increased core inlet flow is limited mainly by the Doppler effect, but automatically triggered operational measures also can contribute to the mitigation of the turbine trip. One very important method is the short-term fast reduction of the recirculation pump speed which is initiated e. g. by a pressure increase in front of the turbine. The large impacts of the starting time and of the rate of the pump speed reduction on the power progression and hence on the deterioration of CPR is presented. Another important procedure to limit the effects of the transient is the fast shutdown of the reactor that is caused when the reactor power reaches the limit value. It is shown that the SCRAM is not fast enough to reduce the first power maximum, but is able to prevent the appearance of a second — much smaller — maximum that would occur around one second after the first one in the absence of a SCRAM.

Experiments on Gas Mixing and Stratification Driven by Jets and Plumes in Large-Scale, Multi-Compartment Geometries

The present paper discusses experiments carried out to examine mixing of different gases (steam, air) and the evolution their distributions in large-scale, multi compartment geometry imitating nuclear reactor containment compartments. The flow and the mixing process in the experiments are driven by plumes and jets representing source structures with different momentum-to-buoyancy strength. The time evolution of the relevant parameters like gas concentrations, velocities and temperatures are followed using dedicated instrumentation. The data obtained is meant to be used for the validation and development of high-resolution, mainly CFD based, 3D computational tools for nuclear reactor containment safety analysis.

Natural Convection Heat Transfer in a Rectangular Liquid Metal Pool With Bottom Heating and Top Cooling

An experimental study is performed to investigate the natural convection heat transfer characteristics with subcooled coolant to create engineering database for basic applications in a lead alloy cooled reactor. Tests are performed in the ALTOS (Applied Liquid-metal Thermal Operation Study) apparatus as part of MITHOS (Metal Integrated Thermo Hydrodynamic Operation System). A relationship is determined between the Nusselt number Nu and the Rayleigh number Ra in the liquid metal rectangular pool. Results are compared with correlations and experimental data in the literature. Given the similar Ra condition, the present test results for Nu of the liquid metal pool with top subcooling are found to be similar to those predicted by the existing correlations or experiments. The current test results are utilized to develop natural convection heat transfer correlations applicable to low Prandtl number Pr fluids that are heated from below and cooled by the external coolant above. Results from this study are slated to be used in designing BORIS (Battery Optimized Reactor Integral System), a small lead cooled modular fast reactor for deployment at remote sites cycled with MOBIS (Modular Optimized Brayton Integral System) for electricity generation, tied with NAVIS (Naval Application Vessel Integral System) for ship propulsion, joined with THAIS (Thermochemical Hydrogen Acquisition Integral System) for hydrogen production, and coupled with DORIS (Desalination Optimized Reactor Integral System) for seawater desalination. Tests are performed with Wood’s metal (Pb-Bi-Sn-Cd) filling a rectangular pool whose lower surface is heated and upper surface cooled by forced convection of water. The test section is 20 cm long, 11.3 cm high and 15 cm wide. The simulant has a melting temperature of 78°C. The constant temperature and heat flux condition was realized for the bottom heating once the steady state had been met. The test parameters include the heated bottom surface temperature of the liquid metal pool, the input power to the bottom surface of the section, and the coolant temperature.

ASTEC: An Integral Code for Simulation of Severe Light Water Reactor Accidents

The integral code ASTEC (Accident Source Term Evaluation Code) commonly developed by IRSN and GRS is a fast running programme, which allows the calculation of entire sequences of severe accidents (SA) in light water reactors from the initiating event up to the release of fission products into the environment, thereby covering all important in-vessel and containment phenomena. Thus, the main fields of ASTEC application are intended to be accident sequence studies, uncertainty and sensitivity studies, probabilistic safety analysis level 2 studies as well as support to experiments. The modular structure of ASTEC allows running each module independently and separately, e.g. for separate effects analyses, as well as a combination of multiple modules for coupled effects testing and integral analyses. Among activities concentrating on the validation of individual ASTEC modules describing specific phenomena, the applicability to reactor cases marks an important step in the development of the code. Feasibility studies on plant applications have been performed for several reactor types such as the German Konvoi PWR 1300, the French PWR 900, and the Russian VVER-1000 and −440 with sequences like station blackout, small- or medium-break loss-of-coolant accident, and loss-of-feedwater transients. Subject of this paper is a short overview on the ASTEC code system and its current status with view to the application to severe accidents sequences at several PWRs, exemplified by selected calculations.

Experiment and Numerical Simulation of Bubble Behaviors in Argon Gas Injection Into Lead-Bismuth Pool

In a lead-bismuth alloy (45%Pb-55%Bi) cooled direct contact boiling water fast reactor (PBWFR), steam can be produced by direct contact of feed water with primary Pb-Bi coolant in the upper core plenum, and Pb-Bi coolant can be circulated by buoyancy forces of steam bubbles. As a basic study to investigate the two-phase flow characteristics in the chimneys of PBWFR, a two-dimensional two-phase flow was simulated by injecting argon gas into Pb-Bi pool in a rectangular vessel (400mm in length, 1500mm in height, 50mm in width), and bubble behaviors were investigated experimentally. Bubble sizes, bubble rising velocities and void fractions were measured using void probes. Argon gas was injected through five nozzles of 4mm in diameter into Pb-Bi at two locations. The experimental conditions are the pressure of atmospheric pressure, Pb-Bi temperatures of 443K, and the flow rate of injection Ar gas is 10, 20, and 30 NL/min. The measured bubble rising velocities were distributed in the range from 1 to 3 m/s. The average velocity was about 0.6 m/s. The measured bubble chord lengths were distributed from 1mm up to 30mm. The average chord length was about 7mm. An analysis was performed by two-dimensional and two-fluid model. The experimental results were compared with the analytical results to evaluate the validity of the analytical model. Although large diameter bubbles were observed in the experiment, the drag force model for spherical bubbles performed better for simulation of the experimental result because of high surface tension force of Pb-Bi.

Large Break LOCA Accident Management Strategies for Accidents With Large Containment Leaks

The goal of this work is the investigation of the influence of different accident management strategies on the thermal-hydraulics in the containment during a Large Break Loss of Coolant Accident with a large containment leak from the beginning of the accident. The increasing relevance of terrorism suggests a closer look at this kind of severe accidents. Normally the course of severe accidents and their associated phenomena are investigated with the assumption of an intact containment from the beginning of the accident. This intact containment has the ability to retain a large part of the radioactive inventory. In these cases there is only a release via a very small leakage due to the untightness of the containment up to cavity bottom melt through. This paper represents the last part of a comprehensive study on the influence of accident management strategies on the source term of VVER-1000 reactors. Basically two different accident sequences were investigated: the “Station Blackout”-sequence and the “Large Break LOCA”. In a first step the source term calculations were performed assuming an intact containment from the beginning of the accident and no accident management action. In a further step the influence of different accident management strategies was studied. The last part of the project was a repetition of the calculations with the assumption of a damaged containment from the beginning of the accident. This paper concentrates on the last step in the case of a Large Break LOCA. To be able to compare the results with calculations performed years ago the calculations were performed using the Source Term Code Package (STCP), hydrogen explosions are not considered. In this study four different scenarios have been investigated. The main parameter was the switch on time of the spray systems. One of the results is the influence of different accident management strategies on the source term. In the comparison with the sequence with intact containment it was demonstrated that the accident management measures have quite lower consequences. In addition it was shown that in the case of a “Large Break LOCA”-sequence the intact containment retains the nuclides up to a factor of 20 000. This is much more than in the case of a “Station Blackout”-sequence. Within the frame of the study 17 source terms have been generated to evaluate in detail accident management strategies for VVER-1000 reactors.

Uncertainty Analysis for a Depressurised Loss of Forced Cooling Event of the PBMR Reactor

This paper presents an uncertainty analysis for a Depressurised Loss of Forced Cooling (DLOFC) event that was performed with the systems CFD (Computational Fluid Dynamics) code Flownex for the PBMR reactor. An uncertainty analysis was performed to determine the variation in maximum fuel, core barrel and reactor pressure vessel (RPV) temperature due to variations in model input parameters. Some of the input parameters that were varied are: thermo-physical properties of helium and the various solid materials, decay heat, neutron and gamma heating, pebble bed pressure loss, pebble bed Nusselt number and pebble bed bypass flows. The Flownex model of the PBMR reactor is a 2-dimensional axi-symmetrical model. It is simplified in terms of geometry and some other input values. However, it is believed that the model adequately indicates the effect of changes in certain input parameters on the fuel temperature and other components during a DLOFC event. Firstly, a sensitivity study was performed where input variables were varied individually according to predefined uncertainty ranges and the results were sorted according to the effect on maximum fuel temperature. In the sensitivity study, only seven variables had a significant effect on the maximum fuel temperature (greater that 5°C). The most significant are power distribution profile, decay heat, reflector properties and effective pebble bed conductivity. Secondly, Monte Carlo analyses were performed in which twenty variables were varied simultaneously within predefined uncertainty ranges. For a one-tailed 95% confidence level, the conservatism that should be added to the best estimate calculation of the maximum fuel temperature for a DLOFC was determined as 53°C. This value will probably increase after some model refinements in the future. Flownex was found to be a valuable tool for uncertainly analyses, facilitating both sensitivity studies and Monte Carlo analyses.

Experiments on Damping of Acoustic Waves in Water-Filled Pipes

A simple, fundamental experimental study was conducted to better understand acoustic wave propagation is fluid-filled pipes. Three experiments were undertaken: the first with zero flow and a closed outlet end, the second with turbulent flow and an open outlet end and the third with zero flow and an open outlet end. The intent was to obtain data for model comparison and to determine the effect of turbulent flow on the system response. New insights are obtained and reported.

Influence of Test Tube Material on Subcooled Flow Boiling Critical Heat Flux in Short Vertical Tube

The steady state subcooled flow boiling critical heat flux (CHF) for the flow velocities (u = 4.0 to 13.3 m/s), the inlet subcoolings (ΔTsub,in = 48.6 to 154.7 K), the inlet pressure (Pin = 735.2 to 969.0 kPa) and the increasing heat input (Q0exp(t/τ), τ = 10, 20 and 33.3 s) are systematically measured with the experimental water loop. The 304 Stainless Steel (SUS304) test tubes of inner diameters (d = 6 mm), heated lengths (L = 66 mm) and L/d = 11 with the inner surface of rough finished (Surface roughness, Ra = 3.18 μm), the Cupro Nickel (Cu-Ni 30%) test tubes of d = 6 mm, L = 60 mm and L/d = 10 with Ra = 0.18 μm and the Platinum (Pt) test tubes of d = 3 and 6 mm, L = 66.5 and 69.6 mm, and L/d = 22.2 and 11.6 respectively with Ra = 0.45 μm are used in this work. The CHF data for the SUS304, Cu-Ni 30% and Pt test tubes were compared with SUS304 ones for the wide ranges of d and L/d previously obtained and the values calculated by the authors’ published steady state CHF correlations against outlet and inlet subcoolings. The influence of the test tube material on CHF is investigated into details and the dominant mechanism of subcooled flow boiling critical heat flux is discussed.

Measurements of Sheath Temperature Profiles in Bruce LVRF Bundles Under Post-Dryout Heat Transfer Conditions in Freon

An experimental program has been completed to study the behaviour of sheath wall temperatures in the Bruce Power Station Low Void Reactivity Fuel (shortened hereafter to Bruce LVRF) bundles under post-dryout (PDO) heat-transfer conditions. The experiment was conducted with an electrically heated simulator of a string of nine Bruce LVRF bundles, installed in the MR-3 Freon heat transfer loop at the Chalk River Laboratories (CRL), Atomic Energy of Canada Limited (AECL). The loop used Freon R-134a as a coolant to simulate typical flow conditions in CANDU® nuclear power stations. The simulator had an axially uniform heat flux profile. Two radial heat flux profiles were tested: a fresh Bruce LVRF profile and a fresh natural uranium (NU) profile. For a given set of flow conditions, the channel power was set above the critical power to achieve dryout, while heater-element wall temperatures were recorded at various overpower levels using sliding thermocouples. The maximum experimental overpower achieved was 64%. For the conditions tested, the results showed that initial dryout occurred at an inner-ring element at low flows and an outer-ring element facing internal subchannels at high flows. Dry-patches (regions of dryout) spread with increasing channel power; maximum wall temperatures were observed at the downstream end of the simulator, and immediately upstream of the mid-bundle spacer plane. In general, maximum wall temperatures were observed at the outer-ring elements facing the internal subchannels. The maximum water-equivalent temperature obtained in the test, at an overpower level of 64%, was significantly below the acceptable maximum temperature, indicating that the integrity of the Bruce LVRF will be maintained at PDO conditions. Therefore, the Bruce LVRF exhibits good PDO heat transfer performance.