Reactivity insertion limits in a typical pool-type research reactor cooled by natural circulation

The course of reactivity insertion in a pool type research reactor, with scram disabled under natural circulation condition is numerically investigated. The analyses were performed by a coupled kinetic–thermal–hydraulic computer code developed specifically for this task. The 10-MW IAEA MTR research reactor was subjected to unprotected reactivity insertion (step and ramp) for both low and high-enriched fuel with continuous reactivity feedback due to coolant and fuel temperature effects. In general, it was found that the power, core mass flow rate and clad temperature under fully established natural circulation are higher for high-enriched fuel than for low enriched fuel. This is unlike the case of decay heat removal, where equal clad temperatures are reported for both fuels. The analysis of reactivity represented by the maximum insertion of positive reactivity ($0.73) demonstrated the high inherent safety features of MTR-type research reactor. Even in the case of total excess reactivity without scram, the high reactivity feedbacks of fuel and moderator temperatures limit the power excursion and avoid consequently escalation of clad temperature to the level of onset of nucleate boiling and sub-cooled void formation. The code can also be modified to provide an accurate capability for the analyses of research reactor transients under forced convection.     

Estimation of fuel and clad temperature of a research reactor during dry period of de-fuelling operation

Heat transfer and fluid flow studies related to spent fuel bundle of a research reactor in fuelling machine has been carried out. When the fuel is in reactor core, the heat generated in the fuel bundle is removed by heavy water under normal reactor operation. However, during the de-fuelling operation, the fuel bundle is exposed to air for some period called dry period. During this period, the decay heat from fuel bundle has to be removed by air flow. This flow of air is induced by natural convection only. In this period, the temperatures of fuel and clad rise. If clad temperature rises beyond a certain limit, structural failure may occur. This failure can result into release of fission products from fuel rod. Hence the temperature of clad has to be within specified limit under all conditions. The objective of this study is to estimate the clad temperature rise during the dry period.In the CFD simulation, the turbulent natural convection flow over fuel and radiation heat transfer are accounted. Standard k-? model for turbulence, Boussinesq approximation for computing the natural convection flow and IMMERSOL model for radiation are used.The steady state and transient CFD simulation of flow and heat is performed, using the CFD code PHOENICS. The steady state analysis provides the maximum temperature the clad will attain if fuel bundle is left exposed to air for sufficiently long time. For safe operation, the clad temperature should be limited to a specified value. From steady state CFD analysis, it is found that steady state clad temperature for various decay powers is higher than the limiting value. Hence transient analysis is also performed. In the transient analysis, the variation of clad temperature with time is predicted for various decay powers. Safe dry time, i.e. the time required for clad to reach the limiting value, is predicted for various decay powers. Determination of safe dry time helps in deciding the time available to the operator to drop the bundle in light water pool for storage. The analysis is found useful in optimizing the de-fuelling process.     

Transient experiments on fast reactor core thermal-hydraulics and its numerical analysis: Inter-subassembly heat transfer and inter-wrapper flow under natural circulation conditions

Sodium experiments were conducted on core thermal-hydraulics simulating a scram transient of a large scale fast breeder reactor using the test facility PLANDTL-DHX with seven fuel subassemblies. The influence of inter-subassembly heat transfer on temperature distribution in the subassembly was revealed via measurements. The flow in the gap between neighboring subassemblies called inter-wrapper flow (IWF) was also studied in relation to its capability of cooling the subassemblies. A computational model is presented for predicting the transient without IWF. The multi-dimensional numerical analysis model employs an empirical correlation to simulate mixing effects between adjacent subchannels. It was shown that the present computational method could evaluate the transient behavior of thermal-hydraulics in the subassemblies accurately from forced to natural circulation accompanied by inter-subassembly heat transfer and flow redistribution in the subassembly. The cooling effects of IWF on the fuel subassemblies were found in spite of natural circulation flow reduction in the primary loop attributable to temperature decreases in the upper non-heated section in the core. The inter-wrapper flow can effectively cool the core under extreme conditions of low flow rates through the core.     

