MONTE CARLO PERTURBATION ANALYSIS OF FUEL TEMPERATURE VARIATIONS IN THE MCNP MODEL OF THE ANNULAR CORE RESEARCH REACTOR.
Abstract not provided.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Abstract not provided.
Two of the major challenges facing the nuclear industry are the large amounts of nuclear waste produced by the current light water reactors and the risk of nuclear proliferation. Most proposed solutions for the waste problem involve fast breeder reactors, but that doesn't help with the nuclear proliferation problem. It is generally understood that a uranium fueled (U-238 + U-235) thermal reactor cannot be a breeder because there aren't enough neutrons. That analysis does not include the effects of fast fission of U-238. This is not surprising because the percentage of fast fissions in existing thermal reactor designs is low, around 3% (in Candu reactors) or less. Table I shows the neutrons involved in 1,000 fissions in a thermal reactor for a conversion ratio of 1.000 when the fuel is mostly U-238 and the fissile content is low enough so that most fast fissions are of U-238. The top section of the table is for U-235 in U-238, the bottom section is for Pu-239 in U-238. The bottom line in each section shows the number of neutrons available for all other losses (escape from the reactor, absorption in non fuel components, and absorption in xenon, samarium, and other fission products). If the percentage of fast fissions of U-238 is 10 or 15 % a thermal breeder may be feasible. This is largely independent of the percentage of fissile content of the fuel. If the reactor starting fissile content is mostly U-235, as the U-235 is replaced with Pu-239 the number of available neutrons increases even if the conversion ratio is only 1.000. The percentage of fast fissions in a reactor is a question of reactor geometry, as the minimum dimension of the fuel elements increases the average distance a fast neutron travels in the fuel before exiting increases and with it the probability of fast fission of U-238.
At the Paul Scherrer Institut, the CASMO/SIMULATE code system is used as reference scheme for core analyses of the Swiss reactors. Besides, Monte Carlo (MC) based neutron transport codes are in parallel being established both for out-of-core applications and full core analyses, in order to provide higher fidelity solutions to problems beyond the modelling capabilities of the standard core modelling methods. In this perspective, a necessity is to be able to conduct the MC calculations at any operating conditions, requiring thus proper initialisation of the thermal-hydraulic feedback as well as at any burnup, prompting the need for a detailed modelling of all material zones in the irradiated fuel assemblies. To tackle the later challenge, the development of the COMPLINK tool, aimed at initialising large-scale MC models of burned core configurations, using as basis upstream 2-D assembly/ 3-D core analyses results, was initiated and is presented in this paper. A main feature is the capability to extend in an automatic manner, a user given MC model in order to represent material data with a higher resolution, keeping track of the location of such materials also within nested repeated structures. This capability prevents the user from the burden of developing specific computational framework for a target application. Since the same functionalities are required to analyse out-of-core spent fuel managements systems, a first verification of COMPLINK is illustrated on the basis of a criticality safety evaluation (CSE) with MCNP6.1 of a Swiss PWR wet storage pool loaded with fresh and burned fuel assemblies. Starting from CASMO-5 depletion calculations at each axial layer of the depleted fuel assemblies and reconstructed based on realistic irradiation histories, COMPLINK is then applied to load the isotopic compositions into the material zones of the corresponding MCNP6.1 assembly models. On this basis, the CSE analyses are conducted with focus on investigating axial and radial modelling effects. All in all, the obtained results are found to be well in line with the impact on CSE expected at the considered burnup levels, providing thus confidence in the correct implementation of the COMPLINK functionalities and confirming its effectiveness as well as efficiency for the initialization of all material zones in large-scale MC models of complex and strongly heterogeneous burned configurations. (authors)
Seismic isolators and damping devices have been used for 30+ years to reduce seismic forces and deployed in 10,000+ structures, including bridges, buildings, and mission-critical infrastructure. Studies have shown that seismic isolation of nuclear power plants can reduce the seismic risk by several orders of magnitude and substantially reduce the overnight capital cost. Advanced reactor developers are considering seismic isolation as an integral design feature in their reactor design. It can either be deployed at the building level, isolating the entire reactor building, or at the component level, isolating individual pieces of safety-class equipment, such as a reactor vessel or a steam generator. For equipment isolation, located close to the reactor vessel, the seismic isolators and/or damping devices may be exposed to gamma and neutron radiation, with a possible effect on their mechanical properties which would affect the response of the isolated equipment under earthquake shaking. The Department of Energy, as part of its Nuclear Energy University Program (NEUP), in collaboration with the Idaho National Laboratory (INL), is funding a project at the University at Buffalo (UB) to investigate the effect of gamma radiation on the mechanical properties of various seismic protective devices. (When isolating equipment, the seismic protective devices will be located outside the reactor vessel where the neutron exposure will be insignificant.
An MCNP model of a co-axial HPGe detector, defined by the specification sheet and without parameter optimization, was used to assess practicality for efficiency characterization and complex spectra simulation. Simulated efficiencies, with and without Gaussian energy broadening, were compared against an experimental calibration for a standard counting geometry. Results show that this minimal detector model predicted efficiencies within 6% of experimental values at energies above 100 keV, while including energy broadening offered little improvement. Furthermore, this work shows that extensive detector characterization and optimization of parameters is not absolutely necessary, thereby reducing the experimental burden associated with detector modeling for many applications.
The IBR-2M is a fast research reactor that operates in supercritical condition for ∼ 800 μs every 200 ms. Two reflector parts in nickel rotate in opposite directions generating 1.8 GWth peak power when they align with the fuel zone changing the reactor status from deep subcritical to supercritical. The reactor core uses high-enriched plutonium fuel and is cooled by sodium. This reactor has been modeled by MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry, whereas the SERPENT model employs the Stereolithographic (STL) geometry representation that can be used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The CUBIT program was also used for a three-dimensional visualization of the Monte Carlo models. SERPENT and MCNP models use the same geometry, material specifications, and nuclear data. The latter are based on the ENDF/B-7.0 library. SERPENT and MCNP using same geometry and same material specifications produce similar k{sub eff} values within 120 pcm.
IBR-2M is a fast research reactor located in Dubna (Russia), which operates in supercritical condition for ~800 μs every 200 ms. Two reflector parts are rotating in opposite directions generating 1.8 GWth peak power when the two parts align with the fuel zone. The reactor core uses plutonium fuel and sodium coolant. The control rods are located in stationary reflector zones. This report describes the computational modeling of the IBR-2M reactor. Two independent computational models were generated using MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry based on macro bodies, e.g. boxes, cylinders, hexagons, spheres, cones, etc., and surfaces. The SERPENT model is based on the STereo Lithographic (STL) geometry representation that is used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The water cooling system of the stationary and the rotating reflectors has a very complicated geometry, involving serpentines and irregular surfaces profiles. The materials in these areas have been homogenized using the GIMP image editor according to the material distributions obtained from the colors of the design drawings. The MCNP and the SERPENT inputs consist of ~1500 and ~1750 lines, respectively. The calculated effective multiplication factors from these two Monte Carlo computer programs with different models agree within 50 pcm.
Purpose: To develop a Monte Carlo N-Particle (MCNP) model for the validation of a prototype backscatter x-ray (BSX) imager, and optimization of BSX technology for medical applications, including selective object-plane imaging. Methods: BSX is an emerging technology that represents an alternative to conventional computed tomography (CT) and projective digital radiography (DR). It employs detectors located on the same side as the incident x-ray source, making use of backscatter and avoiding ring geometry to enclose the imaging object. Current BSX imagers suffer from low spatial resolution. A MCNP model was designed to replicate a BSX prototype used for flaw detection in industrial materials. This prototype consisted of a 1.5mm diameter 60kVp pencil beam surrounded by a ring of four 5.0cm diameter NaI scintillation detectors. The imaging phantom consisted of a 2.9cm thick aluminum plate with five 0.6cm diameter holes drilled halfway. The experimental image was created using a raster scanning motion (in 1.5mm increments). Results: A qualitative comparison between the physical and simulated images showed very good agreement with 1.5mm spatial resolution in plane perpendicular to incident x-ray beam. The MCNP model developed the concept of radiography by selective plane detection (RSPD) for BSX, whereby specific object planes can be imaged by varying kVp. 10keV increments in mean x-ray energy yielded 4mm thick slice resolution in the phantom. Image resolution in the MCNP model can be further increased by increasing the number of detectors, and decreasing raster step size. Conclusion: MCNP modelling was used to validate a prototype BSX imager and introduce the RSPD concept, allowing for selective object-plane imaging. There was very good visual agreement between the experimental and MCNP imaging. Beyond optimizing system parameters for the existing prototype, new geometries can be investigated for volumetric image acquisition in medical applications. This material is based upon work supported under an Integrated University Program Graduate Fellowship sponsored by the Department of Energy Office of Nuclear Energy.
Modeling and simulation is a cost- and time-saving method to study conceptual designs of radiation detection systems. MCNP radiation transport code has been an indispensable tool in development of radiation detection systems, and the PINS systems have also utilized MCNP modeling and simulation in the system development. MCNP models of the PINS systems have been used to understand detector responses and to fill the gap in experimental data. Recently, the MIND backpack system was announced to be developed, and its prototype design was required to conduct laboratory tests and demonstrations. A shadow shield is the component to be designed by Idaho National Laboratory in order to maximize the performance of a backpack system while meeting the sponsor’s requirements. MCNP modeling and simulation was utilized to find the best specification of a shadow shield. This study was performed to down-select a shadow shield design for the MIND system. A total of 24 designs for a shadow shield were identified to build MCNP models, and their simulation results were used to prioritize each design. As a quantitative performance metric, an equation of Figure-of-Merit was introduced in the decision making analysis. In addition, the Analytic Hierarchy Process was adopted to prioritize and recommend best performing shadow shield designs.
The University of Wisconsin Nuclear Reactor converted from HEU to LEU fuel in 2009. As part of the conversion process, an MCNP model was developed to predict flux distributions and depletion in the core. This core model included an unforeseen excess reactivity bias, resulting in an overestimate of core reactivity, and requiring the core configuration to be shuffled in 2011. The original model was constructed by hand and in order to account for radial and axial burnup in the fuel, and to include sufficient detail in the geometric features, the input file was over 17,000 lines long. Given the challenges in assessing the quality of such a long and complicated MCNP model, this project was designed to build an MCNP model based on a minimal set of input quantities in a modular and repeatable fashion. A Python code was implemented to write the MCNP input file from fundamental reactor data, both geometric and material. The recent conversion to LEU fuel provides a unique opportunity to characterize fuel composition quite accurately in the core. The goal of this project is to produce an accurate, maintainable model that can be used with future reactor physics experiments.
Neutronic analyses of the Oregon State TRIGA Reactor (OSTR) have historically been performed using a thoroughly validated MCNP model. This model was originally created in 1997 in support of boron neutron capture therapy (BNCT) research, when the core still utilized highly-enriched uranium (HEU) fuel. The OSTR core was converted to low-enriched uranium (LEU) in the fall of 2008. To more accurately represent the converted state of the core, the MCNP model was modified to incorporate fuel material information from the fuel manufacturer (AREVA). In previous versions of the model, the fuel elements were assumed to have a uniform isotopic composition. Chemical analysis from AREVA provided far more resolved material information, allowing for increased fidelity in the representation of the fuel. The model was also updated to ENDF/B-VII.1 cross section libraries as the previous model utilized ENDF/B-VI. The primary impetus for the update of the MCNP model was the desire to have a more accurate assessment of criticality conditions of the OSTR. The OSTR utilizes three core configurations which are defined by the contents of the B1 grid position. If fuel resides in B1, this is designated as the 'normal' core. Otherwise, there are two aluminum in-core irradiation tubes (ICITs) that can be placed in B1. There is an unlined tube, designated simply as the ICIT, and there is a cadmium-lined tube, designated as the CLICIT. Table 1 contains a summary of the results of a series of MCNP criticality calculations which used the original 'fresh fuel' isotopics from the fall of 2008 in conjunction with critical rod height data from December 2013. Calculations were performed at 15 W, which is where the daily core excess measurements are taken, and at 1 MW, which is the current full-power level at the OSTR. The k-effective values indicate that the fresh fuel material cards are highly inaccurate and over-predict the criticality by at least $1.20 in every core configuration (based on a β{sub eff} value of 0.0075), thus it was desired to perform a depletion calculation to obtain more accurate fuel information. (authors)
The Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL) is one of the premier irradiation facilities in the world. The ATR hosts experiments from a variety of laboratories, companies, and disciplines. Because of its flexibility in providing customized irradiation conditions to experiments, such as coolant chemistry, temperature, and reactor power, analysis for experiments, both from a programmatic and safety perspective, can be a challenging endeavor. Because of these experiment variables, there is always an effort to improve upon the existing analysis methodology in order to provide the sponsor of any given experiment with the most accurate predictions possible. The paper discusses the effort to improve upon the methodology for loading the ATR driver core for experiment analysis. Experiment analysis at INL is generally broken up into three separate disciplines: neutronics, thermal-hydraulics, and structural. The focus of this paper is neutronics analysis. After the experiment has a preliminary design, the neutronics analyst will create a model in the Monte Carlo N-Particle (MCNP) transport code to calculate fuel and material heating rates from MCNP tallies and provide those to the thermal-hydraulics analyst, and calculate reactivity worths, perturbation on the axial fission profile of the ATR driver core by the experiment, void worths, and experiment backup worths. Also, there are usually several programmatic parameters that the experiment sponsor requests that the neutronics analyst calculates. These can include heat generation rates, burnup and displacements per atom (DPA). In order to accurately analyze these parameters, an important input into the physics analyst's MCNP model is the loading of the ATR driver core. The ATR driver core consists of 40, 19-plate highly-enriched uranium (HEU) fuel elements, arranged in a serpentine shape, as pictured in Figure 1. The ATR driver core is responsible for the vast majority of the thermal power of the reactor and the nuclear interactions with the experiments being irradiated in ATR. A unique feature of ATR is the 'power tilt' across the reactor. A power tilt means that the Northwest, Northeast, Southwest, and Southeast Lobes of ATR frequently operate at different powers, providing different irradiation conditions for the experiments in and around the lobes of the reactor. In order to achieve lobe power splits, different fuel elements for each of the 40 positions are chosen by ATR Reactor Engineering so that lobe power splits can be maintained during an operating cycle. A typical fuel loading for an ATR cycle consists of 14 to 18 new fuel elements and 2 to 26 depleted fuel elements. These loadings are typically chosen from three different types of fuel elements. All three ATR fuel element types share a common design, spelled out in Chapter 4 of the ATR Safety Analysis Report (SAR). The fuel element meat is HEU UAlx dispersion fuel sandwiched between aluminum cladding with a length of 121.92 cm (4 ft), with 19 plates per fuel element. The three variations of this common design are designated by the first two letters in their serial numbers: XA, YA, and NB. XA elements are most commonly used. These elements have nominally 1075 g of uranium-235 and 0.66 g of boron-10. The NB elements have the same uranium-235 loading as the XA elements, but no boron in any plates. These elements are used when excess reactivity is needed. The YA elements have the same loading as the XA elements, except that the outermost plate (Plate 19) has no boron or uranium. These elements are used to mitigate aged reflector damage. In order to accurately model the ATR driver core, all three types of fuel elements must be accounted for in the MCNP model. (authors)
This work was performed to support the feasibility study on the potential conversion of the Transient Reactor Test Facility (TREAT) at Idaho National Laboratory from the use of high-enriched uranium (HEU) fuel to the use of low-enriched uranium (LEU) fuel. The analyses were performed by the GTRI Reactor Conversion staff at the Argonne National Laboratory (ANL). The objective of this study was to validate the MCNP model of the TREAT reactor with the well-documented measurements which were taken during the start-up and early operation of TREAT. Furthermore, the effect of carbon graphitization was also addressed. The graphitization level was assumed to be 100% (ANL/GTRI/TM-13/4). For this purpose, a set of experiments was chosen to validate the TREAT MCNP model, involving the approach to criticality procedure, in-core neutron flux measurements with foils, and isothermal temperature coefficient and temperature distribution measurements. The results of this study extended the knowledge base for the TREAT MCNP calculations and established the credibility of the MCNP model to be used in the core conversion feasibility analysis.
Radiation damage estimation is an important component of the post irradiation analysis of the Advanced Graphite Creep (AGC) experiment. It depends primarily on the fast fluence, which is determined using well established methods of spectral adjustment. These are based on best estimates from models such as Monte Carlo N-Particle (MCNP), input cross-sections, and measured activities from flux wires in the experiment. Each of these parameters can propagate uncertainties which will affect the uncertainty in the calculated dose levels for AGC, or any experiment irradiated within a reactor. While the methods of propagating uncertainty are well-established, the final uncertainty estimates they provide are only as good as the estimates of uncertainty in the inputs on which they are based. The purpose of this work is to outline some deficiencies in the ways these input uncertainties are presently estimated, and to outline a methodology by which they can be improved. The fast fluence and radiation damage received by graphite specimens irradiated in the Advanced Graphite Creep (AGC) experiments is presently estimated using spectral adjustment methods that are based on both flux wire activity measurements, and MCNP model predictions. This work describes an ongoing effort to quantify and propagate uncertainties in inputs to the spectral adjustment process, and thereby quantify the resultant error in radiation damage (dpa) estimates. The effort is multi-faceted, and we consider the impacts of both the set of flux wires selected, and the counting process. An expanded set of flux wires is identified that provides a more comprehensive data set on the fast spectrum. To address the counting process itself, a series of round-robin measurements in several reactor metrology laboratories across the Department of Energy (DOE) complex and nuclear industry are being undertaken to refine the American Society for Testing and Materials (ASTM) standards for flux wire measurements. To address the contribution of uncertainty in the MCNP model predictions, an uncertainty quantification (UQ) tool has been developed that statistically samples the model input parameters, runs a series of cases, and assimilates the results to provide an overall uncertainty. The impact of the MCNP UQ tool results is demonstrated by re-analyzing previous AGC flux wire and irradiation data. While the expanded flux wire set obviously cannot be added to these experiments retroactively, plans for future graphite irradiations are outlined.