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)
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
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
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
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)
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
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.
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)
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
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
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)
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
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
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
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.
When the residence time of nuclear fuel rods exceeds a given threshold value, several properties of the pellet material suffer changes and hence the posterior behaviour of the rod is significantly altered. Structural modifications start at the pellet periphery, which is usually referred to as rim zone. It is presently believed that these changes are a consequence of the localized absorption of epithermal neutrons by 238U, which effective cross section presents resonant peaks. Due to the chain of nuclear reactions that take place, several Pu isotopes are born especially at the rim. In particular, the fissile character of 239Pu and 241Pu is the cause of the increased number of fission events that occur in the pellet periphery. For this reason, the power generation rate and the burnup adopt a non uniform distribution in the pellet, reaching at the rim values two or three times higher than the average [1]. The rim zone starts to form for a burnup threshold value of about 50-60 MWd/kgHM and its width increases as the irradiation progresses. The microstructure of this zone is characterized by the presence of small grains, with a typical size of 200 nm, and large pores, of some μm. Even though the rim zone is very thin, it has a significant effect on the mechanical integrity of the pellet, particularly when it makes contact with the cladding, and on the temperature distribution in the whole pellet, because of its low thermal conductivity [1,2]. The numerical codes designed to simulate fuel behaviour under irradiation must include the phenomena associated to high burnup if they aim at extending the prediction range, and this is the purpose with our DIONISIO code. But a detailed analysis of the phenomena that take place in this region demands the use of neutronic codes that solve the Boltzmann transport equations [3] in a number of energy intervals (groups), including adequate considerations in the region of the resonant absorption peaks of 238U. These cell codes predict
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)
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)
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.
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.
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
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.)
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)
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
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
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%.
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 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)
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
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.)
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.
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
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)
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
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)
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)
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)
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)
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
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
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)
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
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
Implementation of burnup credit in spent fuel management systems
Improved calculational methods allow one to take credit for the reactivity reduction associated with fuel burnup. This means reducing the analysis conservatism while maintaining an adequate safety margin. The motivation for using burnup credit in criticality safety applications is based on economic considerations and additional benefits contributing to public health and safety and resource conservation. Interest in the implementation of burnup credit has been shown by many countries. In 1997, the International Atomic Energy Agency (IAEA) started a task to monitor the implementation of burnup credit in spent fuel management systems, to provide a forum to exchange information, to discuss the matter and to gather and disseminate information on the status of national practices of burnup credit implementation in the Member States. The task addresses current and future aspects of burnup credit. This task was continued during the following years. (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)
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)
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)
CARMEN-SYSTEM, Programs System for Thermal Neutron Diffusion and Burnup with Feedback
1 - Description of problem or function: CARMEN is a system of programs developed for the neutronic calculation of PWR cycles. It includes the whole chain of analysis from cell calculations to core calculations with burnup. The core calculations are based on diffusion theory with cross sections depending on the relevant space-dependent feedback effects which are present at each moment along the cycles. The diffusion calculations are in one, two or three dimensions and in two energy groups. The feedback effects which are treated locally are: burnup, water density, power density and fission products. In order to study in detail these parameters the core should be divided into as many zones as different cross section sets are expected to be required in order to reproduce reality correctly. A relevant difference in any feedback parameter between zones produces different cross section sets for the corresponding zones. CARMEN is also capable to perform the following calculations: - Multiplication factor by burnup step with fixed boron concentration - Buckling and control rod insertion - Buckling search by burnup step - Boron search by burnup step - Control rod insertion search by burnup step. 2 - Method of solution: The cell code (LEOPARD-TRACA) generates the fuel assembly cross sections versus burnup. This is the basic library to be used in the CARMEN code proper. With a planar distribution guess for power density, water density and fluxes, the macroscopic cross sections by zone are calculated by CARMEN, and then a diffusion calculation is done in the whole geometry. With the distribution of power density, heat accumulated in the coolant and the thermal and fast fluxes determined in the diffusion calculation, CARMEN calculates the values of the most relevant parameters that influence the macroscopic cross sections by zone: burnup, water density, effective fuel temperature and fission product concentrations. If these parameters by zone are different from the reference
Validation of BGCore System for Burnup Calculations
BGCore is a software package for comprehensive computer simulation of nuclear reactor systems and their fuel cycles. BGCore interfaces the Monte Carlo particles transport code MCNP4C with a SARAF module - an independently developed code for calculating fuel composition during irradiation and spent fuel emissions following discharge. In BGCore system, depletion coupling methodology is based on the multi-group approach that significantly reduces computation time and allows tracking of large number of nuclides during calculations. The objective of this study is validation of the BGCore system against well established and verified, state of the art computer codes for thermal and fast spectrum lattices
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
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.
A burn-up module coupling to an AMPX system
The Reactors and Neutrons Division of the Bariloche Atomic Center uses the AMPX system for the study of high conversion reactors (HCR). Such system allows to make neutronic calculations from the nuclear data library (ENDF/B-IV). The Nuclear Engineering career of the Balseiro Institute developed and implemented a burn-up module at a μ-cell level (BUM: Burn-up Module) which agrees with the requirement to be coupled to the AMPX system. (Author)
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.
Sophistication of burnup analysis system for fast reactor
Improvement on prediction accuracy for neutronics property of fast reactor cores is one of the most important study domains in terms of both achievement of high economical plant efficiency based on reasonably advanced designs and increased reliability and safety margins. In former study, considerable improvement on prediction accuracy in neutronics design has been achieved in the development of the unified constants library as a fruit of a series of critical experiments such as JUPITER in application of the reactor constant adjustments. For design of fast reactor cores, however, improvement of not only static properties but also burnup properties is very important. For such purpose, it is necessary to improve the prediction accuracy on burnup properties using actual burnup data of 'JOYO' and 'MONJU', experimental and prototype fast reactors. Recently, study on effective burnup method for minor actinides becomes important theme. However, there is a problem that analysis work tends to become inefficient for lack of functionality suitable for analysis of composition change due to burnup since the conventional analysis system is targeted to critical assembly systems. Therefore development of burnup analysis system for fast reactors with modularity and flexibility is being done that would contribute to actual core design work and improvement of prediction accuracy. In the previous research, we have developed a prototype system which has functions of performing core and burnup calculations using given constant files (PDS files) and information based on simple and easy user input data. It has also functions of fuel shuffling which is indispensable for production systems. In the present study, we implemented functions for cell calculations and burnup calculations. With this, whole steps in analysis can be carried out with only this system. In addition, we modified the specification of user input to improve the convenience of this system. Since implementations being done so
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)
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)
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)
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)
Automated system for determining the burnup of spent nuclear fuel
Mokritskii V. A.
2014-12-01
Full Text Available The authors analyze their experience in application of semi-conductor detectors and development of a breadboard model of the monitoring system for spent nuclear fuel (SNF. Such system should use CdZnTe-detectors in which one-charging gathering conditions are realized. The proposed technique of real time SNF control during reloading technological operations is based on the obtained research results. Methods for determining the burnup of spent nuclear fuel based on measuring the characteristics of intrinsic radiation are covered in many papers, but those metods do not usually take into account that the nuclear fuel used during the operation has varying degrees of initial enrichment, or a new kind of fuel may be used. Besides, the known methods often do not fit well into the existing technology of fuel loading operations and are not suitable for operational control. Nuclear fuel monitoring (including burnup determination system in this research is based on the measurement of the spectrum of natural gamma-radiation of irradiated fuel assemblies (IFA, as from the point of view of minimizing the time spent, the measurement of IFA gamma spectra directly during fuel loading is optimal. It is the overload time that is regulated rather strictly, and burnup control operations should be coordinated with the schedule of the fuel loading. Therefore, the real time working capacity of the system should be chosen as the basic criterion when constructing the structure of such burnup control systems.
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)
Sophistication of burnup analysis system for fast reactor (2)
Improvement on prediction accuracy for neutronics characteristics of fast reactor cores is one of the most important study domains in terms of both achievement of high economical plant efficiency based on reasonably advanced designs and increased reliability and safety margins. In former study, considerable improvement on prediction accuracy in neutronics design has been achieved in the development of the unified cross-section set as a fruit of a series of critical experiments such as JUPITER in application of the reactor constant adjustments. For design of fast reactor cores improvement of not only static characteristics but also burnup characteristics is very important. For such purpose, it is necessary to improve the prediction accuracy on burnup characteristics using actual burnup data of 'JOYO' and 'MONJU', experimental and prototype fast reactors. Recently, study on effective burnup method for minor actinides becomes important theme. However, there is a problem that analysis work tends to become inefficient for lack of functionality suitable for analysis of composition change due to burnup since the conventional analysis system is targeted to critical assembly systems. Therefore development of burnup analysis system for fast reactors with modularity and flexibility is being done that would contribute to actual core design work and improvement of prediction accuracy. In the previous study, we have developed a prototype system which has functions of performing core and burnup calculations using given constant files (PDS files) and information based on simple and easy user input data. It has also functions of fuel shuffling which is indispensable for power reactor analysis systems. In the present study, by extending the prototype system, features for handling of control rods and energy collapse of group constants have been designed and implemented. Computational results from the present analysis system are stored into restart files which can be accessible by
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)
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)
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.
New Burnup Calculation System for Fusion-Fission Hybrid System
Investigation of nuclear waste incineration has positively been carried out worldwide from the standpoint of environmental issues. Some candidates such as ADS, FBR are under discussion for possible incineration technology. Fusion reactor is one of such technologies, because it supplies a neutron-rich and volumetric irradiation field, and in addition the energy is higher than nuclear reactor. However, it is still hard to realize fusion reactor right now, as well known. An idea of combination of fusion and fission concepts, so-called fusion-fission hybrid system, was thus proposed for the nuclear waste incineration. Even for a relatively lower plasma condition, neutrons can be well multiplied by fission in the nuclear fuel, tritium is thus bred so as to attain its self-sufficiency, enough energy multiplication is then expected and moreover nuclear waste incineration is possible. In the present study, to realize it as soon as possible with the presently proven technology, i.e., using ITER model with the achieved plasma condition of JT60 in JAEA, Japan, a new calculation system for fusion-fission hybrid reactor including transport by MCNP and burnup by ORIGEN has been developed for the precise prediction of the neutronics performance. The author's group already has such a calculation system developed by them. But it had a problem that the cross section libraries in ORIGEN did not have a cross section library, which is suitable specifically for fusion-fission hybrid reactors. So far, those for FBR were approximately used instead in the analysis. In the present study, exact derivation of the collapsed cross section for ORIGEN has been investigated, which means it is directly evaluated from calculated track length by MCNP and point-wise nuclear data in the evaluated nuclear data file like JENDL-3.3. The system realizes several-cycle calculation one time, each of which consists of MCNP criticality calculation, MCNP fixed source calculation with a 3-dimensional precise
Marshall, William BJ J [ORNL; Ade, Brian J [ORNL; Bowman, Stephen M [ORNL; Gauld, Ian C [ORNL; Ilas, Germina [ORNL; Mertyurek, Ugur [ORNL; Radulescu, Georgeta [ORNL
2015-01-01
Oak Ridge National Laboratory and the United States Nuclear Regulatory Commission have initiated a multiyear project to investigate application of burnup credit for boiling-water reactor (BWR) fuel in storage and transportation casks. This project includes two phases. The first phase (1) investigates applicability of peak reactivity methods currently used in spent fuel pools (SFPs) to storage and transportation systems and (2) evaluates validation of both reactivity (k_{eff}) calculations and burnup credit nuclide concentrations within these methods. The second phase will focus on extending burnup credit beyond peak reactivity. This paper documents the first phase, including an analysis of lattice design parameters and depletion effects, as well as both validation components. Initial efforts related to extended burnup credit are discussed in a companion paper. Peak reactivity analyses have been used in criticality analyses for licensing of BWR fuel in SFPs over the last 20 years. These analyses typically combine credit for the gadolinium burnable absorber present in the fuel with a modest amount of burnup credit. Gadolinium burnable absorbers are used in BWR assemblies to control core reactivity. The burnable absorber significantly reduces assembly reactivity at beginning of life, potentially leading to significant increases in assembly reactivity for burnups less than 15–20 GWd/MTU. The reactivity of each fuel lattice is dependent on gadolinium loading. The number of gadolinium-bearing fuel pins lowers initial lattice reactivity, but it has a small impact on the burnup and reactivity of the peak. The gadolinium concentration in each pin has a small impact on initial lattice reactivity but a significant effect on the reactivity of the peak and the burnup at which the peak occurs. The importance of the lattice parameters and depletion conditions are primarily determined by their impact on the gadolinium depletion. Criticality code validation for BWR burnup
Use of burnup credit in criticality safety design analysis of spent fuel storage systems
temperature and density, presence of soluble boron in the core (PWR), use of fixed neutron absorbers (control rods, burnable poison rods, axial power shaping rods), use of integral burnable absorbers (gadolinium or erbium bearing fuel rods, IFBA rods). It will be shown how a bounding approach can be obtained for the impact of these parameters on the reactivity of the storage system. The criticality calculation procedure consists in the following main steps: Isotopic selection and validation; Validation of the criticality calculation code applied; Sensitivity studies on the reactivity effects of axial and horizontal burnup profiles of fuel assemblies; Determination of the criticality acceptance criterion (maximum allowable neutron multiplication factor including the impacts of all the mechanical and calculational uncertainties) and determination of the loading curve. The fundamentals of isotopic selection will be defined, and a survey of the benchmark experiments available for isotopic validation and validation of the criticality calculation code applied will be given. Since the parameters and conditions characterizing the benchmark experiments are usually different from the parameters and conditions describing the spent fuel storage system of interest, a method of checking the applicability of such experiments to the storage system will be briefly described. This method bases the applicability on the similarity of sensitivity coefficients which are defined for the underlying nuclear data characterizing the isotopic compositions and their effect on the spent fuel reactivity. The fact that the axial burnup distribution in a fuel assembly is non-uniform must be considered in the analysis of the storage system. The difference between the system's neutron multiplication factor obtained by using an axially varying burnup profile and the system's neutron multiplication factor obtained by assuming a uniform distribution of the averaged burnup of this profile is known as the 'end
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)
Garcia-Herranz, N.; Cabellos, O. [Madrid Polytechnic Univ., Dept. of Nuclear Engineering (Spain); Cabellos, O.; Sanz, J. [Madrid Polytechnic Univ., 2 Instituto de Fusion Nuclear (Spain); Sanz, J. [Univ. Nacional Educacion a Distancia, Dept. of Power Engineering, Madrid (Spain)
2005-07-01
We present a new code system which combines the Monte Carlo neutron transport code MCNP-4C and the inventory code ACAB as a suitable tool for high burnup calculations. Our main goal is to show that the system, by means of ACAB capabilities, enables us to assess the impact of neutron cross section uncertainties on the inventory and other inventory-related responses in high burnup applications. The potential impact of nuclear data uncertainties on some response parameters may be large, but only very few codes exist which can treat this effect. In fact, some of the most reported effective code systems in dealing with high burnup problems, such as CASMO-4, MCODE and MONTEBURNS, lack this capability. As first step, the potential of our system, ruling out the uncertainty capability, has been compared with that of those code systems, using a well referenced high burnup pin-cell benchmark exercise. It is proved that the inclusion of ACAB in the system allows to obtain results at least as reliable as those obtained using other inventory codes, such as ORIGEN2. Later on, the uncertainty analysis methodology implemented in ACAB, including both the sensitivity-uncertainty method and the uncertainty analysis by the Monte Carlo technique, is applied to this benchmark problem. We estimate the errors due to activation cross section uncertainties in the prediction of the isotopic content up to the high-burnup spent fuel regime. The most relevant uncertainties are remarked, and some of the most contributing cross sections to those uncertainties are identified. For instance, the most critical reaction for Am{sup 242m} is Am{sup 241}(n,{gamma}-m). At 100 MWd/kg, the cross-section uncertainty of this reaction induces an error of 6.63% on the Am{sup 242m} concentration.The uncertainties in the inventory of fission products reach up to 30%.
We present a new code system which combines the Monte Carlo neutron transport code MCNP-4C and the inventory code ACAB as a suitable tool for high burnup calculations. Our main goal is to show that the system, by means of ACAB capabilities, enables us to assess the impact of neutron cross section uncertainties on the inventory and other inventory-related responses in high burnup applications. The potential impact of nuclear data uncertainties on some response parameters may be large, but only very few codes exist which can treat this effect. In fact, some of the most reported effective code systems in dealing with high burnup problems, such as CASMO-4, MCODE and MONTEBURNS, lack this capability. As first step, the potential of our system, ruling out the uncertainty capability, has been compared with that of those code systems, using a well referenced high burnup pin-cell benchmark exercise. It is proved that the inclusion of ACAB in the system allows to obtain results at least as reliable as those obtained using other inventory codes, such as ORIGEN2. Later on, the uncertainty analysis methodology implemented in ACAB, including both the sensitivity-uncertainty method and the uncertainty analysis by the Monte Carlo technique, is applied to this benchmark problem. We estimate the errors due to activation cross section uncertainties in the prediction of the isotopic content up to the high-burnup spent fuel regime. The most relevant uncertainties are remarked, and some of the most contributing cross sections to those uncertainties are identified. For instance, the most critical reaction for Am242m is Am241(n,γ-m). At 100 MWd/kg, the cross-section uncertainty of this reaction induces an error of 6.63% on the Am242m concentration.The uncertainties in the inventory of fission products reach up to 30%
The Fork+ burnup measurement system: Design and first measurement campaign
Previous work with the original Fork detector showed that burnup as determined by reactor records could be accurately allocated to spent nuclear fuel assemblies. The original Fork detector, designed by Los Alamos National Laboratory, used an ion chamber to measure gross gamma count and a fission chamber to measure neutrons from an activation source, 244Cm. In its review of the draft Topical Report on Burnup Credit, the US Nuclear Regulatory Commission indicated it felt uncomfortable with a measurement system that depended on reactor records for calibration. The Fork+ system was developed at Sandia National Laboratories under the sponsorship of the Electric Power Research Institute with the aim of providing this independent measurement capability. The initial Fork+ prototype was used in a measurement campaign at the Maine Yankee reactor. The campaign confirmed the applicability of the sensor approach in the Fork+ system and the efficiency of the hand-portable Fork+ prototype in making fuel assembly measurements. It also indicated potential design modifications that will be necessary before the Fork+ can be used effectively on high-burnup spent fuel
Accelerator-driven transmutation reactor analysis code system (ATRAS)
Sasa, Toshinobu; Tsujimoto, Kazufumi; Takizuka, Takakazu; Takano, Hideki [Japan Atomic Energy Research Inst., Tokai, Ibaraki (Japan). Tokai Research Establishment
1999-03-01
JAERI is proceeding a design study of the hybrid type minor actinide transmutation system which mainly consist of an intense proton accelerator and a fast subcritical core. Neutronics and burnup characteristics of the accelerator-driven system is important from a view point of the maintenance of subcriticality and energy balance during the system operation. To determine those characteristics accurately, it is necessary to involve reactions at high-energy region, which are not treated on ordinary reactor analysis codes. The authors developed a code system named ATRAS to analyze the neutronics and burnup characteristics of accelerator-driven subcritical reactor systems. ATRAS has a function of burnup analysis taking account of the effect of spallation neutron source. ATRAS consists of a spallation analysis code, a neutron transport codes and a burnup analysis code. Utility programs for fuel exchange, pre-processing and post-processing are also incorporated. (author)
Accelerator-driven transmutation reactor analysis code system (ATRAS)
JAERI is proceeding a design study of the hybrid type minor actinide transmutation system which mainly consist of an intense proton accelerator and a fast subcritical core. Neutronics and burnup characteristics of the accelerator-driven system is important from a view point of the maintenance of subcriticality and energy balance during the system operation. To determine those characteristics accurately, it is necessary to involve reactions at high-energy region, which are not treated on ordinary reactor analysis codes. The authors developed a code system named ATRAS to analyze the neutronics and burnup characteristics of accelerator-driven subcritical reactor systems. ATRAS has a function of burnup analysis taking account of the effect of spallation neutron source. ATRAS consists of a spallation analysis code, a neutron transport codes and a burnup analysis code. Utility programs for fuel exchange, pre-processing and post-processing are also incorporated. (author)
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.
An economic evaluation of a storage system for casks with burnup credit
It is generally recognized that casks designed with burnup credit are more economical than those without burnup credit. To estimate how much more economical they are, we made conceptual designs of transport/storage casks with and without burnup credit for PWR and BWR fuels of various uranium enrichment. The casks were designed to contain the maximum number of fuel assemblies under the necessary weight and dimensional limitations as well as the criticality and shielding criteria. The results showed that approximately 8 % to 44 % more fuel assemblies could be contained in casks with burnup credit. We then evaluated the economy of cask storage system incorporating the cask designs obtained above both with and without burnup credit. The results showed that the cost of storing casks with burnup credit is approximately 7 % to 30 % less expensive than storing casks without burnup credit. (J.P.N.)
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...
The fuel burn-up is an important parameter needed to be monitored and determined during a reactor operation and fuel management. The fuel burn-up can be calculated using computer codes and experimentally measured. This work presents the theory and experimental method applied to determine the burn-up of the irradiated and 36% enriched VVR-M2 fuel type assemblies of Dalat reactor. The method is based on measurement of Cs-137 absolute specific activity using gamma spectrometer. Designed measuring system consists of a collimator tube, high purity Germanium detector (HPGe) and associated electronics modules and online computer data acquisition system. The obtained results of measurement are comparable with theoretically calculated results. (author)
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
Application of SCALE4.4 system for burnup credit criticality analysis of PWR spent fuel
An investigation on the application of burnup credit for a PWR spent fuel storage pool has been carried out with the use of the SCALE 4.4 computer code system consisting of SAS2H and CSAS6 modules in association with 44-group SCALE cross-section library. Prior to the application of the computer code system, a series of bench markings have been performed in comparison with available data. A benchmarking of the SAS2h module has been done for experimental concentration data of 54 PWR spent fuel and then correction factors with a 95% probability at a 95% confidence level have been determined on the basis of the calculated and measured concentrations of 38 nuclides. After that, the bias which might have resulted from the use of the CSAS6 module has been calculated for 46 criticality experimental data of UO2 fuel and MOX fuel assemblies. The calculation bias with one-sided tolerance limit factor (2.086) corresponding to a 95% probability at a 95% confidence level has consequently been obtained to be 0.00834. Burnup credit criticality analysis has been done for the PWR spent fuel storage pool by means of the benchmarked or validated code system. It is revealed that the minimum burnup for safe storage is 7560 MWd/tU in 5 wt% enriched fuel if both actinides and fission products in spent fuel are taken into account. However, the minimum value required seems to be 9,565 MWd/tU in the same enriched fuel provided that only the actinides are taken into consideration. (author)
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
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
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
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)
Full text: IAEA uses today FORK Detector in attended and unattended mode for the verification of spent fuels. This system uses a neutron fission chamber and an ionisation chamber to combine total neutron counting and total gamma counting. CANBERRA now proposes the new SMOPY system, which enhances performance as it combines a fission chamber and a CZT gamma spectrometer. This new measurement capability associated with a depletion code embedded in the interpretation software of the system allows a complete identification of the burn-up of any type of fuel (also MOX for example). A first prototype of the SMOPY system was developed in collaboration between AREVA NC CEA and CANBERRA for safeguards but also burn-up credit applications. This prototype has been already used by the IAEA. CANBERRA has now completed the industrialization of this system adding new functionalities. The system allows also axial scanning of the fuel assembly instead of a single point measurement. Two types of interpretation of the measurement have been developed. The first one requires the irradiation history to determine very precisely the burn-up of the fuel assembly, thus allow to verify the operator's declaration. The second method is less precise but doesn't require any data of the fuel to determine the cooling time, burn-up, and the fuel type (MOX or LEU). The new CANBERRA industrialized SMOPY system will allow new possibilities of IAEA verifications and will also permit to address new scenarios of IAEA safeguards activities. (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)
For all the physical components that comprise a nuclear system there is an uncertainty. Assessing the impact of uncertainties in the simulation of fissionable material systems is essential for a best estimate calculation that has been replacing the conservative model calculations as the computational power increases. The propagation of uncertainty in a simulation using a Monte Carlo code by sampling the input parameters is recent because of the huge computational effort required. In this work a sample space of MCNPX calculations was used to propagate the uncertainty. The sample size was optimized using the Wilks formula for a 95. percentile and a two-sided statistical tolerance interval of 95%. Uncertainties in input parameters of the reactor considered included geometry dimensions and densities. It was showed the capacity of the sampling-based method for burnup when the calculations sample size is optimized and many parameter uncertainties are investigated together, in the same input. Particularly it was shown that during the burnup, the variances when considering all the parameters uncertainties is equivalent to the sum of variances if the parameter uncertainties are sampled separately
Development of burnup analysis system for fast reactor (3) (Contract research)
Improvement of the prediction accuracy for neutronics property of fast reactor cores is one of the most important study domains in terms of both achievement of high economical plant efficiency based on reasonably advanced designs and increased reliability and safety margins. In the previous study, considerable improvement of the prediction accuracy in neutronics design has been achieved in the development of the unified constants library as a fruit of a series of critical experiments such as JUPITER in application of the reactor constant adjustments. For design of fast reactor cores, however, improvement of not only static properties but also burnup properties is very important. For such purposes, it is necessary to improve the prediction accuracy of burnup properties using actual burnup data of 'JOYO' and 'MONJU', experimental and prototype fast reactors. Recently, study on effective burnup method for minor actinides becomes important theme. However, there is a problem that analysis work tends to become inefficient for lack of functionality suitable for analysis of composition change due to burnup since the conventional analysis system is targeted to critical assembly systems. Therefore development of burnup analysis system for fast reactors with modularity and flexibility is being done that would contribute to actual core design work and improvement of the prediction accuracy. In the previous study on 'Development of Burnup Analysis System for Fast Reactors (2)' in FY2006, design and implementation of models for detailed geometry of assembly, fuel loading pattern and so on, accompanied with specification and implementation of input file handling to construct data models. In the present study, a prototype system has been implemented in which functionalities are embedded for calculation of macroscopic cross section, core calculation and burnup calculation applying the fruits of the study 'Development of a Framework for the Neutronics Analysis System for Next
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
The purpose of this Technical Committee Meeting was to explore the status of international activities related to the use of burnup credit for spent fuel applications. This was the second major meeting on the issues of burnup credit for spent fuel management systems held since the IAEA began to monitor the uses of burnup credit in spent fuel management systems in 1997. Burnup credit for wet and dry storage systems is needed in many Member States to allow for increased initial fuel enrichment, and to increase the storage capacity and thus to avoid the need for extensive modifications of the spent fuel management systems involved. This document contains 31 individual papers presented at the Meeting; each of the papers was indexed separately
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)
Development of burnup analysis system for fast reactors (2) (Contract research)
Improvement on prediction accuracy for neutronics property of fast reactor cores is one of the most important study domains in terms of both achievement of high economical plant efficiency based on reasonably advanced designs and increased reliability and safety margins. In former study, considerable improvement on prediction accuracy in neutronics design has been achieved in the development of the unified constants library as a fruit of a series of critical experiments such as JUPITER in application of the reactor constant adjustments. For design of fast reactor cores, however, improvement of not only static properties but also burnup properties is very important. For such purpose, it is necessary to improve the prediction accuracy on burnup properties using actual burnup data of 'JOYO' and 'MONJU', experimental and prototype fast reactors. Recently, study on effective burnup method for minor actinides becomes an important theme. However, there is a problem that analysis work tends to become inefficient for lack of functionality suitable for analysis of composition change due to burnup since the conventional analysis system is targeted to critical assembly systems. Therefore development of burnup analysis system for fast reactors with modularity and flexibility can contribute actual core design work and improvement of prediction accuracy. In the previous study on 'Development of Burnup Analysis System (for Fast Reactors)' in FY2005, basic design was conducted to define each component in the system (input, solver, edit) and how to drive them. In this study, detailed design of the system and implementation of the I/O component were conducted according to the results in the basic design followed by proto-typing implementation. (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)
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
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)
Burnup verification at Arkansas Nuclear One-unit 1 using the Fork measurement system
The Fork measurement system, designed at Los Alamos National Laboratory for the International Atomic Energy Agency safeguards program, has been used for several years to examine spent fuel assemblies at nuclear reactors around the world. The objective of the test program described here is to demonstrate the ability of the Fork system to verify the records for assembly burnup at U.S. nuclear utilities. The measurements described here were performed at Arkansas Nuclear One, operated by Energy Operations, Inc. The Fork system was used to examine 34 assemblies in the storage pool of Arkansas Nuclear One-Unit 1. The correlation between the neutron measurements and the reactor records produced an average random deviation in burnup of 3.0% from the calibration, which translates into an average variation of 2.2% in the reactor records for burnup. The system proved to be compatible with utility operations
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)
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)
Sensitivity and uncertainty analysis of burnup reactivity for an accelerator-driven system
A burnup calculation is carried out for an accelerator-driven system (ADS) proposed by the Japan Atomic Energy Agency (JAEA) with the fourth version of JENDL, JENDL-4.0 and the previous one, JENDL-3.3. Considerable differences are seen in burnup reactivity between the nuclear data libraries for an initial phase (first burnup cycle) and an equilibrium phase (tenth burnup cycle). The differences in these values are investigated using two methods: a method by replacing a nuclear data library by nuclide and a sensitivity analysis technique. Among many contributors to them for the both phases, we identify major ones; (1) the initial phase: fission cross section and fission neutron multiplicity of 238Pu, capture cross section of 241Am, and (2) the equilibrium phase: capture cross section of 244Cm and 241Am, and inelastic scattering cross section of 206,207Pb. The uncertainty analysis shows that uncertainties in the burnup reactivity deduced from the JENDL-4.0 covariance data are comparable in magnitude to the differences between the nuclear data libraries, and are dominated by nuclear data parameters of 238Pu. Finally, we show the necessity of uncertainty evaluation of the branching ratio of 241Am capture reaction. (author)
The burn-up code OREST has a spectrum code assigned to it, which determines the neutron spectrum in the actual fuel element mixture at the start and during burn-up and carries out the resonance treatment for the most important uranium and transuranic element isotopes. The reliability of the OREST system is shown for UO2 burn-up in PWR's. Post-irradiation analyses of five UO2 fuel elements of KWO with an initial enrichment of 3.13% by weight of U235 and a mean burn-up of 28.4 GWd/tV are used for comparison. The reliability of OREST information for UO2 fuel in PWR's is proved by the good agreement between experiment and calculation, also compared with KfK's results. (orig./HP)
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)
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 System of Computerized Analysis for Licensing at Atomic industry (SCALA) is a Russian analogue of the well-known SCALE system. For criticality evaluations the ABBN-93 system is used with TWODANT and with joined American KENO and Russian MMK Monte-Carlo code MMKKENO. Using the same cross sections and input models, all these codes give results that coincide within the statistical uncertainties (for Monte-Carlo codes). Validation of criticality calculations using SCALA was performed using data presented in the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Another task of the work was to test the burnup capability of SCALA system in complex geometry in compare with other codes. Benchmark models of VVER type reactor assemblies with UO2 and MOX fuel including the cases with burnable gadolinium absorbers were calculated. KENO-VI and MMK codes were used for power distribution calculations, ORIGEN code was used for the isotopic kinetics calculations. (authors)
BOLD/VENTURE-4, Reactor Analysis System with Sensitivity and Burnup
cross section variation or correlation on nuclide concentrations is provided, but a temperature dependence is coded. Steady state condition with continuous fueling is established by a global iterative scheme that applies the criticality search scheme in the neutronics and models fuel movement directly in the exposure code. Time-dependent sensitivity data applies the forward march, reverse importance approach. The codes do not process data from the user input data stream allowing flexible task assignment along selected calculational paths. Multigroup cross section data are produced locally using the PSR-0063/AMPX II or CCC-0450/SCALE-2 code systems to produce resonance shielding (NITAWL) and cell weighted (XSDRN) microscopic cross sections. Locally, each code is compiled and loaded, and only one version is maintained in a quality assurance state in load module form. An on-line catalog procedure, installed with system support, provides job control instructions with nominal default of space allocation to files. Executing the catalog procedure makes the driver memory resident. The first user input data line must be the control module name used for the run. VENTURE-PC: The VENTURE module applies the finite-difference diffusion or a simple P1 approximation. VENTURE uses an outer-inner iteration scheme with several different data handling methods. Over-relaxation is applied to the inner and outer iterations, and succeeding flux iterates may be accelerated with the Chebychev process. - The BURNER code (module EXPOSURE) uses a difference formulation based on average generation rates; or a matrix exponential formulation to approximate the solution of the coupled burnup differential equations; or an explicit solution for simply coupled nuclide chains. Space dependence is included by working with zone averaged fluxes
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)
SRAC95; general purpose neutronics code system
Okumura, Keisuke; Tsuchihashi, Keichiro [Japan Atomic Energy Research Inst., Tokai, Ibaraki (Japan). Tokai Research Establishment; Kaneko, Kunio
1996-03-01
SRAC is a general purpose neutronics code system applicable to core analyses of various types of reactors. Since the publication of JAERI-1302 for the revised SRAC in 1986, a number of additions and modifications have been made for nuclear data libraries and programs. Thus, the new version SRAC95 has been completed. The system consists of six kinds of nuclear data libraries(ENDF/B-IV, -V, -VI, JENDL-2, -3.1, -3.2), five modular codes integrated into SRAC95; collision probability calculation module (PIJ) for 16 types of lattice geometries, Sn transport calculation modules(ANISN, TWOTRAN), diffusion calculation modules(TUD, CITATION) and two optional codes for fuel assembly and core burn-up calculations(newly developed ASMBURN, revised COREBN). In this version, many new functions and data are implemented to support nuclear design studies of advanced reactors, especially for burn-up calculations. SRAC95 is available not only on conventional IBM-compatible computers but also on scalar or vector computers with the UNIX operating system. This report is the SRAC95 users manual which contains general description, contents of revisions, input data requirements, detail information on usage, sample input data and list of available libraries. (author).
SRAC95; general purpose neutronics code system
SRAC is a general purpose neutronics code system applicable to core analyses of various types of reactors. Since the publication of JAERI-1302 for the revised SRAC in 1986, a number of additions and modifications have been made for nuclear data libraries and programs. Thus, the new version SRAC95 has been completed. The system consists of six kinds of nuclear data libraries(ENDF/B-IV, -V, -VI, JENDL-2, -3.1, -3.2), five modular codes integrated into SRAC95; collision probability calculation module (PIJ) for 16 types of lattice geometries, Sn transport calculation modules(ANISN, TWOTRAN), diffusion calculation modules(TUD, CITATION) and two optional codes for fuel assembly and core burn-up calculations(newly developed ASMBURN, revised COREBN). In this version, many new functions and data are implemented to support nuclear design studies of advanced reactors, especially for burn-up calculations. SRAC95 is available not only on conventional IBM-compatible computers but also on scalar or vector computers with the UNIX operating system. This report is the SRAC95 users manual which contains general description, contents of revisions, input data requirements, detail information on usage, sample input data and list of available libraries. (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)
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)
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)
Development of the next generation reactor analysis code system, MARBLE
A next generation reactor analysis code system, MARBLE, has been developed. MARBLE is a successor of the fast reactor neutronics analysis code systems, JOINT-FR and SAGEP-FR (conventional systems), which were developed for so-called JUPITER standard analysis methods. MARBLE has the equivalent analysis capability to the conventional system because MARBLE can utilize sub-codes included in the conventional system without any change. On the other hand, burnup analysis functionality for power reactors is improved compared with the conventional system by introducing models on fuel exchange treatment and control rod operation and so on. In addition, MARBLE has newly developed solvers and some new features of burnup calculation by the Krylov sub-space method and nuclear design accuracy evaluation by the extended bias factor method. In the development of MARBLE, the object oriented technology was adopted from the view-point of improvement of the software quality such as flexibility, expansibility, facilitation of the verification by the modularization and assistance of co-development. And, software structure called the two-layer system consisting of scripting language and system development language was applied. As a result, MARBLE is not an independent analysis code system which simply receives input and returns output, but an assembly of components for building an analysis code system (i.e. framework). Furthermore, MARBLE provides some pre-built analysis code systems such as the fast reactor neutronics analysis code system. SCHEME, which corresponds to the conventional code and the fast reactor burnup analysis code system, ORPHEUS. (author)
The FEDC Tokamak Systems Code calculates tokamak performance, cost, and configuration as a function of plasma engineering parameters. This version of the code models experimental tokamaks. It does not currently consider tokamak configurations that generate electrical power or incorporate breeding blankets. The code has a modular (or subroutine) structure to allow independent modeling for each major tokamak component or system. A primary benefit of modularization is that a component module may be updated without disturbing the remainder of the systems code as long as the imput to or output from the module remains unchanged
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)
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)
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)
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
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)
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)
SRAC2006: A comprehensive neutronics calculation code system
The SRAC is a code system applicable to neutronics analysis of a variety of reactor types. Since the publication of the second version of the users manual (JAERI-1302) in 1986 for the SRAC system, a number of additions and modifications to the functions and the library data have been made to establish a comprehensive neutronics code system. The current system includes major neutron data libraries (JENDL-3.3, JENDL-3.2, ENDF/B-VII, ENDF/B-VI.8, JEFF-3.1, JEF-2.2, etc.), and integrates five elementary codes for neutron transport and diffusion calculation; PIJ based on the collision probability method applicable to 16 kind of lattice models, SN transport codes ANISN(1D) and TWOTRN(2D), diffusion codes TUD(1D) and CITATION(multi-D). The system also includes an auxiliary code COREBN for multi-dimensional core burn-up calculation. (author)
The advantages of mobile CdZnTe (CZT) detector for nuclear safeguard applications of spent fuel burnup inspection in assembly storage pond are compactness, low cost and ease of operations. In this work, a mobile detection system shield with tungsten alloy was designed and then performed on-site. Net count rate of the 662 keV line of 137Cs was produced linearly with burnup as experimental data simulations shows, in which the deviation from linearity is smaller than 9%. As a result, the feasibility of the method using CZT detector to monitor spent nuclear fuel assembly burnup in a fuel pond was validated. The results calculated with Monte Carlo procedure Geant4 can provide a theoretical guide for the further burnup measurement. (author)
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)
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
This study investigates some antecipated transients without scram for a pressurized water cooled reactor, using coupling of the containment CORAN code to the ALMOD code system, under severe random conditions. This coupling has the objective of including containment model as part of an unified code system. These severe conditions include failure of reactor scram, following a station black-out and emergency power initiation for the burn-up status at the beginning and end of the cycle. Furthermore, for the burn-up status at the end of the cycle, a failure in the closure of the pressurizer relief valve was also investigated. (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 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
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)
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)
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)
FORTRAN code-evaluation system
Capps, J. D.; Kleir, R.
1977-01-01
Automated code evaluation system can be used to detect coding errors and unsound coding practices in any ANSI FORTRAN IV source code before they can cause execution-time malfunctions. System concentrates on acceptable FORTRAN code features which are likely to produce undesirable results.
Coupling CFD code with system code and neutron kinetic code
Vyskocil, Ladislav, E-mail: Ladislav.Vyskocil@ujv.cz; Macek, Jiri
2014-11-15
Highlights: • Coupling interface between CFD code Fluent and system code Athlet was created. • Athlet code is internally coupled with neutron kinetic code Dyn3D. • Explicit coupling of overlapped computational domains was used. • A coupled system of Athlet/Dyn3D+Fluent codes was successfully tested on a real case. - Abstract: The aim of this work was to develop the coupling interface between CFD code Fluent and system code Athlet internally coupled with neutron kinetic code Dyn3D. The coupling interface is intended for simulation of complex transients such as Main Steam Line Break scenarios, which cannot be modeled separately first by system and neutron kinetic code and then by CFD code, because of the feedback between the codes. In the first part of this article, the coupling method is described. Explicit coupling of overlapped computational domains is used in this work. The second part of the article presents a demonstration simulation performed by the coupled system of Athlet/Dyn3D and Fluent. The “Opening a Steam Dump to the Atmosphere” test carried out at the Temelin NPP (VVER-1000) was simulated by the coupled system. In this simulation, the primary and secondary circuits were modeled by Athlet, mixing in downcomer and lower plenum was simulated by Fluent and heat generation in the core was calculated by Dyn3D. The results of the simulation with Athlet/Dyn3D+Fluent were compared with the experimental data and the results from a calculation performed with Athlet/Dyn3D without Fluent.
Coupling CFD code with system code and neutron kinetic code
Highlights: • Coupling interface between CFD code Fluent and system code Athlet was created. • Athlet code is internally coupled with neutron kinetic code Dyn3D. • Explicit coupling of overlapped computational domains was used. • A coupled system of Athlet/Dyn3D+Fluent codes was successfully tested on a real case. - Abstract: The aim of this work was to develop the coupling interface between CFD code Fluent and system code Athlet internally coupled with neutron kinetic code Dyn3D. The coupling interface is intended for simulation of complex transients such as Main Steam Line Break scenarios, which cannot be modeled separately first by system and neutron kinetic code and then by CFD code, because of the feedback between the codes. In the first part of this article, the coupling method is described. Explicit coupling of overlapped computational domains is used in this work. The second part of the article presents a demonstration simulation performed by the coupled system of Athlet/Dyn3D and Fluent. The “Opening a Steam Dump to the Atmosphere” test carried out at the Temelin NPP (VVER-1000) was simulated by the coupled system. In this simulation, the primary and secondary circuits were modeled by Athlet, mixing in downcomer and lower plenum was simulated by Fluent and heat generation in the core was calculated by Dyn3D. The results of the simulation with Athlet/Dyn3D+Fluent were compared with the experimental data and the results from a calculation performed with Athlet/Dyn3D without Fluent
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.
Accumulative Landings System Code Tables
National Oceanic and Atmospheric Administration, Department of Commerce — Code Tables Used In Landings System. These tables assign meanings to the codes that appear in the data tables. Code tables exist for species, gear, state, county,...
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)
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
Gao, Wen
2015-01-01
This comprehensive and accessible text/reference presents an overview of the state of the art in video coding technology. Specifically, the book introduces the tools of the AVS2 standard, describing how AVS2 can help to achieve a significant improvement in coding efficiency for future video networks and applications by incorporating smarter coding tools such as scene video coding. Topics and features: introduces the basic concepts in video coding, and presents a short history of video coding technology and standards; reviews the coding framework, main coding tools, and syntax structure of AV
A guide to the AUS modular neutronics code system
A general description is given of the AUS modular neutronics code system, which may be used for calculations of a very wide range of fission reactors, fusion blankets and other neutron applications. The present system has cross-section libraries derived from ENDF/B-IV and includes modules which provide for lattice calculations, one-dimensional transport calculations, and one, two, and three-dimensional diffusion calculations, burnup calculations and the flexible editing of results. Details of all system aspects of AUS are provided but the major individual modules are only outlined. Sufficient information is given to enable other modules to be added to the system
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
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 the VTT's reactor physics code system
At VTT Energy several international reactor physics codes and nuclear data libraries are used in a variety of applications. The codes and libraries are under constant development and every now and then new updated versions are released, which are taken in use as soon as they have been validated at VTT Energy. The primary aim of the validation is to ensure that the code works properly, and that it can be used correctly. Moreover, the applicability of the codes and libraries are studied in order to establish their advantages and weak points. The capability of generating program-specific nuclear data for different reactor physics codes starting from the same evaluated data is sometimes of great benefit. VTT Energy has acquired a nuclear data processing system based on the NJOY-94.105 and TRANSX-2.15 processing codes. The validity of the processing system has been demonstrated by generating pointwise (MCNP) and groupwise (ANISN) temperature-dependent cross section sets for the benchmark calculations of the Doppler coefficient of reactivity. At VTT Energy the KENO-VI three-dimensional Monte Carlo code is used in criticality safety analyses. The KENO-VI code and the 44GROUPNDF5 data library have been validated at VTT Energy against the ZR-6 and LR-0 critical experiments. Burnup Credit refers to the reduction in reactivity of burned nuclear fuel due to the change in composition during irradiation. VTT Energy has participated in the calculational VVER-440 burnup credit benchmark in order to validate criticality safety calculation tools. (orig.)
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)
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
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)
The effects of fuel cycle parameters on nuclear waste environmental impact are analyzed for an advanced system that includes a Fusion-Fission Hybrid reactor. The system aims at reduction of the long-term radiotoxicity of waste by transmuting highly radiotoxic transuranics. However, the radiological risk of the system is measured by annual doses to the public, which are controlled by reactor operations, fuel cycle processes, waste treatment processes, and design of geological repositories. In this study, the waste environmental impact for a fuel cycle with a Fusion-Fission transmutation is analyzed as a function of three different parameters: burn-up, recovery efficiency and waste form durability for two different geological repositories, one with low actinide solubility and the other with high solubility. It is found that burn-up and recovery efficiency effects on environmental impact strongly depend on repository conditions, while the most influential parameter is found to be the durability of the waste form. (author)
WWER expert system for fuel failure analysis using the RTOP-CA code
The computer expert system for fuel failure analysis of WWER during operation is presented. The diagnostics is based on the measurement of specific activity of reference nuclides in reactor primary coolant and application of a computer code for the data interpretation. The data analysis includes an evaluation of tramp uranium mass in reactor core, detection of failures by iodine and caesium spikes, evaluation of burnup of defective fuel. Evaluation of defective fuel burnup was carried out by applying the relation of caesium nuclides activity in spikes and relations of activities of gaseous fission products for steady state operational conditions. The method of burnup evaluation of defective fuel by use of fission gas activity is presented in details. The neural-network analysis is performed for determination of failed fuel rod number and defect size. Results of the expert system application are illustrated for several fuel campaigns on operating WWER NPPs. (authors)
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.
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
Verification of DeCART/CAPP code system for VHTR by HTTR core analysis
The DeCART/CAPP code system has been developed and verified against the numerical benchmark calculations for an HTTR. The reference calculations have been carried out by the Monte Carlo McCARD code in which a double heterogeneity model was used. Verification results show that the DeCART/CAPP code system gives less negative MTC and RTC than the McCARD code, and thus the DeCART code overestimates the multiplication factors at states with a high moderator and reflector temperature. However, the DeCART/CAPP code system predicts more negative FTC than McCARD code does. In the depletion calculation for the HTTR single cell and single block, the error of DeCART/CAPP code system increases with the burnup (authors)
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.
A nuclear data library production system for advanced lattice codes
A nuclear data library production system has been written to recover and format nuclear data required in the recently released advanced lattice code DRAGON Version-4. State-of-the-art resonance self-shielding calculations require information that goes beyond WIMS-D type models. Dilution-dependent cross sections are required for all resonant reactions and for more than 10 specific dilutions. Ultra-fine multigroup cross section data is also required in the resolved energy domain. Another important aspect of advanced lattice codes is the explicit treatment of most neutron induced reactions in the burnup calculation. We need to perform power normalization due to energy from various neutron induced and decay reactions. Even though the decay energy contributes very little relative to the neutron induced reactions, the information will be very useful for post irradiation behavior of fuel. All this information is collected using DRAGR, a new post-processing module in NJOY99, and formatted in a single direct access hierarchical database. Burnup data is also recovered and the short-life isotopic data is automatically lumped. Moreover, PyNjoy, an object-oriented script, was developed to automate the recovery of Endf/B evaluations, the building, and the management of the database. This system was developed under the Lesser General Public License and is openly available. (author)
Recently several code systems have been developed and utilized at Japan Atomic Energy Research Institute. It is not easy to construct and maintain a code system, but the fact is not well known in the researchers. For this reason, it will be very useful to publish informations about design concepts, characteristics, necessary computer features and amounts of invested manpower for the developments of some recent JAERI code systems. In this report, a general view of required manpower on unification of nuclear codes is discussed and four code systems, i.e., SRAC for thermal reactor analysis, TRITON for Tokamak MHD analysis, SPEEDI for emergent environmental dose prediction and RADHEAT for radiation shielding analysis are presented. They are described from aspects of (1) purpose and schedule of development, (2) outline of system, (3) results of benchmark tests, (4) utilized computer features, (5) invested manpower, and (6) desirable computer features. Finally common aspects of four code systems from viewpoint of necessary computer hardwares and softwares are discussed for future development of code systems. (author)
Improvement of JRR-4 core management code system
Izumo, H.; Watanabe, S.; Nagatomi, H.; Hori, N. [Department of Research Reactor, Tokai Research Establishment, Japan Atomic Energy Institute, Tokai, Ibaraki (Japan)
2000-10-01
In the modification of JRR-4, the fuel was changed from 93% high enrichment uranium aluminized fuel to 20% low enriched uranium silicide fuel in conformity with the framework of reduced enrichment program on JAERI research reactors. As changing of this, JRR-4 core management code system which estimates excess reactivity of core, fuel burn-up and so on, was improved too. It had been difficult for users to operate the former code system because its input-output form was text-form. But, in the new code system (COMMAS-JRR), users are able to operate the code system without using difficult text-form input. The estimation results of excess reactivity of JRR-4 LEU fuel core were showed very good agreements with the measured value. It is the strong points of this new code system to be operated simply by using the windows form pictures act on a personal workstation equip with the graphical-user-interface (GUI), and to estimate accurately the specific characteristics of the LEU core. (author)
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
ESCADRE and ICARE code systems
The French sever accident code development program is following two parallel approaches: the first one is dealing with ''integral codes'' which are designed for giving immediate engineer answers, the second one is following a more mechanistic way in order to have the capability of detailed analysis of experiments, in order to get a better understanding of the scaling problem and reach a better confidence in plant calculations. In the first approach a complete system has been developed and is being used for practical cases: this is the ESCADRE system. In the second approach, a set of codes dealing first with primary circuit is being developed: a mechanistic core degradation code, ICARE, has been issued and is being coupled with the advanced thermalhydraulic code CATHARE. Fission product codes have been also coupled to CATHARE. The ''integral'' ESCADRE system and the mechanistic ICARE and associated codes are described. Their main characteristics are reviewed and the status of their development and assessment given. Future studies are finally discussed. 36 refs, 4 figs, 1 tab
SPECTRAL AMPLITUDE CODING OCDMA SYSTEMS USING ENHANCED DOUBLE WEIGHT CODE
F.N. HASOON
2006-12-01
Full Text Available A new code structure for spectral amplitude coding optical code division multiple access systems based on double weight (DW code families is proposed. The DW has a fixed weight of two. Enhanced double-weight (EDW code is another variation of a DW code family that can has a variable weight greater than one. The EDW code possesses ideal cross-correlation properties and exists for every natural number n. A much better performance can be provided by using the EDW code compared to the existing code such as Hadamard and Modified Frequency-Hopping (MFH codes. It has been observed that theoretical analysis and simulation for EDW is much better performance compared to Hadamard and Modified Frequency-Hopping (MFH codes.
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
Barr, W.L.; Bathke, C.G.; Brooks, J.N.; Bulmer, R.H.; Busigin, A.; DuBois, P.F.; Fenstermacher, M.E.; Fink, J.; Finn, P.A.; Galambos, J.D.; Gohar, Y.; Gorker, G.E.; Haines, J.R.; Hassanein, A.M.; Hicks, D.R.; Ho, S.K.; Kalsi, S.S.; Kalyanam, K.M.; Kerns, J.A.; Lee, J.D.; Miller, J.R.; Miller, R.L.; Myall, J.O.; Peng, Y-K.M.; Perkins, L.J.; Spampinato, P.T.; Strickler, D.J.; Thomson, S.L.; Wagner, C.E.; Willms, R.S.; Reid, R.L. (ed.)
1988-04-01
A tokamak systems code capable of modeling experimental test reactors has been developed and is described in this document. The code, named TETRA (for Tokamak Engineering Test Reactor Analysis), consists of a series of modules, each describing a tokamak system or component, controlled by an optimizer/driver. This code development was a national effort in that the modules were contributed by members of the fusion community and integrated into a code by the Fusion Engineering Design Center. The code has been checked out on the Cray computers at the National Magnetic Fusion Energy Computing Center and has satisfactorily simulated the Tokamak Ignition/Burn Experimental Reactor II (TIBER) design. A feature of this code is the ability to perform optimization studies through the use of a numerical software package, which iterates prescribed variables to satisfy a set of prescribed equations or constraints. This code will be used to perform sensitivity studies for the proposed International Thermonuclear Experimental Reactor (ITER). 22 figs., 29 tabs.
A tokamak systems code capable of modeling experimental test reactors has been developed and is described in this document. The code, named TETRA (for Tokamak Engineering Test Reactor Analysis), consists of a series of modules, each describing a tokamak system or component, controlled by an optimizer/driver. This code development was a national effort in that the modules were contributed by members of the fusion community and integrated into a code by the Fusion Engineering Design Center. The code has been checked out on the Cray computers at the National Magnetic Fusion Energy Computing Center and has satisfactorily simulated the Tokamak Ignition/Burn Experimental Reactor II (TIBER) design. A feature of this code is the ability to perform optimization studies through the use of a numerical software package, which iterates prescribed variables to satisfy a set of prescribed equations or constraints. This code will be used to perform sensitivity studies for the proposed International Thermonuclear Experimental Reactor (ITER). 22 figs., 29 tabs
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.)
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.
The Carmen code (theory and user's manual) is described. This code for assembly and core calculations uses diffusion theory (Citation), with feedback in the cross sections by zone due to the effects of burnup, water density, fuel temperature, Xenon and Samarium. The burnup calculation of a full cycle is solved in only an execution of Carmen, and in a reduced computer time. (auth.)
Yun, Hyung Ju; Kim, Do Yeon; Park, Kwang Heon; Hong, Ser Gi [Dept. of Nuclear Engineering, Kyung Hee University, Seoul (Korea, Republic of)
2016-06-15
Nuclear criticality safety analyses (NCSAs) considering burnup credit were performed for the GBC-32 cask. The used nuclear fuel assemblies (UNFAs) discharged from Hanbit Nuclear Power Plant Unit 3 Cycle 6 were loaded into the cask. Their axial burnup distributions and average discharge burnups were evaluated using the DeCART and Multi-purpose Analyzer for Static and Transient Effects of Reactors (MASTER) codes, and NCSAs were performed using SCALE 6.1/STandardized Analysis of Reactivity for Burnup Credit using SCALE (STARBUCS) and Monte Carlo N-Particle transport code, version 6 (MCNP 6). The axial burnup distributions were determined for 20 UNFAs with various initial enrichments and burnups, which were applied to the criticality analysis for the cask system. The UNFAs for 20- and 30-year cooling times were assumed to be stored in the cask. The criticality analyses indicated that keff values for UNFAs with nonuniform axial burnup distributions were larger than those with a uniform distribution, that is, the end effects were positive but much smaller than those with the reference distribution. The axial burnup distributions for 20 UNFAs had shapes that were more symmetrical with a less steep gradient in the upper region than the reference ones of the United States Department of Energy. These differences in the axial burnup distributions resulted in a significant reduction in end effects compared with the reference.
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
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
System Design Description for the TMAD Code
This document serves as the System Design Description (SDD) for the TMAD Code System, which includes the TMAD code and the LIBMAKR code. The SDD provides a detailed description of the theory behind the code, and the implementation of that theory. It is essential for anyone who is attempting to review or modify the code or who otherwise needs to understand the internal workings of the code. In addition, this document includes, in Appendix A, the System Requirements Specification for the TMAD System
The new methodology developed in this work has the following purposes: a) to implement a burnup capability into the HAMMER-TECHNION/9/computer code by using the CINDER-2/10/computer code to perform the transmutation analysis for the actinides and fission products; b) to implement a reduced version of the CINDER-2 fission product chain structure to treat explicity nearly 99% of all original CINDER-2 fission product absorption in a typical PWR unit cell; c) to treat the effect of the fission product neutron absorption in an unit cell in a multigroup basis; d) to develop a tentative validation procedure for the ENOF/C-V stable and long-lived fission product nuclear data based on the available experimental data/11-14/. The analysis will be performed by using the reduce chain in the coupled system CINDER-2 to generate the time dependent effective four group cross sections for actinides and fission products and CINDER-2 to perform the complete transmutation analysis with its built-in chain structure. (author)
Further validation of the KARATE-440 code system
In the last years several projects aiming at introduction of new VVER-440 fuel types and resulting in more economic fuel cycles were initiated. Increased average enrichment, modification of the lattice pitch and fuels diameter, profiled enrichment, application of burnable absorber, modification of the absorber assembly coupler part could lead to higher burnup and maximum allowed reactor power. The above fuel modifications and the upgraded regimes requiring more accurate calculations have necessitated the further development and validation of the KARATE code system: application of new, more accurate nuclear data, corresponding renewal of the multigroup libraries and the parametrized few group constants. For the further validation the operational data of PAKS (Hungary) NPP and zero reactor measurements were used. The global calculations, where burnup dependent node-wise power distributions, critical boron concentrations, reactivity coefficients of the core etc. are determined, have been validated against operational data of PAKS NPP. Measured critical boron concentrations, reactivity coefficients, radial temperature and axial self power detector signal distributions were used for the comparison for 14 new cycles of each four units. The fine mesh diffusion calculations were benchmarked by using the measurements of ZR-ZR-6 critical facility lattices containing Gd burnable poison (Authors)
Su, Bingjing; Hawari, Ayman, I.
2004-03-30
Modular Pebble Bed Reactor (MPBR) is a high temperature gas-cooled nuclear power reactor currently under study as a next generation reactor system. In addition to its inherently safe design, a unique feature of this reactor is its multi-pass fuel circulation in which the fuel pebbles are randomly loaded and continuously cycled through the core until they reach their prescribed End-of-Life burn-up limit. Unlike the situation with a conventional light water reactor, depending solely on computational methods to perform in-core fuel management for MPBR will be highly inaccurate. An on-line measurement system is needed to accurately assess whether a given pebble has reached its End-of-Life burn-up limit and thereby provide an on-line, automated go/no-go decision on fuel disposition on a pebble-by-pebble basis. This project investigated approaches to analyzing fuel pebbles in real time using gamma spectroscopy and possibly using passive neutron counting of spontaneous fission neutrons to provide the speed, accuracy, and burn-up range required for burnup determination of MPBR. It involved all phases necessary to develop and construct a burn-up monitor, including a review of the design requirements of the system, identification of detection methodologies, modeling and development of potential designs, and finally, the construction and testing of an operational detector system. Based upon the research work performed in this project, the following conclusions are made. In terms of using gamma spectrometry, two possible approaches were identified for burnup assay. The first approach is based on the measurement of the absolute activity of Cs-137. However, due to spectral interference and the need for absolute calibration of the spectrometer, the uncertainty in burnup determination using this approach was found to range from {approx} {+-}40% at beginning of life to {approx} {+-}10% at the discharge burnup. An alternative approach is to use a relative burnup indicator. In this
Modular Pebble Bed Reactor (MPBR) is a high temperature gas-cooled nuclear power reactor currently under study as a next generation reactor system. In addition to its inherently safe design, a unique feature of this reactor is its multi-pass fuel circulation in which the fuel pebbles are randomly loaded and continuously cycled through the core until they reach their prescribed End-of-Life burn-up limit. Unlike the situation with a conventional light water reactor, depending solely on computational methods to perform in-core fuel management for MPBR will be highly inaccurate. An on-line measurement system is needed to accurately assess whether a given pebble has reached its End-of-Life burn-up limit and thereby provide an on-line, automated go/no-go decision on fuel disposition on a pebble-by-pebble basis. This project investigated approaches to analyzing fuel pebbles in real time using gamma spectroscopy and possibly using passive neutron counting of spontaneous fission neutrons to provide the speed, accuracy, and burn-up range required for burnup determination of MPBR. It involved all phases necessary to develop and construct a burn-up monitor, including a review of the design requirements of the system, identification of detection methodologies, modeling and development of potential designs, and finally, the construction and testing of an operational detector system. Based upon the research work performed in this project, the following conclusions are made. In terms of using gamma spectrometry, two possible approaches were identified for burnup assay. The first approach is based on the measurement of the absolute activity of Cs-137. However, due to spectral interference and the need for absolute calibration of the spectrometer, the uncertainty in burnup determination using this approach was found to range from ∼ ±40% at beginning of life to ∼ ±10% at the discharge burnup. An alternative approach is to use a relative burnup indicator. In this case, a self
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)
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%.
Modular ORIGEN-S for multi-physics code systems
The ORIGEN-S code in the SCALE 6.0 nuclear analysis code suite is a well-validated tool to calculate the time-dependent concentrations of nuclides due to isotopic depletion, decay, and transmutation for many systems in a wide range of time scales. Application areas include nuclear reactor and spent fuel storage analyses, burnup credit evaluations, decay heat calculations, and environmental assessments. Although simple to use within the SCALE 6.0 code system, especially with the ORIGEN-ARP graphical user interface, it is generally complex to use as a component within an externally developed code suite because of its tight coupling within the infrastructure of the larger SCALE 6.0 system. The ORIGEN2 code, which has been widely integrated within other simulation suites, is no longer maintained by Oak Ridge National Laboratory (ORNL), has obsolete data, and has a relatively small validation database. Therefore, a modular version of the SCALE/ORIGEN-S code was developed to simplify its integration with other software packages to allow multi-physics nuclear code systems to easily incorporate the well-validated isotopic depletion, decay, and transmutation capability to perform realistic nuclear reactor and fuel simulations. SCALE/ORIGEN-S was extensively restructured to develop a modular version that allows direct access to the matrix solvers embedded in the code. Problem initialization and the solver were segregated to provide a simple application program interface and fewer input/output operations for the multi-physics nuclear code systems. Furthermore, new interfaces were implemented to access and modify the ORIGEN-S input variables and nuclear cross-section data through external drivers. Three example drivers were implemented, in the C, C++, and Fortran 90 programming languages, to demonstrate the modular use of the new capability. This modular version of SCALE/ORIGEN-S has been embedded within several multi-physics software development projects at ORNL, including
Some antecipated transients without scram (ATWS) for a pressurized water cooled reactor, model KWU 1300 MWe, are studied using coupling of the containment code CORAN to the system model code ALMOD, under severe random conditions. This coupling has the objective of including containment model as part of a unified code system. These severe conditions include failure of reactor scram, following a station black-out and emergency power initiation for the burn-up status at the beginning and end of the cycle. Furthermore, for the burn-up status at the end of the cycle a failure in the closure of the pressurizer relief valve was also investigated. For the beginning of the cycle, the containment participates actively during the transient. It is noted that the effect of the burn-up in the fuel is to reduce the seriousness of these transients. On the other hand, the failure in the closure of the pressurized relief valve makes this transients more severe. Moreover, the containment safety or radiological public safety is not affected in any of the cases. (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
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).
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.
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
Computer access security code system
Collins, Earl R., Jr. (Inventor)
1990-01-01
A security code system for controlling access to computer and computer-controlled entry situations comprises a plurality of subsets of alpha-numeric characters disposed in random order in matrices of at least two dimensions forming theoretical rectangles, cubes, etc., such that when access is desired, at least one pair of previously unused character subsets not found in the same row or column of the matrix is chosen at random and transmitted by the computer. The proper response to gain access is transmittal of subsets which complete the rectangle, and/or a parallelepiped whose opposite corners were defined by first groups of code. Once used, subsets are not used again to absolutely defeat unauthorized access by eavesdropping, and the like.
A comparative study of MONTEBURNS and MCNPX 2.6.0 codes in ADS simulations
The possible use of the MONTEBURNS and MCNPX 2.6.0 codes in Accelerator-driven systems (ADSs) simulations for fuel evolution description is discussed. ADSs are investigated for fuel breeding and long-lived fission product transmutation so simulations of fuel evolution have a great relevance. The burnup/depletion capability is present in both studied codes. MONTEBURNS code links Monte Carlo N-Particle Transport Code (MCNP) to the radioactive decay burnup code ORIGEN2, whereas MCNPX depletion/ burnup capability is a linked process involving steady-state flux calculations by MCNPX and nuclide depletion calculations by CINDER90. A lead-cooled accelerator-driven system fueled with thorium was simulated and the results obtained using MONTEBURNS code and the results from MCNPX 2.6.0 code were compared. The system criticality and the variation of the actinide inventory during the burnup were evaluated and the results indicate a similar behavior between the results of each code. (author)
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)
Studies on future application of burnup credit in Hungary
This paper describes the present status of the fuel storage and the possible future applications of burnup credit in wet and dry storage systems in Hungary. It gives a survey of the activities planned in AERI concerning the burnup credit. Some part of these investigations dealing with the influence of the axial changing of the assembly burnup are given in more details. (author)
Development of Three-dimensional Reactor Analysis Code System for Accelerator-Driven System, ADS3D
To investigate an Accelerator-Driven System (ADS) with sub-criticality control mechanism such as control rods or burnable poison, the ADS3D code has been developed on MARBLE which is a next generation reactor analysis code system developed by JAEA. In the past neutronics calculation for the ADS, JAEA employed RZ calculation models to realize efficient investigations. However, it was very difficult to model sub-criticality control mechanisms in RZ calculation models. The ADS3D code system is able to calculate the transportation of protons and neutrons, the burn-up calculation and the fuel exchange in three-dimensional calculation models. It means this code system can treat ADS concepts with sub-criticality control mechanism and makes it possible to investigate a new concept of ADS. (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
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)
Benchmark calculation with MOSRA-SRAC for burnup of a BWR fuel assembly
The Japan Atomic Energy Agency has developed the Modular Reactor Analysis Code System MOSRA to improve the applicability of neutronic characteristics modeling. The cell calculation module MOSRA-SRAC is based on the collision probability method and is one of the core modules of the MOSRA system. To test the module on a real-world problem, it was combined with the benchmark program 'Burnup Credit Criticality Benchmark Phase IIIC.' In this program participants are requested to submit the neutronic characteristics of burnup calculations for a BWR fuel assembly containing fuel rods poisoned with gadolinium (Gd2O3), which is similar to the fuel assembly at TEPCO's Fukushima Daiichi Nuclear Power Station. Because of certain restrictions of the MOSRA-SRAC burnup calculations part of the geometry model was homogenized. In order to verify the validity of MOSRA-SRAC, including the effects of the homogenization, the calculated burnup dependent infinite multiplication factor and the nuclide compositions were compared with those obtained with the burnup calculation code MVP-BURN which had already been validated for many benchmark problems. As a result of the comparisons, the applicability of MOSRA-SRAC module for the BWR assembly has been verified. Furthermore, it can be shown that the effects of the homogenization are smaller than the effects due to the calculation method for both multiplication factor and compositions. (author)
The Electron-Gamma Shower (EGS) code system is a general purpose package for the Monte Carlo simulation of the coupled transport of electrons and photons in an arbitrary geometry for particles with energies above a few keV up to several hundred GeV (depending on the atomic numbers of the target materials). This report introduces a new, enhanced version called EGS5. In addition to explaining and documenting the various enhancements and changes to the previous version (EGS4), this document includes several introductory and advanced tutorials on the use of EGS5, and also contains the EGS5 User Manual. Our intention has been to make this document wholly self-contained so that the user need not refer to the original EGS4 manual (SLAC-265) in order to use the code. To this end, we have taken the liberty of incorporating into Chapter 2 of this report those portions of Chapter 2 of SLAC-265 which describe physics models of EGS4 retained by EGS5, thereby documenting all the physics contained in EGS5. (author)
DRAGON 3.05D, Reactor Cell Calculation System with Burnup
1 - Description of program or function: The computer code DRAGON contains a collection of models that can simulate the neutron behavior of a unit cell or a fuel assembly in a nuclear reactor. It includes all of the functions that characterize a lattice cell code, namely: the interpolation of microscopic cross sections supplied by means of standard libraries; resonance self-shielding calculations in multidimensional geometries; multigroup and multidimensional neutron flux calculations that can take into account neutron leakage; transport-transport or transport-diffusion equivalence calculations as well as editing of condensed and homogenized nuclear properties for reactor calculations; and finally isotopic depletion calculations. 2 - Methods: The code DRAGON contains a multigroup flux solver conceived that can use a various algorithms to solve the neutron transport equation for the spatial and angular distribution of the flux. Each of these algorithms is presented in the form of a one-group solution procedure where the contributions from other energy groups are considered as sources. The current release of DRAGON contains five such algorithms. The JPM option that solves the integral transport equation using the J+- method, (interface current method applied to homogeneous blocks); the SYBIL option that solves the integral transport equation using the collision probability method for simple one dimensional (1-D) or two dimensional (2-D) geometries and the interface current method for 2-D Cartesian or hexagonal assemblies; the EXCELL/NXT option to solve the integral transport equation using the collision probability method for more general 2-D geometries and for three dimensional (3-D) assemblies; the MOCC option to solve the transport equation using the method of cyclic characteristics in 2-D Cartesian, and finally the MCU option to solve the transport equation using the method of characteristics (non cyclic) for 3-D Cartesian geometries. The execution of DRAGON is
Expansion of the CHR bone code system
This report describes the coding system used in the Center for Human Radiobiology (CHR) to identify individual bones and portions of bones of a complete skeletal system. It includes illustrations of various bones and bone segments with their respective code numbers. Codes are also presented for bone groups and for nonbone materials
Summary of high burnup fuel issues and NRC`s plan of action
Meyer, R.O.
1997-01-01
For the past two years the Office of Nuclear Regulatory Research has concentrated mostly on the so-called reactivity-initiated accidents -- the RIAs -- in this session of the Water Reactor Safety Information Meeting, but this year there is a more varied agenda. RIAs are, of course, not the only events of interest for reactor safety that are affected by extended burnup operation. Their has now been enough time to consider a range of technical issues that arise at high burnup, and a list of such issues being addressed in their research program is given here. (1) High burnup capability of the steady-state code (FRAPCON) used for licensing audit calculations. (2) General capability (including high burnup) of the transient code (FRAPTRAN) used for special studies. (3) Adequacy at high burnup of fuel damage criteria used in regulation for reactivity accidents. (4) Adequacy at high burnup of models and fuel related criteria used in regulation for loss-of-coolant accidents (LOCAs). (5) Effect of high burnup on fuel system damage during normal operation, including control rod insertion problems. A distinction is made between technical issues, which may or may not have direct licensing impacts, and licensing issues. The RIAs became a licensing issue when the French test in CABRI showed that cladding failures could occur at fuel enthalpies much lower than a value currently used in licensing. Fuel assembly distortion became a licensing issue when control rod insertion was affected in some operating plants. In this presentation, these technical issues will be described and the NRC`s plan of action to address them will be discussed.
Tandem Mirror Reactor Systems Code (Version I)
A computer code was developed to model a Tandem Mirror Reactor. Ths is the first Tandem Mirror Reactor model to couple, in detail, the highly linked physics, magnetics, and neutronic analysis into a single code. This report describes the code architecture, provides a summary description of the modules comprising the code, and includes an example execution of the Tandem Mirror Reactor Systems Code. Results from this code for two sensitivity studies are also included. These studies are: (1) to determine the impact of center cell plasma radius, length, and ion temperature on reactor cost and performance at constant fusion power; and (2) to determine the impact of reactor power level on cost
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
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
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)
Actinide-only burnup credit methodology for PWR spent nuclear fuel
A conservative methodology is presented that would allow taking credit for burnup in the criticality safety analysis of spent nuclear fuel (SNF) packages. The method is based on the assumption that the isotopic concentration in the SNF and cross sections of each isotope for which credit is taken must be supported by validation experiments. The method allows credit for 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. 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 spent nuclear fuel 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 by use of UO2 and UO2/Puo2 critical experiments. 3. Establish conditions for the SNF (depletion analysis) and package (criticality analysis) which bounds keff. 4. Use the validated codes and bounding conditions to generate package loading criteria (burnup credit loading curves). 5. Verify by measurement that SNF assemblies meet the package loading criteria and confirm proper assembly selection prior to loading. (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)
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
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.
Burnup credit implementation in spent fuel management
The criticality safety analysis of spent fuel management systems has traditionally assumed that the fuel is fresh. This results in significant conservatism in the calculated value of the system's reactivity. The concept of allowing reactivity credit for spent fuel offers economic incentives. Burnup Credit (BUC) could reduce mass limitation during dissolution of highly enriched PWR assemblies at the La Hague reprocessing plant. Furthermore, accounting for burnup credit enables the operator to avoid the use of Gd soluble poison in the dissolver for MOX assemblies. Analyses performed by DOE and its contractors have indicated that using BUC to maximize spent fuel transportation cask capacities is a justifiable concept that would result in public risk benefits and cost savings while fully maintaining criticality safety margins. In order to allow for Fission Products and Actinides in Criticality-Safety analyses, an extensive BUC experimental programme has been developed in France in the framework of the CEA-COGEMA collaboration. 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. Independent measurement systems, e.g. gamma spectrum detection systems, are needed to perform a true independent measurement of assembly burnup, without reliance on reactor records, using the gamma emission signatures fission products (mainly Cesium isotopes). (author)
Reconstruction of pin burnup characteristics from nodal calculations in hexagonal geometry
A reconstruction method has been developed for recovering pin burnup characteristics from fuel cycle calculations performed in hexagonal-z geometry using the nodal diffusion option of the DIF3D/REBUS-3 code system. Intra-modal distributions of group fluxes, nuclide densities, power density, burnup, and fluence are efficiently computed using polynomial shapes constrained to satisfy nodal information. The accuracy of the method has been tested by performing several numerical benchmark calculations and by comparing predicted local burnups to values measured for experimental assemblies in EBR-11. The results indicate that the reconstruction methods are quite accurate, yielding maximum errors in power and nuclide densities that are less than 2% for driver assemblies and typically less than 5% for blanket assemblies. 14 refs., 2 figs., 5 tabs
Graphical abstract: - Highlights: • Continuous Energy Monte-Carlo burnup code. • Instabilities of depletion calculation in loosely coupled system. • Advanced step model for burnup calculations. • Xenon profile oscillation in thermal reactor. • Parametrical study of instabilities. - Abstract: In this paper we use the Continuous Energy Monte-Carlo tool to expose the problem of burnup instabilities occurring in 1D and 2D systems based on PWR geometry. The intensity of power profile oscillations is studied as a function of geometry properties and time step length. We compare two step models for depletion procedure: classic staircase step model and stochastic implicit Euler method, that belongs to the family of predictor–corrector schemes. What is more, we consider the usage of better neutron source intensity value than beginning-of-step approximation. Required methodology was implemented into MCB5 simulation code. The practical conclusions about depletion calculations were formulated and the efficiency of advanced step model was confirmed
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)
AUS98 - The 1998 version of the AUS modular neutronic code system
Robinson, G.S.; Harrington, B.V
1998-07-01
AUS is a neutronics code system which may be used for calculations of a wide range of fission reactors, fusion blankets and other neutron applications. The present version, AUS98, has a nuclear cross section library based on ENDF/B-VI and includes modules which provide for reactor lattice calculations, one-dimensional transport calculations, multi-dimensional diffusion calculations, cell and whole reactor burnup calculations, and flexible editing of results. Calculations of multi-region resonance shielding, coupled neutron and photon transport, energy deposition, fission product inventory and neutron diffusion are combined within the one code system. The major changes from the previous AUS publications are the inclusion of a cross-section library based on ENDF/B-VI, the addition of the MICBURN module for controlling whole reactor burnup calculations, and changes to the system as a consequence of moving from IBM main-frame computers to UNIX workstations This report gives details of all system aspects of AUS and all modules except the POW3D multi-dimensional diffusion module refs., tabs.
AUS98 - The 1998 version of the AUS modular neutronic code system
AUS is a neutronics code system which may be used for calculations of a wide range of fission reactors, fusion blankets and other neutron applications. The present version, AUS98, has a nuclear cross section library based on ENDF/B-VI and includes modules which provide for reactor lattice calculations, one-dimensional transport calculations, multi-dimensional diffusion calculations, cell and whole reactor burnup calculations, and flexible editing of results. Calculations of multi-region resonance shielding, coupled neutron and photon transport, energy deposition, fission product inventory and neutron diffusion are combined within the one code system. The major changes from the previous AUS publications are the inclusion of a cross-section library based on ENDF/B-VI, the addition of the MICBURN module for controlling whole reactor burnup calculations, and changes to the system as a consequence of moving from IBM main-frame computers to UNIX workstations This report gives details of all system aspects of AUS and all modules except the POW3D multi-dimensional diffusion module
Neutronic and burnup studies of accelerator-driven systems dedicated to nuclear waste transmutation
Tucek, Kamil
2004-01-01
Partitioning and transmutation of plutonium, americium, and curium is inevitable if the radiotoxic inventory of spent nuclear fuel is to be reduced by more than a factor of 100. But, admixing minor actinides into the fuel severely degrades system safety parameters, particularly coolant void reactivity, Doppler effect, and (effective) delayed neutron fractions. The incineration process is therefore envisioned to be carried out in dedicated, accelerator-driven sub-critical reactors (ADS). Howev...
Study on the criticality safety evaluation method for burnup credit in JAERI
In relation to burnup credit, three tasks have been carried out at the Japan Atomic Energy Research Institute (JAERI) for establishing the evaluation method of criticality safety for a spent-fuel system, such as storage age ponds and transport casks. The first task is to prepare a benchmark database of criticality experiments and nuclide compositions of spent fuels. The database of nuclide composition is formed by data treasured at JAERI and data collected from the literature. For the database of criticality experiments, the effective multiplication factor of a spent-fuel assembly has been measured at JAERI and data collected from the literature. For the database of criticality experiments, the effective multiplication factor of a spent-fuel assembly has been measured at JAERI. The next task is to develop computer codes. The burnup and criticality codes have been developed and validated by analyzing a large number of benchmarks stored in the aforementioned database. The last task needed to establish the methodology in order to confirm the subcriticality of a spent-fuel system applying burnup credit is described. A reference fuel assembly is introduced so that the criticality of a system can be evaluated by using it, instead of modeling all fuel assemblies explicitly. To determine the nuclide composition of a spent fuel, a simple method is studied utilizing a large number of nuclide composition data stored in the database. Further, the effects of the axial burnup profile and calculation errors are discussed, and the remaining tasks are identified
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
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)
Hirayama, Hideo; Namito, Yoshihito; /KEK, Tsukuba; Bielajew, Alex F.; Wilderman, Scott J.; U., Michigan; Nelson, Walter R.; /SLAC
2005-12-20
In the nineteen years since EGS4 was released, it has been used in a wide variety of applications, particularly in medical physics, radiation measurement studies, and industrial development. Every new user and every new application bring new challenges for Monte Carlo code designers, and code refinements and bug fixes eventually result in a code that becomes difficult to maintain. Several of the code modifications represented significant advances in electron and photon transport physics, and required a more substantial invocation than code patching. Moreover, the arcane MORTRAN3[48] computer language of EGS4, was highest on the complaint list of the users of EGS4. The size of the EGS4 user base is difficult to measure, as there never existed a formal user registration process. However, some idea of the numbers may be gleaned from the number of EGS4 manuals that were produced and distributed at SLAC: almost three thousand. Consequently, the EGS5 project was undertaken. It was decided to employ the FORTRAN 77 compiler, yet include as much as possible, the structural beauty and power of MORTRAN3. This report consists of four chapters and several appendices. Chapter 1 is an introduction to EGS5 and to this report in general. We suggest that you read it. Chapter 2 is a major update of similar chapters in the old EGS4 report[126] (SLAC-265) and the old EGS3 report[61] (SLAC-210), in which all the details of the old physics (i.e., models which were carried over from EGS4) and the new physics are gathered together. The descriptions of the new physics are extensive, and not for the faint of heart. Detailed knowledge of the contents of Chapter 2 is not essential in order to use EGS, but sophisticated users should be aware of its contents. In particular, details of the restrictions on the range of applicability of EGS are dispersed throughout the chapter. First-time users of EGS should skip Chapter 2 and come back to it later if necessary. With the release of the EGS4 version
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)
V.S.O.P. (99/05) computer code system
V.S.O.P. is a computer code system for the comprehensive numerical simulation of the physics of thermal reactors. It implies the setup of the reactor and of the fuel element, processing of cross sections, neutron spectrum evaluation, neutron diffusion calculation in two or three dimensions, fuel burnup, fuel shuffling, reactor control, thermal hydraulics and fuel cycle costs. The thermal hydraulics part (steady state and time-dependent) is restricted to HTRs and to two spatial dimensions. The code can simulate the reactor operation from the initial core towards the equilibrium core. V.S.O.P.(99 / 05) represents the further development of V.S.O.P. (99). Compared to its precursor, the code system has been improved in many details. Major improvements and extensions have been included concerning the neutron spectrum calculation, the 3-d neutron diffusion options, and the thermal hydraulic section with respect to 'multi-pass'-fuelled pebblebed cores. This latest code version was developed and tested under the WINDOWS-XP - operating system. The storage requirement for the executables and the basic libraries associated with the code amounts to about 15 MB. Another 5 MB are required - if desired - for storage of the source code (∼65000 Fortran statements). (orig.)
V.S.O.P. (99/05) computer code system
Ruetten, H.J.; Haas, K.A.; Brockmann, H.; Scherer, W.
2005-11-01
V.S.O.P. is a computer code system for the comprehensive numerical simulation of the physics of thermal reactors. It implies the setup of the reactor and of the fuel element, processing of cross sections, neutron spectrum evaluation, neutron diffusion calculation in two or three dimensions, fuel burnup, fuel shuffling, reactor control, thermal hydraulics and fuel cycle costs. The thermal hydraulics part (steady state and time-dependent) is restricted to HTRs and to two spatial dimensions. The code can simulate the reactor operation from the initial core towards the equilibrium core. V.S.O.P.(99 / 05) represents the further development of V.S.O.P. (99). Compared to its precursor, the code system has been improved in many details. Major improvements and extensions have been included concerning the neutron spectrum calculation, the 3-d neutron diffusion options, and the thermal hydraulic section with respect to 'multi-pass'-fuelled pebblebed cores. This latest code version was developed and tested under the WINDOWS-XP - operating system. The storage requirement for the executables and the basic libraries associated with the code amounts to about 15 MB. Another 5 MB are required - if desired - for storage of the source code ({approx}65000 Fortran statements). (orig.)
ETF system code: composition and applications
A computer code has been developed for application to ETF tokamak system and conceptual design studies. The code determines cost, performance, configuration, and technology requirements as a function of tokamak parameters. The ETF code is structured in a modular fashion in order to allow independent modeling of each major tokamak component. The primary benefit of modularization is that it allows updating of a component module, such as the TF coil module, without disturbing the remainder of the system code as long as the input/output to the modules remains unchanged. The modules may be run independently to perform specific design studies, such as determining the effect of allowable strain on TF coil structural requirements, or the modules may be executed together as a system to determine global effects, such as defining the impact of aspect ratio on the entire tokamak system
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)
A combined 1D/3D fuel burnup analysis of generation IV light water reactor IRIS
A combined 1D/3D methodology for the fuel burnup analysis of generation IV light water reactors with thin boron coating that covers the fuel rods is described in this paper. This methodology is founded on three approximations. The first approximation assumes that the problem of fuel depletion in the entire 3D core can be resolved into two independent problems. One is a 3D Monte Carlo evolution of power distribution in large volumes (nodes) with the KENO-V.a code, and the other is a transport method evolution of burnup dependent fuel composition in 1D Wigner-Seitz cell for each node independently. With the second approximation, the time-dependent fuel composition in the node (e.g., in the fuel assembly) is calculated by using a 1D fuel depletion analysis with the SAS2H control module from the SCALE-4.4a code system. The third approximation involves smearing the boron coating with the clad (by volume homogenization). The proposed SAS2H/KENO-V.a methodology is verified for the case of 2D x-y model of IRIS 15x15 fuel assembly (with a reflective boundary condition) by using two well benchmarked code systems. The first one is MOCUP, a coupled MCNP-4C and ORIGEN2.1 utility code, and the second is KENO-V.a/ORIGEN2.1 code system recently developed by authors of this paper. It has been found that the proposed SAS2H/KENO-V.a methodology gives a satisfactory accuracy for keff and nuclide composition. Finally, this methodology was applied for 3D burnup analysis of IRIS-1000 benchmark≠44 core. Detailed keff and power density evolution with burnup are reported. (author)
Key issues in nuclear fuel cycle concerning high burn-up strategy
In the present high burn-up strategy in Japan, the economic efficiency and reduction of the spent nuclear fuel have been in progress. On the other hand, in the further progress of the strategy, several issues may appear. The amount and activity of nuclides, heat generation, and radiation for a fuel pin in the typical 17x17 PWR assembly were calculated as functions of burn-up and cooling time, using the SWAT code system. Waste loading in glass waste forms from spent UO2 fuel and MOX fuel were discussed, assuming the number of glass canisters of 150 liter per THM is 1.25 at 45 GWd/THM. The number of glass canisters per GWd is almost constant in the range of burn-up up to 70 GWd/THM. The amount of molybdate from Pu-239 fissions linearly increases as a function of burn-up similarly like increase from U-235 fissions. The current vitrification technology may not face serious situation to be required substantial reduction in waste loading relating to molybdate up to 70 GWd/THM. The initial cooling period prior to vitrification, the waste loading and the interim storage period prior to final disposal are major factors which determine the way of storage and final disposal. The higher burn-up above 45 GWd/THM may require pretreatment of HLLW or substantial reduction in waste loading to retain the integrity of the ceramic melter for e.g. five years. Further promotion of high burn-up strategy should be consistent with nuclear fuel cycle including waste management. A potential approach, a conceptual new reprocessing system for thermal reactors is described. (author)
Code Formal Verification of Operation System
Yu Zhang
2010-12-01
Full Text Available with the increasing pressure on non-function attributes (security, safety and reliability requirements of an operation system, high–confidence operation system is becoming more important. Formal verification is the only known way to guarantee that a system is free of programming errors. We research on formal verification of operation system kernel in system code level and take theorem proving and model checking as the main technical methods to resolve the key techniques of verifying operation system kernel in C code level. We present a case study to the verification of real-world C systems code derived from an implementation of μC/OS – II in the end.
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
Comparisons of the predicted and measured isotopic composition for high burnup PWR spent fuels
Comparisons between the calculated and measured isotopic composition for high burnup Korean PWR spent fuel samples were carried out. Spent fuel samples used in this study were obtained from commercial Korean PWRs, Ulchin unit 2 and Yonggwang unit 1. A radiochemical analysis of the spent fuel samples was performed to determine the isotopic compositions of U, Pu, and Nd. The depletion calculations which were carried out using the SAS2H control module in Version 5.1 of the SCALE code system were compared with the results of the radiochemical analyses. The results derived from the measured and calculated concentrations for each isotope of the corresponding samples were generally consistent with the earlier studies and the results were different within a few percent. The validity of the SAS2H control module in Version 5.1 of the SCALE code system could be confirmed in a high burnup spent fuel above 45 GWd/MTU
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)
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)
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)