A preliminary study of target multiplicity for ADSRs

The Accelerator Driven Subcritical Reactor (ADSR) concept is based on the coupling of a particle accelerator to a subcritical reactor core by means of a neutron spallation target interface. This paper investigates the benefits of multiple spallation targets in ADSRs. The motivation behind this is, firstly, to improve the overall reliability of the accelerator-reactor system, and, secondly, to evaluate other potential advantages such as lower beam power requirements. The results show that a system containing two or three spallation targets, coupled to independent accelerators, offers better neutronic performance. This is demonstrated through the increased effective multiplication factor (keff) in the two- and three-target configurations and a more uniform neutron flux distribution. A multiple-target ADSR also proves effective in mitigating the impact of frequent beam interruptions, a pressing issue that needs to be addressed for the ADSR concept to advance. Assuming no simultaneous beam shutdowns, the two- and three-target configurations reduce the risk of fuel cladding failure due to thermal cyclic fatigue. © 2013 Elsevier B.V. All rights reserved.

Complex chemistry DNS of n-heptane spray autoignition at high pressure and intermediate temperature conditions

Direct Numerical Simulations (DNS) of turbulent n-heptane sprays autoigniting at high pressure (P=24bar) and intermediate air temperature (Tair=1000K) have been performed to investigate the physical mechanisms present under conditions where low-temperature chemistry is expected to be important. The initial turbulence in the carrier gas, the global equivalence ratio in the spray region, and the initial droplet size distribution of the spray were varied. Results show that spray ignition exhibits a spotty nature, with several kernels developing independently in those regions where the mixture fraction is close to its most reactive value ξMR (as determined from homogeneous reactor calculations) and the scalar dissipation rate is low. Turbulence reduces the ignition delay time as it promotes mixing between air and the fuel vapor, eventually resulting in lower values of scalar dissipation. High values of the global equivalence ratio are responsible for a larger number of ignition kernels, due to the higher probability of finding regions where ξ=ξMR. Spray polydispersity results in the occurrence of ignition over a wider range of mixture fraction values. This is a consequence of the inhomogeneities in the mixing field that characterize these sprays, where poorly mixed rich spots are seen to alternate with leaner ones which are well-mixed. The DNS simulations presented in this work have also been used to assess the applicability of the Conditional Moment Closure (CMC) method to the simulation of spray combustion. CMC is found to be a valid method for capturing spray autoignition, although care should be taken in the modelling of the unclosed terms appearing in the CMC equations. © 2013 The Combustion Institute.

A fast detailed-chemistry modelling approach for simulating the SI-HCCI transition

An established Stochastic Reactor Model (SRM) is used to simulate the transition from Spark Ignition (SI) to Homogeneous Charge Compression Ignition (HCCI) combustion mode in a four cylinder in-line four-stroke naturally aspirated direct injection SI engine with cam profile switching. The SRM is coupled with GT-Power, a one-dimensional engine simulation tool used for modelling engine breathing during the open valve portion of the engine cycle, enabling multi-cycle simulations. The mode change is achieved by switching the cam profiles and phasing, resulting in a Negative Valve Overlap (NVO), opening the throttle, advancing the spark timing and reducing the fuel mass as well as using a pilot injection. A proven technique for tabulating the model is used to create look-up tables in both SI and HCCI modes. In HCCI mode several tables are required, including tables for the first NVO, transient valve timing NVO, transient valve timing HCCI and steady valve timing HCCI and NVO. This results in the ability to simulate the transition with detailed chemistry in very short computation times. The tables are then used to optimise the transition with the goal of reducing NO x emissions and fluctuations in IMEP. Copyright © 2010 SAE International.

A detailed chemistry simulation of the SI-HCCI transition

