Review of Axial Burnup Distribution Considerations for Burnup Credit Calculations
Wagner, J.C.; DeHart, M.D.
2000-03-01
This report attempts to summarize and consolidate the existing knowledge on axial burnup distribution issues that are important to burnup credit criticality safety calculations. Recently released Nuclear Regulatory Commission (NRC) staff guidance permits limited burnup credit, and thus, has prompted resolution of the axial burnup distribution issue. The reactivity difference between the neutron multiplication factor (keff) calculated with explicit representation of the axial burnup distribution and keff calculated assuming a uniform axial burnup is referred to as the ``end effect.'' This end effect is shown to be dependent on many factors, including the axial-burnup profile, total accumulated burnup, cooling time, initial enrichment, assembly design, and the isotopics considered (i.e., actinide-only or actinides plus fission products). Axial modeling studies, efforts related to the development of axial-profile databases, and the determination of bounding axial profiles are also discussed. Finally, areas that could benefit from further efforts are identified.
Detailed Burnup Calculations for Testing Nuclear Data
NASA Astrophysics Data System (ADS)
Leszczynski, F.
2005-05-01
-section data for burnup calculations, using some of the main available evaluated nuclear data files (ENDF-B-VI-Rel.8, JEFF-3.0, JENDL-3.3), on an isotope-by-isotope basis as much as possible. The selected experimental burnup benchmarks are reference cases for LWR and HWR reactors, with analysis of isotopic composition as a function of burnup. For LWR (H2O-moderated uranium oxide lattices) four benchmarks are included: ATM-104 NEA Burnup credit criticality benchmark; Yankee-Rowe Core V; H.B.Robinson Unit 2 and Turkey Point Unit 3. For HWR (D2O-moderated uranium oxide cluster lattices), three benchmarks were selected: NPD-19-rod Fuel Clusters; Pickering-28-rod Fuel Clusters; and Bruce-37-rod Fuel Clusters. The isotopes with experimental concentration data included in these benchmarks are: Se-79, Sr90, Tc99, Ru106, Sn126, Sb125,1129, Cs133-137, Nd143, 145, Sm149-150, 152, Eul53-155, U234-235, 238, Np237, Pu238-242, Am241-243, and Cm242-248. Results and analysis of differences between calculated and measured absolute and/or relative concentrations of these isotopes for the seven benchmarks are included in this work.
Gauntt, Randall O.; Powers, Dana Auburn; Ashbaugh, Scott G.; Leonard, Mark Thomas; Longmire, Pamela
2010-04-01
In this study, risk-significant pressurized-water reactor severe accident sequences are examined using MELCOR 1.8.5 to explore the range of fission product releases to the reactor containment building. Advances in the understanding of fission product release and transport behavior and severe accident progression are used to render best estimate analyses of selected accident sequences. Particular emphasis is placed on estimating the effects of high fuel burnup in contrast with low burnup on fission product releases to the containment. Supporting this emphasis, recent data available on fission product release from high-burnup (HBU) fuel from the French VERCOR project are used in this study. The results of these analyses are treated as samples from a population of accident sequences in order to employ approximate order statistics characterization of the results. These trends and tendencies are then compared to the NUREG-1465 alternative source term prescription used today for regulatory applications. In general, greater differences are observed between the state-of-the-art calculations for either HBU or low-burnup (LBU) fuel and the NUREG-1465 containment release fractions than exist between HBU and LBU release fractions. Current analyses suggest that retention of fission products within the vessel and the reactor coolant system (RCS) are greater than contemplated in the NUREG-1465 prescription, and that, overall, release fractions to the containment are therefore lower across the board in the present analyses than suggested in NUREG-1465. The decreased volatility of Cs2MoO4 compared to CsI or CsOH increases the predicted RCS retention of cesium, and as a result, cesium and iodine do not follow identical behaviors with respect to distribution among vessel, RCS, and containment. With respect to the regulatory alternative source term, greater differences are observed between the NUREG-1465 prescription and both HBU and LBU predictions than exist between HBU and LBU
Burnup calculation methodology in the serpent 2 Monte Carlo code
Leppaenen, J.; Isotalo, A.
2012-07-01
This paper presents two topics related to the burnup calculation capabilities in the Serpent 2 Monte Carlo code: advanced time-integration methods and improved memory management, accomplished by the use of different optimization modes. The development of the introduced methods is an important part of re-writing the Serpent source code, carried out for the purpose of extending the burnup calculation capabilities from 2D assembly-level calculations to large 3D reactor-scale problems. The progress is demonstrated by repeating a PWR test case, originally carried out in 2009 for the validation of the newly-implemented burnup calculation routines in Serpent 1. (authors)
High-burnup core design using minor actinide-containing metal fuel
Ohta, Hirokazu; Ogata, Takanari; Obara, T.
2013-07-01
A neutronic design study of metal fuel fast reactor (FR) cores is conducted on the basis of an innovative fuel design concept to achieve an extremely high burnup and realize an efficient fuel cycle system. Since it is expected that the burnup reactivity swing will become extremely large in an unprecedented high burnup core, minor actinides (MAs) from light water reactors (LWRs) are added to fresh fuel to improve the core internal conversion. Core neutronic analysis revealed that high burnups of about 200 MWd/kg for a small-scale core and about 300 MWd/kg for a large-scale core can be attained while suppressing the burnup reactivity swing to almost the same level as that of conventional cores with normal burnup. An actinide burnup analysis has shown that the MA consumption ratio is improved to about 60% and that the accumulated MAs originating from LWRs can be efficiently consumed by the high-burnup metal fuel FR. (authors)
PRESSURIZED WATER REACTOR CORE WITH PLUTONIUM BURNUP
Puechl, K.H.
1963-09-24
A pressurized water reactor is described having a core containing Pu/sup 240/ in which the effective microscopic neutronabsorption cross section of Pu/sup 240/ in unconverted condition decreases as the time of operation of the reactor increases, in order to compensate for loss of reactivity resulting from fission product buildup during reactor operation. This means serves to improve the efficiency of the reactor operation by reducing power losses resulting from control rods and burnable poisons. (AEC)
MONTE-CARLO BURNUP CALCULATION UNCERTAINTY QUANTIFICATION AND PROPAGATION DETERMINATION
Sternat, M.; Nichols, T.
2011-06-09
Reactor burnup or depletion codes are used thoroughly in the fields of nuclear forensics and nuclear safeguards. Two common codes include MONTEBURNS and MCNPX/CINDER. These are Monte-Carlo depletion routines utilizing MCNP for neutron transport calculations and either ORIGEN or CINDER for burnup calculations. Uncertainties exist in the MCNP steps, but this information is not passed to the depletion calculations or saved. To quantify this transport uncertainty and determine how it propagates between burnup steps, a statistical analysis of multiple repeated depletion runs is performed. The reactor model chosen is the Oak Ridge Research Reactor (ORR) in a single assembly, infinite lattice configuration. This model was burned for a 150 day cycle broken down into three steps. The output isotopics as well as effective multiplication factor (k-effective) were tabulated and histograms were created at each burnup step using the Scott Method to determine the bin width. The distributions for each code are a statistical benchmark and comparisons made. It was expected that the gram quantities and k-effective histograms would produce normally distributed results since they were produced from a Monte-Carlo routine, but some of the results appear to not. Statistical analyses are performed using the {chi}{sup 2} test against a normal distribution for the k-effective results and several isotopes including {sup 134}Cs, {sup 137}Cs, {sup 235}U, {sup 238}U, {sup 237}Np, {sup 238}Pu, {sup 239}Pu, and {sup 240}Pu.
MONTE-CARLO BURNUP CALCULATION UNCERTAINTY QUANTIFICATION AND PROPAGATION DETERMINATION
Nichols, T.; Sternat, M.; Charlton, W.
2011-05-08
MONTEBURNS is a Monte-Carlo depletion routine utilizing MCNP and ORIGEN 2.2. Uncertainties exist in the MCNP transport calculation, but this information is not passed to the depletion calculation in ORIGEN or saved. To quantify this transport uncertainty and determine how it propagates between burnup steps, a statistical analysis of a multiple repeated depletion runs is performed. The reactor model chosen is the Oak Ridge Research Reactor (ORR) in a single assembly, infinite lattice configuration. This model was burned for a 25.5 day cycle broken down into three steps. The output isotopics as well as effective multiplication factor (k-effective) were tabulated and histograms were created at each burnup step using the Scott Method to determine the bin width. It was expected that the gram quantities and k-effective histograms would produce normally distributed results since they were produced from a Monte-Carlo routine, but some of results do not. The standard deviation at each burnup step was consistent between fission product isotopes as expected, while the uranium isotopes created some unique results. The variation in the quantity of uranium was small enough that, from the reaction rate MCNP tally, round off error occurred producing a set of repeated results with slight variation. Statistical analyses were performed using the {chi}{sup 2} test against a normal distribution for several isotopes and the k-effective results. While the isotopes failed to reject the null hypothesis of being normally distributed, the {chi}{sup 2} statistic grew through the steps in the k-effective test. The null hypothesis was rejected in the later steps. These results suggest, for a high accuracy solution, MCNP cell material quantities less than 100 grams and greater kcode parameters are needed to minimize uncertainty propagation and minimize round off effects.
Burnup concept for a long-life fast reactor core using MCNPX.
Holschuh, Thomas Vernon,; Lewis, Tom Goslee,; Parma, Edward J.,
2013-02-01
This report describes a reactor design with a burnup concept for a long-life fast reactor core that was evaluated using Monte Carlo N-Particle eXtended (MCNPX). The current trend in advanced reactor design is the concept of a small modular reactor (SMR). However, very few of the SMR designs attempt to substantially increase the lifetime of a reactor core, especially without zone loading, fuel reshuffling, or other artificial mechanisms in the core that %E2%80%9Cflatten%E2%80%9D the power profile, including non-uniform cooling, non-uniform moderation, or strategic poison placement. Historically, the limitations of computing capabilities have prevented acceptable margins in the temporal component of the spatial excess reactivity in a reactor design, due primarily to the error in burnup calculations. This research was performed as an initial scoping analysis into the concept of a long-life fast reactor. It can be shown that a long-life fast reactor concept can be modeled using MCNPX to predict burnup and neutronics behavior. The inherent characteristic of this conceptual design is to minimize the change in reactivity over the lifetime of the reactor. This allows the reactor to operate substantially longer at full power than traditional Light Water Reactors (LWRs) or other SMR designs. For the purpose of this study, a single core design was investigated: a relatively small reactor core, yielding a medium amount of power (~200 to 400 MWth). The results of this scoping analysis were successful in providing a preliminary reactor design involving metal U-235/U-238 fuel with HT-9 fuel cladding and sodium coolant at a 20% volume fraction.
Oberle, P.; Broeders, C. H. M.; Dagan, R.
2006-07-01
The increasing tendency towards fuel lifetime extension in thermal nuclear reactors motivated validation work for available evaluation tools for nuclear fuel burnup calculations. In this study two deterministic codes with different transport solvers and one Monte Carlo method are investigated. The code system KAPROS/KARBUS uses the classical deterministic First Collision Probability method utilizing a cylinderized Wigner-Seitz cell. In the SCALES.0/TRITON/NEWT code the Extended Step Characteristic method is applied. In a first step the two deterministic codes are compared with experimental results from the KWO-Isotope Correlation Experiment up to 30 MWD/kg HM burnup, published in 1981. Two pin cell calculations are analyzed by comparison of calculated and experimental results for important heavy isotope vectors. The results are very satisfactory. Subsequently, further validation at higher burnup (< 80 MWD/kg HM) is provided by comparison of the two deterministic codes and the Monte Carlo based burnup code MONTEBURNS for PWR UO{sub 2} fuel assembly calculations. Possible reasons for differences in the results are analyzed and discussed. Especially the influence of cross section data and processing is presented. (authors)
Accident source terms for boiling water reactors with high burnup cores.
Gauntt, Randall O.; Powers, Dana Auburn; Leonard, Mark Thomas
2007-11-01
The primary objective of this report is to provide the technical basis for development of recommendations for updates to the NUREG-1465 Source Term for BWRs that will extend its applicability to accidents involving high burnup (HBU) cores. However, a secondary objective is to re-examine the fundamental characteristics of the prescription for fission product release to containment described by NUREG-1465. This secondary objective is motivated by an interest to understand the extent to which research into the release and behaviors of radionuclides under accident conditions has altered best-estimate calculations of the integral response of BWRs to severe core damage sequences and the resulting radiological source terms to containment. This report, therefore, documents specific results of fission product source term analyses that will form the basis for the HBU supplement to NUREG-1465. However, commentary is also provided on observed differences between the composite results of the source term calculations performed here and those reflected NUREG-1465 itself.
Addressing Fission Product Validation in MCNP Burnup Credit Criticality Calculations
Mueller, Don; Bowen, Douglas G; Marshall, William BJ J
2015-01-01
The US Nuclear Regulatory Commission (NRC) Division of Spent Fuel Storage and Transportation issued Interim Staff Guidance (ISG) 8, Revision 3 in September 2012. This ISG provides guidance for NRC staff members’ review of burnup credit (BUC) analyses supporting transport and dry storage of pressurized water reactor spent nuclear fuel (SNF) in casks. The ISG includes guidance for addressing validation of criticality (k_{eff}) calculations crediting the presence of a limited set of fission products and minor actinides (FP&MAs). Based on previous work documented in NRC Regulatory Guide (NUREG) Contractor Report (CR)-7109, the ISG recommends that NRC staff members accept the use of either 1.5 or 3% of the FP&MA worth—in addition to bias and bias uncertainty resulting from validation of k_{eff} calculations for the major actinides in SNF—to conservatively account for the bias and bias uncertainty associated with the specified unvalidated FP&MAs. The ISG recommends (1) use of 1.5% of the FP&MA worth if a modern version of SCALE and its nuclear data are used and (2) 3% of the FP&MA worth for well qualified, industry standard code systems other than SCALE with the Evaluated Nuclear Data Files, Part B (ENDF/B),-V, ENDF/B-VI, or ENDF/B-VII cross sections libraries. The work presented in this paper provides a basis for extending the use of the 1.5% of the FP&MA worth bias to BUC criticality calculations performed using the Monte Carlo N-Particle (MCNP) code. The extended use of the 1.5% FP&MA worth bias is shown to be acceptable by comparison of FP&MA worths calculated using SCALE and MCNP with ENDF/B-V, -VI, and -VII–based nuclear data. The comparison supports use of the 1.5% FP&MA worth bias when the MCNP code is used for criticality calculations, provided that the cask design is similar to the hypothetical generic BUC-32 cask model and that the credited FP&MA worth is no more than 0.1 Δk_{eff} (ISG-8, Rev. 3, Recommendation 4).
Calculations on fission gas behaviour in the high burnup structure
NASA Astrophysics Data System (ADS)
Blair, P.; Romano, A.; Hellwig, Ch.; Chawla, R.
2006-05-01
The behaviour of fission gas in high burnup fuel during steady-state and transient conditions is of special interest for safety reasons. Despite this, mechanistic models that reflect the fission gas transport processes and reliably predict the evolution of the remaining fission gas in the high burnup structure (HBS) are largely missing today. We start to address this problem by developing a one-dimensional, mass balance model and apply it to LWR UO 2 fuel at the moderate temperatures found in the rim region. We examine the quantity of gas remaining in the HBS fuel matrix at steady state and compare it with experimental values. We find that the current model reproduces the 0.2 wt% observed xenon concentration under certain conditions, viz. fast grain boundary diffusion and an effective volume diffusion coefficient. A sensitivity analysis is also conducted for the model parameters, the relative importance for which is not well established a priori.
Fuel rod and core materials investigations related to LWR extended burnup operation
NASA Astrophysics Data System (ADS)
Kolstad, Erik; Vitanza, Carlo
1992-06-01
The paper deals with tests and recent measurements related to extended burnup fuel performance and describes test facilities and results in the areas of waterside cladding corrosion and irradiation-assisted stress corrosion cracking (IASCC). Fuel temperature data suggest a gradual degradation of UO 2 thermal conductivity with exposure in the range 6-8% per 10 MWd/kgUO 2 at temperatures below 700°C. The effect on the fuel microstructure of interlinkage and resintering phenomena is shown by measuring the surface-to-volume ( S/ V) ratio of the fuel. Changes in S/V with burnup are correlated to power rating and fuel operating temperature. No evidence was found of enhanced fission gas release during load-follow operation in the burnup range 25-45 MWd/kgUO 2. The effect of high lithium concentration (high pH) on the corrosion behaviour of pre-irradiated high burnup Zircaloy-4 fuel rods subjected either to nucleate boiling or to one-phase cooling conditions was studied. The oxide thickness growth rates measured at an average burnup up to 40 MWd/kgUO 2 are consistent with literature data and show no evidence of corrosion enhancement due to the high lithium content and little effect of cooling regime. A test facility for exploring the effects of environmental variables on IASCC behaviour of in-core structural materials is described.
Study on core radius minimization for long life Pb-Bi cooled CANDLE burnup scheme based fast reactor
NASA Astrophysics Data System (ADS)
Afifah, Maryam; Miura, Ryosuke; Su'ud, Zaki; Takaki, Naoyuki; Sekimoto, H.
2015-09-01
Fast Breeder Reactor had been interested to be developed over the world because it inexhaustible source energy, one of those is CANDLE reactor which is have strategy in burn-up scheme, need not control roads for control burn-up, have a constant core characteristics during energy production and don't need fuel shuffling. The calculation was made by basic reactor analysis which use Sodium coolant geometry core parameter as a reference core to study on minimum core reactor radius of CANDLE for long life Pb-Bi cooled, also want to perform pure coolant effect comparison between LBE and sodium in a same geometry design. The result show that the minimum core radius of Lead Bismuth cooled CANDLE is 100 cm and 500 MWth thermal output. Lead-Bismuth coolant for CANDLE reactor enable to reduce much reactor size and have a better void coefficient than Sodium cooled as the most coolant for FBR, then we will have a good point in safety analysis.
Study on core radius minimization for long life Pb-Bi cooled CANDLE burnup scheme based fast reactor
Afifah, Maryam Su’ud, Zaki; Miura, Ryosuke; Takaki, Naoyuki; Sekimoto, H.
2015-09-30
Fast Breeder Reactor had been interested to be developed over the world because it inexhaustible source energy, one of those is CANDLE reactor which is have strategy in burn-up scheme, need not control roads for control burn-up, have a constant core characteristics during energy production and don’t need fuel shuffling. The calculation was made by basic reactor analysis which use Sodium coolant geometry core parameter as a reference core to study on minimum core reactor radius of CANDLE for long life Pb-Bi cooled, also want to perform pure coolant effect comparison between LBE and sodium in a same geometry design. The result show that the minimum core radius of Lead Bismuth cooled CANDLE is 100 cm and 500 MWth thermal output. Lead-Bismuth coolant for CANDLE reactor enable to reduce much reactor size and have a better void coefficient than Sodium cooled as the most coolant for FBR, then we will have a good point in safety analysis.
Energy Science and Technology Software Center (ESTSC)
2010-02-01
Neutron transport, calculation of multiplication factor and neutron fluxes in 2-D configurations: cell calculations, 2-D diffusion and transport, and burnup. Preparation of a cross section library for the code BOXER from a basic library in ENDF/B format (ETOBOX).
Spent fuel pool storage calculations using the ISOCRIT burnup credit tool
Kucukboyaci, Vefa; Marshall, William BJ J
2012-01-01
In order to conservatively apply burnup credit in spent fuel pool criticality safety analyses, Westinghouse has developed a software tool, ISOCRIT, for generating depletion isotopics. This tool is used to create isotopics data based on specific reactor input parameters, such as design basis assembly type; bounding power/burnup profiles; reactor specific moderator temperature profiles; pellet percent theoretical density; burnable absorbers, axial blanket regions, and bounding ppm boron concentration. ISOCRIT generates burnup dependent isotopics using PARAGON; Westinghouse's state-of-the-art and licensed lattice physics code. Generation of isotopics and passing the data to the subsequent 3D KENO calculations are performed in an automated fashion, thus reducing the chance for human error. Furthermore, ISOCRIT provides the means for responding to any customer request regarding re-analysis due to changed parameters (e.g., power uprate, exit temperature changes, etc.) with a quick turnaround.
Karpushkin, T. Yu.
2012-12-15
A technique to calculate the burnup of materials of cells and fuel assemblies using the matrices of first-flight neutron collision probabilities rebuilt at a given burnup step is presented. A method to rebuild and correct first collision probability matrices using average chords prior to the first neutron collision, which are calculated with the help of geometric modules of constructed stochastic neutron trajectories, is described. Results of calculation of the infinite multiplication factor for elementary cells with a modified material composition compared to the reference one as well as calculation of material burnup in the cells and fuel assemblies of a VVER-1000 are presented.
OECD/NEA burnup credit calculational criticality benchmark Phase I-B results
DeHart, M.D.; Parks, C.V.; Brady, M.C.
1996-06-01
In most countries, criticality analysis of LWR fuel stored in racks and casks has assumed that the fuel is fresh with the maximum allowable initial enrichment. This assumption has led to the design of widely spaced and/or highly poisoned storage and transport arrays. If credit is assumed for fuel burnup, initial enrichment limitations can be raised in existing systems, and more compact and economical arrays can be designed. Such reliance on the reduced reactivity of spent fuel for criticality control is referred to as burnup credit. The Burnup Credit Working Group, formed under the auspices of the Nuclear Energy Agency of the Organization for Economic Cooperation and Development, has established a set of well-defined calculational benchmarks designed to study significant aspects of burnup credit computational methods. These benchmarks are intended to provide a means for the intercomparison of computer codes, methods, and data applied in spent fuel analysis. The benchmarks have been divided into multiple phases, each phase focusing on a particular feature of burnup credit analysis. This report summarizes the results and findings of the Phase I-B benchmark, which was proposed to provide a comparison of the ability of different code systems and data libraries to perform depletion analysis for the prediction of spent fuel isotopic concentrations. Results included here represent 21 different sets of calculations submitted by 16 different organizations worldwide and are based on a limited set of nuclides determined to have the most important effect on the neutron multiplication factor of light-water-reactor spent fuel. A comparison of all sets of results demonstrates that most methods agree to within 10% in the ability to estimate the spent fuel concentrations of most actinides. All methods agree within 11% about the average for all fission products studied. Most deviations are less than 10%, and many are less than 5%. The exceptions are Sm 149, Sm 151, and Gd 155.
OECD/NEA Burnup Credit Calculational Criticality Benchmark Phase I-B Results
DeHart, M.D.
1993-01-01
Burnup credit is an ongoing technical concern for many countries that operate commercial nuclear power reactors. In a multinational cooperative effort to resolve burnup credit issues, a Burnup Credit Working Group has been formed under the auspices of the Nuclear Energy Agency of the Organization for Economic Cooperation and Development. This working group has established a set of well-defined calculational benchmarks designed to study significant aspects of burnup credit computational methods. These benchmarks are intended to provide a means for the intercomparison of computer codes, methods, and data applied in spent fuel analysis. The benchmarks have been divided into multiple phases, each phase focusing on a particular feature of burnup credit analysis. This report summarizes the results and findings of the Phase I-B benchmark, which was proposed to provide a comparison of the ability of different code systems and data libraries to perform depletion analysis for the prediction of spent fuel isotopic concentrations. Results included here represent 21 different sets of calculations submitted by 16 different organizations worldwide, and are based on a limited set of nuclides determined to have the most important effect on the neutron multiplication factor of light-water-reactor spent fuel. A comparison of all sets of results demonstrates that most methods are in agreement to within 10% in the ability to estimate the spent fuel concentrations of most actinides. All methods are within 11% agreement about the average for all fission products studied. Furthermore, most deviations are less than 10%, and many are less than 5%. The exceptions are {sup 149}Sm, {sup 151}Sm, and {sup 155}Gd.
Yoshikawa, T.; Iwasaki, T.; Wada, K.; Suyama, K.
2006-07-01
To examine the procedures of the reprocessing, the vitrification and the geologic disposal, precise burn-up calculation for high burn-up and MOX fuels has been performed for not only PWR but also BWR by using SWAT and SWAT2 codes which are the integrated bum-up calculation code systems combined with the bum-up calculation code, ORIGEN2, and the transport calculation code, SRAC (the collision probability method) or MVP (the continuous energy Monte Carlo method), respectively. The calculation results shows that all of the evaluated items (heat generation and concentrations of Mo and Pt) largely increase and those significantly effect to the current procedures of the vitrification and the geologic disposal. The calculation result by SWAT2 confirms that the bundle calculation is required for BWR to be discussed about those effects in details, especially for the MOX fuel. (authors)
Performance upgrades to the MCNP6 burnup capability for large scale depletion calculations
Fensin, M. L.; Galloway, J. D.; James, M. R.
2015-04-11
The first MCNP based inline Monte Carlo depletion capability was officially released from the Radiation Safety Information and Computational Center as MCNPX 2.6.0. With the merger of MCNPX and MCNP5, MCNP6 combined the capability of both simulation tools, as well as providing new advanced technology, in a single radiation transport code. The new MCNP6 depletion capability was first showcased at the International Congress for Advancements in Nuclear Power Plants (ICAPP) meeting in 2012. At that conference the new capabilities addressed included the combined distributive and shared memory parallel architecture for the burnup capability, improved memory management, physics enhancements, and newmore » predictability as compared to the H.B Robinson Benchmark. At Los Alamos National Laboratory, a special purpose cluster named “tebow,” was constructed such to maximize available RAM per CPU, as well as leveraging swap space with solid state hard drives, to allow larger scale depletion calculations (allowing for significantly more burnable regions than previously examined). As the MCNP6 burnup capability was scaled to larger numbers of burnable regions, a noticeable slowdown was realized.This paper details two specific computational performance strategies for improving calculation speedup: (1) retrieving cross sections during transport; and (2) tallying mechanisms specific to burnup in MCNP. To combat this slowdown new performance upgrades were developed and integrated into MCNP6 1.2.« less
Performance upgrades to the MCNP6 burnup capability for large scale depletion calculations
Fensin, M. L.; Galloway, J. D.; James, M. R.
2015-04-11
The first MCNP based inline Monte Carlo depletion capability was officially released from the Radiation Safety Information and Computational Center as MCNPX 2.6.0. With the merger of MCNPX and MCNP5, MCNP6 combined the capability of both simulation tools, as well as providing new advanced technology, in a single radiation transport code. The new MCNP6 depletion capability was first showcased at the International Congress for Advancements in Nuclear Power Plants (ICAPP) meeting in 2012. At that conference the new capabilities addressed included the combined distributive and shared memory parallel architecture for the burnup capability, improved memory management, physics enhancements, and new predictability as compared to the H.B Robinson Benchmark. At Los Alamos National Laboratory, a special purpose cluster named “tebow,” was constructed such to maximize available RAM per CPU, as well as leveraging swap space with solid state hard drives, to allow larger scale depletion calculations (allowing for significantly more burnable regions than previously examined). As the MCNP6 burnup capability was scaled to larger numbers of burnable regions, a noticeable slowdown was realized.This paper details two specific computational performance strategies for improving calculation speedup: (1) retrieving cross sections during transport; and (2) tallying mechanisms specific to burnup in MCNP. To combat this slowdown new performance upgrades were developed and integrated into MCNP6 1.2.
Advanced Fuel Cycle Initiative - Projected Linear Heat Generation Rate and Burnup Calculations
Richard G. Ambrosek; Gray S. Chang; Debbie J. Utterbeck
2005-02-01
This report provides documentation of the physics analysis performed to determine the linear heat generation rate (LHGR) and burnup calculations for the Advanced Fuel Cycle Initiative (AFCI) tests, AFC-1D, AFC-1H, and AFC-1G. The AFC-1D and AFC-1H tests consists of low-fertile metallic fuel compositions and the AFC-1G test consists of non-fertile and low-fertile nitride compositions. These tests will be irradiated in the East Flux Trap (EFT) positions E1, E2, and E3, respectively, during Advanced Test Reactor (ATR) Cycle 135B.
CANDLE: The New Burnup Strategy
Sekimoto, Hiroshi; Ryu, Kouichi; Yoshimura, Yoshikane
2001-11-15
The new burnup strategy CANDLE (Constant Axial shape of Neutron flux, nuclide densities and power shape During Life of Energy production) is proposed. With this burnup strategy, distributions of fuel nuclide densities, neutron flux, and power density move with the same constant speed and without any change in their shapes. The excess reactivity is constant during the burnup. Therefore, any control mechanisms for the burnup are not required. Calculation procedures are presented to find these shapes and the speed of the burning region with the neutron multiplication factor of a reactor employing this burnup strategy.To demonstrate the CANDLE burnup strategy, it is applied to a fast reactor with excellent neutron economy. Only the initially built reactor requires some fissile material such as plutonium or enriched uranium for the nuclear ignition region of its core, but only natural uranium or depleted uranium is required for the other region. Succeeding reactors require only natural or depleted uranium since the burning region of the previous reactor can be utilized for the ignition region. The life of a reactor can be made longer by elongating the core height. The drift speed of the burning region for the presented fast reactor design is {approx}4 cm/yr, which is a preferable value for designing a long-life reactor. The burnup of spent fuel is {approx}40%. It is equivalent to 40% utilization of natural uranium without reprocessing and enrichment.
Calculated Neutron and Gamma-ray Spectra across the Prismatic Very High Temperature Reactor Core
James W. Sterbentz
2008-05-01
Neutron and gamma-ray flux spectra are calculated using the MCNP5 computer code and a one-sixth core model of a prismatic Very High Temperature Reactor based on the General Atomics Gas Turbine-Modular Helium Reactor. Spectra are calculated in the five inner reflector graphite block rings, three annular active core fuel rings, three outer graphite reflector block rings, and the core barrel. The neutron spectra are block and fuel pin averages and are calculated as a function of temperature and burnup. Also provided are the total, fast, and thermal radial profile fluxes and core barrel dpa rates.
Evaluation of fission product worth margins in PWR spent nuclear fuel burnup credit calculations.
Blomquist, R.N.; Finck, P.J.; Jammes, C.; Stenberg, C.G.
1999-02-17
Current criticality safety calculations for the transportation of irradiated LWR fuel make the very conservative assumption that the fuel is fresh. This results in a very substantial overprediction of the actual k{sub eff} of the transportation casks; in certain cases, this decreases the amount of spent fuel which can be loaded in a cask, and increases the cost of transporting the spent fuel to the repository. Accounting for the change of reactivity due to fuel depletion is usually referred to as ''burnup credit.'' The US DOE is currently funding a program aimed at establishing an actinide only burnup credit methodology (in this case, the calculated reactivity takes into account the buildup or depletion of a limited number of actinides). This work is undergoing NRC review. While this methodology is being validated on a significant experimental basis, it implicitly relies on additional margins: in particular, the absorption of neutrons by certain actinides and by all fission products is not taken into account. This provides an important additional margin and helps guarantee that the methodology is conservative provided these neglected absorption are known with reasonable accuracy. This report establishes the accuracy of fission product absorption rate calculations: (1) the analysis of European fission product worth experiments demonstrates that fission product cross-sections available in the US provide very good predictions of fission product worth; (2) this is confirmed by a direct comparison of European and US cross section evaluations; (3) accuracy of Spent Nuclear Fuel (SNF) fission product content predictions is established in a recent ORNL report where several SNF isotopic assays are analyzed; and (4) these data are then combined to establish in a conservative manner the fraction of the predicted total fission product absorption which can be guaranteed based on available experimental data.
Accuracy considerations for Chebyshev rational approximation method (CRAM) in Burnup calculations
Pusa, M.
2013-07-01
The burnup equations can in principle be solved by computing the exponential of the burnup matrix. However, due to the difficult numerical characteristics of burnup matrices, the problem is extremely stiff and the matrix exponential solution has previously been considered infeasible for an entire burnup system containing over a thousand nuclides. It was recently discovered by the author that the eigenvalues of burnup matrices are generally located near the negative real axis, which prompted introducing the Chebyshev rational approximation method (CRAM) for solving the burnup equations. CRAM can be characterized as the best rational approximation on the negative real axis and it has been shown to be capable of simultaneously solving an entire burnup system both accurately and efficiently. In this paper, the accuracy of CRAM is further studied in the context of burnup equations. The approximation error is analyzed based on the eigenvalue decomposition of the burnup matrix. It is deduced that the relative accuracy of CRAM may be compromised if a nuclide concentration diminishes significantly during the considered time step. Numerical results are presented for two test cases, the first one representing a small burnup system with 36 nuclides and the second one a full a decay system with 1531 nuclides. (authors)
HTGR Reactor Physics and Burnup Calculations Using the Serpent Monte Carlo Code
Leppanen, Jaakko; DeHart, Mark D
2009-01-01
One of the main advantages of the continuous-energy Monte Carlo method is its versatility and the capability to model any fuel or reactor configuration without major approximations. This capability becomes particularly valuable in studies involving innovative reactor designs and next-generation systems, which often lie beyond the capabilities of deterministic LWR transport codes. In this study, a conceptual prismatic HTGR fuel assembly was modeled using the Serpent Monte Carlo reactor physics burnup calculation code, under development at VTT Technical Research Centre of Finland since 2004. A new explicit particle fuel model was developed to account for the heterogeneity effects. The results are compared to other Monte Carlo and deterministic transport codes and the study also serves as a test case for the modules and methods in SCALE 6.
Code System to Calculate Mixed Cores in TRIGA Mark II Research Reactor.
Energy Science and Technology Software Center (ESTSC)
2001-08-29
Version 00 TRIGLAV is a computer program for reactor calculations of mixed cores in a TRIGA Mark II research reactor. It can be applied for fuel element burn-up calculations, for power and flux distributions calculations and for reactivity predictions. The TRIGLAV program requires the WIMS-D4 program with the original WIMS cross-section library extended for TRIGA reactor specific nuclides. This package includes the code TRIGAC, which is a new version of TRIGAP.
Radaev, A. I. Schurovskaya, M. V.
2015-12-15
The choice of the spatial nodalization for the calculation of the power density and burnup distribution in a research reactor core with fuel assemblies of the IRT-3M and VVR-KN type using the program based on the Monte Carlo code is described. The influence of the spatial nodalization on the results of calculating basic neutronic characteristics and calculation time is investigated.
NASA Astrophysics Data System (ADS)
Aufiero, M.; Cammi, A.; Fiorina, C.; Leppänen, J.; Luzzi, L.; Ricotti, M. E.
2013-10-01
In this work, the Monte Carlo burn-up code SERPENT-2 has been extended and employed to study the material isotopic evolution of the Molten Salt Fast Reactor (MSFR). This promising GEN-IV nuclear reactor concept features peculiar characteristics such as the on-line fuel reprocessing, which prevents the use of commonly available burn-up codes. Besides, the presence of circulating nuclear fuel and radioactive streams from the core to the reprocessing plant requires a precise knowledge of the fuel isotopic composition during the plant operation. The developed extension of SERPENT-2 directly takes into account the effects of on-line fuel reprocessing on burn-up calculations and features a reactivity control algorithm. It is here assessed against a dedicated version of the deterministic ERANOS-based EQL3D procedure (PSI-Switzerland) and adopted to analyze the MSFR fuel salt isotopic evolution. Particular attention is devoted to study the effects of reprocessing time constants and efficiencies on the conversion ratio and the molar concentration of elements relevant for solubility issues (e.g., trivalent actinides and lanthanides). Quantities of interest for fuel handling and safety issues are investigated, including decay heat and activities of hazardous isotopes (neutron and high energy gamma emitters) in the core and in the reprocessing stream. The radiotoxicity generation is also analyzed for the MSFR nominal conditions. The production of helium and the depletion in tungsten content due to nuclear reactions are calculated for the nickel-based alloy selected as reactor structural material of the MSFR. These preliminary evaluations can be helpful in studying the radiation damage of both the primary salt container and the axial reflectors.
Propagation of Nuclear Data Uncertainties for ELECTRA Burn-up Calculations
NASA Astrophysics Data System (ADS)
Sjöstrand, H.; Alhassan, E.; Duan, J.; Gustavsson, C.; Koning, A. J.; Pomp, S.; Rochman, D.; Österlund, M.
2014-04-01
The European Lead-Cooled Training Reactor (ELECTRA) has been proposed as a training reactor for fast systems within the Swedish nuclear program. It is a low-power fast reactor cooled by pure liquid lead. In this work, we propagate the uncertainties in 239Pu transport data to uncertainties in the fuel inventory of ELECTRA during the reactor lifetime using the Total Monte Carlo approach (TMC). Within the TENDL project, nuclear models input parameters were randomized within their uncertainties and 740 239Pu nuclear data libraries were generated. These libraries are used as inputs to reactor codes, in our case SERPENT, to perform uncertainty analysis of nuclear reactor inventory during burn-up. The uncertainty in the inventory determines uncertainties in: the long-term radio-toxicity, the decay heat, the evolution of reactivity parameters, gas pressure and volatile fission product content. In this work, a methodology called fast TMC is utilized, which reduces the overall calculation time. The uncertainty of some minor actinides were observed to be rather large and therefore their impact on multiple recycling should be investigated further. It was also found that, criticality benchmarks can be used to reduce inventory uncertainties due to nuclear data. Further studies are needed to include fission yield uncertainties, more isotopes, and a larger set of benchmarks.
IN-CORE FUEL MANAGEMENT: PWR Core Calculations Using MCRAC
NASA Astrophysics Data System (ADS)
PetroviĆ, B. G.
1991-01-01
The following sections are included: * INTRODUCTION * IN-CORE FUEL MANAGEMENT CALCULATIONS * In-Core Fuel Management * Methodological Problems of In-Core Fuel Management * In-Core Fuel Management Analytical Tools * PENN STATE FUEL MANAGEMENT PACKAGE * Penn State Fuel Management Package (PFMP) * Assembly Data Description (ADD) * Linking PSU-LEOPARD and MCRAC: An Example * MULTICYCLE REACTOR ANALYSIS CODE (MCRAC) * Main Features and Options of MCRAC code * Core geometry * Diffusion equations * 1.5-group model * Multicycle neutronic analysis * Multicycle cost analysis * Criticality search * Power-dependent xenon feedback calculations * Control rod and burnable absorber simulation * Search for LP with flat BOC power distribution * Artificial ADD option * Variable dimensioning technique * RBI version of MCRAC code * Programming changes in PC version * Fuel interchange option * MCRAC Input/Output * General input description * Sample input * Sample output * EXPERIENCE WITH MCRAC CODE * CONCLUSIONS * REFERENCES
Holly R. Trellue
1998-12-01
Monteburns is a fully automated tool that links the Monte Carlo transport code MCNP with the radioactive decay and burnup code 0RIGEN2. Monteburns produces many criticality and burnup computational parameters based on material feed/removal specifications, power(s), and time intervals. This code processes input from the user indicating the system geometry, initial material compositions, feed/removal, and other code-specific parameters. Results from MCNP, 0RIGEN2, and other calculations are then output successively as the code runs. The principle function of monteburns is to first transfer one-group cross sections and fluxes from MCNP to 0RIGEN2, and then transfer the resulting material compositions (after irradiation and/or decay) from 0RIGEN2 back to MCNP in a repeated, cyclic fashion. The main requirement of the code is that the user have a working MCNP input file and other input parameters; all interaction with 0RIGEN2 and other calculations are performed by monteburns. This report presents the results obtained from the benchmarking of monteburns to measured and previously obtained data from traditional Light Water Reactor systems. The majority of the differences seen between the two were less than five percent. These were primarily a result of variances in cross sections between MCNP, cross section libraries used by other codes, and observed values. With this understanding, this code can now be used with confidence for burnup calculations in three-dimensional systems. It was designed for use in the Accelerator Transmutation of Waste project at Los Alamos National Laboratory but is also being applied to the analysis of isotopic production/destruction of transuranic actinides in a reactor system. The code has now been shown to sufficiently support these calculations.
Gohar, Y.; Zhong, Z.; Talamo, A.; Nuclear Engineering Division
2009-06-09
Argonne National Laboratory (ANL) of USA and Kharkov Institute of Physics and Technology (KIPT) of Ukraine have been collaborating on the conceptual design development of an electron accelerator driven subcritical (ADS) facility, using the KIPT electron accelerator. The neutron source of the subcritical assembly is generated from the interaction of 100 KW electron beam with a natural uranium target. The electron beam has a uniform spatial distribution and electron energy in the range of 100 to 200 MeV. The main functions of the subcritical assembly are the production of medical isotopes and the support of the Ukraine nuclear power industry. Neutron physics experiments and material structure analyses are planned using this facility. With the 100 KW electron beam power, the total thermal power of the facility is {approx}375 kW including the fission power of {approx}260 kW. The burnup of the fissile materials and the buildup of fission products reduce continuously the reactivity during the operation, which reduces the neutron flux level and consequently the facility performance. To preserve the neutron flux level during the operation, fuel assemblies should be added after long operating periods to compensate for the lost reactivity. This process requires accurate prediction of the fuel burnup, the decay behavior of the fission produces, and the introduced reactivity from adding fresh fuel assemblies. The recent developments of the Monte Carlo computer codes, the high speed capability of the computer processors, and the parallel computation techniques made it possible to perform three-dimensional detailed burnup simulations. A full detailed three-dimensional geometrical model is used for the burnup simulations with continuous energy nuclear data libraries for the transport calculations and 63-multigroup or one group cross sections libraries for the depletion calculations. Monte Carlo Computer code MCNPX and MCB are utilized for this study. MCNPX transports the
Nejat, S.M.R. . Dept. of Engineering Physics.)
1993-08-01
The 5-MW highly enriched uranium (HEU)-fueled Tehran Research Reactor is considered for conversion to high-density, low-enriched uranium (LEU) fuel. A preliminary neutronics calculation is performed as part of the conversion goal. In this study, two cores are considered: the HEU reference core and a proposed LEU core similar to the reference core, and a proposed LEU core similar to the reference core, using standardized U[sub 3]Si[sub 2] plates with the option of different [sup 235]U loadings. Various possibilities are investigated for the conversion of HEU to LEU fuel elements with 20% enriched [sup 235]U loadings of 207 to 297 g [sup 235]U/element. For the same equilibrium cycle length, the fuels are compared for flux, power distribution, burnup, and reactivity.
NASA Astrophysics Data System (ADS)
Borodkin, Pavel; Borodkin, Gennady; Khrennikov, Nikolay
2016-02-01
Paper describes the new approach of fitting axial fuel burn-up patterns in peripheral fuel assemblies of VVER-1000 type reactors, on the base of ex-core neutron leakage measurements, neutron-physical calculations and in-core SPND measured data. The developed approach uses results of new ex-vessel measurements on different power units through different reactor cycles and their uncertainties to clear the influence of a fitted fuel burn-up profile to the RPV neutron fluence calculations. The new methodology may be recommended to be included in the routine fluence calculations used in RPV lifetime management and may be taken into account during VVER-1000 core burn-up pattern correction.
Borodkin, P.G.; Borodkin, G.I.; Khrennikov, N.N.
2011-07-01
This paper deals with calculated and semi-analytical evaluations of VVER-1000 reactor core neutron source distributions and their influence on measurements and calculations of the integral through-vessel neutron leakage. Neutron activation measurements analyzed in the paper were carried out in an ex-vessel air cavity at different nuclear power plant units with VVER-1000 during different fuel cycles. The time-integrated neutron source distributions used for DORT calculations were prepared via two different approaches based on (a) calculated fuel burnup (standard routine procedure) and (b) in-core measurements by means of self-powered detectors (SPDs) and thermocouples (TCs) (new approach). Considering that fuel burnup distributions in operating VVER may be evaluated now by the use of analytical methods (calculations) only, it is necessary to develop new approaches for the testing and correction of calculated evaluations of a neutron source. The results presented in this paper allow one to consider the reverse task of the alternative estimation of fuel burnup distributions. The proposed approach is based on the adjustment (fitting) of time-integrated neutron source distributions, and thus fuel burnup patterns, in some part of the reactor core, taking into account neutron leakage measurements, neutron-physical calculations, and in-core SPD and TC measurement data. (authors)
ORPHEE research reactor: 3D core depletion calculation using Monte-Carlo code TRIPOLI-4®
NASA Astrophysics Data System (ADS)
Damian, F.; Brun, E.
2014-06-01
ORPHEE is a research reactor located at CEA Saclay. It aims at producing neutron beams for experiments. This is a pool-type reactor (heavy water), and the core is cooled by light water. Its thermal power is 14 MW. ORPHEE core is 90 cm height and has a cross section of 27x27 cm2. It is loaded with eight fuel assemblies characterized by a various number of fuel plates. The fuel plate is composed of aluminium and High Enriched Uranium (HEU). It is a once through core with a fuel cycle length of approximately 100 Equivalent Full Power Days (EFPD) and with a maximum burnup of 40%. Various analyses under progress at CEA concern the determination of the core neutronic parameters during irradiation. Taking into consideration the geometrical complexity of the core and the quasi absence of thermal feedback for nominal operation, the 3D core depletion calculations are performed using the Monte-Carlo code TRIPOLI-4® [1,2,3]. A preliminary validation of the depletion calculation was performed on a 2D core configuration by comparison with the deterministic transport code APOLLO2 [4]. The analysis showed the reliability of TRIPOLI-4® to calculate a complex core configuration using a large number of depleting regions with a high level of confidence.
Testing actinide fission yield treatment in CINDER90 for use in MCNP6 burnup calculations
Fensin, Michael Lorne; Umbel, Marissa
2015-09-18
Most of the development of the MCNPX/6 burnup capability focused on features that were applied to the Boltzman transport or used to prepare coefficients for use in CINDER90, with little change to CINDER90 or the CINDER90 data. Though a scheme exists for best solving the coupled Boltzman and Bateman equations, the most significant approximation is that the employed nuclear data are correct and complete. Thus, the CINDER90 library file contains 60 different actinide fission yields encompassing 36 fissionable actinides (thermal, fast, high energy and spontaneous fission). Fission reaction data exists for more than 60 actinides and as a result, fissionmore » yield data must be approximated for actinides that do not possess fission yield information. Several types of approximations are used for estimating fission yields for actinides which do not possess explicit fission yield data. The objective of this study is to test whether or not certain approximations of fission yield selection have any impact on predictability of major actinides and fission products. Further we assess which other fission products, available in MCNP6 Tier 3, result in the largest difference in production. Because the CINDER90 library file is in ASCII format and therefore easily amendable, we assess reasons for choosing, as well as compare actinide and major fission product prediction for the H. B. Robinson benchmark for, three separate fission yield selection methods: (1) the current CINDER90 library file method (Base); (2) the element method (Element); and (3) the isobar method (Isobar). Results show that the three methods tested result in similar prediction of major actinides, Tc-99 and Cs-137; however, certain fission products resulted in significantly different production depending on the method of choice.« less
Testing actinide fission yield treatment in CINDER90 for use in MCNP6 burnup calculations
Fensin, Michael Lorne; Umbel, Marissa
2015-09-18
Most of the development of the MCNPX/6 burnup capability focused on features that were applied to the Boltzman transport or used to prepare coefficients for use in CINDER90, with little change to CINDER90 or the CINDER90 data. Though a scheme exists for best solving the coupled Boltzman and Bateman equations, the most significant approximation is that the employed nuclear data are correct and complete. Thus, the CINDER90 library file contains 60 different actinide fission yields encompassing 36 fissionable actinides (thermal, fast, high energy and spontaneous fission). Fission reaction data exists for more than 60 actinides and as a result, fission yield data must be approximated for actinides that do not possess fission yield information. Several types of approximations are used for estimating fission yields for actinides which do not possess explicit fission yield data. The objective of this study is to test whether or not certain approximations of fission yield selection have any impact on predictability of major actinides and fission products. Further we assess which other fission products, available in MCNP6 Tier 3, result in the largest difference in production. Because the CINDER90 library file is in ASCII format and therefore easily amendable, we assess reasons for choosing, as well as compare actinide and major fission product prediction for the H. B. Robinson benchmark for, three separate fission yield selection methods: (1) the current CINDER90 library file method (Base); (2) the element method (Element); and (3) the isobar method (Isobar). Results show that the three methods tested result in similar prediction of major actinides, Tc-99 and Cs-137; however, certain fission products resulted in significantly different production depending on the method of choice.
Power excursion analysis for BWR`s at high burnup
Diamond, D.J.; Neymoith, L.; Kohut, P.
1996-03-01
A study has been undertaken to determine the fuel enthalpy during a rod drop accident and during two thermal-hydraulic transients. The objective was to understand the consequences to high burnup fuel and the sources of uncertainty in the calculations. The analysis was done with RAMONA-4B, a computer code that models the neutron kinetics throughout the core along with the thermal-hydraulics in the core, vessel, and steamline. The results showed that the maximum fuel enthalpy in high burnup fuel will be affected by core design, initial conditions, and modeling assumptions. The important parameters in each of these categories are discussed in the paper.
Moderator poison design and burn-up calculations at the SNS
NASA Astrophysics Data System (ADS)
Lu, W.; Ferguson, P. D.; Iverson, E. B.; Gallmeier, F. X.; Popova, I.
2008-06-01
The spallation neutron source (SNS) at Oak Ridge National Laboratory was commissioned in April 2006. At the nominal operating power (1.4 MW), it will have thermal neutron fluxes approximately an order of magnitude greater than any existing pulsed spallation source. It thus brings a serious challenge to the lifetime of the moderator poison sheets. The SNS moderators are integrated with the inner reflector plug (IRP) at a cost of ˜$2 million a piece. A replacement of the inner reflector plug presents a significant drawback to the facility due to the activation and the operation cost. Although there are a lot of factors limiting the lifetime of the inner reflector plug, like radiation damage to the structural material and helium production of beryllium, the bottle-neck is the lifetime of the moderator poison sheets. Increasing the thickness of the poison sheet extends the lifetime but would sacrifice the neutronic performance of the moderators. A compromise is accepted at the current SNS target system which uses thick Gd poison sheets at a projected lifetime of 6 MW-years of operation. The calculations in this paper reveal that Cd may be a better poison material from the perspective of lifetime and neutronic performance. In replacing Gd, the inner reflector plug could reach a lifetime of 8 MW-years with ˜5% higher peak neutron fluxes at almost no loss of energy resolution.
SOURCE OF BURNUP VALUES FOR COMMERCIAL SPENT NUCLEAR FUEL ASSEMBLIES
BSC
2004-12-01
Waste packages are loaded with commercial spent nuclear fuel (SNF) that satisfies the minimum burnup requirements of a criticality loading curve. The burnup value assigned by the originating nuclear utility to each SNF assembly (assigned burnup) is used to load waste packages in compliance with a criticality loading curve. The burnup provided by a nuclear utility has uncertainties, so conservative calculation methods are used to characterize those uncertainties for incorporation into the criticality loading curves. Procedural safety controls ensure that the correct assembly is loaded into each waste package to prevent a misload that could create a condition affecting the safety margins. Probabilistic analyses show that procedural safety controls can minimize the chance of a misload but can not completely eliminate the possibility. Physical measurements of burnup with instrumentation in the surface facility are not necessary due to the conservative calculation methods used to produce the criticality loading curves. The reactor records assigned burnup of a commercial SNF assembly contains about two percent uncertainty, which is increased to five-percent to ensure conservatism. This five-percent uncertainty is accommodated by adjusting the criticality loading curve. Also, the record keeping methods of nuclear utilities are not uniform and the level of detail required by the NRC has varied over the last several decades. Thus, some SNF assemblies may have assigned burnups that are averages for a batch of assemblies with similar characteristics. Utilities typically have access to more detailed core-follow records that allow the batch average burnup to be changed to an assembly specific burnup. Alternatively, an additional safety margin is incorporated into the criticality loading curve to accommodate SNF assemblies with batch average burnups or greater uncertainties due to the methodology used by the nuclear utility. The utility records provide the assembly identifier
Core Physics and Kinetics Calculations for the Fissioning Plasma Core Reactor
NASA Technical Reports Server (NTRS)
Butler, C.; Albright, D.
2007-01-01
Highly efficient, compact nuclear reactors would provide high specific impulse spacecraft propulsion. This analysis and numerical simulation effort has focused on the technical feasibility issues related to the nuclear design characteristics of a novel reactor design. The Fissioning Plasma Core Reactor (FPCR) is a shockwave-driven gaseous-core nuclear reactor, which uses Magneto Hydrodynamic effects to generate electric power to be used for propulsion. The nuclear design of the system depends on two major calculations: core physics calculations and kinetics calculations. Presently, core physics calculations have concentrated on the use of the MCNP4C code. However, initial results from other codes such as COMBINE/VENTURE and SCALE4a. are also shown. Several significant modifications were made to the ISR-developed QCALC1 kinetics analysis code. These modifications include testing the state of the core materials, an improvement to the calculation of the material properties of the core, the addition of an adiabatic core temperature model and improvement of the first order reactivity correction model. The accuracy of these modifications has been verified, and the accuracy of the point-core kinetics model used by the QCALC1 code has also been validated. Previously calculated kinetics results for the FPCR were described in the ISR report, "QCALC1: A code for FPCR Kinetics Model Feasibility Analysis" dated June 1, 2002.
Calculation on flux-MMF relationship of orthogonal-core
Tajima, K.; Kaga, A.; Anazawa, Y. ); Ichinokura, O. )
1993-03-01
Orthogonal-cores have various potential applications, for instance in parametric transformers and dc-ac converters. The operating characteristics of the devices can be calculated on the basis of the measured relationship of flux to MMF of the orthogonal-core. To achieve optimal design of the applied device, the relationship of flux to MMF must be determined; however, this involves solving a three dimensional nonlinear problem. In this paper, the authors calculate the flux-MMF relationship based on a magnetic circuit model for the orthogonal-core. The computed results agree well with experiment. The method of this study is shown to be valid for calculation of characteristics and useful for optimal design of application devices.
Bai, D.; Levine, S.L. ); Luoma, J.; Mahgerefteh, M. )
1992-01-01
The Three Mile Island unit 1 core reloads have been designed using fast but accurate scoping codes, PSUI-LEOPARD and ADMARC. PSUI-LEOPARD has been normalized to EPRI-CPM2 results and used to calculate the two-group constants, whereas ADMARC is a modern two-dimensional, two-group diffusion theory nodal code. Problems in accuracy were encountered for cycles 8 and higher as the core lifetime was increased beyond 500 effective full-power days. This is because the heavier loaded cores in both {sup 235}U and {sup 10}B have harder neutron spectra, which produces a change in the transport effect in the baffle reflector region, and the burnable poison (BP) simulations were not accurate enough for the cores containing the increased amount of {sup 10}B required in the BP rods. In the authors study, a technique has been developed to take into account the change in the transport effect in the baffle region by modifying the fast neutron diffusion coefficient as a function of cycle length and core exposure or burnup. A more accurate BP simulation method is also developed, using integral transport theory and CPM2 data, to calculate the BP contribution to the equivalent fuel assembly (supercell) two-group constants. The net result is that the accuracy of the scoping codes is as good as that produced by CASMO/SIMULATE or CPM2/SIMULATE when comparing with measured data.
PWR AXIAL BURNUP PROFILE ANALYSIS
J.M. Acaglione
2003-09-17
The purpose of this activity is to develop a representative ''limiting'' axial burnup profile for pressurized water reactors (PWRs), which would encompass the isotopic axial variations caused by different assembly irradiation histories, and produce conservative isotopics with respect to criticality. The effect that the low burnup regions near the ends of spent fuel have on system reactivity is termed the ''end-effect''. This calculation will quantify the end-effects associated with Pressurized Water Reactor (PWR) fuel assemblies emplaced in a hypothetical 21 PWR waste package. The scope of this calculation covers an initial enrichment range of 3.0 through 5.0 wt% U-235 and a burnup range of 10 through 50 GWd/MTU. This activity supports the validation of the process for ensuring conservative generation of spent fuel isotopics with respect to criticality safety applications, and the use of burnup credit for commercial spent nuclear fuel. The intended use of these results will be in the development of PWR waste package loading curves, and applications involving burnup credit. Limitations of this evaluation are that the limiting profiles are only confirmed for use with the B&W 15 x 15 fuel assembly design. However, this assembly design is considered bounding of all other typical commercial PWR fuel assembly designs. This calculation is subject to the Quality Assurance Requirements and Description (QARD) because this activity supports investigations of items or barriers on the Q-list (YMP 2001).
Berna, G.A.; Beyer, G.A.; Davis, K.L.; Lanning, D.D.
1997-12-01
FRAPCON-3 is a FORTRAN IV computer code that calculates the steady-state response of light water reactor fuel rods during long-term burnup. The code calculates the temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include (1) heat conduction through the fuel and cladding, (2) cladding elastic and plastic deformation, (3) fuel-cladding mechanical interaction, (4) fission gas release, (5) fuel rod internal gas pressure, (6) heat transfer between fuel and cladding, (7) cladding oxidation, and (8) heat transfer from cladding to coolant. The code contains necessary material properties, water properties, and heat-transfer correlations. The codes` integral predictions of mechanical behavior have not been assessed against a data base, e.g., cladding strain or failure data. Therefore, it is recommended that the code not be used for analyses of cladding stress or strain. FRAPCON-3 is programmed for use on both mainframe computers and UNIX-based workstations such as DEC 5000 or SUN Sparcstation 10. It is also programmed for personal computers with FORTRAN compiler software and at least 8 to 10 megabytes of random access memory (RAM). The FRAPCON-3 code is designed to generate initial conditions for transient fuel rod analysis by the FRAPTRAN computer code (formerly named FRAP-T6).
Mueller, Don E.; Marshall, William J.; Wagner, John C.; Bowen, Douglas G.
2015-09-01
The U.S. Nuclear Regulatory Commission (NRC) Division of Spent Fuel Storage and Transportation recently issued Interim Staff Guidance (ISG) 8, Revision 3. This ISG provides guidance for burnup credit (BUC) analyses supporting transport and storage of PWR pressurized water reactor (PWR) fuel in casks. Revision 3 includes guidance for addressing validation of criticality (k_{eff}) calculations crediting the presence of a limited set of fission products and minor actinides (FP&MA). Based on previous work documented in NUREG/CR-7109, recommendation 4 of ISG-8, Rev. 3, includes a recommendation to use 1.5 or 3% of the FP&MA worth to conservatively cover the bias due to the specified FP&MAs. This bias is supplementary to the bias and bias uncertainty resulting from validation of k_{eff} calculations for the major actinides in SNF and does not address extension to actinides and fission products beyond those identified herein. The work described in this report involves comparison of FP&MA worths calculated using SCALE and MCNP with ENDF/B-V, -VI, and -VII based nuclear data and supports use of the 1.5% FP&MA worth bias when either SCALE or MCNP codes are used for criticality calculations, provided the other conditions of the recommendation 4 are met. The method used in this report may also be applied to demonstrate the applicability of the 1.5% FP&MA worth bias to other codes using ENDF/B V, VI or VII based nuclear data. The method involves use of the applicant s computational method to generate FP&MA worths for a reference SNF cask model using specified spent fuel compositions. The applicant s FP&MA worths are then compared to reference values provided in this report. The applicants FP&MA worths should not exceed the reference results by more than 1.5% of the reference FP&MA worths.
Recent Developments in No-Core Shell-Model Calculations
Navratil, P; Quaglioni, S; Stetcu, I; Barrett, B R
2009-03-20
We present an overview of recent results and developments of the no-core shell model (NCSM), an ab initio approach to the nuclear many-body problem for light nuclei. In this aproach, we start from realistic two-nucleon or two- plus three-nucleon interactions. Many-body calculations are performed using a finite harmonic-oscillator (HO) basis. To facilitate convergence for realistic inter-nucleon interactions that generate strong short-range correlations, we derive effective interactions by unitary transformations that are tailored to the HO basis truncation. For soft realistic interactions this might not be necessary. If that is the case, the NCSM calculations are variational. In either case, the ab initio NCSM preserves translational invariance of the nuclear many-body problem. In this review, we, in particular, highlight results obtained with the chiral two- plus three-nucleon interactions. We discuss efforts to extend the applicability of the NCSM to heavier nuclei and larger model spaces using importance-truncation schemes and/or use of effective interactions with a core. We outline an extension of the ab initio NCSM to the description of nuclear reactions by the resonating group method technique. A future direction of the approach, the ab initio NCSM with continuum, which will provide a complete description of nuclei as open systems with coupling of bound and continuum states is given in the concluding part of the review.
Topical report on actinide-only burnup credit for PWR spent nuclear fuel packages. Revision 1
None, None
1997-04-01
A methodology for performing and applying nuclear criticality safety calculations, for PWR spent nuclear fuel (SNF) packages with actinide-only burnup credit, is described. The changes in the U-234, U-235, U-236, U-238, Pu-238, Pu-239, Pu-240, Pu-241, Pu-242, and Am-241 concentration with burnup are used in burnup credit criticality analyses. No credit for fission product neutron absorbers is taken. The methodology consists of five major steps. (1) Validate a computer code system to calculate isotopic concentrations of SNF created during burnup in the reactor core and subsequent decay. A set of chemical assay benchmarks is presented for this purpose as well as a method for assessing the calculational bias and uncertainty, and conservative correction factors for each isotope. (2) Validate a computer code system to predict the subcritical multiplication factor, k{sub eff}, of a spent nuclear fuel package. Fifty-seven UO{sub 2}, UO{sub 2}/Gd{sub 2}O{sub 3}, and UO{sub 2}/PuO{sub 2} critical experiments have been selected to cover anticipated conditions of SNF. The method uses an upper safety limit on k{sub eff} (which can be a function of the trending parameters) such that the biased k{sub eff}, when increased for the uncertainty is less than 0.95. (3) Establish bounding conditions for the isotopic concentration and criticality calculations. Three bounding axial profiles have been established to assure the ''end effect'' is accounted for conservatively. (4) Use the validated codes and bounding conditions to generate package loading criteria (burnup credit loading curves). Burnup credit loading curves show the minimum burnup required for a given initial enrichment. The utility burnup record is compared to this requirement after the utility accounts for the uncertainty in its record. Separate curves may be generated for each assembly design, various minimum cooling times and burnable absorber histories. (5) Verify that SNF assemblies meet the package loading criteria
Effect of Fuel Fraction on Small Modified CANDLE Burn-up Based Gas Cooled Fast Reactors
Ariani, Menik; Su'ud, Zaki; Waris, Abdul; Asiah, Nur; Shafii, M. Ali; Khairurrijal
2010-12-23
A conceptual design study of Gas Cooled Fast Reactors with Modified CANDLE Burn-up has been performed. The objective of this research is to get optimal design parameters of such type reactors. The parameters of nuclear design including the critical condition, conversion ratio, and burn-up level were compared. These parameters are calculated by variation in the fuel fraction 47.5% up to 70%. Two dimensional full core multi groups diffusion calculations was performed by CITATION code. Group constant preparations are performed by using SRAC code system with JENDL-3.2 nuclear data library. In this design the reactor cores with cylindrical cell two dimensional R-Z core models are subdivided into several parts with the same volume in the axial directions. The placement of fuel in core arranged so that the result of plutonium from natural uranium can be utilized optimally for 10 years reactor operation. Modified CANDLE burn-up was established successfully in a core radial width 1.4 m. Total thermal power output for reference core is 550 MW. Study on the effect of fuel to coolant ratio shows that effective multiplication factor (k{sub eff}) is in almost linear relations with the change of the fuel volume to coolant ratio.
Effect of Fuel Fraction on Small Modified CANDLE Burn-up Based Gas Cooled Fast Reactors
NASA Astrophysics Data System (ADS)
Ariani, Menik; Su'ud, Zaki; Waris, Abdul; Khairurrijal, Asiah, Nur; Shafii, M. Ali
2010-12-01
A conceptual design study of Gas Cooled Fast Reactors with Modified CANDLE Burn-up has been performed. The objective of this research is to get optimal design parameters of such type reactors. The parameters of nuclear design including the critical condition, conversion ratio, and burn-up level were compared. These parameters are calculated by variation in the fuel fraction 47.5% up to 70%. Two dimensional full core multi groups diffusion calculations was performed by CITATION code. Group constant preparations are performed by using SRAC code system with JENDL-3.2 nuclear data library. In this design the reactor cores with cylindrical cell two dimensional R-Z core models are subdivided into several parts with the same volume in the axial directions. The placement of fuel in core arranged so that the result of plutonium from natural uranium can be utilized optimally for 10 years reactor operation. Modified CANDLE burn-up was established successfully in a core radial width 1.4 m. Total thermal power output for reference core is 550 MW. Study on the effect of fuel to coolant ratio shows that effective multiplication factor (keff) is in almost linear relations with the change of the fuel volume to coolant ratio.
NASA Astrophysics Data System (ADS)
Liu, Xiaojing; Wang, Youhua; Zhu, Jianguo; Guo, Youguang; Lei, Gang; Liu, Chengcheng
2016-05-01
Amorphous and nanocrystalline alloys are now widely used for the cores of high-frequency transformers, and Litz-wire is commonly used as the windings, while it is difficult to calculate the resistance accurately. In order to design a high-frequency transformer, it is important to accurately calculate the core loss and copper loss. To calculate the core loss accurately, the additional core loss by the effect of end stripe should be considered. It is difficult to simulate the whole stripes in the core due to the limit of computation, so a scale down model with 5 stripes of amorphous alloy is simulated by the 2D finite element method (FEM). An analytical model is presented to calculate the copper loss in the Litz-wire, and the results are compared with the calculations by FEM.
Parish, T.A.
1995-03-02
This project involved performing the nuclear design and safety analyses needed to modify the license issued by the Nuclear Regulatory Commission to allow operation of the Texas A& M University Nuclear Science Center Reactor (NSCR) with a core containing low enrichment uranium (LEU) fuel. The specific type of LEU fuel to be considered was the TRIGA 20-20 fuel produced by General Atomic. Computer codes for the neutronic analyses were provided by Argonne National Laboratory (ANL) and the assistance of William Woodruff of ANL in helping the NSCR staff to learn the proper use of the codes is gratefully acknowledged. The codes applied in the LEU analyses were WIMSd4/m, DIF3D, NCTRIGA and PARET. These codes allowed full three dimensional, temperature and burnup dependent calculations modelling the NSCR core to be performed for the first time. In addition, temperature coefficients of reactivity and pulsing calculations were carried out in-house, whereas in the past this modelling had been performed at General Atomic. In order to benchmark the newly acquired codes, modelling of the current NSCR core with highly enriched uranium fuel was also carried out. Calculated results were compared to both earlier licensing calculations and experimental data and the new methods were found to achieve excellent agreement with both. Therefore, even if an LEU core is never loaded at the NSCR, this project has resulted in a significant improvement in the nuclear safety analysis capabilities established and maintained at the NSCR.
Preliminary Study of Burnup Characteristics for a Simplified Small Pebble Bed Reactor
NASA Astrophysics Data System (ADS)
Irwanto, Dwi; Kato, Yukikata; Yamanaka, Ichiro; Obara, Toru
2010-06-01
Simplification of the pebble bed reactor by removing the unloading device from the system was peformed. For this reactor design, a suitable fuel-loading scheme is the Peu à Peu (little by little) fueling scheme. In the Peu à Peu modus, there is no unloading device; as such, the fuels are never discharged and remain at the bottom of the core during reactor operation. This means that the burnup cycle and reactivity is controlled by the addition of fuel. The objectives of the the present study were to find a means of carrying out the exact calculations needed to analyze the Peu à Peu fuel-loading scheme and to optimize the fuel composition, and fuel-loading scheme to achieve better burnup characteristics. The Monte Carlo method is used to perform calculations with high accuracy. Before the calculation of the whole core, the analysis for the infinite geometry was performed. The power generated per mass consumed for each combination of the uranium enrichment and packing fraction was analyzed from the parametric survey. By using the optimal value obtained, a whole-core calculation for the small 20 MWth reactor was performed and the criticality and burnup of this design was analyzed.
Preliminary Study of Burnup Characteristics for a Simplified Small Pebble Bed Reactor
Irwanto, Dwi; Kato, Yukikata; Obara, Toru; Yamanaka, Ichiro
2010-06-22
Simplification of the pebble bed reactor by removing the unloading device from the system was peformed. For this reactor design, a suitable fuel-loading scheme is the Peu a Peu (little by little) fueling scheme. In the Peu a Peu modus, there is no unloading device; as such, the fuels are never discharged and remain at the bottom of the core during reactor operation. This means that the burnup cycle and reactivity is controlled by the addition of fuel. The objectives of the the present study were to find a means of carrying out the exact calculations needed to analyze the Peu a Peu fuel-loading scheme and to optimize the fuel composition, and fuel-loading scheme to achieve better burnup characteristics. The Monte Carlo method is used to perform calculations with high accuracy. Before the calculation of the whole core, the analysis for the infinite geometry was performed. The power generated per mass consumed for each combination of the uranium enrichment and packing fraction was analyzed from the parametric survey. By using the optimal value obtained, a whole-core calculation for the small 20 MWth reactor was performed and the criticality and burnup of this design was analyzed.
Marshall, William BJ J; Ade, Brian J; Bowman, Stephen M; Gauld, Ian C; Ilas, Germina; Mertyurek, Ugur; Radulescu, Georgeta
2015-01-01
Oak Ridge National Laboratory and the United States Nuclear Regulatory Commission have initiated a multiyear project to investigate application of burnup credit for boiling-water reactor (BWR) fuel in storage and transportation casks. This project includes two phases. The first phase (1) investigates applicability of peak reactivity methods currently used in spent fuel pools (SFPs) to storage and transportation systems and (2) evaluates validation of both reactivity (k_{eff}) calculations and burnup credit nuclide concentrations within these methods. The second phase will focus on extending burnup credit beyond peak reactivity. This paper documents the first phase, including an analysis of lattice design parameters and depletion effects, as well as both validation components. Initial efforts related to extended burnup credit are discussed in a companion paper. Peak reactivity analyses have been used in criticality analyses for licensing of BWR fuel in SFPs over the last 20 years. These analyses typically combine credit for the gadolinium burnable absorber present in the fuel with a modest amount of burnup credit. Gadolinium burnable absorbers are used in BWR assemblies to control core reactivity. The burnable absorber significantly reduces assembly reactivity at beginning of life, potentially leading to significant increases in assembly reactivity for burnups less than 15–20 GWd/MTU. The reactivity of each fuel lattice is dependent on gadolinium loading. The number of gadolinium-bearing fuel pins lowers initial lattice reactivity, but it has a small impact on the burnup and reactivity of the peak. The gadolinium concentration in each pin has a small impact on initial lattice reactivity but a significant effect on the reactivity of the peak and the burnup at which the peak occurs. The importance of the lattice parameters and depletion conditions are primarily determined by their impact on the gadolinium depletion. Criticality code validation for BWR burnup
Advances in core loss calculations for magnetic materials
NASA Technical Reports Server (NTRS)
Triner, J. E.
1982-01-01
A new analytical technique which predicts the basic magnetic properties under various operating conditions encountered in state-of-the-art dc-ac/dc converters is discussed. Using a new flux-controlled core excitation circuit, magnetic core characteristics were developed for constant values of ramp flux (square wave voltage excitation) and frequency. From this empirical data, a mathematical loss characteristics equation is developed to analytically predict the specific core loss of several magnetic materials under various waveform excitation conditions. In addition, these characteristics show the circuit designer for the first time the direct functional relatonships between induction level and specific core loss as a function of the two key dc-dc converter operating parameters of input voltage and duty cycle.
Calculation methods for core distortions and mechanical behavior
Sutherland, W.H.
1984-09-01
This paper describes ABADAN, a general purpose, nonlinear, multi-dimensional finite element structural analyses computer code developed for the express purpose of solving large nonlinear problems as typified by the Liquid Metal Fast Breeder Reactor (LMFBR) Core Restraint System design problem. All of the structural modeling features inherent in a general purpose finite element code and required to adequately model an LMFBR core restraint system are demonstrated. Typical results for a radial row and a sixty degree sector model of FFTF are presented. The sixty degree sector results are interpreted in terms of the design criteria that the core restraint system must satisfy. Extensions and adaptations of these modeling techniques to different core restraint design concepts can be readily achieved. 27 figures.
Yasunori Ohoka; Ismile; Hiroshi Sekimoto
2004-07-01
The CANDLE burnup strategy is a new reactor burnup concept, where the distributions of fuel nuclide densities, neutron flux, and power density move with the same constant speed along the core axis from bottom to top or from top to bottom of the core and without any change in their shapes. It can be applied easily to the block-type high temperature gas cooled reactor using an appropriate burnable poison mixed with uranium oxide fuel. In the present study, the burnup distribution and the temperature distribution in the core are investigated and their effects on the CANDLE burnup core characteristics are studied. In this study, the natural gadolinium is used as the burnable poison. With the fuel enrichment of 15%, the natural gadolinium concentration of 3.0% and the fuel pin pitch of 6.6 cm, the CANDLE burnup is realized with the burning region moving speed of 29 cm/year and the axial half width of power density distribution of 1.5 m for uniform group constant case at 900 K. When the effect of nuclide change by burnup is considered, the burning region speed becomes 25 cm/year and the axial half-width of power density distribution becomes 1.25 m. When the temperature distributions effect is considered, the effects on the core characteristics are smaller than the burnup distribution effect. The maximum fuel temperature of the parallel flow case is higher than the counter flow case. (authors)
One method for calculating flux-MMF relationship of orthogonal-core
Tajima, Katsubumi; Kaga, Akio; Anazawa, Yoshihisa . Dept. of Electrical and Electronics Engineering); Ichinokura, Osamu . Dept. of Electrical Engineering)
1993-11-01
The orthogonal-core has various applications, e.g. as a variable inductor, a parametric transformer and a DC-AC converter. This paper describes one method for calculating the flux-MMF relationship of an orthogonal-core. The calculation is based on a 3-dimensional magnetic circuit model of the orthogonal-core. The model is derived by dividing the orthogonal-core, inclusive of the surrounding region, into elements comprising a 3-dimensional magnetic circuit. Using this model, the authors can compute the flux-MMF relationship of the orthogonal core with arbitrary dimensions from the B-H characteristic of the core material. The calculation method presented here is useful for optimum design of devices using an orthogonal-core.
Naruk, S.J.
1987-07-01
Minimum offset of 7 km across the Pinaleno Mountains metamorphic core complex is calculated by integrating the shear strains across the exposed width of the mylonite zone. The calculated displacement equals the offset on the associated detachment fault, estimated from offset marker beds. The method of determining displacement by strain integration may be directly applicable to many other metamorphic core complexes.
NASA Astrophysics Data System (ADS)
Sloma, Tanya Noel
When representing the behavior of commercial spent nuclear fuel (SNF), credit is sought for the reduced reactivity associated with the net depletion of fissile isotopes and the creation of neutron-absorbing isotopes, a process that begins when a commercial nuclear reactor is first operated at power. Burnup credit accounts for the reduced reactivity potential of a fuel assembly and varies with the fuel burnup, cooling time, and the initial enrichment of fissile material in the fuel. With regard to long-term SNF disposal and transportation, tremendous benefits, such as increased capacity, flexibility of design and system operations, and reduced overall costs, provide an incentive to seek burnup credit for criticality safety evaluations. The Nuclear Regulatory Commission issued Interim Staff Guidance 8, Revision 2 in 2002, endorsing burnup credit of actinide composition changes only; credit due to actinides encompasses approximately 30% of exiting pressurized water reactor SNF inventory and could potentially be increased to 90% if fission product credit were accepted. However, one significant issue for utilizing full burnup credit, compensating for actinide and fission product composition changes, is establishing a set of depletion parameters that produce an adequately conservative representation of the fuel's isotopic inventory. Depletion parameters can have a significant effect on the isotopic inventory of the fuel, and thus the residual reactivity. This research seeks to quantify the reactivity impact on a system from dominant depletion parameters (i.e., fuel temperature, moderator density, burnable poison rod, burnable poison rod history, and soluble boron concentration). Bounding depletion parameters were developed by statistical evaluation of a database containing reactor operating histories. The database was generated from summary reports of commercial reactor criticality data. Through depletion calculations, utilizing the SCALE 6 code package, several light
Nuclide analysis in high burnup fuel samples irradiated in Vandellós 2
NASA Astrophysics Data System (ADS)
Zwicky, H. U.; Low, J.; Granfors, M.; Alejano, C.; Conde, J. M.; Casado, C.; Sabater, J.; Lloret, M.; Quecedo, M.; Gago, J. A.
2010-07-01
In the framework of a high burnup fuel demonstration programme, rods with an enrichment of 4.5% 235U were operated to a rod average burnup of about 70 MWd/kgU in the Spanish Vandellós 2 pressurised water reactor. The rods were sent to hot cells and used for different research projects. This paper describes the isotopic composition measurements performed on samples of those rods, and the analysis of the measurement results based on comparison against calculated values. The fraction and composition of fission gases released to the rod free volume was determined for two of the rods. About 8% of Kr and Xe produced by fission were released. From the isotopic composition of the gases, it could be concluded that the gases were not preferentially released from the peripheral part of the fuel column. Local burnup and isotopic content of gamma emitting nuclides were determined by quantitatively evaluating axial gamma scans of the full rods. Nine samples were cut at different axial levels from three of the rods and analysed in two campaigns. More than 50 isotopes of 16 different elements were assessed, most of them by Inductively Coupled Plasma Mass Spectrometry after separation with High Performance Liquid Chromatography. In general, these over 400 data points gave a consistent picture of the isotopic content of irradiated fuel as a function of burnup. Only in a few cases, the analysis provided unexpected results that seem to be wrong, in most cases due to unidentified reasons. Sample burnup analysis was performed by comparing experimental isotopic abundances of uranium and plutonium composition as well as neodymium isotopic concentrations with corresponding CASMO based data. The results were in agreement with values derived independently from gamma scanning and from core design data and plant operating records. Measured isotope abundances were finally assessed using the industry standard SAS2H sequence of the SCALE code system. This exercise showed good agreement between
Whole-Core Heterogeneous Transport Calculations and Their Comparison with Diffusion Results
Nam Zin Cho; Gil Soo Lee; Ser Gi Hong; Chang Keun Jo; Kyung Taek Lee
2000-11-12
Recently the method of characteristics (MOC) has been considered as an effective methodology in lattice calculations. This method gives accurate solutions in complex geometries and strong absorber problems. With increasingly more heterogeneous reactor cores such as a mixed-oxide (MOX) fuel-loaded core or a burnable absorber-loaded core, the limitations due to homogenization and diffusion theory are evident, and the need for whole-core heterogeneous transport calculations is becoming greater. The CRX code based on the MOC is extended to treat whole-core heterogeneous calculation. Since the heterogeneous transport calculation for such a large-scale problem requires large computer memory, a modular ray tracing in which all lattice cells have the same ray distribution for each direction was used to reduce the computer memory requirement. In this scheme, the ray tracing is performed only on different types of cells. Therefore, this ray tracing scheme can significantly reduce the time in tracing along neutron paths and the computer memory for storing track lengths. Also, a parallelization scheme in angular domain rather than in spatial domain and the coarse mesh/coarse group rebalance (CMR/CGR) method in inner and outer iterations were implemented for further reduction of the computer time. To show the effectiveness of the extended CRX code, it is applied to heterogeneous calculation of a benchmark problem core (i.e., 10 x 10 whole core). The results of the transport calculations by CRX are compared with those of TWODANT and with those of the diffusion nodal codes AFEN and NEM (nodal expansion method). Unless the usual homogenization based on single-assembly calculation is drastically improved, the nodal methods would have to be superseded by whole-core heterogeneous calculation methods. It would be feasible to perform whole-core heterogeneous transport calculations routinely if the MOC implemented in the CRX code is enhanced further by more effective acceleration
Inherent safety of minimum-burnup breed and burn reactors
Qvist, S.; Reenspan, E.
2012-07-01
Reactors that aim to sustain the breed and burn (B and B) mode of operation at minimum discharge burnup require excellent neutron economy, Minimum-burnup B and B cores are generally large and feature low neutron leakage probability and a hard neutron spectrum. While highly promising fuel cycles can be achieved with such designs, the very same features are pushing the limits of the core's ability to passively respond safely to unprotected accidents. Low leakage minimum-burnup sodium-cooled B and B cores have a large positive coolant void-worth and coolant temperature reactivity coefficient. In this study, the applicability of major approaches for fast reactor void-worth reduction is evaluated specifically for B and B cores. The design, shuffling scheme and performance of a new metallic-fueled, sodium-cooled minimum burnup B and B core, used as basis for the void-worth reduction analysis, is presented. The analysis shows that reactivity control systems based on passive {sup 6}Li injection during temperature excursions are the only option able to provide negative void-worth without significantly increasing the minimum burnup required for sustaining the B and B mode of operation. A new type of lithium expansion module (LEM) system was developed specifically for B and B cores and its effect on core performance is presented. (authors)
/sup 18/O as a core plus two valence neutrons: A three-body Faddeev calculation
Ueta, K.; Miyake, H.; Mizukami, A.
1983-01-01
The nucleus /sup 18/O is studied assuming a three-body model: two neutrons outside an inert core of /sup 16/O: and solving the Faddeev equations. The calculated spectrum is in good agreement with experiment.
NASA Astrophysics Data System (ADS)
Quinet, Pascal
2014-09-01
A detailed investigation of the atomic structure and radiative parameters involving the lowest states within the 6p4, 6p36d, 6p37s, 6p37p and 6p37d configurations of neutral polonium is reported in the present paper. Using different physical models based on the pseudo-relativistic Hartree-Fock approach, the influence of intravalence, core-valence and core-core electron correlation on the atomic parameters is discussed in detail. This work allowed us to fix the spectroscopic designation of some experimental level energy values and to provide for the first time a set of reliable oscillator strengths corresponding to 31 Po I spectral lines in the wavelength region from 175 to 987 nm.
Local Burn-Up Effects in the NBSR Fuel Element
Brown N. R.; Hanson A.; Diamond, D.
2013-01-31
This study addresses the over-prediction of local power when the burn-up distribution in each half-element of the NBSR is assumed to be uniform. A single-element model was utilized to quantify the impact of axial and plate-wise burn-up on the power distribution within the NBSR fuel elements for both high-enriched uranium (HEU) and low-enriched uranium (LEU) fuel. To validate this approach, key parameters in the single-element model were compared to parameters from an equilibrium core model, including neutron energy spectrum, power distribution, and integral U-235 vector. The power distribution changes significantly when incorporating local burn-up effects and has lower power peaking relative to the uniform burn-up case. In the uniform burn-up case, the axial relative power peaking is over-predicted by as much as 59% in the HEU single-element and 46% in the LEU single-element with uniform burn-up. In the uniform burn-up case, the plate-wise power peaking is over-predicted by as much as 23% in the HEU single-element and 18% in the LEU single-element. The degree of over-prediction increases as a function of burn-up cycle, with the greatest over-prediction at the end of Cycle 8. The thermal flux peak is always in the mid-plane gap; this causes the local cumulative burn-up near the mid-plane gap to be significantly higher than the fuel element average. Uniform burn-up distribution throughout a half-element also causes a bias in fuel element reactivity worth, due primarily to the neutronic importance of the fissile inventory in the mid-plane gap region.
NASA Astrophysics Data System (ADS)
Sambuu, Odmaa; Nanzad, Norov
2009-03-01
A 300 kW small medical reactor was designed to be used for boron neutron capture therapy (BNCT) at KAIST in 1996 [1]. In this paper, analysis for the core life cycle of the original design of the BNCT facility and modifications of the fuel assembly configuration and enrichment to get a proper life cycle were performed and a criticality, neutron flux distribution and fuel burnup calculations were carried out.
Sambuu, Odmaa; Nanzad, Norov
2009-03-31
A 300 kW small medical reactor was designed to be used for boron neutron capture therapy (BNCT) at KAIST in 1996. In this paper, analysis for the core life cycle of the original design of the BNCT facility and modifications of the fuel assembly configuration and enrichment to get a proper life cycle were performed and a criticality, neutron flux distribution and fuel burnup calculations were carried out.
NASA Astrophysics Data System (ADS)
Commerçon, B.; Hennebelle, P.; Levrier, F.; Launhardt, R.; Henning, Th.
2012-03-01
I will present radiation-magneto-hydrodynamics calculations of low-mass and massive dense core collapse, focusing on the first collapse and the first hydrostatic core (first Larson core) formation. The influence of magnetic field and initial mass on the fragmentation properties will be investigated. In the first part reporting low mass dense core collapse calculations, synthetic observations of spectral energy distributions will be derived, as well as classical observational quantities such as bolometric temperature and luminosity. I will show how the dust continuum can help to target first hydrostatic cores and to state about the nature of VeLLOs. Last, I will present synthetic ALMA observation predictions of first hydrostatic cores which may give an answer, if not definitive, to the fragmentation issue at the early Class 0 stage. In the second part, I will report the results of radiation-magneto-hydrodynamics calculations in the context of high mass star formation, using for the first time a self-consistent model for photon emission (i.e. via thermal emission and in radiative shocks) and with the high resolution necessary to resolve properly magnetic braking effects and radiative shocks on scales <100 AU (Commercon, Hennebelle & Henning ApJL 2011). In this study, we investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores (M=100 M_⊙). We identify a new mechanism that inhibits initial fragmentation of massive dense cores, where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation, while moderately magnetized massive dense cores are more appropriate to form OB associations or small star clusters. Finally we will also present synthetic observations of these
Core Calculation of 1 MWatt PUSPATI TRIGA Reactor (RTP) using Monte Carlo MVP Code System
Karim, Julia Abdul
2008-05-20
The Monte Carlo MVP code system was adopted for the Reaktor TRIGA PUSAPTI (RTP) core calculation. The code was developed by a group of researcher of Japan Atomic Energy Agency (JAEA) first in 1994. MVP is a general multi-purpose Monte Carlo code for neutron and photon transport calculation and able to estimate an accurate simulation problems. The code calculation is based on the continuous energy method. This code is capable of adopting an accurate physics model, geometry description and variance reduction technique faster than conventional method as compared to the conventional scalar method. This code could achieve higher computational speed by several factors on the vector super-computer. In this calculation, RTP core was modeled as close as possible to the real core and results of keff flux, fission densities and others were obtained.
J. W. Sterbentz
1999-08-01
Six uranium isotopes and fourteen fission product isotopes were calculated on a mass basis at end-of-life (EOL) conditions for three fuel rods from different Light Water Breeder Reactor (LWBR) measurements. The three fuel rods evaluated here were taken from an LWBR seed module, a standard blanket module, and a reflector (Type IV) module. The calculated results were derived using a depletion methodology previously employed to evaluate many of the radionuclide inventories for spent nuclear fuels at the Idaho National Engineering and Environmental Laboratory. The primary goal of the calculational task was to further support the validation of this particular calculational methodology and its application to diverse reactor types and fuels. Result comparisons between the calculated and measured mass concentrations in the three rods indicate good agreement for the three major uranium isotopes (U-233, U-234, U-235) with differences of less than 20%. For the seed and standard blanket rod, the U-233 and U-234 differences were within 5% of the measured values (these two isotopes alone represent greater than 97% of the EOL total uranium mass). For the major krypton and xenon fission product isotopes, differences of less than 20% and less than 30% were observed, respectively. In general, good agreement was obtained for nearly all the measured isotopes. For these isotopes exhibiting significant differences, possible explanations are discussed in terms of measurement uncertainty, complex transmutations, etc.
Interaction of loading pattern and nuclear data uncertainties in reactor core calculations
Klein, M.; Gallner, L.; Krzykacz-Hausmann, B.; Pautz, A.; Velkov, K.; Zwermann, W.
2012-07-01
Along with best-estimate calculations for design and safety analysis, understanding uncertainties is important to determine appropriate design margins. In this framework, nuclear data uncertainties and their propagation to full core calculations are a critical issue. To deal with this task, different error propagation techniques, deterministic and stochastic are currently developed to evaluate the uncertainties in the output quantities. Among these is the sampling based uncertainty and sensitivity software XSUSA which is able to quantify the influence of nuclear data covariance on reactor core calculations. In the present work, this software is used to investigate systematically the uncertainties in the power distributions of two PWR core loadings specified in the OECD UAM-Benchmark suite. With help of a statistical sensitivity analysis, the main contributors to the uncertainty are determined. Using this information a method is studied with which loading patterns of reactor cores can be optimized with regard to minimizing power distribution uncertainties. It is shown that this technique is able to halve the calculation uncertainties of a MOX/UOX core configuration. (authors)
Hartini, Entin Andiwijayakusuma, Dinan
2014-09-30
This research was carried out on the development of code for uncertainty analysis is based on a statistical approach for assessing the uncertainty input parameters. In the butn-up calculation of fuel, uncertainty analysis performed for input parameters fuel density, coolant density and fuel temperature. This calculation is performed during irradiation using Monte Carlo N-Particle Transport. The Uncertainty method based on the probabilities density function. Development code is made in python script to do coupling with MCNPX for criticality and burn-up calculations. Simulation is done by modeling the geometry of PWR terrace, with MCNPX on the power 54 MW with fuel type UO2 pellets. The calculation is done by using the data library continuous energy cross-sections ENDF / B-VI. MCNPX requires nuclear data in ACE format. Development of interfaces for obtaining nuclear data in the form of ACE format of ENDF through special process NJOY calculation to temperature changes in a certain range.
Accounting for strong localized heterogeneities and local transport effect in core calculations
Ruggieri, J.M.; Doriath, J.Y.; Finck, P.J.; Boyer, R.
1996-09-01
Two methods based on the variational nodal transport method have been developed to account for localized heterogeneities and local transport effects in full core calculations. A local mesh refinement technique relies on using the projected partial ingoing surface currents produced during coarse-mesh iterations as boundary conditions for fine-mesh calculations embedded within the coarse-mesh calculations. The outgoing fine-mesh partial currents are averaged to serve in the coarse-mesh iterations. Then, a mixed transport-diffusion method using two levels of angular approximations for the surface partial currents depending on the node considered has been implemented to account for local transport effects in full core diffusion calculations. These methods have been tested for a model of the Superphenix complementary shutdown rods.
Neutronic calculations for the conversion to LEU of a research reactor core
Varvayanni, M.; Catsaros, N.; Stakakis, E.; Grigoriadis, D.
2008-07-15
For a five-year transitional period the Greek Research Reactor (GRR-1) was operating with a mixed core, containing both Low Enrichment (LEU) and High Enrichment (HEU) Uranium MTR- type fuel assemblies. The neutronic study of the GRR-1 conversion to LEU has been performed using a code system comprising the core-analysis code CITATION-LDI2 and the cell-calculation modules XSDRNPM and NITAWL-II of the SCALE code. A conceptual LEU core configuration was defined and analyzed with respect to the three dimensional multi-group neutron fluxes, the power distribution, the control-rod worth and the compliance with pre-defined Operation Limiting Conditions. Perturbation calculations and reactivity feedback computations were also carried out to provide input to a subsequent thermal-hydraulic study. (author)
Electronic Structure Calculations and Adaptation Scheme in Multi-core Computing Environments
Seshagiri, Lakshminarasimhan; Sosonkina, Masha; Zhang, Zhao
2009-05-20
Multi-core processing environments have become the norm in the generic computing environment and are being considered for adding an extra dimension to the execution of any application. The T2 Niagara processor is a very unique environment where it consists of eight cores having a capability of running eight threads simultaneously in each of the cores. Applications like General Atomic and Molecular Electronic Structure (GAMESS), used for ab-initio molecular quantum chemistry calculations, can be good indicators of the performance of such machines and would be a guideline for both hardware designers and application programmers. In this paper we try to benchmark the GAMESS performance on a T2 Niagara processor for a couple of molecules. We also show the suitability of using a middleware based adaptation algorithm on GAMESS on such a multi-core environment.
NASA Technical Reports Server (NTRS)
Huo, Winifred M.; Kim, Yong-Ki
1999-01-01
Based on the Binary-Encounter-Bethe (BEB) model, the advantage of using relativistic effective core potentials (RECP) in the calculation of total ionization cross sections of heavy atoms or molecules containing heavy atoms is discussed. Numerical examples for Ar, Kr, Xe, and WF6 are presented.
High-resolution Valence and Core Excitation Spectra via First-Principles Calculations and Experiment
NASA Astrophysics Data System (ADS)
Shirley, Eric; Fossard, F.; Gilmore, K.; Hug, G.; Kas, J. J.; Rehr, J. J.; Vila, F.
We calculate the optical and C K-edge near edge spectra of crystalline and molecular C60 measured with high-resolution electron energy-loss spectroscopy. The calculations are carried out using at least three different methods: Bethe-Salpeter calculations using the NIST Bethe-Salpeter Equation solver (NBSE) in the valence and OCEAN (Obtaining Core Excitation with Ab initio methods and NBSE) suite [Gilmore et al., Comp. Phys. Comm., (2015)]; excited-core-hole calculations using XCH [D. Prendergast and G. Galli, Phys. Rev. Lett. 96, 215502 (2006)]; and constrained occupancy using StoBe (Stockholm-Berlin core-excitation code) [StoBe-deMon version 3.0, K. Hermann et al. (2009)]. They include self-energy effects, lifetime-damping, and Debye-Waller effects. A comparison of spectral features to those observed illustrates the sensitivity of certain features to computation details (e.g., self-energy corrections and core-hole screening). This may point to limitations of various approximations, e.g. in conventional BSE paradigm and/or the incomplete treatment of vibrational effects. Supported in part by DOE BES Grant DE-FG03-97ER45623 (JJR, JJK, FV).
Jung, Y. S.; Lee, U. C.; Joo, H. G.
2012-07-01
The traditional subgroup method which has been applied for lattice transport calculations has an inherent limitation for non-uniform temperature distributions. As a measure to incorporate temperature dependence into the subgroup formulation, the subgroup level and number density adjustment method have been proposed. In this paper, the temperature dependent subgroup formulations employed for reflecting the non-uniform temperature effects on the resonance spatial self-shielding are examined for the whole core transport calculation with the thermal feedback. For 2D pin-cell problem with non-uniform temperature profiles, the inherent limitation of conventional subgroup method is confirmed. And the improvement in terms of reactivity is observed with the proposed adjustment scheme. For the real PWR core calculation with thermal feedback in the hot-full-power condition, the noticeable correction for the fuel temperature coefficient by about 10% more negative is obtained with the correction schemes. (authors)
NASA Astrophysics Data System (ADS)
Fensin, Michael Lorne
Monte Carlo-linked depletion methods have gained recent interest due to the ability to more accurately model complex 3-dimesional geometries and better track the evolution of temporal nuclide inventory by simulating the actual physical process utilizing continuous energy coefficients. The integration of CINDER90 into the MCNPX Monte Carlo radiation transport code provides a high-fidelity completely self-contained Monte-Carlo-linked depletion capability in a well established, widely accepted Monte Carlo radiation transport code that is compatible with most nuclear criticality (KCODE) particle tracking features in MCNPX. MCNPX depletion tracks all necessary reaction rates and follows as many isotopes as cross section data permits in order to achieve a highly accurate temporal nuclide inventory solution. This work chronicles relevant nuclear history, surveys current methodologies of depletion theory, details the methodology in applied MCNPX and provides benchmark results for three independent OECD/NEA benchmarks. Relevant nuclear history, from the Oklo reactor two billion years ago to the current major United States nuclear fuel cycle development programs, is addressed in order to supply the motivation for the development of this technology. A survey of current reaction rate and temporal nuclide inventory techniques is then provided to offer justification for the depletion strategy applied within MCNPX. The MCNPX depletion strategy is then dissected and each code feature is detailed chronicling the methodology development from the original linking of MONTEBURNS and MCNP to the most recent public release of the integrated capability (MCNPX 2.6.F). Calculation results of the OECD/NEA Phase IB benchmark, H. B. Robinson benchmark and OECD/NEA Phase IVB are then provided. The acceptable results of these calculations offer sufficient confidence in the predictive capability of the MCNPX depletion method. This capability sets up a significant foundation, in a well established
Calculation of ex-core physical quantities using the 3D importance functions
NASA Astrophysics Data System (ADS)
Trakas, Christos; De Laubiere, Xavier
2014-06-01
Diverse physical quantities are calculated in engineering studies with penalizing hypotheses to assure the required operation margins for each reactor. Today, these physical quantities are obtained by direct calculations from deterministic or Monte Carlo codes. The related states are critical or sub-critical. The current physical quantities are for example: the SRD counting rates (source range detector) in the sub-critical state, the IRD (intermediary range detector) and PRD (power range detector) counting rates (neutron particles only), the deposited energy in the reflector (neutron + photon particles), the fluence or the DPA (displacement per atom) in the reactor vessel (neutron particles only). The reliability of the proposed methodology is tested in the EPR reactor. The main advantage of the new methodology is the simplicity to obtain the physical quantities by an easy matrix calculation importance linked to nuclear power sources for all the cycles of the reactor. This method also allows to by-pass the direct calculations of the physical quantity of irradiated cores by Monte Carlo Codes, these calculations being impossible today (too many isotopic concentrations / MCNP5 limit). This paper presents the first feasibility study for the physical quantities calculation outside of the core by the importance method instead of the direct calculations used currently by AREVA.
Fuel-Cycle of 'CANDLE' Burnup with Depleted Uranium
Hiroshi, Sekimoto
2006-07-01
A new reactor burnup strategy CANDLE (Constant Axial shape of Neutron flux, nuclide densities and power shape During Life of Energy producing reactor) was proposed, where shapes of neutron flux, nuclide densities and power density distributions remain constant but move upward (or downward) along its core axis. This burnup strategy can derive many merits, especially from safety point of view. The change of excess reactivity along burnup is theoretically zero, and the core characteristics, such as power feedback coefficients and power peaking factor, are not changed along burnup. Application of this burnup strategy to neutron rich fast reactors makes excellent performances. Only natural or depleted uranium is required for the replacing fuels. About 40 % of natural or depleted uranium undergoes fission without the conventional reprocessing and enrichment. If the LWR produced energy of X Joules, the CANDLE reactor can produce about 50 X Joules from the depleted uranium left at the enrichment facility for the LWR fuel. If we can say LWRs have produced energy sufficient for full 20 years, we can produce the energy for 1000 years by using the CANDLE reactors with depleted uranium. We need not mine any uranium ore, and do not need reprocessing facility. The burnup of spent fuel becomes 10 times. Therefore, the spent fuel amount per produced energy is also reduced to one-tenth. The equilibrium core contains a lot of instable materials such as higher actinides and fission products, the enough amounts of which can not be obtained easily. The construction of the initial core is a difficult problem. However, by using enriched uranium substituted for actinides in the equilibrium core, we can construct the initial core whose power profile is similar to the equilibrium one and will reach the equilibrium state without any big change during transient. At present we do not have any material standing for such a high burnup. However, the CANDLE burnup can be realized by employing
3D Neutron Transport PWR Full-core Calculation with RMC code
NASA Astrophysics Data System (ADS)
Qiu, Yishu; She, Ding; Fan, Xiao; Wang, Kan; Li, Zeguang; Liang, Jingang; Leroyer, Hadrien
2014-06-01
Nowadays, there are more and more interests in the use of Monte Carlo codes to calculate the detailed power density distributions in full-core reactors. With the Inspur TS1000 HPC Server of Tsinghua University, several calculations have been done based on the EDF 3D Neutron Transport PWR Full-core benchmark through large-scale parallelism. To investigate and compare the results of the deterministic method and Monte Carlo method, EDF R&D and Department of Engineering Physics of Tsinghua University are having a collaboration to make code to code verification. So in this paper, two codes are used. One is the code COCAGNE developed by the EDF R&D, a deterministic core code, and the other is the Monte Carlo code RMC developed by Department of Engineering Physics in Tsinghua University. First, the full-core model is described and a 26-group calculation was performed by these two codes using the same 26-group cross-section library provided by EDF R&D. Then the parallel and tally performance of RMC is discussed. RMC employs a novel algorithm which can cut down most of the communications. It can be seen clearly that the speedup ratio almost linearly increases with the nodes. Furthermore the cell-mapping method applied by RMC consumes little time to tally even millions of cells. The results of the codes COCAGNE and RMC are compared in three ways. The results of these two codes agree well with each other. It can be concluded that both COCAGNE and RMC are able to provide 3D-transport solutions associated with detailed power density distributions calculation in PWR full-core reactors. Finally, to investigate how many histories are needed to obtain a given standard deviation for a full 3D solution, the non-symmetrized condensed 2-group fluxes of RMC are discussed.
Lee, C.; Yang, W. S.
2013-07-01
An improved resonance self-shielding method has been developed to accurately estimate the effective multigroup cross sections for heterogeneous fast reactor assembly and core calculations. In the method, the heterogeneity effect is considered by the use of isotopic escape cross sections while the resonance interference effect is accounted for through the narrow resonance approximation or slowing-down calculations for specific compositions. The isotopic escape cross sections are calculated by solving fixed-source transport equations with the method of characteristics for the whole problem domain. This method requires no pre-calculated resonance integral tables or parameters that are typically necessary in the subgroup method. Preliminary results for multi pin-cell fast reactor problems show that the escape cross sections estimated from the explicit-geometry fixed source calculations produce more accurate eigenvalue and self-shielded effective cross sections than those from conventional one-dimensional geometry models. (authors)
Monte-Carlo Continuous Energy Burnup Code System.
Energy Science and Technology Software Center (ESTSC)
2007-08-31
Version 00 MCB is a Monte Carlo Continuous Energy Burnup Code for a general-purpose use to calculate a nuclide density time evolution with burnup or decay. It includes eigenvalue calculations of critical and subcritical systems as well as neutron transport calculations in fixed source mode or k-code mode to obtain reaction rates and energy deposition that are necessary for burnup calculations. The MCB-1C patch file and data packages as distributed by the NEADB are verymore » well organized and are being made available through RSICC as received. The RSICC package includes the MCB-1C patch and MCB data libraries. Installation of MCB requires MCNP4C source code and utility programs, which are not included in this MCB distribution. They were provided with the now obsolete CCC-700/MCNP-4C package.« less
Extended step characteristic model for quarter-core gamma heating calculations
DeHart, M.D.; Webb, R.L. )
1993-01-01
Discrete ordinates codes are seldom used in lattice or core calculation, because of their limitation to simple geometries, which can be represented using an orthogonal mesh in a given coordinate system. Rough geometric approximations are often applied to obtain an estimate for a flux distribution. However, other methods, such as integral transport or Monte Carlo approaches, are generally more suited to irregular geometries. Each of these methods has its own weaknesses: integral transport methods are limited to problems in which the angular variation of the flux is isotropic or linearly anisotropic; Monte Carlo methods can be time consuming. The extended step characteristic (ESC) method has been developed to apply the discrete ordinates approximation to complicated geometries for which other methods provide less satisfactory solutions. The CENTAUR code has been developed to solve the two-dimensional transport equation using the ESC approach. This paper presents results of CENTAUR calculations for a quarter-core gamma redistribution problem for the Savannah River site (SRS) K reactor, under drained tank conditions following a postulated double-ended guillotine break loss-of-coolant accident. The calculations were used to confirm TWOTRAN calculations, which were based on a coarse approximation of the core geometry. A comparison of the results serves to demonstrate the capabilities and efficiency of the ESC approach.
TORT/MCNP coupling method for the calculation of neutron flux around a core of BWR.
Kurosawa, Masahiko
2005-01-01
For the analysis of BWR neutronics performance, accurate data are required for neutron flux distribution over the In-Reactor Pressure Vessel equipments taking into account the detailed geometrical arrangement. The TORT code can calculate neutron flux around a core of BWR in a three-dimensional geometry model, but has difficulties in fine geometrical modelling and lacks huge computer resource. On the other hand, the MCNP code enables the calculation of the neutron flux with a detailed geometry model, but requires very long sampling time to give enough number of particles. Therefore, a TORT/MCNP coupling method has been developed to eliminate the two problems mentioned above in each code. In this method, the TORT code calculates angular flux distribution on the core surface and the MCNP code calculates neutron spectrum at the points of interest using the flux distribution. The coupling method will be used as the DOT-DOMINO-MORSE code system. This TORT/MCNP coupling method was applied to calculate the neutron flux at points where induced radioactivity data were measured for 54Mn and 60Co and the radioactivity calculations based on the neutron flux obtained from the above method were compared with the measured data. PMID:16604689
Liquid iron-sulfur alloys at outer core conditions by first-principles calculations
NASA Astrophysics Data System (ADS)
Umemoto, Koichiro; Hirose, Kei; Imada, Saori; Nakajima, Yoichi; Komabayashi, Tetsuya; Tsutsui, Satoshi; Baron, Alfred Q. R.
2014-10-01
We perform first-principles calculations to investigate liquid iron-sulfur alloys (Fe, Fe56S8, Fe52S12, and Fe48S16) under high-pressure and high-temperature (150-300 GPa and 4000-6000 K) conditions corresponding to the Earth's outer core. Considering only the density profile, the best match with the preliminary reference Earth model is by liquid Fe-14 wt % S (Fe50S14), assuming sulfur is the only light element. However, its bulk sound velocity is too high, in particular in the deep outer core, suggesting that another light component such as oxygen is required. An experimental check using inelastic X-ray scattering shows good agreement with the calculations. In addition, a present study demonstrates that the Birch's law does not hold for liquid iron-sulfur alloy, consistent with a previous report on pure liquid iron.
Piezoelectric constants for ZnO calculated using classical polarizable core-shell potentials
NASA Astrophysics Data System (ADS)
Dai, Shuangxing; Dunn, Martin L.; Park, Harold S.
2010-11-01
We demonstrate the feasibility of using classical atomistic simulations, i.e. molecular dynamics and molecular statics, to study the piezoelectric properties of ZnO using core-shell interatomic potentials. We accomplish this by reporting the piezoelectric constants for ZnO as calculated using two different classical interatomic core-shell potentials: that originally proposed by Binks and Grimes (1994 Solid State Commun. 89 921-4), and that proposed by Nyberg et al (1996 J. Phys. Chem. 100 9054-63). We demonstrate that the classical core-shell potentials are able to qualitatively reproduce the piezoelectric constants as compared to benchmark ab initio calculations. We further demonstrate that while the presence of the shell is required to capture the electron polarization effects that control the clamped ion part of the piezoelectric constant, the major shortcoming of the classical potentials is a significant underprediction of the clamped ion term as compared to previous ab initio results. However, the present results suggest that overall, these classical core-shell potentials are sufficiently accurate to be utilized for large scale atomistic simulations of the piezoelectric response of ZnO nanostructures.
NASA Astrophysics Data System (ADS)
Antsiferova, E. V.; Bogdanov, V. V.; Derebenko, E. V.; Lagutina, A. V.; Khmelnikov, E. A.
2006-08-01
The up-to-date development of the armored vehicles conditions complication of armor constructions and increased slope of shell armored plates. Combined strikers (C/S) can be used to destroy armored vehicles. We can increase total weight of the core part to increase the striker's power. However, the increase of core part diameter is limited by body dimensions. Thus, we can increase core part weight by increasing its length. Because of C/S interaction with the barriers at large deviation angles, C/S's mechanical trajectory sparks in the barrier. This results in bending stress which occurs in the core part. Because of large deviation angles, the impact of the side surface of oblong core part against the cavity edge occurs. This increases the probability of core part destruction. The calculation technique for oblong core part penetration into different types of barriers is presented. The large number of factors can be calculated using this technique. It is assumed that the core part is destroyed when the tail part impacts against the cavity in the section where specific impact energy exceeds the critical value. Impact elasticity and destruction at bending stress were selected to be destruction criteria. The following core part destruction scenarios were investigated and calculated: (i) core head part is slightly destroyed but tail part of cylindrical shape penetrates deeper; (ii) core tail part is slightly destroyed but head part penetrates deeper, mass loss is taken into account; and (iii) after the impact, the core part is splitted up into two parts, then both of them penetrate into the barrier, one part is of ogival shape, the other is of cylindrical one. This calculation technique was applied to computational program, then critical angles at which core part side surface is still in contact with cavity surface, and the angles at which core part destruction occurs were calculated. Depths of core part penetration for different destruction scenarios were calculated.
Efficient implementation of core-excitation Bethe-Salpeter equation calculations
NASA Astrophysics Data System (ADS)
Gilmore, K.; Vinson, John; Shirley, E. L.; Prendergast, D.; Pemmaraju, C. D.; Kas, J. J.; Vila, F. D.; Rehr, J. J.
2015-12-01
We present an efficient implementation of the Bethe-Salpeter equation (BSE) method for obtaining core-level spectra including X-ray absorption (XAS), X-ray emission (XES), and both resonant and non-resonant inelastic X-ray scattering spectra (N/RIXS). Calculations are based on density functional theory (DFT) electronic structures generated either by ABINIT or QuantumESPRESSO, both plane-wave basis, pseudopotential codes. This electronic structure is improved through the inclusion of a GW self energy. The projector augmented wave technique is used to evaluate transition matrix elements between core-level and band states. Final two-particle scattering states are obtained with the NIST core-level BSE solver (NBSE). We have previously reported this implementation, which we refer to as OCEAN (Obtaining Core Excitations from Ab initio electronic structure and NBSE) (Vinson et al., 2011). Here, we present additional efficiencies that enable us to evaluate spectra for systems ten times larger than previously possible; containing up to a few thousand electrons. These improvements include the implementation of optimal basis functions that reduce the cost of the initial DFT calculations, more complete parallelization of the screening calculation and of the action of the BSE Hamiltonian, and various memory reductions. Scaling is demonstrated on supercells of SrTiO3 and example spectra for the organic light emitting molecule Tris-(8-hydroxyquinoline)aluminum (Alq3) are presented. The ability to perform large-scale spectral calculations is particularly advantageous for investigating dilute or non-periodic systems such as doped materials, amorphous systems, or complex nano-structures.
Calculation of the reactivity feedback due to core-assembly bowing in LMFBRs
Not Available
1983-01-01
The nonuniformity of the temperature distribution in an LMFBR leads to differential thermal expansion of the walls of an assembly hexcan. These thermal expansion differentials cause the hexcan to distort or bow. Consequentially, the assembly experiences a spatial displacement, which results in a change in reactivity for the core. A computational model to calculate the reactivity feedback due to material displacements induced by assembly bowing effects has been developed.
Vitruk, S.G.; Korsun, A.S.; Ushakov, P.A.
1995-09-01
The multilevel mathematical model of neutron thermal hydrodynamic processes in a passive safety core without assemblies duct walls and appropriate computer code SKETCH, consisted of thermal hydrodynamic module THEHYCO-3DT and neutron one, are described. A new effective discretization technique for energy, momentum and mass conservation equations is applied in hexagonal - z geometry. The model adequacy and applicability are presented. The results of the calculations show that the model and the computer code could be used in conceptual design of advanced reactors.
Designing Critical Experiments in Support of Full Burnup Credit
Mueller, Don; Roberts, Jeremy A
2008-01-01
Burnup credit is the process of accounting for the negative reactivity due to fuel burnup and generation of parasitic absorbers over fuel assembly lifetime. For years, the fresh fuel assumption was used as a simple bound in criticality work for used fuel storage and transportation. More recently, major actinides have been included [1]. However, even this yields a highly conservative estimate in criticality calculations. Because of the numerous economical benefits including all available negative reactivity (i.e., full burnup credit) could provide [2], it is advantageous to work toward full burnup credit. Unfortunately, comparatively little work has been done to include non-major actinides and other fission products (FP) in burnup credit analyses due in part to insufficient experimental data for validation of codes and nuclear data. The Burnup Credit Criticality Experiment (BUCCX) at Sandia National Laboratory was a set of experiments with {sup 103}Rh that have relevance for burnup credit [3]. This work uses TSUNAMI-3D to investigate and adjust a BUCCX model to match isotope-specific, energy-dependent k{sub eff} sensitivity profiles to those of a representative high-capacity cask model (GBC-32) [4] for each FP of interest. The isotopes considered are {sup 149}Sm, {sup 143}Nd, {sup 103}Rh, {sup 133}Cs, {sup 155}Gd, {sup 152}Sm, {sup 99}Tc, {sup 145}Nd, {sup 153}Eu, {sup 147}Sm, {sup 109}Ag, {sup 95}Mo, {sup 150}Sm, {sup 101}Ru, and {sup 151}Eu. The goal is to understand the biases and bias uncertainties inherent in nuclear data, and ultimately, to apply these in support of full burnup credit.
Issues related to criticality safety analysis for burnup credit applications
DeHart, M.D.; Parks, C.V.
1995-12-01
Spent fuel transportation and storage cask designs based on a burnup credit approach must consider issues that are not relevant in casks designed under a fresh fuel loading assumption. Parametric analyses are required to characterize the importance of fuel assembly and fuel cycle parameters on spent fuel composition and reactivity. Numerical models are evaluated to determine the sensitivity of criticality safety calculations to modeling assumptions. This paper discusses the results of studies to determine the effect of two important modeling assumptions on the criticality analysis of pressurized-water reactor (PWR) spent fuel: (1) the effect of assumed burnup history (i.e., specific power during and time-dependent variations in operational power) during depletion calculations, and (2) the effect of axial burnup distributions on the neutron multiplication factor calculated for a three-dimensional (3-D) conceptual cask design.
Yamamoto, T.; Suzuki, M.; Ando, Y.
2012-07-01
After the severe core damage of Fukushima Dai-Ichi Nuclear Power Station, radioactive material leaked from the reactor buildings. As part of monitoring of radioactivity in the site, measurements of radioactivity in soils at three fixed points have been performed for {sup 134}Cs and {sup 137}Cs with gamma-ray spectrometry and for Pu, Pu, and {sup 240}Pu with {alpha}-ray spectrometry. Correlations of radioactivity ratios of {sup 134}Cs to {sup 137}Cs, and {sup 238}Pu to the sum of {sup 239}Pu and {sup 240}Pu with fuel burnup were studied by using theoretical burnup calculations and measurements on isotopic inventories, and compared with the Cs and Pu radioactivity rations in the soils. The comparison indicated that the burnup of the fuel sourcing the radioactivity was from 18 to 38 GWd/t, which corresponded to that of the fuel in the highest power and, therefore, the highest decay heat in operating high-burnup fueled BWR cores. (authors)
NASA Astrophysics Data System (ADS)
Sboev, A. G.; Ilyashenko, A. S.; Vetrova, O. A.
1997-02-01
The method of bucking evaluation, realized in the MOnte Carlo code MCS, is described. This method was applied for calculational analysis of well known light water experiments TRX-1 and TRX-2. The analysis of this comparison shows, that there is no coincidence between Monte Carlo calculations, obtained by different ways: the MCS calculations with given experimental bucklings; the MCS calculations with given bucklings evaluated on base of full core MCS direct simulations; the full core MCNP and MCS direct simulations; the MCNP and MCS calculations, where the results of cell calculations are corrected by the coefficients taking into the account the leakage from the core. Also the buckling values evaluated by full core MCS calculations have differed from experimental ones, especially in the case of TRX-1, when this difference has corresponded to 0.5 percent increase of Keff value.
Ab initio no core configuration interaction calculations in the natural orbital basis
NASA Astrophysics Data System (ADS)
Constantinou, Chrysovalantis; Caprio, Mark A.; Vary, James P.; Maris, Pieter
2015-10-01
The natural orbital basis has been successfully used in the past in atomic and molecular structure calculations. The natural orbitals used in those calculations are calculated by diagonalizing the electron one-body density matrix. Here we develop natural orbitals for nuclear no-core configuration interaction (NCCI) calculations. A NCCI calculation using an initial single particle basis, such as the harmonic oscillator basis, must first be performed in order to obtain a one-body density matrix. The eigenvectors of the one-body density matrix are the natural orbitals, and the corresponding eigenvalues are the occupations of these natural orbitals in the nuclear wave function. According to these occupancies, the most important natural orbitals, in the sense of the most occupied, can then be selected and used in a NCCI calculation. We discuss ab initio nuclear NCCI calculations for light nuclei and assess their ability to provide faster convergence. Supported by the US DOE (under Grants DE-FG02-95ER-40934, DESC0008485 SciDAC/NUCLEI, and DE-FG02-87ER40371), and the US NSF (under Grant 0904782). Computational resources provided by NERSC (supported by US DOE Contract DE-AC02-05CH11231), and NDCRC.
Calculation of Design Parameters for an Equilibrium LEU Core in the NBSR
Hanson, A.L.; Diamond, D.
2011-09-30
A plan is being developed for the conversion of the NIST research reactor (NBSR) from high-enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. Previously, the design of the LEU fuel had been determined in order to provide the users of the NBSR with the same cycle length as exists for the current HEU fueled reactor. The fuel composition at different points within an equilibrium fuel cycle had also been determined. In the present study, neutronics parameters have been calculated for these times in the fuel cycle for both the existing HEU and the proposed LEU equilibrium cores. The results showed differences between the HEU and LEU cores that would not lead to any significant changes in the safety analysis for the converted core. In general the changes were reasonable except that the figure-of-merit for neutrons that can be used by experimentalists shows there will be a 10% reduction in performance. The calculations included kinetics parameters, reactivity coefficients, reactivity worths of control elements and abnormal configurations, and power distributions.
Liquid iron-hydrogen alloys at outer core conditions by first-principles calculations
NASA Astrophysics Data System (ADS)
Umemoto, Koichiro; Hirose, Kei
2015-09-01
We examined the density, bulk sound (compressional) velocity, and Grüneisen parameter of liquid pure Fe, Fe100H28 (0.50 wt % H), Fe88H40 (0.81 wt % H), and Fe76H52 (1.22 wt % H) at Earth's outer core pressure and temperature (P-T) conditions (~100 to 350 GPa, 4000 to 7000 K) based on first-principles molecular dynamics calculations. The results demonstrate that the thermodynamic Grüneisen parameter of liquid iron alloy decreases with increasing pressure, temperature, and hydrogen concentration, indicating a relatively small temperature gradient in the outer core when hydrogen is present. Along such temperature profile, both the density and compressional velocity of liquid iron containing ~1 wt % hydrogen match seismological observations. It suggests that hydrogen could be a primary light element in the core, although the shear velocity of the inner core is not reconciled with solid Fe-H alloy and thus requires another impurity element.
Full Core 3-D Simulation of a Partial MOX LWR Core
S. Bays; W. Skerjanc; M. Pope
2009-05-01
A comparative analysis and comparison of results obtained between 2-D lattice calculations and 3-D full core nodal calculations, in the frame of MOX fuel design, was conducted. This study revealed a set of advantages and disadvantages, with respect to each method, which can be used to guide the level of accuracy desired for future fuel and fuel cycle calculations. For the purpose of isotopic generation for fuel cycle analyses, the approach of using a 2-D lattice code (i.e., fuel assembly in infinite lattice) gave reasonable predictions of uranium and plutonium isotope concentrations at the predicted 3-D core simulation batch average discharge burnup. However, it was found that the 2-D lattice calculation can under-predict the power of pins located along a shared edge between MOX and UO2 by as much as 20%. In this analysis, this error did not occur in the peak pin. However, this was a coincidence and does not rule out the possibility that the peak pin could occur in a lattice position with high calculation uncertainty in future un-optimized studies. Another important consideration in realistic fuel design is the prediction of the peak axial burnup and neutron fluence. The use of 3-D core simulation gave peak burnup conditions, at the pellet level, to be approximately 1.4 times greater than what can be predicted using back-of-the-envelope assumptions of average specific power and irradiation time.
Calculation of scattering characteristic of complex target on multi-core platform
NASA Astrophysics Data System (ADS)
Guo, Xing; Wu, Zhensen; Linghu, Longxiang
2013-09-01
The scattering characteristic of complex target from terrestrial and celestial background radiation has been widely used in such engineering fields as remote sensing, feature extraction, tracking and recognition of target thus having been an attractive field for many scientists for decades. In our method, the model of target is constructed using 3DMAX and the surface is divided into triangle facets firstly. Bidirectional Reflectance Distribution Function (BRDF) is introduced and MODTRAN is applied to calculate background radiation for a given time at a given place. Finally the scattering of each facet is added up to get the scattering of the target. As the background radiance comes in all directions and in a wide spectrum and the complex target always consists of thousands of facets, in general it takes hours to complete the calculation. Consequently this limits its use in the real time applications. Recent years have seen the continual development of multi-core CPU. As a result parallel programming on multi-cores has been more and more popular. In this paper, the openMP, Intel CILK ++, Intel Threading Building Blocks (TBB) are used separately to leverage the processing power of multi-cores processors. Our experiments are conducted on a DELL desktop based on an Intel I7- 2600K CPU running at 3.40 GHz with 8 cores and 16.0 GB RAM. The Intel Composer 2013 is employed to build the program. Also in OpenMP implementation, gcc is used. The results demonstrate that highest speedups for three parallel models are 5.06X, 5.02X, 5.15X respectively.
No Core CI calculations for light nuclei with chiral 2- and 3-body forces
NASA Astrophysics Data System (ADS)
Maris, Pieter; Metin Aktulga, H.; Binder, Sven; Calci, Angelo; Çatalyürek, Ümit V.; Langhammer, Joachim; Ng, Esmond; Saule, Erik; Roth, Robert; Vary, James P.; Yang, Chao
2013-08-01
The atomic nucleus is a self-bound system of strongly interacting nucleons. In No-Core Configuration Interaction calculations, the nuclear wavefunction is expanded in Slater determinants of single-nucleon wavefunctions (Configurations), and the many-body Schrödinger equation becomes a large sparse matrix problem. The challenge is to reach numerical convergence to within quantified numerical uncertainties for physical observables using finite truncations of the infinite-dimensional basis space. We discuss strategies for constructing and solving the resulting large sparse matrices for a set of low-lying eigenvalues and eigenvectors on current multicore computer architectures. Several of these strategies have been implemented in the code MFDn, a hybrid MPI/OpenMP Fortran code for ab initio nuclear structure calculations that scales well to over 200,000 cores. We discuss how the similarity renormalization group can be used to improve the numerical convergence. We present results for excitation energies and other selected observables for 8Be and 12C using realistic 2- and 3-body forces obtained from chiral perturbation theory. Finally, we demonstrate that collective phenomena such as rotational band structures can emerge from these microscopic calculations.
High Burnup Fuel Behavior Modeling
Jahingir, M.; Rand, R.; Stachowski, R.; Miles, B.; Kusagaya, K.
2007-07-01
This paper discusses the development and qualification of the PRIME03 code to address high burnup mechanisms and to improve uranium utilization in current and new reactor designs. Materials properties and behavioral models have been updated from previous thermal-mechanical codes to reflect the effects of burnup on fuel pellet thermal conductivity, Zircaloy creep, fuel pellet relocation, and fission gas release. These new models are based on results of in-pool and post irradiation examination (PIE) of commercial boiling water reactor (BWR) fuel rods at high burnup and results from international experimental programs. The new models incorporated into PRIME03 also address specific high burnup effects associated with formation of pellet rim porosity at high exposure. The PRIME03 code is qualified by comparison of predicted and measured fuel performance parameters for a large number of high, low, and moderate burnup test and commercial reactor rod. The extensive experimental qualification of the PRIME03 prediction capabilities confirms that it is a reliable best-estimate predictor of fuel rod thermal-mechanical performance over a wide range of design and operating conditions. (authors)
Statistical error propagation in ab initio no-core full configuration calculations of light nuclei
NASA Astrophysics Data System (ADS)
Navarro Pérez, R.; Amaro, J. E.; Ruiz Arriola, E.; Maris, P.; Vary, J. P.
2015-12-01
We propagate the statistical uncertainty of experimental N N scattering data into the binding energy of 3H and 4He. We also study the sensitivity of the magnetic moment and proton radius of the 3H to changes in the N N interaction. The calculations are made with the no-core full configuration method in a sufficiently large harmonic oscillator basis. For those light nuclei we obtain Δ Estat(3H) =0.015 MeV and Δ Estat(4He) =0.055 MeV .
Effect of Control Blade History, and Axial Coolant Density and Burnup Profiles on BWR Burnup Credit
Ade, Brian J; Marshall, William BJ J; Martinez-Gonzalez, Jesus S
2015-05-01
Oak Ridge National Laboratory (ORNL) and the US Nuclear Regulatory Commission (NRC) have initiated a multiyear project to investigate the application of burnup credit (BUC) for boiling water reactor (BWR) fuel in storage and transportation systems (often referred to as casks) and spent fuel pools (SFPs). This work is divided into two main phases. The first phase investigated the applicability of peak reactivity methods currently used in SFPs to transportation and storage casks and the validation of reactivity calculations and spent fuel compositions within these methods. The second phase focuses on extending BUC beyond peak reactivity. This paper documents the analysis of the effects of control blade insertion history, and moderator density and burnup axial profiles for extended BWR BUC.
NASA Astrophysics Data System (ADS)
Caprio, Mark A.; Maris, Pieter; Vary, James P.
2014-03-01
The emergence of rotational bands has recently been observed in no-core configuration interaction (NCCI) calculations for p-shell nuclei, as evidenced by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. Yrast and low-lying excited bands are found. The results demonstrate the possibility of well-developed rotational structure in NCCI calculations, using realistic nucleon-nucleon interactions, and within finite, computationally-accessible configuration spaces. This talk will focus on results for rotation in both the even-mass and odd-mass Be isotopes (7 <= A <= 12). Supported by US DOE (DE-FG02-95ER-40934, DESC0008485 SciDAC/NUCLEI, DE-FG02-87ER40371), US NSF (0904782), and Research Corporation for Science Advancement (Cottrell Scholar Award). Computational resources provided by NERSC (US DOE DE-AC02-05CH11231).
Whole-core neutron transport calculations without fuel-coolant homogenization
Smith, M. A.; Tsoulfanidis, N.; Lewis, E. E.; Palmiotti, G.; Taiwo, T. A.
2000-02-10
The variational nodal method implemented in the VARIANT code is generalized to perform full core transport calculations without spatial homogenization of cross sections at either the fuel-pin cell or fuel assembly level. The node size is chosen to correspond to one fuel-pin cell in the radial plane. Each node is divided into triangular finite subelements, with the interior spatial flux distribution represented by piecewise linear trial functions. The step change in the cross sections at the fuel-coolant interface can thus be represented explicitly in global calculations while retaining the fill spherical harmonics capability of VARIANT. The resulting method is applied to a two-dimensional seven-group representation of a LWR containing MOX fuel assemblies. Comparisons are made of the accuracy of various space-angle approximations and of the corresponding CPU times.
Energy-Deposition and Damage Calculations in Core-Vessel Inserts at the Spallation Neutron Source
Murphy, B.D.
2002-06-25
Heat-deposition and damage calculations are described for core-vessel inserts in the target area of the Spallation Neutron Source. Two separate designs for these inserts (or neutron beam tubes) were studied; a single-unit insert and a multi-unit insert. The single unit contains a neutron guide; the multi unit does not. Both units are constructed of stainless steel. For the single unit, separate studies were carried out with the guide composed of stainless steel, glass, and aluminum. Results are also reported for an aluminum window on the front of the insert, a layer of nickel on the guide, a cadmium shield surrounding the guide, and a stainless steel plug in the beam-tube opening. The locations of both inserts were the most forward positions to be occupied by each design respectively thus ensuring that the calculations are conservative.
Strategies for Application of Isotopic Uncertainties in Burnup Credit
Gauld, I.C.
2002-12-23
Uncertainties in the predicted isotopic concentrations in spent nuclear fuel represent one of the largest sources of overall uncertainty in criticality calculations that use burnup credit. The methods used to propagate the uncertainties in the calculated nuclide concentrations to the uncertainty in the predicted neutron multiplication factor (k{sub eff}) of the system can have a significant effect on the uncertainty in the safety margin in criticality calculations and ultimately affect the potential capacity of spent fuel transport and storage casks employing burnup credit. Methods that can provide a more accurate and realistic estimate of the uncertainty may enable increased spent fuel cask capacity and fewer casks needing to be transported, thereby reducing regulatory burden on licensee while maintaining safety for transporting spent fuel. This report surveys several different best-estimate strategies for considering the effects of nuclide uncertainties in burnup-credit analyses. The potential benefits of these strategies are illustrated for a prototypical burnup-credit cask design. The subcritical margin estimated using best-estimate methods is discussed in comparison to the margin estimated using conventional bounding methods of uncertainty propagation. To quantify the comparison, each of the strategies for estimating uncertainty has been performed using a common database of spent fuel isotopic assay measurements for pressurized-light-water reactor fuels and predicted nuclide concentrations obtained using the current version of the SCALE code system. The experimental database applied in this study has been significantly expanded to include new high-enrichment and high-burnup spent fuel assay data recently published for a wide range of important burnup-credit actinides and fission products. Expanded rare earth fission-product measurements performed at the Khlopin Radium Institute in Russia that contain the only known publicly-available measurement for {sup 103
Hybrid parallel code acceleration methods in full-core reactor physics calculations
Courau, T.; Plagne, L.; Ponicot, A.; Sjoden, G.
2012-07-01
When dealing with nuclear reactor calculation schemes, the need for three dimensional (3D) transport-based reference solutions is essential for both validation and optimization purposes. Considering a benchmark problem, this work investigates the potential of discrete ordinates (Sn) transport methods applied to 3D pressurized water reactor (PWR) full-core calculations. First, the benchmark problem is described. It involves a pin-by-pin description of a 3D PWR first core, and uses a 8-group cross-section library prepared with the DRAGON cell code. Then, a convergence analysis is performed using the PENTRAN parallel Sn Cartesian code. It discusses the spatial refinement and the associated angular quadrature required to properly describe the problem physics. It also shows that initializing the Sn solution with the EDF SPN solver COCAGNE reduces the number of iterations required to converge by nearly a factor of 6. Using a best estimate model, PENTRAN results are then compared to multigroup Monte Carlo results obtained with the MCNP5 code. Good consistency is observed between the two methods (Sn and Monte Carlo), with discrepancies that are less than 25 pcm for the k{sub eff}, and less than 2.1% and 1.6% for the flux at the pin-cell level and for the pin-power distribution, respectively. (authors)
NASA Technical Reports Server (NTRS)
Tuma, Margaret L.; Weisshaar, Andreas; Li, Jian; Beheim, Glenn
1995-01-01
To determine the feasibility of coupling the output of a single-mode optical fiber into a single-mode rib waveguide in a temperature varying environment, a theoretical calculation of the coupling efficiency between the two was investigated. Due to the complex geometry of the rib guide, there is no analytical solution to the wave equation for the guided modes, thus, approximation and/or numerical techniques must be utilized to determine the field patterns of the guide. In this study, three solution methods were used for both the fiber and guide fields; the effective-index method (EIM), Marcatili's approximation, and a Fourier method. These methods were utilized independently to calculate the electric field profile of each component at two temperatures, 20 C and 300 C, representing a nominal and high temperature. Using the electric field profile calculated from each method, the theoretical coupling efficiency between an elliptical-core optical fiber and a rib waveguide was calculated using the overlap integral and the results were compared. It was determined that a high coupling efficiency can be achieved when the two components are aligned. The coupling efficiency was more sensitive to alignment offsets in the y direction than the x, due to the elliptical modal field profile of both components. Changes in the coupling efficiency over temperature were found to be minimal.
NASA Astrophysics Data System (ADS)
Espel, Federico Puente
The main objective of this PhD research is to develop a high accuracy modeling tool using a Monte Carlo based coupled system. The presented research comprises the development of models to include the thermal-hydraulic feedback to the Monte Carlo method and speed-up mechanisms to accelerate the Monte Carlo criticality calculation. Presently, deterministic codes based on the diffusion approximation of the Boltzmann transport equation, coupled with channel-based (or sub-channel based) thermal-hydraulic codes, carry out the three-dimensional (3-D) reactor core calculations of the Light Water Reactors (LWRs). These deterministic codes utilize nuclear homogenized data (normally over large spatial zones, consisting of fuel assembly or parts of fuel assembly, and in the best case, over small spatial zones, consisting of pin cell), which is functionalized in terms of thermal-hydraulic feedback parameters (in the form of off-line pre-generated cross-section libraries). High accuracy modeling is required for advanced nuclear reactor core designs that present increased geometry complexity and material heterogeneity. Such high-fidelity methods take advantage of the recent progress in computation technology and coupled neutron transport solutions with thermal-hydraulic feedback models on pin or even on sub-pin level (in terms of spatial scale). The continuous energy Monte Carlo method is well suited for solving such core environments with the detailed representation of the complicated 3-D problem. The major advantages of the Monte Carlo method over the deterministic methods are the continuous energy treatment and the exact 3-D geometry modeling. However, the Monte Carlo method involves vast computational time. The interest in Monte Carlo methods has increased thanks to the improvements of the capabilities of high performance computers. Coupled Monte-Carlo calculations can serve as reference solutions for verifying high-fidelity coupled deterministic neutron transport methods
Recent development of VTT`s calculation system for VVERs
Kaloinen, E.; Antila, M.; Kyrki-Rajamaeki, R.
1994-12-31
A comprehensive code system for reactor physics and dynamics calculations in hexagonal core geometry has been developed at VTT ENERGY. The main modules of the system consist of the CASMO-HEX assembly burnup code and the three-dimensional codes HEXBU-3D for burnup simulation of VVER cores and HEXTRAN for transient and accident analysis of the core and coolant circuit. Figure 1 presents a schematic view of the calculational steps from a nuclear data library to static and dynamic simulations of the reactor. The code system has recently undergone several improvements and modifications, particularly to match the modern requirements for safety analyses. This paper gives a survey of the code system in its present state.
Thermal-hydraulic calculations for the conversion to LEU of a research reactor core
Grigoriadis, D.; Varvayanni, M.; Catsaros, N.; Stakakis, E.
2008-07-15
The thermal-hydraulic analysis performed for the needs of the conversion of the open pool 5MW Greek Research Reactor (GRR-1) to a pure Low Enrichment (LEU) configuration is presented. The methodology was based on a complete set of neutronic calculations performed for the new core configuration, in compliance with pre-defined Operation Limiting Conditions. The hottest channel analysis approach was adopted, and peaking factors were used to account for fabrication or measuring uncertainties. Calculations were carried out using the numerical codes NATCON, PLTEMP and PARET provided by Argonne National Laboratory (ANL). Two main different classes of conditions were considered, namely i) steady state normal operating conditions and ii) transient cases related to accidental events including reactivity feedback effects. For steady state operating conditions the behaviour of the new configuration was examined both for forced and natural convection cooling modes. Transient calculations considered several initiating events including reactivity insertion accidents (slow or fast reactivity insertion) and total or partial loss-of-flow accidents, i.e. in accordance to guidelines provided by the IAEA for research Reactors. (author)
Santra, Robin; Christ, Kevin V.; Greene, Chris H.
2004-04-01
The first three electronically excited states in the alkaline-earth-metal atoms magnesium, calcium, and strontium comprise the (nsnp){sup 3}P{sub J}{sup o}(J=0,1,2) fine-structure manifold. All three states are metastable and are of interest for optical atomic clocks as well as for cold-collision physics. An efficient technique--based on a physically motivated potential that models the presence of the ionic core--is employed to solve the Schroedinger equation for the two-electron valence shell. In this way, radiative lifetimes, laser-induced clock shifts, and long-range interaction parameters are calculated for metastable Mg, Ca, and Sr.
NASA Astrophysics Data System (ADS)
Li, S.-Y.; Niklasson, G. A.; Granqvist, C. G.
2011-06-01
Composites including VO2-based thermochromic nanoparticles are able to combine high luminous transmittance Tlum with a significant modulation of the solar energy transmittance ΔTsol at a "critical" temperature in the vicinity of room temperature. Thus nanothermochromics is of much interest for energy efficient fenestration and offers advantages over thermochromic VO2-based thin films. This paper presents calculations based on effective medium theory applied to dilute suspensions of core-shell nanoparticles and demonstrates that, in particular, moderately thin-walled hollow spherical VO2 nanoshells can give significantly higher values of ΔTsol than solid nanoparticles at the expense of a somewhat lowered Tlum. This paper is a sequel to a recent publication [S.-Y. Li, G. A. Niklasson, and C. G. Granqvist, J. Appl. Phys. 108, 063525 (2010)].
Ab Initio No-Core Shell Model Calculations Using Realistic Two- and Three-Body Interactions
Navratil, P; Ormand, W E; Forssen, C; Caurier, E
2004-11-30
There has been significant progress in the ab initio approaches to the structure of light nuclei. One such method is the ab initio no-core shell model (NCSM). Starting from realistic two- and three-nucleon interactions this method can predict low-lying levels in p-shell nuclei. In this contribution, we present a brief overview of the NCSM with examples of recent applications. We highlight our study of the parity inversion in {sup 11}Be, for which calculations were performed in basis spaces up to 9{Dirac_h}{Omega} (dimensions reaching 7 x 10{sup 8}). We also present our latest results for the p-shell nuclei using the Tucson-Melbourne TM three-nucleon interaction with several proposed parameter sets.
Analysis on burnup step effect for evaluating reactor criticality and fuel breeding ratio
Saputra, Geby; Purnama, Aditya Rizki; Permana, Sidik; Suzuki, Mitsutoshi
2014-09-30
Criticality condition of the reactors is one of the important factors for evaluating reactor operation and nuclear fuel breeding ratio is another factor to show nuclear fuel sustainability. This study analyzes the effect of burnup steps and cycle operation step for evaluating the criticality condition of the reactor as well as the performance of nuclear fuel breeding or breeding ratio (BR). Burnup step is performed based on a day step analysis which is varied from 10 days up to 800 days and for cycle operation from 1 cycle up to 8 cycles reactor operations. In addition, calculation efficiency based on the variation of computer processors to run the analysis in term of time (time efficiency in the calculation) have been also investigated. Optimization method for reactor design analysis which is used a large fast breeder reactor type as a reference case was performed by adopting an established reactor design code of JOINT-FR. The results show a criticality condition becomes higher for smaller burnup step (day) and for breeding ratio becomes less for smaller burnup step (day). Some nuclides contribute to make better criticality when smaller burnup step due to individul nuclide half-live. Calculation time for different burnup step shows a correlation with the time consuming requirement for more details step calculation, although the consuming time is not directly equivalent with the how many time the burnup time step is divided.
Applicability of 3D Monte Carlo simulations for local values calculations in a PWR core
NASA Astrophysics Data System (ADS)
Bernard, Franck; Cochet, Bertrand; Jinaphanh, Alexis; Jacquet, Olivier
2014-06-01
As technical support of the French Nuclear Safety Authority, IRSN has been developing the MORET Monte Carlo code for many years in the framework of criticality safety assessment and is now working to extend its application to reactor physics. For that purpose, beside the validation for criticality safety (more than 2000 benchmarks from the ICSBEP Handbook have been modeled and analyzed), a complementary validation phase for reactor physics has been started, with benchmarks from IRPHEP Handbook and others. In particular, to evaluate the applicability of MORET and other Monte Carlo codes for local flux or power density calculations in large power reactors, it has been decided to contribute to the "Monte Carlo Performance Benchmark" (hosted by OECD/NEA). The aim of this benchmark is to monitor, in forthcoming decades, the performance progress of detailed Monte Carlo full core calculations. More precisely, it measures their advancement towards achieving high statistical accuracy in reasonable computation time for local power at fuel pellet level. A full PWR reactor core is modeled to compute local power densities for more than 6 million fuel regions. This paper presents results obtained at IRSN for this benchmark with MORET and comparisons with MCNP. The number of fuel elements is so large that source convergence as well as statistical convergence issues could cause large errors in local tallies, especially in peripheral zones. Various sampling or tracking methods have been implemented in MORET, and their operational effects on such a complex case have been studied. Beyond convergence issues, to compute local values in so many fuel regions could cause prohibitive slowing down of neutron tracking. To avoid this, energy grid unification and tallies preparation before tracking have been implemented, tested and proved to be successful. In this particular case, IRSN obtained promising results with MORET compared to MCNP, in terms of local power densities, standard
Non-Invasive Characterization of Burnup for PWR Spent Fuel Rods with Burnups > 80 GWd/t
Caruso, S.; Murphy, M.; Jatuff, F.; Chawla, R.
2006-07-01
High-resolution gamma spectroscopy has been employed for the measurement of {sup 134}Cs/{sup 137}Cs, {sup 154}Eu/{sup 137}Cs and {sup 134}Cs/{sup 154}Eu gamma intensity ratios from spent fuel with the purpose of deriving pin-averaged single-ratio burnup indicators for high and ultra-high burnups. Two UO{sub 2} pressurised water reactor (PWR) fuel rod segments with record burnup levels >80 GWd/t have been experimentally characterised. Additionally, pin cell depletion calculations have been performed for each sample with the deterministic code CASMO-4, using both its JEF2.2- and its ENDF/B-IV-based libraries, for three different descriptions of the fuel rod irradiation histories, in order to test the sensitivity of the results to neutron cross sections and to the depletion model employed. Measured and calculated ratios have then been compared. It is shown that the {sup 134}Cs/{sup 137}Cs ratio, frequently used as burnup monitor, is considerably less accurate for values exceeding 50 GWd/t; discrepancies of up to {approx}25% are found between measured and calculated values. The ratios built with the {sup 154}Eu concentration show much larger discrepancies, essentially because this isotope is rather poorly predicted as revealed by the use of different basic cross section data. (authors)
Varivtsev, A. V. Zhemkov, I. Yu.
2014-12-15
The application of the improved method for calculating the radiation heat generation in the elements of an experimental device located at the periphery of the BOR-60 reactor core results in a significant reduction in the discrepancies between the calculated and the experimental data. This allows us to conclude that the improved method has an advantage over the one used earlier.
Jung, Y. S.; Joo, H. G.; Yoon, J. I.
2013-07-01
The nTRACER direct whole core transport code employing the planar MOC solution based 3-D calculation method, the subgroup method for resonance treatment, the Krylov matrix exponential method for depletion, and a subchannel thermal/hydraulic calculation solver was developed for practical high-fidelity simulation of power reactors. Its accuracy and performance is verified by comparing with the measurement data obtained for three pressurized water reactor cores. It is demonstrated that accurate and detailed multi-physic simulation of power reactors is practically realizable without any prior calculations or adjustments. (authors)
NASA Astrophysics Data System (ADS)
Zhivkov, P.; Furman, W.; Stoyanov, Ch
2014-09-01
The main characteristics of the neutron field formed within the massive (512 kg) natural uranium target assembly (TA) QUINTA irradiated by deuteron beam of JINR Nuclotron with energies 1,2,4, and 8 GeV as well as the spatial distributions and the integral numbers of (n,f), (n,γ) and (n,xn)- reactions were calculated and compared with experimental data [1] . The MCNPX 27e code with ISABEL/ABLA/FLUKA and INCL4/ABLA models of intra-nuclear cascade (INC) and experimental cross-sections of the corresponding reactions were used. Special attention was paid to the elucidation of the role of charged particles (protons and pions) in the fission of natural uranium of TA QUINTA. Extensive calculations have been done for quasi-infinite (with very small neutron leakage) depleted uranium TA BURAN having mass about 20 t which are intended to be used in experiments at Nuclotron in 2014-2016. As in the case of TA QUINTA which really models the central zone of TA BURAN the total numbers of fissions, produced 239Pu nuclei and total neutron multiplicities are predicted to be proportional to proton or deuteron energy up to 12 GeV. But obtained values of beam power gain are practically constant in studied incident energy range and are approximately four. These values are in contradiction with the experimental result [2] obtained for the depleted uranium core weighting three tons at incident proton energy 0.66 GeV.
Cluster form factor calculation in the ab initio no-core shell model
Navratil, Petr
2004-11-01
We derive expressions for cluster overlap integrals or channel cluster form factors for ab initio no-core shell model (NCSM) wave functions. These are used to obtain the spectroscopic factors and can serve as a starting point for the description of low-energy nuclear reactions. We consider the composite system and the target nucleus to be described in the Slater determinant (SD) harmonic oscillator (HO) basis while the projectile eigenstate to be expanded in the Jacobi coordinate HO basis. This is the most practical case. The spurious center of mass components present in the SD bases are removed exactly. The calculated cluster overlap integrals are translationally invariant. As an illustration, we present results of cluster form factor calculations for <{sup 5}He vertical bar{sup 4}He+n>, <{sup 5}He vertical bar{sup 3}H+d>, <{sup 6}Li vertical bar{sup 4}He+d>, <{sup 6}Be vertical bar{sup 3}He+{sup 3}He>, <{sup 7}Li vertical bar{sup 4}He+{sup 3}H>, <{sup 7}Li vertical bar{sup 6}Li+n>, <{sup 8}Be vertical bar{sup 6}Li+d>, <{sup 8}Be vertical bar{sup 7}Li+p>, <{sup 9}Li vertical bar{sup 8}Li+n>, and <{sup 13}C vertical bar{sup 12}C+n>, with all the nuclei described by multi-({Dirac_h}/2{pi}){omega} NCSM wave functions.
Gross Gamma Dose Rate Measurements for TRIGA Spent Nuclear Fuel Burnup Validation
Winston, P.L.; Sterbentz, J.W.
2002-07-01
Gross gamma-ray dose rates from six spent TRIGA fuel elements were measured and compared to calculated values as a means to validate the reported element burnups. A newly installed and functional gamma-ray detection subsystem of the In-Cell Examination System was used to perform the measurements and is described in some detail. The analytical methodology used to calculate the corresponding dose rates is presented along with the calculated values. Comparison of the measured and calculated dose rates for the TRIGA fuel elements indicates good agreement (less than a factor of 2 difference). The intent of the subsystem is to measure the gross gamma dose rate and correlate the measurement to a calculated dose rate based on the element s known burnup and other pertinent spent fuel information. Although validation of the TRIGA elements' burnup is of primary concern in this paper, the measurement and calculational techniques can be used to either validate an element's reported burnup or provide a burnup estimate for an element with an unknown burnup. (authors)
Gross Gamma Dose Rate Measurements for TRIGA Spent Nuclear Fuel Burnup Validation
Winston, Philip Lon; Sterbentz, James William
2001-04-01
Gross gamma-ray dose rates from six spent TRIGA fuel elements were measured and compared to calculated values as a means to validate the reported element burnups. A newly installed and functional gamma-ray detection subsystem of the In-Cell Examination System was used to perform the measurements and is described in some detail. The analytical methodology used to calculate the corresponding dose rates is presented along with the calculated values. Comparison of the measured and calculated dose rates for the TRIGA fuel elements indicates good agreement (less than a factor of 2 difference). The intent of the subsystem is to measure the gross gamma dose rate and correlate the measurement to a calculated dose rate based on the element s known burnup and other pertinent spent fuel information. Although validation of the TRIGA elements’ burnup is of primary concern in this paper, the measurement and calculational techniques can be used to either validate an element’s reported burnup or provide a burnup estimate for an element with an unknown burnup.
Raman spectra calculations for Si-Ge core-shell nanocrystals using ab initio real-space methods
NASA Astrophysics Data System (ADS)
Bobbitt, N. Scott; Chelikowsky, James R.
We use a real-space pseudopotential method within density functional theory to calculate Raman spectra for Si-Ge core-shell nanocrystals. We examine the lattice strain induced by the interface of the core and the shell. We calculate how this strain affects the vibrational modes and Raman spectra. We also find that the relative size of the Si and Ge peaks in the Raman spectrum is proportional to the size of the Si core and Ge shell regions, which suggests that Raman spectroscopy can be used to experimentally determine the relative size of the core and the outer shell in these nanocrystals. This work is supported by the DOE under Grant Number DE-FG02-06ER46286. Computations were performed on machines at TACC and NERSC.
Review of Technical Studies in the United States in Support of Burnup Credit Regulatory Guidance
Wagner, John C; Parks, Cecil V; Mueller, Don; Gauld, Ian C
2010-01-01
Taking credit for the reduction in reactivity associated with fuel depletion can enable more cost-effective, higher-density storage, transport, disposal, and reprocessing of spent nuclear fuel (SNF) while maintaining sufficient subcritical margin to establish an adequate safety basis. Consequently, there continues to be considerable interest in the United States (U.S.), as well as internationally, in the increased use of burnup credit in SNF operations, particularly related to storage, transport, and disposal of commercial SNF. This interest has motivated numerous technical studies related to the application of burnup credit, both domestically and internationally, as well as the design of SNF storage, transport and disposal systems that rely on burnup credit for maintaining subcriticality. Responding to industry requests and needs, the U.S. Nuclear Regulatory Commission (NRC) initiated a burnup credit research program in 1999, with support from the Oak Ridge National Laboratory (ORNL), to develop regulatory guidance and the supporting technical bases for allowing and expanding the use of burnup credit in pressurized-water reactor SNF storage and transport applications. Although this NRC research program has not been continuous since its inception, considerable progress has been achieved in many key areas in terms of increased understanding of relevant phenomena and issues, availability of relevant information and data, and subsequently updated regulatory guidance for expanded use of burnup credit. This paper reviews technical studies performed by ORNL for the U.S. NRC burnup credit research program. Examples of topics include reactivity effects associated with reactor operating characteristics, fuel assembly characteristics, burnable absorbers, control rods, spatial burnup distributions, cooling time, and assembly misloading; methods and data for validation of isotopic composition predictions; methods and data for validation of criticality calculations; and
Phenomena and Parameters Important to Burnup Credit
Parks, C.V.
2001-01-10
Since the mid-1980s, a significant number of studies have been directed at understanding the phenomena and parameters important to implementation of burnup credit in out-of-reactor applications involving pressurized-water-reactor (PWR) spent fuel. The efforts directed at burnup credit involving boiling-water-reactor (BWR) spent fuel have been more limited. This paper reviews the knowledge and experience gained from work performed in the US and other countries in the study of burnup credit. Relevant physics and analysis phenomenon are identified, and an assessment of their importance to burnup credit implementation for transport and dry cask storage is given.
NASA Astrophysics Data System (ADS)
Maris, P.; Caprio, M. A.; Vary, J. P.
2015-01-01
The emergence of rotational bands is observed in no-core configuration interaction (NCCI) calculations for the Be isotopes (7 ≤A ≤12 ), as evidenced by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. Yrast and low-lying excited bands are found. The results indicate well-developed rotational structure in NCCI calculations, using the JISP16 realistic nucleon-nucleon interaction within finite, computationally accessible configuration spaces.
Technical Development on Burn-up Credit for Spent LWR Fuel
Gauld, I.C.
2001-12-26
Technical development on burn-up credit for spent LWR fuels had been performed at JAERI since 1990 under the contract with Science and Technology Agency of Japan entitled ''Technical Development on Criticality Safety Management for Spent LWR Fuels.'' Main purposes of this work are to obtain the experimental data on criticality properties and isotopic compositions of spent LWR fuels and to verify burnup and criticality calculation codes. In this work three major experiments of exponential experiments for spent fuel assemblies to obtain criticality data, non-destructive gamma-ray measurement of spent fuel rods for evaluating axial burn-up profiles, and destructive analyses of spent fuel samples for determining precise burn-up and isotopic compositions were carried out. The measured data obtained were used for validating calculation codes as well as an examination of criticality safety analyses. Details of the work are described in this report.
NASA Astrophysics Data System (ADS)
Serrano-Purroy, D.; Clarens, F.; González-Robles, E.; Glatz, J. P.; Wegen, D. H.; de Pablo, J.; Casas, I.; Giménez, J.; Martínez-Esparza, A.
2012-08-01
Two weak points in Performance Assessment (PA) exercises regarding the alteration of Spent Nuclear Fuel (SNF) are the contribution of the so-called Instant Release Fraction (IRF) and the effect of High Burn-Up Structure (HBS). This manuscript focuses on the effect of HBS in matrix (long term) and instant release of a Pressurised Water Reactor (PWR) SNF irradiated in a commercial reactor with a mean Burn-Up (BU) of 60 GWd/tU. In order to study the HBS contribution, two samples from different radial positions have been prepared. One from the centre of the SNF, labelled CORE, and one from the periphery, enriched with HBS and labelled OUT. Static leaching experiments have been carried out with two synthetic leaching solutions: bicarbonate (BIC) and Bentonitic Granitic Groundwater (BGW), and in all cases under oxidising conditions. IRF values have been calculated from the determined Fraction of Inventory in Aqueous Phase (FIAP). In all studied cases, some radionuclides (RN): Rb, Sr and Cs, have shown higher release rates than uranium, especially at the beginning of the experiment, and have been considered as IRF. Redox sensitive RN like Mo and Tc have been found to dissolve slightly faster than uranium and further studies might be needed to confirm if they can also be considered part of the IRF. Most of the remaining studied RN, mainly actinides and lanthanides, have been found to dissolve congruently with the uranium matrix. Finally, Zr, Ru and Rh presented lower release rates than the matrix. Higher matrix release has been determined for CORE than for OUT samples showing that the formation of HBS might have a protective effect against the oxidative corrosion of the SNF. On the contrary, no significant differences have been observed between the two studied leaching solutions (BIC and BGW). Two different IRF contributions have been determined. One corresponding to the fraction of inventory segregated in the external open grain boundaries, directly available to water and
C 1s and N 1s core excitation of aniline: Experiment by electron impact and ab initio calculations
Duflot, D.; Flament, J.-P.; Giuliani, A.; Heinesch, J.; Grogna, M.; Hubin-Franskin, M.-J.
2007-05-15
Core shell excitation spectra of aniline at the carbon and nitrogen 1s edges have been obtained by inner-shell electron energy-loss spectroscopy recorded under scattering conditions where electric dipolar conditions dominate, with higher resolution than in the previous studies. They are interpreted with the aid of ab initio configuration interaction calculations. The spectrum at the C 1s edge is dominated by an intense {pi}{sup *} band. The calculated chemical shift due to the different chemical environment at the carbon 1s edge calculated is in agreement with the experimental observations within a few tenths of an eV. The transition energies of the most intense bands in the C 1s excitation spectrum are discussed at different levels of calculations. In the nitrogen 1s excitation spectrum the most intense bands are due to Rydberg-valence transitions involving the {sigma}{sup *}-type molecular orbitals, in agreement with the experiment. This assignment is different from that of extended Hueckel molecular orbital calculations. The geometries of the core excited states have been calculated and compared to their equivalent core molecules and benzene.
Calculation of the elastic properties of a triangular cell core for lightweight composite mirrors
NASA Astrophysics Data System (ADS)
Penado, F. Ernesto; Clark, James H., III; Walton, Joshua P.; Romeo, Robert C.; Martin, Robert N.
2007-09-01
The use of composite materials in the fabrication of optical telescope mirrors offers many advantages over conventional methods, including lightweight, portability and the potential for lower manufacturing costs. In the construction of the substrate for these mirrors, sandwich construction offers the advantage of even lower weight and higher stiffness. Generally, an aluminum or Nomex honeycomb core is used in composite applications requiring sandwich construction. However, the use of a composite core offers the potential for increased stiffness and strength, low thermal distortion compatible with that of the facesheets, the absence of galvanic corrosion and the ability to readily modify the core properties. In order to design, analyze and optimize these mirrors, knowledge of the mechanical properties of the core is essential. In this paper, the mechanical properties of a composite triangular cell core (often referred to as isogrid) are determined using finite element analysis of a representative unit cell. The core studied offers many advantages over conventional cores including increased thermal and dimensional stability, as well as low weight. Results are provided for the engineering elastic moduli of cores made of high stiffness composite material as a function of the ply layup and cell size. Finally, in order to illustrate the use of these properties in a typical application, a 1.4-m diameter composite mirror is analyzed using the finite element method, and the resulting stiffness and natural frequencies are presented.
van Wüllen, Christoph
2012-03-21
State-of-the art effective core potentials (ECPs) that replace electrons of inner atomic cores involve non-local potentials. If such an effective core potential is added to the Hamiltonian of a system in a magnetic field, the resulting Hamiltonian is not gauge invariant. This means, magnetic properties such as magnetisabilities and magnetic shieldings (or magnetic susceptibilities and nuclear magnetic resonance chemical shifts) calculated with different gauge origins are different even for exact solutions of the Schrödinger equation. It is possible to restore gauge invariance of the Hamiltonian by adding magnetic field dependent terms arising from the effective core potential. Numerical calculations on atomic and diatomic model systems (potassium mono-cation and potassium dimer) clearly demonstrate that the standard effective core potential Hamiltonian violates gauge invariance, and this affects the calculation of magnetisabilities more strongly than the calculation of magnetic shieldings. The modified magnetic field dependent effective core potential Hamiltonian is gauge invariant, and therefore it is the correct starting point for distributed gauge origin methods. The formalism for gauge including atomic orbitals (GIAO) and individual gauge for localized orbitals methods is worked out. ECP GIAO results for the potassium dimer are presented. The new method performs much better than a previous ECP GIAO implementation that did not account for the non-locality of the potential. For magnetic shieldings, deviations are clearly seen, but they amount to few ppm only. For magnetisabilities, our new ECP GIAO implementation is a major improvement, as demonstrated by the comparison of all-electron and ECP results. PMID:22443751
Accident source terms for light-water nuclear power plants using high-burnup or MOX fuel.
Salay, Michael; Gauntt, Randall O.; Lee, Richard Y.; Powers, Dana Auburn; Leonard, Mark Thomas
2011-01-01
Representative accident source terms patterned after the NUREG-1465 Source Term have been developed for high burnup fuel in BWRs and PWRs and for MOX fuel in a PWR with an ice-condenser containment. These source terms have been derived using nonparametric order statistics to develop distributions for the timing of radionuclide release during four accident phases and for release fractions of nine chemical classes of radionuclides as calculated with the MELCOR 1.8.5 accident analysis computer code. The accident phases are those defined in the NUREG-1465 Source Term - gap release, in-vessel release, ex-vessel release, and late in-vessel release. Important differences among the accident source terms derived here and the NUREG-1465 Source Term are not attributable to either fuel burnup or use of MOX fuel. Rather, differences among the source terms are due predominantly to improved understanding of the physics of core meltdown accidents. Heat losses from the degrading reactor core prolong the process of in-vessel release of radionuclides. Improved understanding of the chemistries of tellurium and cesium under reactor accidents changes the predicted behavior characteristics of these radioactive elements relative to what was assumed in the derivation of the NUREG-1465 Source Term. An additional radionuclide chemical class has been defined to account for release of cesium as cesium molybdate which enhances molybdenum release relative to other metallic fission products.
Mai, V. T.; Fujii, T.; Wada, K.; Kitada, T.; Takaki, N.; Yamaguchi, A.; Watanabe, H.; Unesaki, H.
2012-07-01
Considering the importance of thorium data and concerning about the accuracy of Th-232 cross section library, a series of experiments of thorium critical core carried out at KUCA facility of Kyoto Univ. Research Reactor Inst. have been analyzed. The core was composed of pure thorium plates and 93% enriched uranium plates, solid polyethylene moderator with hydro to U-235 ratio of 140 and Th-232 to U-235 ratio of 15.2. Calculations of the effective multiplication factor, control rod worth, reactivity worth of Th plates have been conducted by MVP code using JENDL-4.0 library [1]. At the experiment site, after achieving the critical state with 51 fuel rods inserted inside the reactor, the measurements of the reactivity worth of control rod and thorium sample are carried out. By comparing with the experimental data, the calculation overestimates the effective multiplication factor about 0.90%. Reactivity worth of the control rods evaluation using MVP is acceptable with the maximum discrepancy about the statistical error of the measured data. The calculated results agree to the measurement ones within the difference range of 3.1% for the reactivity worth of one Th plate. From this investigation, further experiments and research on Th-232 cross section library need to be conducted to provide more reliable data for thorium based fuel core design and safety calculation. (authors)
Nuclide Importance and the Steady-State Burnup Equation
Sekimoto, Hiroshi; Nemoto, Atsushi
2000-05-15
Conventional methods for evaluating some characteristic values of nuclides relating to burnup in a given neutron spectrum are reviewed in a mathematically systematic way, and a new method based on the importance theory is proposed. In this method, these characteristic values of a nuclide are equivalent to the importances of the nuclide. By solving the equation adjoint to the steady-state burnup equation with a properly chosen source term, the importances for all nuclides are obtained simultaneously.The fission number importance, net neutron importance, fission neutron importance, and absorbed neutron importance are evaluated and discussed. The net neutron importance is a measure directly estimating neutron economy, and it can be evaluated simply by calculating the fission neutron importance minus the absorbed neutron importance, where only the absorbed neutron importance depends on the fission product. The fission neutron importance and absorbed neutron importance are analyzed separately, and detailed discussions of the fission product effects are given for the absorbed neutron importance.
Recommendations for Addressing Axial Burnup in the PWR Burnup Credit Analyses
Wagner, J.C.
2002-10-23
This report presents studies performed to support the development of a technically justifiable approach for addressing the axial-burnup distribution in pressurized-water reactor (PWR) burnup-credit criticality safety analyses. The effect of the axial-burnup distribution on reactivity and proposed approaches for addressing the axial-burnup distribution are briefly reviewed. A publicly available database of profiles is examined in detail to identify profiles that maximize the neutron multiplication factor, k{sub eff}, assess its adequacy for PWR burnup credit analyses, and investigate the existence of trends with fuel type and/or reactor operations. A statistical evaluation of the k{sub eff} values associated with the profiles in the axial-burnup-profile database was performed, and the most reactive (bounding) profiles were identified as statistical outliers. The impact of these bounding profiles on k{sub eff} is quantified for a high-density burnup credit cask. Analyses are also presented to quantify the potential reactivity consequence of loading assemblies with axial-burnup profiles that are not bounded by the database. The report concludes with a discussion on the issues for consideration and recommendations for addressing axial burnup in criticality safety analyses using burnup credit for dry cask storage and transportation.
A study on fast reactor core mechanics by an ex-reactor test and comparisons with calculations
Tottori, Shoji; Kawanaka, Ikunori; Nakagawa, Masatoshi; Arie, Kazuo; Itoh, Kunihiro; Ohya, Takeaki; Motomiya, Takeo; Adachi, Hironori
1996-07-01
This paper presents and discusses the results of core bowing experiments performed with an ex-reactor rig holding a half hexagon array of 22 sub-assemblies (S/As) simulating the Japanese DFBR conditions and the comparisons of the measured results with calculations by individually developed codes--ARKAS, RAINBOW, SANBOW. The main conclusions of this study are (1) interwrapper loads and S/A displacements within the array were measured at selected positions for a series of five tests simulating the DFBR core bowing modes, (2) the overall comparison between the non-friction calculation and measurement showed good agreement for loads, displacements and their directions, and (3) validation of the friction algorithm has also been carried out and further improvement of the agreement was obtained.
Chi, C.-C.; Hsiao, C.-H.; Ouyang, Chuenhou; Skoropata, E.; Lierop, J. van
2015-05-07
Significant efforts towards understanding bi-magnetic core-shell nanoparticles are underway currently as they provide a pathway towards properties unavailable with single-phased systems. Recently, we have demonstrated that the magnetism of γ-Fe2O3/CoO core-shell nanoparticles, in particular, at high temperatures, originates essentially from an interfacial doped iron-oxide layer that is formed by the migration of Co{sup 2+} from the CoO shell into the surface layers of the γ-Fe2O3 core [Skoropata et al., Phys. Rev. B 89, 024410 (2014)]. To examine directly the nature of the intermixed layer, we have used high-resolution transmission electron microscopy (HRTEM) and first-principles calculations to examine the impact of the core-shell intermixing at the atomic level. By analyzing the HRTEM images and energy dispersive spectra, the level and nature of intermixing was confirmed, mainly as doping of Co into the octahedral site vacancies of γ-Fe2O3. The average Co doping depths for different processing temperatures (150 °C and 235 °C) were 0.56 nm and 0.78 nm (determined to within 5% through simulation), respectively, establishing that the amount of core-shell intermixing can be altered purposefully with an appropriate change in synthesis conditions. Through first-principles calculations, we find that the intermixing phase of γ-Fe2O3 with Co doping is ferromagnetic, with even higher magnetization as compared to that of pure γ-Fe2O3. In addition, we show that Co doping into different octahedral sites can cause different magnetizations. This was reflected in a change in overall nanoparticle magnetization, where we observed a 25% reduction in magnetization for the 235 °C versus the 150 °C sample, despite a thicker intermixed layer.
Code System for 2-Group, 3D Neutronic Kinetics Calculations Coupled to Core Thermal Hydraulics.
Energy Science and Technology Software Center (ESTSC)
2000-05-12
Version 00 QUARK is a combined computer program comprising a revised version of the QUANDRY three-dimensional, two-group neutron kinetics code and an upgraded version of the COBRA transient core analysis code (COBRA-EN). Starting from either a critical steady-state (k-effective or critical dilute Boron problem) or a subcritical steady-state (fixed source problem) in a PWR plant, the code allows one to simulate the neutronic and thermal-hydraulic core transient response to reactivity accidents initiated both inside themore » vessel (such as a control rod ejection) and outside the vessel (such as the sudden change of the Boron concentration in the coolant). QUARK output can be used as input to PSR-470/NORMA-FP to perform a subchannel analysis from converged coarse-mesh nodal solutions.« less
Simulating the Dynamics of Earth's Core: Using NCCS Supercomputers Speeds Calculations
NASA Technical Reports Server (NTRS)
2002-01-01
If one wanted to study Earth's core directly, one would have to drill through about 1,800 miles of solid rock to reach liquid core-keeping the tunnel from collapsing under pressures that are more than 1 million atmospheres and then sink an instrument package to the bottom that could operate at 8,000 F with 10,000 tons of force crushing every square inch of its surface. Even then, several of these tunnels would probably be needed to obtain enough data. Faced with difficult or impossible tasks such as these, scientists use other available sources of information - such as seismology, mineralogy, geomagnetism, geodesy, and, above all, physical principles - to derive a model of the core and, study it by running computer simulations. One NASA researcher is doing just that on NCCS computers. Physicist and applied mathematician Weijia Kuang, of the Space Geodesy Branch, and his collaborators at Goddard have what he calls the,"second - ever" working, usable, self-consistent, fully dynamic, three-dimensional geodynamic model (see "The Geodynamic Theory"). Kuang runs his model simulations on the supercomputers at the NCCS. He and Jeremy Bloxham, of Harvard University, developed the original version, written in Fortran 77, in 1996.
Guerin, P.; Baudron, A. M.; Lautard, J. J.
2006-07-01
This paper describes a new technique for determining the pin power in heterogeneous core calculations. It is based on a domain decomposition with overlapping sub-domains and a component mode synthesis technique for the global flux determination. Local basis functions are used to span a discrete space that allows fundamental global mode approximation through a Galerkin technique. Two approaches are given to obtain these local basis functions: in the first one (Component Mode Synthesis method), the first few spatial eigenfunctions are computed on each sub-domain, using periodic boundary conditions. In the second one (Factorized Component Mode Synthesis method), only the fundamental mode is computed, and we use a factorization principle for the flux in order to replace the higher order Eigenmodes. These different local spatial functions are extended to the global domain by defining them as zero outside the sub-domain. These methods are well-fitted for heterogeneous core calculations because the spatial interface modes are taken into account in the domain decomposition. Although these methods could be applied to higher order angular approximations - particularly easily to a SPN approximation - the numerical results we provide are obtained using a diffusion model. We show the methods' accuracy for reactor cores loaded with UOX and MOX assemblies, for which standard reconstruction techniques are known to perform poorly. Furthermore, we show that our methods are highly and easily parallelizable. (authors)
Three-Dimensional Full Core Power Calculations for Pressurized Water Reactors
Evans, Thomas M; Davidson, Gregory G; Slaybaugh, Rachel N
2010-01-01
We have implemented a new multilevel parallel decomposition in the Denovo discrete ordinates radiation transport code. In concert with Krylov subspace iterative solvers, the multilevel decomposition allows concurrency over energy in addition to space-angle. The original space-angle partitioning in Denovo placed an eective limit on the scalability of the transport solver that was highly dependent on the problem size. The added phase-space concurrency combined with the high-performance Krylov solvers has enabled weak scaling to 100K cores on the Jaguar XT5 supercomputer. Furthermore, the multilevel decomposition provides enough concurrency to scale to exascale computing and beyond.
Weak-coupling calculations in the /sup 208/Pb core region
McGrory, J.B.
1981-01-01
The structure of Tl and Hg isotopes near the /sup 208/Pb core is discussed in terms of a weak-coupling model where low-lying proton-hole states are coupled to low-lying neutron-hole states. The even Pb isotopes are first discussed in terms of a generalized seniority-2 approximation which is shown to be very accurate. The weak-coupling results are not satisfactory, and it is suggested that the defeat may be in the n-p interaction.
No-Core Shell Model Calculations in Light Nuclei with Three-Nucleon Forces
Barrett, B R; Vary, J P; Nogga, A; Navratil, P; Ormand, W E
2004-01-08
The ab initio No-Core Shell Model (NCSM) has recently been expanded to include nucleon-nucleon (NN) and three-nucleon (3N) interactions at the three-body cluster level. Here it is used to predict binding energies and spectra of p-shell nuclei based on realistic NN and 3N interactions. It is shown that 3N force (3NF) properties can be studied in these nuclear systems. First results show that interactions based on chiral perturbation theory lead to a realistic description of {sup 6}Li.
MCWO - Linking MCNP And ORIGEN2 For Fuel Burnup Analysis
Gray S Chang
2005-04-01
The UNIX BASH (Bourne Again Shell) script MCWO has been developed at the Idaho National Engineering and Environment Laboratory (INEEL) to couple the Monte Carlo transport code MCNP with the depletion and buildup code ORIGEN2. MCWO is a fully automated tool that links the Monte Carlo transport code MCNP with the radioactive decay and burnup code ORIGEN2. MCWO can handle a large number of fuel burnup and material loading specifications, Advanced Test Reactor (ATR) powers, and irradiation time intervals. The program processes input from the user that specifies the system geometry, initial material compositions, feed/removal specifications, and other code-specific parameters. Calculated results from MCNP, ORIGEN2, and data process module calculations are then output successively as the code runs. The principal function of MCWO is to transfer one-group cross-section and flux values from MCNP to ORIGEN2, and then transfer the resulting material compositions (after irradiation and/or decay) from ORIGEN2 back to MCNP in a repeated, cyclic fashion. The basic requirement of the code is that the user have a working MCNP input file and other input parameters; all interaction with ORIGEN2 and other calculations are performed by UNIX BASH script MCWO. This paper presents the MCWO-calculated results of the RERTR-1 and -2 , and the Weapons-Grade Mixed Oxiide fuel (Wg-MOX) fuel experiments in ATR and compares the MCWO-calculated results with the measured data.
NASA Astrophysics Data System (ADS)
Zhao, Zhongxiang
Modular Pebble Bed Reactor (MPBR) is a high temperature gas-cooled nuclear power reactor. This project investigated the feasibility of using the passive neutron counting and active neutron/gamma counting for the on line fuel burnup measurement for MPBR. To investigate whether there is a correlation between neutron emission and fuel burnup, the MPBR fuel depletion was simulated under different irradiation conditions by ORIGEN2. It was found that the neutron emission from an irradiated pebble increases with burnup super-linearly and reaches to 104 neutron/sec/pebble at the discharge burnup. The photon emission from an irradiated pebble was found to be in the order of 1013 photon/sec/pebble at all burnup levels. Analysis shows that the neutron emission rate of an irradiated pebble is sensitive to its burnup history and the spectral-averaged one-group cross sections used in the depletion calculations, which consequently leads to large uncertainty in the correlation between neutron emission and burnup. At low burnup levels, the uncertainty in the neutron emission/burnup correlation is too high and the neutron emission rate is too low so that it is impossible to determine a pebble's burnup by on-line neutron counting at low burnup levels. At high burnup levels, the uncertainty in the neutron emission rate becomes less but is still large in quantity. However, considering the super-linear feature of the correlation, the uncertainty in burnup determination was found to be ˜7% at the discharge burnup, which is acceptable. Therefore, total neutron emission rate of a pebble can be used as a burnup indicator to determine whether a pebble should be discharged or not. The feasibility of using passive neutron counting methods for the on-line burnup measurement was investigated by using a general Monte Carlo code, MCNP, to assess the detectability of the neutron emission and the capability to discriminate gamma noise by commonly used neutron detectors. It was found that both He-3
Some Nuclear Calculations of U-235-D2O Gaseous-Core Cavity Reactors
NASA Technical Reports Server (NTRS)
Ragsdale, Robert G.; Hyland, Robert E.
1961-01-01
The results of a multigroup, diffusion theory study of spherical gaseous-core cavity reactors are presented in this report. The reactor cavity of gaseous U235 is enclosed by a region of hydrogen gas and is separated from an external D2O moderator-reflector by a zirconium structural shell. Some cylindrical reactors are also investigated. A parametric study of spherical reactors indicates that, for the range of variables studied, critical mass increases as: (1) Fuel region is compressed within the reactor cavity, (2) moderator thickness is decreased, (3) structural shell thickness is increased, and (4) moderator temperature is increased. A buckling analogy is used to estimate the critical mass of fully reflected cylindrical reactors from spherical results without fuel compression. For a reactor cavity of a 120-centimeter radius uniformly filled with fuel, no structural shell, a moderator temperature of 70 F, and a moderator thickness of 100 centimeters, the critical mass of a spherical reactor is 3.1 kilograms while that of a cylinder with a length-to-diameter ratio of 1.0 (L/D = 1) is approximately 3.8 kilograms and, with L/D = 2, 5.9 kilograms. For the range of variables considered for U235-D2O gaseous-core cavity reactors, the systems are characterized by 95 to 99 percent thermal absorptions, with the flux reaching a maximum in the moderator about 10 to 15 centimeters from the reactor cavity.
González-Cataldo, F.; Wilson, Hugh F.; Militzer, B.
2014-05-20
By combining density functional molecular dynamics simulations with a thermodynamic integration technique, we determine the free energy of metallic hydrogen and silica, SiO{sub 2}, at megabar pressures and thousands of degrees Kelvin. Our ab initio solubility calculations show that silica dissolves into fluid hydrogen above 5000 K for pressures from 10 and 40 Mbars, which has implications for the evolution of rocky cores in giant gas planets like Jupiter, Saturn, and a substantial fraction of known extrasolar planets. Our findings underline the necessity of considering the erosion and redistribution of core materials in giant planet evolution models, but they also demonstrate that hot metallic hydrogen is a good solvent at megabar pressures, which has implications for high-pressure experiments.
Development of burnup dependent fuel rod model in COBRA-TF
NASA Astrophysics Data System (ADS)
Yilmaz, Mine Ozdemir
The purpose of this research was to develop a burnup dependent fuel thermal conductivity model within Pennsylvania State University, Reactor Dynamics and Fuel Management Group (RDFMG) version of the subchannel thermal-hydraulics code COBRA-TF (CTF). The model takes into account first, the degradation of fuel thermal conductivity with high burnup; and second, the fuel thermal conductivity dependence on the Gadolinium content for both UO2 and MOX fuel rods. The modified Nuclear Fuel Industries (NFI) model for UO2 fuel rods and Duriez/Modified NFI Model for MOX fuel rods were incorporated into CTF and fuel centerline predictions were compared against Halden experimental test data and FRAPCON-3.4 predictions to validate the burnup dependent fuel thermal conductivity model in CTF. Experimental test cases from Halden reactor fuel rods for UO2 fuel rods at Beginning of Life (BOL), through lifetime without Gd2O3 and through lifetime with Gd 2O3 and a MOX fuel rod were simulated with CTF. Since test fuel rod and FRAPCON-3.4 results were based on single rod measurements, CTF was run for a single fuel rod surrounded with a single channel configuration. Input decks for CTF were developed for one fuel rod located at the center of a subchannel (rod-centered subchannel approach). Fuel centerline temperatures predicted by CTF were compared against the measurements from Halden experimental test data and the predictions from FRAPCON-3.4. After implementing the new fuel thermal conductivity model in CTF and validating the model with experimental data, CTF model was applied to steady state and transient calculations. 4x4 PWR fuel bundle configuration from Purdue MOX benchmark was used to apply the new model for steady state and transient calculations. First, one of each high burnup UO2 and MOX fuel rods from 4x4 matrix were selected to carry out single fuel rod calculations and fuel centerline temperatures predicted by CTF/TORT-TD were compared against CTF /TORT-TD /FRAPTRAN
NASA Astrophysics Data System (ADS)
Mejdoubi, Abdelilah; Brosseau, Christian
2007-11-01
Herein, we report finite-element calculations of the effective (relative) permittivity of composite materials consisting of inclusions and inclusion arrays with a core-shell structure embedded in a surrounding host. The material making up the core of the two-dimensional structures, or cross sections of infinite three-dimensional objects (parallel, infinitely long, and identical cylinders) where the properties and characteristics are invariant along the perpendicular cross sectional plane, is assumed to have a negative real part of the permittivity, while the coating material (annular shell) is considered to be lossless. While strictly valid only in a dc situation, our analysis can be extended to treat electric fields that oscillate with time, provided that the wavelengths and attenuation lengths associated with the fields are much larger than the microstructure dimension in order that the homogeneous (effective-medium) representation of the composite structure makes sense. While one may identify features of the electrostatic resonance (ER) which are common to core-shell structures characterized by permittivities with real parts of opposite signs, it appears that the predicted ER positions are sensitive to the shell thickness and can be tuned through varying this geometric parameter. For example, we observe that the ER is broadened and shifted as the loss and the shell thickness are increased, respectively. We also argue that such core shell may also be valuable in controlling ER characteristics via polarization in an external electric field. In addition, by considering calculations of the electric field distribution, we find that the ER results in very strong and local-field enhancements into small parts of the shell perimeter. Our findings open up possibilities for the development of hybrid structures that could exploit the ER features for a particular application.
Depletion analysis of the UMLRR reactor core using MCNP6
NASA Astrophysics Data System (ADS)
Odera, Dim Udochukwu
Accurate knowledge of the neutron flux and temporal nuclide inventory in reactor physics calculations is necessary for a variety of application in nuclear engineering such as criticality safety, safeguards, and spent fuel storage. The Monte Carlo N- Particle (MCNP6) code with integrated buildup depletion code (CINDER90) provides a high-fidelity tool that can be used to perform 3D, full core simulation to evaluate fissile material utilization, and nuclide inventory calculations as a function of burnup. The University of Massachusetts Lowell Research Reactor (UMLRR) reactor has been modeled with the deterministic based code, VENTURE and with an older version of MCNP (MCNP5). The MIT developed MCODE (MCNP ORIGEN DEPLETION CODE) was used previously to perform some limited depletion calculations. This work chronicles the use of MCNP6, released in June 2013, to perform coupled neutronics and depletion calculation. The results are compared to previously benchmarked results. Furthermore, the code is used to determine the ratio of fission products 134Cs and 137Cs (burnup indicators), and the resultant ratio is compared to the burnup of the UMLRR.
60Co as AN On-Line Burnup Indicator for Multi-Pass Pebble Bed Reactors
NASA Astrophysics Data System (ADS)
Hawari, Ayman I.; Chen, Jianwei
2003-06-01
Multi-pass pebble bed reactor concepts are characterized by circulating fuel systems that cycle the pebbles in and out of the core until the burnup limit is reached. Currently modular designs of such reactors, with nominal powers of approximately 300 MW-thermal, are under consideration for deployment internationally. A concern of the proposed designs is the ability to perform online measurements of the fuel burnup to determine whether a pebble has reached its end-of-life burnup limit (~ 80,000 MWD/MTU). In this work, computational simulations were performed to assess the utilization of a passive gamma ray spectrometric approach to perform this task. However, in addition to using the inherent signatures of the irradiated fuel, the use of the 59Co(n,γ)60Co reaction as a burnup indicator is considered. The results show that the activity ratio of 134Cs/60Co can provide an indicator that is accurate to within 5% at burnup greater than 20,000 MWD/MTU as the power is varied between 50% and 200% of the reactor's thermal power.
Exarhos, C.A.
1985-06-20
The US Department of Energy, Consumers Power Company, Exxon Nuclear Company, and General Public Utilities Nuclear Corporation have participated since 1979 in a cooperative Extended Burnup Demonstration Program. Under the program, standard ENC-fabricated reload fuel in the Big Rock Point and Oyster Creek reactor cores has been irradiated to discharge burnups at or beyond 35,000 MWD/MTU, one to two cycles beyond its originally projected exposure life. The program provides for examination of the fuel at poolside before and after each extended burnup cycle as well as for limited destructive hot cell examination. The 1984 progress report covers work performed under the EBD program between April 1983 and March 1984. Major milestones reached during the period include completion of a hot cell examination on four high burnup rods from Big Rock Point and of a poolside on the Oyster Creek EBD fuel at discharge. The hot cell examination of four rods at burnups to 37.2 GWD/MTU confirmed poolside measurements on the same fuel, showing the urania and gadolinia-bearing fuel rods to be in excellent condition. No major cladding degradation, pellet restructuring, or pellet-clad interaction was found in any of the samples examined. The Oyster Creek fuel, examined at an assembly average exposure of 34.5 GWD/MTU, showed good performance with regard to both diametral creepdown and clad oxide accumulation.
Investigation of Burnup Credit Modeling Issues Associated with BWR Fuel
Wagner, J.C.
2000-10-12
/or decreasing the neutron absorber concentration. However, regulations associated with permanent disposal require consideration of scenarios and/or package conditions that are not relevant or credible for storage or transportation, and as a result, necessitate credit for burnup in BWR fuel to maintain capacity objectives. Burnup credit relies on depletion calculations to provide a conservative estimate of spent fuel contents and subsequent criticality calculations to assess the value of k{sub eff} for a spent fuel cask or a fuel configuration under a variety of postulated conditions. Therefore, validation is necessary to quantify biases and uncertainties between analytic predictions and measured isotopics. However, the design and operational aspects of BWRs result in a more heterogeneous and time-varying reactor configuration than those of PWRs. Thus, BWR spent fuel analyses and validation efforts are significantly more complicated than those of their PWR counterparts. BWR spent fuel assemblies are manufactured with variable enrichments, both radially and axially, are exposed to time- and spatially-varying void distributions, contain integral burnable absorber rods, and are subject to partial control-blade insertion during operation. The latter is especially true in older fuel assemblies. Away-from-reactor depletion tools used for characterization of spent fuel have typically been developed and validated for more homogeneous PWR fuel assemblies without integral burnable absorber rods, and thus must be reassessed for BWR configurations to determine a conservative methodology for estimating the isotopic content of spent BWR fuel. This report examines the use of SAS2H8 for calculating spent BWR fuel isotopics for burnup-credit criticality safety analyses and assesses the adequacy of SAS2H for this task. The effects of SAS2H modeling assumptions on calculated spent BWR fuel isotopics and the effects of depletion assumptions on calculated k{sub inf} values are investigated. Detailed
Burnup Credit Approach Used in the Yucca Mountain License Application
Scaglione, John M; Wagner, John C
2010-01-01
The United States Department of Energy has submitted a license application (LA) for construction authorization of a deep geologic repository at Yucca Mountain, Nevada. The license application is currently under review by the United States Nuclear Regulatory Commission (NRC). This paper will describe the methodology and approach used in the LA to address the issue of criticality and the role of burnup credit during the postclosure period. The most significant and effective measures for prevention of criticality in the repository include multiple redundant barriers that act to isolate fissionable material from water (which can act as a moderator, corrosive agent, and transporter of fissile material); inherent geometry of waste package internals and waste forms; presence of fixed neutron absorbers in waste package internals; and fuel burnup for commercial spent nuclear fuel. A probabilistic approach has been used to screen criticality from the total system performance assessment. Within the probabilistic approach, criticality is considered an event, and the total probability of a criticality event occurring within 10,000 years of disposal is calculated and compared against the regulatory criterion. The total probability of criticality includes contributions associated with both internal (within waste packages) and external (external to waste packages) criticality for each of the initiating events that could lead to waste package breach. The occurrence of and conditions necessary for criticality in the repository have been thoroughly evaluated using a comprehensive range of parameter distributions. A simplified design-basis modeling approach has been used to evaluate the probability of criticality by using numerous significant and conservative assumptions. Burnup credit is used only for evaluations of in-package configurations and uses a combination of conservative and bounding modeling approximations to ensure conservatism. This paper will review the NRC regulatory
NASA Astrophysics Data System (ADS)
Shen, Zhichao; Ni, Sidao; Wu, Wenbo; Sun, Daoyuan
2016-04-01
The core-mantle boundary (CMB) topography plays a key role in constraining geodynamic modeling and core-mantle coupling. It's effective to resolve the intermediate lateral scale topography (hundreds of km) with short period core reflected seismic phases (ScP) due to their small Fresnel-zones at short epicentral distances (<3336 km (30°)). We developed a method based on the ray theory and representation theorem to calculate short period ScP synthetics for intermediate lateral scale CMB topography. The CMB topography we introduced here is axisymmetric and specified with two parameters: H (height) and L (diameter, or lateral length scale). Our numerical computation shows that a bump (H > 0) and dip (H < 0) model would cause defocusing/weakening and focusing/amplifying effects on ScP amplitude. Moreover, the effect of frequency and combination of L and H are quantified with the amplification coefficients. Then we applied this method to estimate a possible CMB topography beneath northeastern Japan, and a CMB model with L = 140 km, H = 1.2 km overall matches the observed pattern of 2D PcP/ScP amplitude ratios. However, it is difficult to totally rule out other factors that may also affect PcP/ScP pattern because of limitation of ray-based algorithms we used here. A hybrid method combining ray theory and numerical method is promising for studying complicated 3D structure and CMB topography in the future.
NASA Astrophysics Data System (ADS)
Correa, Alfredo; Schleife, Andre; Kanai, Yosuke
2014-03-01
In order to understand the interaction of projectile atoms with targets under particle radiation in materials, e.g. in space applications or nuclear reactors, it is critical to investigate electronic and ionic contributions to stopping power. The goal of such efforts is detailed understanding of radiation damages as well as fundamental effects such as ion-electron interaction. While ionic stopping has been successfully modeled by molecular dynamics in the past, only recently a computational framework came within reach that is capable of accurately describing electronic stopping from first principles. Using our large-scale implementation of real-time time-dependent density functional theory in non-adiabatic Ehrenfest molecular dynamics, we are able to gain deep insight into electronic stopping for systems with hundreds of atoms and thousands of electrons, taking into account their quantum-mechanical electron-electron interaction. We discuss distinct contributions of valence and core electrons of aluminum target atoms to electronic stopping, and we study their importance for different projectile (hydrogen and helium atoms) velocities. There is striking influence of the stopping geometry especially for fast projectiles, and we find excellent agreement with experiment. Prepared by LLNL under Contract DE-AC52-07NA27344.
Ab initio no core calculations of light nuclei and preludes to Hamiltonian quantum field theory
Vary, J. P.; Maris, P.; Honkanen, H.; Li, J.; Shirokov, A. M.; Brodsky, S. J.; Harindranath, A.
2009-12-17
Recent advances in ab initio quantum many-body methods and growth in computer power now enable highly precise calculations of nuclear structure. The precision has attained a level sufficient to make clear statements on the nature of 3-body forces in nuclear physics. Total binding energies, spin-dependent structure effects, and electroweak properties of light nuclei play major roles in pinpointing properties of the underlying strong interaction. Eventually, we anticipate a theory bridge with immense predictive power from QCD through nuclear forces to nuclear structure and nuclear reactions. Light front Hamiltonian quantum field theory offers an attractive pathway and we outline key elements.
Ab initio no core calculations of light nuclei and preludes to Hamiltonian quantum field theory
Vary, J.P.; Maris, P.; Shirokov, A.M.; Honkanen, H.; li, J.; Brodsky, S.J.; Harindranath, A.; Teramond, G.F.de; /Costa Rica U.
2009-08-03
Recent advances in ab initio quantum many-body methods and growth in computer power now enable highly precise calculations of nuclear structure. The precision has attained a level sufficient to make clear statements on the nature of 3-body forces in nuclear physics. Total binding energies, spin-dependent structure effects, and electroweak properties of light nuclei play major roles in pinpointing properties of the underlying strong interaction. Eventually,we anticipate a theory bridge with immense predictive power from QCD through nuclear forces to nuclear structure and nuclear reactions. Light front Hamiltonian quantum field theory offers an attractive pathway and we outline key elements.
NASA Astrophysics Data System (ADS)
Raza, Zamaan; Shulumba, Nina; Caffrey, Nuala M.; Dubrovinsky, Leonid; Abrikosov, Igor A.
2015-06-01
A recently discovered phase of orthorhombic iron carbide o-Fe7C3 [Prescher et al., Nat. Geosci. 8, 220 (2015), 10.1038/ngeo2370] is assessed as a potentially important phase for interpretation of the properties of the Earth's core. In this paper, we carry out first-principles calculations on o-Fe7C3 , finding properties to be in broad agreement with recent experiments, including a high Poisson's ratio (0.38). Our enthalpy calculations suggest that o-Fe7C3 is more stable than Eckstrom-Adcock hexagonal iron carbide (h-Fe7C3 ) below approximately 100 GPa. However, at 150 GPa, the two phases are essentially degenerate in terms of Gibbs free energy, and further increasing the pressure towards Earth's core conditions stabilizes h-Fe7C3 with respect to the orthorhombic phase. Increasing the temperature tends to stabilize the hexagonal phase at 360 GPa, but this trend may change beyond the limit of the quasiharmonic approximation.
Gauntt, Randall O.; Ross, Kyle W.; Smith, James Dean; Longmire, Pamela
2010-04-01
The Oak Ridge National Laboratory computer code, ORIGEN2.2 (CCC-371, 2002), was used to obtain the elemental composition of irradiated low-enriched uranium (LEU)/mixed-oxide (MOX) pressurized-water reactor fuel assemblies. Described in this report are the input parameters for the ORIGEN2.2 calculations. The rationale for performing the ORIGEN2.2 calculation was to generate inventories to be used to populate MELCOR radionuclide classes. Therefore the ORIGEN2.2 output was subsequently manipulated. The procedures performed in this data reduction process are also described herein. A listing of the ORIGEN2.2 input deck for two-cycle MOX is provided in the appendix. The final output from this data reduction process was three tables containing the radionuclide inventories for LEU/MOX in elemental form. Masses, thermal powers, and activities were reported for each category.
Fuel Modelling at Extended Burnup: IAEA Coordinated Research Project FUMEX-II
Killeen, J.C.; Turnbull, J.A.; Sartori, E.
2007-07-01
The International Atomic Energy Agency sponsored a Coordinated Research Project on Fuel Modelling at Extended Burnup (FUMEX-II). Eighteen fuel modelling groups participated with the intention of improving their capabilities to understand and predict the behaviour of water reactor fuel at high burnups. The exercise was carried out in coordination with the OECD/NEA. The participants used a mixture of data derived from actual irradiation histories of high burnup experimental fuel and commercial irradiations where post-irradiation examination measurements are available, combined with idealised power histories intended to represent possible future extended dwell commercial irradiations and test code capabilities at high burnup. All participants have been asked to model nine priority cases out of some 27 cases made available to them for the exercise from the IAEA/OECD International Fuel Performance Experimental Database. Calculations carried out by the participants, particularly for the idealised cases, have shown how varying modelling assumptions affect the high burnup predictions, and have led to an understanding of the requirements of future high burnup experimental data to help discriminate between modelling assumptions. This understanding is important in trying to model transient and fault behaviour at high burnup. It is important to recognise that the code predictions presented here should not be taken to indicate that some codes do not perform well. The codes have been designed for different applications and have differing assumptions and validation ranges; for example codes intended to predict Candu fuel operation with thin wall collapsible cladding do not need the clad creep and gap conductivity modelling found in PWR codes. Therefore, when a case is based on Candu technology or PWR technology, it is to be expected that the codes may not agree. However, it is the very differences in such behaviour that is useful in helping to understand the effects of such
Summary of high burnup fuel issues and NRC`s plan of action
Meyer, R.O.
1997-01-01
For the past two years the Office of Nuclear Regulatory Research has concentrated mostly on the so-called reactivity-initiated accidents -- the RIAs -- in this session of the Water Reactor Safety Information Meeting, but this year there is a more varied agenda. RIAs are, of course, not the only events of interest for reactor safety that are affected by extended burnup operation. Their has now been enough time to consider a range of technical issues that arise at high burnup, and a list of such issues being addressed in their research program is given here. (1) High burnup capability of the steady-state code (FRAPCON) used for licensing audit calculations. (2) General capability (including high burnup) of the transient code (FRAPTRAN) used for special studies. (3) Adequacy at high burnup of fuel damage criteria used in regulation for reactivity accidents. (4) Adequacy at high burnup of models and fuel related criteria used in regulation for loss-of-coolant accidents (LOCAs). (5) Effect of high burnup on fuel system damage during normal operation, including control rod insertion problems. A distinction is made between technical issues, which may or may not have direct licensing impacts, and licensing issues. The RIAs became a licensing issue when the French test in CABRI showed that cladding failures could occur at fuel enthalpies much lower than a value currently used in licensing. Fuel assembly distortion became a licensing issue when control rod insertion was affected in some operating plants. In this presentation, these technical issues will be described and the NRC`s plan of action to address them will be discussed.
Coolant Density and Control Blade History Effects in Extended BWR Burnup Credit
Ade, Brian J; Marshall, William BJ J; Bowman, Stephen M; Gauld, Ian C; Ilas, Germina; Martinez-Gonzalez, Jesus S
2015-01-01
Oak Ridge National Laboratory and the US Nuclear Regulatory Commission have initiated a multiyear project to investigate the application of burnup credit (BUC) for boiling water reactor (BWR) fuel in storage and transportation casks. This project includes two phases. The first phase investigates the applicability of peak reactivity methods currently used for spent fuel pools to spent fuel storage and transportation casks and the validation of reactivity (k_{eff}) calculations and predicted spent fuel compositions. The second phase focuses on extending BUC beyond peak reactivity. This paper documents work performed to date investigating some aspects of extended BUC. (The technical basis for application of peak reactivity methods to BWR fuel in storage and transportation systems is presented in a companion paper.) Two reactor operating parameters are being evaluated to establish an adequate basis for extended BWR BUC: (1) the effect of axial void profile and (2) the effect of control blade utilization during operation. A detailed analysis of core simulator data for one cycle of a modern operating BWR plant was performed to determine the range of void profiles and the variability of the profile experienced during irradiation. Although a single cycle does not provide complete data, the data obtained are sufficient to determine the primary effects and to identify conservative modeling approaches. These data were used in a study of the effect of axial void profile. The first stage of the study was determination of the necessary moderator density temporal fidelity in depletion modeling. After the required temporal fidelity was established, multiple void profiles were used to examine the effect on cask reactivity. The results of these studies are being used to develop recommendations for conservatively modeling the void profile effects for BWR depletion calculations. The second operational parameter studied was control blade history. Control blades are inserted in
NASA Astrophysics Data System (ADS)
Marchand, P.; Masson, J.; Chabrier, G.; Hennebelle, P.; Commerçon, B.; Vaytet, N.
2016-07-01
We develop a detailed chemical network relevant to calculate the conditions that are characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of potassium, sodium, and hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to nH = 1012 cm-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process. The multi-dimensional multi-species equilibrium abundance table and a copy of the code are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/592/A18
NASA Astrophysics Data System (ADS)
Wenzel, Jan; Holzer, Andre; Wormit, Michael; Dreuw, Andreas
2015-06-01
The extended second order algebraic-diagrammatic construction (ADC(2)-x) scheme for the polarization operator in combination with core-valence separation (CVS) approximation is well known to be a powerful quantum chemical method for the calculation of core-excited states and the description of X-ray absorption spectra. For the first time, the implementation and results of the third order approach CVS-ADC(3) are reported. Therefore, the CVS approximation has been applied to the ADC(3) working equations and the resulting terms have been implemented efficiently in the adcman program. By treating the α and β spins separately from each other, the unrestricted variant CVS-UADC(3) for the treatment of open-shell systems has been implemented as well. The performance and accuracy of the CVS-ADC(3) method are demonstrated with respect to a set of small and middle-sized organic molecules. Therefore, the results obtained at the CVS-ADC(3) level are compared with CVS-ADC(2)-x values as well as experimental data by calculating complete basis set limits. The influence of basis sets is further investigated by employing a large set of different basis sets. Besides the accuracy of core-excitation energies and oscillator strengths, the importance of cartesian basis functions and the treatment of orbital relaxation effects are analyzed in this work as well as computational timings. It turns out that at the CVS-ADC(3) level, the results are not further improved compared to CVS-ADC(2)-x and experimental data, because the fortuitous error compensation inherent in the CVS-ADC(2)-x approach is broken. While CVS-ADC(3) overestimates the core excitation energies on average by 0.61% ± 0.31%, CVS-ADC(2)-x provides an averaged underestimation of -0.22% ± 0.12%. Eventually, the best agreement with experiments can be achieved using the CVS-ADC(2)-x method in combination with a diffuse cartesian basis set at least at the triple-ζ level.
Wenzel, Jan Holzer, Andre; Wormit, Michael; Dreuw, Andreas
2015-06-07
The extended second order algebraic-diagrammatic construction (ADC(2)-x) scheme for the polarization operator in combination with core-valence separation (CVS) approximation is well known to be a powerful quantum chemical method for the calculation of core-excited states and the description of X-ray absorption spectra. For the first time, the implementation and results of the third order approach CVS-ADC(3) are reported. Therefore, the CVS approximation has been applied to the ADC(3) working equations and the resulting terms have been implemented efficiently in the adcman program. By treating the α and β spins separately from each other, the unrestricted variant CVS-UADC(3) for the treatment of open-shell systems has been implemented as well. The performance and accuracy of the CVS-ADC(3) method are demonstrated with respect to a set of small and middle-sized organic molecules. Therefore, the results obtained at the CVS-ADC(3) level are compared with CVS-ADC(2)-x values as well as experimental data by calculating complete basis set limits. The influence of basis sets is further investigated by employing a large set of different basis sets. Besides the accuracy of core-excitation energies and oscillator strengths, the importance of cartesian basis functions and the treatment of orbital relaxation effects are analyzed in this work as well as computational timings. It turns out that at the CVS-ADC(3) level, the results are not further improved compared to CVS-ADC(2)-x and experimental data, because the fortuitous error compensation inherent in the CVS-ADC(2)-x approach is broken. While CVS-ADC(3) overestimates the core excitation energies on average by 0.61% ± 0.31%, CVS-ADC(2)-x provides an averaged underestimation of −0.22% ± 0.12%. Eventually, the best agreement with experiments can be achieved using the CVS-ADC(2)-x method in combination with a diffuse cartesian basis set at least at the triple-ζ level.
Wenzel, Jan; Holzer, Andre; Wormit, Michael; Dreuw, Andreas
2015-06-01
The extended second order algebraic-diagrammatic construction (ADC(2)-x) scheme for the polarization operator in combination with core-valence separation (CVS) approximation is well known to be a powerful quantum chemical method for the calculation of core-excited states and the description of X-ray absorption spectra. For the first time, the implementation and results of the third order approach CVS-ADC(3) are reported. Therefore, the CVS approximation has been applied to the ADC(3) working equations and the resulting terms have been implemented efficiently in the adcman program. By treating the α and β spins separately from each other, the unrestricted variant CVS-UADC(3) for the treatment of open-shell systems has been implemented as well. The performance and accuracy of the CVS-ADC(3) method are demonstrated with respect to a set of small and middle-sized organic molecules. Therefore, the results obtained at the CVS-ADC(3) level are compared with CVS-ADC(2)-x values as well as experimental data by calculating complete basis set limits. The influence of basis sets is further investigated by employing a large set of different basis sets. Besides the accuracy of core-excitation energies and oscillator strengths, the importance of cartesian basis functions and the treatment of orbital relaxation effects are analyzed in this work as well as computational timings. It turns out that at the CVS-ADC(3) level, the results are not further improved compared to CVS-ADC(2)-x and experimental data, because the fortuitous error compensation inherent in the CVS-ADC(2)-x approach is broken. While CVS-ADC(3) overestimates the core excitation energies on average by 0.61% ± 0.31%, CVS-ADC(2)-x provides an averaged underestimation of -0.22% ± 0.12%. Eventually, the best agreement with experiments can be achieved using the CVS-ADC(2)-x method in combination with a diffuse cartesian basis set at least at the triple-ζ level. PMID:26049476
ISOTOPIC MODEL FOR COMMERCIAL SNF BURNUP CREDIT
A.H. Wells
2004-11-17
The purpose of this report is to demonstrate a process for selecting bounding depletion parameters, show that they are conservative for pressurized water reactor (PWR) and boiling water reactor (BWR) spent nuclear fuel (SNF), and establish the range of burnup for which the parameters are conservative. The general range of applicability is for commercial light water reactor (LWR) SNF with initial enrichments between 2.0 and 5.0 weight percent {sup 235}U and burnups between 10 and 50 gigawatt-day per metric ton of uranium (GWd/MTU).
DeHart, M.D.
1996-05-01
Spent fuel transportation and storage cask designs based on a burnup credit approach must consider issues that are not relevant in casks designed under a fresh-fuel loading assumption. For example, the spent fuel composition must be adequately characterized and the criticality analysis model can be complicated by the need to consider axial burnup variations. Parametric analyses are needed to characterize the importance of fuel assembly and fuel cycle parameters on spent fuel composition and reactivity. Numerical models must be evaluated to determine the sensitivity of criticality safety calculations to modeling assumptions. The purpose of this report is to describe analyses and evaluations performed in order to demonstrate the effect physical parameters and modeling assumptions have on the criticality analysis of spent fuel. The analyses in this report include determination and ranking of the most important actinides and fission products; study of the effect of various depletion scenarios on subsequent criticality calculations; establishment of trends in neutron multiplication as a function of fuel enrichment, burnup, cooling time- and a parametric and modeling evaluation of three-dimensional effects (e.g., axially varying burnup and temperature/density effects) in a conceptual cask design. The sensitivity and parametric evaluations were performed with the consideration of two different burnup credit approaches: (1) only actinides in the fuel are considered in the criticality analysis, and (2) both actinides and fission products are considered. Calculations described in this report were performed using the criticality and depletion sequences available in the SCALE code system and the SCALE 27-group burnup library. Although the results described herein do not constitute a validation of SCALE for use in spent fuel analysis, independent validation efforts have been completed and are described in other reports.
Tikhomirov, A. V.; Ponomarenko, G. L.
2012-07-01
An additional verification of bundled software (BS) SAPFIR-95 and amp;RC [1] and code KORSAR/GP [2] was performed. Both software products were developed in A.P. Alexandrov NITI and certified by ROSTEKHNADZOR of RF for numeric simulation of stationary, transitional and emergency conditions of VVER reactors. A benchmark model for neutronics calculations was created within the limits of this work. The cold subcritical state of VVER - 1000 reactor stationary fuelling was simulated on the basis of FA with an increased height of the fuel column (TVS-2M) considering detailed presentation of radial and front neutron reflectors. A case of passing of pure condensate slug through the core in initially deep subcritical state during start of the first RCP set after refueling was considered as an examined condition of reactor operation. A relatively small size of the slug, its spatial position near the reflectors (lower and lateral), as well as failure of the inserted control rods of the control and protection system (CPS CR) to reach the lower limit of the fuel column stipulate for methodical complexity of a correct calculation of the neutron multiplication constant (K{sub eff}) using engineering codes. Code RC was used as a test program in the process of reactor calculated 3-D modeling. Code MCNP5 [3] was used as the precision program, which solves the equation of neutrons transfer by Monte-Carlo method and which was developed in the US (Los-Alamos). As a result of comparative calculations dependency of K{sub eff} on two parameters was evaluated - boron acid concentration (Cb) and CPS CR position. Reactivity effect was evaluated, which is implemented as a result of failure of all CPS control rods to reach the lower fuel limit calculated using the engineering codes mentioned above. (authors)
Calculating recombination rates and biexciton binding/antibinding in core-shell dots and nano-rods
NASA Astrophysics Data System (ADS)
Shumway, John
2014-03-01
Predicting radiative lifetimes and photoluminescence (PL) emission energies from electron-hole recombination in nano structures is complicated by correlation. Quantum correlations--particularly the attraction between the recombining electron and hole--reduce the PL emission energy but also modify the wave functions, enhancing recombination rates. Interactions with spectator particles can also affect energies and lifetimes, though sometimes the sign of these changes is non-intuitive. Path-integral quantum Monte Carlo (PI-QMC) is a wave-function free computational quantum approach that can easily handle interactions between several electrons and holes in a nanostructure. We present an application to core-shell dots and nano-rods, where proper treatment of correlation is necessary to understand the binding/antibinding transition in the biexciton. The imaginary-time paths provide further insights into the properties of the electron-hole states. We show how changing the topology of the paths can be used to calculate recombination rates and give insights into the recombination process. Fluctuations in the paths are used to calculate responses to electric and magnetic fields. These calculations are performed with the open source pi-qmc code available on GitHub and as a community resource on the nanoHUB. Work supported by NSF OCI 1148502.
Arinkin, F.; Chakrov, P.; Chekushina, L.; Gizatulin,, Sh.; Koltochnik, S.; Hanan, N.; Garner, P.; Nuclear Engineering Division; Kazakhstan Ministry of Energy and Mineral Resources
2010-03-01
In 2010 life test of three LEU (19.7%) lead test assemblies (LTA) is expected in the existing WWR-K reactor core with regular WWR-C-type fuel assemblies and a smaller core with a beryllium insert. Preliminary analysis of test safety is to be carried out. It implies reconstruction of the reactor core history for last three years, including burnup calculation for each regular fuel assembly (FA), as well as calculation of characteristics of the test core. For the planned configuration of the test core a number of characteristics have been calculated. The obtained data will be used as input for calculations on LTA test core steady-state thermal hydraulics and on transient analysis.
Model for evolution of grain size in the rim region of high burnup UO2 fuel
NASA Astrophysics Data System (ADS)
Xiao, Hongxing; Long, Chongsheng; Chen, Hongsheng
2016-04-01
The restructuring process of the high burnup structure (HBS) formation in UO2 fuel results in sub-micron size grains that accelerate the fission gas swelling, which will raise some concern over the safety of extended the nuclear fuel operation life in the reactor. A mechanistic and engineering model for evolution of grain size in the rim region of high burnup UO2 fuel based on the experimental observations of the HBS in the literature is presented. The model takes into account dislocations evolution under irradiation and the grain subdivision occur successively at increasing local burnup. It is assumed that the original driving force for subdivision of grain in the HBS of UO2 fuel is the production and accumulation of dislocation loops during irradiation. The dislocation loops can also be annealed through thermal diffusion when the temperature is high enough. The capability of this model is validated by the comparison with the experimental data of temperature threshold of subdivision, dislocation density and sub-grain size as a function of local burnup. It is shown that the calculated results of the dislocation density and subdivided grain size as a function of local burnup are in good agreement with the experimental results.
NFCSim: A Dynamic Fuel Burnup and Fuel Cycle Simulation Tool
Schneider, Erich A.; Bathke, Charles G.; James, Michael R.
2005-07-15
NFCSim is an event-driven, time-dependent simulation code modeling the flow of materials through the nuclear fuel cycle. NFCSim tracks mass flow at the level of discrete reactor fuel charges/discharges and logs the history of nuclear material as it progresses through a detailed series of processes and facilities, generating life-cycle material balances for any number of reactors. NFCSim is an ideal tool for analysis - of the economics, sustainability, or proliferation resistance - of nonequilibrium, interacting, or evolving reactor fleets. The software couples with a criticality and burnup engine, LACE (Los Alamos Criticality Engine). LACE implements a piecewise-linear, reactor-specific reactivity model for its criticality calculations. This model constructs fluence-dependent reactivity traces for any facility; it is designed to address nuclear economies in which either a steady state is never obtained or is a poor approximation. LACE operates in transient and equilibrium fuel management regimes at the refueling batch level, derives reactor- and cycle-dependent initial fuel compositions, and invokes ORIGEN2.x to carry out burnup calculations.
Choi, Sunghwan; Kwon, Oh-Kyoung; Kim, Jaewook; Kim, Woo Youn
2016-09-15
We investigated the performance of heterogeneous computing with graphics processing units (GPUs) and many integrated core (MIC) with 20 CPU cores (20×CPU). As a practical example toward large scale electronic structure calculations using grid-based methods, we evaluated the Hartree potentials of silver nanoparticles with various sizes (3.1, 3.7, 4.9, 6.1, and 6.9 nm) via a direct integral method supported by the sinc basis set. The so-called work stealing scheduler was used for efficient heterogeneous computing via the balanced dynamic distribution of workloads between all processors on a given architecture without any prior information on their individual performances. 20×CPU + 1GPU was up to ∼1.5 and ∼3.1 times faster than 1GPU and 20×CPU, respectively. 20×CPU + 2GPU was ∼4.3 times faster than 20×CPU. The performance enhancement by CPU + MIC was considerably lower than expected because of the large initialization overhead of MIC, although its theoretical performance is similar with that of CPU + GPU. © 2016 Wiley Periodicals, Inc. PMID:27431905
Wang, T K; Peir, J J
2000-01-01
The purpose of this work is to establish a method for evaluating the burn-up values of the rod-type TRIGA spent fuel by using gamma-ray spectrometry of the short-lived fission products 97Zr/97Nb, 132I, and 140La. Fuel irradiation history is not needed in this method. Short-lived fission-product activities were established by reirradiating the spent fuels in a nuclear reactor. Based on the measured activities, 235U burn-up values can be deduced by iterative calculations. The complication caused by 239Pu production and fission is also discussed in detail. The burn-up values obtained by this method are in good agreement with those deduced from the conventional method based on long-lived fission products 137Cs, 134Cs/137Cs ratio and 106Ru/137Cs ratio. PMID:10670930
Sensitivity and Uncertainty Analysis to Burn-up Estimates on ADS Using ACAB Code
Cabellos, O; Sanz, J; Rodriguez, A; Gonzalez, E; Embid, M; Alvarez, F; Reyes, S
2005-02-11
Within the scope of the Accelerator Driven System (ADS) concept for nuclear waste management applications, the burnup uncertainty estimates due to uncertainty in the activation cross sections (XSs) are important regarding both the safety and the efficiency of the waste burning process. We have applied both sensitivity analysis and Monte Carlo methodology to actinides burnup calculations in a lead-bismuth cooled subcritical ADS. The sensitivity analysis is used to identify the reaction XSs and the dominant chains that contribute most significantly to the uncertainty. The Monte Carlo methodology gives the burnup uncertainty estimates due to the synergetic/global effect of the complete set of XS uncertainties. These uncertainty estimates are valuable to assess the need of any experimental or systematic reevaluation of some uncertainty XSs for ADS.
Sensitivity and Uncertainty Analysis to Burnup Estimates on ADS using the ACAB Code
Cabellos, O.; Sanz, J.; Rodriguez, A.; Gonzalez, E.; Embid, M.; Alvarez, F.; Reyes, S.
2005-05-24
Within the scope of the Accelerator Driven System (ADS) concept for nuclear waste management applications, the burnup uncertainty estimates due to uncertainty in the activation cross sections (XSs) are important regarding both the safety and the efficiency of the waste burning process. We have applied both sensitivity analysis and Monte Carlo methodology to actinides burnup calculations in a lead-bismuth cooled subcritical ADS. The sensitivity analysis is used to identify the reaction XSs and the dominant chains that contribute most significantly to the uncertainty. The Monte Carlo methodology gives the burnup uncertainty estimates due to the synergetic/global effect of the complete set of XS uncertainties. These uncertainty estimates are valuable to assess the need of any experimental or systematic re-evaluation of some uncertainty XSs for ADS.
Stout, R.B.; Merckx, K.R.; Holm, J.S.
1981-01-01
This study calculates the reduced uranium requirements and the economic incentives for increasing the burnup of current design LWR fuels from the current range of 25 to 35 MWD/Kg to a range of 45 to 55 MWD/Kg. The changes in fuel management strategies which may be required to accommodate these high burnup fuels and longer fuel cycles are discussed. The material behavior problems which may present obstacles to achieving high burnup or to license fuel are identified and discussed. These problems are presented in terms of integral fuel response and the informational needs for commercial and licensing acceptance. Research and development programs are outlined which are aimed at achieving a licensing position and commercial acceptance of high burnup fuels.
Antineutrinos for Reactor Safeguards: Effect of Fuel Loading and Burnup on the Signal
NASA Astrophysics Data System (ADS)
Erickson, Anna; Bernstein, Adam; Bowden, Nathaniel
2014-02-01
Various types of nuclear reactor related information, including relative power level and fuel evolution parameters, can be inferred remotely using antineutrino detectors. We show that it is possible to verify assembly-level burnup using information derived from an antineutrino detector if the nominal reactor fuel loading is known. Alternatively, if the core power is measured using an independent method, for example, a thermal hydraulic element, and the nominal core behavior is known, the antineutrino detector has a capability to determine previously unknown MOX loading in the core.
A Modal Expansion Equilibrium Cycle Perturbation Method for Optimizing High Burnup Fast Reactors
NASA Astrophysics Data System (ADS)
Touran, Nicholas W.
This dissertation develops a simulation tool capable of optimizing advanced nuclear reactors considering the multiobjective nature of their design. An Enhanced Equilibrium Cycle (EEC) method based on the classic equilibrium method is developed to evaluate the response of the equilibrium cycle to changes in the core design. Advances are made in the consideration of burnup-dependent cross sections and dynamic fuel performance (fission gas release, fuel growth, and bond squeeze-out) to allow accuracy in high-burnup reactors such as the Traveling Wave Reactor. EEC is accelerated for design changes near a reference state through a new modal expansion perturbation method that expands arbitrary flux perturbations on a basis of λ-eigenmodes. A code is developed to solve the 3-D, multigroup diffusion equation with an Arnoldi-based solver that determines hundreds of the reference flux harmonics and later uses these harmonics to determine expansion coefficients required to approximate the perturbed flux. The harmonics are only required for the reference state, and many substantial and localized perturbations from this state are shown to be well-approximated with efficient expressions after the reference calculation is performed. The modal expansion method is coupled to EEC to produce the later-in-time response of each design perturbation. Because the code determines the perturbed flux explicitly, a wide variety of core performance metrics may be monitored by working within a recently-developed data management system called the ARMI. Through ARMI, the response of each design perturbation may be evaluated not only for the flux and reactivity, but also for reactivity coefficients, thermal hydraulics parameters, economics, and transient performance. Considering the parameters available, an automated optimization framework is designed and implemented. A non-parametric surrogate model using the Alternating Conditional Expectation (ACE) algorithm is trained with many design
An Automated, Multi-Step Monte Carlo Burnup Code System.
TRELLUE, HOLLY R.
2003-07-14
Version 02 MONTEBURNS Version 2 calculates coupled neutronic/isotopic results for nuclear systems and produces a large number of criticality and burnup results based on various material feed/removal specifications, power(s), and time intervals. MONTEBURNS is a fully automated tool that links the LANL MCNP Monte Carlo transport code with a radioactive decay and burnup code. Highlights on changes to Version 2 are listed in the transmittal letter. Along with other minor improvements in MONTEBURNS Version 2, the option was added to use CINDER90 instead of ORIGEN2 as the depletion/decay part of the system. CINDER90 is a multi-group depletion code developed at LANL and is not currently available from RSICC. This MONTEBURNS release was tested with various combinations of CCC-715/MCNPX 2.4.0, CCC-710/MCNP5, CCC-700/MCNP4C, CCC-371/ORIGEN2.2, ORIGEN2.1 and CINDER90. Perl is required software and is not included in this distribution. MCNP, ORIGEN2, and CINDER90 are not included.
An Automated, Multi-Step Monte Carlo Burnup Code System.
Energy Science and Technology Software Center (ESTSC)
2003-07-14
Version 02 MONTEBURNS Version 2 calculates coupled neutronic/isotopic results for nuclear systems and produces a large number of criticality and burnup results based on various material feed/removal specifications, power(s), and time intervals. MONTEBURNS is a fully automated tool that links the LANL MCNP Monte Carlo transport code with a radioactive decay and burnup code. Highlights on changes to Version 2 are listed in the transmittal letter. Along with other minor improvements in MONTEBURNS Version 2,more » the option was added to use CINDER90 instead of ORIGEN2 as the depletion/decay part of the system. CINDER90 is a multi-group depletion code developed at LANL and is not currently available from RSICC. This MONTEBURNS release was tested with various combinations of CCC-715/MCNPX 2.4.0, CCC-710/MCNP5, CCC-700/MCNP4C, CCC-371/ORIGEN2.2, ORIGEN2.1 and CINDER90. Perl is required software and is not included in this distribution. MCNP, ORIGEN2, and CINDER90 are not included.« less
Extended Burnup Credit for BWR Spent Nuclear Fuel in Storage and Transportation Systems
Ade, Brian J; Bowman, Stephen M; Gauld, Ian C; Ilas, Germina; Martinez, J. S.
2015-01-01
[Full Text] Oak Ridge National Laboratory and the United States Nuclear Regulatory Commission have initiated a multiyear project to investigate the application of burnup credit (BUC) for boiling-water reactor (BWR) fuel in storage and transportation casks. This project includes two phases. The first phase investigates the applicability of peak reactivity methods currently used for spent fuel pools to spent fuel storage and transportation casks and the validation of reactivity (k_{eff}) calculations and depleted fuel compositions. The second phase focuses on extending BUC beyond peak reactivity. This paper documents work performed to date, investigating some aspects of extended BUC, and it also describes the plan to complete the evaluations. The technical basis for application of peak reactivity methods to BWR fuel in storage and transportation systems is presented in a companion paper. Two reactor operating parameters are being evaluated to establish an adequate basis for extended BWR BUC, including investigation of the axial void profile effect and the effect of control blade utilization during operation. A detailed analysis of core simulator data for one cycle of an operating BWR plant was performed to determine the range of void profiles and the variability of the profile experienced during irradiation. While a single cycle does not provide complete data, the data obtained are sufficient to use to determine the primary effects and identify conservative modeling approaches. Using data resulting from a single cycle, the axial void profile is studied by first determining the temporal fidelity necessary in depletion modeling, and then using multiple void profiles to examine the effect of the void profile on cask reactivity. The results of these studies are being used to develop recommendations for conservatively modeling the void profile effects for BWR depletion calculations. The second operational parameter studied is control blade exposure. Control blades
Revised Burnup Code System SWAT: Description and Validation Using Postirradiation Examination Data
Suyama, Kenya; Mochizuki, Hiroki; Kiyosumi, Takehide
2002-05-15
The burnup code system Step-Wise Burnup Analysis Code System (SWAT) is revised for use in a burnup credit analysis. An important feature of the revised SWAT is that its functions are achieved by calling validated neutronics codes without any changes to the original codes. This feature is realized with a system function of the operating system, which allows the revised SWAT to be independent of the development status of each code.A package of the revised SWAT contains the latest libraries based on JENDL-3.2 and the second version of the JNDC FP library. These libraries allow us to analyze burnup problems, such as an analysis of postirradiation examination (PIE), using the latest evaluated data of not only cross sections but also fission yield and decay constants.Another function of the revised SWAT is a library generator for the ORIGEN2 code, which is one of the most reliable burnup codes. ORIGEN2 users can obtain almost the same results with the revised SWAT using the library prepared by this function.The validation of the revised SWAT is conducted by calculation of the Organization for Economic Cooperation and Development/Nuclear Energy Agency burnup credit criticality safety benchmark Phase I-B and analyses of PIE data for spent fuel from Takahama Unit 3. The analysis of PIE data shows that the revised SWAT can predict the isotopic composition of main uranium and plutonium with a deviation of 5% from experimental results taken from UO{sub 2} fuels of 17 x 17 fuel assemblies. Many results of fission products including samarium are within a deviation of 10%. This means that the revised SWAT has high reliability to predict the isotopic composition for pressurized water reactor spent fuel.
Parametric Study of the Effect of Burnable Poison Rods for PWR Burnup Credit
Wagner, J.C.
2001-09-28
The Interim Staff Guidance on burnup credit (ISG-8) issued by the United States Nuclear Regulatory Commission's (U.S. NRC) Spent Fuel Project Office recommends restricting the use of burnup credit to assemblies that have not used burnable absorbers. This recommended restriction eliminates a large portion of the currently discharged spent fuel assemblies from cask loading, and thus severely limits the practical usefulness of burnup credit. In the absence of readily available information on burnable poison rod (BPR) design specifications and usage in U.S. pressurized-water-reactors (PWRs), and the subsequent reactivity effect of BPR exposure on discharged spent nuclear fuel (SNF), NRC staff has indicated a need for additional information in these areas. In response, this report presents a parametric study of the effect of BPR exposure on the reactivity of SNF for various BPR designs, fuel enrichments, and exposure conditions, and documents BPR design specifications. Trends in the reactivity effects of BPRs are established with infinite pin-cell and assembly array calculations with the SCALE and HELIOS code packages, respectively. Subsequently, the reactivity effects of BPRs for typical initial enrichment and burnup combinations are quantified based on three-dimensional (3-D) KENO V.a Monte Carlo calculations with a realistic rail-type cask designed for burnup credit. The calculations demonstrate that the positive reactivity effect due to BPR exposure increases nearly linearly with burnup and is dependent on the number, poison loading, and design of the BPRs and the initial fuel enrichment. Expected typical reactivity increases, based on one-cycle BPR exposure, were found to be less than 1% {Delta}k. Based on the presented analysis, guidance is offered on an appropriate approach for calculating bounding SNF isotopic data for assemblies exposed to BPRs. Although the analyses do not address the issue of validation of depletion methods for assembly designs with BPRs
High burnup effects in WWER fuel rods
Smirnov, V.; Smirnov, A.
1996-03-01
Since 1987 at the Research Institute of Atomic Reactors, the examinations of the WWER spent fuel assemblies has been carried out. These investigations are aimed to gain information on WWER spent fuel conditions in order to validate the fuel assemblies use during the 3 and 4 year fuel cycle in the WWER-440 and WWER-1000 units. At present time, the aim is to reach an average fuel burnup of 55 MWd/kgU. According to this aim, a new investigation program on the WWER spent fuel elements is started. The main objectives of this program are to study the high burnup effects and their influence on the WWER fuel properties. This paper presented the main statistical values of the WWER-440 and WWER-1000 reactors` fuel assemblies and their fragment parameters. Average burnup of fuel in the investigated fuel assemblies was in the range of 13 to 49.7 MWd/kgU. In this case, the numer of fuel cycles was from 1 to 4 during operation of the fuel assemblies.
Effect of Control Blade History, and Axial Coolant Density and Burnup Profiles on BWR Burnup Credit
Marshall, William BJ J
2016-01-01
A technical basis for peak reactivity boiling water reactor (BWR) burnup credit (BUC) methods was recently generated, and the technical basis for extended BWR BUC is now being developed. In this paper, a number of effects related to extended BWR BUC are analyzed, including three major operational effects in BWRs: the coolant density axial distribution, the use of control blades during operation, and the axial burnup profile. Specifically, uniform axial moderator density profiles are analyzed and compared to previous results and an additional temporal fidelity study combing moderator density profiles for three different fuel assemblies is presented. Realistic control blade histories and cask criticality results are compared to previously generated constructed control blade histories. Finally, a preliminary study of the axial burnup profile is provided.
Lashkari, A; Khalafi, H; Kazeminejad, H
2013-05-01
In this work, kinetic parameters of Tehran research reactor (TRR) mixed cores have been calculated. The mixed core configurations are made by replacement of the low enriched uranium control fuel elements with highly enriched uranium control fuel elements in the reference core. The MTR_PC package, a nuclear reactor analysis tool, is used to perform the analysis. Simulations were carried out to compute effective delayed neutron fraction and prompt neutron lifetime. Calculation of kinetic parameters is necessary for reactivity and power excursion transient analysis. The results of this research show that effective delayed neutron fraction decreases and prompt neutron lifetime increases with the fuels burn-up. Also, by increasing the number of highly enriched uranium control fuel elements in the reference core, the prompt neutron lifetime increases, but effective delayed neutron fraction does not show any considerable change. PMID:24976672
Assessment of Reactivity Margins and Loading Curves for PWR Burnup Credit Cask Designs
Wagner, J.C.
2002-12-17
This report presents studies to assess reactivity margins and loading curves for pressurized water reactor (PWR) burnup-credit criticality safety evaluations. The studies are based on a generic high-density 32-assembly cask and systematically vary individual calculational (depletion and criticality) assumptions to demonstrate the impact on the predicted effective neutron multiplication factor, k{sub eff}, and burnup-credit loading curves. The purpose of this report is to provide a greater understanding of the importance of input parameter variations and quantify the impact of calculational assumptions on the outcome of a burnup-credit evaluation. This study should provide guidance to regulators and industry on the technical areas where improved information will most enhance the estimation of accurate subcritical margins. Based on these studies, areas where future work may provide the most benefit are identified. The report also includes an evaluation of the degree of burnup credit needed for high-density casks to transport the current spent nuclear fuel inventory. By comparing PWR discharge data to actinide-only based loading curves and determining the number of assemblies that meet the loading criteria, this evaluation finds that additional negative reactivity (through either increased credit for fuel burnup or cask design/utilization modifications) is necessary to accommodate the majority of current spent fuel assemblies in high-capacity casks. Assemblies that are not acceptable for loading in the prototypic high-capacity cask may be stored or transported by other means (e.g., lower capacity casks that utilize flux traps and/or increased fixed poison concentrations or high-capacity casks with design/utilization modifications).
Benefits of the delta K of depletion benchmarks for burnup credit validation
Lancaster, D.; Machiels, A.
2012-07-01
Pressurized Water Reactor (PWR) burnup credit validation is demonstrated using the benchmarks for quantifying fuel reactivity decrements, published as 'Benchmarks for Quantifying Fuel Reactivity Depletion Uncertainty,' EPRI Report 1022909 (August 2011). This demonstration uses the depletion module TRITON available in the SCALE 6.1 code system followed by criticality calculations using KENO-Va. The difference between the predicted depletion reactivity and the benchmark's depletion reactivity is a bias for the criticality calculations. The uncertainty in the benchmarks is the depletion reactivity uncertainty. This depletion bias and uncertainty is used with the bias and uncertainty from fresh UO{sub 2} critical experiments to determine the criticality safety limits on the neutron multiplication factor, k{sub eff}. The analysis shows that SCALE 6.1 with the ENDF/B-VII 238-group cross section library supports the use of a depletion bias of only 0.0015 in delta k if cooling is ignored and 0.0025 if cooling is credited. The uncertainty in the depletion bias is 0.0064. Reliance on the ENDF/B V cross section library produces much larger disagreement with the benchmarks. The analysis covers numerous combinations of depletion and criticality options. In all cases, the historical uncertainty of 5% of the delta k of depletion ('Kopp memo') was shown to be conservative for fuel with more than 30 GWD/MTU burnup. Since this historically assumed burnup uncertainty is not a function of burnup, the Kopp memo's recommended bias and uncertainty may be exceeded at low burnups, but its absolute magnitude is small. (authors)
Radiochemical Assays of Irradiated VVER-440 Fuel for Use in Spent Fuel Burnup Credit Activities
Jardine, L J
2005-04-25
The objective of this spent fuel burnup credit work was to study and describe a VVER-440 reactor spent fuel assembly (FA) initial state before irradiation, its operational irradiation history and the resulting radionuclide distribution in the fuel assembly after irradiation. This work includes the following stages: (1) to pick out and select a specific spent (irradiated) FA for examination; (2) to describe the FA initial state before irradiation; (3) to describe the irradiation history, including thermal calculations; (4) to examine the burnup distribution of select radionuclides along the FA height and cross-section; (5) to examine the radionuclide distributions; (6) to determine the Kr-85 release into the plenum; (7) to select and prepare FA rod specimens for destructive examinations; (8) to determine the radionuclide compositions, isotope masses and burnup in the rod specimens; and (9) to analyze, document and process the results. The specific workscope included the destructive assay (DA) of spent fuel assembly rod segments with an {approx}38.5 MWd/KgU burnup from a single VVER-440 fuel assembly from the Novovorenezh reactor in Russia. Based on irradiation history criteria, four rods from the fuel assembly were selected and removed from the assembly for examination. Next, 8 sections were cut from the four rods and sent for destructive analysis of radionuclides by radiochemical analyses. The results were documented in a series of seven reports over a period of {approx}1 1/2 years.
DANDE: a linked code system for core neutronics/depletion analysis
LaBauve, R.J.; England, T.R.; George, D.C.; MacFarlane, R.E.; Wilson, W.B.
1985-06-01
This report describes DANDE - a modular neutronics, depletion code system for reactor analysis. It consists of nuclear data processing, core physics, and fuel depletion modules, and allows one to use diffusion and transport methods interchangeably in core neutronics calculations. This latter capability is especially important in the design of small modular cores. Additional unique features include the capability of updating the nuclear data file during a calculation; a detailed treatment of depletion, burnable poisons as well as fuel; and the ability to make geometric changes such as control rod repositioning and fuel relocation in the course of a calculation. The detailed treatment of reactor fuel burnup, fission-product creation and decay, as well as inventories of higher-order actinides is a necessity when predicting the behavior of reactor fuel under increased burn conditions. The operation of the code system is made clear in this report by following a sample problem.
Malouch, F.
2011-07-01
The control of the temperature in material samples irradiated in a material testing reactor requires the knowledge of the nuclear heating caused by the energy deposition by neutrons and photons interacting in the irradiation device structures. Thus, a neutron-photonic three-dimensional calculation scheme has been developed to evaluate the nuclear heating in experimental devices irradiated in the core of the OSIRIS MTR reactor (CEA/Saclay Center). The aim is to obtain a predictive tool for the nuclear heating estimation in irradiation devices. This calculation scheme is mainly based on the TRIPOLI-4 three-dimensional continuous-energy Monte Carlo transport code, developed by CEA (Saclay Center). An experimental validation has been carried out on the basis of nuclear heating measurements performed in the OSIRIS core. After an overview of the experimental devices irradiated in the OSIRIS reactor, we present the calculation scheme and the first results of the experimental validation. (authors)
NASA Technical Reports Server (NTRS)
Tuma, Margaret L.; Beheim, Glenn
1995-01-01
The effective-index method and Marcatili's technique were utilized independently to calculate the electric field profile of a rib channel waveguide. Using the electric field profile calculated from each method, the theoretical coupling efficiency between a single-mode optical fiber and a rib waveguide was calculated using the overlap integral. Perfect alignment was assumed and the coupling efficiency calculated. The coupling efficiency calculation was then repeated for a range of transverse offsets.
Determination of deuterium-tritium critical burn-up parameter by four temperature theory
NASA Astrophysics Data System (ADS)
Nazirzadeh, M.; Ghasemizad, A.; Khanbabei, B.
2015-12-01
Conditions for thermonuclear burn-up of an equimolar mixture of deuterium-tritium in non-equilibrium plasma have been investigated by four temperature theory. The photon distribution shape significantly affects the nature of thermonuclear burn. In three temperature model, the photon distribution is Planckian but in four temperature theory the photon distribution has a pure Planck form below a certain cut-off energy and then for photon energy above this cut-off energy makes a transition to Bose-Einstein distribution with a finite chemical potential. The objective was to develop four temperature theory in a plasma to calculate the critical burn up parameter which depends upon initial density, the plasma components initial temperatures, and hot spot size. All the obtained results from four temperature theory model are compared with 3 temperature model. It is shown that the values of critical burn-up parameter calculated by four temperature theory are smaller than those of three temperature model.
Determination of deuterium–tritium critical burn-up parameter by four temperature theory
Nazirzadeh, M.; Ghasemizad, A.; Khanbabei, B.
2015-12-15
Conditions for thermonuclear burn-up of an equimolar mixture of deuterium-tritium in non-equilibrium plasma have been investigated by four temperature theory. The photon distribution shape significantly affects the nature of thermonuclear burn. In three temperature model, the photon distribution is Planckian but in four temperature theory the photon distribution has a pure Planck form below a certain cut-off energy and then for photon energy above this cut-off energy makes a transition to Bose-Einstein distribution with a finite chemical potential. The objective was to develop four temperature theory in a plasma to calculate the critical burn up parameter which depends upon initial density, the plasma components initial temperatures, and hot spot size. All the obtained results from four temperature theory model are compared with 3 temperature model. It is shown that the values of critical burn-up parameter calculated by four temperature theory are smaller than those of three temperature model.
Microstructural modeling of thermal conductivity of high burn-up mixed oxide fuel
NASA Astrophysics Data System (ADS)
Teague, Melissa; Tonks, Michael; Novascone, Stephen; Hayes, Steven
2014-01-01
Predicting the thermal conductivity of oxide fuels as a function of burn-up and temperature is fundamental to the efficient and safe operation of nuclear reactors. However, modeling the thermal conductivity of fuel is greatly complicated by the radially inhomogeneous nature of irradiated fuel in both composition and microstructure. In this work, radially and temperature-dependent models for effective thermal conductivity were developed utilizing optical micrographs of high burn-up mixed oxide fuel. The micrographs were employed to create finite element meshes with the OOF2 software. The meshes were then used to calculate the effective thermal conductivity of the microstructures using the BISON [1] fuel performance code. The new thermal conductivity models were used to calculate thermal profiles at end of life for the fuel pellets. These results were compared to thermal conductivity models from the literature, and comparison between the new finite element-based thermal conductivity model and the Duriez-Lucuta model was favorable.
Microstructural Modeling of Thermal Conductivity of High Burn-up Mixed Oxide Fuel
Melissa Teague; Michael Tonks; Stephen Novascone; Steven Hayes
2014-01-01
Predicting the thermal conductivity of oxide fuels as a function of burn-up and temperature is fundamental to the efficient and safe operation of nuclear reactors. However, modeling the thermal conductivity of fuel is greatly complicated by the radially inhomogeneous nature of irradiated fuel in both composition and microstructure. In this work, radially and temperature-dependent models for effective thermal conductivity were developed utilizing optical micrographs of high burn-up mixed oxide fuel. The micrographs were employed to create finite element meshes with the OOF2 software. The meshes were then used to calculate the effective thermal conductivity of the microstructures using the BISON fuel performance code. The new thermal conductivity models were used to calculate thermal profiles at end of life for the fuel pellets. These results were compared to thermal conductivity models from the literature, and comparison between the new finite element-based thermal conductivity model and the Duriez–Lucuta model was favorable.
Gauld, I. C.; Parks, C. V.
2000-12-11
This report has been prepared to review the technical issues important to the prediction of isotopic compositions and source terms for high-burnup, light-water-reactor (LWR) fuel as utilized in the licensing of spent fuel transport and storage systems. The current trend towards higher initial ^{235}U enrichments, more complex assembly designs, and more efficient fuel management schemes has resulted in higher spent fuel burnups than seen in the past. This trend has led to a situation where high-burnup assemblies from operating LWRs now extend beyond the area where available experimental data can be used to validate the computational methods employed to calculate spent fuel inventories and source terms. This report provides a brief review of currently available validation data, including isotopic assays, decay heat measurements, and shielded dose-rate measurements. Potential new sources of experimental data available in the near term are identified. A review of the background issues important to isotopic predictions and some of the perceived technical challenges that high-burnup fuel presents to the current computational methods are discussed. Based on the review, the phenomena that need to be investigated further and the technical issues that require resolution are presented. The methods and data development that may be required to address the possible shortcomings of physics and depletion methods in the high-burnup and high-enrichment regime are also discussed. Finally, a sensitivity analysis methodology is presented. This methodology is currently being investigated at the Oak Ridge National Laboratory as a computational tool to better understand the changing relative significance of the underlying nuclear data in the different enrichment and burnup regimes and to identify the processes that are dominant in the high-burnup regime. The potential application of the sensitivity analysis methodology to help establish a range of applicability for experimental
Mueller, Don; Rearden, Bradley T; Reed, Davis Allan
2010-01-01
One of the challenges associated with implementation of burnup credit is the validation of criticality calculations used in the safety evaluation; in particular the availability and use of applicable critical experiment data. The purpose of the validation is to quantify the relationship between reality and calculated results. Validation and determination of bias and bias uncertainty require the identification of sets of critical experiments that are similar to the criticality safety models. A principal challenge for crediting fission products (FP) in a burnup credit safety evaluation is the limited availability of relevant FP critical experiments for bias and bias uncertainty determination. This paper provides an evaluation of the available critical experiments that include FPs, along with bounding, burnup-dependent estimates of FP biases generated by combining energy dependent sensitivity data for a typical burnup credit application with the nuclear data uncertainty information distributed with SCALE 6. A method for determining separate bias and bias uncertainty values for individual FPs and illustrative results is presented. Finally, a FP bias calculation method based on data adjustment techniques and reactivity sensitivity coefficients calculated with the SCALE sensitivity/uncertainty tools and some typical results is presented. Using the methods described in this paper, the cross-section bias for a representative high-capacity spent fuel cask associated with the ENDF/B-VII nuclear data for 16 most important stable or near stable FPs is predicted to be no greater than 2% of the total worth of the 16 FPs, or less than 0.13 % k/k.
Saavedra, Steven F; Charlton, William S; Solodov, Alexander A; Ehinger, Michael H
2010-01-01
Globally, there exists a long history in reprocessing in evaluation of the shipper/receiver difference (SRD) on spent nuclear fuel (SNF) received and processed. Typically, the declared shipper s values for uranium and plutonium in SNF (based on calculations involving the initial manufacturer s data and reactor operating history) are used as the input quantities to the head-end process of the facility. Problems have been encountered when comparing these values with measured results of the input accountability tank contents. A typical comparison yields a systematic bias indicated as a loss of 5 7 percent of the plutonium (Pu) and approximately 1 percent for the uranium (U). Studies suggest that such deviation can be attributed to the non-linear nature of the axial burnup values of the SNF. Oak Ridge National Laboratory and Texas A&M University are co-investigating the development of a new method, via Nondestructive Assay (NDA) techniques, to improve the accuracy in burnup and Pu content quantification. Two major components have been identified to achieve this objective. The first component calculates a measurement-based burnup profile along the axis of a fuel rod. Gamma-ray data is collected at numerous locations along the axis of the fuel rod using a High Purity Germanium (HPGe) detector designed for a wide range of gamma-ray energies. Using two fission products, 137Cs and 134Cs, the burnup is calculated at each measurement location and a profile created along the axis of the rod based on the individual measurement locations. The second component measures the U/Pu ratio using an HPGe detector configured for relatively low-energy gamma-rays including x-rays. Fluorescence x-rays from U and Pu are measured and compared to the U/Pu ratio determined from a destructive analysis of the sample. This will be used to establish a relationship between the measured and actual values. This relationship will be combined with the burnup analysis results to establish a relationship
Analysis of burnup and isotopic compositions of BWR 9 x 9 UO{sub 2} fuel assemblies
Suzuki, M.; Yamamoto, T.; Ando, Y.; Nakajima, T.
2012-07-01
In order to extend isotopic composition data focusing on fission product nuclides, measurements are progressing using facilities of JAEA for five samples taken from high burnup BWR 9 x 9 UO{sub 2} fuel assemblies. Neutronics analysis with an infinite assembly model was applied to the preliminary measurement data using a continuous-energy Monte Carlo burnup calculation code MVP-BURN with nuclear libraries based on JENDL-3.3 and JENDL-4.0. The burnups of the samples were determined to be 28.0, 39.3, 56.6, 68.1, and 64.0 GWd/t by the Nd-148 method. They were compared with those calculated using node-average irradiation histories of power and in-channel void fractions which were taken from the plant data. The comparison results showed that the deviations of the calculated burnups from the measurements were -4 to 3%. It was confirmed that adopting the nuclear data library based on JENDL-4.0 reduced the deviations of the calculated isotopic compositions from the measurements for {sup 238}Pu, {sup 144}Nd, {sup 145}Nd, {sup 146}Nd, {sup 148}Nd, {sup 134}Cs, {sup 154}Eu, {sup 152}Sm, {sup 154}Gd, and {sup 157}Gd. On the other hand, the effect of the revision in the nuclear. data library on the neutronics analysis was not significant for major U and Pu isotopes. (authors)
Gritzay, Olena; Kalchenko, Oleksandr; Klimova, Nataliya; Razbudey, Volodymyr; Sanzhur, Andriy; Binney, Stephen
2005-05-24
The presented results show our consecutive steps in developing a neutron source with parameters required by Boron Neutron Capture Therapy (BNCT) at the Kyiv Research Reactor (KRR). The main goal of this work was to analyze the influence of installation of different types of uranium converters close to the reactor core on neutron beam characteristics and on level of reactor safety. The general Monte Carlo radiation transport code MCNP, version 4B, has been used for these calculations.
Uncertainties in the effects of burnup and their impact on criticality safety licensing criteria
Carlson, R.W.; Fisher, L.E.
1990-07-13
Current criteria for criticality safety for spent fuel shipping and storage casks are conservative because no credit is permitted for the effects of burnup of the fuel inside the cask. Cask designs that will transport and store large numbers of fuel assemblies (20 or more) must devote a substantial part of their payload to criticality control measures if they are to meet this criteria. The Department of Energy is developing the data necessary to support safety analyses that incorporate the effects of burnup for the next generation of spent fuel shipping casks. The efforts described here are devoted to the development of acceptance criteria that will be the basis for accepting safety analyses. Preliminary estimates of the uncertainties of the effects of burnup have been developed to provide a basis for the consideration of critically safety criteria. The criticality safety margins in a spent fuel shipping or storage cask are dominated by the portions of a fuel assembly that are in low power regions of a reactor core, and the reactor operating conditions are very different from spent fuel storage or transport cask conditions. Consequently, the experience that has been gathered during years of reactor operation does not apply directly to the prediction of criticality safety margins for spent fuel shipping or storage casks. The preliminary estimates of the uncertainties presented in this paper must be refined by both analytical and empirical studies that address both the magnitude of the uncertainties and their interdependence. 9 refs., 5 figs.
Depletion calculations for the McClellan Nuclear Radiation Center.
Klann, R. T.; Newell, D. L.
1997-12-08
Depletion calculations have been performed for the McClellan reactor history from January 1990 through August 1996. A database has been generated for continuing use by operations personnel which contains the isotopic inventory for all fuel elements and fuel-followed control rods maintained at McClellan. The calculations are based on the three-dimensional diffusion theory code REBUS-3 which is available through the Radiation Safety Information Computational Center (RSICC). Burnup-dependent cross-sections were developed at zero power temperatures and full power temperatures using the WIMS code (also available through RSICC). WIMS is based on discretized transport theory to calculate the neutron flux as a function of energy and position in a one-dimensional cell. Based on the initial depletion calculations, a method was developed to allow operations personnel to perform depletion calculations and update the database with a minimal amount of effort. Depletion estimates and calculations can be performed by simply entering the core loading configuration, the position of the control rods at the start and end of cycle, the reactor power level, the duration of the reactor cycle, and the time since the last reactor cycle. The depletion and buildup of isotopes of interest (heavy metal isotopes, erbium isotopes, and fission product poisons) are calculated for all fuel elements and fuel-followed control rods in the MNRC inventory. The reactivity loss from burnup and buildup of fission product poisons and the peak xenon buildup after shutdown are also calculated. The reactivity loss from going from cold zero power to hot full power can also be calculated by using the temperature-dependent, burnup-dependent cross-sections. By calculating all of these reactivity effects, operations personnel are able to estimate the total excess reactivity necessary to run the reactor for the given cycle. This method has also been used to estimate the worth of individual control rods. Using this
Interaction of dislocations in UO2 during high burn-up structure formation
NASA Astrophysics Data System (ADS)
Baranov, V. G.; Lunev, A. V.; Tenishev, A. V.; Khlunov, A. V.
2014-01-01
Dislocation dynamics is used to investigate the distribution of dislocations in oxide nuclear fuel under irradiation using the values of dislocation density from experiments. A model is constructed to account for the effects of irradiation on dislocation movement and for the brittle behavior of the material. Results show that the ground state of interacting dislocations in UO2 during irradiation is a periodic structure with spacing between walls equal to 100-300 nm at experimental dislocation densities. These regions adorned by dislocation walls are called sub-grains and represent the result of polygonization. The threshold of polygonization is shown to depend on the fluctuations of the stress field produced by interaction of many dislocations. These fluctuations reach a critical value when a critical dislocation density is reached (˜4 × 1014 m-2). The calculated value matches experimental data on dislocation density measurement of irradiated uranium dioxide at burn-up corresponding to the formation of high burn-up structure.
Benchmark of Atucha-2 PHWR RELAP5-3D control rod model by Monte Carlo MCNP5 core calculation
Pecchia, M.; D'Auria, F.; Mazzantini, O.
2012-07-01
Atucha-2 is a Siemens-designed PHWR reactor under construction in the Republic of Argentina. Its geometrical complexity and peculiarities require the adoption of advanced Monte Carlo codes for performing realistic neutronic simulations. Therefore core models of Atucha-2 PHWR were developed using MCNP5. In this work a methodology was set up to collect the flux in the hexagonal mesh by which the Atucha-2 core is represented. The scope of this activity is to evaluate the effect of obliquely inserted control rod on neutron flux in order to validate the RELAP5-3D{sup C}/NESTLE three dimensional neutron kinetic coupled thermal-hydraulic model, applied by GRNSPG/UNIPI for performing selected transients of Chapter 15 FSAR of Atucha-2. (authors)
Zhang, D.; Rahnema, F.
2013-07-01
The coarse mesh transport method (COMET) is a highly accurate and efficient computational tool which predicts whole-core neutronics behaviors for heterogeneous reactor cores via a pre-computed eigenvalue-dependent response coefficient (function) library. Recently, a high order perturbation method was developed to significantly improve the efficiency of the library generation method. In that work, the method's accuracy and efficiency was tested in a small PWR benchmark problem. This paper extends the application of the perturbation method to include problems typical of the other water reactor cores such as BWR and CANDU bundles. It is found that the response coefficients predicted by the perturbation method for typical BWR bundles agree very well with those directly computed by the Monte Carlo method. The average and maximum relative errors in the surface-to-surface response coefficients are 0.02%-0.05% and 0.06%-0.25%, respectively. For CANDU bundles, the corresponding quantities are 0.01%-0.05% and 0.04% -0.15%. It is concluded that the perturbation method is highly accurate and efficient with a wide range of applicability. (authors)
Methodology for the Weapons-Grade MOX Fuel Burnup Analysis in the Advanced Test Reactor
G. S. Chang
2005-08-01
A UNIX BASH (Bourne Again SHell) script CMO has been written and validated at the Idaho National Laboratory (INL) to couple the Monte Carlo transport code MCNP with the depletion and buildup code ORIGEN-2 (CMO). The new Monte Carlo burnup analysis methodology in this paper consists of MCNP coupling through CMO with ORIGEN-2, and is therefore called the MCWO. MCWO is a fully automated tool that links the Monte Carlo transport code MCNP with the radioactive decay and burnup code ORIGEN-2. MCWO is capable of handling a large number of fuel burnup and material loading specifications, Advanced Test Reactor (ATR) lobe powers, and irradiation time intervals. MCWO processes user input that specifies the system geometry, initial material compositions, feed/removal specifications, and other code-specific parameters. Calculated results from MCNP, ORIGEN-2, and data process module calculations are output in succession as MCWO executes. The principal function of MCWO is to transfer one-group cross-section and flux values from MCNP to ORIGEN-2, and then transfer the resulting material compositions (after irradiation and/or decay) from ORIGEN-2 back to MCNP in a repeated, cyclic fashion. The basic requirements of MCWO are a working MCNP input file and some additional input parameters; all interaction with ORIGEN-2 as well as other calculations are performed by CMO. This paper presents the MCWO-calculated results for the Reduced Enrichment Research and Test Reactor (RERTR) experiments RERTR-1 and RERTR-2 as well as the Weapons-Grade Mixed Oxide (WG-MOX) fuel testing in ATR. Calculations performed for the WG-MOX test irradiation, which is managed by the Oak Ridge National Laboratory (ORNL), supports the DOE Fissile Materials Disposition Program (FMDP). The MCWO-calculated results are compared with measured data.
NASA Astrophysics Data System (ADS)
Long, Andrew W.; Wong, Bryan M.
2012-09-01
We present a new pseudospectral approach for incorporating many-body, nonlocal exact exchange interactions to understand the formation of electron gases in core-shell nanowires. Our approach is efficiently implemented in the open-source software package PAMELA (Pseudospectral Analysis Method with Exchange & Local Approximations) that can calculate electronic energies, densities, wavefunctions, and band-bending diagrams within a self-consistent Schrödinger-Poisson formalism. The implementation of both local and nonlocal electronic effects using pseudospectral methods is key to PAMELA's efficiency, resulting in significantly reduced computational effort compared to finite-element methods. In contrast to the new nonlocal exchange formalism implemented in this work, we find that the simple, conventional Schrödinger-Poisson approaches commonly used in the literature (1) considerably overestimate the number of occupied electron levels, (2) overdelocalize electrons in nanowires, and (3) significantly underestimate the relative energy separation between electronic subbands. In addition, we perform several calculations in the high-doping regime that show a critical tunneling depth exists in these nanosystems where tunneling from the core-shell interface to the nanowire edge becomes the dominant mechanism of electron gas formation. Finally, in order to present a general-purpose set of tools that both experimentalists and theorists can easily use to predict electron gas formation in core-shell nanowires, we document and provide our efficient and user-friendly PAMELA source code that is freely available at http://alum.mit.edu/www/usagi.
NASA Astrophysics Data System (ADS)
Ohno, M.; Decleva, P.
1993-05-01
The carbon and oxygen 1s core excitation spectra of free CO and NiCO are calculated by ab initio 1h1p/1h1p and 2h2p/2h2p configuration interaction (CI) method using an extended basis set. We employed the ground state as well as core-hole relaxed orbitals. For free CO, we obtain a reasonably good description of the electron energy loss spectroscopy (EELS) spectra. The present interpretation of the spectra agrees with others. For NiCO, we obtain a reasonably good description of the near edge x-ray absorption fine structure (NEXAFS) spectra of the CO/Ni(100) system and that of the electron energy loss spectroscopy (EELS) spectra of the gas phase Ni(CO)4 . We show the existence of the Rydberg-derived additional excited states in the NEXAFS spectra of the chemisorbed molecule and give an interpretation of these states. The disappearance of the giant shake-up satellite in the NEXAFS spectra of the adsorbate is explained in terms of the hindrance of the cooperative core-hole screening mechanism in the π* resonantly excited state. The core-hole screening mechanism in the σ* resonantly excited state is also investigated.
Wagner, J.C.
2001-08-02
This report proposes and documents a computational benchmark problem for the estimation of the additional reactivity margin available in spent nuclear fuel (SNF) from fission products and minor actinides in a burnup-credit storage/transport environment, relative to SNF compositions containing only the major actinides. The benchmark problem/configuration is a generic burnup credit cask designed to hold 32 pressurized water reactor (PWR) assemblies. The purpose of this computational benchmark is to provide a reference configuration for the estimation of the additional reactivity margin, which is encouraged in the U.S. Nuclear Regulatory Commission (NRC) guidance for partial burnup credit (ISG8), and document reference estimations of the additional reactivity margin as a function of initial enrichment, burnup, and cooling time. Consequently, the geometry and material specifications are provided in sufficient detail to enable independent evaluations. Estimates of additional reactivity margin for this reference configuration may be compared to those of similar burnup-credit casks to provide an indication of the validity of design-specific estimates of fission-product margin. The reference solutions were generated with the SAS2H-depletion and CSAS25-criticality sequences of the SCALE 4.4a package. Although the SAS2H and CSAS25 sequences have been extensively validated elsewhere, the reference solutions are not directly or indirectly based on experimental results. Consequently, this computational benchmark cannot be used to satisfy the ANS 8.1 requirements for validation of calculational methods and is not intended to be used to establish biases for burnup credit analyses.
NASA Astrophysics Data System (ADS)
Anoufa, M.; Kiat, J. M.; Kornev, I.; Bogicevic, C.
2013-10-01
In this paper, we want to emphasize the fact that many experimental properties of ceramics can be explained by the existence of a core-shell structure of the grains, particularly at small sizes. In this framework, we have studied BaTiO3 (BT) ceramics constituted of core-shell nanoparticles, nanowires, or nanoplanes by using ab initio derived effective Hamiltonian calculations whose application range is for large values of shell thickness and low values of shell permittivity. Many differences and new features compared to the situation of nanodots are induced by the core-shell structure. For instance, phase sequences are different; there is also a coexistence of vortices found by Naumov, Bellaiche, and Fu [I. I. Naumov, L. Bellaiche, and H. Fu, Nature (London)10.1038/nature03107 432, 737 (2004)] in the case of isolated dots with a homogeneous polarization, a transition from cubic paraelectric phase towards nonpolar rhombohedral phase, anomalies in dielectric permittivity associated with the onset of toroidal moments, etc. Afterwards, we compare these results with those obtained by the Landau theory of core-shell ceramics we have recently published. However, the ab initio calculations fail to capture the physics at small shell thickness and/or high shell permittivity, whereas the Landau theory fails to predict the peculiar properties of the phases in which vortices exist. Therefore, in a tentative way to build a global theory, we have constructed a Landau potential using both the polarization and the toroidal moment as competing order parameters, which allows us to propose a phase diagram, whatever the thickness and permittivity of the shell are.
Mitenkova, E. F.; Novikov, N. V.; Blokhin, A. I.
2012-07-01
The different (U-Pu) fuel compositions are considered for next generation of sodium fast breeder reactors. The considerable discrepancies in axial and radial neutron spectra for hybrid reactor systems compared to the cores with UO{sub 2} fuel cause increasing uncertainty of generating the group nuclear constants in those reactor systems. The calculation results of BFS-62-3A critical assembly which is considered as full-scale model of BN-600 hybrid core with steel reflector specify quite different spectra in local areas. For those systems the MCNP 5 calculations demonstrate significant sensitivity of effective multiplication factor K{sub eff} and spectral indices to nuclear data libraries. For {sup 235}U, {sup 238}U, {sup 239}Pu the results of calculated radial fission rate distributions against the reconstructed ones are analyzed. Comparative analysis of spectral indices, neutron spectra and radial fission rate distributions are performed using the different versions of ENDF/B, JENDL-3.3, JENDL-4, JEFF-3.1.1 libraries and BROND-3 for Fe, Cr isotopes. For analyzing the fission rate sensitivity to the plutonium presence in the fuel {sup 239}Pu is substituted for {sup 235}U (enrichment 90%) in the FA areas containing the plutonium. For {sup 235}U, {sup 238}U, {sup 239}Pu radial fission rate distributions the explanation of pick values discrepancies is based on the group fission constants analyses and possible underestimation of some features at the experimental data recovery method (Westcott factors, temperature dependence). (authors)
NASA Astrophysics Data System (ADS)
Moore, Keith; McLaughlin, Brendan M.; Lane, Ian C.
2016-04-01
BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio 2Σ+ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron effective core-potential and even-tempered augmented polarized core-valence basis sets (aug-pCVnZ-PP, n = 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X2Σ+ state (extrapolated to the CBS limit) is 16 895.12 cm-1 (2.094 eV), which agrees within 0.1% of a revised experimental value of <16 910.6 cm-1, while the calculated re is within 0.03 pm of the experimental result.
NASA Astrophysics Data System (ADS)
Müller, Wolfgang; Meyer, Wilfried
1984-04-01
Extensive all-electron SCF and valence CI calculations are presented for alkali dimer systems with consideration of intershell correlation effects by use of an effective core polarization potential (CPP), which contains only a single adjustable atomic parameter. High accuracy is obtained for the ground-state spectroscopic constants of the studied molecules. The maximum deviations from accurate experimental data are as follows: 1% or 0.03 Å for Re, 2% or 100 cm-1 for De, 0.5% or 1 cm-1 for ωe, and 0.2% or 100 cm-1 for ionization energies. For experimentally uncertain or unknown values reliable predictions can thus be made. The calculated dipole moments for LiK and NaK agree with experiment to within 0.1%, but for LiNa we obtain a deviation of 8% or 0.036 D. An analysis of molecular core polarization contributions reveals the reasons for some systematic defects in previous pseudopotential calculations.
Navarro, J.; Aryaeinejad, R.; Nigg, D.W.
2011-07-01
The goal of this work was to perform a feasibility study and establish measurement techniques to determine the burnup of the Advanced Test Reactor (ATR) fuels at the Idaho National Laboratory (INL). Three different detectors of high purity germanium (HPGe), lanthanum bromide (LaBr{sub 3}), and high pressure xenon (HPXe) in two detection system configurations of below and above the water pool were used in this study. The last two detectors were used for the first time in fuel burnup measurements. The results showed that a better quality spectra can be achieved with the above the water pool configuration. Both short and long cooling time fuels were investigated in order to determine which measurement technique, absolute or fission product ratio, is better suited in each scenario and also to establish what type of detector should be used in each case for the best burnup measurement. The burnup and cooling time calibrations were established using experimental absolute activities or isotopic ratios and ORIGEN burnup calculations. A method was developed to do burnup and cooling time calibrations using fission isotopes activities without the need to know the exact geometry. (authors)
Kishi, Hiroki; Miyazawa, Miki; Matsushima, Naoki; Yamauchi, Jun
2014-02-21
We investigate the X-ray photoelectron spectroscopy (XPS) binding energies of As 3d in Si for various defects in neutral and charged states by first-principles calculation. It is found that the complexes of a substitutional As and a vacancy in charged and neutral states explain the experimentally observed unknown peak very well.
Evaluation of burnup credit for fuel storage analysis -- Experience in Spain
Conde, J.M.; Recio, M.
1995-04-01
Several Spanish light water reactor commercial nuclear power plants are close to maximum spent-fuel pool storage capacity. The utilities are working on the implementation of state-of-the-art methods to increase the storage capacity, including both changes in the pool design (recracking) and the implementation of new analysis approaches with reduced conservation (burnup credit). Burnup credit criticality safety analyses have been approved for two pressurized water reactor plants (four units) and one boiling water reactor (BWR); an other BWR storage analysis is being developed at this moment. The elimination of the ``fresh fuel assumption`` increases the complexity of the criticality analysis to be performed, sometimes putting into question the capability of the analytic tools to properly describe this new situation and increasing the scope of the scenarios to be analyzed. From a regulatory perspective, the reactivity reduction associated with burnup of the fuel can be given credit only if the exposure of each fuel bundle can be known with enough accuracy. Subcriticality of spent-fuel storage depends mainly on the initial fuel enrichment, storage geometry, fuel exposure history, and cooling time. The last two aspects introduced new uncertainties in the criticality analysis that should be quantified in an adequate way. In addition, each and every fuel bundle has its own specific exposure history, so that strong assumptions and simplified calculational schemes have to be developed to undertake the analysis. The Consejo de Seguridad Nuclear (CSN), Spanish regulatory authority on the matter of nuclear safety and radiation protection, plays an active role in the development of analysis methods to support burnup credit, making proposals that may be beneficial in terms of risk and cost while keeping the widest safety margins possible.
Hanson A. L.; Diamond D.
2014-06-30
A plan is being developed for the conversion of the NIST research reactor (NBSR) from high-enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. The LEU fuel may be a monolithic foil (LEUm) of U10Mo (10% molybdenum by weight in an alloy with uranium) or a dispersion of U7Mo in aluminum (LEUd). A previous report provided neutronic calculations for the LEUm fuel and this report presents the neutronics parameters for the LEUd fuel. The neutronics parameters for the LEUd fuel are compared to those previously obtained for the present HEU fuel and the proposed LEUm fuel. The results show no significant differences between the LEUm and the LEUd other than the LEUd fuel requires slightly less uranium than the LEUm fuel due to less molybdenum being present. The calculations include kinetics parameters, reactivity coefficients, reactivity worths of control elements and abnormal configurations, and power distributions under normal operation and with misloaded fuel elements.
PWR cores with silicon carbide cladding
Dobisesky, J. P.; Carpenter, D.; Pilat, E.; Kazimi, M. S.
2012-07-01
The feasibility of using silicon carbide rather than Zircaloy cladding, to reach higher power levels and higher discharge burnups in PWRs has been evaluated. A preliminary fuel design using fuel rods with the same dimensions as in the Westinghouse Robust Fuel Assembly but with fuel pellets having 10 vol% central void has been adopted to mitigate the higher fuel temperatures that occur due to the lower thermal conductivity of the silicon carbide and to the persistence of the open clad-pellet gap over most of the fuel life. With this modified fuel design, it is possible to achieve 18 month cycles that meet present-day operating constraints on peaking factor, boron concentration, reactivity coefficients and shutdown margin, while allowing batch average discharge burnups up to 80 MWD/kgU and peak rod burnups up to 100 MWD/kgU. Power uprates of 10% and possibly 20% also appear feasible. For non-uprated cores, the silicon carbide-clad fuel has a clear advantage that increases with increasing discharge burnup. Even for comparable discharge burnups, there is a savings in enriched uranium. Control rod configuration modifications may be required to meet the shutdown margin criterion for the 20% up-rate. Silicon carbide's ability to sustain higher burnups than Zircaloy also allows the design of a licensable two year cycle with only 96 fresh assemblies, avoiding the enriched uranium penalty incurred with use of larger batch sizes due to their excessive leakage. (authors)
Dieudonne, C.; Dumonteil, E.; Malvagi, F.; Diop, C. M.
2013-07-01
For several years, Monte Carlo burnup/depletion codes have appeared, which couple a Monte Carlo code to simulate the neutron transport to a deterministic method that computes the medium depletion due to the neutron flux. Solving Boltzmann and Bateman equations in such a way allows to track fine 3 dimensional effects and to get rid of multi-group hypotheses done by deterministic solvers. The counterpart is the prohibitive calculation time due to the time-expensive Monte Carlo solver called at each time step. Therefore, great improvements in term of calculation time could be expected if one could get rid of Monte Carlo transport sequences. For example, it may seem interesting to run an initial Monte Carlo simulation only once, for the first time/burnup step, and then to use the concentration perturbation capability of the Monte Carlo code to replace the other time/burnup steps (the different burnup steps are seen like perturbations of the concentrations of the initial burnup step). This paper presents some advantages and limitations of this technique and preliminary results in terms of speed up and figure of merit. Finally, we will detail different possible calculation scheme based on that method. (authors)
Burn-up and neutron economy of accelerator-driven reactor
Takahashi, H.; Yang, W.; An, Y.; Yamazaki, Y.
1997-07-01
It is desirable to have only a small reactivity change in the large burn-up of a solid fuel fast reactor, so that the number of replacements or shuffling of the fuel can be reduced, and plant factor accordingly increased. Also, this reduces the number of control rods needed for the change in burn-up reactivity. In subcritical operation, power controlled by beam power is suggested, but this practice is not as economical as the use of control rods and makes more careful operation of the accelerator is required due to changes in the wake field. In subcritical operation, even a slightly subcritical one, the safety problems associated with a hard neutron spectrum can be alleviated. Neutron leakage from a flattened core, which is needed for operation of the critical fast reactor can be lessen by using the non flat core which has good neutron economy. For generating nuclear energy, it is essential to have a high neutron economy, although breeding the fuel is not welcomed in the present political climate, as is needed for transmuting long lived fission products. In contrast to the breeder, the accelerator driven reactor can separate the energy production from fuel production and processing. Thus, it is suited for non-proliferation of nuclear material by prohibiting the processing and production of fuel in the unrestricted area so this can be only done in international controlled areas which are restricted and remote.
Terry, William Knox; Gougar, Hans D; Ougouag, Abderrafi Mohammed-El-Ami
2002-07-01
A new deterministic method has been developed for the neutronics analysis of a pebble-bed reactor (PBR). The method accounts for the flow of pebbles explicitly and couples the flow to the neutronics. The method allows modeling of once-through cycles as well as cycles in which pebbles are recirculated through the core an arbitrary number of times. This new work is distinguished from older methods by the systematically semi-analytical approach it takes. In particular, whereas older methods use the finite-difference approach (or an equivalent one) for the discretization and the solution of the burnup equation, the present work integrates the relevant differential equation analytically in discrete and complementary sub-domains of the reactor. Like some of the finite-difference codes, the new method obtains the asymptotic fuel-loading pattern directly, without modeling any intermediate loading pattern. This is a significant advantage for the design and optimization of the asymptotic fuel-loading pattern. The new method is capable of modeling directly both the once-through-then-out fuel cycle and the pebble recirculating fuel cycle. Although it currently includes a finite-difference neutronics solver, the new method has been implemented into a modular code that incorporates the framework for the future coupling to an efficient solver such as a nodal method and to modern cross section preparation capabilities. In its current state, the deterministic method presented here is capable of quick and efficient design and optimization calculations for the in-core PBR fuel cycle. The method can also be used as a practical "scoping" tool. It could, for example, be applied to determine the potential of the PBR for resisting nuclear-weapons proliferation and to optimize proliferation-resistant features. However, the purpose of this paper is to show that the method itself is viable. Refinements to the code are under way, with the objective of producing a powerful reactor physics
Code System to Calculate Cross Sections for PWR Fuel Assembly Calculations.
Energy Science and Technology Software Center (ESTSC)
1994-11-15
Version 00 The MARIA System calculates cross sections for PWR fuel assembly calculations. It generates the cross sections library for the diffusion calculations with burnup and feedback effects (CARMEN System, NEA 0649 and RSIC CCC-487) and the k(infinite) and M**2 parameters for the nodal calculations (SIMULA, NEA 0768). MARIA includes three modules. PRELIM generates the input data for the fuel assembly calculation module, for all fuel assembly types in the core and at any conditionmore » of power rate and temperature. WIMS-TRACA is a modified version of the fuel assembly calculation program WIMS-D/4 (NEA 0329 and RSIC CCC-576), which generates the collapsed cross sections versus burn up needed by the CARMEN code (reference cell, boron, xenon, samarium, and light water). POSWIM calculates the transport corrections to the diffusion constant of the absorber materials generated by WIMS-TRACA, to be used directly in the diffusion code when rods or burnable absorber rods are present.« less
NASA Astrophysics Data System (ADS)
Wurster, James
2016-09-01
In this paper, we introduce Nicil: Non-Ideal magnetohydrodynamics Coefficients and Ionisation Library. Nicil is a stand-alone Fortran90 module that calculates the ionisation values and the coefficients of the non-ideal magnetohydrodynamics terms of Ohmic resistivity, the Hall effect, and ambipolar diffusion. The module is fully parameterised such that the user can decide which processes to include and decide upon the values of the free parameters, making this a versatile and customisable code. The module includes both cosmic ray and thermal ionisation; the former includes two ion species and three species of dust grains (positively charged, negatively charged, and neutral), and the latter includes five elements which can be doubly ionised. We demonstrate tests of the module, and then describe how to implement it into an existing numerical code.
Radulescu, Georgeta; Gauld, Ian C.; Ilas, Germina; Wagner, John C.
2014-11-01
This paper describes a depletion code validation approach for criticality safety analysis using burnup credit for actinide and fission product nuclides in spent nuclear fuel (SNF) compositions. The technical basis for determining the uncertainties in the calculated nuclide concentrations is comparison of calculations to available measurements obtained from destructive radiochemical assay of SNF samples. Probability distributions developed for the uncertainties in the calculated nuclide concentrations were applied to the SNF compositions of a criticality safety analysis model by the use of a Monte Carlo uncertainty sampling method to determine bias and bias uncertainty in effective neutron multiplication factor. Application of the Monte Carlo uncertainty sampling approach is demonstrated for representative criticality safety analysis models of pressurized water reactor spent fuel pool storage racks and transportation packages using burnup-dependent nuclide concentrations calculated with SCALE 6.1 and the ENDF/B-VII nuclear data. Furthermore, the validation approach and results support a recent revision of the U.S. Nuclear Regulatory Commission Interim Staff Guidance 8.
Radulescu, Georgeta; Gauld, Ian C.; Ilas, Germina; Wagner, John C.
2014-11-01
This paper describes a depletion code validation approach for criticality safety analysis using burnup credit for actinide and fission product nuclides in spent nuclear fuel (SNF) compositions. The technical basis for determining the uncertainties in the calculated nuclide concentrations is comparison of calculations to available measurements obtained from destructive radiochemical assay of SNF samples. Probability distributions developed for the uncertainties in the calculated nuclide concentrations were applied to the SNF compositions of a criticality safety analysis model by the use of a Monte Carlo uncertainty sampling method to determine bias and bias uncertainty in effective neutron multiplication factor. Application ofmore » the Monte Carlo uncertainty sampling approach is demonstrated for representative criticality safety analysis models of pressurized water reactor spent fuel pool storage racks and transportation packages using burnup-dependent nuclide concentrations calculated with SCALE 6.1 and the ENDF/B-VII nuclear data. Furthermore, the validation approach and results support a recent revision of the U.S. Nuclear Regulatory Commission Interim Staff Guidance 8.« less
Applications of a Monte Carlo whole-core microscopic depletion method
Hutton, J.L.; Butement, A.W.; Watt, S.; Shadbolt, R.D.
1995-12-31
In the WIMS-6 (Ref. 1) reactor physics program scheme a three-dimensional microscopic depletion method has been developed using Monte Carlo fluxes. Together with microscopic cross sections, these give nuclide reaction rates, which are used to solve nuclide depletion equations for each region. An extension of the method, enabling rapid whole-core calculations, has been implemented in the long-established companion code MONK5W. Predictions at successive depletion time steps are based on a calculational route where both geometry and cross sections are accurately represented, providing a reliable and independent approach for benchmarking other methods. Newly developed tracking and storage procedures in MONK5W enable whole core burnup modeling on a desktop computer. Theory and applications are presented in this paper.
Dawahra, S; Khattab, K; Saba, G
2015-07-01
Comparative studies for the conversion of the fuel from HEU to LEU in the Miniature Neutron Source Reactor (MNSR) have been performed using the MCNP4C and GETERA codes. The precise calculations of (135)Xe and (149)Sm concentrations and reactivities were carried out and compared during the MNSR operation time and after shutdown for the existing HEU fuel (UAl4-Al, 90% enriched) and the potential LEU fuels (U3Si2-Al, U3Si-Al, U9Mo-Al, 19.75% enriched and UO2, 12.6% enriched) in this paper using the MCNP4C and GETERA codes. It was found that the (135)Xe and (149)Sm reactivities did not reach their equilibrium reactivities during the daily operating time of the reactor. The (149)Sm reactivities could be neglected compared to (135)Xe reactivities during the reactor operating time and after shutdown. The calculations for the UAl4-Al produced the highest (135)Xe reactivity in all the studied fuel group during the reactor operation (0.39 mk) and after the reactor shutdown (0.735 mk), It followed by U3Si-Al (0.34 mk, 0.653 mk), U3Si2-Al (0.33 mk, 0.634 mk), U9Mo-Al (0.3 mk, 0.568 mk) and UO2 (0.24 mk, 0.448 mk) fuels, respectively. Finally, the results showed that the UO2 was the best candidate for fuel conversion to LEU in the MNSR since it gave the lowest (135)Xe reactivity during the reactor operation and after shutdown. PMID:25816783
Thermal Behavior of Advanced UO{sub 2} Fuel at High Burnup
Muller, E.; Lambert, T.; Silberstein, K.; Therache, B.
2007-07-01
To improve the fuel performance, advanced UO{sub 2} products are developed to reduce significantly Pellet-Cladding Interaction and Fission Gas Release to increase high burnup safety margins on Light Water Reactors. To achieve the expected improvements, doping elements are currently used, to produce large grain viscoplastic UO{sub 2} fuel microstructures. In that scope, AREVA NP is conducting the qualification of a new UO{sub 2} fuel pellet obtained by optimum chromium oxide doping. To assess the fuel thermal performance, especially the fuel conductivity degradation with increasing burnup and also the kinetics of fission gas release under transient operating conditions, an instrumented in-pile experiment, called REMORA, has been developed by the CEA. One segment base irradiated for five cycles in a French EDF commercial PWR ({approx} 62 GWd/tM) was consequently re-instrumented with a fuel centerline thermocouple and an advanced pressure sensor. The design of this specific sensor is based on the counter-pressure principle and avoids any drift phenomenon due to nuclear irradiation. This rodlet was then irradiated in the GRIFFONOS rig of the Osiris experimental reactor at CEA Saclay. This device, located in the periphery of the core, is designed to perform test under conditions close to those prevailing in French PWR reactor. Power variations are carried out by translating the device relatively to the core. Self - powered neutron detectors are positioned in the loop in order to monitor the power the whole time of the irradiation. The re-irradiation of the REMORA experiment consisted of a stepped ramp to power in order to point out a potential degradation of the fuel thermal conductivity with increasing burnup. During the first part of the irradiation, most of the measurements were performed at low power in order to take into account the irradiation effects on UO{sub 2} thermal conductivity at high burnup in low range of temperature. The second part of the irradiation
Carlsson, Johan; Tucek, Kamil; Wider, Hartmut
2006-07-01
This paper concerns two independent safety investigations on critical and sub-critical heavy liquid metal cooled fast reactors using simple flow paths. The first investigation applies to locating the steam generators in the risers instead of the down-comers of a simple flow path designed sub-critical reactor of 600 MW{sub th} power. This was compared to a similar design, but with the steam generators located in the downcomers. The transients investigated were Total-Loss-of-Power and unprotected Loss-Of-Flow. It is shown that this reactor peaks at 1041 K after 29 hours during a Total-Loss-Of-Power accident. The difference between locating the steam generators in the risers and the downcomers is insignificant for this accident type. During an unprotected Loss-Of-Flow accident at full power, the core outlet temperature stabilizes at 1010 K, which is 337 K above nominal outlet temperature. The second investigation concerns a 1426 MW{sub th} critical reactor where the influence of the core height versus the core outlet temperature is studied during an unprotected Loss-Of-Flow and Total-Loss-Of-Power accident. A pancake type core geometry of 1.0 m height and 5.8 m diameter, is compared to a compact core of 2 m height and 4.5 m diameter. Moderators, like BeO and hydrides, and their influence on safety coefficients and burnup swings are also presented. Both cores incinerate transuranics from spent LWR fuel with minor actinide fraction of 5%. We show that LFRs can be designed both to breed and burn transuranics from LWRs. It is shown that the hydrides lead to the most favorable reactivity feedbacks, but the poorest reactivity swing. The computational fluid dynamics code STAR-CD was used for all thermal hydraulic calculations, and the MCNP and MCB for neutronics, and burn-up calculations. (authors)
Cerne, S.P.; Hermann, O.W.; Westfall, R.M.
1987-10-01
This study presents the reactivity loss of spent PWR fuel due to burnup in terms of the infinite lattice multiplications factor, k/sub infinity/. Calculations were performed using the SAS2 and CSAS1 control modules of the SCALE system. The k/sub infinity/ values calculated for all combinations of six enrichments, seven burnups, and five cooling times. The results are presented as a primary function of enrichment in both tabular and graphic form. An equation has been developed to estimate the tabulated values of k/sub infinity/'s by specifying enrichment, cooling time, and burnup. Atom densities for fresh fuel, and spent fuel at cooling times of 2, 10, and 20 years are included. 13 refs., 8 figs., 8 tabs.
Modeling Cladding-Coolant Heat Transfer of High-Burnup Fuel During RIA
Wenfeng Liu; Kazimi, Mujid S.
2006-07-01
This paper describes a model for the cladding-coolant heat transfer of high burnup fuel during a Reactivity Initiated Accident (RIA) which is implemented in the fuel performance code FRAPTRAN 1.2. The minimum stable film boiling temperature, affected by the subcooling and the clad oxidation, is modeled by a modified Henry correlation. This accounts for the effects of thermal properties of the cladding surface on the transient temperature drop during liquid-solid contact. The transition boiling regime is described as the interpolation of the heat flux between two anchor points on the boiling curve: the Critical Heat Flux (CHF) and minimum stable film boiling. The CHF correlation is based on the Zuber hydrodynamic model multiplied by a subcooling factor. Frederking correlation is chosen to model the film boiling regime. The heat conduction through the oxide layer of the cladding surface of high burnup fuel is calculated by solving heat conduction equations with thermal properties of zirconia taken from MATPRO. This model is validated in the FRAPTRAN code for test cases of both high burnup and fresh test fuel rods including the burnup level (0--56 MW d/kg), peak fuel enthalpy deposit (70--190 cal/g), degree of subcooling (0--80 deg. C), and extent of oxidation (0--25 micron). The modified code demonstrates the capability of differentiating between the departure from nucleate boiling (DNB) and none-DNB cases. The predicted peak cladding temperature (PCT) and duration of DNB achieves generally good agreement with the experimental data. It is found that the cladding surface oxidation of high burnup fuel causes an early rewetting of cladding or suppresses DNB due to two factors: 1) Thick zirconia layer may delay the heat conducted to the surface while keeping the surface heat transfer in the most effective nucleate boiling regime. 2) The transient liquid-solid contact resulting from vapor breaking down would cause a lower interface temperature for an oxidized surface
Validation of depletion codes for burnup credit evaluation of LWR assemblies
Ranta-aho, A.
2006-07-01
This paper reports the comparison of the CASMO-4E predictions with the radiochemical assay data from assemblies irradiated in Takahama-3 PWR and Fukushima-Daini-2 BWR, and the most recently reported spent fuel data from the VVER-440 assembly irradiated in Novovoronezh 4. Some of the calculations were repeated with the ABURN burnup code, which is a combination of the MCNP4C Monte Carlo code and the ORIGEN2 depletion code. The cross section libraries applied were based on the ENDF/B-VI and the JEF-2.2 data. (authors)
Shamasundar, B.I.; Fehrenbach, M.E.
1981-05-01
The KUGEL computer code is designed to perform thermal/hydraulic analysis and coated-fuel particle performance calculations for axisymmetric pebble bed reactor (PBR) cores. This computer code was developed as part of a Department of Energy (DOE)-funded study designed to verify the published core performance data on PBRs. The KUGEL code is designed to interface directly with the 2DB code, a two-dimensional neutron diffusion code, to obtain distributions of thermal power, fission rate, fuel burnup, and fast neutron fluence, which are needed for thermal/hydraulic and fuel performance calculations. The code is variably dimensioned so that problem size can be easily varied. An interpolation routine allows variable mesh size to be used between the 2DB output and the two-dimensional thermal/hydraulic calculations.
Armstrong, J.; Hamilton, H.; Hyland, B.
2013-07-01
A program is underway at Atomic Energy of Canada Limited (AECL) to develop a new fuel bundle concept to enable greater burnups for PT-HWR (pressure tube heavy water reactor) cores. One option that AECL is investigating is an internally cooled annular fuel (ICAF) element concept. ICAF contains annular cylindrical pellets with cladding on the inner and outer diameters. Coolant flows along the outside of the element and through the centre. With such a concept, the maximum fuel temperature as a function of linear element rating is significantly reduced compared to conventional, solid-rod type fuel. The preliminary ICAF bundle concept considered in this study contains 24 half-metre long internally cooled annular fuel elements and one non-fuelled centre pin. The introduction of the non-fuelled centre pin reduces the coolant void reactivity (CVR), which is the increase in reactivity that occurs on voiding the coolant in accident scenarios. Lattice cell and full core physics calculations of the preliminary ICAF fuel bundle concept have been performed for medium burnups of approximately 18 GWd/tU using WIMS-AECL and reactor fuel simulation program (RFSP). The results will be used to assist in concept configuration optimization. The effects of radial and axial core power distributions, linear element power ratings, refuelling rates and operational power ramps have been analyzed. The results suggest that burnups of greater than 18 GWd/tU can be achieved in current reactor designs. At approximately 18 GWd/tU, expected maximum linear element ratings in a PT-HWR with online-refuelling are approximately 90 kW/m. These conditions would be prohibitive for solid-rod fuel, but may be possible in ICAF fuel given the reduced maximum fuel temperature as a function of linear element rating. (authors)