Analysis of reactivity induced accidents at Pakistan Research Reactor-1

Analysis of reactivity induced accidents in Pakistan Research Reactor-1 (PARR-1) utilizing low enriched uranium (LEU) fuel, has been carried out using standard computer code PARET. The present core comprises of 29 standard and five control fuel elements. Various modes of reactivity insertions have been considered. The events studied include: start-up accident; accidental drop of a fuel element on the core; flooding of a beam tube with water; removal of an in-pile experiment during reactor operation etc. For each of these transients, time histories of reactor power, energy released and clad surface temperature etc. were calculated. The results reveal that the peak clad temperatures remain well below the clad melting temperature during these accidents. It is concluded that the reactor, which is operated safely at a steady-state power level of 10 MW, with coolant flow rate of 950 m³/h, will also be safe against any possible reactivity induced accident and will not result in a fuel failure.     

PASCAR: Long burning small modular reactor based on natural circulation

PASCAR is a 100 MWt/35 MWe lead-bismuth-cooled small modular reactor which requires no on-site refueling and well suits to be used as a distributed power source in either a single unit or a cluster for electricity, heat supply, and desalination. This paper includes both steady-state and transient performance evaluations for neutronics and thermal-hydraulics. Through design optimization studies for minimizing a burn-up reactivity loss, the metallic fuels-loaded core was designed with less than 1$ reactivity swing over 20-year cycle. A radial peaking power location shows the slow inward migration from outer enrichment zones while maintaining peaking factor within 1.35, reducing radiation damage and corrosion duty of high temperature environments. Equipped with coolant flow path large enough to ensure low pressure drop, this reactor is intended to operate by only natural circulation of chemically inert coolant within relatively low temperature range, 320-420 °C. Peak outlet temperature is nearly 450 °C where an Al-containing duplex cladding has sufficient corrosion resistance. Despite of 50% decrease of fuel thermal conductivity after swelling, inherent negative reactivity feedback and passive decay heat removal capability could secure an ample safety margin of peak fuel centerline temperature in tow safety analyses, unprotected transient overpower and unprotected loss of heat sink. The likelihood of loss of coolant, loss of flow, and local blockage is virtually eliminated by employing respectively a double-walled vessel, pump-less cooling, and cross-flow allowed open square assemblies. Simple fabrication, modular construction, and long burning cycle would compensate for economic disadvantages over smaller power and lower temperature than those of conventional fast reactors.     

CFD investigation of flow inversion in typical MTR research reactor undergoing thermal–hydraulic transients

Three dimensional CFD full simulations of the fast loss of flow accident (FLOFA) of the IAEA 10 MW generic MTR research reactor are conducted. In this system the flow is initially downward. The transient scenario starts when the pump coasts down exponentially with a time constant of 1 s. As a result the temperatures of the heating element, the clad, and the coolant rise. When the flow reaches 85% of its nominal value the control rod system scrams and the power drops sharply resulting in the temperatures of the different components to drop. As the coolant flow continues to drop, the decay heat causes the temperatures to increase at a slower rate in the beginning. When the flow becomes laminar, the rate of temperature increase becomes larger and when the pumps completely stop a flow inversion occurs because of natural convection. The temperature will continue to rise at even higher rates until natural convection is established, that is when the temperatures settle off. The interesting 3D patterns of the flow during the inversion process are shown and investigated. The temperature history is also reported and is compared with those estimated by one-dimensional codes. Generally, very good agreement is achieved which provides confidence in the modeling approach.     

Investigation of thermohydraulic parameters during natural convection cooling of TRIGA reactor

