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.     

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.     

Safety analysis of the reference research reactor MTR during reactivity insertion accident using the code MERSAT

The IAEA’s reference research reactor MTR-10 MW has been modeled using the code MERSAT. The developed MERSAT model consists of detailed representation of primary and secondary loops including reactor pool, bypass, main pump, heat exchanger and reactor core with the corresponding neutronics and thermalhydraulic characteristics. Following the successful accomplishment of the steady state operation at nominal power of 10 MW, reactivity insertion accident (RIA) for three different initial reactivity values of $1.5/0.5 s, $1.35/0.5 s and $0.1/1.0 s have been simulated. The predicted peaks of reactor power, hot channel fuel, clad and coolant temperatures demonstrate inherent safety features of the reference MTR reactor. Only in case of the fast RIA of $1.5/0.5 s, where the peak power of 133.66 MW arrived 0.625 s after the start of the transient, the maximum hot channel clad temperature arrives at the condition of subcooled boiling with the subsequent void formation. However, due to the strong negative reactivity feedback effects of coolant and fuel temperatures the void formation persists for a very short time so that thermalhydraulic conditions remained far from exceeding the safety design limits of thermalhydraulic instability and DNB. Finally, the simulation results show good agreement with previous international benchmark analyses accomplished with other qualified channel and thermalhydraulic system codes.     

Reactivity feedbacks of a material test research reactor fueled with various low enriched uranium dispersion fuels

The reactivity feedbacks of a material test research reactor using various low enriched uranium fuels, having same uranium density were calculated. For this purpose, the original aluminide fuel (UAl_x–Al) containing 4.40 gU/cm³ of an MTR was replaced with silicide (U₃Si–Al and U₃Si₂–Al) and oxide (U₃O₈–Al) dispersion fuels having the same uranium density as of the original fuel. Calculations were carried out to find the fuel temperature reactivity feedback, moderator temperature reactivity feedback, moderator density reactivity feedback and moderator void reactivity feedback. Nuclear reactor analysis codes including WIMS-D4 and CITATION were employed to carry out these calculations. It was observed that the magnitudes all the respective reactivity feedbacks from 38 °C to 50 °C and 100 °C, at the beginning of life, of all the fuels were very close to each other. The fuel temperature reactivity feedback of the U₃O₈–Al was about 2% more than the original UAl_x–Al fuel. The magnitudes of the moderator temperature, moderator density and moderator void reactivity feedbacks of all the fuels, showed very minor variations from the original aluminide fuel.     

Assessment of fuel conversion from HEU to LEU in the Syrian MNSR reactor using the MCNP code

Assessment of fuel conversion from high enriched uranium (HEU) to low enriched uranium (LEU) fuel in the Syrian MNSR reactor was conducted in this paper. Three 3-D neutronic models for the Syrian MNSR reactor using the MCNP-4C code were developed to assess the possibility of fuel conversion from 89.87% HEU fuel (UAl₄–Al) to 19.75% LEU fuel (UO₂). The first model showed that 347 fuel rods with HEU fuel were required to obtain a reactor core with 5.17 mk unadjusted excess reactivity. The second model showed that only 200 LEU fuel rods distributed in the reactor core like the David star figure were required to obtain a reactor core with 4.85 mk unadjusted excess reactivity. The control rod worth using the LEU fuel was enhanced. Finally, the third model showed that distribution of 200 LEU fuel rods isotropically in the 10 circles of the reactor core failed to convert the fuel since the calculated core unadjusted excess reactivity for this model was 10.45 mk. This value was far beyond the reactor operation limits and highly exceeded the current MNSR core unadjusted excess reactivity (5.17 mk).     

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.     

Burnup study for Pakistan Research Reactor-1 utilizing high density low enriched uranium fuel

Burnup study for Pakistan Research Reactor-1 (PARR-1), which is a typical swimming pool type MTR utilizing high density low enriched uranium fuel, was performed by using Fuel Cycle Analysis Program (FCAP). Existing equilibrium core of PARR-1, which is relatively economical but provides less neutron fluxes per unit power than the first equilibrium core, was formed by adding five more fuel elements in the first equilibrium core. This study shows that if the fuel loading is increased in the first equilibrium core of PARR-1 by replacing the fuel of density 3.28 gU/cm³ by the fuel of density 4.00 gU/cm³ then the new equilibrium core can provide 10% higher neutron fluxes at the irradiation sites and will also require 1.5 kg less fuel than that required for existing equilibrium core for one-year full power operation at 10 MW. The new core provides neutron fluxes at 13% lower cost and if the size of this core is further reduced by three fuel elements then this core can provide 20% higher thermal neutron flux at the central flux trap at 9% lower cost. A possible use of U-Mo (5 w/o Mo) fuel of density 8.5 gU/cm³ in PARR-1 with an increase in existing water channel width from 2.1 to 2.45 mm (Ann. Nucl. Energy 32(1), 29–62) would provide up to 41% more thermal neutron flux at the central flux trap at 13% lower cost than the existing equilibrium core. The power peaking factors in these cores are similar to the power peaking factors of the existing equilibrium core and these cores are likely to operate within the safety constraints as defined for the existing equilibrium core of PARR-1.     

