with high precision the neutron flux, burnup and concentration of every isotope, fissile, fissionable or fertile, gaseous or solid, all of them as functions of radius and time. But this formidable task is not suitable to be included in a fuel performance code, which must attend the great number of thermomechanical and thermochemical processes within the fuel rod. To accommodate both requirements, a simplified treatment is adopted consisting of restricting the balance equations to more relevant nuclides and reducing the energy spectrum to a single group. The purpose is to obtain empirical expressions to represent, with the higher possible approximation degree, the absorption, capture and fission cross sections of these isotopes as functions of the initial enrichment in 235U, the average burnup and the radial coordinate. The curves obtained with a so drastic simplification demand a careful testing before incorporation in the general fuel behaviour code. This testing is performed via comparison with the reliable reactor codes. The first antecedent in this type of analysis is found in the RADAR model [4] which was validated against the WIMS [5,6] code. The TUBRNP model, included in the TRANSURANUS code [7] and the RAPID model [8] are also based on the same concept. In this work curves fitted for the cross sections of 235U, 236U, 238U, 239Pu, 240Pu, 241Pu and 242Pu are obtained from the predictions of the reactor cell codes HUEMUL [9] and CONDOR [10] for an average burnup ranging from fresh fuel to 120 MWd/kgHM and for an initial enrichment ranging from natural uranium to 12%. The final purpose is to extend the application range of the DIONISIO code [11,12,13] (originally designed to predict the fuel behavior in normal operation conditions) to the high burnup domain. The predictions of DIONISIO were compared with a large number of experimental data, obtaining an excellent agreement
Burnup calculation code system COMRAD96
COMRAD was one of the burnup code system developed by JAERI. COMRAD96 is a transfered version of COMRAD to Engineering Work Station. It is divided to several functional modules, 'Cross Section Treatment', 'Generation and Depletion Calculation', and 'Post Process'. It enables us to analyze a burnup problem considering a change of neutron spectrum using UNITBURN. Also it can display the γ Spectrum on a terminal. This report is the general description and user's manual of COMRAD96. (author)
Integrated burnup calculation code system SWAT
SWAT is an integrated burnup code system developed for analysis of post irradiation examination, transmutation of radioactive waste, and burnup credit problem. It enables us to analyze the burnup problem using neutron spectrum depending on environment of irradiation, combining SRAC which is Japanese standard thermal reactor analysis code system and ORIGEN2 which is burnup code widely used all over the world. SWAT makes effective cross section library based on results by SRAC, and performs the burnup analysis with ORIGEN2 using that library. SRAC and ORIGEN2 can be called as external module. SWAT has original cross section library on based JENDL-3.2 and libraries of fission yield and decay data prepared from JNDC FP Library second version. Using these libraries, user can use latest data in the calculation of SWAT besides the effective cross section prepared by SRAC. Also, User can make original ORIGEN2 library using the output file of SWAT. This report presents concept and user's manual of SWAT. (author)
Systemization of burnup sensitivity analysis code. 2
Towards the practical use of fast reactors, it is a very important subject to improve prediction accuracy for neutronic properties in LMFBR cores from the viewpoint of improvements on plant efficiency with rationally high performance cores and that on reliability and safety margins. A distinct improvement on accuracy in nuclear core design has been accomplished by the development of adjusted nuclear library using the cross-section adjustment method, in which the results of criticality experiments of JUPITER and so on are reflected. In the design of large LMFBR cores, however, it is important to accurately estimate not only neutronic characteristics, for example, reaction rate distribution and control rod worth but also burnup characteristics, for example, burnup reactivity loss, breeding ratio and so on. For this purpose, it is desired to improve prediction accuracy of burnup characteristics using the data widely obtained in actual core such as the experimental fast reactor 'JOYO'. The analysis of burnup characteristics is needed to effectively use burnup characteristics data in the actual cores based on the cross-section adjustment method. So far, a burnup sensitivity analysis code, SAGEP-BURN, has been developed and confirmed its effectiveness. However, there is a problem that analysis sequence become inefficient because of a big burden to users due to complexity of the theory of burnup sensitivity and limitation of the system. It is also desired to rearrange the system for future revision since it is becoming difficult to implement new functions in the existing large system. It is not sufficient to unify each computational component for the following reasons; the computational sequence may be changed for each item being analyzed or for purpose such as interpretation of physical meaning. Therefore, it is needed to systemize the current code for burnup sensitivity analysis with component blocks of functionality that can be divided or constructed on occasion. For
Systemization of burnup sensitivity analysis code
To practical use of fact reactors, it is a very important subject to improve prediction accuracy for neutronic properties in LMFBR cores from the viewpoints of improvements on plant efficiency with rationally high performance cores and that on reliability and safety margins. A distinct improvement on accuracy in nuclear core design has been accomplished by development of adjusted nuclear library using the cross-section adjustment method, in which the results of critical experiments of JUPITER and so on are reflected. In the design of large LMFBR cores, however, it is important to accurately estimate not only neutronic characteristics, for example, reaction rate distribution and control rod worth but also burnup characteristics, for example, burnup reactivity loss, breeding ratio and so on. For this purpose, it is desired to improve prediction accuracy of burnup characteristics using the data widely obtained in actual core such as the experimental fast reactor core 'JOYO'. The analysis of burnup characteristics is needed to effectively use burnup characteristics data in the actual cores based on the cross-section adjustment method. So far, development of a analysis code for burnup sensitivity, SAGEP-BURN, has been done and confirmed its effectiveness. However, there is a problem that analysis sequence become inefficient because of a big burden to user due to complexity of the theory of burnup sensitivity and limitation of the system. It is also desired to rearrange the system for future revision since it is becoming difficult to implement new functionalities in the existing large system. It is not sufficient to unify each computational component for some reasons; computational sequence may be changed for each item being analyzed or for purpose such as interpretation of physical meaning. Therefore it is needed to systemize the current code for burnup sensitivity analysis with component blocks of functionality that can be divided or constructed on occasion. For this
High burnup models in computer code fair
An advanced fuel analysis code FAIR has been developed for analyzing the behavior of fuel rods of water cooled reactors under severe power transients and high burnups. The code is capable of analyzing fuel pins of both collapsible clad, as in PHWR and free standing clad as in LWR. The main emphasis in the development of this code is on evaluating the fuel performance at extended burnups and modelling of the fuel rods for advanced fuel cycles. For this purpose, a number of suitable models have been incorporated in FAIR. For modelling the fission gas release, three different models are implemented, namely Physically based mechanistic model, the standard ANS 5.4 model and the Halden model. Similarly the pellet thermal conductivity can be modelled by the MATPRO equation, the SIMFUEL relation or the Halden equation. The flux distribution across the pellet is modelled by using the model RADAR. For modelling pellet clad interaction (PCMI)/ stress corrosion cracking (SCC) induced failure of sheath, necessary routines are provided in FAIR. The validation of the code FAIR is based on the analysis of fuel rods of EPRI project ''Light water reactor fuel rod modelling code evaluation'' and also the analytical simulation of threshold power ramp criteria of fuel rods of pressurized heavy water reactors. In the present work, a study is carried out by analysing three CRP-FUMEX rods to show the effect of various combinations of fission gas release models and pellet conductivity models, on the fuel analysis parameters. The satisfactory performance of FAIR may be concluded through these case studies. (author). 12 refs, 5 figs
Revised SWAT. The integrated burnup calculation code system
SWAT is an integrated burnup code system developed for analysis of post irradiation examination, transmutation of radioactive waste, and burnup credit problem. This report shows an outline and a user's manual of revised SWAT. This revised SWAT includes expansion of functions, increasing supported machines, and correction of several bugs reported from users of previous SWAT. (author)
Revised SWAT. The integrated burnup calculation code system
Suyama, Kenya; Mochizuki, Hiroki [Department of Fuel Cycle Safety Research, Nuclear Safety Research Center, Tokai Research Establishment, Japan Atomic Energy Research Institute, Tokai, Ibaraki (Japan); Kiyosumi, Takehide [The Japan Research Institute, Ltd., Tokyo (Japan)
2000-07-01
SWAT is an integrated burnup code system developed for analysis of post irradiation examination, transmutation of radioactive waste, and burnup credit problem. This report shows an outline and a user's manual of revised SWAT. This revised SWAT includes expansion of functions, increasing supported machines, and correction of several bugs reported from users of previous SWAT. (author)
Burnup calculation methodology in the serpent 2 Monte Carlo code
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)
Validation of IRBURN calculation code system through burnup benchmark analysis
Assessment of the reactor fuel composition during the irradiation time, fuel management and criticality safety analysis require the utilization of a validated burnup calculation code system. In this work a newly developed burnup calculation code system, IRBURN, is introduced for the estimation and analysis of the fuel burnup in LWR reactors. IRBURN provides the full capabilities of the Monte Carlo neutron and photon transport code MCNP4C as well as the versatile code for calculating the buildup and decay of nuclides in nuclear materials, ORIGEN2.1, along with other data processing and linking subroutines. This code has the capability of using different depletion calculation schemes. The accuracy and precision of the implemented algorithms to estimate the eigenvalue and spent fuel isotope concentrations are demonstrated by validation against reliable benchmark problem analyses. A comparison of IRBURN results with experimental data demonstrates that the code predicts the spent fuel concentrations within 10% accuracy. Furthermore, standard deviations of the average values for isotopic concentrations including IRBURN data decreases considerably in comparison with the same parameter excluding IRBURN results, except for a few sets of isotopes. The eigenvalue comparison between our results and the benchmark problems shows a good prediction of the k-inf values during the entire burnup history with the maximum difference of 1% at 100 MWd/kgU.
Systemization of burnup sensitivity analysis code (2) (Contract research)
Towards the practical use of fast reactors, it is a very important subject to improve prediction accuracy for neutronic properties in LMFBR cores from the viewpoint of improvements on plant economic efficiency with rationally high performance cores and that on reliability and safety margins. A distinct improvement on accuracy in nuclear core design has been accomplished by the development of adjusted nuclear library using the cross-section adjustment method, in which the results of critical experiments of JUPITER and so on are reflected. In the design of large LMFBR cores, however, it is important to accurately estimate not only neutronic characteristics, for example, reaction rate distribution and control rod worth but also burnup characteristics, for example, burnup reactivity loss, breeding ratio and so on. For this purpose, it is desired to improve prediction accuracy of burnup characteristics using the data widely obtained in actual core such as the experimental fast reactor 'JOYO'. The analysis of burnup characteristic is needed to effectively use burnup characteristics data in the actual cores based on the cross-section adjustment method. So far, a burnup sensitivity analysis code, SAGEP-BURN, has been developed and confirmed its effectiveness. However, there is a problem that analysis sequence become inefficient because of a big burden to users due to complexity of the theory of burnup sensitivity and limitation of the system. It is also desired to rearrange the system for future revision since it is becoming difficult to implement new functions in the existing large system. It is not sufficient to unify each computational component for the following reasons: the computational sequence may be changed for each item being analyzed or for purpose such as interpretation of physical meaning. Therefore, it is needed to systemize the current code for burnup sensitivity analysis with component blocks of functionality that can be divided or constructed on occasion
Methods of RECORD, an LWR fuel assembly burnup code
The RECORD computer code is a detailed rector physics code for performing efficient LWR fuel assembly calculations, taking into account most of the features found in BWR and PWR fuel designs. The code calculates neutron spectrum, reaction rates and reactivity as a function of fuel burnup, and it generates the few-group data required for use in full scale core simulation and fuel management calculations. The report describes the methods of the RECORD computer code and the basis for fundamental models selected, and gives a review of code qualifications against measured data. (Auth. /RF)
Evolution of the ELESTRES code for application to extended burnups
The computer code ELESTRES is frequently used at Atomic Energy of Canada Limited to assess the integrity of CANDU fuel under normal operating conditions. The code also provides initial conditions for evaluating fuel behaviour during high-temperature transients. This paper describes recent improvements in the code in the areas of pellet expansion and of fission gas release. Both of these are very important considerations in ensuring fuel integrity at extended burnups. Firstly, in calculations of pellet expansion, the code now accounts for the effect of thermal stresses on the volume of gas bubbles at the boundaries of UO2 grains. This has a major influence on the expansion of the pellet during power-ramps. Secondly, comparisons with data showed that the previous fission gas package significantly underpredicted the fission gas release at high burnups. This package has now been improved via modifications to the following modules: distance between neighbouring bubbles on grain boundaries; diffusivity; and thermal conductivity. The predictions of the revised version of the code show reasonable agreement with measurements of ridge strains and of fission gas release. An illustrative example demonstrates that the code can be used to identify a fuel design that would: reduce the sheath stresses at circumferential ridges by a factor of 2-10; and keep the gas pressure at very high burnups to below the coolant pressure
OTTER 3 - A single channel, axial burnup code
OTTER 3 is a single channel, axial burnup code, written in Fortran for the KDF 9 computer, and suitable for studying fuel management schemes of the continuous charge/discharge type. A general fuel shuffling scheme is allowed, and both unidirectional and bidirectional fuel feed can be studied. A 2-group neutron diffusion code is incorporated, the flux equations being solved by the forward elimination - backward substitution technique for the inner problem and a source iteration technique accelerated by Chebyshev extrapolation for the outer problem. (author)
Burnup calculations using serpent code in accelerator driven thorium reactors
In this study, burnup calculations have been performed for a sodium cooled Accelerator Driven Thorium Reactor (ADTR) using the Serpent 1.1.16 Monte Carlo code. The ADTR has been designed for burning minor actinides, mixed 232Th and mixed 233U fuels. A solid Pb-Bi spallation target in the center of the core is used and sodium as coolant. The system is designed for a heating power of 2 000 MW and for an operation time of 600 days. For burnup calculations the Advanced Matrix Exponential Method CRAM (Chebyshev Rational Approximation Method) and different nuclear data libraries (ENDF7, JEF2.2, JEFF3.1.1) were used. The effective multiplication factor change from 0.93 to 0.97 for different nuclear data libraries during the reactor operation period. (orig.)
Burnup calculations using serpent code in accelerator driven thorium reactors
Korkmaz, M.E.; Agar, O. [Karamanoglu Mehmetbey Univ., Karaman (Turkey). Physics Dept.; Yigit, M. [Aksaray Univ. (Turkey). Physics Dept.
2013-07-15
In this study, burnup calculations have been performed for a sodium cooled Accelerator Driven Thorium Reactor (ADTR) using the Serpent 1.1.16 Monte Carlo code. The ADTR has been designed for burning minor actinides, mixed {sup 232}Th and mixed {sup 233}U fuels. A solid Pb-Bi spallation target in the center of the core is used and sodium as coolant. The system is designed for a heating power of 2 000 MW and for an operation time of 600 days. For burnup calculations the Advanced Matrix Exponential Method CRAM (Chebyshev Rational Approximation Method) and different nuclear data libraries (ENDF7, JEF2.2, JEFF3.1.1) were used. The effective multiplication factor change from 0.93 to 0.97 for different nuclear data libraries during the reactor operation period. (orig.)
BURNCAL: A Nuclear Reactor Burnup Code Using MCNP Tallies
BURNCAL is a Fortran computer code designed to aid in analysis, prediction, and optimization of fuel burnup performance in a nuclear reactor. The code uses output parameters generated by the Monte Carlo neutronics code MCNP to determine the isotopic inventory as a function of time and power density. The code allows for multiple fueled regions to be analyzed. The companion code, RELOAD, can be used to shuffle fueled regions or reload regions with fresh fuel. BURNCAL can be used to study the reactivity effects and isotopic inventory as a function of time for a nuclear reactor system. Neutron transmutation, fission, and radioactive decay are included in the modeling of the production and removal terms for each isotope of interest. For a fueled region, neutron transmutation, fuel depletion, fission-product poisoning, actinide generation, and burnable poison loading and depletion effects are included in the calculation. Fueled and un-fueled regions, such as cladding and moderator, can be analyzed simultaneously. The nuclides analyzed are limited only by the neutron cross section availability in the MCNP cross-section library. BURNCAL is unique in comparison to other burnup codes in that it does not use the calculated neutron flux as input to other computer codes to generate the nuclide mixture for the next time step. Instead, BURNCAL directly uses the neutron absorption tally/reaction information generated by MCNP for each nuclide of interest to determine the nuclide inventory for that region. This allows for the full capabilities of MCNP to be incorporated into the calculation and a more accurate and robust analysis to be performed
Miniature neutron source reactor burnup calculations using IRBURN code system
Highlights: ► Fuel consumption of Iranian MNSR during 15 years of operation has been investigated. ► Calculations have been performed by the IRBURN code. Precision and accuracy of the implemented model has been validated. ► Our study shows the consumption rate of MNSR is about 1%. - Abstract: Fuel consumption of Iranian miniature neutron source reactor (MNSR) during 15 years of operation has been investigated. Reactor core neutronic parameters such as flux and power distributions, control rod worth and effective multiplication factor at BOL and after 15 years of irradiation has been calculated. The Monte Carlo-based depletion code system IRBURN has been used for studying the reactor core neutronic parameters as well as the isotopic inventory of the fuel during burnup. The precision and accuracy of the implemented model has been verified via validation the results for neutronic parameters in the MNSR final safety analysis report. The results show that keff decreases from 1.0034 to 0.9897 and the total U-235 consumption in the core is about 13.669 g after 15 years of operational time. Finally, our studying shows the consumption rate of MNSR is about 1%.
Burnup calculations using the ORIGEN code in the CONKEMO computing system
This article describes the CONKEMO computing system for kinetic multigroup calculations of nuclear reactors and their physical characteristics during burnup. The ORIGEN burnup calculation code has been added to the system. The results of an international benchmark calculation are also presented. (author)
Development and validation of Monte-Carlo burnup calculation code MCNTRANS
A new nuclear fuel burnup calculation code MCNTRANS based on MCNP was introduced in this paper. The neutronics calculation parameter was extracted from the MCNP5 reaction rate tally result, while a graph theory algorithm was implemented to track the burnup chain and the analytic solution of the Bateman equation was given. At the same time, the detailed physical process was considered to improve the accuracy and serviceability of this code, and prediction-correction method was used to allow a large burnup step. The OECD/NEA and JAERI pin cell benchmark problems were used to validate the code MCNTRANS while a reference result was given by other code. It can be concluded that the calculation results of MCNTRANS are generally consistent with the experimental result and that of the other burnup codes, and part of the actinides and fission products calculation result show better accuracy. (authors)
Development of high-burnup fuel analysis code EXBURN-I
A computer code EXBURN-I has been developed which analyses LWR fuel behavior in high-burnup region in normal operation and transient conditions. In the high-burnup region, fuel behavior is affected considerably by such burnup-dependent factors as FP gas release, waterside corrosion of cladding, and pellet property change. To analyze these phenomena, in the present version, the base code FEMAXI-IV has been improved and incorporated such new models as pellet thermal conductivity change, burnup-dependent FP gas release rate, and cladding waterside corrosion. The present report describes the whole structure of the code, adopted models, and material properties, followed by input manual and sample input/output. Verification and further improvement of the code performance by experimental data will be done in the next stage. (author)
Development of an MCNP-tally based burnup code and validation through PWR benchmark exercises
The aim of this study is to evaluate the capabilities of a newly developed burnup code called BUCAL1. The code provides the full capabilities of the Monte Carlo code MCNP5, through the use of the MCNP tally information. BUCAL1 uses the fourth order Runge Kutta method with the predictor-corrector approach as the integration method to determine the fuel composition at a desired burnup step. Validation of BUCAL1 was done by code vs. code comparison. Results of two different kinds of codes are employed. The first one is CASMO-4, a deterministic multi-group two-dimensional transport code. The second kind is MCODE and MOCUP, a link MCNP-ORIGEN codes. These codes use different burnup algorithms to solve the depletion equations system. Eigenvalue and isotope concentrations were compared for two PWR uranium and thorium benchmark exercises at cold (300 K) and hot (900 K) conditions, respectively. The eigenvalue comparison between BUCAL1 and the aforementioned two kinds of codes shows a good prediction of the systems'k-inf values during the entire burnup history, and the maximum difference is within 2%. The differences between the BUCAL1 isotope concentrations and the predictions of CASMO-4, MCODE and MOCUP are generally better, and only for a few sets of isotopes these differences exceed 10%.
Development of 3d reactor burnup code based on Monte Carlo method and exponential Euler method
Burnup analysis plays a key role in fuel breeding, transmutation and post-processing in nuclear reactor. Burnup codes based on one-dimensional and two-dimensional transport method have difficulties in meeting the accuracy requirements. A three-dimensional burnup analysis code based on Monte Carlo method and Exponential Euler method has been developed. The coupling code combines advantage of Monte Carlo method in complex geometry neutron transport calculation and FISPACT in fast and precise inventory calculation, meanwhile resonance Self-shielding effect in inventory calculation can also be considered. The IAEA benchmark text problem has been adopted for code validation. Good agreements were shown in the comparison with other participants' results. (authors)
Revised Burnup Code System SWAT: Description and Validation Using Postirradiation Examination Data
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 UO2 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
Development and validation of a fast reactor core burnup code - FARCOB
Mohanakrishnan, P. [Indira Gandhi Centre for the Atomic Research, CDO, Reactor Physics Division, Kalpakkam, TN 603 102 (India)], E-mail: mohana@igcar.gov.in
2008-02-15
A large fast breeder reactor requires the accurate estimation of power produced in different parts of the reactor core and blanket during any operating condition for a safe and economic operation through out reactor life time. A fast reactor core simulation code FARCOB based on multigroup diffusion theory has been developed in IGCAR for core simulation of PFBR reactor under construction. FARCOB uses centre mesh differencing scheme with triangular meshes in the X-Y plane. Steady state solution results match exactly with those of other reputed codes DIF3D and VENTURE for SNR-300 benchmarks. For burnup simulation, core is divided into radial and axial burnup zones and burnup equations are solved at constant power. Burnable fuel and blanket number densities are found and stored for each mesh, so that the user can alter burnup zones and core geometry after a burnup step. For validation, results of FARCOB has been compared with results of other institutes in two burnup benchmarks (ANL 1000 MWe benchmark and BN-600 hybrid core benchmark). It is found that FARCOB results match well with those of the other institutes.
Siemens PWR burnup credit criticality analysis methodology: Depletion code and verification methods
Application of burnup credit requires knowledge of the reactivity state of the irradiated fuel for which burnup credit is taken. The isotopic inventory of the irradiated fuel has to be calculated, therefore, by means of depletion codes. Siemens performs depletion calculations for PWR fuel burnup credit applications with the aid of the code package SAV. This code package is based on the first principles approach, i.e., avoids cycle or reactor specific fitting or adjustment parameters. This approach requires a general and comprehensive qualification of SAV by comparing experimental with calculational results. In the paper on hand the attention is focused mainly on the evaluation of chemical assay data received from different experimental programmes. (author)
The purpose of this paper is to discuss the theories, techniques and computer codes that are frequently used in numerical reactor criticality and burnup calculations. It is a part of an integrated nuclear reactor calculation scheme conducted by the Reactors Department, Inshas Nuclear Research Centre. The crude part in numerical reactor criticality and burnup calculations includes the determination of neutron flux distribution which can be obtained in principle as a solution of Boltzmann transport equation. Numerical methods used for solving transport equations are discussed. Emphasis are made on numerical techniques based on multigroup diffusion theory. These numerical techniques include nodal, modal, and finite difference ones. The most commonly known computer codes utilizing these techniques are reviewed. Some of the main computer codes that have been already developed at the Reactors Department and related to numerical reactor criticality and burnup calculations have been presented
Validation of a new continuous Monte Carlo burnup code using a Mox fuel assembly
The reactivity of nuclear fuel decreases with irradiation (or burnup) due to the transformation of heavy nuclides and the formation of fission products. Burnup credit studies aim at accounting for fuel irradiation in criticality studies of the nuclear fuel cycle (transport, storage, etc...). The principal objective of this study is to evaluate the potential capabilities of a newly developed burnup code called 'BUCAL1'. BUCAL1 differs in comparison with other burnup codes as it does not use the calculated neutron flux as input to other computer codes to generate the nuclide inventory for the next time step. Instead, BUCAL1 directly uses the neutron reaction tally information generated by MCNP for each nuclide of interest to determine the new nuclides inventory. This allows the full capabilities of MCNP to be incorporated into the calculation and a more accurate and robust analysis to be performed. Validation of BUCAL1 was processed by code-to-code comparisons using predictions of several codes from the NEA/OCED. Infinite multiplication factors (k∞) and important fission product and actinide concentrations were compared for a MOX core benchmark exercise. Results of calculations are analysed and discussed.
New high burnup fuel models for NRC`s licensing audit code, FRAPCON
Lanning, D.D.; Beyer, C.E.; Painter, C.L. [Pacific Northwest Laboratory, Richland, WA (United States)
1996-03-01
Fuel behavior models have recently been updated within the U.S. Nuclear Regulatory Commission steady-state FRAPCON code used for auditing of fuel vendor/utility-codes and analyses. These modeling updates have concentrated on providing a best estimate prediction of steady-state fuel behavior up to the maximum burnup level s of current data (60 to 65 GWd/MTU rod-average). A decade has passed since these models were last updated. Currently, some U.S. utilities and fuel vendors are requesting approval for rod-average burnups greater than 60 GWd/MTU; however, until these recent updates the NRC did not have valid fuel performance models at these higher burnup levels. Pacific Northwest Laboratory (PNL) has reviewed 15 separate effects models within the FRAPCON fuel performance code (References 1 and 2) and identified nine models that needed updating for improved prediction of fuel behavior at high burnup levels. The six separate effects models not updated were the cladding thermal properties, cladding thermal expansion, cladding creepdown, fuel specific heat, fuel thermal expansion and open gap conductance. Comparison of these models to the currently available data indicates that these models still adequately predict the data within data uncertainties. The nine models identified as needing improvement for predicting high-burnup behavior are fission gas release (FGR), fuel thermal conductivity (accounting for both high burnup effects and burnable poison additions), fuel swelling, fuel relocation, radial power distribution, fuel-cladding contact gap conductance, cladding corrosion, cladding mechanical properties and cladding axial growth. Each of the updated models will be described in the following sections and the model predictions will be compared to currently available high burnup data.
Development of burnup calculation function in reactor Monte Carlo code RMC
This paper presents the burnup calculation capability of RMC, which is a new Monte Carlo (MC) neutron transport code developed by Reactor Engineering Analysis Laboratory (REAL) in Tsinghua University of China. Unlike most of existing MC depletion codes which explicitly couple the depletion module, RMC incorporates ORIGEN 2.1 in an implicit way. Different burn step strategies, including the middle-of-step approximation and the predictor-corrector method, are adopted by RMC to assure the accuracy under large burnup step size. RMC employs a spectrum-based method of tallying one-group cross section, which can considerably saves computational time with negligible accuracy loss. According to the validation results of benchmarks and examples, it is proved that the burnup function of RMC performs quite well in accuracy and efficiency. (authors)
Burnup calculation capability in the PSG2 / Serpent Monte Carlo reactor physics code
The PSG continuous-energy Monte Carlo reactor physics code has been developed at VTT Technical Research Centre of Finland since 2004. The code is mainly intended for group constant generation for coupled reactor simulator calculations and other tasks traditionally handled using deterministic lattices physics codes. The name was recently changed from acronym PSG to 'Serpent', and the capabilities have been extended by implementing built-in burnup calculation routines that enable the code to be used for fuel cycle studies and the modelling of irradiated fuels. This paper presents the methodology used for burnup calculation. Serpent has two fundamentally different options for solving the Bateman depletion equations: 1) the Transmutation Trajectory Analysis method (TTA), based on the analytical solution of linearized depletion chains and 2) the Chebyshev Rational Approximation Method (CRAM), an advanced matrix exponential solution developed at VTT. The first validation results are compared to deterministic CASMO-4E calculations. It is also shown that the overall running time in Monte Carlo burnup calculation can be significantly reduced using specialized calculation techniques, and that the continuous-energy Monte Carlo method is becoming a viable alternative to deterministic assembly burnup codes. (authors)
DELIGHT-6: one dimensional lattice burn-up code for high temperature gas-cooled reactors
The code, DELIGHT-6, performs multi-group neutron spectrum calculation and provides few-group constans for succeeding core calculations. The main objective of the code is to serve as the lattice burn-up code for the core of a very high temperature gas-cooled reactor. The fuel rods of the reactor contain many coated fuel particles resulting double heterogeneous arrangement. The main calculational schema of DELIGHT-6 code is as follows; (1) Energy range for fast neutrons covers from 10 MeV to 2.38 eV and is divided into 61 fine groups. The thermal neutrons covers the rest of the energy range from 2.38 eV to 0 eV. Thermal spectrum is calculated by P1 or P0 approximation with 50 fine groups. (2) To treat resonance absorption, IR method is employed. (3) Zero and one dimensional models are available for the fuel lattice geometry and used for criticality and burn-up calculations. Collision probability method is adopted for the calculation of one dimensional model. (4) Shielding factor of burnable poison is calculated by collision probability method. (5) Other functions of the code are; 1. Spatial shielding factor calculation of 240Pu, 2. Calculation of neutron streaming effect caused by a gap or a hole in the fuel lattice, 3. Calculation of neutron flux distribution in the fuel lattice by diffusion theory, 4. Calculation of Xe and Sm absorption cross sections with burn-up. (6) Cross section library in both fast and thermal energy range is compiled from ENDF/B-4 except burn-up data of Xm, Sm and pseudo FPs which are supplied by ENDF/B-3. (7) The code provides the macroscopic group constants of fuel lattice with burn-up in CITATION input format. (jin)
Accuracy assessment of a new Monte Carlo based burnup computer code
Highlights: ► A new burnup code called BUCAL1 was developed. ► BUCAL1 uses the MCNP tallies directly in the calculation of the isotopic inventories. ► Validation of BUCAL1 was done by code to code comparison using VVER-1000 LEU Benchmark Assembly. ► Differences from BM value were found to be ± 600 pcm for k∞ and ±6% for the isotopic compositions. ► The effect on reactivity due to the burnup of Gd isotopes is well reproduced by BUCAL1. - Abstract: This study aims to test for the suitability and accuracy of a new home-made Monte Carlo burnup code, called BUCAL1, by investigating and predicting the neutronic behavior of a “VVER-1000 LEU Assembly Computational Benchmark”, at lattice level. BUCAL1 uses MCNP tally information directly in the computation; this approach allows performing straightforward and accurate calculation without having to use the calculated group fluxes to perform transmutation analysis in a separate code. ENDF/B-VII evaluated nuclear data library was used in these calculations. Processing of the data library is performed using recent updates of NJOY99 system. Code to code comparisons with the reported Nuclear OECD/NEA results are presented and analyzed.
Point reactivity burnup code DELIGHT-4 for high temperature, gas-cooled reactor cells
The code DELIGHT-4 has been developed for analizing burnup characteristics of the graphite moderated reactor cells and producing the few-group constants. Calculation models for the code are as follows: (1) The number of neutron energy groups is 61 for fast neutrons (10 MeV -- 2.38 eV) and 50 for thermal neutrons (2.38 eV -- 0 eV). (2) The doubly space-heterogeneous effect of fuel (dispersion of coated fuel particles in fuel compacts and regular array of fuel rods in graphite blocks) is considered in the calculation of resonance absorption. (3) The double heterogenity of burnable poison (dispersion of absorber grains in rods) can be considered. (4) The chemical binding effect of graphite is introduced in the scattering of thermal neutrons. (5) The calculations of criticality and burnup are by a few-energy-group models (up to 10 groups for both fast and thermal neutrons), and nuclide chains of thorium-uranium and uranium-plutonium are used for burnup calculation. (6) Neutron streaming effect through holes and gaps in cells can be considered in criticality calculation. (7) The flux distribution in cells can be calculated. The cell-averaged few group constants can be produced in card form for 1-D transport approximation code SLALOM, 2-D S sub( n) code TWOTRAN, 1-D diffusion code BRIQUET, 2-D diffusion code ZADOC-3 and 3-D diffusion code CITATION-DEGA. (author)
Recent developments of the TRANSURANUS code with emphasis on high burnup phenomena
TRANSURANUS is a computer program for the thermal and mechanical analysis of fuel rods in nuclear reactors, which is developed at the Institute for Transuranium Elements. The code is in use in several European organisations, both in research and industry. In the paper the recent developments are summarised: the burnup degradation of the fuel thermal conductivity as well as the effects of gadolinium on the radial power distribution and thermal conductivity. Fission gas release from the High Burnup Structure is discussed. Finally, a new numerical method is outlined that is able to treat the highly non-linear mechanical equations in transients (RIAs and LOCAs). (author)
The burnup analysis program based on three dimensional discrete ordinates (SN) neutron/photon transport method has been developed by the FDS team, China, to aid in analysis, prediction, and optimization of fuel burnup performance in a nuclear reactor. The program uses output parameters generated by three-dimensional SN trans- port code to determine the isotopic inventory and anisotropic flux distribution as a function of time. For a fueled region, neutron transmutation, fuel depletion, fission-product poisoning, actinide generation, and burnable poison loading and depletion effects are included in the calculation. The IAEA benchmark test problem has been correctly calculated and analyzed to validate the system. (authors)
Burnup calculations of TR-2 Research Reactor with Monteburns Monte Carlo Code
Full text: In this study, some neutronic calculations of first and second core cycles of 5 MW pool type TR-2 Research Reactor have been performed using Multi-Step Monte Carlo Burnup Code System MONTEBURNS and the results were compared with the values of experiments and other codes. Time dependent keff distribution and burnup ratios belong to first and second core cycles of TR-2 Research Reactor were compared and quite good consistence in the results were observed. After modeling the first and second core cycles of TR-2 with MCNP5 Monte Carlo code, MCNP5 used in MONTEBURNS code has been parallelized in 8 HP ProLiant BL680C G5 systems with 4 quad-core Intel Xeon E7330 CPU, utilizing the MPI parallel protocol and simulations were performed on the 128 cores Linux parallel computing machine system. The computation time was reduced by parallelization of MONTEBURNS which uses MCNP in many steps. (authors)
SWAT, Step-Wise Burnup Analysis Code System to Combine SRAC-95 Cell Calculation Code and ORIGEN2
1 - Description of program or function: SWAT evaluates isotopic composition of spent nuclear fuel, especially for burnup credit issues by driving codes SRAC95 and ORIGEN2.1 or ORIGEN2. SWAT is an automated driver code system. At the initial development phase, it was constructed by combining source programs of SRAC and ORIGEN2. To overcome the problem associated with code updates, SWAT chose to use system function of UNIX operating system to execute SRAC95 and ORIGEN2. So that, SWAT is independent of development and modification of SRAC95 and ORIGEN2.1. In SWAT, ORIGEN2(82) or ORIGEN2.1 is used for burnup calculations using the matrix exponential method. An updated decay library is included in the distribution. SWAT uses SRAC95 for neutron spectrum and effective cross section calculation in 107 groups, using the collision probability method for given geometry and isotopic composition. One or two dimensional cell geometries are supported in SRAC95. NEA-1698/02: The main purpose of new package is to run SWAT on several machines not supported in previous package (IA64 under Linux, Windows with cygwin and Sun,...) and several commercial FORTRAN compiler (Intel, PGI, Fujitsu). 2 - Methods: In calculating the problem-dependent cross section in SWAT, the total burnup history is divided into 'burnup steps'. Power, boric acid concentration, temperature of each region, and void ratio of coolant are given as history data. For each burnup step, the neutron spectrum and effective cross section are evaluated by SRAC95 using the information given in previous burnup calculation and cell geometry information. The user can select geometry options for the collision probability method in SRAC95. 3 - Restrictions on the complexity of the problem: Resonance absorption calculation with ultra-fine group cross section can not be directly applicable for 2D geometry
Modification of the code SCTEMP and RIA transient analysis at high burnup
The code SCTEMP has been modified for RIA analysis of high burnup fuel. New models were introduced for this purpose, including rim type radial power profile, thermal conductivity degradation, and heat transfer coefficient covering different regimes. Thus, the thermal response during a RIA transient at high burnup can be simulated. An analytical exercise was made with parameter variations of pulse time amplitude, radial power profile, thermal conductivity degradation, energy deposition and heat transfer boundary conditions. In order to keep this exercise within a realistic frame, the fuel used as reference is a high burnup fuel which was tested at Halden and for which the fuel thermal characteristics are known based on in-pile fuel temperature measurements. The calculation results are analysed and discussed in terms of their significance for fuel behaviour under the transient and for design of possible RIA experiments. Pulse time amplitude effects and energy deposition effects of high burnup fuel were evaluated quantitatively. The rim of high burnup fuel is assumed to play a very important role for RIA transients. (author)
Using SERPENT Monte Carlo and Burnup code to model Traveling Wave Reactors - TWR
This paper is mainly devoted to the proof-of-principle implementation of the SERPENT code for the simulation of traveling wave reactors. Traveling wave reactors are both fast reactors and nuclear burning wave reactors in which the breeding and burning of nuclear fuel appear almost simultaneously. SERPENT is a neutron transport code whose last official update package is SERPENT 1.1.19 and whose SERPENT 2 version is currently in progress. The investigation of SERPENT 1.1.19 and of SERPENT 2 codes for multiprocessor tasks with long burnup steps was performed. It appears that SERPENT 2 has eliminated parallelization problems efficiently. Methods to remove the influence of the ignition zone were considered, and neutron transport simulations with various fragmentations of the burnup zone were performed. (authors)
There have been two versions of SWAT depending on details of its development history: the revised SWAT that uses the deterministic calculation code SRAC as a neutron transportation solver, and the SWAT3.1 that uses the continuous energy Monte Carlo code MVP or MCNP5 for the same purpose. It takes several hours, however, to execute one calculation by the continuous energy Monte Carlo code even on the super computer of the Japan Atomic Energy Agency. Moreover, two-dimensional burnup calculation is not practical using the revised SWAT because it has problems on production of effective cross section data and applying them to arbitrary fuel geometry when a calculation model has multiple burnup zones. Therefore, SWAT4.0 has been developed by adding, to SWAT3.1, a function to utilize the deterministic code SARC2006, which has shorter calculation time, as an outer module of neutron transportation solver for burnup calculation. SWAT4.0 has been enabled to execute two-dimensional burnup calculation by providing an input data template of SRAC2006 to SWAT4.0 input data, and updating atomic number densities of burnup zones in each burnup step. This report describes outline, input data instruction, and examples of calculations of SWAT4.0. (author)
In the design of the incore thermionic reactor system developed under the Advanced Thermionic Initiative (ATI), the fuel is highly enriched uranium dioxide and the moderating medium is zirconium hydride. The traditional burnup and fuel depletion analysis codes have been found to be inadequate for these calculations, largely because of the material and geometry modeled and because the neutron spectra assumed for the codes such as LEOPARD and ORIGEN do not even closely fit that for a small, thermal reactor using ZrH as moderator. More sophisticated codes such as the transport lattice type code WIMS often lack some materials, such as ZrH. Thus a new method which could accurately calculate the neutron spectrum and the appropriate reaction rates within the fuel element is needed. The method developed utilizes and interconnects the accuracy of the Monte Carlo Neutron/Photon (MCNP) method to calculate reaction rates for the important isotopes, and a time dependent depletion routine to calculate the temporal effects on isotope concentrations. This effort required the modification of MCNP itself to perform the additional task of accomplishing burnup calculations. The modified version called, MCNPBURN, evolved to be a general dual purpose code which can be used for standard calculations as well as for burn-up
Horhoianu, G.; Moscalu, D.R.; Popescu, I.A
1998-06-01
Extending burnup is a practical way to improve the economics of water-reactor operation, via enhanced fuel utilisation and reduced spent fuel volume. A dedicated fuel behaviour modeling computer code (entitled ROFEM-1B) has been developed in order to analyse high burnup fuel performance. The code was benchmarked on an experimental data base which include a significant number of irradiation experiments performed in TRIGA-INR Pitesti research reactor. Five fuel rod behaviour during irradiation up to 50 MWd kg{sup -1}UO{sub 2}{sup -1} burnup have been analysed by the code in the framework of the first phase of the international FUMEX code exercise co-ordinated by IAEA Vienna. The input experimental data package has been prepared by IFE-OECD Halden, Norway laboratory. In the second phase of the FUMEX exercise the participants have analysed eight simplified theoretical cases. The paper presents and discuss the results obtained with ROFEM-1B in this exercise and the comparison between code predictions and experimental data.
Development of a fuel rod thermal-mechanical analysis code for high burnup fuel
The thermal-mechanical analysis code for high burnup BWR fuel rod has been developed by NFI. The irradiation data accumulated up to the assembly burnup of 55 GWd/t in commercial BWRs were adopted for the modeling. In the code, pellet thermal conductivity degradation with burnup progress was considered. Effects of the soluble FPs, irradiation defects and porosity increase due to RIM effect were taken into the model. In addition to the pellet thermal conductivity degradation, the pellet swelling due to the RIM porosity was studied. The modeling for the high burnup effects was also carried out for (U, Gd)O2 and MOX fuel. The thermal conductivities of all pellet types, UO2, (U, Gd)O2 and (U, Pu)O2 pellets, are expressed by the same form of equation with individual coefficient γ in the code. The pellet center temperature was calculated using this modeling code, and compared with measured values for the code verification. The pellet center temperature calculated using the thermal conductivity degradation model agreed well with the measured values within ±150 deg. C. The influence of rim porosity on pellet center temperature is small, and the temperature increase in only 30 deg. C at 75 GWd/t and 200 W/cm. The pellet center temperature of MOX fuel was also calculated, and it was found that the pellet center temperature of MOX fuel with 10wt% PuO2 is about 60 deg. C higher than UO2 fuel at 75 GWd/t and 200 W/cm. (author)
Progress on burnup calculation methods coupling Monte Carlo and depletion codes
Leszczynski, Francisco [Comision Nacional de Energia Atomica, San Carlos de Bariloche, RN (Argentina). Centro Atomico Bariloche]. E-mail: lesinki@cab.cnea.gob.ar
2005-07-01
Several methods of burnup calculations coupling Monte Carlo and depletion codes that were investigated and applied for the author last years are described. here. Some benchmark results and future possibilities are analyzed also. The methods are: depletion calculations at cell level with WIMS or other cell codes, and use of the resulting concentrations of fission products, poisons and actinides on Monte Carlo calculation for fixed burnup distributions obtained from diffusion codes; same as the first but using a method o coupling Monte Carlo (MCNP) and a depletion code (ORIGEN) at a cell level for obtaining the concentrations of nuclides, to be used on full reactor calculation with Monte Carlo code; and full calculation of the system with Monte Carlo and depletion codes, on several steps. All these methods were used for different problems for research reactors and some comparisons with experimental results of regular lattices were performed. On this work, a resume of all these works is presented and discussion of advantages and problems found are included. Also, a brief description of the methods adopted and MCQ system for coupling MCNP and ORIGEN codes is included. (author)
First steps towards a validation of the new burnup and depletion code TNT
In the frame of the fusion of the core design calculation capabilities, represented by V.S.O.P., and the accident calculation capabilities, represented by MGT(-3D), the successor of the TINTE code, difficulties were observed in defining an interface between a program backbone and the ORIGEN code respectively the ORIGENJUEL code. The estimation of the effort of refactoring the ORIGEN code or to write a new burnup code from scratch, led to the decision that it would be more efficient writing a new code, which could benefit from existing programming and software engineering tools from the computer code side and which can use the latest knowledge of nuclear reactions, e.g. consider all documented reaction channels. Therefore a new code with an object-oriented approach was developed at IEK-6. Object-oriented programming is currently state of the art and provides mostly an improved extensibility and maintainability. The new code was named TNT which stands for Topological Nuclide Transformation, since the code makes use of the real topology of the nuclear reactions. Here we want to present some first validation results from code to code benchmarks with the codes ORIGEN V2.2 and FISPACT2005 and whenever possible analytical results also used for the comparison. The 2 reference codes were chosen due to their high reputation in the field of fission reactor analysis (ORIGEN) and fusion facilities (FISPACT). (orig.)
An iterative method for burn-up determination in the non-destructive analysis of irradiated fuel elements using the ORIGEN computer code is presented. On the bases of data obtained from ORIGEN code the calibration coefficient for the neutron flux is determined as a function of one fission product activity while the burn-up is determined as a function of the calibration coefficient for a given irradiation history. These functions are used for determining the burn-up of nuclear fuel elements measured by gamma-scanning. The method is tested for fuel elements irradiated in a TRIGA reactor facility. (Author)
Development of a MCNP–ORIGEN burn-up calculation code system and its accuracy assessment
Highlights: • MCNP and ORIGEN are coupled to perform nuclides depletion and decay calculation. • Coupled system MCORE uses “modified predictor corrector” approach. • MCORE can use different depletion schemes and simulate fuel shuffling. • MCORE is assessed by a “VVER-1000 LEU Assembly Computational Benchmark”. • MCORE is also assessed by a fast reactor benchmark problem. - Abstract: An MCNP–ORIGEN burn-up calculation code system, named MCORE (MCNP and ORIGEN burn-up Evaluation code), is developed in this work. MCORE makes use of the Monte Carlo neutron and photon transport code MCNP4C and nuclides depletion and decay calculation code ORIGEN2.1. MCNP and ORIGEN are coupled by data processing and linking subroutines. In MCORE, a so called “modified predictor corrector” approach is used. MCORE provides the capability of using different depletion calculation schemes and simulating fuel shuffling. Total nuclide density changes in active cells are considered in MCORE. The validity and applicability of the developed code are tested by investigating and predicting the neutronic and isotopic behavior of a “VVER-1000 LEU Assembly Computational Benchmark” at lattice level and a “Physics of Plutonium Recycling” fast reactor at core level (OECD-NEA). The comparison results show that the MCORE code predicts the nuclide composition within 5% accuracy and k∞ within 800 pcm at the end of the burn-up for LEU assembly (40 MWD/kg HM). For a fast reactor, the results obtained by MCORE are in the range of reported results except for 243Am. In general, MCORE results show a good agreement with the benchmark values
DELIGHT-6(revised): one dimensional lattice burnup code for high temperature gas cooled reactors
The code, DELIGHT-6, performs the multi-group neutron spectrum calculation and provides the few-group constants for burnup calculations of a high temperature gas-cooled reactor core, whose fuel elements containing many coated fuel particles are arranged in double heterogeneity. The main revisions in the DELIGHT-6 (Revised) are as follows; (1)The option of a sphere fuel cell calculation is added for the core design of pebble bed type high temperature gas-cooled reactor. (2)The yield and decay constants of fission products for burnup calculation is revised. (3)The following auxiliary functions are added; (i) Automatic calculation of averaged atom number density in the fuel region, (ii) Estimation of local neutron flux distribution (disadvantage factor), (iii) Preparation of the data for the fine mesh core calculation. (author)
MODRIB - a zero dimensional code for criticality and burn-up of LWR's
The computer program MODRIB is a zero-dimensional code for calculating criticality and burn-up of light water reactors (LWR's). It is a version of an Italian code RIBOT-2 with an updated cross-section data library. The nuclear constants of MODRIB-code are calculated with a two group scheme (fast and thermal), where the fast group is an average of three fast groups. The code requires as input data essential extensive reactor parameters such as fuel rod radius, clad thickness, fuel enrichment, lattice pitch, water density and temperature etc. A summary of the physical model description and the input-output procedures are given in this report. Selected results of two sample problems are also given for the purpose of checking the validity and reliability of the code. The first is BWR and the second is PWR. The calculation time for a criticality problem with burn-up is about 8 seconds for the first time step and about 3 seconds for each subsequent time step on the ICL-1906 computer facility. The requirements on the memory size is less than 32 K-word. (author)
ALEPH 1.1.2: A Monte Carlo burn-up code
In the last 40 years, Monte Carlo particle transport has been applied to a multitude of problems such as shielding and medical applications, to various types of nuclear reactors, . . . The success of the Monte Carlo method is mainly based on its broad application area, on its ability to handle nuclear data not only in its most basic but also most complex form (namely continuous energy cross sections, complex interaction laws, detailed energy-angle correlations, multi-particle physics, . . . ), on its capability of modeling geometries from simple 1D to complex 3D, . . . There is also a current trend in Monte Carlo applications toward high detail 3D calculations (for instance voxel-based medical applications), something for which deterministic codes are neither suited nor performant as to computational time and precision. Apart from all these fields where Monte Carlo particle transport has been applied successfully, there is at least one area where Monte Carlo has had limited success, namely burn-up and activation calculations where the time parameter is added to the problem. The concept of Monte Carlo burn-up consists of coupling a Monte Carlo code to a burn-up module to improve the accuracy of depletion and activation calculations. For every time step the Monte Carlo code will provide reaction rates to the burn-up module which will return new material compositions to the Monte Carlo code. So if static Monte Carlo particle transport is slow, then Monte Carlo particle transport with burn-up will be even slower as calculations have to be performed for every time step in the problem. The computational issues to perform accurate Monte Carlo calculations are however continuously reduced due to improvements made in the basic Monte Carlo algorithms, due to the development of variance reduction techniques and due to developments in computer architecture (more powerful processors, the so-called brute force approach through parallel processors and networked systems
Zagar, T.; Ravnik, M.; Persic, A. (J.Stefan Institute, Ljubljana (Slovenia))
1999-12-15
Results of fuel element burn-up determination by measurement and calculation are given. Fuel element burn-up was calculated with two different programs TRIGLAV and TRIGAC using different models. New TRIGLAV code is based on cylindrical, two-dimensional geometry with four group diffusion approximation. TRIGAC program uses one-dimensional cylindrical geometry with twogroup diffusion approximation. Fuel element burn-up was measured with reactivity method. In this paper comparison and analysis of these three methods is presented. Results calculated with TRIGLAV show considerably better alignment with measured values than results calculated with TRIGAC. Some two-dimensional effects in fuel element burn-up can be observed, for instance smaller standard fuel element burn-up in mixed core rings and control rod influence on nearby fuel elements. (orig.)
ELESTRES 2.1 computer code for high burnup CANDU fuel performance analysis
The ELESTRES (ELEment Simulation and sTRESses) computer code models the thermal, mechanical and micro structural behaviours of CANDU® fuel element under normal operating conditions. The main purpose of the code is to calculate fuel temperatures, fission gas release, internal gas pressure, fuel pellet deformation, and fuel sheath strains in fuel element design analysis and assessments. It is also used to provide initial conditions for evaluating fuel behaviour during high temperature transients. ELESTRES 2.1 was developed for high burnup fuel application, based on an industry standard tool version of the code, through the implementation or modification to code models such as fission gas release, fuel pellet densification, flux depression (radial power distribution in the fuel pellet), fuel pellet thermal conductivity, fuel sheath creep, fuel sheath yield strength, fuel sheath oxidation, two dimensional heat transfer between the fuel pellet and the fuel sheath; and an automatic finite element meshing capability to handle various fuel pellet shapes. The ELESTRES 2.1 code design and development was planned, implemented, verified, validated, and documented in accordance with the AECL software quality assurance program, which meets the requirements of the Canadian Standards Association standard for software quality assurance CSA N286.7-99. This paper presents an overview of the ELESTRES 2.1 code with descriptions of the code's theoretical background, solution methodologies, application range, input data, and interface with other analytical tools. Code verification and validation results, which are also discussed in the paper, have confirmed that ELESTRES 2.1 is capable of modelling important fuel phenomena and the code can be used in the design assessment and the verification of high burnup fuels. (author)
In order to account for the reactivity-reducing effect of burn-up in the criticality safety analysis for systems with irradiated nuclear fuel (''burnup credit''), numerical methods to determine the enrichment and burnup dependent nuclide inventory (''burnup code'') and its resulting multiplication factor keff (''criticality code'') are applied. To allow for reliable conclusions, for both calculation systems the systematic deviations of the calculation results from the respective true values, the bias and its uncertainty, are being quantified by calculation and analysis of a sufficient number of suitable experiments. This quantification is specific for the application case under scope and is also called validation. GRS has developed a methodology to validate a calculation system for the application of burnup credit in the criticality safety analysis for irradiated fuel assemblies from pressurized water reactors. This methodology was demonstrated by applying the GRS home-built KENOREST burnup code and the criticality calculation sequence CSAS5 from SCALE code package. It comprises a bounding approach and alternatively a stochastic, which both have been exemplarily demonstrated by use of a generic spent fuel pool rack and a generic dry storage cask, respectively. Based on publicly available post irradiation examination and criticality experiments, currently the isotopes of uranium and plutonium elements can be regarded for.
A multi-platform linking code for fuel burnup and radiotoxicity analysis
A linking code between ORIGEN2.1 and MCNP has been developed at the Departamento de Engenharia Nuclear/UFMG to calculate coupled neutronic/isotopic results for nuclear systems and to produce a large number of criticality, burnup and radiotoxicity results. In its previous version, it evaluated the isotopic composition evolution in a Heat Pipe Power System model as well as the radiotoxicity and radioactivity during lifetime cycles. In the new version, the code presents features such as multi-platform execution and automatic results analysis. Improvements made in the code allow it to perform simulations in a simpler and faster way without compromising accuracy. Initially, the code generates a new input for MCNP based on the decisions of the user. After that, MCNP is run and data, such as recoverable energy per prompt fission neutron, reaction rates and keff, are automatically extracted from the output and used to calculate neutron flux and cross sections. These data are then used to construct new ORIGEN inputs, one for each cell in the core. Each new input is run on ORIGEN and generates outputs that represent the complete isotopic composition of the core on that time step. The results show good agreement between GB (Coupled Neutronic/Isotopic code) and Monteburns (Automated, Multi-Step Monte Carlo Burnup Code System), developed by the Los Alamos National Laboratory. - Highlights: • Multi-platform execution and automatic results analysis are implemented. • Improvements allow simulations in a simpler and faster way. • The results show good agreement between GB and Monteburns codes
LOLA-SYSTEM, JEN-UPM PWR Fuel Management System Burnup Code System
1 - Description of program or function: The LOLA-SYSTEM is a part of the JEN-UPM code package for PWR fuel management, scope or design calculations. It is a code package for core burnup calculations using nodal theory based on a FLARE type code. The LOLA-SYSTEM includes four modules: the first one (MELON-3) generates the constants of the K-inf and M2 correlations to be input into SIMULA-3. It needs the K-inf and M2 fuel assembly values at different conditions of moderator temperature, Boron concentration, burnup, etc., which are provided by MARIA fuel assembly calculations. The main module (SIMULA-3) is the core burnup calculation code in three dimensions and one group of energy. It normally uses a geometrical representation of one node per fuel assembly or per quarter of fuel assembly. It has included a thermal hydraulic feedback on flow and voids and criticality searches on boron concentration and control rods insertion. The CONCON code makes the calculation of the albedo, transport factors, K-inf and M2 correction factors to be input into SIMULA-3. The calculation is made in the XY transversal plane. The CONAXI code is similar to CONCON, but in the axial direction. 2 - Method of solution: MELON-3 makes a mean squares fit of K-inf and M2 values at different conditions in order to determine the constants of the feedback correlations. SIMULA-3 uses a modified one-group nodal theory, with a new transport kernel that provides the same node interface leakages as a fine mesh diffusion calculation. CONCON and CONAXI determine the transport and correction factors, as well as the albedo, to be input into SIMULA-3. They are determined by a method of leakages equivalent to the detailed diffusion calculation of CARMEN or VENTURE; these factors also include the heterogeneity effects inside the node. 3 - Restrictions on the complexity of the problem: Number of axial nodes less than or equal 34. Number of material types less than or equal 30. Number of fuel assembly types less
High burnup fuel simulation during a reactivity initiated transient using the FRAPTRAN code
CIP01 test simulated a RIA (Reactivity Initiated Accident) transient at the experimental CABRI reactor under sodium cooling conditions on November 29th 2001. The father rod was fabricated by ENUSA and irradiated during five cycles in Vandellos 2 NPP and it reached a maximum burnup at EOL (End of Life) of 75 GWd/tU. The test rodlet was refabricated in Studsvik from the span 5 of the father rod. The cladding material was ZIRLO with a mean oxide layer of 75 μm. This note describes the modelling of the CIP01 test carried out in the CABRI reactor under RIA conditions by the use of the transient FRAPTRAN code (1.2 version). The CIP01 test consisted in a power pulse of 32.4 ms with a deposited energy of 98.2 cal/g at PPN (Peak Power Node). The code results will be presented and compared to the experimental data available. In addition, the FRAPTRAN and the SCANAIR code estimations will be also contrasted. This work is framed within the CSN-CIEMAT agreement on 'Thermo-Mechanical Behaviour of the Nuclear Fuel at High Burnup'. (Author)
Burnup span sensitivity analysis of different burnup coupling schemes
Highlights: ► The objective of this work is the burnup span sensitivity analysis of different coupling schemes. ► Three kinds of schemes have been implemented in a new MCNP–ORIGEN linkage program. ► Two kinds of schemes are based predictor–corrector technique and the third is based on Euler explicit method. ► The analysis showed that the predictor–corrector approach better accounts for nonlinear behavior of burnup. ► It is sufficiently good to use the Euler method at small spans but for large spans use of second order scheme is mandatory. - Abstract: The analysis of core composition changes is complicated by the fact that the time and spatial variations in isotopic composition depend on the neutron flux distribution and vice versa. Fortunately, changes in core composition occur relatively slowly and hence the burnup analysis can be performed by dividing the burnup period into some burnup spans and assuming that the averaged flux and cross sections are constant during each burn up span. The burnup span sensitivity analysis attempts to find how much the burnup spans could be increased without any significant change in results. This goal has been achieved by developing a new MCNP–ORIGEN linkage program named MOBC (MCNP–ORIGEN Burnup Calculation). Three kinds of coupling scheme have been implemented in MOBC. Two of these are based on second order predictor–corrector technique and enable us to choose larger time steps, whilst the third one is based on Euler explicit first order method and is faster than the other two. The validity of the developed program has been evaluated by the code vs. code comparison technique. Two different types of codes are employed. The first one is based on deterministic two dimensional transport method, like CASMO-4 and HELIOS codes, and the second one is based on Monte Carlo method, like MCODE code. Only one coupling technique is employed in each of these state of the art codes, while the MOBC excels in its ability to
Development of a code and models for high burnup fuel performance analysis
First the high burnup LWR fuel behavior is discussed and necessary models for the analysis are reviewed. These aspects of behavior are the changes of power history due to the higher enrichment, the temperature feedback due to fission gas release and resultant degradation of gap conductance, axial fission gas transport in fuel free volume, fuel conductivity degradation due to fission product solution and modification of fuel micro-structure. Models developed for these phenomena, modifications in the code, and the benchmark results mainly based on Risoe fission gas project is presented. Finally the rim effect which is observe only around the fuel periphery will be discussed focusing into the fuel conductivity degradation and swelling due to the porosity development. (author). 18 refs, 13 figs, 3 tabs
The predictions of the ENIGMA code have been compared with data from high burn-up fuel experiments from the Halden and RISO reactors. The experiments modelled were IFA-504 and IFA-558 from Halden and the test II-5 from the RISO power burnup test series. The code has well modelled the fuel thermal performance and has provided a good measure of iodine release from pre-interlinked fuel. After interlinkage the iodine predictions remain a good fit for one experiment, but there is significant overprediction for a second experiment (IFA-558). Stable fission gas release is also well modelled and the predictions are within the expected uncertainly band throughout the burn-up range. This report presents code predictions for stable fission gas release to 40GWd/tU, iodine release measurements to 50GWd/tU and thermal performance (fuel centre temperature) to 55GWd/tU. Fuel ratings of up to 38kW/m were modelled at the high burn-up levels. The code is shown to accurately or conservatively predict all these parameters. (author). 1 ref., 6 figs
Mathematic preprocessor for RELAP5 code using Microsoft Excel
Computational program are used for thermal hydraulic analysis of accidents and transients conditions in nuclear power plants. The RELAP5 code has been developed to simulate accidents and transients conditions, performing a best estimate analysis, in Pressurized Water Reactors (PWR) and auxiliary systems. The RELAP5 code, which has been used as a toll for licensing nuclear facilities in Brazil, is the objective of the study performed in this work. The main problem in using the RELAP5 code is the huge amount of information necessary to model the nuclear reactor and thus to simulate thermal-hydraulic accidents. Moreover, the RELAP5 code input data requires a large amount of mathematical operations to calculate the geometry of the plant components. Therefore, in order to make easier the data input for the RELAP5 code a friendly preprocessor has been developed. The preprocessor accepts basic information about the geometry of the plant components and performs all the calculations needed for the RELAP5 input. This preprocessor has been developed based on the MS-Excel software. (author)
The high burnup database for CANDU fuel includes several cases from both power station and experimental reactor irradiations, with achieved burnups of up to 800 MW.h/kgU. The power history for each of these cases is different, encompassing low steady-state, declining, and power-ramps. This variety offers a good opportunity to check the models of fuel behaviour, and to identify areas for improvement. The main parameters for comparing calculated versus measured data are the fission gas release and the sheath hoop strain. Good agreement of calculated values of these two parameters with experimental data indicates that the global behaviour of the fuel element is adequately simulated by our codes. The ELESIM computer code was used as the simulation tool. The models for fission gas release, swelling and for fuel pellet expansion were thoroughly analysed. Changes were proposed for both models. The fuel pellet expansion model was modified to account for gaseous swelling, which becomes very important at high burnups. As well, the mathematics of the fission gas release model was upgraded for the diffusional release of fission gas atoms to the grain boundaries. A revised version of the ELESIM computer code was used to simulate the cases from the high burnup database. Satisfactory agreement was found for most cases. The discrepancies are discussed in view of alternative mechanisms that can operate and be enhanced at high burnup. These include stoichiometry changes with burnup that affects fission gas release, and also outer pellet rim fission gas release by a grain boundary diffusion process. The main conclusion of this study is that the revised version of the ELESIM code is able to simulate with reasonable accuracy high burnup as well as low burnup CANDU fuel. This includes irradiations of steady-state, declining, or ramped fuel power histories with a prolonged hold at high power. However, future improvements to ELESIM are needed to model fuel power histories with short dwell
A PWR PCI failure criterion to burnups of 60 GW·d/t using the ENIGMA code
A fuel performance modelling code (ENIGMA) has been used to analyse the empirical PCI failure criterion in terms of a clad failure stress as a function of burnup and fast neutron dose. The Studsvik database has been analysed. Results indicate a rising and then saturating failure stress with burnup and fast neutron dose. Using the PCI failure limits, equivalent to 95/95 confidence limits, an ENIGMA stress-based methodology is used to derive PWR PCI failure limits up to 60 GW·d/t U using a conservative assumption that the failure stress does not increase at high burnup and neutron dose. In addition experimental ramp data on gadolinia-doped fuel rods do not indicate any increased susceptibility to PCI failure implying that the UO2 criterion can be used for gadolinia doped fuel. (author)
By coupling a three-dimensional Discrete Ordinate (SN) code with an Exponential Euler Method inventory code, a three-dimensional code system was developed for accurate and rapid fuel burnup and material activation analyses. Resonance self-shielding effect was corrected dynamically in each step and each fuel zone. The update of compositions and nuclide densities in transport calculation was treated by setting nuclide density limitation, keeping important nuclides and using pseudo fission products (FPeff). The validation was described, which illustrates the robustness and reliability of this code system. (authors)
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.
The program system KENOREST version 1998 will be presented, which is a useful tool for burnup and reactivity calculations for LWR fuel. The three-dimensional Monte Carlo code KENO-V.a is coupled with the one-dimensional GRS burnup program system OREST-98. The objective is to achieve a better modelling of plutonium and actinide build-up or burnout for advanced heterogeneous fuel assembly designs. Further objectives are directed to reliable calculations of the pin power distributions and of reactor safety parameters including axial and radial rod temperatures for fuel assemblies of modern design. The stand-alone-code KENO-V.a version is used without any changes in the program source. The OREST-98 system was developed to handle multirod problems and additional burnup dependent moderator conditions which can be applied to stretch-out simulations in the reactor. A new interface module RESPEFF between KENO and OREST transforms the 2-d or 3-d KENO flux results to the one-dimensional lattice code OREST in a fully automated manner to maintain reaction rate balance between the codes. First results for assembly multiplication factors, isotope inventories are compared with OECD results. (author)
Highlights: • Development of a neural network model to predict the requiered plutonium content. • The accuracy of this model is very good (0.37% of error). • Development of a neural network model to predict evolution of average cross sections. • Predictions allow calculating fuel depletion quickly and with a very good accuracy. • This approach has been applied to the PWR MOX case in a dynamic fuel cycle code. - Abstract: Dynamic fuel cycle simulation tools calculate nuclei inventories and mass flows evolution in an entire fuel cycle, from the mine to the final disposal. Usually, the fuel depletion in reactor is handled by a fuel loading model and a mean cross section predictor. In the case of a PWR–MOX, a fuel loading model provides from a plutonium stock the plutonium fraction in the fresh fuel needed to reach a specific burnup. A mean cross section predictor aims to assess isotopic cross sections required for building Bateman equations for any fresh fuel composition with a sufficient accuracy and a reasonable computing time. This paper presents a methodology based on neural networks for building a fuel loading model and a cross section predictor for a PWR reactor loaded with MOX fuel. The mean error of the plutonium content prediction from the fuel loading model is 0.37%. Furthermore, the mean cross section predictor allows completion of the fuel depletion calculation in less than one minute with excellent accuracy. A maximum deviation of 3% on main nuclei is obtained at the end of cycle between inventories calculated from neural networks and from the reference coupled neutron transport/fuel depletion calculation
Behler, Matthias; Hannstein, Volker; Kilger, Robert; Moser, Franz-Eberhard; Pfeiffer, Arndt; Stuke, Maik
2014-06-15
In order to account for the reactivity-reducing effect of burn-up in the criticality safety analysis for systems with irradiated nuclear fuel (''burnup credit''), numerical methods to determine the enrichment and burnup dependent nuclide inventory (''burnup code'') and its resulting multiplication factor k{sub eff} (''criticality code'') are applied. To allow for reliable conclusions, for both calculation systems the systematic deviations of the calculation results from the respective true values, the bias and its uncertainty, are being quantified by calculation and analysis of a sufficient number of suitable experiments. This quantification is specific for the application case under scope and is also called validation. GRS has developed a methodology to validate a calculation system for the application of burnup credit in the criticality safety analysis for irradiated fuel assemblies from pressurized water reactors. This methodology was demonstrated by applying the GRS home-built KENOREST burnup code and the criticality calculation sequence CSAS5 from SCALE code package. It comprises a bounding approach and alternatively a stochastic, which both have been exemplarily demonstrated by use of a generic spent fuel pool rack and a generic dry storage cask, respectively. Based on publicly available post irradiation examination and criticality experiments, currently the isotopes of uranium and plutonium elements can be regarded for.
The effect of fuel burn-up on the radioactivation behavior of cladding hull materials was investigated using the ORIGEN-S code for various materials of Zircaloy-4, Zirlo, HANA-4, and HANA-6 and for various fuel burn-ups of 30, 45, 60, and 75 GWD/MTU. The Zircaloy-4 material is the only one that does not contain Nb as an alloy constituent, and it was revealed that 125Sb, 125mTe, and 55Fe are the major sources of radioactivity. On the other hand, 93mNb was identified as the most radioactive nuclide for the other materials although minor radioactive nuclides varied owing to their different initial constituents. The radioactivity of 94Nb was of particular focus owing to its acceptance limit against a Korean intermediate-/low-level waste repository. The radioactivation calculation results revealed that only Zircaloy-4 is acceptable for the Korean repository, while the other materials required at least 4,900 of Nb decontamination factor owing to the high radioactivity of 94Nb regardless of the fuel burn-up. A discussion was also made on the feasibility of Zr recovery methods (chlorination and electrorefining) for selective recovery of Zr so that it can be disposed of in the Korean repository. (author)
The double-heterogeneity characterising pebble-bed high temperature reactors (HTRs) makes Monte Carlo based calculation tools the most suitable for detailed core analyses. These codes can be successfully used to predict the isotopic evolution during irradiation of the fuel of this kind of cores. At the moment, there are many computational systems based on MCNP that are available for performing depletion calculation. All these systems use MCNP to supply problem dependent fluxes and/or microscopic cross sections to the depletion module. This latter then calculates the isotopic evolution of the fuel resolving Bateman's equations. In this paper, a comparative analysis of three different MCNP-based depletion codes is performed: Montburns2.0, MCNPX2.6.0 and BGCore. Monteburns code can be considered as the reference code for HTR calculations, since it has been already verified during HTR-N and HTR-N1 EU project. All calculations have been performed on a reference model representing an infinite lattice of thorium-plutonium fuelled pebbles. The evolution of k-inf as a function of burnup has been compared, as well as the inventory of the important actinides. The k-inf comparison among the codes shows a good agreement during the entire burnup history with the maximum difference lower than 1%. The actinide inventory prediction agrees well. However significant discrepancy in Am and Cm concentrations calculated by MCNPX as compared to those of Monteburns and BGCore has been observed. This is mainly due to different Am-241 (n,γ) branching ratio utilized by the codes. The important advantage of BGCore is its significantly lower execution time required to perform considered depletion calculations. While providing reasonably accurate results BGCore runs depletion problem about two times faster than Monteburns and two to five times faster than MCNPX.
ABRAC: A microcomputer-based Fortran code for multi-cyle burnup
Pressurized-water reactors have reactor physics and fuel management characteristics which are very amenable to simplified analysis. Given models which account for the dominant features of core and fuel performance, it is possible to rapidly perform relatively accurate scoping studies of many years of reactor operation in just a few hours on a modern (386-class) microcomputer. Models are described for burnup-dependent cross-section generation, for burnup of fuel under irradiation, and for computation of radial power distributions in hexagonal geometry assuming hexagonal fuel assemblies. Comparisons with more elaborate methods are given in order to validate the models, and to quantify the accuracy of the results. 16 refs., 5 figs., 5 tabs
Burnup simulations of different fuel grades using the MCNPX Monte Carlo code
Asah-Opoku Fiifi
2014-01-01
Full Text Available Global energy problems range from the increasing cost of fuel to the unequal distribution of energy resources and the potential climate change resulting from the burning of fossil fuels. A sustainable nuclear energy would augment the current world energy supply and serve as a reliable future energy source. This research focuses on Monte Carlo simulations of pressurized water reactor systems. Three different fuel grades - mixed oxide fuel (MOX, uranium oxide fuel (UOX, and commercially enriched uranium or uranium metal (CEU - are used in this simulation and their impact on the effective multiplication factor (Keff and, hence, criticality and total radioactivity of the reactor core after fuel burnup analyzed. The effect of different clad materials on Keff is also studied. Burnup calculation results indicate a buildup of plutonium isotopes in UOX and CEU, as opposed to a decline in plutonium radioisotopes for MOX fuel burnup time. For MOX fuel, a decrease of 31.9% of the fissile plutonium isotope is observed, while for UOX and CEU, fissile plutonium isotopes increased by 82.3% and 83.8%, respectively. Keff results show zircaloy as a much more effective clad material in comparison to zirconium and stainless steel.
PLUTON is a three-group neutronic code analyzing, as functions of time and burnup, the change of radial profiles, together with average values, of power density, burnup, concentration of trans-uranium elements, plutonium buildup, depletion of fissile elements, and fission product generation in water reactor fuel rod with standard UO2, UO2-Gd2O3, inhomogeneous MOX, and UO2-ThO2. The PLUTON code, which has been designed to be run on Windows PC, has adopted a theoretical shape function of neutron attenuation in pellet, which enables users to perform a very fast and accurate calculation easily. The present code includes the irradiation conditions of the Halden Reactor which gives verification data for the code. The total list of trans-uranium elements included in the calculations consists of 92U233-239, 93Np237-239, 94Pu238-243, 95Am241-244 (including isomers), and 96Cm242-245. Poisoning fission products are represented by 54Xe131,133,135, 48Cd113, 62Sm149,151,152, 64Gd154-160, 63Eu153,155, 36Kr83,85, 42Mo95, 43Tc99, 45Rh103, 47Ag109, 53I127,129,131, 55Cs133, 57La139, 59Pr141, 60Nd143-150, 61Pm147. Fission gases and volatiles included in the code are 36Kr83-86, 54Xe129-136, 52Te125-130, 53I127-131, 55Cs133-137, and 56Ba135-140. Verification has been performed up to 83 GWd/tU, and a satisfactory agreement has been obtained. (author)
Implementation and adaption of the Computer Code ECOSYS/EXCEL for Austria as OECOSYS/EXCEL
During 1989, under contract to the Austrian Chamber of the Federal Chancellor, department VII, the radioecological forecast model OECOSYS was implemented by the Austrian Research Centre Seibersdorf on a VAX computer using VAX Fortran. OECOSYS allows the prediction of the consequences after a large scale contamination event. During 1992, under contract to the Austrian Federal Ministry of Health, Sports and Consumer Protection, department III OECOSYS - in the version of 1989 - was implemented on PC's in Seibersdorf and the Ministry using OS/2 and Microsoft -Fortran. In March 1993, the Ministry ordered an update which had become necessary and the evaluation of two exercise scenarios. Since that time the prognosis model with its auxiliary program and communication facilities is kept on stand-by and yearly exercises are performed to maintain its readiness. The current report describes the implementation and adaption to Austrian conditions of the newly available EXCEL version of the German ECOSYS prognosis model as OECOSYS. (author)
The aesthetics of code : on excellence in instrumental action
Pineiro, Erik
2003-01-01
Software systems form an essential part of Western society,serving as tools to uphold institutions, processes andservices. It is understandable, therefore, that the mostfundamental aspects of programs are their function and utility.But they are not, however, the only things programmers areconcerned with when writing them. On the contrary, programmers also discuss about many otheraspects of software, including the beauty of code. Theydistinguish between different programming styles and express...
Dieudonne, C.; Dumonteil, E.; Malvagi, F.; Diop, C. M. [Commissariat a l' Energie Atomique et aux Energies Alternatives CEA, Service d' Etude des Reacteurs et de Mathematiques Appliquees, DEN/DANS/DM2S/SERMA/LTSD, F91191 Gif-sur-Yvette cedex (France)
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)
This paper describes results of testing of the TRANSURANUS burn-up model (TUBRNP routine) for Gd-doped WWER-1000 fuel pin based on results of HELIOS code. The testing covers the analysis of different types of nuclear fuel rods from a neutronic point of view that one can encounter in the VVER-1000 reactor core. The HELIOS computations simulate the assembly geometry, and combine 4 different 235U enrichment configurations with 4 different Gd2O3-concentrations. For each of these combinations the radial distribution of the concentrations of 155Gd and 157Gd compute in one Gd-doped rod. Based on these results the recommendations on using cross section of Gd in TRANSURANUS TUBRNP model were proposed. (author)
Highlights: → The specifications required for the analyses of the destructive assay data taken from irradiated fuel in Ohi-1 and Ohi-2 PWRs were documented in this paper. → These data were analyzed using the SWAT2.1 code, and the calculation results showed good agreement with experimental results. → These destructive assay data are suitable for the benchmarking of the burnup calculation code systems. - Abstract: The isotopic composition of spent nuclear fuels is vital data for studies on the nuclear fuel cycle and reactor physics. The Japan Atomic Energy Agency (JAEA) has been active in obtaining such data for pressurized water reactor (PWR) and boiling water reactor (BWR) fuels, and some data has already been published. These data have been registered with the international Spent Fuel Isotopic Composition Database (SFCOMPO) and widely used as international benchmarks for burnup calculation codes and libraries. In this paper, Assay Data of Spent Nuclear Fuel from two fuel assemblies irradiated in the Ohi-1 and Ohi-2 PWRs in Japan are shown. The destructive assay data from Ohi-2 have already been published. However, these data were not suitable for the benchmarking of calculation codes and libraries because several important specifications and data were not included. This paper summarizes the details of destructive assay data and specifications required for analyses of isotopic composition from Ohi-1 and Ohi-2. For precise burnup analyses, the burnup values of destructive assay samples were re-evaluated in this study. These destructive assay data were analyzed using the SWAT2.1 code, and the calculation results showed good agreement with experimental results. This indicates that the quality of destructive assay data from Ohi-1 and Ohi-2 PWRs is high, and that these destructive assay data are suitable for the benchmarking of burnup calculation code systems.
Lemehov, Sergei E; Suzuki, Motoe [Japan Atomic Energy Research Inst., Tokai, Ibaraki (Japan). Tokai Research Establishment
2001-08-01
PLUTON is a three-group neutronic code analyzing, as functions of time and burnup, the change of radial profiles, together with average values, of power density, burnup, concentration of trans-uranium elements, plutonium buildup, depletion of fissile elements, and fission product generation in water reactor fuel rod with standard UO{sub 2}, UO{sub 2}-Gd{sub 2}O{sub 3}, inhomogeneous MOX, and UO{sub 2}-ThO{sub 2}. The PLUTON code, which has been designed to be run on Windows PC, has adopted a theoretical shape function of neutron attenuation in pellet, which enables users to perform a very fast and accurate calculation easily. The present code includes the irradiation conditions of the Halden Reactor which gives verification data for the code. The total list of trans-uranium elements included in the calculations consists of {sub 92}U{sup 233-239}, {sub 93}Np{sup 237-239}, {sub 94}Pu{sup 238-243}, {sub 95}Am{sup 241-244} (including isomers), and {sub 96}Cm{sup 242-245}. Poisoning fission products are represented by {sub 54}Xe{sup 131,133,135}, {sub 48}Cd{sup 113}, {sub 62}Sm{sup 149,151,152}, {sub 64}Gd{sup 154-160}, {sub 63}Eu{sup 153,155}, {sub 36}Kr{sup 83,85}, {sub 42}Mo{sup 95}, {sub 43}Tc{sup 99}, {sub 45}Rh{sup 103}, {sub 47}Ag{sup 109}, {sub 53}I{sup 127,129,131}, {sub 55}Cs{sup 133}, {sub 57}La{sup 139}, {sub 59}Pr{sup 141}, {sub 60}Nd{sup 143-150}, {sub 61}Pm{sup 147}. Fission gases and volatiles included in the code are {sub 36}Kr{sup 83-86}, {sub 54}Xe{sup 129-136}, {sub 52}Te{sup 125-130}, {sub 53}I{sup 127-131}, {sub 55}Cs{sup 133-137}, and {sub 56}Ba{sup 135-140}. Verification has been performed up to 83 GWd/tU, and a satisfactory agreement has been obtained. (author)
The reactivity effect of the asymmetry of axial burnup profile in burnup credit criticality safety is studied for a realistic PWR spent fuel transport cask proposed in the current OECD/NEA Phase II-C benchmark problem. The axial burnup profiles are simulated in 21 material zones based on in-core flux measurements varying from strong asymmetry to more or less no asymmetry. Criticality calculations in a 3-D model have been performed using the continuous energy Monte Carlo code MCNP-4B2 and the nuclear data library JENDL-3.2. Calculation conditions are determined with consideration of the axial fission source convergence. Calculations are carried out not only for cases proposed in the benchmark but also for additional cases assuming symmetric burnup profile. The actinide-only approach supposed for first domestic introduction of burnup credit into criticality evaluation is also considered in addition to the actinide plus fission product approach adopted in the benchmark. The calculated results show that keff and the end effect increase almost linearly with increasing burnup axial offset that is defined as one of typical parameters showing the intensity of axial burnup asymmetry. The end effect is more sensitive to the asymmetry of burnup profile for the higher burnup. For an axially distributed burnup, the axial fission source distribution becomes strongly asymmetric as its peak shifts toward the top end of the fuel's active zone where the local burnup is less than that of the bottom end. The peak of fission source distribution becomes higher with the increase of either the asymmetry of burnup profile or the assembly-averaged burnup. The conservatism of the assumption of uniform axial burnup based on the actinide-only approach is estimated quantitatively in comparison with the keff result calculated with experiment-based strongest asymmetric axial burnup profile with the actinide plus fission product approach. (author)
Borio Di Tigliole, A.; Bruni, J.; Panza, F. [Dept. of Nuclear and Theoretical Physics, Univ. of Pavia, 27100 Pavia (Italy); Italian National Inst. of Nuclear Physics INFN, Section of Pavia, Via A. Bassi, 6, 27100 Pavia (Italy); Alloni, D.; Cagnazzo, M.; Magrotti, G.; Manera, S.; Prata, M.; Salvini, A. [Italian National Inst. of Nuclear Physics INFN, Section of Pavia, Via A. Bassi, 6, 27100 Pavia (Italy); Applied Nuclear Energy Laboratory LENA, Univ. of Pavia, Via Aselli, 41, 27100 Pavia (Italy); Chiesa, D.; Clemenza, M.; Pattavina, L.; Previtali, E.; Sisti, M. [Physics Dept. G. Occhialini, Univ. of Milano Bicocca, 20126 Milano (Italy); Italian National Inst. of Nuclear Physics INFN, Section of Milano Bicocca, P.zza della Scienza, 3, 20126 Milano (Italy); Cammi, A. [Italian National Inst. of Nuclear Physics INFN, Section of Milano Bicocca, P.zza della Scienza, 3, 20126 Milano (Italy); Dept. of Energy Enrico Fermi Centre for Nuclear Studies CeSNEF, Polytechnic Univ. of Milan, Via U. Bassi, 34/3, 20100 Milano (Italy)
2012-07-01
Aim of this work was to perform a rough preliminary evaluation of the burn-up of the fuel of TRIGA Mark II research reactor of the Applied Nuclear Energy Laboratory (LENA) of the Univ. of Pavia. In order to achieve this goal a computation of the neutron flux density in each fuel element was performed by means of Monte Carlo code MCNP (Version 4C). The results of the simulations were used to calculate the effective cross sections (fission and capture) inside fuel and, at the end, to evaluate the burn-up and the uranium consumption in each fuel element. The evaluation, showed a fair agreement with the computation for fuel burn-up based on the total energy released during reactor operation. (authors)
Aim of this work was to perform a rough preliminary evaluation of the burn-up of the fuel of TRIGA Mark II research reactor of the Applied Nuclear Energy Laboratory (LENA) of the Univ. of Pavia. In order to achieve this goal a computation of the neutron flux density in each fuel element was performed by means of Monte Carlo code MCNP (Version 4C). The results of the simulations were used to calculate the effective cross sections (fission and capture) inside fuel and, at the end, to evaluate the burn-up and the uranium consumption in each fuel element. The evaluation, showed a fair agreement with the computation for fuel burn-up based on the total energy released during reactor operation. (authors)
Ramamoorthy, Karthikeyan
The main aim of this research is the development and validation of computational schemes for advanced lattice codes. The advanced lattice code which forms the primary part of this research is "DRAGON Version4". The code has unique features like self shielding calculation with capabilities to represent distributed and mutual resonance shielding effects, leakage models with space-dependent isotropic or anisotropic streaming effect, availability of the method of characteristics (MOC), burnup calculation with reaction-detailed energy production etc. Qualified reactor physics codes are essential for the study of all existing and envisaged designs of nuclear reactors. Any new design would require a thorough analysis of all the safety parameters and burnup dependent behaviour. Any reactor physics calculation requires the estimation of neutron fluxes in various regions of the problem domain. The calculation goes through several levels before the desired solution is obtained. Each level of the lattice calculation has its own significance and any compromise at any step will lead to poor final result. The various levels include choice of nuclear data library and energy group boundaries into which the multigroup library is cast; self shielding of nuclear data depending on the heterogeneous geometry and composition; tracking of geometry, keeping error in volume and surface to an acceptable minimum; generation of regionwise and groupwise collision probabilities or MOC-related information and their subsequent normalization thereof, solution of transport equation using the previously generated groupwise information and obtaining the fluxes and reaction rates in various regions of the lattice; depletion of fuel and of other materials based on normalization with constant power or constant flux. Of the above mentioned levels, the present research will mainly focus on two aspects, namely self shielding and depletion. The behaviour of the system is determined by composition of resonant
Application of two-dimensional burnup computer codes to the operation of nuclear power plants
The needs for three-dimensional computer calculations of the power density distribution in WWER type reactors are outlined. In most cases, however, two-dimensional calculations provide sufficiently exact results and result in a decrease in computer costs. The application, performance and computer codes of two-dimensional calculations are dealt with. (author)
Allowing Monte Carlo (MC) codes to perform fuel cycle calculations requires coupling to a point depletion solver. In order to perform depletion calculations, one-group (1-g) cross sections must be provided in advance. This paper focuses on generating accurate 1-g cross section values that are necessary for evaluation of nuclide densities as a function of burnup. The proposed method is an alternative to the conventional direct reaction rate tally approach, which requires extensive computational efforts. The method presented here is based on the multi-group (MG) approach, in which pre-generated MG sets are collapsed with MC calculated flux. In our previous studies, we showed that generating accurate 1-g cross sections requires their tabulation against the background cross-section (σ0) to account for the self-shielding effect. However, in previous studies, the model that was used to calculate σ0 was simplified by fixing Bell and Dancoff factors. This work demonstrates that 1-g values calculated under the previous simplified model may not agree with the tallied values. Therefore, the original background cross section model was extended by implicitly accounting for the Dancoff and bell factors. The method developed here reconstructs the correct value of σ0 by utilizing statistical data generated within the MC transport calculation by default. The proposed method was implemented into BGCore code system. The 1-g cross section values generated by BGCore were compared with those tallied directly from the MCNP code. Very good agreement (<0.05%) in the 1-g cross values was observed. The method dose not carry any additional computational burden and it is universally applicable to the analysis of thermal as well as fast reactor systems. (author)
Inclusion of historical dependences of fuel burn-up into MOBY-DICK Code
The paper briefly describes inclusion of historical dependences of cross sections into MOBY-DICK code. Changes in program and its library are specified, especially from the point of view of programs user. Preliminary testing on assembly level is described and also testing on core level for 'transient' loading pattern calculations is demonstrated on examples of 18th and 19th cycles of the Dukovany NPP Unit III. Some features are addressed in the end of the paper (Authors)
Incorporation of the variation in conductivity with burnup in the stability of code predictive LAPUR
In the field of nuclear safety, the analysis of the stability of boiling water reactors is one of the biggest challenges for researchers. LAPUR code that allows to obtain the parameters of stability of the plant (Decay rate and frequency), being one of the programs used by IBERDROLA can be used for these calculations. With the collaboration of the research group TIN of the Polytechnic University of Valencia, a model of loss of conductivity of uranium has joined with the burned LAPUR. This update allows you to play the phenomenon in a more realistic way. This improvement has been validated and verified contrasting results with reference values.
Highlights: • The inventory of the radioactive nuclides was calculated using the GETERA code. • The consumptions of 235U for the HEU, MEU and LEU were: 2360, 2334 and 2320 g. • The amounts of 239Pu produced were: 67.07, 157.86 g for the MEU and LEU. • The core radioactivity for the MEU and LEU were: 8.84 × 1016 and 9.31 × 1016 Bq. - Abstract: Efforts have been made recently to study the possibility of core conversion of the 10 MW MTR type research reactor from the HEU to LEU fuels due to the proliferation issue. Since the inventory of the reactor core is a required parameter to study the atmospheric dispersion calculation for a postulated accidental airborne radionuclide release from the reactor, the inventory of the radioactive nuclides accumulated in the UAlx–Al fuels: HEU (93% 235U), MEU (45% 235U) and LEU (20% 235U) after 200 days of the reactor operating time was calculated using the GETERA code. The result showed, after 200 days of the reactor operation time (35% burnup), that the total consumptions of 235U and 238U for the HEU, MEU and LEU fuels were: 2360, 2334 and 2320 g for the 235U and 13, 105 and 238 g for the 238U, respectively. The amounts of 239Pu produced in the core were: 67.07, 157.86 g for the MEU and LEU fuels, respectively, compared with 7.95 g. The total core radioactivity after 200 days for the MEU and LEU cores were: 8.84 × 1016 and 9.31 × 1016 Bq, respectively, compared with 8.63 × 1016 Bq for the HEU core
Gholamzadeh Zohreh; Hossein Feghhi Seyed Amir; Soltani Leila; Rezazadeh Marzieh; Tenreiro Claudio; Joharifard Mahdi
2014-01-01
Decrease of the economically accessible uranium resources and the inherent proliferation resistance of thorium fuel motivate its application in nuclear power systems. Estimation of the nuclear reactor’s neutronic parameters during different operational situations is of key importance for the safe operation of nuclear reactors. In the present research, thorium oxide fuel burn-up calculations for a demonstrative model of a heavy water- -cooled reactor have been performed using MCNPX 2.6 code. N...
Gholamzadeh Zohreh
2014-12-01
Full Text Available Decrease of the economically accessible uranium resources and the inherent proliferation resistance of thorium fuel motivate its application in nuclear power systems. Estimation of the nuclear reactor’s neutronic parameters during different operational situations is of key importance for the safe operation of nuclear reactors. In the present research, thorium oxide fuel burn-up calculations for a demonstrative model of a heavy water- -cooled reactor have been performed using MCNPX 2.6 code. Neutronic parameters for three different thorium fuel matrices loaded separately in the modelled thermal core have been investigated. 233U, 235U and 239Pu isotopes have been used as fissile element in the thorium oxide fuel, separately. Burn-up of three different fuels has been calculated at 1 MW constant power. 135X and 149Sm concentration variations have been studied in the modelled core during 165 days burn-up. Burn-up of thorium oxide enriched with 233U resulted in the least 149Sm and 135Xe productions and net fissile production of 233U after 165 days. The negative fuel, coolant and void reactivity of the used fuel assures safe operation of the modelled thermal core containing (233U-Th O2 matrix. Furthermore, utilisation of thorium breeder fuel demonstrates several advantages, such as good neutronic economy, 233U production and less production of long-lived α emitter high radiotoxic wastes in biological internal exposure point of view
Analysis of the burnup of the control rods with the COREMASTER-Presto code
An evaluation of the capacity of the COREMASTER-Presto code, to evaluate generically the burnt of the control bars in the Laguna Verde reactors plant (CLV) is made. It was found that the code only reports burnt values of the control rods in MWD/TM, in spite of having with a second order polynomial model, for the conversion to remainder of the Boron-10 (B-10). It was observed that said model is adequate only for burnt smaller to 45,000 MWD/TM. To evaluate the burnt of the control rods it was reproduced the balance cycle of 18 months for the CLV, executing Cm-Presto during 13 consecutive cycles. First without rod burnt, taking this as the base case. Later on, cases with 1, 2 and up to 13 cycles with rod burnt were generated. When comparing results it was observed that the control rods pattern it loses reactivity lineally with the burnt one. By each 10 G Wd/T of burnt of the nucleus it is decreased the reactivity of the pattern rods ∼ 1 pcm in hot condition and of ∼ 20 pcm in cold condition. When burning three cycles those rods more burnt reached the 13,900 MWD/TM, equivalent to 36% of B-10 reduction, near value to 34% proposed by aging in the one lost study of B-10. It was observed that Cm-Presto it doesn't burn the superior node of the control rods when these are completely extracted. A one big lost of B-10, of the order of 50%, it represents only a decrease of 11% of the reactivity value of the rod. One can affirm that even when it is strongly decreased the content of B-10, the rod is continue considering as a black absorber, that is to say, thermal neutron that enters in the neutron rod that is absorbed. (Author)
Des macros Excel pour exporter et importer des modules de code VBA
Duclos, Jean Baptiste
2008-01-01
Cet article présente un outil Excel qui permet de gérer facilement l’exportation et l’importation de macros VBA, en sauvegardant les programmes des macros dans des fichiers textes. Les macros d’importation et d’exportation de code peuvent être facilement utilisées dans d’autres programmes.
SFR whole core burnup calculations with TRIPOLI-4 Monte Carlo code
Under the Working Party on Scientific Issues of Reactor Systems (WPRS) of the OECD/NEA, an international collaboration benchmark was recently established on the neutronic analysis of four Sodium-cooled Fast Reactor (SFR) concepts of the Generation- IV nuclear energy systems. As the whole core Monte Carlo depletion calculation is one of the essential challenges of current reactor physics studies, the continuous-energy TRIPOLI-4 Monte Carlo transport code was firstly used in this study to perform whole core 3D neutronic calculations for these four SFR cores. Two medium size (1000 MWt) and two large size (3600 MWt) SFR of GEN-IV systems were analyzed. The medium size SFR concepts are from the Advanced Burner Reactors (ABR). The large size SFR concepts are from the self-breeding reactors. The TRIPOLI-4 depletion calculations were made with MOX and metallic U-Pu-Zr fuels for the ABR cores and with MOX and Carbide (U,Pu)C fuels for the self-breeding cores. The whole core reactor physics parameters calculations were performed for the BOEC and EOEC (Beginning and End of Equilibrium Cycle) conditions. This paper summarizes the TRIPOLI-4 calculation results of Keff, βeff, sodium void worth, Doppler constant, control rod worth, and core power distributions for the BOEC and EOEC conditions. The one-cycle depletion calculation results of the core inventory of U and TRU (Pu, Am, Cm, and Np) are also analyzed, after 328.5 days depletion irradiation for the ABR cores, 410 days for the large MOX core, and 500 days for the large carbide core. (author)
Results of the RTOP-CA code calculations for experiments in the research MIR reactor are presented. The MIR-reactor tests were made to study the activity release from defective WWER fuel at high burnup (∼60 MWd/kgU). The RTOP-CA calculations are compared to experimental data on radial distributions of burnup as well as radial profiles of Pu and Xe concentrations in fuel pellets. The RTOP-CA predictions are also compared to the data on activity release (radionuclides of I, Cs, Xe and Kr) from the test fuel rod with an artificial defect in cladding. Additional calculations were performed for WWER-1000 fuel of an advanced design. In these calculation series the effect of design innovations on activity release from defective fuel rods was estimated. It is demonstrated that in case of a failure the new generation of WWER fuel shows lower levels of activity release into primary coolant. (authors)
A fuel performance code FEMAXI-6 has been developed for the analysis of LWR fuel rod behaviors in normal operation and transient (not accident) conditions. The code uses FEM for mechanical analysis, and has incorporated thermal and mechanical models of phenomena anticipated in high burn-up fuel rods, such as fuel thermal conductivity degradation and pellet-clad bonding. In the present study, PCMI induced by swelling in a high burn-up BWR type fuel rod has been analyzed by the FEMAXI-6 code. During a power ramp for the high burn-up fuel, instantaneous pellet swelling can significantly exceed the level that is predicted by a 'steady-rate' swelling model, causing a large circumferential strain in cladding. This phenomenon has been simulated by a new swelling model to take into account the fission gas bubble growth, and we found that the new model can give satisfactory predictions on cladding diametral expansion in comparison with measurements in test rod. The bubble growth model assumes an equilibrium between bubble size and external pressure on the bubble, and simultaneous solution is obtained with both bubble size determination and diffusion equation of fission gas atoms. In addition, a pellet-clad bonding model which assumes firm mechanical coupling between pellet outer surface and cladding inner surface predicts an elevated tensile stress in the axial direction of cladding during ramp, indicating the generation of bi-axial stress state in the cladding. These analyses by the FEMAXI-6 code enable us to predict the magnitude of mechanical loading on cladding during transient and also serve for failure investigation. Clearly, prediction by code calculation depends on the creep and stress-strain properties of highly irradiated cladding. (author)
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)
This paper is focused on the description and validation of the new dynamic model aimed at the processes taking place on grain boundaries in polycrystalline Light Water Reactor fuel based on the use of UO2, which has been recently developed for the START-3 code. The analysis embraces such processes as formation of fine surface clusters and larger intergranular pores, equi-axial grain growth, direct release and percolation of fission gas to the open surfaces. This model is also closely linked to intragranular behaviour of fission gas and essentially overlaps modelling of High Burnup Structure, as well as high temperature- and power transient-assisted processes. The model embodies some of the state-of-the-art approaches to numerical description of the processes taking place on grain boundaries, incorporating considerations of diffusion theory with respect to analysis of the dynamics of grain boundary pore growth/shrinkage caused by self-diffusion of the fuel material around them. Besides, it widely uses the elements of probability theory thereby accounting for stochastic nature of the analyzed phenomena. The several examples of model validation, illustrating credibility of pertinent results as applied to a wide enough range of application, including high-temperature out-of-pile annealing, High Burnup Structure Effects and transient behaviour of Light Water Reactor fuel, are also presented. As the validation shows, the developed model can be accepted as an important element of overall dynamic modelling with a view to justification of reliability of high burnup Light Water Reactor fuel, and safety analysis, as well
Validation of lattice code 'EXCEL' with TIC experiments on uniform and regularly perturbed lattices
Temporary International Collective (TIC) was established in 1972 by an agreement among seven countries, namely, Bulgaria, Czechoslovakia, Germany, Hungary, Poland, Romania and Union of Soviet Socialist Republics. The main objective of TIC was to provide the experimental data for the reactor physics analysis of water cooled and water moderated power reactors (WWER). Extensive experimental work for different core configurations was carried out by TIC countries to investigate the physics behaviour of WWER lattices and the results were published in TIC volumes. Two VVER-1000 MWe reactors are currently in an advanced stage of construction and due for commissioning in Kudankulam, Tamil Nadu, India. Indigenous development of in-core fuel management computer codes for the analysis of hexagonal lattice cores is also in an advanced stage to address various design, operation and safety issues of VVER type cores. The validation of the above TIC lattice experiments will help in the identification of deficiencies in reactor physics design computational codes and the associated nuclear data libraries. In this paper, TIC experiments on uniform and regularly perturbed lattices have been analyzed as part of the validation of indigenous computer codes, EXCEL, TRIHEX-FA and HEXPIN developed at Light Water Reactors Physics Section, B.A.R.C. Neutron-nuclear multi-group cross-section libraries in WIMS/D format in 69/172 energy groups have been released by IAEA at the conclusion of WIMS library update project (WLUP). In the present study we have used libraries based on ENDF/B-6, ENDF/B-7, JEFF3.1 and JENDL3.2 evaluated nuclear datasets. The results of the theoretical analyses bring out the performance of the code system and various cross-section libraries.
The bundle-duct interaction analysis code ''BAMBOO'' has been developed for the purpose of predicting deformation of a wire-wrapped fuel pin bundle of a fast breeder reactor (FBR). The BAMBOO code calculates helical bowing and oval-distortion of all the fuel pins in a fuel subassembly. We developed deformation models in order to precisely analyze the irradiation induced deformation by the code: a model to analyze fuel pin self-bowing induced by circumferential gradient of void swelling as well as thermal expansion, and a model to analyze dispersion of the orderly arrangement of a fuel pin bundle. We made deformation analyses of high burn-up fuel subassemblies in Phenix reactor and compared the calculated results with the post irradiation examination data of these subassemblies for the verification of these models. From the comparison we confirmed that the calculated values of the oval-distortion and bowing reasonably agreed with the PIE results if these models were used in the analysis of the code. (author)
The deep burn fuel cycle for the incineration of military plutonium in the GT-MHR is studied using the Monte-Carlo burnup code. The irradiation is DF is so rich in fissile isotopes that the TF cannot guarantee a negative reactive feedback, and the presence of erbium as burnable poison is absolutely necessary for the reactivity safety reasons. At beginning of life (BOL) the fuel composed of DF, consisting of fresh military plutonium, after an irradiation period of three years the fuel is reprocessed into post driver fuel (PDF). The mass flow of the GT-MHR fuelled by military plutonium at the equilibrium of the fuel composition shows that 66% of 239Pu is burned in three years and 92% in six years. (authors)
The isotopic composition calculations were performed for 26 spent fuel samples from the Obrigheim PWR reactor and 55 spent fuel samples from 7 PWR reactors using the SAS2H module of the SCALE4.4 code system with 27, 44 and 238 group cross-section libraries and the SWAT code system with the 107 group cross-section library. For the analyses of samples from the Obrigheim PWR reactor, geometrical models were constructed for each of SCALE4.4/SAS2H and SWAT. For the analyses of samples from 7 PWR reactors, the geometrical model already adopted in the SCALE/SAS2H was directly converted to the model of SWAT. The four kinds of calculation results were compared with the measured data. For convenience, the ratio of the measured to calculated values was used as a parameter. When the ratio is less than unity, the calculation overestimates the measurement, and the ratio becomes closer to unity, they have a better agreement. For many important nuclides for burnup credit criticality safety evaluation, the four methods applied in this study showed good coincidence with measurements in general. More precise observations showed, however: (1) Less unity ratios were found for Pu-239 and -241 for selected 16 samples out of the 26 samples from the Obrigheim reactor (10 samples were deselected because their burnups were measured with Cs-137 non-destructive method, less reliable than Nd-148 method the rest 16 samples were measured with); (2) Larger than unity ratios were found for Am-241 and Cm-242 for both the 16 and 55 samples; (3) Larger than unity ratios were found for Sm-149 for the 55 samples; (4) SWAT was generally accompanied by larger ratios than those of SAS2H with some exceptions. Based on the measured-to-calculated ratios for 71 samples of a combined set in which 16 selected samples and 55 samples were included, the correction factors that should be multiplied to the calculated isotopic compositions were generated for a conservative estimate of the neutron multiplication factor
Reactivity effects of nonuniform axial burnup distributions on spent fuel
When conducting future criticality safety analyses on spent fuel shipping casks, burnup credit may play a significant role in determining the number of fuel assemblies that can be safely loaded into each cask. An important area in burnup credit analysis is the burnup variation along the length of the fuel assembly, which is determined by the location of the assembly in the reactor core and its residence time. A study of the effects of axial burnup distributions on reactivity has been conducted, using data from existing power plant fuel. Utilizing a one-dimensional, two-group diffusion code, named REALAX, the reactivity effects of axial burnup profiles have been calculated for various PWR fuel assemblies. The reactivity effects calculated by the code are defined in terms of k for the axially dependent burnup distribution minus k for a uniform axial burnup distribution at the assembly average burnup divided by k for a uniform axial burnup distribution at the assembly average burnup. Criticality safety specialists can take advantage of the quick-running code to determine axial effects of different assembly burnup profiles. In general, the positive reactivity effects of axial burnup distributions increase as burnup increases, though they do not increase faster than the overall decrease in reactivity due to burnup
Reactivity effects of nonuniform axial burnup distributions on spent fuel
Leary, R.W. II; Parish, T.A. [Texas A & M Univ., College Station, TX (United States)
1995-12-01
When conducting future criticality safety analyses on spent fuel shipping casks, burnup credit may play a significant role in determining the number of fuel assemblies that can be safely loaded into each cask. An important area in burnup credit analysis is the burnup variation along the length of the fuel assembly, which is determined by the location of the assembly in the reactor core and its residence time. A study of the effects of axial burnup distributions on reactivity has been conducted, using data from existing power plant fuel. Utilizing a one-dimensional, two-group diffusion code, named REALAX, the reactivity effects of axial burnup profiles have been calculated for various PWR fuel assemblies. The reactivity effects calculated by the code are defined in terms of k for the axially dependent burnup distribution minus k for a uniform axial burnup distribution at the assembly average burnup divided by k for a uniform axial burnup distribution at the assembly average burnup. Criticality safety specialists can take advantage of the quick-running code to determine axial effects of different assembly burnup profiles. In general, the positive reactivity effects of axial burnup distributions increase as burnup increases, though they do not increase faster than the overall decrease in reactivity due to burnup.
Supercell burnup model for the physics design of BWR fuel assemblies
A code called SUPERB has been developed for the BWR fuel assembly burnup analyses using supercell model. Each of the characteristic heterogeneities of a BWR fuel assembly like water gap, poisoned pins, control blade etc., is treated by invoking appropriate supercell concept. The burnup model of SUPERB is so devised as to strike a balance between accuracy and speed. This is achieved by building isotopic densities in each fuel pin separately while the depletion equations are solved only in a few groups of pins or burnup zones and the multigroup neutron spectra are differentiated in fewer group of pincell types. Multiple fuel ring burnup is considered only for Gd isotopes. A special empirical formula allows the microscopic cross section of Gd isotopes to be varied even during burnup integration. The supercell model has been tested against Monte Carlo results for the fresh cold clean Tarapur fuel assembly with two Gd fuel pins. The burnup model of SUPERB has been validated against one of the most sophisticated codes LWR-WIMS for a benchmark problem involving all the complexities of a BWR fuel assembly. The agreement of SUPERB results with both Monte Carlo and LWR-WIMS results is found to be excellent. (auth.)
Paladino, Patricia Andrea
2006-07-01
Computational program are used for thermal hydraulic analysis of accidents and transients conditions in nuclear power plants. The RELAP5 code has been developed to simulate accidents and transients conditions, performing a best estimate analysis, in Pressurized Water Reactors (PWR) and auxiliary systems. The RELAP5 code, which has been used as a toll for licensing nuclear facilities in Brazil, is the objective of the study performed in this work. The main problem in using the RELAP5 code is the huge amount of information necessary to model the nuclear reactor and thus to simulate thermal-hydraulic accidents. Moreover, the RELAP5 code input data requires a large amount of mathematical operations to calculate the geometry of the plant components. Therefore, in order to make easier the data input for the RELAP5 code a friendly preprocessor has been developed. The preprocessor accepts basic information about the geometry of the plant components and performs all the calculations needed for the RELAP5 input. This preprocessor has been developed based on the MS-Excel software. (author)
Burnup study of 18 months and 16/20 months cycle AP1000 cores using CASMO4E and SIMULATE-3 codes
AP1000 reactor is an advanced pressurized water reactor equipped with passive safety systems. AP1000 reactor core is designed for 18 month cycle length and can also be used for 16/20 month alternate cycle lengths to meet energy requirements during high demand periods. The purpose of this study is to analyze the feasibility of AP1000 core for both 18 and 16/20 alternate cycle lengths by using CASMO4E and SIMULATE-3 code package. For this purpose, burnup analysis of both the schemes is carried out from initial core loading through optimized transition cores to equilibrium core. The study is performed by modeling three dimensional full core in SIMULATE-3 with each fuel assembly divided into 40 axial and 4 radial quadrant nodes. Once and twice burned fuel reloading from one cycle to the next and removal of burnable poison rods from the core after first cycle options are used in these codes. The results of this study indicate that both the cycle schemes can be utilized by varying the core loading pattern. Moreover, reactivity coefficients, total power peaking factors and enthalpy rise factors are calculated which indicate that the AP1000 core provide adequate safety margins in both the cycle schemes. (author)
Aiming at validation of depletion codes against WWER-440 spent fuel data some calculations of isotope composition of WWER-440 spent fuel assembly have been carried out by the NESSEL-NUKO code system. The initial data and data for the comparisons were taken from the ISTS burn up credit project data, recently published in the ISTC report 'Radiochemical Assays of Irradiated WWER-440 Fuel for Use in Spent Fuel Burnup Credit Activities. The specific work scope included the destructive assay (DA) of spent fuel assembly rod segments with an - -38.5 MWd/KgU burn up from a single WWER-440 fuel assembly from the Novovorenezh reactor in Russia (Authors)
Activity ratio measurement and burnup analysis for high burnup PWR fuels
Applying burnup credit to spent fuel transportation and storage system is beneficial. To take burnup credit to criticality safety design for a spent fuel transportation cask and storage rack, the burnup of target fuel assembly based on core management data must be confirmed by experimental methods. Activity ratio method, in which measured the ratio of the activity of a nuclide to that of another, is one of the ways to confirm burnup history. However, there is no public data of gamma-ray spectrum from high burnup fuels and validation of depletion calculation codes is not sufficient in the evaluation of the burnup or nuclide inventories. In this study, applicability evaluation of activity ratio method was carried out for high burnup fuel samples taken from PWR lead use assembly. In the gamma-ray measurement experiments, energy spectrum was taken in the Reactor Fuel Examination Facility (RFEF) of Japan Atomic Energy Agency (JAEA), and 134Cs/137Cs and 154Eu/137Cs activity ratio were obtained. With the MVP-BURN code, the activity ratios were calculated by depletion calculation tracing the operation history. As a result, 134Cs/137Cs and 154Eu/137Cs activity ratios for UO2 fuel samples show good agreements and the activity ratio method has good applicability to high burnup fuels. 154Eu/134Cs activity ratio for Gd2O3+UO2 fuels also shows good agreements between calculation results and experimental results as well as the activity ratio for UO2 fuels. It also becomes clear that it is necessary to pay attention to not only burnup but also axial burnup distribution history when confirming the burnup of UO2+Gd2O3 fuel with 134Cs/137Cs activity ratios. (author)
Analysis of high burnup fuel safety issues
Lee, Chan Bock; Kim, D. H.; Bang, J. G.; Kim, Y. M.; Yang, Y. S.; Jung, Y. H.; Jeong, Y. H.; Nam, C.; Baik, J. H.; Song, K. W.; Kim, K. S
2000-12-01
Safety issues in steady state and transient behavior of high burnup LWR fuel above 50 - 60 MWD/kgU were analyzed. Effects of burnup extension upon fuel performance parameters was reviewed, and validity of both the fuel safety criteria and the performance analysis models which were based upon the lower burnup fuel test results was analyzed. It was found that further tests would be necessary in such areas as fuel failure and dispersion for RIA, and high temperature cladding corrosion and mechanical deformation for LOCA. Since domestic fuels have been irradiated in PWR up to burnup higher than 55 MWD/kgU-rod. avg., it can be said that Korea is in the same situation as the other countries in the high burnup fuel safety issues. Therefore, necessary research areas to be performed in Korea were derived. Considering that post-irradiation examination(PIE) for the domestic fuel of burnup higher than 30 MWD/kgU has not been done so far at all, it is primarily necessary to perform PIE for high burnup fuel, and then simulation tests for RIA and LOCA could be performed by using high burnup fuel specimens. For the areas which can not be performed in Korea, international cooperation will be helpful to obtain the test results. With those data base, safety of high burnup domestic fuels will be confirmed, current fuel safety criteria will be re-evaluated, and finally transient high burnup fuel behavior analysis technology will be developed through the fuel performance analysis code development.
Analysis of high burnup fuel safety issues
Safety issues in steady state and transient behavior of high burnup LWR fuel above 50 - 60 MWD/kgU were analyzed. Effects of burnup extension upon fuel performance parameters was reviewed, and validity of both the fuel safety criteria and the performance analysis models which were based upon the lower burnup fuel test results was analyzed. It was found that further tests would be necessary in such areas as fuel failure and dispersion for RIA, and high temperature cladding corrosion and mechanical deformation for LOCA. Since domestic fuels have been irradiated in PWR up to burnup higher than 55 MWD/kgU-rod. avg., it can be said that Korea is in the same situation as the other countries in the high burnup fuel safety issues. Therefore, necessary research areas to be performed in Korea were derived. Considering that post-irradiation examination(PIE) for the domestic fuel of burnup higher than 30 MWD/kgU has not been done so far at all, it is primarily necessary to perform PIE for high burnup fuel, and then simulation tests for RIA and LOCA could be performed by using high burnup fuel specimens. For the areas which can not be performed in Korea, international cooperation will be helpful to obtain the test results. With those data base, safety of high burnup domestic fuels will be confirmed, current fuel safety criteria will be re-evaluated, and finally transient high burnup fuel behavior analysis technology will be developed through the fuel performance analysis code development
Burnup performances of boron nitride and boron coated nuclear fuels
The nuclear fuels of urania (UOV) and 5% and 10% gadolinia (Gd2O3) containing UO2 previously produced by sol-gel technique were coated with first boron nitride (BN) then boron (B) thin layer by chemical vapor deposition (CVD) and also by plasma enhanced chemical vapor deposition (PECVD) techniques to increase the fuel cycle length and to improve the physical properties. From the cross-sectional view of BN and B layers taken from scanning electron microscope (SEM), the excellent adherence of BN onto fuel and B onto BN layer was observed in both cases. The behavior of fuel burnup, depletion of BN and B, the effect of coating thickness and also Gd2O3 content on the burnup performances of the fuels were identified by using the code WIMS-D/4 for Pressurized Water Reactor (PWR) and Boiling Water Reactor (BWR) cores. The optimum thickness ratio of B to BN was found as 4 and their thicknesses were chosen as 40 mm and 10 mm respectively in both reactor types to get extended cycle length. The assemblies consisting of fuels with 5% Gd2O3 and also coated with 10 mm BN and 40 mm B layers were determined as candidates for getting higher burnup in both types of reactors
REBUS: A burnup credit experimental programme
An international programme called REBUS (REactivity tests for a direct evaluation of the Burn-Up credit on Selected irradiated LWR fuel bundles) for the investigation of the burn-up credit has been initiated by the Belgian Nuclear Research Center SCK-CEN and Belgonucleaire. At present it is sponsored by USNRC, EdF from France and VGB, representing German nuclear utilities. The programme aims to establish a neutronic benchmark for reactor physics codes. This benchmark would qualify the codes to perform calculations of the burn-up credit. The benchmark exercise will investigate the following fuel types with associated burn-up. 1. Reference absorber test bundle, 2. Fresh commercial PWR UO2 fuel, 3. Irradiated commercial PWR UO2 fuel (50 GWd/tM), 4. Fresh PWR UO2 fuel, 5. Irradiated PWR UO2 fuel (30 GWd/tM). Reactivity effects will be measured in the critical facility VENUS. The accumulated burn-up of all rods will be measured non-destructively by gamma-spectrometry. Some rods will be analyzed destructively with respect to accumulated burn-up, actinides content and TOP-18 fission products (i.e. those non-gaseous fission products that have most implications on the reactivity). The experimental implementation of the programme will start in 2000. (author)
Theory analysis and simple calculation of travelling wave burnup scheme
Travelling wave burnup scheme is a new burnup scheme that breeds fuel locally just before it burns. Based on the preliminary theory analysis, the physical imagine was found. Through the calculation of a R-z cylinder travelling wave reactor core with ERANOS code system, the basic physical characteristics of this new burnup scheme were concluded. The results show that travelling wave reactor is feasible in physics, and there are some good features in the reactor physics. (authors)
Highlights: • A successful validation of the burn-up simulation system EVOLCODE is presented here. • A Sensitivity/Uncertainty model was applied for uncertainty propagation/assessment. • Cross sections are for most cases the main contributors to inventory uncertainties. • The improved model helps to explain some simulation-experiment discrepancies. • Some hints for the improvement of basic data libraries are provided. - Abstract: A validation of the burn-up simulation system EVOLCODE 2.0 is presented here, involving the experimental measurement of U and Pu isotopes and some fission fragments production ratios after a burn-up of around 30 GWd/tU in a Pressurized Light Water Reactor (PWR). This work provides an in-depth analysis of the validation results, including the possible sources of the uncertainties. An uncertainty analysis based on the sensitivity methodology has been also performed, providing the uncertainties in the isotopic content propagated from the cross sections uncertainties. An improvement of the classical Sensitivity/Uncertainty (S/U) model has been developed to take into account the implicit dependence of the neutron flux normalization, that is, the effect of the constant power of the reactor. The improved S/U methodology, neglected in this kind of studies, has proven to be an important contribution to the explanation of some simulation-experiment discrepancies for which, in general, the cross section uncertainties are, for the most relevant actinides, an important contributor to the simulation uncertainties, of the same order of magnitude and sometimes even larger than the experimental uncertainties and the experiment-simulation differences. Additionally, some hints for the improvement of the JEFF3.1.1 fission yield library and for the correction of some errata in the experimental data are presented
COGEMA/TRANSNUCLEAIRE's experience with burnup credit
Facing a continuous increase in the fuel enrichments, COGEMA and TRANSNUCLEAIRE have implemented step by step a burnup credit programme to improve the capacity of their equipment without major physical modification. Many authorizations have been granted by the French competent authority in wet storage, reprocessing and transport since 1981. As concerns transport, numerous authorizations have been validated by foreign competent authorities. Up to now, those authorizations are restricted to PWR Fuel type assemblies made of enriched uranium. The characterization of the irradiated fuel and the reactivity of the systems are evaluated by calculations performed with well qualified French codes developed by the CEA (French Atomic Energy Commission): CESAR as a depletion code and APPOLO-MORET as a criticality code. The authorizations are based on the assurance that the burnup considered is met on the least irradiated part of the fuel assemblies. Besides, the most reactive configuration is calculated and the burnup credit is restricted to major actinides only. This conservative approach allows not to take credit for any axial profile. On the operational side, the procedures have been reevaluated to avoid misloadings and a burnup verification is made before transport, storage and reprocessing. Depending on the level of burnup credit, it consists of a qualitative (go/no-go) verification or of a quantitative measurement. Thus the use of burnup credit is now a common practice in France and Germany and new improvements are still in progress: extended qualifications of the codes are made to enable the use of six selected fission products in the criticality evaluations. (author)
The gas tag burnup changes investigated were limited to the three tags (Kr-78/Kr-80, Xe-126/Xe-129 and Kr-82/Kr-80) currently accepted as being the most desirable. Control rod tag burnup was significantly greater than fuel rod tag burnup. This occurs because control rods stay in the reactor longer and occupy positions of greater low-energy flux. Thus, minimum tag spacings were set by the control rods as 1.079 for Kr-78/Kr-80, 1.189 for Xe-126/Xe-129 and 1.134 for Kr-82/Kr-80
An advanced model GRSWEL-A for fission gas behavior and micro-structural evolutions in Light Water Reactor (LWR) fuels was developed for and embedded in the START-3 fuel performance code. This paper represents the physical basis and verification of the model with emphasis on analysis of High Burn-up Structure (HBS), which is generally ascribed to a so-called rim-layer of high burn-up fuel pellets. Specifically, the issues of microscopic polygonization, loss of matrix fission gas, growth of fuel porosity and fission gas release are highlighted. The effects of HBS on total fission gas release, temperature distribution in the pellet, pellet swelling and permanent strain of the cladding are considered in the appropriate section of the paper by means of comparative and sensitivity analysis with the use of both modeling and available experimental data. In all the cases, an accounting for the present effects is found to be an important integral part of thorough analysis of LWR fuel behavior. Aside from the description of current capabilities of modeling, some priority directions of further improvement are outlined. (author)
Automated generation of burnup chain for reactor analysis applications
This paper presents the development of an automated generation of a new burnup chain for reactor analysis applications. The JENDL FP Decay Data File 2011 and Fission Yields Data File 2011 were used as the data sources. The nuclides in the new chain are determined by restrictions of the half-life and cumulative yield of fission products or from a given list. Then, decay modes, branching ratios and fission yields are recalculated taking into account intermediate reactions. The new burnup chain is output according to the format for the SRAC code system. Verification was performed to evaluate the accuracy of the new burnup chain. The results show that the new burnup chain reproduces well the results of a reference one with 193 fission products used in SRAC. Further development and applications are being planned with the burnup chain code. (author)
Investigation of burnup credit implementation for BWR fuel
Burnup Credit allows considering the reactivity decrease due to fuel irradiation in criticality studies for the nuclear fuel cycle. Its implementation requires to carefully analyze the validity of the assumptions made to define the axial profile of the burnup and void fraction (for BWR), to determine the composition of the irradiated fuel and to compute the criticality simulation. In the framework of Burnup Credit implementation for BWR fuel, this paper proposes to investigate part of these items. The studies presented in this paper concern: the influence of the burnup and of the void fraction on BWR spent fuel content and on the effective multiplication factor of an infinite array of BWR assemblies. A code-to-code comparison for BWR fuel depletion calculations relevant to Burnup Credit is also performed. (authors)
High Burnup Fuel Performance and Safety Research
Bang, Je Keun; Lee, Chan Bok; Kim, Dae Ho (and others)
2007-03-15
The worldwide trend of nuclear fuel development is to develop a high burnup and high performance nuclear fuel with high economies and safety. Because the fuel performance evaluation code, INFRA, has a patent, and the superiority for prediction of fuel performance was proven through the IAEA CRP FUMEX-II program, the INFRA code can be utilized with commercial purpose in the industry. The INFRA code was provided and utilized usefully in the universities and relevant institutes domesticallly and it has been used as a reference code in the industry for the development of the intrinsic fuel rod design code.
Fission gas release modelling at high burnup
A large quantity of experimental data on fission gas release is now available in the public domain. It covers a wide variety of fuel types and burnups of up to more than 70 GWd/tU. This data, together with gas release measurements from British Energy's AGRs, has been used to build a comprehensive validation database for the fuel performance code ENIGMA. Validation of ENIGMA version 5.11 against this database has identified a requirement for model development to improve predictions at high burnup. A modified gas release model has been produced and tested. (author)
Burnup and plutonium distribution of WWER-440 fuel pin at extended burnup
The formation of rim region in LWR UO2 based nuclear fuel at high burnup is a common observation. This region has very high porosity due to excessive gas release. Such a region is also characterized by a significantly high plutonium concentration and high local burnup compared to the internal fuel region. Spatial distribution of these parameters has been incorporated with fuel behavior and performance analysis codes by using mostly empirical relations. Variation of these parameters depends on the neutron flux as well as neutron energy spectrum. Detailed neutronics analysis is necessary for the accurate prediction of these parameters. This study is performed by MCNP4B Monte Carlo code for the calculation of local neutron flux, ORIGEN2 for burnup and depletion calculations, and MONTEBURNS for coupling these codes. For the analysis, a typical WWER-440 fuel pin and surrounding water moderator are considered in a hexagonal pin cell. Fuel pin is divided into a number of radial segments. A relatively small mesh size is used at the region near the surface to reveal the rim effect. The variation of plutonium and local burnup are obtained for high burnup. Results are compared with existing experimental observations for WWER-440 fuel and other theoretical predictions
Review of Axial Burnup Distribution Considerations for Burnup Credit Calculations
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
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.
The status of development of burnup credit for criticality safety analyses in Spain is described in this paper. Ongoing activities in the country in this field, both national and international, are resumed. Burnup credit is currently being applied to wet storage of PWR fuel, and credit to integral burnable absorbers is given for BWR fuel storage. It is envisaged to apply burnup credit techniques to the new generation of transport casks now in the design phase. The analysis methodologies submitted for the analyses of PWR and BWR fuel wet storage are outlined. Analytical activities in the country are described, as well as international collaborations in this field. Perspectives for future research and development of new applications are finally resumed. (author)
Detailed Burnup Calculations for Testing Nuclear Data
Leszczynski, F.
2005-05-01
A general method (MCQ) has been developed by introducing a microscopic burnup scheme that uses the Monte Carlo calculated fluxes and microscopic reaction rates of a complex system and a depletion code for burnup calculations as a basis for solving nuclide material balance equations for each spatial region in which the system is divided. Continuous energy-dependent cross-section libraries and full 3D geometry of the system can be input for the calculations. The resulting predictions for the system at successive burnup time steps are thus based on a calculation route where both geometry and cross sections are accurately represented, without geometry simplifications and with continuous energy data, providing an independent approach for benchmarking other methods and nuclear data of actinides, fission products, and other burnable absorbers. The main advantage of this method over the classical deterministic methods currently used is that the MCQ System is a direct 3D method without the limitations and errors introduced on the homogenization of geometry and condensation of energy of deterministic methods. The Monte Carlo and burnup codes adopted until now are the widely used MCNP and ORIGEN codes, but other codes can be used also. For using this method, there is need of a well-known set of nuclear data for isotopes involved in burnup chains, including burnable poisons, fission products, and actinides. For fixing the data to be included in this set, a study of the present status of nuclear data is performed, as part of the development of the MCQ method. This study begins with a review of the available cross-section data of isotopes involved in burnup chains for power and research nuclear reactors. The main data needs for burnup calculations are neutron cross sections, decay constants, branching ratios, fission energy, and yields. The present work includes results of selected experimental benchmarks and conclusions about the sensitivity of different sets of cross
Results of the isotopic concentrations of VVER calculational burnup credit benchmark No. 2(CB2)
Results of the nuclide concentrations are presented of VVER Burnup Credit Benchmark No. 2(CB2) that were performed in The Nuclear Technology Center of Cuba with available codes and libraries. The CB2 benchmark specification as the second phase of the VVER burnup credit benchmark is summarized. The CB2 benchmark focused on VVER burnup credit study proposed on the 97' AER Symposium. The obtained results are isotopic concentrations of spent fuel as a function of the burnup and cooling time. The depletion point 'ORIGEN2' code and other codes were used for the calculation of the spent fuel concentration. (author)
Burn-up measurements coupling gamma spectrometry and neutron measurement
The need to apply for burn-up credit arises with the increase of the initial enrichment of nuclear fuel. When burn-up credit is used in criticality safety studies, it is often necessary to confirm it by measurement. For the last 10 years, CANBERRA has manufactured the PYTHON system for such measurements. However, the method used in the PYTHON itself uses certain reactor data to arrive at burn-up estimates. Based on R and D led by CEA and COGEMA in the framework of burn-up measurement for burn-up credit and safeguards applications, CANBERRA is developing the next generation of burn-up measurement device. This new product, named SMOPY, is able to measure burn-up of any kind of irradiated fuel assembly with a combination of gamma spectrometry and passive neutron measurements. The measurement data is used as input to the CESAR depletion code, which has been developed and qualified by CEA and COGEMA for burn-up credit determinations. In this paper, we explain the complementary nature of the gamma and neutron measurements. In addition, we draw on our previous experience from PYTHON system and from COGEMA La Hague to show what types of evaluations are required to qualify the SMOPY system, to estimate its uncertainties, and to detect discrepancies in the fuel data given by the reactor plant to characterize the irradiated fuel assembly. (authors)
Tang, Guoping; Mayes, Melanie A.; Parker, Jack C.; Jardine, Philip M.
2010-09-01
We implemented the widely used CXTFIT code in Excel to provide flexibility and added sensitivity and uncertainty analysis functions to improve transport parameter estimation and to facilitate model discrimination for multi-tracer experiments on structured soils. Analytical solutions for one-dimensional equilibrium and nonequilibrium convection dispersion equations were coded as VBA functions so that they could be used as ordinary math functions in Excel for forward predictions. Macros with user-friendly interfaces were developed for optimization, sensitivity analysis, uncertainty analysis, error propagation, response surface calculation, and Monte Carlo analysis. As a result, any parameter with transformations (e.g., dimensionless, log-transformed, species-dependent reactions, etc.) could be estimated with uncertainty and sensitivity quantification for multiple tracer data at multiple locations and times. Prior information and observation errors could be incorporated into the weighted nonlinear least squares method with a penalty function. Users are able to change selected parameter values and view the results via embedded graphics, resulting in a flexible tool applicable to modeling transport processes and to teaching students about parameter estimation. The code was verified by comparing to a number of benchmarks with CXTFIT 2.0. It was applied to improve parameter estimation for four typical tracer experiment data sets in the literature using multi-model evaluation and comparison. Additional examples were included to illustrate the flexibilities and advantages of CXTFIT/Excel. The VBA macros were designed for general purpose and could be used for any parameter estimation/model calibration when the forward solution is implemented in Excel. A step-by-step tutorial, example Excel files and the code are provided as supplemental material.
Effect of Self-Shielding on Burn-Up Calculation of ETRR-2 Reactor
There exist two approaches for burn-up calculation. The first on is to use cell parameters generated using cell calculation code at different degrees of burn-up. The other is to use microscopic cross sections with self-shielding in order to compensate for the variation of spectrum at different degree of burn-up. The effect of using different forms of self-shielding factors on burn-up calculation for ETRR-2 reactor has been determined. The results of the two approaches are inter-compared up to 50% burn-up
A Burnup Analysis of PBMR-400MWth Reactor Core
The purpose of this study is to analyze the burnup characteristics of 400MWth PBMR using Monte Carlo method. In the world, the deterministic method is widely used to model such that system but it still has a disadvantage which is not flexible in simulating the burnup cycle. Although this method applies some techniques to increase the accuracy of calculation results but it is necessary to model this system by a suitable computer code that can verify and validate the results of the deterministic method. A method which uses a Monte Carlo technique for simulating the burnup cycle was performed. A reactor physics computer code uses in this method is MONTEBURN 2.0 which accurately and efficiently computes the neutronic and material properties of the fuel cycle. MONTEBURN is a fully automated tool that links the MCNP Monte Carlo transport code with a radioactive decay and burnup code ORIGEN. In this model, the calculations are based on a detailed core modeling using MCNP. The fuel pebble is thoroughly modeled by introducing unit cell modeling for the graphite matrix and fuel kernels in the pebble. For the burnup model, a start-up core was studied with considering the movement of pebbles. By shifting down one layer at each discrete time step and inserting fresh fuel from the top, this cyclic calculation is continued until equilibrium burnup cycle is achieved. In this study, the time dependence of multiplication factor keff, the spatial dependence of flux profile, power distribution, burnup, and inventory of isotopes in the start up process are analyzed. The results will provide the basis data of the burnup process and be also utilized as the verified data to validate a compute code for PBMR core analysis which will be developed in near future
Analysis of burnup credit on spent fuel storage
Chemical analyses were carried out on high burnup UO2 (65 GWd/t) and MOX (45 GWd/t) spent fuel pins. Measured data of the composition of nuclides from 234U to 242Pu were used for evaluation of ORIGEN-2/82 code. Criticality calculations were executed for the casks which were being designed to store 52 BWR or 21 PWR spent fuel assemblies. The reactivity biases were evaluated for (1) axial and horizontal profiles of burnup, and void history (BWR), (2) operational histories such as control rod insertion history, BPR insertion history and others, and (3) calculational accuracy of ORIGEN-2/82 code on the composition of nuclides. Present evaluation shows that introduction of burnup credit has a substantial merit in criticality safety analysis of the cask, even if these reactivity biases are considered. The concept of equivalent uniform burnup was adapted for present reactivity bias evaluation and showed a possibility of simplifying the reactivity bias evaluation in burnup credit. Finally, adapting procedures of burnup credit such as the burnup meter were evaluated. (author)
Evaluation of RSG-GAS Core Management Based on Burnup Calculation
Evaluation of RSG-GAS Core Management Based on Burnup Calculation. Presently, U3Si2-Al dispersion fuel is used in RSG-GAS core and had passed the 60th core. At the beginning of each cycle the 5/1 fuel reshuffling pattern is used. Since 52nd core, operators did not use the core fuel management computer code provided by vendor for this activity. They use the manually calculation using excel software as the solving. To know the accuracy of the calculation, core calculation was carried out using two kinds of 2 dimension diffusion codes Batan-2DIFF and SRAC. The beginning of cycle burn-up fraction data were calculated start from 51st to 60th using Batan-EQUIL and SRAC COREBN. The analysis results showed that there is a disparity in reactivity values of the two calculation method. The 60th core critical position resulted from Batan-2DIFF calculation provide the reduction of positive reactivity 1.84 % Δk/k, while the manually calculation results give the increase of positive reactivity 2.19 % Δk/k. The minimum shutdown margin for stuck rod condition for manual and Batan-3DIFF calculation are -3.35 % Δk/k dan -1.13 % Δk/k respectively, it means that both values met the safety criteria, i.e <-0.5 % Δk/k. Excel program can be used for burn-up calculation, but it is needed to provide core management code to reach higher accuracy. (author)
Fuel burnup characteristics for the NRU research reactor
The driver fuel of the NRU research reactor at AECL, Chalk River is a low enriched uranium (LEU) fuel alloy of Al-61 wt% U3Si, consisting of particles of U3Si dispersed in a continuous aluminum matrix, with 19.8% U235 in uranium. This paper describes the burnup characteristics for this type of fuel in NRU, including the determination of fuel depletion using the neutronic simulation code TRIAD, comparisons between simulated and measured burnup values, and the regulatory licensing operational average fuel burnup limit. (author)
Fuel burnup characteristics for the NRU research reactor
Leung, T.C., E-mail: leungt@aecl.ca [Atomic Energy of Canada Limited, Chalk River, Ontario (Canada)
2013-07-01
The driver fuel of the NRU research reactor at AECL, Chalk River is a low enriched uranium (LEU) fuel alloy of Al-61 wt% U{sub 3}Si, consisting of particles of U{sub 3}Si dispersed in a continuous aluminum matrix, with 19.8% U235 in uranium. This paper describes the burnup characteristics for this type of fuel in NRU, including the determination of fuel depletion using the neutronic simulation code TRIAD, comparisons between simulated and measured burnup values, and the regulatory licensing operational average fuel burnup limit. (author)
Calculation study of TNPS spent fuel pool using burnup credit
Exampled by the spent fuel pool of TNPS which is consist of 2 × 5 fuel storage racks, the spent fuel high-density storage based on burnup credit (BUC) and related criticality safety issues were studied. The MONK9A code was used to analyze keff, of different enrichment fuels at different burnups. A reference loading curve was proposed in accordance with the system keff's changing with the burnup of different initially enriched nuclear fuels. The capacity of the spent fuel pool increases by 31% compared with the one that does not consider BUC. (authors)
TRIGA criticality experiment for testing burn-up calculations
Persic, Andreja; Ravnik, Matjaz; Zagar, Tomaz [Jozef Stefan Institute, Reactor Physics Division, Ljubljana (Slovenia)
1999-07-01
A criticality experiment with partly burned TRIGA fuel is described. 20 wt % enriched standard TRIGA fuel elements initially containing 12 wt % U are used. Their average burn-up is 1.4 MWd. Fuel element burn-up is calculated in 2-D four group diffusion approximation using TRIGLAV code. The burn-up of several fuel elements is also measured by reactivity method. The excess reactivity of several critical and subcritical core configurations is measured. Two core configurations contain the same fuel elements in the same arrangement as were used in the fresh TRIGA fuel criticality experiment performed in 1991. The results of the experiment may be applied for testing the computer codes used for fuel burn-up calculations. (author)
Program package for 2D burnup calculation
The program package for 2 dimension burnup calculation was developed for TRIGA Mark III reactor. The package consists of 3 modules: PRESIX, SIXTUS-2, and BURN; 1 library, and 2 input files. PRESIX module prepared cross sections for diffusion calculation. SIXTUS-2 module, a two dimensional diffusion code in hexagonal geometry, calculates keff, neutron fluxes and power distributions. BURN module performs the burnup of fuel elements and stored the result in the ELEM.DAT file. PRESIX.LIB is two group cross section library for major reactor core components prepared using WIMS-D4 code. PRES.INP, the first input file, reads information on reactor power and core loading pattern. ELEM.DAT, the second input file, is prepared for specific TRIGA reactor and dependent on operation history. To verify the reactor model and computational methods, the calculated excess reactivities were compared to the measurement. The results are in good agreement. (author)
In this report an analysis of WWER-440 fuel of the second generation supplied by Russian JSC TVEL for high burnup fuel cycle is presented. The certificated code START-3 is applied to modeling of fuel rod operation parameters. Reliability of high-burnup fuel on the base of 5-6 year operation is demonstrated. Special attention is paid to aspects for attainment of burnup 70 MWd/kgU, including experimental and fuel modeling support and fuel operation experience
Burnup effects of MOX fuel pincells in PWR - OECD/NEA burnup credit benchmark analysis -
The burnup effects were analyzed for various cases of MOX fuel pincells of fresh and irradiated fuels by using the HELIOS, MCNP-4/B, CRX and CDP computer codes. The investigated parameters were burnup, cooling time and combinations of nuclides in the fuel region. The fuel compositions for each case were provided by BNFL (British Nuclear Fuel Limited) as a part of the problem specification so that the results could be focused on the calculation of the neutron multiplication factor. The results of the analysis show that the largest saving effect of the neutron multiplication factor due to burnup credit is 30 %. This is mainly due to the consideration of actinides and fission products in the criticality analysis
Burnup-dependent cross section data for research reactors
Studies currently in progress consider research and test reactors which commonly have burnups of 50 atom percent in 235-U and may reach as high a 70 atom percent. At these levels of burnup changes in cross-section data with burnup become significant. Some preliminary studies of these effects lead to the development of a modified version of REBUS-2 which supports changes in cross-section data with burnup. This version of REBUS-2 allows for changes in the cross-section data only at each time sub-interval in the problem, and these cross-section changes for capture and fission are based on a least squares polynomial fit as a function of burnup. In this paper an attempt is made to evaluate the importance of burnup dependent data for the various isotopes and/or groups, and to assess the accuracy of this method by comparing the REBUS-2 results with results obtained from PDQ-7. The 10 MW IAEA benchmark problem has been selected for this study. A description of the reactor and the XY model can be found in the IAEA Guidebook. The EPRI-CELL4 code was used to generate burnup dependent cross section data for use with both REBUS-2 and PDQ-7. Cross-section data were generated at 10 time steps to a burnup of approximately 50 atom percent in 235-U. The agreement between the PDQ-7 results and the REBUS-2 results with fitted burnup dependent cross-section data are quite good. Burnup dependent cross sections are essential for accurate estimates of cycle lengths and reactivities, and low order polynomial fits of capture and fission data for selected isotopes and energy groups can provide this capability
Light a CANDLE. An innovative burnup strategy of nuclear reactors
CANDLE is a new burnup strategy for nuclear reactors, which stands for Constant Axial Shape of Neutron Flux, Nuclide Densities and Power Shape During Life of Energy Production. When this candle-like burnup strategy is adopted, although the fuel is fixed in a reactor core, the burning region moves, at a speed proportionate to the power output, along the direction of the core axis without changing the spatial distribution of the number density of the nuclides, neutron flux, and power density. Excess reactivity is not necessary for burnup and the shape of the power distribution and core characteristics do not change with the progress of burnup. It is not necessary to use control rods for the control of the burnup. This booklet described the concept of the CANDLE burnup strategy with basic explanations of excess neutrons and its specific application to a high-temperature gas-cooled reactor and a fast reactor with excellent neutron economy. Supplementary issues concerning the initial core and high burnup were also referred. (T. Tanaka)
Burnup credit methodology validation against WWER experimental data
A methodology for criticality safety analyses with burnup credit application has been developed for WWER spent fuel management facilities. This methodology is based on two worldwide used code systems: SCALE 4.4 for depletion and criticality calculations and NESSEL-NUKO - for depletion calculations. The methodology is in process of extensive validation for WWER applications. The depletion code systems NESSEL-NUKO and SCALE4.4 (control module SAS2H) have been validated on the basis of comparison with the calculated results obtained by other depletion codes for the CB2 Calculational Burnup Credit Benchmark. The validation of these code systems for WWER-440 and WWER-1000 spent fuel assembly depletion analysis based on comparisons with appropriate experimental data commenced last year. In this paper some results from burnup methodology validation against measured nuclide concentration given in the ISTC project 2670 for WWER-440 and from ORNL publication for WWER-1000 are presented. (authors)
Study on burn-up credit and minor actinide in post-irradiation analysis
Accuracy of burnup calculation for actinide is very important as to the study of burn-up credit. For minor-actinides such as Am243 and Cm244, however, typical burnup calculation codes are not accurate enough. The accuracy for both nuclides was studied by using the SWAT code. The study showed that the C/E values of both nuclides could be improved at the same time by changing the cross section of Pu242. A study of burnup calculation related to the cross section of Pu242 should be performed to improve the accuracy for both nuclides. (author)
Effect of fuel burnup history on neutronic characteristics of WWER-1000 core
The paper analyzes fuel burnup history effect on neutronic characteristics of WWER-1000 core with use of the DYN3D codes. The DYN3D code employs the local Pu-239 concentration as an indicator of burnup spectral history. The calculations have been performed for the first four fuel loadings of Khmelnitsky NPP unit 2 and stationary fuel loading with TVSA. The effect of fuel burnup history is shown both on macro-characteristics on the reactor core and on local values of burnup and power
The purpose of this report is to describe an updated set of reactor models for pressurized-water reactors (PWRs) and boiling-water reactors (BWRs) operating on uranium fuel cycles and the methods used to generate the information for these models. Since new fuel cycle schemes and reactor core designs are introduced from time to time by reactor manufacturers and fuel vendors, an effort has been made to update these reactor models periodically and to expand the data bases used by the ORIGEN2 computer code. In addition, more sophisticated computational techniques than previously available were used to calculate the resulting reactor model cross-section libraries. The PWR models were based on a Westinghouse design, while the BWR models were based on a General Electric BWR/6 design. The specific reactor types considered in this report are as follows (see Glossary for the definition of these and other terms): (1) PWR-US, (2) PWR-UE, (3) BWR-US, (4) BWR-USO, and (5) BWR-UE. Each reactor model includes a unique data library that may be used to simulate the buildup and deletion of isotopes in nuclear materials using the ORIGEN2 computer code. 33 refs., 44 tabs
This project was performed by KAERI in the frame of construction of the international cooperative basis on the nuclear energy. This was supported from MOST under the title of 'Establishment of Technical Collaboration basis between Korea and France for the development of severe accident assessment computer code under high burn up condition'. The current operating NPP are converting the burned fuel to the wasted fuel after burn up of 40 GWD/MTU. But in Korea, burn up of more than 60 GWD/MTU will be expected because of the high fuel efficiency but also cost saving for storing the wasted fuel safely. The domestic research for the purpose of developing the fuel and the cladding that can be used under the high burn up condition up to 100 GWD/MTU is in progress now. But the current computer code adopts the model and the data that are valid only up to the 40 GWD/MTU at most. Therefore the current model could not take into account the phenomena that may cause differences in the fission product release behavior or in the core damage process due to the high burn up operation (more than 40 GWD/MTU). To evaluate the safety of the NPP with the high burn up fuel, the improvement of current severe accident code against the high burn up condition is an important research item. Also it should start without any delay. Therefore, in this study, an expert group was constructed to establish the research basis for the severe accident under high burn up conditions. From this expert group, the research items regarding the high burn up condition were selected and identified through discussion and technical seminars. Based on these selected items, the meeting between IRSN and KAERI to find out the cooperative research items on the severe accident under the high burn up condition was held in the IRSN headquater in Paris. After the meeting, KAERI and IRSN agreed to cooperate with each other on the selected items, and to co-host the international seminar, and to develop the model and to
The codes WIMSD-4 and BORGLES - part of the MTR-PC code package- have been applied to prepare the microscopic cross section library for the main elements of MNSR core for 6 neutron energy groups. The generated library was utilized from the 3D code CITATION to perform the calculation of fuel burn up and depletion including the identification of main fission products and its effects on the multiplication factor. In this regard some modifications have been introduced to the subroutine NUCY in CITATION to incorporate estimating the concentration of the related actinides and fission products. The burn up results indicated that the core life time of MNSR is being mainly estimated by Sm-149 following by Gd-157 and Cd-113. The accumulation of these actinides during 100 continuous operation days caused a reduction of ca. 2 mk for the excess reactivity. This result seems to be in good agreement with the available empirical value of 1.8 mk which relates to the whole discontinuous operation period of the reactor since its start and up to now. The calculation procedure simulates the sporadic operation with an adequate continuous operation period. This approximation is valid for the long lived actinides that mainly dictate the core life time. However, it is an overestimation for the concentration of short lived radioactive products like Xe-135. In the framework of this analysis the possibility of replacement of current MNSR fuel through low enriched fuels has been explored for two the fuel types U02-Mg and U3Si-Al. The results indicate that the first type (UO2-Mg) realize the criticality conditions with low enrichment of 20%, whereas the second type (U3Si-Al) required increasing the uranium enrichment up to 33%. For both fuel types the contribution of plutonium isotopes on the criticality has been also evaluated. Additionally, the influence of mixing burnable absorbers (Gd-113, Cd- 113) with the fresh fuels was investigated to identify their long-term control effect on the
Hernandez, J.L.; Alonso, G.; Perusquia, R.; Montes, J.L.; Hernandez, H. [ININ, 52045 Ocoyoacac, Estado de Mexico (Mexico)]. e-mail: jlhm@nuclear.inin-mx
2003-07-01
An evaluation of the capacity of the COREMASTER-Presto code, to evaluate generically the burnt of the control bars in the Laguna Verde reactors plant (CLV) is made. It was found that the code only reports burnt values of the control rods in MWD/TM, in spite of having with a second order polynomial model, for the conversion to remainder of the Boron-10 (B-10). It was observed that said model is adequate only for burnt smaller to 45,000 MWD/TM. To evaluate the burnt of the control rods it was reproduced the balance cycle of 18 months for the CLV, executing Cm-Presto during 13 consecutive cycles. First without rod burnt, taking this as the base case. Later on, cases with 1, 2 and up to 13 cycles with rod burnt were generated. When comparing results it was observed that the control rods pattern it loses reactivity lineally with the burnt one. By each 10 G Wd/T of burnt of the nucleus it is decreased the reactivity of the pattern rods {approx} 1 pcm in hot condition and of {approx} 20 pcm in cold condition. When burning three cycles those rods more burnt reached the 13,900 MWD/TM, equivalent to 36% of B-10 reduction, near value to 34% proposed by aging in the one lost study of B-10. It was observed that Cm-Presto it doesn't burn the superior node of the control rods when these are completely extracted. A one big lost of B-10, of the order of 50%, it represents only a decrease of 11% of the reactivity value of the rod. One can affirm that even when it is strongly decreased the content of B-10, the rod is continue considering as a black absorber, that is to say, thermal neutron that enters in the neutron rod that is absorbed. (Author)
Prediction of fission gas release at high burn-up
Reliable design of LWR fuel rods requires the fission gas release to be predicted as accurately as possible. Indeed that physical phenomenon governs both the fuel temperatures and the inner gas pressure. Fission gas release data have been reviewed by the NRC and it has been concluded that a fission gas release enhancement occurs at burn-up above 20 GWd/tM. To correct deficient fission gas release models which do not include burn-up dependence, the NRC developed an empirical correction method to describe burn-up enhancement effect. BELGONUCLEAIRE has developed its own fission gas release model which is utilized in licensing calculation through the COMETHE code. Fission gas release predictions at high burn-up are confronted to the experimental data as well as to the predictions of the NRC correlation. The physics of the fission gas release phenomenon is discussed
Transnucleaire's experience with burnup credit in transport operations
Facing a continued increase in fuel enrichment values, Transnucleaire has progressively implemented a burnup credit programme in order to maintain or, where possible, to improve the capacity of its transport packagings without physical modification. Many package design approvals, based on a notion of burnup credit, have been granted by the French competent authority for transport since the early eighties, and many of these approvals have been validated by foreign competent authorities. Up to now, these approvals are restricted to fuel assemblies made of enriched uranium and irradiated in pressurized water reactors (PWR). The characterization of the irradiated fuel and the reactivity of the package are evaluated by calculation, performed using qualified French codes developed by the CEA (Commisariat a l'Energie Atomique/French Atomic Energy Commission): CESAR as a depletion code and APOLO-MORET as a criticality code. The approvals are based on the hypothesis that the burnup considered is that applied on the least irradiated region of the fuel assemblies, the conservative approach being not to take credit for any axial profile of burnup along the fuel assembly. The most reactive configuration is calculated and the burnup credit is also restricted to major actinides only. On the operational side and in compliance with regulatory requirements, verification is made before transport, in order to meet safety objectives as required by the transport regulations. Besides a review of documentation related to the irradiation history of each fuel assembly, it consists of either a qualitative (go/no-go) verification or of a quantitative measurement, depending on the level of burnup credit. Thus the use of burnup credit is now a common practice with Transnucleaire's packages, particularly in France and Germany. New improvements are still in progress and qualifications of the calculation code are now well advanced, which will allow in the near future the use of six selected
Burnup dependent core neutronic analysis for PBMR
The strategy for core neutronics modeling is based on SCALE4.4 code KENOV.a module that uses Monte Carlo calculational methods. The calculations are based on detailed unit cell and detailed core modeling. The fuel pebble is thoroughly modeled by introducing unit cell modeling for the graphite matrix and the fuel kernels in the pebble. The core is then modeled by placing these pebbles randomly throughout the core, yet not loosing track of any one of them. For the burnup model, a cyclic manner is adopted by coupling the KENOV.a and ORIGEN-S modules. Shifting down one slice at each discrete time step, and inserting fresh fuel from the top, this cyclic calculation model continues until equilibrium burnup cycle is achieved. (author)
Burn-up credit analysis of RBMK-1000 an WWER-1000 spent nuclear fuel accounting only for actinides is carried out and a method is proposed for actinide burn-up credit. Two burn-up credit approaches are analyzed, which consider a system without and with the distribution of isotopes along the height of the fuel assembly. Calculations are performed using SCALE and MCNP computer codes
Fuel burnup analysis for the Moroccan TRIGA research reactor
Highlights: ► A fuel burnup analysis of the 2 MW TRIGA MARK II Moroccan research reactor was established. ► Burnup calculations were done by means of the in-house developed burnup code BUCAL1. ► BUCAL1 uses the MCNP tallies directly in the calculation of the isotopic inventories. ► The reactor life time was found to be 3360 MW h considering full power operating conditions. ► Power factors and fluxes of the in-core irradiation positions are strongly affected by burnup. -- Abstract: The fundamental advantage and main reason to use Monte Carlo methods for burnup calculations is the possibility to generate extremely accurate burnup dependent one group cross-sections and neutron fluxes for arbitrary core and fuel geometries. Yet, a set of values determined for a material at a given position and time remains accurate only in a local region, in which neutron spectrum and flux vary weakly — and only for a limited period of time, during which changes of the local isotopic composition are minor. This paper presents the approach of fuel burnup evaluation used at the Moroccan TRIGA MARK II research reactor. The approach is essentially based upon the utilization of BUCAL1, an in-house developed burnup code. BUCAL1 is a FORTRAN computer code designed to aid in analysis, prediction, and optimization of fuel burnup performance in nuclear reactors. The code was developed to incorporate the neutron absorption reaction tally information generated directly by MCNP5 code in the calculation of fissioned or neutron-transmuted isotopes for multi-fueled regions. The fuel cycle length and changes in several core parameters such as: core excess reactivity, control rods position, fluxes at the irradiation positions, axial and radial power factors and other parameters are estimated. Besides, this study gives valuable insight into the behavior of the reactor and will ensure better utilization and operation of the reactor during its life-time and it will allow the establishment of
Effect of Core Configurations on Burn-Up Calculations For MTR Type Reactors
Three-dimensional burn-up calculations of MTR-type research reactor were performed using different patterns of control rods , to examine their effect on power density and neutron flux distributions throughout the entire core and on the local burn-up distribution. Calculations were performed using the computer codes' package MTRPC system, using the cell calculation transport code WIMS-D4 and the core calculation diffusion code CITVAP. A depletion study was done and the effects on the reactor fuel were studied, then an empirical formula was generated for every fuel element type, to correlate irradiation to burn-up percentage. Keywords: Neutronic Calculations, Burn-Up, MTR-Type Research Reactors, MTRPC Package, Empirical Formula For Fuel Burn-Up.
Taking burnup credit for interim storage and transportation system for BWR fuels
In the back-end issues of nuclear fuel cycle, selection of reprocessing or one-through is a big issue. For both of the cases, a reasonable interim storage and transportation system is required. This study proposes an advanced practical monitoring and evaluation system. The system features the followings: (l) Storage racks and transportation casks taking credit for burnup. (2) A burnup estimation system using a compact monitor with Cd- Te detectors and fission chambers. (3) A neutron emission-rate evaluation methodology, especially important for high burnup MOX fuels. (4) A nuclear materials management system for safeguards. Current storage system and transport casks are designed on the basis of a fresh fuel assumption. The assumption is too conservative. Taking burnup credit gives a reasonable design while keeping conservatism. In order to establish a reasonable burnup credit design system, a calculation system has been developed for determining isotope compositions, burnup, and criticality. The calculation system consists of some modules such as TGBLA, ORIGEN, CITATION, MCNP and KENO. The TGBLA code is a fuel design code for LWR fuels developed in TOSHIBA Corporation. The code takes operational history such as, power density, void fraction into account. This code is applied to the back-end issues for a more accurate design of a storage and a transportation system. The ORIGEN code is well-known one-point isotope depletion code. In the calculation system, the code calculates isotope compositions using libraries generated from the TGBLA code. The CITATION code, the MCNP code, and the KENO code are three dimensional diffusion code, continuous energy Monte Carlo code, discrete energy Monte Carlo code, respectively. Those codes calculate k- effective of the storage and transportation systems using isotope compositions generated from the ORIGEN code. The CITATION code and the KENO code are usually used for practical designs. The MCNP code is used for reference
Methods used in burn-up determination of the irradiated fuel rods at TRIGA reactor
A short presentation of the methods used at INR TRIGA reactor for the burn-up determination is given together with some considerations on ORIGEN 2 computer code used for calculating fission products activities and nuclide concentration. Burn-up is determined by gamma spectroscopy and thermal power monitoring. (Author)
Application of Candle burnup to small fast reactor
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 to an axial direction. An equilibrium state was obtained for a large fast reactor (core radius is 2 m and reflector thickness is 0.5 m) successfully by using a newly developed direct analysis code. However, it is difficult to apply this burnup strategy to small reactors, since its neutron leakage becomes large and neutron economy becomes worse. Fuel enrichment should be increased in order to sustain the criticality. However, higher enrichment of fresh fuel makes the CANDLE burnup difficult. We try to find some small reactor designs, which can realize the CANDLE burnup. We have successfully find a design, which is not the CANDLE burnup in the strict meaning, but satisfies qualitatively its characteristics mentioned at the top of this abstract. In the final paper, the general description of CANDLE burnup and some results on the obtained small fast reactor design are presented.(author)
Development of a Burnup Program based on the Krylov Subspace Method
The depletion calculation of the DeCART code has been performed by the support of the ORIGEN code. Recently, a burnup program based on the Krylov subspace method is developed and implemented to the DeCART code. Numerical solution for the burnup equation by the Krylov subspace method is well described. Therefore, this paper describes the Krylov subspace method for a burnup equation briefly in Section 2, and focuses on the DeCART solution for a pin cell problem by comparing it with the HELIOS solution
Development of a Burnup Program based on the Krylov Subspace Method
Cho, Jin-Young; Shim, Hyung-Jin; Kim, Kang-Seog; Song, Jae-Seung; Lee, Chung-Chan [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)
2007-10-15
The depletion calculation of the DeCART code has been performed by the support of the ORIGEN code. Recently, a burnup program based on the Krylov subspace method is developed and implemented to the DeCART code. Numerical solution for the burnup equation by the Krylov subspace method is well described. Therefore, this paper describes the Krylov subspace method for a burnup equation briefly in Section 2, and focuses on the DeCART solution for a pin cell problem by comparing it with the HELIOS solution.
A burnup credit calculation methodology for PWR spent fuel transportation
A burnup credit calculation methodology for PWR spent fuel transportation has been developed and validated in CEA/Saclay. To perform the calculation, the spent fuel composition are first determined by the PEPIN-2 depletion analysis. Secondly the most important actinides and fission product poisons are automatically selected in PEPIN-2 according to the reactivity worth and the burnup for critically consideration. Then the 3D Monte Carlo critically code TRIMARAN-2 is used to examine the subcriticality. All the resonance self-shielded cross sections used in this calculation system are prepared with the APOLLO-2 lattice cell code. The burnup credit calculation methodology and related PWR spent fuel transportation benchmark results are reported and discussed. (authors)
Burnup instabilities in the full-core HTR model simulation
Highlights: • We performed full-core burnup calculation coupled with Monte Carlo code. • Depletion instabilities have been detected for HTR system at high burnup. • We assess the stability of time step models in application to core calculation. • Discussion of the modeling factors related to burnup core simulation is presented. - Abstract: The phenomenon of numerical instabilities present in the Monte Carlo burnup calculations has been shown and explained by many authors using models of LWR, often simplified. Some theoretical considerations about origins of oscillations are very general, however it may be difficult to apply it easily to other models as a prediction of stability. Physics of HTR core differs significantly from the properties of light water system and the reliable extrapolation of the current numerical results is not possible. Moreover, most of the works concerning HTR burnup calculations put no emphasis on the spatial stability of the simulation and apply very long time steps. The awareness in this field of research seems to be not sufficient. In this paper, we focus on the demonstration of depletion instabilities in the simulations of HTR core dedicated for deep burnup of plutonium and minor actinides. We apply various methodology of time step implemented in advanced Continuous Energy Monte Carlo burnup code MCB version 5. Stability analysis is very rare for the full core calculations and the awareness of the oscillation’s problem is obligatory for the reliable modeling of a fuel cycle. In the summary of this work we systematize and discuss factors related to the stability of depletion and review available solutions
Increased burnup of fuel elements
The specialists' group for fuel elements of the Kerntechnische Gesellschaft e.V. held a meeting on ''Increased Burnup of Fuel Elements'' on 9th and 10th of November 1982 at the GKSS Research Center Geesthacht. Most papers dealt with the problems of burnup increase of fuel elements for light water reactors with respect to fuel manufacturing, power plant operation and reprocessing. Review papers were given on the burnup limits for high temperature gas cooled reactors and sodium fast breeder reactors. The meeting ended with a presentation of the technical equipment of the hot laboratory of the GKSS and the programs which are in progress there. (orig.)
Dependence of control rod worth on fuel burnup
Savva, P., E-mail: savvapan@ipta.demokritos.g [NCSR ' DEMOKRITOS' , PoB 60228, 15310 Aghia Paraskevi (Greece); Varvayanni, M., E-mail: melina@ipta.demokritos.g [NCSR ' DEMOKRITOS' , PoB 60228, 15310 Aghia Paraskevi (Greece); Catsaros, N., E-mail: nicos@ipta.demokritos.g [NCSR ' DEMOKRITOS' , PoB 60228, 15310 Aghia Paraskevi (Greece)
2011-02-15
Research highlights: Diffusion and MC calculations for rod worth dependence on burnup and Xe in reactors. One-step rod withdrawal/insertion are used for rod worth estimation. The study showed that when Xe is present the rods worth is significantly reduced. Rod worth variation with burnup depends on rod position in core. Rod worth obtained with MC code is higher than that obtained from deterministic. - Abstract: One important parameter in the design and the analysis of a nuclear reactor core is the reactivity worth of the control rods, i.e. their efficiency to absorb excess reactivity. The control rod worth is affected by parameters such as the fuel burnup in the rod vicinity, the Xe concentration in the core, the operational time of the rod and its position in the core. In the present work, two different computational approaches, a deterministic and a stochastic one, were used for the determination of the rods worth dependence on the fuel burnup level and the Xe concentration level in a conceptual, symmetric reactor core, based on the MTR fuel assemblies used in the Greek Research Reactor (GRR-1). For the deterministic approach the neutronics code system composed by the SCALE modules NITAWL and XSDRN and the diffusion code CITATION was used, while for the stochastic one the Monte Carlo code TRIPOLI was applied. The study showed that when Xe is present in the core, the rods worth is significantly reduced, while the rod worth variation with increasing burnup depends on the rods position in the core grid. The rod worth obtained with the use of the Monte Carlo code is higher than the one obtained from the deterministic code.
Analysis of burnup credit on spent fuel transport / storage casks - estimation of reactivity bias
Chemical analyses of high burnup UO2 (65 GWd/t) and MOX (45 GWd/t) spent fuel pins were carried out. Measured data of nuclides' composition from U234 to P 242 were used for evaluation of ORIGEN-2/82 code and a nuclear fuel design code (NULIF). Critically calculations were executed for transport and storage casks for 52 BWR or 21 PWR spent fuel assemblies. The reactivity biases were evaluated for axial and horizontal profiles of burnup, and historical void fraction (BWR), operational histories such as control rod insertion history, BPR insertion history and others, and calculational accuracy of ORIGEN-2/82 on nuclides' composition. This study shows that introduction of burnup credit has a large merit in criticality safety analysis of casks, even if these reactivity biases are considered. The concept of equivalent uniform burnup was adapted for the present reactivity bias evaluation and showed the possibility of simplifying the reactivity bias evaluation in burnup credit. (authors)
High burnup experience in PWRs
The purpose of this paper is to summarize the high burnup experience of Westinghouse PWR fuel. The emphasis is on two regions of commercial PWR fuel that attained region average burnups greater than 36,000 MWD/MTU. One region operated under load follow conditions. The other region operated at base load conditions with a high average linear heat rating. Coolant activity data and post irradiation data were obtained. The post-irradiation data consisted of visual examinations, crud sampling, rod-to-rod dimensional changes, fuel column length changes, rod and assembly growth, assembly bow, fuel rod profilometry, grid spring relaxation, and fuel assembly sipping tests. The data showed that the fuel operated reliably to this burnup. Plans for irradiation to higher burnups are also discussed
One of the major problems in burnup studies is the reasonably fast and accurate calculation of the space-and-energy dependent neutron flux and reaction rates for realistic power reactor fuel geometries and compositions, and its optimal integration in the global reactor calculations. The scope of the present research was to develop improved methods trying to satisfy the above requirements. In the epithermal region, simple and efficient approximation is proposed which allows the analytical solution for the space dependence of the spherical harmonics flux moments, and hence the derivation of the recurrence relations between he flux moments at successive lethargy pivotal points. A new matrix formalism to invert the coefficient matrix of band structure resulted in a reduce computer time and memory demands. The research on epithermal region is finalized in computing programme SPLET, which calculates the space-lethargy distribution of the spherical harmonics neutron flux moments, and the related integral quantities as reaction rates and resonance integrals. For partial verification of the above methods a Monte Carlo procedure was developed. Using point-wise representation of variables, a flexible and fast convergent integral transport method SEPT i developed. Expanding the neutron source and flux in finite series of arbitrary polynomials, the space-and-energy dependent integral transport equation is transformed into a general linear algebraic form, which is solved numerically. A simple and efficient procedure for deriving multipoint equations and constructing matrix is proposed and examined, and no unwanted oscillations were noticed. The energy point method was combined with the spherical harmonics method as well. A multi zone few-group program SPECTAR for global reactor calculations was developed. For testing, the flux distribution, neutron leakage and effective multiplication factor for the PWR reactor of the power station San Onofre were calculated. In order to verify
External linking scripts between Monte Carlo transport codes and burnup codes, and complete integration of burnup capability into Monte Carlo transport codes, have been or are currently being developed. Monte Carlo linked burnup methodologies may serve as an excellent benchmark for new deterministic burnup codes used for advanced systems; however, there are some instances where deterministic methodologies break down (i.e., heavily angularly biased systems containing exotic materials without proper group structure) and Monte Carlo burn up may serve as an actual design tool. Therefore, researchers are also developing these capabilities in order to examine complex, three-dimensional exotic material systems that do not contain benchmark data. Providing a reference scheme implies being able to associate statistical errors to any neutronic value of interest like k(eff), reaction rates, fluxes, etc. Usually in Monte Carlo, standard deviations are associated with a particular value by performing different independent and identical simulations (also referred to as 'cycles', 'batches', or 'replicas'), but this is only valid if the calculation itself is not biased. And, as will be shown in this paper, there is a bias in the methodology that consists of coupling transport and depletion codes because Bateman equations are not linear functions of the fluxes or of the reaction rates (those quantities being always measured with an uncertainty). Therefore, we have to quantify and correct this bias. This will be achieved by deriving an unbiased minimum variance estimator of a matrix exponential function of a normal mean. The result is then used to propose a reference scheme to solve Boltzmann/Bateman coupled equations, thanks to Monte Carlo transport codes. Numerical tests will be performed with an ad hoc Monte Carlo code on a very simple depletion case and will be compared to the theoretical results obtained with the reference scheme. Finally, the statistical error propagation
CB2 result evaluation (VVER-440 burnup credit benchmark)
The second portion of the four-piece international calculational benchmark on the VVER burnup credit (CB2) prepared in the collaboration with the OECD/NEA/NSC Burnup Credit Criticality Benchmarks Working Group and proposed to the AER research community has been evaluated. The evaluated results of calculations performed by analysts from Cuba, the Czech Republic, Finland, Germany, Russia, Slovakia and the United Kingdom are presented. The goal of this study is to compare isotopic concentrations calculated by the participants using various codes and libraries for depletion of the VVER-440 fuel pin cell. No measured values were available for the comparison. (author)
Power excursion analysis for BWR`s at high burnup
Diamond, D.J.; Neymoith, L.; Kohut, P. [Brookhaven National Lab., Upton, NY (United States)
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.
2008-01-01
Citigroup,one of the World top 500 companies,has now settled in Excel Center,Financial Street. The opening ceremony of Excel Center and the entry ceremony of Citigroup in the center were held on March 31.Government leaders of Xicheng District,the Excel CEO and the heads of Asia-Pacific Region leaders of Citibank all participated in the ceremony.
Belo, Thiago F.; Fiel, Joao Claudio B., E-mail: thiagofbelo@hotmail.com [Instituto Militar de Engenharia (IME), Rio de Janeiro, RJ (Brazil)
2015-07-01
Nuclear reactor core analysis involves neutronic modeling and the calculations require problem dependent nuclear data generated with few neutron energy groups, as for instance the neutron cross sections. The methods used to obtain these problem-dependent cross sections, in the reactor calculations, generally uses nuclear computer codes that require a large processing time and computational memory, making the process computationally very expensive. Presently, analysis of the macroscopic cross section, as a function of nuclear parameters, has shown a very distinct behavior that cannot be represented by simply using linear interpolation. Indeed, a polynomial representation is more adequate for the data parameterization. To provide the cross sections of rapidly and without the dependence of complex systems calculations, this work developed a set of parameterized cross sections, based on the Tchebychev polynomials, by fitting the cross sections as a function of nuclear parameters, which include fuel temperature, moderator temperature and density, soluble boron concentration, uranium enrichment, and the burn-up. In this study is evaluated the problem-dependent about fission, scattering, total, nu-fission, capture, transport and absorption cross sections for a typical PWR fuel element reactor, considering burn-up cycle. The analysis was carried out with the SCALE 6.1 code package. The results of comparison with direct calculations with the SCALE code system and also the test using project parameters, such as the temperature coefficient of reactivity and fast fission factor, show excellent agreements. The differences between the cross-section parameterization methodology and the direct calculations based on the SCALE code system are less than 0.03 percent. (author)
The commercial and technological impact of high burnup
Deregulation of electricity markets is driving prices downward. Consequently utilities continue to demand the minimization of electrical production costs. Fuel cycle cost savings are valued as a strong contributor, although directly representing only about one third of electricity generating costs. Burnups consistent with the current enrichment limit of 5 w/0 will be required. Significant progress has already been achieved by Siemens in meeting the demands of utilities for increased fuel burnup. The technological challenges imposed are mainly related to corrosion and hydrogen pickup of the clad, the properties of the fuel and the dimensional changes of the structure. Clad materials with increased corrosion resistance have been developed. The high burnup behaviour of the fuel has been extensively investigated and the decrease of thermal conductivity, the rim effect and the increase of fission gas release can be described, with good accuracy, in fuel rod computer codes. Advanced statistical design methods have been developed and introduced. In summary, most of the questions about the fuel operational behaviour and reliability in the high burnup range have been solved or the solutions are visible. The main licensing challenges for high burnup fuel are currently seen for accident condition analyses, especially for RIA and LOCA. (author)
IFPE/US-PWR-16 X 16 Lead Test Assembly Extended Burnup Demonstration Program
power and thermal hydraulic data as well the fuel design, geometry, and composition data is contained in tables and figures within the main body of the report. The axial power shapes, fuel rod power and burnup histories, and PIE data are attached to the report as appendices. Additionally, Excel spreadsheets containing the axial power shapes and fuel rod power and burnup histories are included in the transmittal to ease use in modeling codes. Section and Description: Appendix 1 Axial power shapes; Appendix 2 Fuel rod power and burnup histories; Appendix 3 Internal void volume and fission gas release measurements; Appendix 4 Diametral strain measurements; Appendix 5 Fuel rod gamma activity measurements; Appendix 6 Fuel burnup analyses; Appendix 7 Fuel density measurements; Appendix 8 Fuel grain and pore size distributions; Appendix 9 Cladding oxide thickness measurements; Appendix 10 Hydrogen concentration and absorption measurements; Appendix 11 Cladding mechanical property testing
PWR AXIAL BURNUP PROFILE ANALYSIS
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 andW 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)
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).
Results of the isotopic concentrations of WWER calculation Burnup Credit Benchmark NO.2 (CB2)
The purpose of this document is to present the results of the nuclide concentrations of the WWER Burnup Credit Benchmark NO.2 (CB2) that were performed in The Nuclear Technology Center of Cuba with available codes and libraries. The CB2 benchmark specification as the second phase of the WWER burnup credit benchmark is summarized in [1]. The CB2 benchmark focused on WWER burnup credit study proposed on the 97' Atomic Energy Research symposium [2]. The obtained results are isotopic concentrations of spent fuel as a function of the burnup and cooling time. The depletion point 'ORIGEN2'[3] code was used for the calculation of the spent fuel concentration. This work also comprises the results obtained by other codes [4]. (Author)
The impact of burn-up credit in criticality studies
Nowadays optimization goes with everything. So French engineering firms try to demonstrate that fuel transport casks and storage pools are able to receive assemblies with higher 235U initial enrichments. Fuel Burnup distribution contributes to demonstrate it. This instruction has to elaborate a way to take credit of burnup effects on criticality safety designs. The calculation codes used are CESAR 4.21-APOLLO 1-MORET III. The assembly studied (UO2) is irradiated in a French Pressurized Water Reactor like EDF nuclear power reactor: PWR 1300 MWe, 17 x 17 array. Its initial enrichment in 235U equals 4.5%. The studies exposed in this report have evaluated the effects of: i) the 15 fission products considered in Burnup Credit (95Mo, 99Tc, 101Ru, 103Rh, 109Ag, 133Cs, 143Nd, 145Nd, 147Sm, 149Sm, 150Sm, 151Sm, 152Sm, 153Eu, 155Gd), ii) the calculated abundances corrected or not by fixed factors, iii) the choice of one cross sections library used by CESAR 4.21, iu) the zone number elected in the axial burnup distribution zoning, u) the kind of cut applied on (regular/optimized). Two axial distribution profiles are studied: one with 44 GWd/t average burnup, the other with 20 GWd/t average burnup. The second one considers a shallow control rods insertion in the upper limit of the assembly. The results show a margin in reactivity about 0.045 with consideration of the 6 most absorbent fission products (103Rh, 133Cs, 143Nd, 149Sm, 152Sm, 155Gd), and about 0.06 for all Burnup Credit fission products whole. Those results have been calculated with an average burnup of 44 GWj/t. In a conservative approach, corrective factors must be apply on the abundance of some fission products. The cross sections library used by CESAR 4.21 (BBL 4) is sufficient and gives satisfactory results. The zoning of the assembly axial distribution burnup in 9 regular zones grants a satisfying calculation time/result precision compromise. (author)
Verschuuren, Gerard M
2013-01-01
Covering a variety of Excel simulations, from gambling to genetics, this introduction is for people interested in modeling future events, without the cost of an expensive textbook. The simulations covered offer a fun alternative to the usual Excel topics and include situations such as roulette, password cracking, sex determination, population growth, and traffic patterns, among many others.
Tyrrell, Sidney
2006-01-01
In this brief article, the author instructs teachers on how to produce an interactive spreadsheet from scratch in about 20 minutes and en route equip themselves and their students, with handy "Excel" skills. The aim is to introduce the basics of "Excel," plus some fun bits, speedily and with a purpose; producing something that is useful in its own…
CANDU lattice uncertainties during burnup
Uncertainties associated with fundamental nuclear data accompany evaluated nuclear data libraries in the form of covariance matrices. As nuclear data are important parameters in reactor physics calculations, any associated uncertainty causes a loss of confidence in the calculation results. The quantification of output uncertainties is necessary to adequately establish safety margins of nuclear facilities. In this work, microscopic cross-section has been propagated through lattice burnup calculations applied to a generic CANDU® model. It was found that substantial uncertainty emerges during burnup even when fission yield fraction and decay rate uncertainties are neglected. (author)
Fast reactor 3D core and burnup analysis using VESTA
Luciano, N.; Shamblin, J.; Maldonado, I. [Nuclear Engineering Dept., Univ. of Tennessee, Knoxville, TN 37996-2300 (United States)
2012-07-01
Burnup analyses using the VESTA code have been performed on a MOX-fuelled fast reactor model as specified by an IAEA computational benchmark. VESTA is a relatively new code that has been used for burnup credit calculations and thermal reactor models, but not typically for fast reactor applications. The detailed input and results of the IAEA benchmark provides an opportunity to gauge the use of VESTA in a fast reactor application. VESTA employs an ultra-fine multi-group binning approach that accelerates Monte Carlo burnup calculations. Using VESTA to compute the end of cycle (EOC) power fractions by enrichment zone showed agreement with the published values within 5%. When comparing the ultra-fine multi-group binning approach to the tally-based approach, EOC isotopic masses also agree within 5%. Using the ultra-fine multi-group binning approach, we obtain a wall-time speedup factor of 35 when compared to the tally-based approach for computing a k{sub eff} eigenvalue with burnup problem. The authors conclude the use of VESTA's ultra-fine multi-group binning approach with Monte Carlo transport performs accurate depletion calculations for this fast reactor benchmark. (authors)
Extension of the TRANSURANUS burn-up model
The validation range of the model in the TRANSURANUS fuel performance code for calculating the radial power density and burn-up in UO2 fuel has been extended from 64 MWd/kgHM up to 102 MWd/kgHM, thereby improving also its precision. In addition, the first verification of calculations with post-irradiation examination data is reported for LWR-MOX fuel with a rod average burn-up up to 45 MWd/kgHM. The extension covers the inclusion of new isotopes in order to account for the production of 238Pu. The corresponding one-group cross-sections used in the equations rely on results obtained with ALEPH, a new Monte Carlo burn-up code. The experimental verification is based on electron probe microanalysis (EPMA) and on secondary ion mass spectrometry (SIMS) as well as radiochemical data of fuel irradiated in commercial power plants. The deviations are quantified in terms of frequency distributions of the relative errors. The relative errors on the burn-up distributions in both fuel types remain below 12%, corresponding to the experimental scatter
The criticality safety analysis of spent fuel systems has traditionally assumed that the fuel is fresh. This results in significant conservatism in the calculated value of the system's reactivity. Improved calculational methods allows one to take credit for the reactivity reduction associated with fuel burnup, hence reducing the analysis conservatism while maintaining an adequate criticality safety margin. Motivation for using burnup credit in criticality safety applications is generally based on economic considerations. Although economics may be a primary factor in deciding to use burnup credit, other benefits may be realized. Many of the additional benefits of burnup credit that are not strictly economic, may be considered to contribute to public health and safety, and resource conservation and environmental quality. Interest in the implementation of burnup credit has been shown by many countries. A summary of the information gathered by the IAEA about ongoing activities and regulatory status of burnup credit in different countries is included. Burnup credit implementation introduces new parameters and effects that should be addressed in the criticality analysis (e.g., axial and radial burnup shapes, fuel irradiation history, and others). Analysis of these parameters introduces new variations as well as the uncertainties, that should be considered in the safety assessment of the system. Also, the need arises to validate the isotopic composition that results from a depletion calculation, as well as to extend the current validation range of criticality codes to cover spent fuel. The use of burnup credit implies a verification of the fuel burnup before loading for transport, storage, disposal, or reprocessing each assembly, to make sure that the burnup level achieved complies with the criteria established. Methods and procedures used in different countries are described in this report
Burnup determination of water reactor fuel
The present meeting was scheduled by the International Atomic Energy Agency in consultation with the Members of the International Working Group on Water Reactor Fuel Performance and Technology. The meeting was hosted by the Commission of the European Communities, at the Transuranium Research Laboratory, Joint Research Centre Karlsruhe, in the Federal Republic of Germany. This subject was dealt with for the first time by the IAEA. It was found to correspond adequately to this type of Specialist Meeting and to be suitable at a moment when the extension of burnup constitutes a major technical and economical issue in fuel technology. It was stressed that analysis of highly burnt fuels, mixed oxides and burnable absorber bearing fuels required extension of the experimental data base, to comply with the increasing demand for an improved fuel management, including better qualification of reactor physics codes. Twenty-seven participants from eleven countries plus two international organizations attended the Meeting. Twelve papers were given during three technical sessions, followed by a panel discussion which allowed to formulate the conclusions of the meeting and recommendations to the Agency. In addition, participants were invited to give an outline of their national programmes, related to Burnup Determination of Water Reactor Fuel. A separate abstract was prepared for each of these 12 papers. Refs, figs and tabs
Power excursion analysis for high burnup cores
A study was undertaken of power excursions in high burnup cores. There were three objectives in this study. One was to identify boiling water reactor (BWR) and pressurized water reactor (PWR) transients in which there is significant energy deposition in the fuel. Another was to analyze the response of BWRs to the rod drop accident (RDA) and other transients in which there is a power excursion. The last objective was to investigate the sources of uncertainty in the RDA analysis. In a boiling water reactor, the events identified as having significant energy deposition in the fuel were a rod drop accident, a recirculation flow control failure, and the overpressure events; in a pressurized water reactor, they were a rod ejection accident and boron dilution events. The RDA analysis was done with RAMONA-4B, a computer code that models the space- dependent neutron kinetics throughout the core along with the thermal hydraulics in the core, vessel, and steamline. The results showed that the calculated maximum fuel enthalpy in high burnup fuel will be affected by core design, initial conditions, and modeling assumptions. The important uncertainties in each of these categories are discussed in the report
RAPID program to predict radial power and burnup distribution of UO{sub 2} fuel
Lee, Chan Bock; Song, Jae Sung; Bang, Je Gun; Kim, Dae Ho [Korea Atomic Energy Research Institute, Taejon (Korea)
1999-02-01
Due to the radial variation of the neutron flux and its energy spectrum inside UO{sub 2} fuel, the fission density and fissile isotope production rates are varied radially in the pellet, and it becomes necessary to know the accurate radial power and burnup variation to predict the high burnup fuel behavior such as rim effects. Therefore, to predict the radial distribution of power, burnup and fissionable nuclide densities in the pellet with the burnup and U-235 enrichment, RAPID(RAdial power and burnup Prediction by following fissile Isotope Distribution in the pellet) program was developed. It considers the specific radial variation of the neutron reaction of the nuclides while the constant radial variation of neutron reaction except neutron absorption of U-238 regardless of the nuclides, the burnup and U-235 enrichment is assumed in TUBRNP model which is recognized as the one of the most reliable models. Therefore, it is expected that RAPID may be more accurate than TUBRNP, specially at high burnup region. RAPID is based upon and validated by the detailed reactor physics code, HELIOS which is one of few codes that can calculates the radial variations of the nuclides inside the pellet. Comparison of RAPID prediction with the measured data of the irradiated fuels showed very good agreement. RAPID can be used to calculate the local variations of the fissionable nuclide concentrations as well as the local power and burnup inside that pellet as a function of the burnup up to 10 w/o U-235 enrichment and 150 MWD/kgU burnup under the LWR environment. (author). 8 refs., 50 figs., 1 tab.
Determination of axial profit performed burnup credit by SCALE 4.3-system
SCALE 4.3 is a modular code system designed for realizing standard computational analysis for licensing evaluation. Since now, spent fuel storage pools criticality analysis have been done considering this fuel as fresh, with its maximum enrichment. With burnup credit we can obtain cheaper and compact configurations. The procedure for calculating a spent fuel storage consists of a burnup calculation plus a criticality calculation. We can perform a conservative approximation for the burnup calculations using 1-D results, but, besides the geometry configurations for the 3-D criticality calculation. we need an appropriate approximation to model the burnup axial variation. We assume that for a burnup profile set, the most conservative profile is between the lower and the upper range of this profile, set. We consider only combinations of the maximum and minimum burnup in each axial region, for each burnup range. This gives an estimation of the different burnup shapes effect and the general characteristics of the most conservative shape. (Author) 6 refs
Adeyemi, Oluseyi
2011-01-01
Sourcing Excellence is one of the key performance indicators (KPIs) in this world of ever changing sourcing strategies. Manufacturing companies need to access and diagnose the reliability and competencies of existing suppliers in order to coordinate and develop them. This would help in managing...... excellence in manufacturing companies. In addition, an advanced approach to sourcing involves integrating and coordinating suppliers, common materials, processes, designs and technologies across worldwide buying, design and operating locations. However, since most manufacturing companies do not have well...... integrating the models of JIT sourcing and global sourcing would also add value and enhance excellent sourcing practices. These and other possibilities of sourcing would be considered in this research....
Burnup credit calculations on long-term disposal
One of the considered options for handling of irradiated nuclear fuel is the final disposal in some kind of repository. This necessitates the long-term investigation of subcriticality, heat production, public dose etc. NEA WPNCS Burnup Credit Expert Group defined a new benchmark to test the codes and data used for such problems. The effect of cooling time should be investigated. This implies that the decay data and not the cross sections influence the results. Composition of 4.5 % UO2 fuel with 50 MWd/kgU is given at the assembly removal from the core. Change of composition should be evaluated for 30 values of cooling time up to 1 million years. Keff should be evaluated with these compositions for a container housing 21 fuel assemblies. Initial concentration of 115 isotopes is given. For criticality calculations the usual 'burnup credit set' is used (14 actinides and 15 fission products). Results for additional isotopes is not presented now. The investigated fuel is 17 x 17 PWR UO2 type, with 25 guide tubes. The selected cooling times covers the time intervals of the usual handling procedures around the reactors (few years storing in storage pool, transport), interim storage (hundred years), and the long time scale of disposal up to 1 million years. Results: 1) For major actinides, ORIGEN and MULTICELL based keff results are practically identical up to 1000 years, far beyond the cooling times it was intended. 2) For actinides and fission products, the agreement is excellent up to 100 years, which covers the interim storage. 3) The difference of keff results about 0.02 at 1000 years. The reason is mainly the presence of Np-237, not considered in the previous case. It is produced from Am-241 by α-decay (432 years). Compositions calculated by ORIGEN and TIBSO results the same keff values for cooling times up to 1 million years. Changes in keff with cooling time have clear physical explanation. Compositions calculated by ORIGEN and MULTICELL results the same keff
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
TRIGA fuel burn-up calculations and its confirmation
The Cesium (Cs-137) isotopic concentration due to irradiation of TRIGA Fuel Elements FE(s) is calculated and measured at the Atominstitute (ATI) of Vienna University of Technology (VUT). The Cs-137 isotope, as proved burn-up indicator, was applied to determine the burn-up of the TRIGA Mark II research reactor FE. This article presents the calculations and measurements of the Cs-137 isotope and its relevant burn-up of six selected Spent Fuel Elements SPE(s). High-resolution gamma-ray spectroscopy based non-destructive method is employed to measure spent fuel parameters. By the employment of this method, the axial distribution of Cesium-137 for six SPE(s) is measured, resulting in the axial burn-up profiles. Knowing the exact irradiation history and material isotopic inventory of an irradiated FE, six SPE(s) are selected for on-site gamma scanning using a special shielded scanning device developed at the ATI. This unique fuel inspection unit allows to scan each millimeter of the FE. For this purpose, each selected FE was transferred to the fuel inspection unit using the standard fuel transfer cask. Each FE was scanned at a scale of 1 cm of its active length and the Cs-137 activity was determined as proved burn-up indicator. The measuring system consists of a high-purity germanium detector (HPGe) together with suitable fast electronics and on-line PC data acquisition module. The absolute activity of each centimeter of the FE was measured and compared with reactor physics calculations. The ORIGEN2, a one-group depletion and radioactive decay computer code, was applied to calculate the activity of the Cs-137 and the burn-up of selected SPE. The deviation between calculations and measurements was in range from 0.82% to 12.64%.
Development of the SUHAM-U code for burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel is described. Developed SUHAM-U code has capacity to calculate burnup in each fuel or poison zone of each cell of VVER-1000 fuel assembly. In so doing Surface Harmonics Method is used for calculation of the detail neutron spectra in fuel assembly at separated burnup values. Verification of SUHAM-U code by burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel has been carried out. Comparisons were carried out with calculations by UNK and RECOL codes. UNK code uses the first collisions probabilities method for solution of the neutron transport equation and RECOL code uses Monte-Carlo method with point-wise continues energy presentation of cross-sections. The main conclusion of all comparisons is the SUHAM-U code calculates the fuel burnup of VVER-1000 fuel assemblies with uranium and MOX fuel with enough high accuracy. Time expenditures are adduced. (authors)
Boyarinov, V. F.; Davidenko, V. D.; Polismakov, A. A.; Tsibulsky, V. F. [Russian Research Center Kurchatov Inst., Nuclear Reactor Inst., 123182, Moscow (Russian Federation)
2006-07-01
Development of the SUHAM-U code for burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel is described. Developed SUHAM-U code has capacity to calculate burnup in each fuel or poison zone of each cell of VVER-1000 fuel assembly. In so doing Surface Harmonics Method is used for calculation of the detail neutron spectra in fuel assembly at separated burnup values. Verification of SUHAM-U code by burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel has been carried out. Comparisons were carried out with calculations by UNK and RECOL codes. UNK code uses the first collisions probabilities method for solution of the neutron transport equation and RECOL code uses Monte-Carlo method with point-wise continues energy presentation of cross-sections. The main conclusion of all comparisons is the SUHAM-U code calculates the fuel burnup of VVER-1000 fuel assemblies with uranium and MOX fuel with enough high accuracy. Time expenditures are adduced. (authors)
This paper describes the experience gained in Germany in implementing burnup credit in wet storage and dry transport systems of spent PWR, BWR, and MOX fuel. It gives a survey of the levels of burnup credit presently used, the regulatory status and activities planned, the fuel depletion codes and criticality calculation codes employed, the verification methods used for validating these codes, the modeling assumptions made to ensure that the burnup credit criticality analysis is based on a fuel irradiation history which leads to bounding neutron multiplication factors, and the implementation of procedures used for fuel loading verification. (author)
The influence of pitch, burnup and absorber rods on the spent fuel pool criticality
It has been shown that supercriticality might occur for some postulated accidents for the TRIGA spent fuel pool at ''Josef Stefan'' Institute in Ljubljana, Slovenia. However, in the previous studies, the effect of burnup was not accounted for. In this work the dependence of criticality on fuel burnup, the pitch among the elements and the number of uniformly mixed absorber rods for a square arrangement is presented. The Monte Carlo computer code MCNP4B with ENDF-B/VI library and detailed three dimensional geometry was used. WIMS-D code was used to model the isotopic composition of the fuel for 5, 10, 20 and 30 % burnup without cooling time. The results show, that out of the three studied effects: pitch from contact (3.75cm) up to rack design pitch (8cm), number of absorbers from 0 to 8 and burnup up to 30 %, the pitch has the greatest influence on the multiplication factor keff. In the interval in which the pitch was changed, keff decreased for up to 0.45. The number of absorber rods affects the multiplication factor much less. This effect is bigger for more compact arrangements, e.g. for contact of fuel elements with 8 absorber rods among them, keff values are smaller for almost 0.20 than for arrangement without absorber rods regardless of the burnup. The effect of burnup is the smallest since in no case keff decreases for more than 0.10, even for high burnups of 30 %. (author)
HAMCIND, Cell Burnup with Fission Products Poisoning
1 - Description of program or function: HAMCIND is a cell burnup code based in a coupling between HAMMER-TECHNION and CINDER. The fission product poisoning is taken into account in an explicit fashion. 2 - Method of solution: The nonlinear coupled set of equations for the neutron transport and nuclide transmutation equations and nuclide transmutation equations in a unit cell is solved by HAMCIND in a quasi-static approach. The spectral transport equation is solved by HAMMER-TECHNION at the beginning of each time-step while the nuclide transmutation equations are solved by CINDER for every time-step. The HAMMER-TECHNION spectral calculations are performed taking into account the fission product contribution to the macroscopic cross sections (fast and thermal), in the inelastic scattering matrix and even in the thermal scattering matrices. 3 - Restrictions on the complexity of the problem: Restrictions and/or limitations for HAMCIND depend upon the local operating system
Development and verification of Monte Carlo burnup calculation system
Monte Carlo burnup calculation code system has been developed to evaluate accurate various quantities required in the backend field. From the Actinide Research in a Nuclear Element (ARIANE) program, by using, the measured nuclide compositions of fuel rods in the fuel assemblies irradiated in the commercial Netherlands BWR, the analyses have been performed for the code system verification. The code system developed in this paper has been verified through analysis for MOX and UO2 fuel rods. This system enables to reduce large margin assumed in the present criticality analysis for LWR spent fuels. (J.P.N.)
Burnup credit implementation in WWER spent fuel management systems: Status and future aspects
This paper describes the motivation for possible burnup credit implementation in WWER spent fuel management systems in Bulgaria. The activities being done are described, namely: the development and verification of a 3D few-group diffusion burnup model; the application of the KORIGEN code for evaluation of WWER fuel nuclear inventory during reactor core lifetime and after spent fuel discharge; using the SCALE modular system (PC Version 4.1) for criticality safety analyses of spent fuel storage facilities. Future plans involving such important tasks as validation and verification of computer systems and libraries for WWER burnup credit analysis are shown. (author)
Determination of the burn-up of TRIGA fuel elements by calculation and reactivity experiments
The burnup of 17 fuel elements of the TRIGA Mark-II reactor in Vienna was measured. Different types of fuel elements had been simultaneously used for several years. The measured burnup values are compared with those calculated on the basis of core configuration and reactor operation history records since the beginning of operation. A one-dimensional, two-group diffusion computer code TRIGAP was used for the calculations. Comparison with burnup values determined by γ-scanning is also made. (orig./HP)
Taking burnup credit for interim storage and transportation system for BWR fuels
In order to establish a realistic burnup credit design system, a calculation system has been developed for determining isotope compositions, burnup, and criticality. The calculation system consists of several modules such as TGBLA, ORIGEN, CITATION, MCNP, and KENO. The TGBLA code is a fuel design code for LWR fuels developed in TOSHIBA Corporation. A compact measurement system for a fuel assembly has been being developed to meet requirements for the burnup determination, the neutron emission-rate evaluation, and the nuclear materials management. For a spent MOX fuel, a neutron emission rate measurement method has been being developed. The system consists of Cd-Te detectors and / or fission chambers. Some model calculations were carried out for the latest design BWR fuels. The effect of taking burnup credit for a transportation cask is shown. (authors)
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.
Technical Development on Burn-up Credit for Spent LWR Fuel
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
OREST - The hammer-origen burnup program system
Reliable prediction of the characteristics of irradiated light water reactor fuels (e.g., afterheat power, neutron and gamma radiation sources, final uranium and plutonium contents) is needed for many aspects of the nuclear fuel cycle. Two main problems must be solved: the simulation of all isotopic nuclear reactions and the simulation of neutron fluxes setting the reactions in motion. In state-of-the-art computer techniques, a combination of specialized codes for lattice cell and burnup calculations is preferred to solve these cross-linked problems in time or burnup step approximation. In the program system OREST, developed for official and commercial tasks in the Federal Republic of Germany nuclear fuel cycle, the well-known codes HAMMER and ORIGEN and directly coupled with a fuel rod temperature module
M. H. Altaf; N.H. Badrun
2014-01-01
Burnup dependent steady state thermal hydraulic analysis of TRIGA Mark-II research reactor has been carried out utilizing coupled point kinetics, neutronics and thermal hydraulics code EUREKA-2/RR. From the previous calculations of neutronics parameters including percentage burnup of individual fuel elements performed so far for 700 MWD burnt core of TRIGA reactor showed that the fuel rod predicted as hottest at the beginning of cycle (fresh core) was found to remain as the hottest until 200 ...
Ultrasonic measurement of high burn-up fuel elastic properties
The ultrasonic method developed for the evaluation of high burn-up fuel elastic properties is presented hereafter. The objective of the method is to provide data for fuel thermo-mechanical calculation codes in order to improve industrial nuclear fuel and materials or to design new reactor components. The need for data is especially crucial for high burn-up fuel modelling for which the fuel mechanical properties are essential and for which a wide range of experiments in MTR reactors and high burn-up commercial reactor fuel examinations have been included in programmes worldwide. To contribute to the acquisition of this knowledge the LAIN activity is developing in two directions. First one is development of an ultrasonic focused technique adapted to active materials study. This technique was used few years ago in the EdF laboratory in Chinon to assess the ageing of materials under irradiation. It is now used in a hot cell at ITU Karlsruhe to determine the elastic moduli of high burnup fuels from 0 to 110 GWd/tU. Some of this work is presented here. The second on going programme is related to the qualification of acoustic sensors in nuclear environments, which is of a great interest for all the methods, which work, in a hostile nuclear environment
Prediction of fission gas pressure from high burnup oxide fuel
The ELESIM fuel performance code incorporates a fundamentally based treatment of the relevant physical processes affecting fission gas release. The fission gas release model treats fission gas diffusion, formation and subsequent interlinkage of intergranular bubbles, grain boundary storage of gas, grain growth and fuel swelling. The latter case considers the contributions of thermal expansion, densification, solid fission products, and gas bubbles. The effect of porosity on fuel thermal conductivity is taken into account. Previously we showed predictions of the gas release model agreed well with measured values for oxide fuel with burnups to about 300 MW.h/kg U. The applicability of the model to high burnup fuel is examined using examples from the literature. The fission gas release range considered is about 1-100% for burnups to 1000 MW.h/kg U in thermal reactor fuel and 2400 MW.h/kg U in fast reactor fuel. Predicted and measured releases are shown to be in good agreement, suggesting that the fundamental model is correct. In some models, empirical correction factors are required at high burnup to achieve agreement between predicted and measured release values; no such factor is required in ELESIM. (auth)
Radionuclide Release from High Burnup Fuel
In this paper we investigate the production, evolution and release of radioactive fission products in a light water reactor. The production of the nuclides is determined by the neutronics, their evolution in the fuel by local temperature and by the fuel microstructure and the rate of release is governed by the scenario and the properties of the microstructure where the nuclides reside. The problem combines fields of reactor physics, fuel behaviour analysis and accident analysis. Radionuclide evolution during fuel reactor life is also important for determination of instant release fraction of final repository analysis. The source term problem is investigated by literature study and simulations with reactor physics code Serpent as well as fuel performance code ENIGMA. The capabilities of severe accident management codes MELCOR and ASTEC for describing high burnup structure effects are reviewed. As the problem is multidisciplinary in nature the transfer of information between the codes is studied. While the combining of the different fields as they currently are is challenging, there are some possibilities to synergy. Using reactor physics tools capable of spatial discretization is necessary for determining the HBS inventory. Fuel performance studies can provide insight how the HBS should be modelled in severe accident codes, however the end effect is probably very small considering the energetic nature of the postulated accidents in these scenarios. Nuclide release in severe accidents is affected by fuel oxidation, which is not taken into account by ANSI/ANS-5.4 but could be important in some cases, and as such, following the example of severe accident models would benefit the development of fuel performance code models. (author)
MCWO - Linking MCNP And ORIGEN2 For Fuel Burnup Analysis
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 Oxide fuel (Wg-MOX) fuel experiments in ATR and compares the MCWO-calculated results with the measured data
End effect analysis with various axial burnup distributions in high density spent fuel storage racks
Highlights: • Criticality tests are carried out with various axial burnup distributions of fuel assemblies for spent fuel storage racks. • KENO-Va code system was used to obtain criticalities with 10 axial segments. • ORIGEN-S code system was used to obtain burnup dependent axial compositions. • The criticality and burnup dependent reactivity difference are obtained from the results. • End effect quantifications are satisfactory confirming the previous suggestions. - Abstract: End effect of spent fuel comes from the difference between uniform and actual axial burnup distributions of fuel assemblies. It is significant to control the criticality safety in spent fuel storage and transportation. This work is focused on estimation of end effect in the spent fuel of light water reactor for the spent fuel storage rack region-II. High and low burnups of corresponding different uranium enrichments are taken into consideration to analyze the end effect with different axial burnup distributions such as uniform, MOC and EOC profiles. Two types of fuel assemblies such as CE type and Westinghouse type are considered. The whole calculations have been carried out by using the SCALE6 code including ORIGEN-S and KENO-Va
Barreto, Humberto
2015-01-01
This article is not the usual Excel pedagogy fare in that it does not provide an application or example taught via a spreadsheet. Instead, it briefly reviews the history of spreadsheets in the economics classroom and explores the current environment, with an emphasis on modern learning theory. The conclusion is not surprising: spreadsheets improve…
Measurement techniques for verifying burnup
Ewing, R.I. (Sandia National Lab., Albuquerque, NM (US)); Bierman, S.R. (Pacific Northwest Lab., Richland, WA (US))
1992-05-01
Measurements of the nuclear radiation from spent reactor fuel are being considered to qualify assemblies for loading into casks that will be used to transport spent fuel from utility sites to a federal storage facility. To ensure nuclear criticality safety, the casks are being designed to accept assemblies that meet restrictions as to burnup, initial enrichment and cooling time. This paper reports that measurements could be used to ensure that only fuel assemblies that meet the restrictions are selected for loading.
Measurement techniques for verifying burnup
Measurements of the nuclear radiation from spent reactor fuel are being considered to qualify assemblies for loading into casks that will be used to transport spent fuel from utility sites to a federal storage facility. To ensure nuclear criticality safety, the casks are being designed to accept assemblies that meet restrictions as to burnup, initial enrichment and cooling time. This paper reports that measurements could be used to ensure that only fuel assemblies that meet the restrictions are selected for loading
Burn-up dependent steady-state thermal hydraulic analysis of Pakistan research reactor-1
Muhammad Atta; Iqbal Masood; Mahmood Tayyab
2011-01-01
The burn-up dependent steady-state thermal hydraulic analysis of Pakistan research reactor-1, reference operating core, has been carried out utilizing standard computer codes WIMS/D4, CITATION, and RELAP5/MOD3.4. Reactor codes WIMS/D4 and CITATION have been used for the calculations of neutronic parameters including peaking factors and power profiles at different burn-up considering a xenon free core and also the equilibrium xenon values. RELAP5/MOD3.4 code was utilized for the determin...
Ostertag, Richard
2011-01-01
Focus of this thesis is to introduce methods used within the concept of Operational Excellence, being one of possible approaches to manage business processes, and subsequent application of these methods to accounting operations of a concrete company. Using application of these methods, thesis aims to evaluate quality of the process outputs as well as the utilization rate of available labor resources within analyzed business unit. Concrete suggestions to improve quality of process outputs and ...
Analysis on burnup step effect for evaluating reactor criticality and fuel breeding ratio
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
Consideration of the depletion phenomena and isotopic uncertainties in burnup-credit criticality analysis places an increasing reliance on computational tools and significantly increases the overall complexity of the calculations. An automated analysis and data management capability is essential for practical implementation of large-scale burnup credit analyses that can be performed in a reasonable amount of time. STARBUCS is a new prototypic analysis sequence being developed for the SCALE code system to perform automated criticality calculations of spent fuel systems employing burnup credit. STARBUCS is designed to help analyze the dominant burnup credit phenomena including spatial burnup gradients and isotopic uncertainties. A search capability also allows STARBUCS to iterate to determine the spent fuel parameters (e.g., enrichment and burnup combinations) that result in a desired keff for a storage configuration. Although STARBUCS was developed to address the analysis needs for spent fuel transport and storage systems, it provides sufficient flexibility to allow virtually any configuration of spent fuel to be analyzed, such as storage pools and reprocessing operations. STARBUCS has been used extensively at Oak Ridge National Laboratory (ORNL) to study burnup credit phenomena in support of the NRC Research program
TOPICAL REPORT ON ACTINIDE-ONLY BURNUP CREDIT FOR PWR SPENT NUCLEAR FUEL PACKAGES
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, keff, of a spent nuclear fuel package. Fifty-seven UO2, UO2/Gd2O3, and UO2/PuO2 critical experiments have been selected to cover anticipated conditions of SNF. The method uses an upper safety limit on keff (which can be a function of the trending parameters) such that the biased keff, 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 and confirm proper assembly selection prior to loading
Future disposal burnup credit process and effort
The United States Department of Energy's Office of Civilian Radioactive Waste Management has developed a risk-informed, performance based methodology for disposal criticality analyses. The methodology is documented in the Disposal Criticality Analysis Methodology Topical Report, YMP/TR-004Q (YMP 2000). The methodology includes taking credit for the burnup of irradiated commercial light water reactor fuel in criticality analyses, i.e., burnup credit. This paper summarizes the ongoing and planned future burnup credit activities associated with the methodology. (author)
Burnup credit activities in the United States
This report covers progress in burnup credit activities that have occurred in the United States of America (USA) since the International Atomic Energy Agency's (IAEA's) Advisory Group Meeting (AGM) on Burnup Credit was convened in October 1997. The Proceeding of the AGM were issued in April 1998 (IAEA-TECDOC-1013, April 1998). The three applications of the use of burnup credit that are discussed in this report are spent fuel storage, spent fuel transportation, and spent fuel disposal. (author)
New burnup calculation of TRIGA IPR-R1 reactor
The IPR-R1 TRIGA Mark I research reactor, located at the Nuclear Technology Development Center - CDTN, Belo Horizonte, Brazil, operates since 1960.The reactor is operating for more than fifty years and has a long history of operation. Determining the current composition of the fuel is very important to calculate various parameters. The reactor burnup calculation has been performed before, however, new techniques, methods, software and increase of the processing capacity of the new computers motivates new investigations to be performed. This work presents the evolution of effective multiplication constant and the results of burnup. This new model has a more detailed geometry with the introduction of the new devices, like the control rods and the samarium discs. This increase of materials in the simulation in burnup calculation was very important for results. For these series of simulations a more recently cross section library, ENDF/B-VII, was used. To perform the calculations two Monte Carlo particle transport code were used: Serpent and MCNPX. The results obtained from two codes are presented and compared with previous studies in the literature. (author)
The REBUS experimental programme for burn-up credit
An international programme called REBUS for the investigation of the burn-up credit has been initiated by the Belgian Nuclear Research Centre SCK·CEN and Belgonucleaire with the support of EdF and IRSN from France and VGB, representing German nuclear utilities and NUPEC, representing the Japanese industry. Recently also ORNL from the U.S. jointed the programme. The programme aims to establish a neutronic benchmark for reactor physics codes in order to qualify the codes for calculations of the burn-up credit. The benchmark exercise investigate the following fuel types with associated burn-up: reference fresh 3.3% enriched UO2 fuel, fresh commercial PWR UO2 fuel and irradiated commercial PWR UO2 fuel (54 GWd/tM), fresh PWR MOX fuel and irradiated PWR MOX fuel (20 GWd/tM). The experiments on the three configurations with fresh fuel have been completed. The experiments show a good agreement between calculation and experiments for the different measured parameters: critical water level, reactivity effect of the water level and fission-rate and flux distributions. In 2003 the irradiated BR3 MOX fuel bundle was loaded into the VENUS reactor and the associated experimental programme was carried out. The reactivity measurements in this configuration with irradiated fuel show a good agreement between experimental and preliminary calculated values. (author)
Fuel rod behaviour at high burnup WWER fuel cycles
The modernisation of WWER fuel cycles is carried out on the base of complete modelling and experimental justification of fuel rods up to 70 MWd/kgU. The modelling justification of the reliability of fuel rod and fuel rod with gadolinium is carried out with the use of certified START-3 code. START-3 code has a continuous experimental support. The thermophysical and strength reliability of WWER-440 fuel is justified for fuel rod and pellet burnups 65 MWd/kgU and 74 MWd/U, accordingly. Results of analysis are demonstrated by the example of uranium-gadolinium fuel assemblies of second generation under 5-year cycle with a portion of 6-year assemblies and by the example of successfully completed pilot operation of 5-year cycle fuel assemblies during 6 years at unit 3 of Kolskaja NPP. The thermophysical and strength reliability of WWER-1000 fuel is justified for a fuel rod burnup 66 MWd/kgU by the example of fuel operation under 4-year cycles and 6-year test operation of fuel assemblies at unit 1 of Kalininskaya NPP. By the example of 5-year cycle at Dukovany NPP Unit 2 it was demonstrated that WWER fuel rod of a burnup 58 MWd/kgU ensure reliable operation under load following conditions. The analysis has confirmed sufficient reserves of Russian fuel to implement program of JSC 'TVEL' in order to improve technical and economical parameters of WWER fuel cycles
New burnup calculation of TRIGA IPR-R1 reactor
Meireles, Sincler P. de; Campolina, Daniel de A.M.; Santos, Andre A. Campagnole dos; Menezes, Maria A.B.C.; Mesquita, Amir Z., E-mail: sinclercdtn@hotmail.com.br [Centro de Desenvolvimento da Tecnologia Nuclear (CDTN/CNEN-MG), Belo Horizonte, MG (Brazil)
2015-07-01
The IPR-R1 TRIGA Mark I research reactor, located at the Nuclear Technology Development Center - CDTN, Belo Horizonte, Brazil, operates since 1960.The reactor is operating for more than fifty years and has a long history of operation. Determining the current composition of the fuel is very important to calculate various parameters. The reactor burnup calculation has been performed before, however, new techniques, methods, software and increase of the processing capacity of the new computers motivates new investigations to be performed. This work presents the evolution of effective multiplication constant and the results of burnup. This new model has a more detailed geometry with the introduction of the new devices, like the control rods and the samarium discs. This increase of materials in the simulation in burnup calculation was very important for results. For these series of simulations a more recently cross section library, ENDF/B-VII, was used. To perform the calculations two Monte Carlo particle transport code were used: Serpent and MCNPX. The results obtained from two codes are presented and compared with previous studies in the literature. (author)
Burnup analysis of the VVER-1000 reactor using thorium-based fuel
Korkmaz, Mehmet E.; Agar, Osman; Bueyueker, Eylem [Karamanoglu Mehmetbey Univ., Karaman (Turkey). Faculty of Kamil Ozdag Science
2014-12-15
This paper aims to investigate {sup 232}Th/{sup 233}U fuel cycles in a VVER-1000 reactor through calculation by computer. The 3D core geometry of VVER-1000 system was designed using the Serpent Monte Carlo 1.1.19 Code. The Serpent Code using parallel programming interface (Message Passing Interface-MPI), was run on a workstation with 12-core and 48 GB RAM. {sup 232}Th/{sup 235}U/{sup 238}U oxide mixture was considered as fuel in the core, when the mass fraction of {sup 232}Th was increased as 0.05-0.1-0.2-0.3-0.4 respectively, the mass fraction of {sup 238}U equally was decreased. In the system, the calculations were made for 3 000 MW thermal power. For the burnup analyses, the core is assumed to deplete from initial fresh core up to a burnup of 16 MWd/kgU without refuelling considerations. In the burnup calculations, a burnup interval of 360 effective full power days (EFPDs) was defined. According to burnup, the mass changes of the {sup 232}Th, {sup 233}U, {sup 238}U, {sup 237}Np, {sup 239}Pu, {sup 241}Am and {sup 244}Cm were evaluated, and also flux and criticality of the system were calculated in dependence of the burnup rate.
Burnup analysis of the VVER-1000 reactor using thorium-based fuel
This paper aims to investigate 232Th/233U fuel cycles in a VVER-1000 reactor through calculation by computer. The 3D core geometry of VVER-1000 system was designed using the Serpent Monte Carlo 1.1.19 Code. The Serpent Code using parallel programming interface (Message Passing Interface-MPI), was run on a workstation with 12-core and 48 GB RAM. 232Th/235U/238U oxide mixture was considered as fuel in the core, when the mass fraction of 232Th was increased as 0.05-0.1-0.2-0.3-0.4 respectively, the mass fraction of 238U equally was decreased. In the system, the calculations were made for 3 000 MW thermal power. For the burnup analyses, the core is assumed to deplete from initial fresh core up to a burnup of 16 MWd/kgU without refuelling considerations. In the burnup calculations, a burnup interval of 360 effective full power days (EFPDs) was defined. According to burnup, the mass changes of the 232Th, 233U, 238U, 237Np, 239Pu, 241Am and 244Cm were evaluated, and also flux and criticality of the system were calculated in dependence of the burnup rate.
Actinide-Only Burnup Credit for PWR Spent Nuclear Fuel Packages
The objective of this topical report is to present to the NRC for review and acceptance a methodology for using burnup credit in the design of criticality control systems for PWR spent fuel transportation packages, while maintaining the criticality safety margins and related requirements of 10 CFR Part 71 and 72. The proposed methodology consists of five major steps as summarized below: (1) Validate a computer code system to calculate isotopic concentrations in SNF created during burnup in the reactor core and subsequent decay. (2) Validate a computer code system to predict the subcritical multiplication factor, keff, of a spent nuclear fuel package. (3) Establish bounding conditions for the isotopic concentration and criticality calculations. (4) Use the validated codes and bounding conditions to generate package loading criteria (burnup credit loading curves). and (5) Verify that SNF assemblies meet the package loading criteria and confirm proper fuel assembly selection prior to loading. (This step is required but the details are outside the scope of this topical report.) When reviewed and accepted by the NRC, this topical report will serve as a criterion document for criticality control analysts and will provide steps for the use of actinide-only burnup credit in the design of criticality control systems. The NRC-accepted burnup credit methodology will be used by commercial SNF storage and transportation package designers. Design-specific burnup credit criticality analyses will be defined, developed, and documented in the Safety Analysis Report (SAR) for each specific storage or transportation package that uses burnup credit. These SARs will then be submitted to the NRC for review and approval. This topical report is expected to be referenced in a number of storage and transportation cask applications to be submitted by commercial cask and canister designers to the NRC. Therefore, NRC acceptance of this topical report will result in increased efficiency of the
Topical Report on Actinide-Only Burnup Credit for PWR Spent Nuclear Fuel Packages. Revision 2
None, None
1998-09-01
The objective of this topical report is to present to the NRC for review and acceptance a methodology for using burnup credit in the design of criticality control systems for PWR spent fuel transportation packages, while maintaining the criticality safety margins and related requirements of 10 CFR Part 71 and 72. The proposed methodology consists of five major steps as summarized below: (1) Validate a computer code system to calculate isotopic concentrations in SNF created during burnup in the reactor core and subsequent decay. (2) Validate a computer code system to predict the subcritical multiplication factor, keff, of a spent nuclear fuel package. (3) Establish bounding conditions for the isotopic concentration and criticality calculations. (4) Use the validated codes and bounding conditions to generate package loading criteria (burnup credit loading curves). and (5) Verify that SNF assemblies meet the package loading criteria and confirm proper fuel assembly selection prior to loading. (This step is required but the details are outside the scope of this topical report.) When reviewed and accepted by the NRC, this topical report will serve as a criterion document for criticality control analysts and will provide steps for the use of actinide-only burnup credit in the design of criticality control systems. The NRC-accepted burnup credit methodology will be used by commercial SNF storage and transportation package designers. Design-specific burnup credit criticality analyses will be defined, developed, and documented in the Safety Analysis Report (SAR) for each specific storage or transportation package that uses burnup credit. These SARs will then be submitted to the NRC for review and approval. This topical report is expected to be referenced in a number of storage and transportation cask applications to be submitted by commercial cask and canister designers to the NRC. Therefore, NRC acceptance of this topical report will result in increased efficiency of the
Phenomena and Parameters Important to Burnup Credit
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
Phenomena and parameters important to burnup credit
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 United States 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. (author)
COREBN: A core burn-up calculation module for SRAC2006
COREBN is an auxiliary code of the SRAC system for multi-dimensional core burn-up calculation based on the diffusion theory and interpolation of macroscopic cross-sections tabulated to local parameters such as burn-up degree, moderator temperature and so on. The macroscopic cross-sections are prepared by cell burn-up calculations with the collision probability method of SRAC. SRAC and COREBN have wide applicability for various types of cell and core geometries. They have been used mainly for the purpose of core burn-up management of research reactors in Japan Atomic Energy Agency. The report is a revision of the users manual for the latest version of COREBN served with the SRAC released in 2006. (author)
OREST, LWR Burnup Simulation Using Program HAMMER and ORIGEN
1 - Description of program or function: In OREST, the 1-dimensional lattice code HAMMER and the isotope generation and depletion code ORIGEN are directly coupled for burnup simulation in light-water reactor fuels (GRS recommended). Additionally heavy water and graphite moderated systems can be calculated. New version differs from the previous version in the following features: An 84-group-library LIB84 for up to 200 isotopes is used to update the 3-group -POISON-XS. LIB84 uses the same energy boundaries as THERMOS and HAMLET in . In this way, high flexibility is achieved in very different reactor models. The coupling factor between THERMOS and HAMLET is now directly transferred from HAMMER to THERES and omits the equation 4 (see page 6 of the manual). Sandwich-reactor fuel reactivity and burnup calculations can be started with NGEOM = 1. Thorium graphite reactivity and burnup calculations can be started with NLIBE = 1. High enriched U-235 heavy water moderated reactivity and burnup calculations can be started. HAMLET libraries in for U-235, U-236, U-238, Np-237, Pu-238, Pu-239, Pu-240, Pu-242, Am-241, Am-243 and Zirconium are updated using resonance parameters. NEA-1324/04: A new version of the module hamme97.f has replaced the old one. 2 - Method of solution: For the user-defined irradiation history, an input data processor generates program loops over small burnup steps for the main codes HAMMER and ORIGEN. The user defined assembly description is transformed to an equivalent HAMMER fuel cell. HAMMER solves the integral neutron transport equation in a four-region cylindrical or sandwiched model with reflecting boundaries and runs with fuel power calculated rod temperatures. ORIGEN runs with HAMMER-calculated cross sections and neutron spectra and calculates isotope concentrations during burnup by solving the buildup-, depletion- and decay-chain equations. An output data processor samples the outputs of the program modules and generates tabular works for the
Burn-up dependent steady-state thermal hydraulic analysis of Pakistan research reactor-1
Muhammad Atta
2011-01-01
Full Text Available The burn-up dependent steady-state thermal hydraulic analysis of Pakistan research reactor-1, reference operating core, has been carried out utilizing standard computer codes WIMS/D4, CITATION, and RELAP5/MOD3.4. Reactor codes WIMS/D4 and CITATION have been used for the calculations of neutronic parameters including peaking factors and power profiles at different burn-up considering a xenon free core and also the equilibrium xenon values. RELAP5/MOD3.4 code was utilized for the determination of peak fuel centerline, clad and coolant temperatures to ensure the safety of the reactor throughout the cycle. The calculations reveal that the reactor is safe and no nucleate boiling will commence at any part of the core throughout the cycle and that the safety margin increases with burnup as peaking factors decrease.
Burn-up dependent steady-state thermal hydraulic analysis of Pakistan research reactor-1
The burn-up dependent steady-state thermal hydraulic analysis of Pakistan research reactor-1, reference operating core, has been carried out utilizing standard computer codes WIMS/D4, CITATION, and RELAP5/MOD3.4. Reactor codes WIMS/D4 and CITATION have been used for the calculations of neutronic parameters including peaking factors and power profiles at different burn-up considering a xenon free core and also the equilibrium xenon values. RELAP5/MOD3.4 code was utilized for the determination of peak fuel centerline, clad and coolant temperatures to ensure the safety of the reactor throughout the cycle. The calculations reveal that the reactor is safe and no nucleate boiling will commence at any part of the core throughout the cycle and that the safety margin increases with burnup as peaking factors decrease. (author)
For adopting burnup credit in transport or storage of spent fuel (SF), development of a reliable burnup calculation code is crucial. For this purpose, data of Post Irradiation Examination (PIE) have been extensively analyzed to evaluate accuracy of burnup calculation codes for a 14*14 or 15*15 PWR fuel assembly. This study shows results of analysis of this latest PIE with SWAT and ORIGEN2.1. SWAT is an integrated burnup code system for a 17*17 PWR fuel assembly that has been developed by Tohoku University and JAERI. The results show that SWAT can more precisely predict nuclide composition of latest PWR assembly than ORIGEN2.1. (O.M.)
Roman, Steven
2008-01-01
To achieve the maximum control and flexibility from Microsoft® Excel often requires careful custom programming using the VBA (Visual Basic for Applications) language. Writing Excel Macros with VBA, 2nd Edition offers a solid introduction to writing VBA macros and programs, and will show you how to get more power at the programming level: focusing on programming languages, the Visual Basic Editor, handling code, and the Excel object model.
Simulation of the behaviour of nuclear fuel under high burnup conditions
Highlights: • Increasing the time of nuclear fuel into reactor generates high burnup structure. • We analyze model to simulate high burnup scenarios for UO2 nuclear fuel. • We include these models in the DIONISIO 2.0 code. • Tests of our models are in very good agreement with experimental data. • We extend the range of predictability of our code up to 60 MWd/KgU average. - Abstract: In this paper we summarize all the models included in the latest version of the DIONISIO code related to the high burnup scenario. Due to the extension of nuclear fuels permanence under irradiation, physical and chemical modifications are developed in the fuel material, especially in the external corona of the pellet. The codes devoted to simulation of the rod behaviour under irradiation need to introduce modifications and new models in order to describe those phenomena and be capable to predict the behaviour in all the range of a general pressurized water reactor. A complex group of subroutines has been included in the code in order to predict the radial distribution of power density, burnup, concentration of diverse nuclides and porosity within the pellet. The behaviour of gadolinium as burnable poison also is modelled into the code. The results of some of the simulations performed with DIONISIO are presented to show the good agreement with the data selected for the FUMEX I/II/III exercises, compiled in the NEA data bank
Syncrude Canada Ltd., operator of the oil sands mine and processing plant near Fort McMurray, Alberta, produces 11% of Canada's crude oil and is the country's largest private-sector employer of native Canadians. Syncrude has the goal of employing about 10% native Canadians, which is about the percentage of natives in the regional population. Examples are presented of successful native employment and entrepreneurship at Syncrude. Doreen Janvier, once employed at Syncrude's mine wash bays, was challenged to form her own company to contract out labor services. Her company, DJM Enterprises, now has a 2-year contract to operate three highly sophisticated wash bays used to clean mining equipment, and is looking to bid on other labor contracts. Mabel Laviolette serves as liaison between the oil containment and recovery team, who recover oil skimmed off Syncrude's tailings basin, and the area manager. The team approach and the seasonal nature of the employment fit in well with native cultural patterns. The excellence of native teamwork is also illustrated in the mine rescue team, one unit of which is entirely native Canadian. Part of Syncrude's aboriginal policy is to encourage development of aboriginal enterprises, such as native-owned Clearwater Welding and Fabricating Ltd., which has held welding and fabricating contracts with most major companies in the region and is a major supplier of skilled tradesmen to Syncrude. Syncrude also provides employment and training, encourages natives to continue their education, and promotes local community development. 4 figs
Tsoupikova, Daria
2007-02-01
This paper describes the research and development of a virtual reality visualization project "Passing excellence" about the world famous architectural ensemble "Kizhi". The Kizhi Pogost is located on an island in Lake Onega in northern Karelia in Russia. It is an authentic museum of an ancient wood building tradition which presents a unique artistic achievement. This ensemble preserves a concentration of masterpieces of the Russian heritage and is included in the List of Most Endangered Sites of the World Monuments Watch protected by World Heritage List of UNESCO. The project strives to create a unique virtual observation of the dynamics of the architectural changes of the museum area beginning from the 15th Century up to the 21st Century. The visualization is being created to restore the original architecture of Kizhi island based on the detailed photographs, architectural and geometric measurements, textural data, video surveys and resources from the Kizhi State Open-Air Museum archives. The project is being developed using Electro, an application development environment for the tiled display high-resolution graphics visualization system and can be shown on the virtual reality systems such as the GeoWall TM and the C-Wall.
Plutonium and Minor Actinides Recycling in Standard BWR using Equilibrium Burnup Model
Abdul Waris
2008-03-01
Full Text Available Plutonium (Pu and minor actinides (MA recycling in standard BWR with equilibrium burnup model has been studied. We considered the equilibrium burnup model as a simple time independent burnup method, which can manage all possible produced nuclides in any nuclear system. The equilibrium burnup code was bundled with a SRAC cell-calculation code to become a coupled cell-burnup calculation code system. The results show that the uranium enrichment for the criticality of the reactor, the amount of loaded fuel and the required natural uranium supply per year decrease for the Pu recycling and even much lower for the Pu & MA recycling case compared to those of the standard once-through BWR case. The neutron spectra become harder with the increasing number of recycled heavy nuclides in the reactor core. The total fissile rises from 4.77% of the total nuclides number density in the reactor core for the standard once-through BWR case to 6.64% and 6.72% for the Plutonium recycling case and the Pu & MA recycling case, respectively. The two later data may become the main basis why the required uranium enrichment declines and consequently diminishes the annual loaded fuel and the required natural uranium supply. All these facts demonstrate the advantage of plutonium and minor actinides recycling in BWR.
Results of the isotopic concentrations of VVER calculational burnup credit benchmark no. 2(cb2
The characterization of the irradiated fuel materials is becoming more important with the Increasing use of nuclear energy in the world. The purpose of this document is to present the results of the nuclide concentrations calculated Using Calculation VVER Burnup Credit Benchmark No. 2(CB2). The calculations were Performed in The Nuclear Technology Center of Cuba. The CB2 benchmark specification as the second phase of the VVER burnup credit benchmark is Summarized in [1]. The CB2 benchmark focused on VVER burnup credit study proposed on the 97' AER Symposium [2]. It should provide a comparison of the ability of various code systems And data libraries to predict VVER-440 spent fuel isotopes (isotopic concentrations) using Depletion analysis. This phase of the benchmark calculations is still in progress. CB2 should be finished by summer 1999 and evaluated results could be presented on the next AER Symposium. The obtained results are isotopic concentrations of spent fuel as a function of the burnup and Cooling time. The depletion point ORIGEN2[3] code was used for the calculation of the spent Fuel concentration. The depletion analysis was performed using the VVER-440 irradiated fuel assemblies with in-core Irradiation time of 3 years, burnup of the 30000 mwd/TU, and an after discharge cooling Time of 0 and 1 year. This work also comprises the results obtained by other codes[4].
A simplified burnup calculation strategy with refueling in static molten salt reactor
Molten Salt Reactors, by nature can be refuelled and reprocessed online. Thus, a simulation methodology has to be developed which can consider online refueling and reprocessing aspect of the reactor. To cater such needs a simplified burnup calculation strategy to account for refueling and removal of molten salt fuel at any desired burnup has been identified in static molten salt reactor in batch mode as a first step of way forward. The features of in-house code ITRAN has been explored for such calculations. The code also enables us to estimate the reactivity introduced in the system due to removal of any number of considered nuclides at any burnup. The effect of refueling fresh fuel and removal of burned fuel has been studied in batch mode with in-house code ITRAN. The effect of refueling and burnup on change in reactivity per day has been analyzed. The analysis of removal of 233Pa at a particular burnup has been carried out. The similar analysis has been performed for some other nuclides also. (author)
Issues for effective implementation of burnup credit
In the United States, burnup credit has been used in the criticality safety evaluation for storage pools at pressurized water reactors (PWRs) and considerable work has been performed to lay the foundation for use of burnup credit in dry storage and transport cask applications and permanent disposal applications. Many of the technical issues related to the basic physics phenomena and parameters of importance are similar in each of these applications. However, the nuclear fuel cycle in the United States has never been fully integrated and the implementation of burnup credit to each of these applications is dependent somewhat on the specific safety bases developed over the history of each operational area. This paper will briefly review the implementation status of burnup credit for each application area and explore some of the remaining issues associated with effective implementation of burnup credit. (author)
MTR core loading pattern optimization using burnup dependent group constants
Iqbal Masood
2008-01-01
Full Text Available A diffusion theory based MTR fuel management methodology has been developed for finding superior core loading patterns at any stage for MTR systems, keeping track of burnup of individual fuel assemblies throughout their history. It is based on using burnup dependent group constants obtained by the WIMS-D/4 computer code for standard fuel elements and control fuel elements. This methodology has been implemented in a computer program named BFMTR, which carries out detailed five group diffusion theory calculations using the CITATION code as a subroutine. The core-wide spatial flux and power profiles thus obtained are used for calculating the peak-to-average power and flux-ratios along with the available excess reactivity of the system. The fuel manager can use the BFMTR code for loading pattern optimization for maximizing the excess reactivity, keeping the peak-to-average power as well as flux-ratio within constraints. The results obtained by the BFMTR code have been found to be in good agreement with the corresponding experimental values for the equilibrium core of the Pakistan Research Reactor-1.
High burnup fuel development program in Japan
A step wise burnup extension program has been progressing in Japan to reduce the LWR fuel cycle cost. At present, the maximum assembly burnup limit of BWR 8 Χ 8 type fuel (B. Step II fuel) is 50GWd/t and a limited numbers of 9 Χ 9 type fuel (B. Step III fuel) with 55GWd/t maximum assembly burnup has been licensed by regulatory agencies recently. Though present maximum assembly burnup limit for PWR fuel is 48GWd/t (P. Step I fuel), the licensing work has been progressing for irradiation testing on a limited number of fuel assemblies with extended burnup of up to 55GWd/t (p. Step II fuel) Design of high burnup fuel and fabrication test are carried out by vendors, and subsequent irradiation test of fuel rods is conducted jointly by utilities and vendors to prepare for licensing. It is usual to make an irradiation test for vectarion, using lead use assemblies by government to confirm fuel integrity and reliability and win the public confidence. Nuclear Power Engineering Corporation (NUPE C) is responsible for verification test. The fuel are subjected to post irradiation examination (PIE) and no unfavorable indications of fuel behavior have found both in NUPE C verification test and joint irradiation test by utilities and vendors. Burnup extension is an urgent task for LWR fuel in Japan in order to establish the domestic fuel cycle. It is conducted in joint efforts of industries, government and institutes. However, watching a situation of burnup extension in the world, we are not going ahead of other countries in the achievement of burnup extension. It is due to a conservative policy in the nuclear safety of the country. This is the reason why the burnup extension program in Japan is progressing 'slow and steady' As for the data obtained, no unfavorable indications of fuel behavior have found both in NUPE C verification test and joint irradiation test by utilities and vendors until now
The implementation of burnup credit in VVER-440 spent fuel
The countries using Russian reactors VVER-440 cooperate in reactor physics in Atomic Energy Research (AER). One of topic areas is 'Physical Problems of Spent Fuel, Radwaste and Decommissioning' (Working Group E). In this article, in the first part is an overview about our activity for numerical and experimental verification of codes which participants use for calculation of criticality, isotopic concentration, activity, neutron and gamma sources and shielding is shown. The set of numerical benchmarks (CB1, CB2, CB3 and CB4) is very similar (the same idea, the VVER-440) to the OECD/NEA/NSC Burnup Credit Criticality Benchmarks, Phases 1 and 2. In the second part, verification of the SCALE 4.4 system (only criticality and nuclide concentrations) for VVER-440 fuel is shown. In the third part, dependence of criticality on burnup (only actinides and actinides + fission products) for transport cask C30 with VVER-440 fuel by optimal moderation is shown. In the last part, current status in implementation burnup credit in Slovakia is shown. (author)
Value of 236U to actinide-only burnup credit
The US Department of Energy (DOE) submitted a topical report to the US Nuclear Regulatory Commission (NRC) in May 1995 in order to gain approval of a method for criticality analysis of transport packages that takes account for the change in actinide isotopes with burnup [pressurized water reactors (PWRs) only]. Historically, the NRC has conservatively assumed that the fuel was in its initial conditions (without any burnable absorbers). In order to permit credit for the changes in actinide content, the NRC has required validation of the depletion and criticality codes for spent nuclear fuel, justification of conservative depletion modeling, and finally confirmation measurements before loading. The NRC requested additional information on March 22, 1996. The DOE responded by a revision of the topical report in May 1997. The NRC again responded with another set of requests of additional information in April 1998. In that set of questions, the NRC challenged the use of 236U in burnup credit. Uranium-236 is not found in any significant amount in any available critical experiments. The authors explore the value of 236U to actinide-only burnup credit
BISON, 1-D Burnup and Transport in Slab, Cylindrical, Spherical Geometry
1 - Description of problem or function: BISON-1.5 solves the one- dimensional Boltzmann transport equation for neutron and gamma-rays and transmutation equations for fuel nuclides. 2 - Method of solution: In the transport calculation stage the one- dimensional Boltzmann transport equation is solved by the discrete ordinates method. In the burnup calculation stage, transmutation equations for fuel nuclides are solved by Bateman's method. The neutron flux obtained in the transport calculation stage is used to determine the transmutation rates in the burnup calculation stage. Both stages are repeated in tandem till the end of the burnup cycle. 3 - Restrictions on the complexity of the problem: A 42-group neutron and 21-group gamma-ray cross section library is prepared in the code package. Core storage for array variables is dynamically allocated by the code, so there are no restrictions on the size of each array
The impact of radial and axial burnup profile on the criticality of WWER-440 spent fuel cask is presented in the paper. The calculations are performed based on two AER Benchmark problems for WWER-440 irradiated fuel assembly. The radial zonewise dependent spent fuel inventory has been calculated by the NESSEL - NUKO code system. The axial dependent isotope concentrations have been determined by the modular code system SCALE4.4. For criticality calculations the SCALE4.4 has been applied. Calculations have been carried out for cask with 30 WWER-440 fuel assemblies with initial enrichment 3.6% of 235U and burnup up to 40 MWd/kgU. The influence of radial and axial burnup credit on the cask criticality has been evaluated
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.