Important steady-state thermohydraulic parameters of the TRIGA research reactor operating under natural convection mode of coolant flow were investigated using NCTRIGA computer code. Neutronic parameters used in preparing the input of NCTRIGA were taken from the analysis performed by 3-D Monte Carlo code MCNP4C. Benchmarking of the NCTRIGA calculated results were performed against the experimental data measured by the thermocouples in the instrumented fuel element (IFE) during the steady state operation of the reactor under natural convection mode of coolant flow. Various thermohydraulic parameters like the coolant velocity, flow rate and mass flow rate were generated for the hot channel as well as for the two channels comprising instrumented fuels. Calculated peak fuel temperatures at different power levels were compared with the measured values and also with the calculations performed by PARET code. Axial temperature profile at the fuel centreline, fuel surface and coolant in the hot channel were generated. Fuel surface heat flux, heat transfer coefficient and Reynolds’s number for the hot channel were also calculated. The effect of fuel-cladding gap and the influence of fuel rod spacing were investigated to validate the performance of NCTRIGA code. The investigated results were found to be in good agreement with the experimental values, which indicates that the NCTRIGA code can be used with confidence for TRIGA reactor analysis.     

Development of a coupled neutronic/thermal-hydraulic tool with multi-scale capabilities and applications to HPLWR core analysis

The High Performance Light Water Reactor (HPLWR) is a thermal spectrum nuclear reactor cooled and moderated with light water operated at supercritical pressure. It is an innovative reactor concept, which requires developing and applying advanced analysis tools as described in the paper.The relevant water density reduction associated with the heat-up, together with the multi-pass core design, results in a pronounced coupling between neutronic and thermal-hydraulic analyses, which takes into account the strong natural influence of the in-core distribution of power generation and water properties. The neutron flux gradients within the multi-pass core, together with the pronounced dependence of water properties on the temperature, require to consider a fine spatial resolution in which the individual fuel pins are resolved to provide precise evaluation of the clad temperature, currently considered as one of the crucial design criteria. These goals have been achieved considering an advanced analysis method based on the usage of existing codes which have been coupled with developed interfaces.Initially neutronic and thermal-hydraulic full core calculations have been iterated until a consistent solution is found to determine the steady state full power condition of the HPLWR core. Results of few group neutronic analyses might be less reliable in case of HPLWR 3-pass core than for conventional LWRs because of considerable changes of the neutron spectrum within the core, hence 40 groups transport theory has been preferred to the usual 2 groups diffusion theory. Successively, with the usage of a developed pin-power reconstruction technique capable to account for the innovative fuel assembly design, sub-channel investigations of the individual fuel assemblies have been performed evaluating pin-wise clad temperatures. Obtained results will be discussed giving a detailed insight of the revolutionary HPLWR 3 pass core concept and understanding the physical reasons, which influence the local clad temperatures.The obtained results represent a new quality in core analyses, which takes into full consideration the coupling between neutronics and thermal-hydraulics as well as the spatial effects of the fuel assembly heterogeneity in determining the local pin-power and the associated maximum clad temperature.     

Uncontrolled withdrawal of a control rod without scram

In the present work the thermal-hydraulics of reactivity-induced transients in low enriched uranium (LEU) core of a typical material test research reactor (MTR) are analyzed using the previous program developed by Khater et al. The analysis was done for uncontrolled withdrawal of a control rod with scram-disabled conditions. Initiating reactivity events with and without the influence of reactivity efficiency curve (“S” curve) were considered. The results of the proposed transients are analyzed and compared with each other. In transient without the “S” curve influence, a high primary peak power of 406.18 MW is attained and a clad melt down takes place after 1.85 s. In the transient with the “S” curve influence, a high super prompt-critical situation is produced (1.762$ at 0.895 s) with a very high primary peak power of 801.05 MW at 0.912 s. Also, a fast clad melt down is resulted in the hot channel at 1.088 s and a stable film boiling is established. This study indicates that, compared to the application of linear reactivity curve, the application of the reactivity efficiency curve results in the prediction of higher peaks in power and temperatures (fuel, clad and coolant) with a fast clad melt down.     