A Stochastic Reactor Model (SRM) has been used to simulate the transition from Spark Ignition (SI) mode to Homogeneous Charge Compression Ignition (HCCI) mode in a four cylinder in-line four-stroke naturally aspirated direct injection SI engine with cam profile switching. The SRM is coupled with GT-Power, a one-dimensional engine simulation tool used for modelling engine breathing during the open valve portion of the engine cycle, enabling multi-cycle simulations. The model is initially calibrated in both modes using steady state data from SI and HCCI operation. The mode change is achieved by switching the cam profiles and phasing, resulting in a Negative Valve Overlap (NVO), opening the throttle, advancing the spark timing and reducing the fuel mass as well as utilising a pilot injection. Experimental data is presented along with the simulation results. The model is used to investigate key control parameters and their effects on parameters that are difficult to measure experimentally. The effect of the spark in the first HCCI cycles is found to have a major impact on the stability of the transition. Copyright © 2010 SAE International.

Modeling of SFR cores with Serpent-DYN3D codes sequence

DYN3D reactor dynamics nodal diffusion code was originally developed for the analysis of Light Water Reactors. In this paper, we demonstrate the feasibility of using DYN3D for modeling of fast spectrum reactors. A homogenized cross sections data library was generated using continuous energy Monte-Carlo code Serpent which provides significant modeling flexibility compared with traditional deterministic lattice transport codes and tolerable execution time. A representative sodium cooled fast reactor core was modeled with the Serpent-DYN3D code sequence and the results were compared with those produced by ERANOS code and with a 3D full core Monte-Carlo solution. Very good agreement between the codes was observed for the core integral parameters and power distribution suggesting that the DYN3D code with cross section library generated using Serpent can be reliably used for the analysis of fast reactors. © 2012 Elsevier Ltd. All rights reserved.

Neutronic optimization in high conversion Th-233U fuel assembly with simulated annealing

This paper reports on fuel design optimization of a PWR operating in a self sustainable Th-233U fuel cycle. Monte Carlo simulated annealing method was used in order to identify the fuel assembly configuration with the most attractive breeding performance. In previous studies, it was shown that breeding may be achieved by employing heterogeneous Seed-Blanket fuel geometry. The arrangement of seed and blanket pins within the assemblies may be determined by varying the designed parameters based on basic reactor physics phenomena which affect breeding. However, the amount of free parameters may still prove to be prohibitively large in order to systematically explore the design space for optimal solution. Therefore, the Monte Carlo annealing algorithm for neutronic optimization is applied in order to identify the most favorable design. The objective of simulated annealing optimization is to find a set of design parameters, which maximizes some given performance function (such as relative period of net breeding) under specified constraints (such as fuel cycle length). The first objective of the study was to demonstrate that the simulated annealing optimization algorithm will lead to the same fuel pins arrangement as was obtained in the previous studies which used only basic physics phenomena as guidance for optimization. In the second part of this work, the simulated annealing method was used to optimize fuel pins arrangement in much larger fuel assembly, where the basic physics intuition does not yield clearly optimal configuration. The simulated annealing method was found to be very efficient in selecting the optimal design in both cases. In the future, this method will be used for optimization of fuel assembly design with larger number of free parameters in order to determine the most favorable trade-off between the breeding performance and core average power density.

Comparison of square and hexagonal fuel lattices for high conversion PWRs

This paper reports on an investigation into fuel design choices of a pressurized water reactor operating in a self-sustainable Th- 233U fuel cycle. In order to evaluate feasibility of this concept, two types of fuel assembly lattices were considered: square and hexagonal. The hexagonal lattice may offer some advantages over the square one. For example, the fertile blanket fuel can be packed more tightly reducing the blanket volume fraction in the core and potentially allowing to achieve higher core average power density. The calculations were carried out with Monte-Carlo based BGCore code system and the results were compared to those obtained with Serpent Monte-Carlo code and deterministic transport code BOXER. One of the major design challenges associated with the SB concept is high power peaking due to the high concentration of fissile material in the seed region. The second objective of this work is to estimate the maximum achievable core power density by evaluation of limiting thermal hydraulic parameters. The analysis showed that both fuel assembly designs have a potential of achieving net breeding. Although hexagonal lattice was found to be somewhat more favorable because it allows achieving higher power density, while having breeding performance comparable to the square lattice case. © Carl Hanser Verlag München.

Coupled neutronic thermo-hydraulic analysis of full PWR core with Monte-Carlo based BGCore system