Test irradiations of full-sized U3Si2-Al fuel plates up to very high fission densities

In the course of the licensing procedure of the ‘Forschungsneutronenquelle Heinz Maier-Leibnitz’, i.e. the new 20 MW high-flux research reactor FRM II in Garching near Munich, extensive test irradiations have been performed to qualify the U₃Si₂-Al dispersion fuel with a relatively high density of highly enriched uranium (93 wt% of ²³⁵U) up to very high fission densities. Two of the three FRM II type fuel plates used in the irradiation tests contained U₃Si₂-Al dispersion fuel with HEU densities of 3.0 gU/cm³ or 1.5 gU/cm³ (‘homogeneous plates’) and one plate had two adjacent zones of either density (‘mixed plate’). They were irradiated in the French MTR reactors SILOE and OSIRIS in the years before 2002. The local plate thickness was measured on certain tracks along the plates during interruptions of the irradiation. The maximum fission density obtained in the U₃Si₂ fuel particles was 1.4 × 10²² f/cm³ and 1.1 × 10²² f/cm³ in the 1.5 gU/cm³ and 3.0 gU/cm³ fuel zones, respectively. In the course of the irradiations, the plate thickness increased monotonically and approximately linearly, leading to a maximum plate thickness swelling of 14% and 21% and a corresponding volume increase of the fuel particles of 106% and 81%, respectively. Our results are discussed and compared with the data from the literature.     

Reactivity feedback coefficients of a material test research reactor fueled with high-density U3Si2 dispersion fuels

The reactivity feedback coefficients of a material test research reactor fueled with high-density U₃Si₂ dispersion fuels were calculated. For this purpose, the low-density LEU fuel of an MTR was replaced with high-density U₃Si₂ LEU fuels currently being developed under the RERTR program. Calculations were carried out to find the fuel temperature reactivity coefficient, moderator temperature reactivity coefficient and moderator density reactivity coefficient. Nuclear reactor analysis codes including WIMS-D4 and CITATION were employed to carry out these calculations. It is observed that the average values of fuel temperature reactivity feedback coefficient, moderator temperature reactivity coefficient and moderator density reactivity coefficient from 20 °C to 100 °C, at the beginning of life, followed the relationships (in units of Δk/k × 10⁻⁵ K⁻¹) −2.116 − 0.118 ρ_U, 0.713 − 37.309/ρ_U and −12.765 − 34.309/ρ_U, respectively for 4.0 ≤ ρ_U (g/cm³) ≤ 6.0.     

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.     

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.     

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.     

Simulation of reactivity transients in current MTRs

This work aims at simulation of reactivity induced transients in High Enriched Uranium (HEU) and Low Enriched Uranium (LEU) cores of a typical Material Test research Reactor (MTR) using PARET code. The transient problem was forced through specification of externally inserted reactivity as a function of time. Reactivity insertions are idealized by ramps and steps. Superdelayed-critical transients, superprompt-critical transients and quasistatic transients are selected for the analysis. Ramp and step reactivity functions were employed to simulate these perturbations. The effect of initial power on transient behavior has also been investigated. The low enriched uranium core is analyzed for transients without scram. The magnitudes of maximum reactivity insertions are chosen to be in the range of $0.05 to 2.0 for different reactivity insertion times. Transient simulation with scram reveals that response of both HEU and LEU-cores is similar for selected ‘ramps’ and ‘steps’. The difference is observed in the peak values of power and coolant, clad and fuel temperatures. Trip level is achieved earlier in case of LEU-core. The peak clad temperatures in both LEU and HEU-cores remain below the melting point of aluminum-clad for the selected reactivity insertions. Simulation show that the LEU-core is more sensitive to perturbations at low power as compared to the transients at full power. For reactivity transients at low power level, power rises sharply to a higher peak value. In transients at full power, the peak power barely exceeds the trip level. The power oscillations after the first peak are observed for transients without scram.     

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.     

Use of Thorium-Plutonium MOX in the inner pins of CANDU fuel bundles