Thermal-hydraulic behavior of the NUR nuclear research reactor during a fast loss of flow transient

Nuclear safety analysis remains of crucial importance for both the design and the operation of nuclear reactors. Safety analysis usually entails the simulation of several selected postulated accidents, which can be divided into two main categories, namely reactivity insertion accident (RIA) and loss of flow accident (LOFA). In this paper, thermal-hydraulic simulations of fast LOFA accident were carried out on the new core configuration of the material test research reactor NUR. For this purpose, the nuclear reactor analysis PARET code was used to determine the reactor performance by calculating the reactor power, the reactivity and the temperatures of different components (fuel, clad and coolant) as a function of time. It was observed that during the transient the maximum clad temperature remained well below the critical temperature limit of 110 °C, and the maximum coolant temperature did not exceed the onset of nucleate boiling point of 120 °C. It is concluded that the reactor can be operated at full power level with sufficient safety margins with regard to such kind of transients.     

Sensitivity of reactivity insertion limits with respect to safety parameters in a typical MTR

The sensitivity of various safety parameters, affecting the reactivity insertion limits imposed by clad melting temperature for a typical pool type research reactor, have been investigated in this work. The analysis was done for low enriched uranium (LEU) core with scram disabled conditions. The temperature coefficients of fuel and coolant, void/density coefficient and β_eff were individually varied and the reactor behavior for different ramp reactivity transients was studied. In this work ramp reactivity insertions from 1.6 to 2 $/0.5 s were selected and peak power, maximum fuel, clad and coolant temperatures were determined. Results show that peak power decreases with an increase in the Doppler coefficient of reactivity. However, it rises with an increase in the reactivity insertion. Core remains insensitive to the coolant temperature coefficient of reactivity for ramps in the range of 1.6–1.9/0.5 s. Peak power decreases with an increase in the void coefficient of reactivity (0.1 $/%void to 0.8 $/%void). With a decrease in the void coefficient of reactivity, the maximum fuel and clad temperatures show a non-linear rise. Power and temperature peaks in the transient are sensitive to the values of β_eff. Finally, it can be concluded that LEU is a safe core due to its smaller β_eff, larger Doppler coefficient and void coefficient of reactivity. It is inferred through this work that reactivity insertion limits of LEU core are quite insensitive to β_eff, the Doppler coefficient and the coolant temperature coefficient of reactivity. They are highly sensitive to the change of the void coefficient of reactivity in the core.     

Analysis of loss of flow accident at Pakistan research reactor-1

The main objective of the reactor safety is to keep the reactor core in a condition, which does not permit any release of radioactivity into the environment. In order to ensure this, the reactor must have sufficient safety margins during all possible operational conditions (normal as well as accidental). To accomplish this, a study has been carried out, for the analysis of loss of flow accident (LOFA), which is one of the probable scenarios among other possible events such as reactivity-induced-accidents, loss of coolant accident, etc. The study has been carried out for Pakistan research reactor, PARR-1, which was initially converted from HEU to LEU fuel. It is a swimming pool type reactor using MTR type fuel. Presently, a new core is proposed to be assembled containing LEU and some of the used (less burnt) HEU fuel elements. The accident is assumed when the reactor is running at a steady-state power level of 9.8 MW. Computer code PARET and standard correlations were employed to compute various parameters. Results predict nucleate boiling in the core but the temperatures would remain far below the fuel clad melting point.     

Effects of high density dispersion fuel loading on the dynamics of a low enriched uranium fueled material test research reactor