BGCore reactor analysis system was recently developed at Ben-Gurion University for calculating in-core fuel composition and spent fuel emissions following discharge. It couples the Monte Carlo transport code MCNP with an independently developed burnup and decay module SARAF. Most of the existing MCNP based depletion codes (e.g. MOCUP, Monteburns, MCODE) tally directly the one-group fluxes and reaction rates in order to prepare one-group cross sections necessary for the fuel depletion analysis. BGCore, on the other hand, uses a multi-group (MG) approach for generation of one group cross-sections. This coupling approach significantly reduces the code execution time without compromising the accuracy of the results. Substantial reduction in the BGCore code execution time allows consideration of problems with much higher degree of complexity, such as introduction of thermal hydraulic (TH) feedback into the calculation scheme. Recently, a simplified TH feedback module, THERMO, was developed and integrated into the BGCore system. To demonstrate the capabilities of the upgraded BGCore system, a coupled neutronic TH analysis of a full PWR core was performed. The BGCore results were compared with those of the state of the art 3D deterministic nodal diffusion code DYN3D (Grundmann et al.; 2000). Very good agreement in major core operational parameters including k-eff eigenvalue, axial and radial power profiles, and temperature distributions between the BGCore and DYN3D results was observed. This agreement confirms the consistency of the implementation of the TH feedback module. Although the upgraded BGCore system is capable of performing both, depletion and TH analyses, the calculations in this study were performed for the beginning of cycle state with pre-generated fuel compositions. © 2011 Published by Elsevier B.V.

Detailed design of242mAm breeding in pressurized water reactors

There is growing interest in the use of 242mAm as a nuclear fuel. Because of its very high thermal fission cross section and its large number of neutrons released per fission, it can be used for various unique applications, such as space propulsion, medical applications, and compact energy sources. Since the thermal absorption cross section of 242mAm is very high, the best way to obtain 242mAm is by the capture of fast or epithermal neutrons in 241Am. However, fast spectrum reactors are not readily available. In this paper, we explore the possibility of producing 242mAm in existing pressurized water reactors (PWRs) with minimal interference in reactor performance. As suggested in previous studies on the subject, the 242mAm breeding targets are shielded with strong thermal absorbers in order to suppress the thermal neutron flux that causes 242mAm destruction. Since 242mAm enrichment within the Am target mainly depends on the neutron energy distribution, which in turn depends on the Am target thickness and on the neutron filter cutoff energy (thermal absorber type), this unique Am target design was developed. In our study, Cd, Sm, and Gd were considered as thermal neutron filters, as suggested by Cesana et al. The most favorable results were obtained by irradiating Am targets covered either with Gd or Cd. In these cases, up to 8.65% enrichment of 242mAm is obtained after 4.5 yr (three successive PWR fuel cycles) of irradiation. It was also found that significant quantities [up to 1.3 kg/GW (electric)-yr] of 242mAm can be obtained in PWR reactors without notable interference with reactor performance. However, in order to maintain the original fuel cycle length, the enrichment of the driver (UO2) fuel must be increased by ∼1%, raised from the conventional 4.5 to 5.5%, depending on the thermal neutron filter used. The most important reactivity feedback coefficients for fuel assemblies containing the 242mAm breeding targets were evaluated and found to be close to those of a standard PWR. Another product of neutron capture in the 241Am reaction is 238Pu. It was found that in a typical 1000 MW (electric) PWR core with one-third of the fuel assemblies containing 241Am targets, up to 15.1 kg of 238Pu enriched to 80% can be produced per year.

HTGR fuel element depletion benchmark: Stage three results

Recently, a new numerical benchmark exercise for High Temperature Gas Cooled Reactor (HTGR) fuel depletion was defined. The purpose of this benchmark is to provide a comparison basis for different codes and methods applied to the burnup analysis of HTGRs. The benchmark specifications include three different models: (1) an infinite lattice of tristructural isotropic (TRISO) fuel particles, (2) an infinite lattice of fuel pebbles, and (3) a prismatic fuel including fuel and coolant channels. In this paper, we present the results of the third stage of the benchmark obtained with MCNP based depletion code BGCore and deterministic lattice code HELIOS 1.9. The depletion calculations were performed for three-dimensional model of prismatic fuel with explicitly described TRISO particles as well as for two-dimensional model, in which double heterogeneity of the TRISO particles was eliminated using reactivity equivalent physical transformation (RPT). Generally, good agreement in the results of the calculations obtained using different methods and codes was observed.