This research is focused on using Thorium-Plutonium MOX fuel in the inner fuel pins of the CANDU fuel bundles for plutonium incineration and reduction of uranium demand and to reduce coolant void reactivity. The delayed neutron fraction and the power distribution amongst the fuel elements of the fuel bundle have been considered as main safety parameters.The 700 MWe Advanced CANDU Reactor (ACR-700) was selected as a case study. The inner eight UO₂ fuel pins of the ACR-700 fuel bundle are replaced by Thorium-Plutonium MOX fuel pins in the proposed design with 3% reactor grade PuO₂. This amount represents 23.4 w/o of the fuel in the bundle. The outer two fuel rings (35 pins) enrichment is reduced from 2.1 w/o U-235 to 2 w/o U-235. The simulation using MCNP6 showed that about 27% reduction of uranium demand can be achieved. The proposed fuel bundle eliminate the use of burnable poisons in the central pin that was used for negative coolant void reactivity and more reduction in the coolant void reactivity was achieved (about 3.5 mk less than the reference fuel bundle). The power distribution throughout the fuel bundle is more flat in the proposed fuel bundle. Use of this fuel bundle reduces the delayed neutron fraction from 540 pcm in the reference case to 480 pcm in the proposed case.     

Neutronic analysis for core conversion (HEU–LEU) of Pakistan research reactor-2 (PARR-2)

Neutronic analyses for the core conversion of Pakistan research reactor-2 (PARR-2) from high enriched uranium (HEU) fuel to low enriched uranium (LEU) fuel has been performed. Neutronic model has been verified for 90.2% enriched HEU fuel (UAl₄–Al). For core conversion, UO₂ fuel was chosen as an appropriate fuel option because of higher uranium density. Clad has been changed from aluminum to zircalloy-4. Uranium enrichment of 12.6% has been optimized based on the design basis criterion of excess reactivity 4 mk in miniature neutron source reactor (MNSR). Lattice calculations for cross-section generation have been performed utilizing WIMS while core modeling was carried out employing three dimensions option of CITATION. Calculated neutronic parameters were compared for HEU and LEU fuels. Comparison shows that to get same thermal neutron flux at inner irradiation sites, reactor power has to be increased from 30 to 33 kW for LEU fuel. Reactivity coefficients calculations show that doppler and void coefficient values of LEU fuel are higher while moderator coefficient of HEU fuel is higher. It is concluded that from neutronic point of view LEU fuel UO₂ of 12.6% enrichment with zircalloy-4 clad is suitable to replace the existing HEU fuel provided that dimensions of fuel pin and total number of fuel pins are kept same as for HEU fuel.     

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

The effects of using high density low enriched uranium on the neutronic parameters 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 cross-sections affecting the reactor in different ways, therefore fuels U–Mo (9 w/o) which contain the same elements in same ratio were selected for analysis. Simulations were carried out to calculate core excess reactivity, neutron flux spectrum, prompt neutron generation time, effective delayed neutron fraction and feedback coefficients including Doppler feedback coefficient, and reactivity coefficients for change of water density and temperature. Nuclear reactor analysis codes including WIMS-D4 and CITATION were employed to carry out these calculations. It is observed that the excess reactivity at the beginning of life does not increase as the uranium density of fuel. Both the prompt neutron generation time and the effective delayed neutron fraction decrease as the uranium density increases. The absolute value of Doppler feedback coefficient increases while the absolute values of reactivity coefficients for change of water density and temperature decrease.     

国际RERTR计划发展现状

通过国际间十年的共同努力,降低研究试验堆铀浓缩度(RERTR)计划已经取得了很大的进展,找到了一种性能良好、铀密度高的低浓铀(LEU)燃料——U_3Si_2-Al 弥散型燃料。美国核管会已于1988年7月批准研究堆使用4.3g/cm~3铀密度的 U_3Si_2-Al 燃料。本文简要地介绍了几个主要国家在这方面的进展情况,以及这种燃料元件在制造、性能方面的特点。     

Determination of thermal hydraulic data of GHARR-1 under reactivity insertion transients using the PARET/ANL code

The PARET/ANL code has been adapted by the IAEA for testing transient behaviour in research reactors since it provides a coupled thermal hydrodynamic and point kinetics capability for estimating thermal hydraulic margins. A two-channel power peaking profile of the Ghana Research Reactor-1 (GHARR-1) has been developed for the PARET/ANL (Version 7.3; 2007) using the Monte Carlo N-Particle code (MCNP) to determine the thermal hydraulic data for reactivity insertion transients in the range of 2.0 × 10⁻³ Δk/k to 5.5 × 10⁻³ Δk/k. Peak clad and coolant temperatures ranged from 59.18 °C to 112.36 °C and 42.95 °C to 79.42 °C respectively. Calculated safety margins (DNBR) satisfied the MNSR thermal hydraulic design criteria for which no boiling occurs in the reactor core. The generated thermal hydraulic data demonstrated a high inherent safety feature of GHARR-1 for which the high negative reactivity feedback of the moderator limits power excursion and consequently the escalation of the clad temperature.