The effects of using high density low enriched uranium on the dynamics of a material test research reactor were studied. For this purpose, the low density LEU fuel of an MTR was replaced with high density LEU fuels currently being developed under the RERTR program. Since the alloying elements have different properties affecting the reactor in different ways, fuels U–Mo (9w/o) which contain the same elements in same ratio were selected for analysis. Simulations were carried out to determine the reactor performance under reactivity insertion and loss of flow transients. Nuclear reactor analysis code PARET was employed to carry out these calculations. It is observed that during the fast reactivity insertion transient, the maximum reactor power is achieved and the energy released till the power reaches its maximum increases by 45% and 18.5%, respectively, as uranium density increases from 6.57 gU/cm³ to 8.90 gU/cm³. This results in increased maximum temperatures of fuel, clad and coolant outlet, achieved during the transient, by 27.7 K, 19.7 K and 7.9 K, respectively. The time required to reach the peak power decreases. During the slow reactivity insertion transient, the maximum reactor power achieved increases slightly by 0.3% as uranium density increases from 6.57 gU/cm³ to 8.90 gU/cm³ but the energy generated till the power reaches its maximum decreases by 5.7%. The temperatures of fuel, clad and coolant outlet remain almost the same for all types of fuels. During the loss of flow transients, no appreciable difference in the power and temperature profiles was observed and the graph plots overlapped each other.     

Prediction, analysis and solution of the flow inversion phenomenon in a typical MTR-reactor with upward core cooling

Research reactors of power greater than 20 MW are usually designed to be cooled with upward coolant flow direction inside the reactor core. This is mainly to prevent flow inversion problems following a pump coast down. However, in some designs and under certain operating conditions, flow inversion phenomenon is predicted. In the present work, the best-estimate Material Testing Reactors Thermal-Hydraulic Analysis program (MTRTHA) is used to simulate a typical MTR reactor behavior with upward cooling. The MTRTHA model consists of five interactively coupled submodels for: (a) coolant, (b) fuel plate, (c) chimney, lower plenum, suction box and cold leg, (d) flap valve and (e) natural circulation flow. The model divides the active core into a specified axial regions and the fuel plate into a specified radial zones, then a nodal calculation is performed for both average and hot channels with a chopped cosine shaped heat generation flux. The reactor simulation under loss of off-site power is performed for two cases namely: two-flap valves open and one flap-valve fails to open. The simulation is performed under a hypothetical case of loss of off-site power. Unfortunately, the flow inversion phenomenon is predicted under certain decay heat and/or pool temperature values below the design values. In most cases, the flow inversion phenomenon is accompanied by boiling which is undesirable phenomenon in this type of reactors as it could affect the fuel-clad integrity. The model results for the flow inversion phenomenon prediction are analyzed and a solution of the problem is suggested.     

Adsorptive Characteristics of Fibrous Adsorbent for Uranium

Transient thermal-hydraulic behaviors of the JRR-3 which is an open-pool type research reactor has been analyzed with the THYDE-P1 code. The focal point is the thermal-hydraulic behaviors related to the core flow reversal during the transition from forced circulation downflow to natural circulation upflow. In the case of a loss-of-coolant accident (LOCA), for example, the core flow reversal is expected to occur just after the water pool isolation from the primary cooling loop with a leak. The core flow reversal should cause a sudden increase in fuel temperature and a steep decrease in the departure-from-nucleate-boiling ratio (DNBR) and the phenomenon is, therefore, very important especially for safety design and evaluation of research reactors.

Major purposes of the present work are to clarify physical phenomena during the transient and to identify important parameters affecting the peak fuel temperature and the minimum DNBR. The results calculated with THYDE-P1 assuming the sequences of events of the loss-of-offsite power and LOCA help us to understand the phenomena both qualitatively and quantitatively, with respect to the safety design and evaluation.     

CFD simulation of the IAEA 10 MW generic MTR reactor under loss of flow transient