The potential use of 241Am as proliferation resistant burnable poison in PWRs

This paper discusses the use of 241Am as proliferation resistant burnable poison for light water reactors. Homogeneous addition of small (as little as 0.12%) amounts of 241Am to the conventional light water reactor fuel results in significant increase in 238Pu/Pu ratio in the discharged fuel improving its proliferation resistance. Moreover, 241Am, admixed to the fuel, acts as burnable absorber allowing for substantial reduction in conventional reactivity control means without a notable fuel cycle length penalty. This is possible due to favorable characteristics of 241Am transmutation chain. The fuel cycle length penalty of introducing 241Am into the core is evaluated and discussed, as well as the impact of He production in the fuel pins and degradation of reactivity feedback coefficients. Proliferation resistance and reactivity control features related to the use of 241Am are compared to those of using 237Np, which has also been suggested as an additive to the conventional fuel in order to improve its proliferation resistance. It was found that 241Am admixture is more favorable than 237Np admixture because of the smaller fuel cycle length penalty and higher burnable poison savings. Addition of either 237Np or 241Am would provide substantial but not ultimate protection from misuse of Pu originating in the spent fuel from the commercial power reactors. Therefore, the benefits from application of the concept would have to be carefully evaluated against the additional costs and proliferation risks associated with manufacturing of 237Np or 241Am doped fuel. Although this work concerns specifically with PWRs, the conclusions could also be applied to BWRs and, to some extent, to other thermal spectrum reactor types. © 2009 Elsevier Ltd. All rights reserved.

Decay heat in fast reactors with transuranic fuels

In this work, we performed an evaluation of decay heat power of advanced, fast spectrum, lead and molten salt-cooled reactors, with flexible conversion ratio. The decay heat power was calculated using the BGCore computer code, which explicitly tracks over 1700 isotopes in the fuel throughout its burnup and subsequent decay. In the first stage, the capability of the BGCore code to accurately predict the decay heat power was verified by performing a benchmark calculation for a typical UO2 fuel in a Pressurized Water Reactor environment against the (ANSI/ANS-5.1-2005, "Decay Heat Power in Light Water Reactors," American National Standard) standard. Very good agreement (within 5%) between the two methods was obtained. Once BGCore calculation capabilities were verified, we calculated decay power for fast reactors with different coolants and conversion ratios, for which no standard procedure is currently available. Notable differences were observed for the decay power of the advanced reactor as compared with the conventional UO2 LWR. The importance of the observed differences was demonstrated by performing a simulation of a Station Blackout transient with the RELAP5 computer code for a lead-cooled fast reactor. The simulation was performed twice: using the code-default ANS-79 decay heat curve and using the curve calculated specifically for the studied core by BGCore code. The differences in the decay heat power resulted in failure to meet maximum cladding temperature limit criteria by ∼100 °C in the latter case, while in the transient simulation with the ANS-79 decay heat curve, all safety limits were satisfied. The results of this study show that the design of new reactor safety systems must be based on decay power curves specific to each individual case in order to assure the desired performance of these systems. © 2009 Elsevier B.V. All rights reserved.

Flexible conversion ratio fast reactors: Overview

Conceptual designs of lead-cooled and liquid salt-cooled fast flexible conversion ratio reactors were developed. The performance achievable by the unity conversion ratio cores of these reactors was compared to an existing supercritical carbon dioxide-cooled (S-CO2) fast reactor design and an uprated version of an existing sodium-cooled fast reactor. All concepts have cores rated at 2400 MWt. The cores of the liquid-cooled reactors are placed in a large-pool-type vessel with dual-free level, which also contains four intermediate heat exchangers (IHXs) coupling a primary coolant to a compact and efficient supercritical CO2 Brayton cycle power conversion system. The S-CO2 reactor is directly coupled to the S-CO2 Brayton cycle power conversion system. Decay heat is removed passively using an enhanced reactor vessel auxiliary cooling system (RVACS) and a passive secondary auxiliary cooling system (PSACS). The selection of the water-cooled versus air-cooled heat sink for the PSACS as well as the analysis of the probability that the PSACS may fail to complete its mission was performed using risk-informed methodology. In addition to these features, all reactors were designed to be self-controllable. Further, the liquid-cooled reactors utilized common passive decay heat removal systems whereas the S-CO2 uses reliable battery powered blowers for post-LOCA decay heat removal to provide flow in well defined regimes and to accommodate inadvertent bypass flows. The multiple design limits and challenges which constrained the execution of the four fast reactor concepts are elaborated. These include principally neutronics and materials challenges. The neutronic challenges are the large positive coolant reactivity feedback, small fuel temperature coefficient, small effective delayed neutron fraction, large reactivity swing and the transition between different conversion ratio cores. The burnup, temperature and fluence constraints on fuels, cladding and vessel materials are elaborated for three categories of material - materials currently available, available on a relatively short time scale and available only with significant development effort. The selected fuels are the metallic U-TRU-Zr (10% Zr) for unity conversion ratio and TRU-Zr (75% Zr) for zero conversion ratio. The principal selected cladding and vessel materials are HT-9 and A533 or A508, respectively, for current availability, T-91 and 9Cr-1Mo steel for relatively short-term availability and oxide dispersion strengthened ferritic steel (ODS) available only with significant development. © 2009 Elsevier B.V. All rights reserved.

Implementation of multi-group cross-section methodology in BGCore MC-depletion code

BGCore is a software package for comprehensive computer simulation of nuclear reactor systems and their fuel cycles. The BGCore interfaces Monte Carlo particles transport code MCNP4C with a SARAF module - an independently developed code for calculating in-core fuel composition and spent fuel emissions following discharge. In BGCore system, depletion coupling methodology is based on the multi-group approach that significantly reduces computation time and allows tracking of large number of nuclides during calculations. In this study, burnup calculation capabilities of BGCore system were validated against well established and verified, computer codes for thermal and fast spectrum lattices. Very good agreement in k eigenvalue and nuclide densities prediction was observed for all cases under consideration. In addition, decay heat prediction capabilities of the BGCore system were benchmarked against the most recent edition of ANS Standard methodology for UO2 fuel decay power prediction in LWRs. It was found that the difference between ANS standard data and that predicted by the BGCore does not exceed 5%.

Efficient generation of one-group cross sections for coupled Monte Carlo depletion calculations

Coupled Monte Carlo depletion systems provide a versatile and an accurate tool for analyzing advanced thermal and fast reactor designs for a variety of fuel compositions and geometries. The main drawback of Monte Carlo-based systems is a long calculation time imposing significant restrictions on the complexity and amount of design-oriented calculations. This paper presents an alternative approach to interfacing the Monte Carlo and depletion modules aimed at addressing this problem. The main idea is to calculate the one-group cross sections for all relevant isotopes required by the depletion module in a separate module external to Monte Carlo calculations. Thus, the Monte Carlo module will produce the criticality and neutron spectrum only, without tallying of the individual isotope reaction rates. The onegroup cross section for all isotopes will be generated in a separate module by collapsing a universal multigroup (MG) cross-section library using the Monte Carlo calculated flux. Here, the term "universal" means that a single MG cross-section set will be applicable for all reactor systems and is independent of reactor characteristics such as a neutron spectrum; fuel composition; and fuel cell, assembly, and core geometries. This approach was originally proposed by Haeck et al. and implemented in the ALEPH code. Implementation of the proposed approach to Monte Carlo burnup interfacing was carried out through the BGCORE system. One-group cross sections generated by the BGCORE system were compared with those tallied directly by the MCNP code. Analysis of this comparison was carried out and led to the conclusion that in order to achieve the accuracy required for a reliable core and fuel cycle analysis, accounting for the background cross section (σ0) in the unresolved resonance energy region is essential. An extension of the one-group cross-section generation model was implemented and tested by tabulating and interpolating by a simplified σ0 model. A significant improvement of the one-group cross-section accuracy was demonstrated.