Three-dimensional simulation of the IAEA 10 MW generic reactor under loss of flow transient is introduced using the CFD code, Fluent. The IAEA reactor calculation is a safety-related benchmark problem for an idealized material testing reactor (MTR) pool type specified in order to compare calculational methods used in various research centers. The flow transients considered include fast loss of flow accidents (FLOFA) and slow loss of flow accidents (SLOFA) modeled with exponential flow decay and time constants of 1 and 25 s, respectively. The transients were initiated from a power of 12 MW with a flow trip point at 85% nominal flow and a 200 ms time delay. The simulation shows comparable results as those published by other research groups. However, interesting 3D patterns are shown that are usually lost based on the one-dimensional simulations that other research groups have introduced. In addition, information about the maximum clad surface temperature, the maximum fuel element temperature as well as the location of hot spots in fuel channel is also reported.     

Expansion of the method of correlated temperatures to clad and fuel design analyses

In the thermal design of nuclear reactor cores, specified design limits (temperatures and linear power rating) should not be exceeded by the operating values of certain elements (coolant, clad and fuel). However, a certain number of channels or fuel pins could be permitted to exceed the specified limits without affecting the reactor's safety while still allowing reliable operation. An expansion of the method of correlated temperatures, developed for coolant temperature analysis, was performed to enable clad temperature and fuel centerline melting analyses for reactor core reliability studies. Since generation of random numbers is involved, calculational procedures, tailored to designer needs, were developed in order to reduce computational time. The method is applied to a typical LMFBR core and results are presented for various assumed clad and fuel design limits.     

Thermal hydraulic and safety analysis for core conversion (HEU–LEU) of Syrian Miniature Neutron Source Reactor

The paper presents the behavior and properties analysis of the low enriched uranium fuel compared with the original high enriched uranium fuel. The MNSR reactor core was modeled with both fuel materials and the reactor behavior was studied during the steady state and abnormal conditions. The MERSAT code was used in the analysis. The steady state thermal hydraulic analysis results were compared with that obtained from the experimental results hold during commissioning the Syrian MNSR. Comparison with experimental data shows that the steady-state behavior of the HEU core was accurately predicted by the MERSAT code calculations. The validated model was then used to analyze LEU cores with two proposed UO₂ fuel pin designs. With each LEU core, the steady state and 3.77 mk rod withdrawal transient were run and the results were compared with the available published data in the literatures for the low enriched uranium fuel core. The results reveal that the low enriched uranium fuel showed a good behavior and the peak clad temperatures remain well below the clad melting temperature during reactivity insertion accident.     

Simulation of uncontrolled loss of flow transients of a material test research reactor fuelled with low and high enriched uranium dispersion fuels

The effects of using low and high enrichment uranium fuel on the uncontrolled loss of flow transients in a material test research reactor were studied. For this purpose, simulations were carried out of an MTR fuelled separately with LEU and HEU fuel, to determine the reactor performance under loss of flow transients with totally failed external control systems. The coolant pump was assumed to loose its performance and the coolant flow rate reduced according to the relation m(t)/m₀ = exp(−t/25) to a new stable level. The new reduced flows m/m₀ = 0.2, 0.4, 0.6 and 0.8 were modeled. The nuclear reactor analysis code PARET was employed to carry out these calculations. It was observed that the reactors stabilized at new power levels which were lower than the original power level, with the power of HEU fuelled reactor slightly lesser than that of the LEU fuelled reactor. However, at the start of transient, the LEU fuelled reactor had a lower power level resulting in lower fuel, clad and coolant temperatures than the HEU fuelled reactor.     

液体燃料熔盐堆三维稳态分析程序开发

液体燃料熔盐堆的物理热工特性与固体燃料反应堆有很大的不同,在分析计算中必须考虑燃料流动特性的影响,一般分析固体反应堆的程序均不能直接用于分析液体燃料熔盐堆。根据熔盐堆的流动特性,建立了液体燃料熔盐堆的三维中子动力学模型和流动传热模型,开发了针对液体燃料熔盐堆的三维稳态核热耦合程序,并以此分析了稳态情况下MOSART堆的物理热工特性。结果表明,堆芯流速对快中子和热中子影响较小,对堆芯温度和缓发中子分布影响较大。