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)
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)
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)
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)
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)
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 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)
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.
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
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)
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)
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.
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)
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.
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)
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)
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)
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)
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)
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
Development and verification of Monte Carlo burnup calculation system
Monte Carlo burnup calculation code system has been developed to evaluate accurate various quantities required in the backend field. From the Actinide Research in a Nuclear Element (ARIANE) program, by using, the measured nuclide compositions of fuel rods in the fuel assemblies irradiated in the commercial Netherlands BWR, the analyses have been performed for the code system verification. The code system developed in this paper has been verified through analysis for MOX and UO2 fuel rods. This system enables to reduce large margin assumed in the present criticality analysis for LWR spent fuels. (J.P.N.)
Burnup 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)
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.
Integrated burnup calculation code system SWAT
SWAT is an integrated burnup code system developed for analysis of post irradiation examination, transmutation of radioactive waste, and burnup credit problem. It enables us to analyze the burnup problem using neutron spectrum depending on environment of irradiation, combining SRAC which is Japanese standard thermal reactor analysis code system and ORIGEN2 which is burnup code widely used all over the world. SWAT makes effective cross section library based on results by SRAC, and performs the burnup analysis with ORIGEN2 using that library. SRAC and ORIGEN2 can be called as external module. SWAT has original cross section library on based JENDL-3.2 and libraries of fission yield and decay data prepared from JNDC FP Library second version. Using these libraries, user can use latest data in the calculation of SWAT besides the effective cross section prepared by SRAC. Also, User can make original ORIGEN2 library using the output file of SWAT. This report presents concept and user's manual of SWAT. (author)
Systemization of burnup sensitivity analysis code. 2
Towards the practical use of fast reactors, it is a very important subject to improve prediction accuracy for neutronic properties in LMFBR cores from the viewpoint of improvements on plant efficiency with rationally high performance cores and that on reliability and safety margins. A distinct improvement on accuracy in nuclear core design has been accomplished by the development of adjusted nuclear library using the cross-section adjustment method, in which the results of criticality experiments of JUPITER and so on are reflected. In the design of large LMFBR cores, however, it is important to accurately estimate not only neutronic characteristics, for example, reaction rate distribution and control rod worth but also burnup characteristics, for example, burnup reactivity loss, breeding ratio and so on. For this purpose, it is desired to improve prediction accuracy of burnup characteristics using the data widely obtained in actual core such as the experimental fast reactor 'JOYO'. The analysis of burnup characteristics is needed to effectively use burnup characteristics data in the actual cores based on the cross-section adjustment method. So far, a burnup sensitivity analysis code, SAGEP-BURN, has been developed and confirmed its effectiveness. However, there is a problem that analysis sequence become inefficient because of a big burden to users due to complexity of the theory of burnup sensitivity and limitation of the system. It is also desired to rearrange the system for future revision since it is becoming difficult to implement new functions in the existing large system. It is not sufficient to unify each computational component for the following reasons; the computational sequence may be changed for each item being analyzed or for purpose such as interpretation of physical meaning. Therefore, it is needed to systemize the current code for burnup sensitivity analysis with component blocks of functionality that can be divided or constructed on occasion. For
Systemization of burnup sensitivity analysis code
To practical use of fact reactors, it is a very important subject to improve prediction accuracy for neutronic properties in LMFBR cores from the viewpoints of improvements on plant efficiency with rationally high performance cores and that on reliability and safety margins. A distinct improvement on accuracy in nuclear core design has been accomplished by development of adjusted nuclear library using the cross-section adjustment method, in which the results of critical experiments of JUPITER and so on are reflected. In the design of large LMFBR cores, however, it is important to accurately estimate not only neutronic characteristics, for example, reaction rate distribution and control rod worth but also burnup characteristics, for example, burnup reactivity loss, breeding ratio and so on. For this purpose, it is desired to improve prediction accuracy of burnup characteristics using the data widely obtained in actual core such as the experimental fast reactor core 'JOYO'. The analysis of burnup characteristics is needed to effectively use burnup characteristics data in the actual cores based on the cross-section adjustment method. So far, development of a analysis code for burnup sensitivity, SAGEP-BURN, has been done and confirmed its effectiveness. However, there is a problem that analysis sequence become inefficient because of a big burden to user due to complexity of the theory of burnup sensitivity and limitation of the system. It is also desired to rearrange the system for future revision since it is becoming difficult to implement new functionalities in the existing large system. It is not sufficient to unify each computational component for some reasons; computational sequence may be changed for each item being analyzed or for purpose such as interpretation of physical meaning. Therefore it is needed to systemize the current code for burnup sensitivity analysis with component blocks of functionality that can be divided or constructed on occasion. For this
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.)
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
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
Propagation of statistical and nuclear data uncertainties in Monte Carlo burn-up calculations
Garcia-Herranz, Nuria [Departamento de Ingenieria Nuclear, Universidad Politecnica de Madrid, UPM (Spain)], E-mail: nuria@din.upm.es; Cabellos, Oscar [Departamento de Ingenieria Nuclear, Universidad Politecnica de Madrid, UPM (Spain); Sanz, Javier [Departamento de Ingenieria Energetica, Universidad Nacional de Educacion a Distancia, UNED (Spain); Juan, Jesus [Laboratorio de Estadistica, Universidad Politecnica de Madrid, UPM (Spain); Kuijper, Jim C. [NRG - Fuels, Actinides and Isotopes Group, Petten (Netherlands)
2008-04-15
Two methodologies to propagate the uncertainties on the nuclide inventory in combined Monte Carlo-spectrum and burn-up calculations are presented, based on sensitivity/uncertainty and random sampling techniques (uncertainty Monte Carlo method). Both enable the assessment of the impact of uncertainties in the nuclear data as well as uncertainties due to the statistical nature of the Monte Carlo neutron transport calculation. The methodologies are implemented in our MCNP-ACAB system, which combines the neutron transport code MCNP-4C and the inventory code ACAB. A high burn-up benchmark problem is used to test the MCNP-ACAB performance in inventory predictions, with no uncertainties. A good agreement is found with the results of other participants. This benchmark problem is also used to assess the impact of nuclear data uncertainties and statistical flux errors in high burn-up applications. A detailed calculation is performed to evaluate the effect of cross-section uncertainties in the inventory prediction, taking into account the temporal evolution of the neutron flux level and spectrum. Very large uncertainties are found at the unusually high burn-up of this exercise (800 MWd/kgHM). To compare the impact of the statistical errors in the calculated flux with respect to the cross uncertainties, a simplified problem is considered, taking a constant neutron flux level and spectrum. It is shown that, provided that the flux statistical deviations in the Monte Carlo transport calculation do not exceed a given value, the effect of the flux errors in the calculated isotopic inventory are negligible (even at very high burn-up) compared to the effect of the large cross-section uncertainties available at present in the data files.
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)
Development of an MCNP-tally based burnup code and validation through PWR benchmark exercises
The aim of this study is to evaluate the capabilities of a newly developed burnup code called BUCAL1. The code provides the full capabilities of the Monte Carlo code MCNP5, through the use of the MCNP tally information. BUCAL1 uses the fourth order Runge Kutta method with the predictor-corrector approach as the integration method to determine the fuel composition at a desired burnup step. Validation of BUCAL1 was done by code vs. code comparison. Results of two different kinds of codes are employed. The first one is CASMO-4, a deterministic multi-group two-dimensional transport code. The second kind is MCODE and MOCUP, a link MCNP-ORIGEN codes. These codes use different burnup algorithms to solve the depletion equations system. Eigenvalue and isotope concentrations were compared for two PWR uranium and thorium benchmark exercises at cold (300 K) and hot (900 K) conditions, respectively. The eigenvalue comparison between BUCAL1 and the aforementioned two kinds of codes shows a good prediction of the systems'k-inf values during the entire burnup history, and the maximum difference is within 2%. The differences between the BUCAL1 isotope concentrations and the predictions of CASMO-4, MCODE and MOCUP are generally better, and only for a few sets of isotopes these differences exceed 10%.
Development and validation of burnup function in reactor Monte Carlo 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 middle-of-step approximation and predictor-corrector method, are adopted by RMC to assure accuracy under large step size. RMC employs a spectrum-based method of tallying one-group cross section, which can considerably save computational time with negligible accuracy loss. According to validation results of benchmarks and examples, it is proved that the burnup function of RMC performs quite well in accuracy and efficiency. (author)
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
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%.
SERPENT Monte Carlo reactor physics code
SERPENT is a three-dimensional continuous-energy Monte Carlo reactor physics burnup calculation code, developed at VTT Technical Research Centre of Finland since 2004. The code is specialized in lattice physics applications, but the universe-based geometry description allows transport simulation to be carried out in complicated three-dimensional geometries as well. The suggested applications of SERPENT include generation of homogenized multi-group constants for deterministic reactor simulator calculations, fuel cycle studies involving detailed assembly-level burnup calculations, validation of deterministic lattice transport codes, research reactor applications, educational purposes and demonstration of reactor physics phenomena. The Serpent code has been publicly distributed by the OECD/NEA Data Bank since May 2009 and RSICC in the U. S. since March 2010. The code is being used in some 35 organizations in 20 countries around the world. This paper presents an overview of the methods and capabilities of the Serpent code, with examples in the modelling of WWER-440 reactor physics. (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
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
External linking scripts between Monte Carlo transport codes and burnup codes, and complete integration of burnup capability into Monte Carlo transport codes, have been or are currently being developed. Monte Carlo linked burnup methodologies may serve as an excellent benchmark for new deterministic burnup codes used for advanced systems; however, there are some instances where deterministic methodologies break down (i.e., heavily angularly biased systems containing exotic materials without proper group structure) and Monte Carlo burn up may serve as an actual design tool. Therefore, researchers are also developing these capabilities in order to examine complex, three-dimensional exotic material systems that do not contain benchmark data. Providing a reference scheme implies being able to associate statistical errors to any neutronic value of interest like k(eff), reaction rates, fluxes, etc. Usually in Monte Carlo, standard deviations are associated with a particular value by performing different independent and identical simulations (also referred to as 'cycles', 'batches', or 'replicas'), but this is only valid if the calculation itself is not biased. And, as will be shown in this paper, there is a bias in the methodology that consists of coupling transport and depletion codes because Bateman equations are not linear functions of the fluxes or of the reaction rates (those quantities being always measured with an uncertainty). Therefore, we have to quantify and correct this bias. This will be achieved by deriving an unbiased minimum variance estimator of a matrix exponential function of a normal mean. The result is then used to propose a reference scheme to solve Boltzmann/Bateman coupled equations, thanks to Monte Carlo transport codes. Numerical tests will be performed with an ad hoc Monte Carlo code on a very simple depletion case and will be compared to the theoretical results obtained with the reference scheme. Finally, the statistical error propagation
MCOR - Monte Carlo depletion code for reference LWR calculations
Research highlights: → Introduction of a reference Monte Carlo based depletion code with extended capabilities. → Verification and validation results for MCOR. → Utilization of MCOR for benchmarking deterministic lattice physics (spectral) codes. - Abstract: The MCOR (MCnp-kORigen) code system is a Monte Carlo based depletion system for reference fuel assembly and core calculations. The MCOR code is designed as an interfacing code that provides depletion capability to the LANL Monte Carlo code by coupling two codes: MCNP5 with the AREVA NP depletion code, KORIGEN. The physical quality of both codes is unchanged. The MCOR code system has been maintained and continuously enhanced since it was initially developed and validated. The verification of the coupling was made by evaluating the MCOR code against similar sophisticated code systems like MONTEBURNS, OCTOPUS and TRIPOLI-PEPIN. After its validation, the MCOR code has been further improved with important features. The MCOR code presents several valuable capabilities such as: (a) a predictor-corrector depletion algorithm, (b) utilization of KORIGEN as the depletion module, (c) individual depletion calculation of each burnup zone (no burnup zone grouping is required, which is particularly important for the modeling of gadolinium rings), and (d) on-line burnup cross-section generation by the Monte Carlo calculation for 88 isotopes and usage of the KORIGEN libraries for PWR and BWR typical spectra for the remaining isotopes. Besides the just mentioned capabilities, the MCOR code newest enhancements focus on the possibility of executing the MCNP5 calculation in sequential or parallel mode, a user-friendly automatic re-start capability, a modification of the burnup step size evaluation, and a post-processor and test-matrix, just to name the most important. The article describes the capabilities of the MCOR code system; from its design and development to its latest improvements and further ameliorations
MCOR - Monte Carlo depletion code for reference LWR calculations
Puente Espel, Federico, E-mail: fup104@psu.edu [Department of Mechanical and Nuclear Engineering, Pennsylvania State University (United States); Tippayakul, Chanatip, E-mail: cut110@psu.edu [Department of Mechanical and Nuclear Engineering, Pennsylvania State University (United States); Ivanov, Kostadin, E-mail: kni1@psu.edu [Department of Mechanical and Nuclear Engineering, Pennsylvania State University (United States); Misu, Stefan, E-mail: Stefan.Misu@areva.com [AREVA, AREVA NP GmbH, Erlangen (Germany)
2011-04-15
Research highlights: > Introduction of a reference Monte Carlo based depletion code with extended capabilities. > Verification and validation results for MCOR. > Utilization of MCOR for benchmarking deterministic lattice physics (spectral) codes. - Abstract: The MCOR (MCnp-kORigen) code system is a Monte Carlo based depletion system for reference fuel assembly and core calculations. The MCOR code is designed as an interfacing code that provides depletion capability to the LANL Monte Carlo code by coupling two codes: MCNP5 with the AREVA NP depletion code, KORIGEN. The physical quality of both codes is unchanged. The MCOR code system has been maintained and continuously enhanced since it was initially developed and validated. The verification of the coupling was made by evaluating the MCOR code against similar sophisticated code systems like MONTEBURNS, OCTOPUS and TRIPOLI-PEPIN. After its validation, the MCOR code has been further improved with important features. The MCOR code presents several valuable capabilities such as: (a) a predictor-corrector depletion algorithm, (b) utilization of KORIGEN as the depletion module, (c) individual depletion calculation of each burnup zone (no burnup zone grouping is required, which is particularly important for the modeling of gadolinium rings), and (d) on-line burnup cross-section generation by the Monte Carlo calculation for 88 isotopes and usage of the KORIGEN libraries for PWR and BWR typical spectra for the remaining isotopes. Besides the just mentioned capabilities, the MCOR code newest enhancements focus on the possibility of executing the MCNP5 calculation in sequential or parallel mode, a user-friendly automatic re-start capability, a modification of the burnup step size evaluation, and a post-processor and test-matrix, just to name the most important. The article describes the capabilities of the MCOR code system; from its design and development to its latest improvements and further ameliorations. Additionally
Usage of burnt fuel isotopic compositions from engineering codes in Monte-Carlo code calculations
Aleshin, Sergey S.; Gorodkov, Sergey S.; Shcherenko, Anna I. [Nuclear Research Centre ' ' Kurchatov Institute' ' , Moscow (Russian Federation)
2015-09-15
A burn-up calculation of VVER's cores by Monte-Carlo code is complex process and requires large computational costs. This fact makes Monte-Carlo codes usage complicated for project and operating calculations. Previously prepared isotopic compositions are proposed to use for the Monte-Carlo code (MCU) calculations of different states of VVER's core with burnt fuel. Isotopic compositions are proposed to calculate by an approximation method. The approximation method is based on usage of a spectral functionality and reference isotopic compositions, that are calculated by engineering codes (TVS-M, PERMAK-A). The multiplication factors and power distributions of FA and VVER with infinite height are calculated in this work by the Monte-Carlo code MCU using earlier prepared isotopic compositions. The MCU calculation data were compared with the data which were obtained by engineering codes.
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)
Coupled Monte Carlo and burnup calculations for research and power reactors
Reactor physics calculations require the solution of transport equation under consideration of complex energy structure of cross sections and spatial distribution of materials. Since the material composition changes during operation due to burn-up this changes have to be regarded in detail. An accurate method for reactor physics calculations is the continuous energy Monte Carlo method e.g. represented by the code MCNP(4C). The performance and quality of this code is demonstrated by own and world-wide applications. With the Monte Carlo method coupled with a program for the solution of the equations for nuclide build-up and decay not only the composition of fresh fuel but also the composition of irradiated fuel can be taken into account without external calculation of nuclide composition of irradiated materials. The coupled method represents at time the most accurate solution of criticality problems if the calculated parameters are determined with sufficient low statistical uncertainty. By this method important safety related parameters and reactivity conditions can be analysed for research reactors, standard LWR design as well as for innovative reactor or fuel designs or systems for actinide burning. (author)
Development of a MCNP–ORIGEN burn-up calculation code system and its accuracy assessment
Highlights: • MCNP and ORIGEN are coupled to perform nuclides depletion and decay calculation. • Coupled system MCORE uses “modified predictor corrector” approach. • MCORE can use different depletion schemes and simulate fuel shuffling. • MCORE is assessed by a “VVER-1000 LEU Assembly Computational Benchmark”. • MCORE is also assessed by a fast reactor benchmark problem. - Abstract: An MCNP–ORIGEN burn-up calculation code system, named MCORE (MCNP and ORIGEN burn-up Evaluation code), is developed in this work. MCORE makes use of the Monte Carlo neutron and photon transport code MCNP4C and nuclides depletion and decay calculation code ORIGEN2.1. MCNP and ORIGEN are coupled by data processing and linking subroutines. In MCORE, a so called “modified predictor corrector” approach is used. MCORE provides the capability of using different depletion calculation schemes and simulating fuel shuffling. Total nuclide density changes in active cells are considered in MCORE. The validity and applicability of the developed code are tested by investigating and predicting the neutronic and isotopic behavior of a “VVER-1000 LEU Assembly Computational Benchmark” at lattice level and a “Physics of Plutonium Recycling” fast reactor at core level (OECD-NEA). The comparison results show that the MCORE code predicts the nuclide composition within 5% accuracy and k∞ within 800 pcm at the end of the burn-up for LEU assembly (40 MWD/kg HM). For a fast reactor, the results obtained by MCORE are in the range of reported results except for 243Am. In general, MCORE results show a good agreement with the benchmark values
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
The Monte Carlo code MONK is a general program written to provide a high degree of flexibility to the user. MONK is distinguished by its detailed representation of nuclear data in point form i.e., the cross-section is tabulated at specific energies instead of the more usual group representation. The nuclear data are unadjusted in the point form but recently the code has been modified to accept adjusted group data as used in fast and thermal reactor applications. The various geometrical handling capabilities and importance sampling techniques are described. In addition to the nuclear data aspects, the following features are also described; geometrical handling routines, tracking cycles, neutron source and output facilities. 12 references. (U.S.)
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)
Dufek, Jan; Anglart, Henryk
2013-01-01
Numerically stable Monte Carlo burnup calculations of nuclear fuel cycles are now possible with the previously derived Stochastic Implicit Euler method based coupling scheme. In this paper, we show that this scheme can be easily extended to include the thermal-hydraulic feedback during the Monte Carlo burnup simulations, while preserving its unconditional stability property. At each time step, the implicit solution (for the end-of-step neutron flux, fuel nuclide densities and thermal-hydrauli...
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)
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)
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.
Monte Carlo studies on the burnup measurement for the high temperature gas cooling reactor
Online fuel pebble burnup measurement in a future high temperature gas cooling reactor is proposed for implementation through a high purity germanium (HPGe) gamma spectrometer. By using KORIGEN software and MCNP Monte Carlo simulations, the single pebble gamma radiations to be recorded in the detector are simulated under different irradiation histories. A specially developed algorithm is applied to analyze the generated spectra to reconstruct the gamma activity of the 137Cs monitoring nuclide. It is demonstrated that by taking into account the intense interfering peaks, the 137Cs activity in the spent pebbles can be derived with a standard deviation of 3.0% (1σ). The results support the feasibility of utilizing the HPGe spectrometry in the online determination of the pebble burnup in future modular pebble bed reactors. (authors)
Monte Carlo solver for UWB1 nuclear fuel depletion code
Highlights: • A new Monte Carlo solver was developed in order to speed-up depletion calculations. • For LWR model, UWB1 Monte Carlo solver is on average 10 times faster than MCNP6. • The UWB1 code will allow faster calculation analysis of BA parameters in fuel design. - Abstract: Recent nuclear reactor burnable absorber research tries to introduce new materials in the nuclear fuel. As a part of this effort, a fast computational tool is being developed for the advanced nuclear fuel. The first version of the newly developed UWB1 fast nuclear fuel depletion code significantly reduced calculation time by omitting the solution step for the Boltzmann transport equation. However, estimation of neutron multiplication factor during depletion was not sufficiently calculated. Therefore, at least one transport calculation for fuel depletion is necessary. This paper presents a new Monte Carlo solver that is implemented into the UWB1 code. The UWB1 Monte Carlo solver calculates neutron multiplication factor and neutron flux in the fuel for collapsed cross sections. Accuracy of the solver is supported by using current nuclear data stored in the ENDF/B-VII.1 library. Speed of the solver is the product of development focusing on minimization of CPU utilization at the expense of RAM demands. The UWB1 Monte Carlo solver is approximately 14 times faster than the MCNP6 reference code when one transport equation solution within fuel depletion is compared. Another speed-up can be achieved by employing advanced depletion scheme in the coupled transport and burnup equations. The resulting faster code will be used in optimization studies for ideal burnable absorber material selection where many various materials and concentrations will be evaluated
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)
Monte Carlo uncertainty propagation approaches in ADS burn-up calculations
Highlights: ► Two Monte Carlo uncertainty propagation approaches are compared. ► How to make both approaches equivalent is presented and applied. ► ADS burn-up calculation is selected as the application of approaches. ► The cross-section uncertainties of 239Pu and 241Pu are propagated. ► Cross-correlations appear as a source of differences between approaches. - Abstract: In activation calculations, there are several approaches to quantify uncertainties: deterministic by means of sensitivity analysis, and stochastic by means of Monte Carlo. Here, two different Monte Carlo approaches for nuclear data uncertainty are presented: the first one is the Total Monte Carlo (TMC). The second one is by means of a Monte Carlo sampling of the covariance information included in the nuclear data libraries to propagate these uncertainties throughout the activation calculations. This last approach is what we named Covariance Uncertainty Propagation, CUP. This work presents both approaches and their differences. Also, they are compared by means of an activation calculation, where the cross-section uncertainties of 239Pu and 241Pu are propagated in an ADS activation calculation
Revised Burnup Code System SWAT: Description and Validation Using Postirradiation Examination Data
The burnup code system Step-Wise Burnup Analysis Code System (SWAT) is revised for use in a burnup credit analysis. An important feature of the revised SWAT is that its functions are achieved by calling validated neutronics codes without any changes to the original codes. This feature is realized with a system function of the operating system, which allows the revised SWAT to be independent of the development status of each code.A package of the revised SWAT contains the latest libraries based on JENDL-3.2 and the second version of the JNDC FP library. These libraries allow us to analyze burnup problems, such as an analysis of postirradiation examination (PIE), using the latest evaluated data of not only cross sections but also fission yield and decay constants.Another function of the revised SWAT is a library generator for the ORIGEN2 code, which is one of the most reliable burnup codes. ORIGEN2 users can obtain almost the same results with the revised SWAT using the library prepared by this function.The validation of the revised SWAT is conducted by calculation of the Organization for Economic Cooperation and Development/Nuclear Energy Agency burnup credit criticality safety benchmark Phase I-B and analyses of PIE data for spent fuel from Takahama Unit 3. The analysis of PIE data shows that the revised SWAT can predict the isotopic composition of main uranium and plutonium with a deviation of 5% from experimental results taken from UO2 fuels of 17 x 17 fuel assemblies. Many results of fission products including samarium are within a deviation of 10%. This means that the revised SWAT has high reliability to predict the isotopic composition for pressurized water reactor spent fuel
Development and validation of a fast reactor core burnup code - FARCOB
Mohanakrishnan, P. [Indira Gandhi Centre for the Atomic Research, CDO, Reactor Physics Division, Kalpakkam, TN 603 102 (India)], E-mail: mohana@igcar.gov.in
2008-02-15
A large fast breeder reactor requires the accurate estimation of power produced in different parts of the reactor core and blanket during any operating condition for a safe and economic operation through out reactor life time. A fast reactor core simulation code FARCOB based on multigroup diffusion theory has been developed in IGCAR for core simulation of PFBR reactor under construction. FARCOB uses centre mesh differencing scheme with triangular meshes in the X-Y plane. Steady state solution results match exactly with those of other reputed codes DIF3D and VENTURE for SNR-300 benchmarks. For burnup simulation, core is divided into radial and axial burnup zones and burnup equations are solved at constant power. Burnable fuel and blanket number densities are found and stored for each mesh, so that the user can alter burnup zones and core geometry after a burnup step. For validation, results of FARCOB has been compared with results of other institutes in two burnup benchmarks (ANL 1000 MWe benchmark and BN-600 hybrid core benchmark). It is found that FARCOB results match well with those of the other institutes.
Development of a New Monte Carlo reactor physics code
Monte Carlo neutron transport codes are widely used in various reactor physics applications, traditionally related to criticality safety analyses, radiation shielding problems, detector modelling and validation of deterministic transport codes. The main advantage of the method is the capability to model geometry and interaction physics without major approximations. The disadvantage is that the modelling of complicated systems is very computing-intensive, which restricts the applications to some extent. The importance of Monte Carlo calculation is likely to increase in the future, along with the development in computer capacities and parallel calculation. An interesting near-future application for the Monte Carlo method is the generation of input parameters for deterministic reactor simulator codes. These codes are used in coupled LWR full-core analyses and typically based on few-group nodal diffusion methods. The input data consists of homogenised few-group constants, presently generated using deterministic lattice transport codes. The task is becoming increasingly challenging, along with the development in nuclear technology. Calculations involving high-burnup fuels, advanced MOX technology and next-generation reactor systems are likely to cause problems in the future, if code development cannot keep up with the applications. A potential solution is the use of Monte Carlo based lattice transport codes, which brings all the advantages of the calculation method. So far there has been only a handful of studies on group constant generation using the Monte Carlo method, although the interest has clearly increased during the past few years. The homogenisation of reaction cross sections is simple and straightforward, and it can be carried out using any Monte Carlo code. Some of the parameters, however, require the use of special techniques that are usually not available in general-purpose codes. The main problem is the calculation of neutron diffusion coefficients, which
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
Development of Monte Carlo-based pebble bed reactor fuel management code
Highlights: • A new Monte Carlo-based fuel management code for OTTO cycle pebble bed reactor was developed. • The double-heterogeneity was modeled using statistical method in MVP-BURN code. • The code can perform analysis of equilibrium and non-equilibrium phase. • Code-to-code comparisons for Once-Through-Then-Out case were investigated. • Ability of the code to accommodate the void cavity was confirmed. - Abstract: A fuel management code for pebble bed reactors (PBRs) based on the Monte Carlo method has been developed in this study. The code, named Monte Carlo burnup analysis code for PBR (MCPBR), enables a simulation of the Once-Through-Then-Out (OTTO) cycle of a PBR from the running-in phase to the equilibrium condition. In MCPBR, a burnup calculation based on a continuous-energy Monte Carlo code, MVP-BURN, is coupled with an additional utility code to be able to simulate the OTTO cycle of PBR. MCPBR has several advantages in modeling PBRs, namely its Monte Carlo neutron transport modeling, its capability of explicitly modeling the double heterogeneity of the PBR core, and its ability to model different axial fuel speeds in the PBR core. Analysis at the equilibrium condition of the simplified PBR was used as the validation test of MCPBR. The calculation results of the code were compared with the results of diffusion-based fuel management PBR codes, namely the VSOP and PEBBED codes. Using JENDL-4.0 nuclide library, MCPBR gave a 4.15% and 3.32% lower keff value compared to VSOP and PEBBED, respectively. While using JENDL-3.3, MCPBR gave a 2.22% and 3.11% higher keff value compared to VSOP and PEBBED, respectively. The ability of MCPBR to analyze neutron transport in the top void of the PBR core and its effects was also confirmed
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.
SPQR: a Monte Carlo reactor kinetics code
The SPQR Monte Carlo code has been developed to analyze fast reactor core accident problems where conventional methods are considered inadequate. The code is based on the adiabatic approximation of the quasi-static method. This initial version contains no automatic material motion or feedback. An existing Monte Carlo code is used to calculate the shape functions and the integral quantities needed in the kinetics module. Several sample problems have been devised and analyzed. Due to the large statistical uncertainty associated with the calculation of reactivity in accident simulations, the results, especially at later times, differ greatly from deterministic methods. It was also found that in large uncoupled systems, the Monte Carlo method has difficulty in handling asymmetric perturbations
The Serpent Monte Carlo Code: Status, Development and Applications in 2013
Leppänen, Jaakko; Pusa, Maria; Viitanen, Tuomas; Valtavirta, Ville; Kaltiaisenaho, Toni
2014-06-01
The Serpent Monte Carlo reactor physics burnup calculation code has been developed at VTT Technical Research Centre of Finland since 2004, and is currently used in 100 universities and research organizations around the world. This paper presents the brief history of the project, together with the currently available methods and capabilities and plans for future work. Typical user applications are introduced in the form of a summary review on Serpent-related publications over the past few years.
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)
Coded aperture optimization using Monte Carlo simulations
Coded apertures using Uniformly Redundant Arrays (URA) have been unsuccessfully evaluated for two-dimensional and three-dimensional imaging in Nuclear Medicine. The images reconstructed from coded projections contain artifacts and suffer from poor spatial resolution in the longitudinal direction. We introduce a Maximum-Likelihood Expectation-Maximization (MLEM) algorithm for three-dimensional coded aperture imaging which uses a projection matrix calculated by Monte Carlo simulations. The aim of the algorithm is to reduce artifacts and improve the three-dimensional spatial resolution in the reconstructed images. Firstly, we present the validation of GATE (Geant4 Application for Emission Tomography) for Monte Carlo simulations of a coded mask installed on a clinical gamma camera. The coded mask modelling was validated by comparison between experimental and simulated data in terms of energy spectra, sensitivity and spatial resolution. In the second part of the study, we use the validated model to calculate the projection matrix with Monte Carlo simulations. A three-dimensional thyroid phantom study was performed to compare the performance of the three-dimensional MLEM reconstruction with conventional correlation method. The results indicate that the artifacts are reduced and three-dimensional spatial resolution is improved with the Monte Carlo-based MLEM reconstruction.
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)
Successful vectorization - reactor physics Monte Carlo code
Most particle transport Monte Carlo codes in use today are based on the ''history-based'' algorithm, wherein one particle history at a time is simulated. Unfortunately, the ''history-based'' approach (present in all Monte Carlo codes until recent years) is inherently scalar and cannot be vectorized. In particular, the history-based algorithm cannot take advantage of vector architectures, which characterize the largest and fastest computers at the current time, vector supercomputers such as the Cray X/MP or IBM 3090/600. However, substantial progress has been made in recent years in developing and implementing a vectorized Monte Carlo algorithm. This algorithm follows portions of many particle histories at the same time and forms the basis for all successful vectorized Monte Carlo codes that are in use today. This paper describes the basic vectorized algorithm along with descriptions of several variations that have been developed by different researchers for specific applications. These applications have been mainly in the areas of neutron transport in nuclear reactor and shielding analysis and photon transport in fusion plasmas. The relative merits of the various approach schemes will be discussed and the present status of known vectorization efforts will be summarized along with available timing results, including results from the successful vectorization of 3-D general geometry, continuous energy Monte Carlo. (orig.)
Recent developments of the TRANSURANUS code with emphasis on high burnup phenomena
TRANSURANUS is a computer program for the thermal and mechanical analysis of fuel rods in nuclear reactors, which is developed at the Institute for Transuranium Elements. The code is in use in several European organisations, both in research and industry. In the paper the recent developments are summarised: the burnup degradation of the fuel thermal conductivity as well as the effects of gadolinium on the radial power distribution and thermal conductivity. Fission gas release from the High Burnup Structure is discussed. Finally, a new numerical method is outlined that is able to treat the highly non-linear mechanical equations in transients (RIAs and LOCAs). (author)
The burnup analysis program based on three dimensional discrete ordinates (SN) neutron/photon transport method has been developed by the FDS team, China, to aid in analysis, prediction, and optimization of fuel burnup performance in a nuclear reactor. The program uses output parameters generated by three-dimensional SN trans- port code to determine the isotopic inventory and anisotropic flux distribution as a function of time. For a fueled region, neutron transmutation, fuel depletion, fission-product poisoning, actinide generation, and burnable poison loading and depletion effects are included in the calculation. The IAEA benchmark test problem has been correctly calculated and analyzed to validate the system. (authors)
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.
Parallel processing Monte Carlo radiation transport codes
Issues related to distributed-memory multiprocessing as applied to Monte Carlo radiation transport are discussed. Measurements of communication overhead are presented for the radiation transport code MCNP which employs the communication software package PVM, and average efficiency curves are provided for a homogeneous virtual machine
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
Monte Carlo simulation code modernization
CERN. Geneva
2015-01-01
The continual development of sophisticated transport simulation algorithms allows increasingly accurate description of the effect of the passage of particles through matter. This modelling capability finds applications in a large spectrum of fields from medicine to astrophysics, and of course HEP. These new capabilities however come at the cost of a greater computational intensity of the new models, which has the effect of increasing the demands of computing resources. This is particularly true for HEP, where the demand for more simulation are driven by the need of both more accuracy and more precision, i.e. better models and more events. Usually HEP has relied on the "Moore's law" evolution, but since almost ten years the increase in clock speed has withered and computing capacity comes in the form of hardware architectures of many-core or accelerated processors. To harness these opportunities we need to adapt our code to concurrent programming models taking advantages of both SIMD and SIMT architectures. Th...
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)
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)
Morse Monte Carlo Radiation Transport Code System
Emmett, M.B.
1975-02-01
The report contains sections containing descriptions of the MORSE and PICTURE codes, input descriptions, sample problems, deviations of the physical equations and explanations of the various error messages. The MORSE code is a multipurpose neutron and gamma-ray transport Monte Carlo code. Time dependence for both shielding and criticality problems is provided. General three-dimensional geometry may be used with an albedo option available at any material surface. The PICTURE code provide aid in preparing correct input data for the combinatorial geometry package CG. It provides a printed view of arbitrary two-dimensional slices through the geometry. By inspecting these pictures one may determine if the geometry specified by the input cards is indeed the desired geometry. 23 refs. (WRF)
Monte Carlo calculations of the REBUS critical experiment for validation of burnup credit
The application of burnup credit (BUC) to criticality safety analysis for Spent Nuclear Fuel (SNF) configurations requires the implementation of both estimation of the SNF composition with the aid of depletion calculation tools and estimation of the SNF reactivity with the aid of criticality calculation tools. Amongst the several experimental programs dedicated to the validation of both calculation tools, REBUS is distinguished by a combination of chemical analysis and critical experiment. In addition to detailed assays of irradiated fuel, the reactivity worth of the fuel rods under investigation is measured both before and after irradiation. Since a whole bundle of fuel rods is used in the experiment, the change in reactivity is significant enough to be observable by Monte Carlo calculations. Thus, the calculation tools which see the most widespread use in SNF critical safety applications can be validated directly. Apart from the effective neutron multiplication factor keff, REBUS also provides measurements of the flux and fission rate distributions. While the program comprises investigation of commercial UO2 fuel rods and mixed oxide (MOX) fuel from a research reactor, the presentation will focus on the commercial UO2 fuel with an overview of the experimental setup and first results from the analysis. (author)
First steps towards a validation of the new burnup and depletion code TNT
In the frame of the fusion of the core design calculation capabilities, represented by V.S.O.P., and the accident calculation capabilities, represented by MGT(-3D), the successor of the TINTE code, difficulties were observed in defining an interface between a program backbone and the ORIGEN code respectively the ORIGENJUEL code. The estimation of the effort of refactoring the ORIGEN code or to write a new burnup code from scratch, led to the decision that it would be more efficient writing a new code, which could benefit from existing programming and software engineering tools from the computer code side and which can use the latest knowledge of nuclear reactions, e.g. consider all documented reaction channels. Therefore a new code with an object-oriented approach was developed at IEK-6. Object-oriented programming is currently state of the art and provides mostly an improved extensibility and maintainability. The new code was named TNT which stands for Topological Nuclide Transformation, since the code makes use of the real topology of the nuclear reactions. Here we want to present some first validation results from code to code benchmarks with the codes ORIGEN V2.2 and FISPACT2005 and whenever possible analytical results also used for the comparison. The 2 reference codes were chosen due to their high reputation in the field of fission reactor analysis (ORIGEN) and fusion facilities (FISPACT). (orig.)
The MCNPX Monte Carlo Radiation Transport Code
MCNPX (Monte Carlo N-Particle eXtended) is a general-purpose Monte Carlo radiation transport code with three-dimensional geometry and continuous-energy transport of 34 particles and light ions. It contains flexible source and tally options, interactive graphics, and support for both sequential and multi-processing computer platforms. MCNPX is based on MCNP4c and has been upgraded to most MCNP5 capabilities. MCNP is a highly stable code tracking neutrons, photons and electrons, and using evaluated nuclear data libraries for low-energy interaction probabilities. MCNPX has extended this base to a comprehensive set of particles and light ions, with heavy ion transport in development. Models have been included to calculate interaction probabilities when libraries are not available. Recent additions focus on the time evolution of residual nuclei decay, allowing calculation of transmutation and delayed particle emission. MCNPX is now a code of great dynamic range, and the excellent neutronics capabilities allow new opportunities to simulate devices of interest to experimental particle physics, particularly calorimetry. This paper describes the capabilities of the current MCNPX version 2.6.C, and also discusses ongoing code development
THE MCNPX MONTE CARLO RADIATION TRANSPORT CODE
WATERS, LAURIE S. [Los Alamos National Laboratory; MCKINNEY, GREGG W. [Los Alamos National Laboratory; DURKEE, JOE W. [Los Alamos National Laboratory; FENSIN, MICHAEL L. [Los Alamos National Laboratory; JAMES, MICHAEL R. [Los Alamos National Laboratory; JOHNS, RUSSELL C. [Los Alamos National Laboratory; PELOWITZ, DENISE B. [Los Alamos National Laboratory
2007-01-10
MCNPX (Monte Carlo N-Particle eXtended) is a general-purpose Monte Carlo radiation transport code with three-dimensional geometry and continuous-energy transport of 34 particles and light ions. It contains flexible source and tally options, interactive graphics, and support for both sequential and multi-processing computer platforms. MCNPX is based on MCNP4B, and has been upgraded to most MCNP5 capabilities. MCNP is a highly stable code tracking neutrons, photons and electrons, and using evaluated nuclear data libraries for low-energy interaction probabilities. MCNPX has extended this base to a comprehensive set of particles and light ions, with heavy ion transport in development. Models have been included to calculate interaction probabilities when libraries are not available. Recent additions focus on the time evolution of residual nuclei decay, allowing calculation of transmutation and delayed particle emission. MCNPX is now a code of great dynamic range, and the excellent neutronics capabilities allow new opportunities to simulate devices of interest to experimental particle physics; particularly calorimetry. This paper describes the capabilities of the current MCNPX version 2.6.C, and also discusses ongoing code development.
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)
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)
Criticality benchmarking of ANET Monte Carlo code
In this work the new Monte Carlo code ANET is tested on criticality calculations. ANET is developed based on the high energy physics code GEANT of CERN and aims at progressively satisfying several requirements regarding both simulations of GEN II/III reactors, as well as of innovative nuclear reactor designs such as the Accelerator Driven Systems (ADSs). Here ANET is applied on three different nuclear configurations, including a subcritical assembly, a Material Testing Reactor and the conceptual configuration of an ADS. In the first case, calculation of the effective multiplication factor (keff) are performed for the Training Nuclear Reactor of the Aristotle University of Thessaloniki, while in the second case keff is computed for the fresh fueled core of the Portuguese research reactor (RPJ) just after its conversion to Low Enriched Uranium, considering the control rods at the position that renders the reactor critical. In both cases ANET computations are compared with corresponding results obtained by three different well established codes, including both deterministic (XSDRNPM/CITATION) and Monte Carlo (TRIPOLI, MCNP). In the RPI case, keff computations are also compared with observations during the reactor core commissioning since the control rods are considered at criticality position. The above verification studies show ANET to produce reasonable results since they are satisfactorily compared with other models as well as with observations. For the third case (ADS), preliminary ANET computations of keff for various intensities of the proton beam are presented, showing also a reasonable code performance concerning both the order of magnitude and the relative variation of the computed parameter. (author)
In order to account for the reactivity-reducing effect of burn-up in the criticality safety analysis for systems with irradiated nuclear fuel (''burnup credit''), numerical methods to determine the enrichment and burnup dependent nuclide inventory (''burnup code'') and its resulting multiplication factor keff (''criticality code'') are applied. To allow for reliable conclusions, for both calculation systems the systematic deviations of the calculation results from the respective true values, the bias and its uncertainty, are being quantified by calculation and analysis of a sufficient number of suitable experiments. This quantification is specific for the application case under scope and is also called validation. GRS has developed a methodology to validate a calculation system for the application of burnup credit in the criticality safety analysis for irradiated fuel assemblies from pressurized water reactors. This methodology was demonstrated by applying the GRS home-built KENOREST burnup code and the criticality calculation sequence CSAS5 from SCALE code package. It comprises a bounding approach and alternatively a stochastic, which both have been exemplarily demonstrated by use of a generic spent fuel pool rack and a generic dry storage cask, respectively. Based on publicly available post irradiation examination and criticality experiments, currently the isotopes of uranium and plutonium elements can be regarded for.
LOLA-SYSTEM, JEN-UPM PWR Fuel Management System Burnup Code System
1 - Description of program or function: The LOLA-SYSTEM is a part of the JEN-UPM code package for PWR fuel management, scope or design calculations. It is a code package for core burnup calculations using nodal theory based on a FLARE type code. The LOLA-SYSTEM includes four modules: the first one (MELON-3) generates the constants of the K-inf and M2 correlations to be input into SIMULA-3. It needs the K-inf and M2 fuel assembly values at different conditions of moderator temperature, Boron concentration, burnup, etc., which are provided by MARIA fuel assembly calculations. The main module (SIMULA-3) is the core burnup calculation code in three dimensions and one group of energy. It normally uses a geometrical representation of one node per fuel assembly or per quarter of fuel assembly. It has included a thermal hydraulic feedback on flow and voids and criticality searches on boron concentration and control rods insertion. The CONCON code makes the calculation of the albedo, transport factors, K-inf and M2 correction factors to be input into SIMULA-3. The calculation is made in the XY transversal plane. The CONAXI code is similar to CONCON, but in the axial direction. 2 - Method of solution: MELON-3 makes a mean squares fit of K-inf and M2 values at different conditions in order to determine the constants of the feedback correlations. SIMULA-3 uses a modified one-group nodal theory, with a new transport kernel that provides the same node interface leakages as a fine mesh diffusion calculation. CONCON and CONAXI determine the transport and correction factors, as well as the albedo, to be input into SIMULA-3. They are determined by a method of leakages equivalent to the detailed diffusion calculation of CARMEN or VENTURE; these factors also include the heterogeneity effects inside the node. 3 - Restrictions on the complexity of the problem: Number of axial nodes less than or equal 34. Number of material types less than or equal 30. Number of fuel assembly types less
High burnup fuel simulation during a reactivity initiated transient using the FRAPTRAN code
CIP01 test simulated a RIA (Reactivity Initiated Accident) transient at the experimental CABRI reactor under sodium cooling conditions on November 29th 2001. The father rod was fabricated by ENUSA and irradiated during five cycles in Vandellos 2 NPP and it reached a maximum burnup at EOL (End of Life) of 75 GWd/tU. The test rodlet was refabricated in Studsvik from the span 5 of the father rod. The cladding material was ZIRLO with a mean oxide layer of 75 μm. This note describes the modelling of the CIP01 test carried out in the CABRI reactor under RIA conditions by the use of the transient FRAPTRAN code (1.2 version). The CIP01 test consisted in a power pulse of 32.4 ms with a deposited energy of 98.2 cal/g at PPN (Peak Power Node). The code results will be presented and compared to the experimental data available. In addition, the FRAPTRAN and the SCANAIR code estimations will be also contrasted. This work is framed within the CSN-CIEMAT agreement on 'Thermo-Mechanical Behaviour of the Nuclear Fuel at High Burnup'. (Author)
Startup of “CANDLE” burnup in a Gas-cooled Fast Reactor using Monte Carlo method
Highlights: ► In equilibrium state of a CANDLE core, the burning region contains fission products and actinides. ► These isotopes are not available. The solution is startup of a reactor using easily available materials. ► At the end of core life the fuel for the equilibrium core is produced. ► In this work the startup of a CANDLE-GFR has been evaluated using Monte Carlo technique. ► The results show that the equilibrium state could be achieved after some minor transients. -- Abstract: During the past decade, the CANDLE burnup strategy has been proposed as an innovative fuel cycle and reactor design for complete utilization of uranium resources. In this strategy the shapes of neutron flux, nuclide densities and power density distribution remain constant but the burning region moves in axial direction. The feasibility of this strategy has been demonstrated widely by using the diffusion technique in conjunction with nuclide transmutation equations. On the other hand since the Monte Carlo method provides the exact solution to the neutron transport, the Monte Carlo technique is becoming more widely used in routine burnup calculations. The main objective of this work is startup of CANDLE burnup in a Gas cooled Fast Reactor using a Monte Carlo burnup scheme. In this case only natural or depleted uranium is required for fresh fuel region. However, the construction of the first CANDLE core is faced with a big problem. In equilibrium state the burning region contains a spectrum of fission products as well as higher actinides. These isotopes are not easily available for constructing the initial CANDLE core. The solution is startup of a special reactor using the enriched uranium in starter zone. At the end of core life the fuel for the next core is produced with the composition close to the equilibrium state. An originally MCNP–ORIGEN linkage program named MOBC has been used for criticality and isotopic evaluation of the core. The results of analysis showed that
MORSE Monte Carlo radiation transport code system
This report is an addendum to the MORSE report, ORNL-4972, originally published in 1975. This addendum contains descriptions of several modifications to the MORSE Monte Carlo Code, replacement pages containing corrections, Part II of the report which was previously unpublished, and a new Table of Contents. The modifications include a Klein Nishina estimator for gamma rays. Use of such an estimator required changing the cross section routines to process pair production and Compton scattering cross sections directly from ENDF tapes and writing a new version of subroutine RELCOL. Another modification is the use of free form input for the SAMBO analysis data. This required changing subroutines SCORIN and adding new subroutine RFRE. References are updated, and errors in the original report have been corrected
Burnup span sensitivity analysis of different burnup coupling schemes
Highlights: ► The objective of this work is the burnup span sensitivity analysis of different coupling schemes. ► Three kinds of schemes have been implemented in a new MCNP–ORIGEN linkage program. ► Two kinds of schemes are based predictor–corrector technique and the third is based on Euler explicit method. ► The analysis showed that the predictor–corrector approach better accounts for nonlinear behavior of burnup. ► It is sufficiently good to use the Euler method at small spans but for large spans use of second order scheme is mandatory. - Abstract: The analysis of core composition changes is complicated by the fact that the time and spatial variations in isotopic composition depend on the neutron flux distribution and vice versa. Fortunately, changes in core composition occur relatively slowly and hence the burnup analysis can be performed by dividing the burnup period into some burnup spans and assuming that the averaged flux and cross sections are constant during each burn up span. The burnup span sensitivity analysis attempts to find how much the burnup spans could be increased without any significant change in results. This goal has been achieved by developing a new MCNP–ORIGEN linkage program named MOBC (MCNP–ORIGEN Burnup Calculation). Three kinds of coupling scheme have been implemented in MOBC. Two of these are based on second order predictor–corrector technique and enable us to choose larger time steps, whilst the third one is based on Euler explicit first order method and is faster than the other two. The validity of the developed program has been evaluated by the code vs. code comparison technique. Two different types of codes are employed. The first one is based on deterministic two dimensional transport method, like CASMO-4 and HELIOS codes, and the second one is based on Monte Carlo method, like MCODE code. Only one coupling technique is employed in each of these state of the art codes, while the MOBC excels in its ability to
Development of a code and models for high burnup fuel performance analysis
First the high burnup LWR fuel behavior is discussed and necessary models for the analysis are reviewed. These aspects of behavior are the changes of power history due to the higher enrichment, the temperature feedback due to fission gas release and resultant degradation of gap conductance, axial fission gas transport in fuel free volume, fuel conductivity degradation due to fission product solution and modification of fuel micro-structure. Models developed for these phenomena, modifications in the code, and the benchmark results mainly based on Risoe fission gas project is presented. Finally the rim effect which is observe only around the fuel periphery will be discussed focusing into the fuel conductivity degradation and swelling due to the porosity development. (author). 18 refs, 13 figs, 3 tabs
The predictions of the ENIGMA code have been compared with data from high burn-up fuel experiments from the Halden and RISO reactors. The experiments modelled were IFA-504 and IFA-558 from Halden and the test II-5 from the RISO power burnup test series. The code has well modelled the fuel thermal performance and has provided a good measure of iodine release from pre-interlinked fuel. After interlinkage the iodine predictions remain a good fit for one experiment, but there is significant overprediction for a second experiment (IFA-558). Stable fission gas release is also well modelled and the predictions are within the expected uncertainly band throughout the burn-up range. This report presents code predictions for stable fission gas release to 40GWd/tU, iodine release measurements to 50GWd/tU and thermal performance (fuel centre temperature) to 55GWd/tU. Fuel ratings of up to 38kW/m were modelled at the high burn-up levels. The code is shown to accurately or conservatively predict all these parameters. (author). 1 ref., 6 figs
Criticality qualification of a new Monte Carlo code for reactor core analysis
In order to accurately simulate Accelerator Driven Systems (ADS), the utilization of at least two computational tools is necessary (the thermal-hydraulic problem is not considered in the frame of this work), namely: (a) A High Energy Physics (HEP) code system dealing with the 'Accelerator part' of the installation, i.e. the computation of the spectrum, intensity and spatial distribution of the neutrons source created by (p, n) reactions of a proton beam on a target and (b) a neutronics code system, handling the 'Reactor part' of the installation, i.e. criticality calculations, neutron transport, fuel burn-up and fission products evolution. In the present work, a single computational tool, aiming to analyze an ADS in its integrity and also able to perform core analysis for a conventional fission reactor, is proposed. The code is based on the well qualified HEP code GEANT (version 3), transformed to perform criticality calculations. The performance of the code is tested against two qualified neutronics code systems, the diffusion/transport SCALE-CITATION code system and the Monte Carlo TRIPOLI code, in the case of a research reactor core analysis. A satisfactory agreement was exhibited by the three codes.
On the use of SERPENT Monte Carlo code to generate few group diffusion constants
Piovezan, Pamela, E-mail: pamela.piovezan@ctmsp.mar.mil.b [Centro Tecnologico da Marinha em Sao Paulo (CTMSP), Sao Paulo, SP (Brazil); Carluccio, Thiago; Domingos, Douglas Borges; Rossi, Pedro Russo; Mura, Luiz Felipe, E-mail: fermium@cietec.org.b, E-mail: thiagoc@ipen.b [Fermium Tecnologia Nuclear, Sao Paulo, SP (Brazil); Instituto de Pesquisas Energeticas e Nucleares (IPEN/CNEN-SP), Sao Paulo, SP (Brazil)
2011-07-01
The accuracy of diffusion reactor codes strongly depends on the quality of the groups constants processing. For many years, the generation of such constants was based on 1-D infinity cell transport calculations. Some developments using collision probability or the method of characteristics allow, nowadays, 2-D assembly group constants calculations. However, these 1-D and 2-D codes how some limitations as , for example, on complex geometries and in the neighborhood of heavy absorbers. On the other hand, since Monte Carlos (MC) codes provide accurate neutro flux distributions, the possibility of using these solutions to provide group constants to full-core reactor diffusion simulators has been recently investigated, especially for the cases in which the geometry and reactor types are beyond the capability of the conventional deterministic lattice codes. The two greatest difficulties on the use of MC codes to group constant generation are the computational costs and the methodological incompatibility between analog MC particle transport simulation and deterministic transport methods based in several approximations. The SERPENT code is a 3-D continuous energy MC transport code with built-in burnup capability that was specially optimized to generate these group constants. In this work, we present the preliminary results of using the SERPENT MC code to generate 3-D two-group diffusion constants for a PWR like assembly. These constants were used in the CITATION diffusion code to investigate the effects of the MC group constants determination on the neutron multiplication factor diffusion estimate. (author)
Criticality qualification of a new Monte Carlo code for reactor core analysis
Catsaros, N. [Institute of Nuclear Technology - Radiation Protection, NCSR ' DEMOKRITOS' , P.O. Box 60228, 15310 Aghia Paraskevi (Greece); Gaveau, B. [MAPS, Universite Paris VI, 4 Place Jussieu, 75005 Paris (France); Jaekel, M. [Laboratoire de Physique Theorique, Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris (France); Maillard, J. [MAPS, Universite Paris VI, 4 Place Jussieu, 75005 Paris (France); CNRS-IDRIS, Bt 506, BP167, 91403 Orsay (France); CNRS-IN2P3, 3 rue Michel Ange, 75794 Paris (France); Maurel, G. [Faculte de Medecine, Universite Paris VI, 27 rue de Chaligny, 75012 Paris (France); MAPS, Universite Paris VI, 4 Place Jussieu, 75005 Paris (France); Savva, P., E-mail: savvapan@ipta.demokritos.g [Institute of Nuclear Technology - Radiation Protection, NCSR ' DEMOKRITOS' , P.O. Box 60228, 15310 Aghia Paraskevi (Greece); Silva, J. [MAPS, Universite Paris VI, 4 Place Jussieu, 75005 Paris (France); Varvayanni, M.; Zisis, Th. [Institute of Nuclear Technology - Radiation Protection, NCSR ' DEMOKRITOS' , P.O. Box 60228, 15310 Aghia Paraskevi (Greece)
2009-11-15
In order to accurately simulate Accelerator Driven Systems (ADS), the utilization of at least two computational tools is necessary (the thermal-hydraulic problem is not considered in the frame of this work), namely: (a) A High Energy Physics (HEP) code system dealing with the 'Accelerator part' of the installation, i.e. the computation of the spectrum, intensity and spatial distribution of the neutrons source created by (p, n) reactions of a proton beam on a target and (b) a neutronics code system, handling the 'Reactor part' of the installation, i.e. criticality calculations, neutron transport, fuel burn-up and fission products evolution. In the present work, a single computational tool, aiming to analyze an ADS in its integrity and also able to perform core analysis for a conventional fission reactor, is proposed. The code is based on the well qualified HEP code GEANT (version 3), transformed to perform criticality calculations. The performance of the code is tested against two qualified neutronics code systems, the diffusion/transport SCALE-CITATION code system and the Monte Carlo TRIPOLI code, in the case of a research reactor core analysis. A satisfactory agreement was exhibited by the three codes.
The high burnup database for CANDU fuel includes several cases from both power station and experimental reactor irradiations, with achieved burnups of up to 800 MW.h/kgU. The power history for each of these cases is different, encompassing low steady-state, declining, and power-ramps. This variety offers a good opportunity to check the models of fuel behaviour, and to identify areas for improvement. The main parameters for comparing calculated versus measured data are the fission gas release and the sheath hoop strain. Good agreement of calculated values of these two parameters with experimental data indicates that the global behaviour of the fuel element is adequately simulated by our codes. The ELESIM computer code was used as the simulation tool. The models for fission gas release, swelling and for fuel pellet expansion were thoroughly analysed. Changes were proposed for both models. The fuel pellet expansion model was modified to account for gaseous swelling, which becomes very important at high burnups. As well, the mathematics of the fission gas release model was upgraded for the diffusional release of fission gas atoms to the grain boundaries. A revised version of the ELESIM computer code was used to simulate the cases from the high burnup database. Satisfactory agreement was found for most cases. The discrepancies are discussed in view of alternative mechanisms that can operate and be enhanced at high burnup. These include stoichiometry changes with burnup that affects fission gas release, and also outer pellet rim fission gas release by a grain boundary diffusion process. The main conclusion of this study is that the revised version of the ELESIM code is able to simulate with reasonable accuracy high burnup as well as low burnup CANDU fuel. This includes irradiations of steady-state, declining, or ramped fuel power histories with a prolonged hold at high power. However, future improvements to ELESIM are needed to model fuel power histories with short dwell
A PWR PCI failure criterion to burnups of 60 GW·d/t using the ENIGMA code
A fuel performance modelling code (ENIGMA) has been used to analyse the empirical PCI failure criterion in terms of a clad failure stress as a function of burnup and fast neutron dose. The Studsvik database has been analysed. Results indicate a rising and then saturating failure stress with burnup and fast neutron dose. Using the PCI failure limits, equivalent to 95/95 confidence limits, an ENIGMA stress-based methodology is used to derive PWR PCI failure limits up to 60 GW·d/t U using a conservative assumption that the failure stress does not increase at high burnup and neutron dose. In addition experimental ramp data on gadolinia-doped fuel rods do not indicate any increased susceptibility to PCI failure implying that the UO2 criterion can be used for gadolinia doped fuel. (author)
By coupling a three-dimensional Discrete Ordinate (SN) code with an Exponential Euler Method inventory code, a three-dimensional code system was developed for accurate and rapid fuel burnup and material activation analyses. Resonance self-shielding effect was corrected dynamically in each step and each fuel zone. The update of compositions and nuclide densities in transport calculation was treated by setting nuclide density limitation, keeping important nuclides and using pseudo fission products (FPeff). The validation was described, which illustrates the robustness and reliability of this code system. (authors)
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
Fast code for Monte Carlo simulations
A computer code to generate the dynamic evolution of the Ising model on a square lattice, following the Metropolis algorithm is presented. The computer time consumption is reduced by a factor of 8 when one compares our code with traditional multiple spin codes. The memory allocation size is also reduced by a factor of 4. The code is easily generalizable for other lattices and models. (author)
MORET: Version 4.B. A multigroup Monte Carlo criticality code
MORET 4 is a three dimensional multigroup Monte Carlo code which calculates the effective multiplication factor (keff) of any configurations more or less complex as well as reaction rates in the different volumes of the geometry and the leakage out of the system. MORET 4 is the Monte Carlo code of the APOLLO2-MORET 4 standard route of CRISTAL, the French criticality package. It is the most commonly used Monte Carlo code for French criticality calculations. During the last four years, the MORET 4 team has developed or improved the following major points: modernization of the geometry, implementation of perturbation algorithms, source distribution convergence, statistical detection of stationarity, unbiased variance estimation and creation of pre-processing and post-processing tools. The purpose of this paper is not only to present the new features of MORET but also to detail clearly the physical models and the mathematical methods used in the code. (author)
MOx benchmark calculations by deterministic and Monte Carlo codes
Highlights: ► MOx based depletion calculation. ► Methodology to create continuous energy pseudo cross section for lump of minor fission products. ► Mass inventory comparison between deterministic and Monte Carlo codes. ► Higher deviation was found for several isotopes. - Abstract: A depletion calculation benchmark devoted to MOx fuel is an ongoing objective of the OECD/NEA WPRS following the study of depletion calculation concerning UOx fuels. The objective of the proposed benchmark is to compare existing depletion calculations obtained with various codes and data libraries applied to fuel and back-end cycle configurations. In the present work the deterministic code NEWT/ORIGEN-S of the SCALE6 codes package and the Monte Carlo based code MONTEBURNS2.0 were used to calculate the masses of inventory isotopes. The methodology to apply the MONTEBURNS2.0 to this benchmark is also presented. Then the results from both code were compared.
A New Monte Carlo Neutron Transport Code at UNIST
Monte Carlo neutron transport code named MCS is under development at UNIST for the advanced reactor design and research purpose. This MC code can be used for fixed source calculation and criticality calculation. Continuous energy neutron cross section data and multi-group cross section data can be used for the MC calculation. This paper presents the overview of developed MC code and its calculation results. The real time fixed source calculation ability is also tested in this paper. The calculation results show good agreement with commercial code and experiment. A new Monte Carlo neutron transport code is being developed at UNIST. The MC codes are tested with several benchmark problems: ICSBEP, VENUS-2, and Hoogenboom-Martin benchmark. These benchmarks covers pin geometry to 3-dimensional whole core, and results shows good agreement with reference results
Benchmarking Monte Carlo codes for criticality safety using subcritical measurements
Monte Carlo codes that are used for criticality safety evaluations are typically validated using critical experiments in which the neutron multiplication factor is unity. However, the conditions for most fissile material operations do not coincide to those of the critical experiments. This paper demonstrates that Monte Carlo methods and nuclear data can be validated using subcritical measurements whose conditions may coincide more closely to actual configurations of fissile material. (orig.)
Design of shielding of LILW containers by Monte Carlo codes
Accurate prediction of dose rates from containers with radioactive waste is becoming more important regarding more rigorous regulative in this area. The usual approach to the problem consists in combining numerical and measuring methods. In this paper a Monte Carlo calculations were used for calculating doses from a standard 200 liter drum which contains the intermediate level radioactive waste. Two different Monte Carlo codes were applied and compared, for the same combination of parameters. (author)
Verification of the spectral history correction method with fully coupled Monte Carlo code BGCore
Recently, a new method for accounting for burnup history effects on few-group cross sections was developed and implemented in the reactor dynamic code DYN3D. The method relies on the tracking of the local Pu-239 density which serves as an indicator of burnup spectral history. The validity of the method was demonstrated in PWR and VVER applications. However, the spectrum variation in BWR core is more pronounced due to the stronger coolant density change. Therefore, the purpose of the current work is to further investigate the applicability of the method to BWR analysis. The proposed methodology was verified against recently developed BGCore system, which couples Monte Carlo neutron transport with depletion and thermal hydraulic solvers and thus capable of providing a reference solution for 3D simulations. The results dearly show that neglecting the spectral history effects leads to a very large deviation (e.g. 2000 pcm in reactivity) from fee reference solution. However, a very good agreement between DYN3D and BGCore is observed (on the order of 200 pcm in reactivity), when the. Pu-correction method is applied. (author)
Current status of the PSG Monte Carlo neutron transport code
PSG is a new Monte Carlo neutron transport code, developed at the Technical Research Centre of Finland (VTT). The code is mainly intended for fuel assembly-level reactor physics calculations, such as group constant generation for deterministic reactor simulator codes. This paper presents the current status of the project and the essential capabilities of the code. Although the main application of PSG is in lattice calculations, the geometry is not restricted in two dimensions. This paper presents the validation of PSG against the experimental results of the three-dimensional MOX fuelled VENUS-2 reactor dosimetry benchmark. (authors)
Vectorization techniques for neutron transport Monte Carlo codes
Four Monte Carlo codes, KENO IV, MORSE-DD, MCNP and VIM, have been vectorized already at JAERI Computing Center aiming at an increase in clculation performance, and speed-up ratios of vectorized codes to the original ones were found to be low values between 1.3 and 1.5. In this report the vectorization processes for these four codes are reviewed comprehensively, and methods of analysis for vectorization, modification of control structures of codes and debugging techniques are discussed. The reason for low speed-up ratios is also discussed. (author)
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 Monte Carlo code for ion beam therapy
Anaïs Schaeffer
2012-01-01
Initially developed for applications in detector and accelerator physics, the modern Fluka Monte Carlo code is now used in many different areas of nuclear science. Over the last 25 years, the code has evolved to include new features, such as ion beam simulations. Given the growing use of these beams in cancer treatment, Fluka simulations are being used to design treatment plans in several hadron-therapy centres in Europe. Fluka calculates the dose distribution for a patient treated at CNAO with proton beams. The colour-bar displays the normalized dose values. Fluka is a Monte Carlo code that very accurately simulates electromagnetic and nuclear interactions in matter. In the 1990s, in collaboration with NASA, the code was developed to predict potential radiation hazards received by space crews during possible future trips to Mars. Over the years, it has become the standard tool to investigate beam-machine interactions, radiation damage and radioprotection issues in the CERN accelerator com...
Development of a New Monte Carlo reactor physics code
Leppänen, Jaakko
2007-01-01
Monte Carlo neutron transport codes are widely used in various reactor physics applications, traditionally related to criticality safety analyses, radiation shielding problems, detector modelling and validation of deterministic transport codes. The main advantage of the method is the capability to model geometry and interaction physics without major approximations. The disadvantage is that the modelling of complicated systems is very computing-intensive, which restricts the applications to so...
The Monte Carlo code TRAMO - Capabilities and instructions for application
The report is intended for readers familiar with the fundamentals of the Monte Carlo method. Those readers might be interested in learning about successful generalisations as well as new ideas for curbing the statistical errors involved. Another intention however is to explain the significant basic features of the multigroup Monte Carlo code TRAMO, including the required input, so that readers will be able to performing the required adjustments to the specific calculation technique and develop their own tools for performing their specific calculations. An indispensable code needed for such TRAMO applications is the TRAWEI Monte Carlo code which calculates he required weightings for applications of the variance reducing Weight Window Method; other codes required are those for generating the neutron cross-section data and the group data. The TRAMO code calculates, with given source distribution of neutrons in multigroup approximation, multigroup flux data, integrated group flux data, and dose values for given partial volumes and surfaces. There are further code versions for calculation of neutron and gamma fluxes, or criticality data, but these are not considered in the report. (orig./CB)
MORSE Monte Carlo radiation transport code system
For a number of years the MORSE user community has requested additional help in setting up problems using various options. The sample problems distributed with MORSE did not fully demonstrate the capability of the code. At Oak Ridge National Laboratory the code originators had a complete set of sample problems, but funds for documenting and distributing them were never available. Recently the number of requests for listings of input data and results for running some particular option the user was trying to implement has increased to the point where it is not feasible to handle them on an individual basis. Consequently it was decided to package a set of sample problems which illustrates more adequately how to run MORSE. This write-up may be added to Part III of the MORSE report. These sample problems include a combined neutron-gamma case, a neutron only case, a gamma only case, an adjoint case, a fission case, a time-dependent fission case, the collision density case, an XCHEKR run and a PICTUR run
Taylor series development in the Monte Carlo code Tripoli-4
Mazzolo, Alain; Zoia, Andrea; Martin, Brunella
2014-06-01
Perturbation methods for one or several variables based on the Taylor series development up to the second order is presented for the collision estimator in the framework of the Monte Carlo code Tripoli-4. Comparisons with the correlated sampling method implemented in Tripoli-4 demonstrate the need of including the cross derivatives in the development.
Behler, Matthias; Hannstein, Volker; Kilger, Robert; Moser, Franz-Eberhard; Pfeiffer, Arndt; Stuke, Maik
2014-06-15
In order to account for the reactivity-reducing effect of burn-up in the criticality safety analysis for systems with irradiated nuclear fuel (''burnup credit''), numerical methods to determine the enrichment and burnup dependent nuclide inventory (''burnup code'') and its resulting multiplication factor k{sub eff} (''criticality code'') are applied. To allow for reliable conclusions, for both calculation systems the systematic deviations of the calculation results from the respective true values, the bias and its uncertainty, are being quantified by calculation and analysis of a sufficient number of suitable experiments. This quantification is specific for the application case under scope and is also called validation. GRS has developed a methodology to validate a calculation system for the application of burnup credit in the criticality safety analysis for irradiated fuel assemblies from pressurized water reactors. This methodology was demonstrated by applying the GRS home-built KENOREST burnup code and the criticality calculation sequence CSAS5 from SCALE code package. It comprises a bounding approach and alternatively a stochastic, which both have been exemplarily demonstrated by use of a generic spent fuel pool rack and a generic dry storage cask, respectively. Based on publicly available post irradiation examination and criticality experiments, currently the isotopes of uranium and plutonium elements can be regarded for.
The effect of fuel burn-up on the radioactivation behavior of cladding hull materials was investigated using the ORIGEN-S code for various materials of Zircaloy-4, Zirlo, HANA-4, and HANA-6 and for various fuel burn-ups of 30, 45, 60, and 75 GWD/MTU. The Zircaloy-4 material is the only one that does not contain Nb as an alloy constituent, and it was revealed that 125Sb, 125mTe, and 55Fe are the major sources of radioactivity. On the other hand, 93mNb was identified as the most radioactive nuclide for the other materials although minor radioactive nuclides varied owing to their different initial constituents. The radioactivity of 94Nb was of particular focus owing to its acceptance limit against a Korean intermediate-/low-level waste repository. The radioactivation calculation results revealed that only Zircaloy-4 is acceptable for the Korean repository, while the other materials required at least 4,900 of Nb decontamination factor owing to the high radioactivity of 94Nb regardless of the fuel burn-up. A discussion was also made on the feasibility of Zr recovery methods (chlorination and electrorefining) for selective recovery of Zr so that it can be disposed of in the Korean repository. (author)
Acceleration of a Monte Carlo radiation transport code
Execution time for the Integrated TIGER Series (ITS) Monte Carlo radiation transport code has been reduced by careful re-coding of computationally intensive subroutines. Three test cases for the TIGER (1-D slab geometry), CYLTRAN (2-D cylindrical geometry), and ACCEPT (3-D arbitrary geometry) codes were identified and used to benchmark and profile program execution. Based upon these results, sixteen top time-consuming subroutines were examined and nine of them modified to accelerate computations with equivalent numerical output to the original. The results obtained via this study indicate that speedup factors of 1.90 for the TIGER code, 1.67 for the CYLTRAN code, and 1.11 for the ACCEPT code are achievable. copyright 1996 American Institute of Physics
Study on random number generator in Monte Carlo code
The Monte Carlo code uses a sequence of pseudo-random numbers with a random number generator (RNG) to simulate particle histories. A pseudo-random number has its own period depending on its generation method and the period is desired to be long enough not to exceed the period during one Monte Carlo calculation to ensure the correctness especially for a standard deviation of results. The linear congruential generator (LCG) is widely used as Monte Carlo RNG and the period of LCG is not so long by considering the increasing rate of simulation histories in a Monte Carlo calculation according to the remarkable enhancement of computer performance. Recently, many kinds of RNG have been developed and some of their features are better than those of LCG. In this study, we investigate the appropriate RNG in a Monte Carlo code as an alternative to LCG especially for the case of enormous histories. It is found that xorshift has desirable features compared with LCG, and xorshift has a larger period, a comparable speed to generate random numbers, a better randomness, and good applicability to parallel calculation. (author)
A semianalytic Monte Carlo code for modelling LIDAR measurements
Palazzi, Elisa; Kostadinov, Ivan; Petritoli, Andrea; Ravegnani, Fabrizio; Bortoli, Daniele; Masieri, Samuele; Premuda, Margherita; Giovanelli, Giorgio
2007-10-01
LIDAR (LIght Detection and Ranging) is an optical active remote sensing technology with many applications in atmospheric physics. Modelling of LIDAR measurements appears useful approach for evaluating the effects of various environmental variables and scenarios as well as of different measurement geometries and instrumental characteristics. In this regard a Monte Carlo simulation model can provide a reliable answer to these important requirements. A semianalytic Monte Carlo code for modelling LIDAR measurements has been developed at ISAC-CNR. The backscattered laser signal detected by the LIDAR system is calculated in the code taking into account the contributions due to the main atmospheric molecular constituents and aerosol particles through processes of single and multiple scattering. The contributions by molecular absorption, ground and clouds reflection are evaluated too. The code can perform simulations of both monostatic and bistatic LIDAR systems. To enhance the efficiency of the Monte Carlo simulation, analytical estimates and expected value calculations are performed. Artificial devices (such as forced collision, local forced collision, splitting and russian roulette) are moreover foreseen by the code, which can enable the user to drastically reduce the variance of the calculation.
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)
TRIPOLI-3: a neutron/photon Monte Carlo transport code
The present version of TRIPOLI-3 solves the transport equation for coupled neutron and gamma ray problems in three dimensional geometries by using the Monte Carlo method. This code is devoted both to shielding and criticality problems. The most important feature for particle transport equation solving is the fine treatment of the physical phenomena and sophisticated biasing technics useful for deep penetrations. The code is used either for shielding design studies or for reference and benchmark to validate cross sections. Neutronic studies are essentially cell or small core calculations and criticality problems. TRIPOLI-3 has been used as reference method, for example, for resonance self shielding qualification. (orig.)
SPAMCART: a code for smoothed particle Monte Carlo radiative transfer
Lomax, O
2016-01-01
We present a code for generating synthetic SEDs and intensity maps from Smoothed Particle Hydrodynamics simulation snapshots. The code is based on the Lucy (1999) Monte Carlo Radiative Transfer method, i.e. it follows discrete luminosity packets, emitted from external and/or embedded sources, as they propagate through a density field, and then uses their trajectories to compute the radiative equilibrium temperature of the ambient dust. The density is not mapped onto a grid, and therefore the calculation is performed at exactly the same resolution as the hydrodynamics. We present two example calculations using this method. First, we demonstrate that the code strictly adheres to Kirchhoff's law of radiation. Second, we present synthetic intensity maps and spectra of an embedded protostellar multiple system. The algorithm uses data structures that are already constructed for other purposes in modern particle codes. It is therefore relatively simple to implement.
ORPHEE research reactor: 3D core depletion calculation using Monte-Carlo code TRIPOLI-4®
Damian, F.; Brun, E.
2014-06-01
ORPHEE is a research reactor located at CEA Saclay. It aims at producing neutron beams for experiments. This is a pool-type reactor (heavy water), and the core is cooled by light water. Its thermal power is 14 MW. ORPHEE core is 90 cm height and has a cross section of 27x27 cm2. It is loaded with eight fuel assemblies characterized by a various number of fuel plates. The fuel plate is composed of aluminium and High Enriched Uranium (HEU). It is a once through core with a fuel cycle length of approximately 100 Equivalent Full Power Days (EFPD) and with a maximum burnup of 40%. Various analyses under progress at CEA concern the determination of the core neutronic parameters during irradiation. Taking into consideration the geometrical complexity of the core and the quasi absence of thermal feedback for nominal operation, the 3D core depletion calculations are performed using the Monte-Carlo code TRIPOLI-4® [1,2,3]. A preliminary validation of the depletion calculation was performed on a 2D core configuration by comparison with the deterministic transport code APOLLO2 [4]. The analysis showed the reliability of TRIPOLI-4® to calculate a complex core configuration using a large number of depleting regions with a high level of confidence.
Highlights: • Coupling of Monte Carlo code Serpent and thermal–hydraulics code RELAP5. • A convergence criterion is developed based on the statistical uncertainty of power. • Correlation between MC statistical uncertainty and coupled error is quantified. • Both UO2 and MOX single assembly models are used in the coupled simulation. • Validation of coupling results with a multi-group transport code DeCART. - Abstract: Coupled multi-physics approach plays an important role in improving computational accuracy. Compared with deterministic neutronics codes, Monte Carlo codes have the advantage of a higher resolution level. In the present paper, a three-dimensional continuous-energy Monte Carlo reactor physics burnup calculation code, Serpent, is coupled with a thermal–hydraulics safety analysis code, RELAP5. The coupled Serpent/RELAP5 code capability is demonstrated by the improved axial power distribution of UO2 and MOX single assembly models, based on the OECD-NEA/NRC PWR MOX/UO2 Core Transient Benchmark. Comparisons of calculation results using the coupled code with those from the deterministic methods, specifically heterogeneous multi-group transport code DeCART, show that the coupling produces more precise results. A new convergence criterion for the coupled simulation is developed based on the statistical uncertainty in power distribution in the Monte Carlo code, rather than ad-hoc criteria used in previous research. The new convergence criterion is shown to be more rigorous, equally convenient to use but requiring a few more coupling steps to converge. Finally, the influence of Monte Carlo statistical uncertainty on the coupled error of power and thermal–hydraulics parameters is quantified. The results are presented such that they can be used to find the statistical uncertainty to use in Monte Carlo in order to achieve a desired precision in coupled simulation
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
Verification of Monte Carlo transport codes by activation experiments
With the increasing energies and intensities of heavy-ion accelerator facilities, the problem of an excessive activation of the accelerator components caused by beam losses becomes more and more important. Numerical experiments using Monte Carlo transport codes are performed in order to assess the levels of activation. The heavy-ion versions of the codes were released approximately a decade ago, therefore the verification is needed to be sure that they give reasonable results. Present work is focused on obtaining the experimental data on activation of the targets by heavy-ion beams. Several experiments were performed at GSI Helmholtzzentrum fuer Schwerionenforschung. The interaction of nitrogen, argon and uranium beams with aluminum targets, as well as interaction of nitrogen and argon beams with copper targets was studied. After the irradiation of the targets by different ion beams from the SIS18 synchrotron at GSI, the γ-spectroscopy analysis was done: the γ-spectra of the residual activity were measured, the radioactive nuclides were identified, their amount and depth distribution were detected. The obtained experimental results were compared with the results of the Monte Carlo simulations using FLUKA, MARS and SHIELD. The discrepancies and agreements between experiment and simulations are pointed out. The origin of discrepancies is discussed. Obtained results allow for a better verification of the Monte Carlo transport codes, and also provide information for their further development. The necessity of the activation studies for accelerator applications is discussed. The limits of applicability of the heavy-ion beam-loss criteria were studied using the FLUKA code. FLUKA-simulations were done to determine the most preferable from the radiation protection point of view materials for use in accelerator components.
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
Nanodosimetric verification in proton therapy: Monte Carlo Codes Comparison
Full text: Nanodosimetry strives to develop a novel dosimetry concept suitable for advanced modalities of cancer radiotherapy, such as proton therapy. This project aims to evaluate the plausibility of the physical models implemented in the Geant4 Very Low Energy (Geant4-DNA) extensions by comparing nanodosimetric quantities calculated with Geant4-DNA and the PTB Monte Carlo track structure code. Nanodosimetric track structure parameters were calculated for cylindrical targets representing DNA and nucleosome segments and converted into the probability of producing a DSB using the model proposed by Garty et al. [1]. Monoenergetic protons and electrons of energies typical for 6-electron spectra were considered as primary particles. Good agreement was found between the two codes for electrons of energies above 200 eV. Below this energy Geant4-DNA produced slightly higher numbers of ionisations in the sensitive volumes and higher probabilities for DSB formation. For protons, Geant4-DNA also gave higher numbers of ionisations and DSB probabilities, particularly in the low energy range, while a satisfactory agreement was found for energies higher than I MeV. Comparing two codes can be useful as any observed divergence in results between the two codes provides valuable information as to where further consideration of the underlying physical models used in each code may be required. Consistently it was seen that the largest difference between the codes was in the low energy ranges for each particle type. (author)
Proton therapy Monte Carlo SRNA-VOX code
Ilić Radovan D.
2012-01-01
Full Text Available The most powerful feature of the Monte Carlo method is the possibility of simulating all individual particle interactions in three dimensions and performing numerical experiments with a preset error. These facts were the motivation behind the development of a general-purpose Monte Carlo SRNA program for proton transport simulation in technical systems described by standard geometrical forms (plane, sphere, cone, cylinder, cube. Some of the possible applications of the SRNA program are: (a a general code for proton transport modeling, (b design of accelerator-driven systems, (c simulation of proton scattering and degrading shapes and composition, (d research on proton detectors; and (e radiation protection at accelerator installations. This wide range of possible applications of the program demands the development of various versions of SRNA-VOX codes for proton transport modeling in voxelized geometries and has, finally, resulted in the ISTAR package for the calculation of deposited energy distribution in patients on the basis of CT data in radiotherapy. All of the said codes are capable of using 3-D proton sources with an arbitrary energy spectrum in an interval of 100 keV to 250 MeV.
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)
Computed radiography simulation using the Monte Carlo code MCNPX
Simulating x-ray images has been of great interest in recent years as it makes possible an analysis of how x-ray images are affected owing to relevant operating parameters. In this paper, a procedure for simulating computed radiographic images using the Monte Carlo code MCNPX is proposed. The sensitivity curve of the BaFBr image plate detector as well as the characteristic noise of a 16-bit computed radiography system were considered during the methodology's development. The results obtained confirm that the proposed procedure for simulating computed radiographic images is satisfactory, as it allows obtaining results comparable with experimental data. (author)
Adjoint Monte Carlo techniques and codes for organ dose calculations
Adjoint Monte Carlo simulations can be effectively used for the estimation of doses in small targets when the sources are extended in large volumes or surfaces. The main features of two computer codes for calculating doses at free points or in organs of an anthropomorphic phantom are described. In the first program (REBEL-3) natural gamma-emitting sources are contained in the walls of a dwelling room; in the second one (POKER-CAMP) the user can specify arbitrary gamma sources with different spatial distributions in the environment: in (or on the surface of) the ground and in the air. 3 figures
A code system has been developed to produce neutron cross section libraries for the MVP continuous energy Monte Carlo code from an evaluated nuclear data library in the ENDF format. The code system consists of 9 computer codes, and can process nuclear data in the latest ENDF-6 format. By using the present system, MVP neutron cross section libraries for important nuclides in reactor core analyses, shielding and fusion neutronics calculations have been prepared from JENDL-3.1, JENDL-3.2, JENDL-FUSION file and ENDF/B-VI data bases. This report describes the format of MVP neutron cross section library, the details of each code in the code system and how to use them. (author)
Highlights: → The specifications required for the analyses of the destructive assay data taken from irradiated fuel in Ohi-1 and Ohi-2 PWRs were documented in this paper. → These data were analyzed using the SWAT2.1 code, and the calculation results showed good agreement with experimental results. → These destructive assay data are suitable for the benchmarking of the burnup calculation code systems. - Abstract: The isotopic composition of spent nuclear fuels is vital data for studies on the nuclear fuel cycle and reactor physics. The Japan Atomic Energy Agency (JAEA) has been active in obtaining such data for pressurized water reactor (PWR) and boiling water reactor (BWR) fuels, and some data has already been published. These data have been registered with the international Spent Fuel Isotopic Composition Database (SFCOMPO) and widely used as international benchmarks for burnup calculation codes and libraries. In this paper, Assay Data of Spent Nuclear Fuel from two fuel assemblies irradiated in the Ohi-1 and Ohi-2 PWRs in Japan are shown. The destructive assay data from Ohi-2 have already been published. However, these data were not suitable for the benchmarking of calculation codes and libraries because several important specifications and data were not included. This paper summarizes the details of destructive assay data and specifications required for analyses of isotopic composition from Ohi-1 and Ohi-2. For precise burnup analyses, the burnup values of destructive assay samples were re-evaluated in this study. These destructive assay data were analyzed using the SWAT2.1 code, and the calculation results showed good agreement with experimental results. This indicates that the quality of destructive assay data from Ohi-1 and Ohi-2 PWRs is high, and that these destructive assay data are suitable for the benchmarking of burnup calculation code systems.
Lemehov, Sergei E; Suzuki, Motoe [Japan Atomic Energy Research Inst., Tokai, Ibaraki (Japan). Tokai Research Establishment
2001-08-01
PLUTON is a three-group neutronic code analyzing, as functions of time and burnup, the change of radial profiles, together with average values, of power density, burnup, concentration of trans-uranium elements, plutonium buildup, depletion of fissile elements, and fission product generation in water reactor fuel rod with standard UO{sub 2}, UO{sub 2}-Gd{sub 2}O{sub 3}, inhomogeneous MOX, and UO{sub 2}-ThO{sub 2}. The PLUTON code, which has been designed to be run on Windows PC, has adopted a theoretical shape function of neutron attenuation in pellet, which enables users to perform a very fast and accurate calculation easily. The present code includes the irradiation conditions of the Halden Reactor which gives verification data for the code. The total list of trans-uranium elements included in the calculations consists of {sub 92}U{sup 233-239}, {sub 93}Np{sup 237-239}, {sub 94}Pu{sup 238-243}, {sub 95}Am{sup 241-244} (including isomers), and {sub 96}Cm{sup 242-245}. Poisoning fission products are represented by {sub 54}Xe{sup 131,133,135}, {sub 48}Cd{sup 113}, {sub 62}Sm{sup 149,151,152}, {sub 64}Gd{sup 154-160}, {sub 63}Eu{sup 153,155}, {sub 36}Kr{sup 83,85}, {sub 42}Mo{sup 95}, {sub 43}Tc{sup 99}, {sub 45}Rh{sup 103}, {sub 47}Ag{sup 109}, {sub 53}I{sup 127,129,131}, {sub 55}Cs{sup 133}, {sub 57}La{sup 139}, {sub 59}Pr{sup 141}, {sub 60}Nd{sup 143-150}, {sub 61}Pm{sup 147}. Fission gases and volatiles included in the code are {sub 36}Kr{sup 83-86}, {sub 54}Xe{sup 129-136}, {sub 52}Te{sup 125-130}, {sub 53}I{sup 127-131}, {sub 55}Cs{sup 133-137}, and {sub 56}Ba{sup 135-140}. Verification has been performed up to 83 GWd/tU, and a satisfactory agreement has been obtained. (author)
Vectorization of continuous energy Monte Carlo code VIM
VIM is a continuous energy Monte Carlo code for criticality calculation. The random walk control system which uses combinatorial geometry system has been vectorized on FACOM VP-100. Vectorization has been done by the event bank method which controls simultaneous multiple particle's random walks, since behavior of neutron is independent. In vectorization of VIM code, we have two problems. One is a large overhead introduced by program modifications for vectorization. Another is a lowering of vector processing efficiency, since the vector length decreases with time according to the absorption and leakage of neutron and cut off of neutron for variance reduction. The average vector length during the random walks has been kept long by utilizing cross section library of single energy band and by reducing the number of the event banks. The performance ratio of vectorized version to the original one is 1.39 for the simple geometry and 1.13 for the complex geometry. (author)
Parallel computing by Monte Carlo codes MVP/GMVP
General-purpose Monte Carlo codes MVP/GMVP are well-vectorized and thus enable us to perform high-speed Monte Carlo calculations. In order to achieve more speedups, we parallelized the codes on the different types of parallel computing platforms or by using a standard parallelization library MPI. The platforms used for benchmark calculations are a distributed-memory vector-parallel computer Fujitsu VPP500, a distributed-memory massively parallel computer Intel paragon and a distributed-memory scalar-parallel computer Hitachi SR2201, IBM SP2. As mentioned generally, linear speedup could be obtained for large-scale problems but parallelization efficiency decreased as the batch size per a processing element(PE) was smaller. It was also found that the statistical uncertainty for assembly powers was less than 0.1% by the PWR full-core calculation with more than 10 million histories and it took about 1.5 hours by massively parallel computing. (author)
Parallelization of a Monte Carlo particle transport simulation code
Hadjidoukas, P.; Bousis, C.; Emfietzoglou, D.
2010-05-01
We have developed a high performance version of the Monte Carlo particle transport simulation code MC4. The original application code, developed in Visual Basic for Applications (VBA) for Microsoft Excel, was first rewritten in the C programming language for improving code portability. Several pseudo-random number generators have been also integrated and studied. The new MC4 version was then parallelized for shared and distributed-memory multiprocessor systems using the Message Passing Interface. Two parallel pseudo-random number generator libraries (SPRNG and DCMT) have been seamlessly integrated. The performance speedup of parallel MC4 has been studied on a variety of parallel computing architectures including an Intel Xeon server with 4 dual-core processors, a Sun cluster consisting of 16 nodes of 2 dual-core AMD Opteron processors and a 200 dual-processor HP cluster. For large problem size, which is limited only by the physical memory of the multiprocessor server, the speedup results are almost linear on all systems. We have validated the parallel implementation against the serial VBA and C implementations using the same random number generator. Our experimental results on the transport and energy loss of electrons in a water medium show that the serial and parallel codes are equivalent in accuracy. The present improvements allow for studying of higher particle energies with the use of more accurate physical models, and improve statistics as more particles tracks can be simulated in low response time.
A comparison between the Monte Carlo radiation transport codes MCNP and MCBEND
Sawamura, Hidenori; Nishimura, Kazuya [Computer Software Development Co., Ltd., Tokyo (Japan)
2001-01-01
In Japan, almost of all radiation analysts are using the MCNP code and MVP code on there studies. But these codes have not had automatic variance reduction. MCBEND code made by UKAEA have automatic variance reduction. And, MCBEND code is user friendly more than other Monte Carlo Radiation Transport Codes. Our company was first introduced MCBEND code in Japan. Therefore, we compared with MCBEND code and MCNP code about functions and production capacity. (author)
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
On the use of the Serpent Monte Carlo code for few-group cross section generation
Research highlights: → B1 methodology was used for generation of leakage-corrected few-group cross sections in the Serpent Monte-Carlo code. → Few-group constants generated by Serpent were compared with those calculated by Helios deterministic lattice transport code. → 3D analysis of a PWR core was performed by a nodal diffusion code DYN3D employing two-group cross section sets generated by Serpent and Helios. → An excellent agreement in the results of 3D core calculations obtained with Helios and Serpent generated cross-section libraries was observed. - Abstract: Serpent is a recently developed 3D continuous-energy Monte Carlo (MC) reactor physics burnup calculation code. Serpent is specifically designed for lattice physics applications including generation of homogenized few-group constants for full-core core simulators. Currently in Serpent, the few-group constants are obtained from the infinite-lattice calculations with zero neutron current at the outer boundary. In this study, in order to account for the non-physical infinite-lattice approximation, B1 methodology, routinely used by deterministic lattice transport codes, was considered for generation of leakage-corrected few-group cross sections in the Serpent code. A preliminary assessment of the applicability of the B1 methodology for generation of few-group constants in the Serpent code was carried out according to the following steps. Initially, the two-group constants generated by Serpent were compared with those calculated by Helios deterministic lattice transport code. Then, a 3D analysis of a Pressurized Water Reactor (PWR) core was performed by the nodal diffusion code DYN3D employing two-group cross section sets generated by Serpent and Helios. At this stage thermal-hydraulic (T-H) feedback was neglected. The DYN3D results were compared with those obtained from the 3D full core Serpent MC calculations. Finally, the full core DYN3D calculations were repeated taking into account T-H feedback and
Monte Carlo simulation of medical linear accelerator using primo code
The use of monte Carlo simulation has become very important in the medical field and especially in calculation in radiotherapy. Various Monte Carlo codes were developed simulating interactions of particles and photons with matter. One of these codes is PRIMO that performs simulation of radiation transport from the primary electron source of a linac to estimate the absorbed dose in a water phantom or computerized tomography (CT). PRIMO is based on Penelope Monte Carlo code. Measurements of 6 MV photon beam PDD and profile were done for Elekta precise linear accelerator at Radiation and Isotopes Center Khartoum using computerized Blue water phantom and CC13 Ionization Chamber. accept Software was used to control the phantom to measure and verify dose distribution. Elektalinac from the list of available linacs in PRIMO was tuned to model Elekta precise linear accelerator. Beam parameter of 6.0 MeV initial electron energy, 0.20 MeV FWHM, and 0.20 cm focal spot FWHM were used, and an error of 4% between calculated and measured curves was found. The buildup region Z max was 1.40 cm and homogenous profile in cross line and in line were acquired. A number of studies were done to verily the model usability one of them is the effect of the number of histories on accuracy of the simulation and the resulted profile for the same beam parameters. The effect was noticeable and inaccuracies in the profile were reduced by increasing the number of histories. Another study was the effect of Side-step errors on the calculated dose which was compared with the measured dose for the same setting.It was in range of 2% for 5 cm shift, but it was higher in the calculated dose because of the small difference between the tuned model and measured dose curves. Future developments include simulating asymmetrical fields, calculating the dose distribution in computerized tomographic (CT) volume, studying the effect of beam modifiers on beam profile for both electron and photon beams.(Author)
KAMCCO, a reactor physics Monte Carlo neutron transport code
KAMCCO is a 3-dimensional reactor Monte Carlo code for fast neutron physics problems. Two options are available for the solution of 1) the inhomogeneous time-dependent neutron transport equation (census time scheme), and 2) the homogeneous static neutron transport equation (generation cycle scheme). The user defines the desired output, e.g. estimates of reaction rates or neutron flux integrated over specified volumes in phase space and time intervals. Such primary quantities can be arbitrarily combined, also ratios of these quantities can be estimated with their errors. The Monte Carlo techniques are mostly analogue (exceptions: Importance sampling for collision processes, ELP/MELP, Russian roulette and splitting). Estimates are obtained from the collision and track length estimators. Elastic scattering takes into account first order anisotropy in the center of mass system. Inelastic scattering is processed via the evaporation model or via the excitation of discrete levels. For the calculation of cross sections, the energy is treated as a continuous variable. They are computed by a) linear interpolation, b) from optionally Doppler broadened single level Breit-Wigner resonances or c) from probability tables (in the region of statistically distributed resonances). (orig.)
Verification of Monte Carlo transport codes FLUKA, Mars and Shield
The present study is a continuation of the project 'Verification of Monte Carlo Transport Codes' which is running at GSI as a part of activation studies of FAIR relevant materials. It includes two parts: verification of stopping modules of FLUKA, MARS and SHIELD-A (with ATIMA stopping module) and verification of their isotope production modules. The first part is based on the measurements of energy deposition function of uranium ions in copper and stainless steel. The irradiation was done at 500 MeV/u and 950 MeV/u, the experiment was held at GSI from September 2004 until May 2005. The second part is based on gamma-activation studies of an aluminium target irradiated with an argon beam of 500 MeV/u in August 2009. Experimental depth profiling of the residual activity of the target is compared with the simulations. (authors)
Monte Carlo Code System Development for Liquid Metal Reactor
Kim, Chang Hyo; Shim, Hyung Jin; Han, Beom Seok; Park, Ho Jin; Park, Dong Gyu [Seoul National University, Seoul (Korea, Republic of)
2007-03-15
We have implemented the composition cell class and the use cell to MCCARD for hierarchy input processing. For the inputs of KALlMER-600 core consisted of 336 assemblies, we require the geometric data of 91,056 pin cells. Using hierarchy input processing, it was observed that the system geometries are correctly handled with the geometric data of total 611 cells; 2 cells for fuel rods, 2 cells for guide holes, 271 translation cells for rods, and 336 translation cells for assemblies. We have developed monte carlo decay-chain models based on decay chain model of REBUS code for liquid metal reactor analysis. Using developed decay-chain models, the depletion analysis calculations have performed for the homogeneous and heterogeneous model of KALlMER-600. The k-effective for the depletion analysis agrees well with that of REBUS code. and the developed decay chain models shows more efficient performance for time and memories, as compared with the existing decay chain model The chi-square criterion has been developed to diagnose the temperature convergence for the MC TjH feedback calculations. From the application results to the KALlMER pin and fuel assembly problem, it is observed that the new criterion works well Wc have applied the high efficiency variance reduction technique by splitting Russian roulette to estimate the PPPF of the KALIMER core at BOC. The PPPF of KALlMER core at BOC is 1.235({+-}0.008). The developed technique shows four time faster calculation, as compared with the existin2 calculation Subject Keywords Monte Carlo
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
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)
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)
A Monte Carlo track structure code for low energy protons
Endo, S; Nikjoo, H; Uehara, S; Hoshi, M; Ishikawa, M; Shizuma, K
2002-01-01
A code is described for simulation of protons (100 eV to 10 MeV) track structure in water vapor. The code simulates molecular interaction by interaction for the transport of primary ions and secondary electrons in the form of ionizations and excitations. When a low velocity ion collides with the atoms or molecules of a target, the ion may also capture or lose electrons. The probabilities for these processes are described by the quantity cross-section. Although proton track simulation at energies above Bragg peak (>0.3 MeV) has been achieved to a high degree of precision, simulations at energies near or below the Bragg peak have only been attempted recently because of the lack of relevant cross-section data. As the hydrogen atom has a different ionization cross-section from that of a proton, charge exchange processes need to be considered in order to calculate stopping power for low energy protons. In this paper, we have used state-of-the-art Monte Carlo track simulation techniques, in conjunction with the pub...
Monte Carlo simulation in UWB1 depletion code
UWB1 depletion code is being developed as a fast computational tool for the study of burnable absorbers in the University of West Bohemia in Pilsen, Czech Republic. In order to achieve higher precision, the newly developed code was extended by adding a Monte Carlo solver. Research of fuel depletion aims at development and introduction of advanced types of burnable absorbers in nuclear fuel. Burnable absorbers (BA) allow the compensation of the initial reactivity excess of nuclear fuel and result in an increase of fuel cycles lengths with higher enriched fuels. The paper describes the depletion calculations of VVER nuclear fuel doped with rare earth oxides as burnable absorber based on performed depletion calculations, rare earth oxides are divided into two equally numerous groups, suitable burnable absorbers and poisoning absorbers. According to residual poisoning and BA reactivity worth, rare earth oxides marked as suitable burnable absorbers are Nd, Sm, Eu, Gd, Dy, Ho and Er, while poisoning absorbers include Sc, La, Lu, Y, Ce, Pr and Tb. The presentation slides have been added to the article
The Monte Carlo code MCSHAPE: Main features and recent developments
MCSHAPE is a general purpose Monte Carlo code developed at the University of Bologna to simulate the diffusion of X- and gamma-ray photons with the special feature of describing the full evolution of the photon polarization state along the interactions with the target. The prevailing photon–matter interactions in the energy range 1–1000 keV, Compton and Rayleigh scattering and photoelectric effect, are considered. All the parameters that characterize the photon transport can be suitably defined: (i) the source intensity, (ii) its full polarization state as a function of energy, (iii) the number of collisions, and (iv) the energy interval and resolution of the simulation. It is possible to visualize the results for selected groups of interactions. MCSHAPE simulates the propagation in heterogeneous media of polarized photons (from synchrotron sources) or of partially polarized sources (from X-ray tubes). In this paper, the main features of MCSHAPE are illustrated with some examples and a comparison with experimental data. - Highlights: • MCSHAPE is an MC code for the simulation of the diffusion of photons in the matter. • It includes the proper description of the evolution of the photon polarization state. • The polarization state is described by means of the Stokes vector, I, Q, U, V. • MCSHAPE includes the computation of the detector influence in the measured spectrum. • MCSHAPE features are illustrated with examples and comparison with experiments
Inclusion of historical dependences of fuel burn-up into MOBY-DICK Code
The paper briefly describes inclusion of historical dependences of cross sections into MOBY-DICK code. Changes in program and its library are specified, especially from the point of view of programs user. Preliminary testing on assembly level is described and also testing on core level for 'transient' loading pattern calculations is demonstrated on examples of 18th and 19th cycles of the Dukovany NPP Unit III. Some features are addressed in the end of the paper (Authors)
Incorporation of the variation in conductivity with burnup in the stability of code predictive LAPUR
In the field of nuclear safety, the analysis of the stability of boiling water reactors is one of the biggest challenges for researchers. LAPUR code that allows to obtain the parameters of stability of the plant (Decay rate and frequency), being one of the programs used by IBERDROLA can be used for these calculations. With the collaboration of the research group TIN of the Polytechnic University of Valencia, a model of loss of conductivity of uranium has joined with the burned LAPUR. This update allows you to play the phenomenon in a more realistic way. This improvement has been validated and verified contrasting results with reference values.
Highlights: • The inventory of the radioactive nuclides was calculated using the GETERA code. • The consumptions of 235U for the HEU, MEU and LEU were: 2360, 2334 and 2320 g. • The amounts of 239Pu produced were: 67.07, 157.86 g for the MEU and LEU. • The core radioactivity for the MEU and LEU were: 8.84 × 1016 and 9.31 × 1016 Bq. - Abstract: Efforts have been made recently to study the possibility of core conversion of the 10 MW MTR type research reactor from the HEU to LEU fuels due to the proliferation issue. Since the inventory of the reactor core is a required parameter to study the atmospheric dispersion calculation for a postulated accidental airborne radionuclide release from the reactor, the inventory of the radioactive nuclides accumulated in the UAlx–Al fuels: HEU (93% 235U), MEU (45% 235U) and LEU (20% 235U) after 200 days of the reactor operating time was calculated using the GETERA code. The result showed, after 200 days of the reactor operation time (35% burnup), that the total consumptions of 235U and 238U for the HEU, MEU and LEU fuels were: 2360, 2334 and 2320 g for the 235U and 13, 105 and 238 g for the 238U, respectively. The amounts of 239Pu produced in the core were: 67.07, 157.86 g for the MEU and LEU fuels, respectively, compared with 7.95 g. The total core radioactivity after 200 days for the MEU and LEU cores were: 8.84 × 1016 and 9.31 × 1016 Bq, respectively, compared with 8.63 × 1016 Bq for the HEU core
Gholamzadeh Zohreh; Hossein Feghhi Seyed Amir; Soltani Leila; Rezazadeh Marzieh; Tenreiro Claudio; Joharifard Mahdi
2014-01-01
Decrease of the economically accessible uranium resources and the inherent proliferation resistance of thorium fuel motivate its application in nuclear power systems. Estimation of the nuclear reactor’s neutronic parameters during different operational situations is of key importance for the safe operation of nuclear reactors. In the present research, thorium oxide fuel burn-up calculations for a demonstrative model of a heavy water- -cooled reactor have been performed using MCNPX 2.6 code. N...
Monte Carlo N Particle code - Dose distribution of clinical electron beams in inhomogeneous phantoms
H A Nedaie; Mosleh-Shirazi, M. A.; Allahverdi, M.
2013-01-01
Electron dose distributions calculated using the currently available analytical methods can be associated with large uncertainties. The Monte Carlo method is the most accurate method for dose calculation in electron beams. Most of the clinical electron beam simulation studies have been performed using non- MCNP [Monte Carlo N Particle] codes. Given the differences between Monte Carlo codes, this work aims to evaluate the accuracy of MCNP4C-simulated electron dose distributions in a homogenous...
Parallel implementation of the Monte Carlo transport code EGS4 on the hypercube
Monte Carlo transport codes are commonly used in the study of particle interactions. The CALOR89 code system is a combination of several Monte Carlo transport and analysis programs. In order to produce good results, a typical Monte Carlo run will have to produce many particle histories. On a single processor computer, the transport calculation can take a huge amount of time. However, if the transport of particles were divided among several processors in a multiprocessor machine, the time can be drastically reduced
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
Method of tallying adjoint fluence and calculating kinetics parameters in Monte Carlo codes
A method of using iterated fission probability to estimate the adjoint fluence during particles simulation, and using it as the weighting function to calculate kinetics parameters βeff and A in Monte Carlo codes, was introduced in this paper. Implements of this method in continuous energy Monte Carlo code MCNP and multi-group Monte Carlo code MCMG are both elaborated. Verification results show that, with regardless additional computing cost, using this method, the adjoint fluence accounted by MCMG matches well with the result computed by ANISN, and the kinetics parameters calculated by MCNP agree very well with benchmarks. This method is proved to be reliable, and the function of calculating kinetics parameters in Monte Carlo codes is carried out effectively, which could be the basement for Monte Carlo codes' utility in the analysis of nuclear reactors' transient behavior. (authors)
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)
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.
Highlights: • Pu-239 based spectral history method was tested on 3D BWR single assembly case. • Burnup of a BWR fuel assembly was performed with the nodal code DYN3D. • Reference solution was obtained by coupled Monte-Carlo thermal-hydraulic code BGCore. • The proposed method accurately reproduces moderator density history effect for BWR test case. - Abstract: This research focuses on the verification of a recently developed methodology accounting for spectral history effects in 3D full core nodal simulations. The traditional deterministic core simulation procedure includes two stages: (1) generation of homogenized macroscopic cross section sets and (2) application of these sets to obtain a full 3D core solution with nodal codes. The standard approach adopts the branch methodology in which the branches represent all expected combinations of operational conditions as a function of burnup (main branch). The main branch is produced for constant, usually averaged, operating conditions (e.g. coolant density). As a result, the spectral history effects that associated with coolant density variation are not taken into account properly. Number of methods to solve this problem (such as micro-depletion and spectral indexes) were developed and implemented in modern nodal codes. Recently, we proposed a new and robust method to account for history effects. The methodology was implemented in DYN3D and involves modification of the few-group cross section sets. The method utilizes the local Pu-239 concentration as an indicator of spectral history. The method was verified for PWR and VVER applications. However, the spectrum variation in BWR core is more pronounced due to the stronger coolant density change. The purpose of the current work is investigating the applicability of the method to BWR analysis. The proposed methodology was verified against recently developed BGCore system, which couples Monte Carlo neutron transport with depletion and thermal-hydraulic solvers and
Burnup dependent core neutronic analysis for PBMR
The strategy for core neutronics modeling is based on SCALE4.4 code KENOV.a module that uses Monte Carlo calculational methods. The calculations are based on detailed unit cell and detailed core modeling. The fuel pebble is thoroughly modeled by introducing unit cell modeling for the graphite matrix and the fuel kernels in the pebble. The core is then modeled by placing these pebbles randomly throughout the core, yet not loosing track of any one of them. For the burnup model, a cyclic manner is adopted by coupling the KENOV.a and ORIGEN-S modules. Shifting down one slice at each discrete time step, and inserting fresh fuel from the top, this cyclic calculation model continues until equilibrium burnup cycle is achieved. (author)
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
Results of the RTOP-CA code calculations for experiments in the research MIR reactor are presented. The MIR-reactor tests were made to study the activity release from defective WWER fuel at high burnup (∼60 MWd/kgU). The RTOP-CA calculations are compared to experimental data on radial distributions of burnup as well as radial profiles of Pu and Xe concentrations in fuel pellets. The RTOP-CA predictions are also compared to the data on activity release (radionuclides of I, Cs, Xe and Kr) from the test fuel rod with an artificial defect in cladding. Additional calculations were performed for WWER-1000 fuel of an advanced design. In these calculation series the effect of design innovations on activity release from defective fuel rods was estimated. It is demonstrated that in case of a failure the new generation of WWER fuel shows lower levels of activity release into primary coolant. (authors)
A fuel performance code FEMAXI-6 has been developed for the analysis of LWR fuel rod behaviors in normal operation and transient (not accident) conditions. The code uses FEM for mechanical analysis, and has incorporated thermal and mechanical models of phenomena anticipated in high burn-up fuel rods, such as fuel thermal conductivity degradation and pellet-clad bonding. In the present study, PCMI induced by swelling in a high burn-up BWR type fuel rod has been analyzed by the FEMAXI-6 code. During a power ramp for the high burn-up fuel, instantaneous pellet swelling can significantly exceed the level that is predicted by a 'steady-rate' swelling model, causing a large circumferential strain in cladding. This phenomenon has been simulated by a new swelling model to take into account the fission gas bubble growth, and we found that the new model can give satisfactory predictions on cladding diametral expansion in comparison with measurements in test rod. The bubble growth model assumes an equilibrium between bubble size and external pressure on the bubble, and simultaneous solution is obtained with both bubble size determination and diffusion equation of fission gas atoms. In addition, a pellet-clad bonding model which assumes firm mechanical coupling between pellet outer surface and cladding inner surface predicts an elevated tensile stress in the axial direction of cladding during ramp, indicating the generation of bi-axial stress state in the cladding. These analyses by the FEMAXI-6 code enable us to predict the magnitude of mechanical loading on cladding during transient and also serve for failure investigation. Clearly, prediction by code calculation depends on the creep and stress-strain properties of highly irradiated cladding. (author)
Fast reactor 3D core and burnup analysis using VESTA
Luciano, N.; Shamblin, J.; Maldonado, I. [Nuclear Engineering Dept., Univ. of Tennessee, Knoxville, TN 37996-2300 (United States)
2012-07-01
Burnup analyses using the VESTA code have been performed on a MOX-fuelled fast reactor model as specified by an IAEA computational benchmark. VESTA is a relatively new code that has been used for burnup credit calculations and thermal reactor models, but not typically for fast reactor applications. The detailed input and results of the IAEA benchmark provides an opportunity to gauge the use of VESTA in a fast reactor application. VESTA employs an ultra-fine multi-group binning approach that accelerates Monte Carlo burnup calculations. Using VESTA to compute the end of cycle (EOC) power fractions by enrichment zone showed agreement with the published values within 5%. When comparing the ultra-fine multi-group binning approach to the tally-based approach, EOC isotopic masses also agree within 5%. Using the ultra-fine multi-group binning approach, we obtain a wall-time speedup factor of 35 when compared to the tally-based approach for computing a k{sub eff} eigenvalue with burnup problem. The authors conclude the use of VESTA's ultra-fine multi-group binning approach with Monte Carlo transport performs accurate depletion calculations for this fast reactor benchmark. (authors)
Verification of the Monte Carlo code RMC with a whole PWR MOX/UO2 core benchmark
Several types of V and V work are being carried out for the Reactor Monte Carlo code RMC, including the heterogeneous whole core configurations. In this paper, a whole PWR MOX/UO2 core benchmark which contains both UO2 and MOX assemblies with different enrichments and various burn-up points is chosen to verify RMC's criticality calculation capability, and the results of RMC and other codes are discussed and compared, such as eigenvalues, assembly power distributions, pin power distributions and so on. The discrepancies in eigenvalues and power distributions are satisfactory, which proves the accuracy of RMC's criticality calculation. Also, the influences of different cross-section libraries are discussed upon the results of RMC. Besides these results, the detailed comparisons between RMC and MCNP with the same ENDF/B-VII.0 cross-section library are carried out in this paper, including the comparisons of control rod worths calculated by both RMC and MCNP. According to the results, RMC and MCNP agree quite well in eigenvalues, power distributions and other results. The discrepancies of eigenvalues and control rod worth are fairly small and the relative differences of assembly and pin power distributions are acceptable. All these results contribute to the conclusion that the criticality calculation performance of RMC is accurate and excellent. (author)
Fuel burnup analysis for the Moroccan TRIGA research reactor
Highlights: ► A fuel burnup analysis of the 2 MW TRIGA MARK II Moroccan research reactor was established. ► Burnup calculations were done by means of the in-house developed burnup code BUCAL1. ► BUCAL1 uses the MCNP tallies directly in the calculation of the isotopic inventories. ► The reactor life time was found to be 3360 MW h considering full power operating conditions. ► Power factors and fluxes of the in-core irradiation positions are strongly affected by burnup. -- Abstract: The fundamental advantage and main reason to use Monte Carlo methods for burnup calculations is the possibility to generate extremely accurate burnup dependent one group cross-sections and neutron fluxes for arbitrary core and fuel geometries. Yet, a set of values determined for a material at a given position and time remains accurate only in a local region, in which neutron spectrum and flux vary weakly — and only for a limited period of time, during which changes of the local isotopic composition are minor. This paper presents the approach of fuel burnup evaluation used at the Moroccan TRIGA MARK II research reactor. The approach is essentially based upon the utilization of BUCAL1, an in-house developed burnup code. BUCAL1 is a FORTRAN computer code designed to aid in analysis, prediction, and optimization of fuel burnup performance in nuclear reactors. The code was developed to incorporate the neutron absorption reaction tally information generated directly by MCNP5 code in the calculation of fissioned or neutron-transmuted isotopes for multi-fueled regions. The fuel cycle length and changes in several core parameters such as: core excess reactivity, control rods position, fluxes at the irradiation positions, axial and radial power factors and other parameters are estimated. Besides, this study gives valuable insight into the behavior of the reactor and will ensure better utilization and operation of the reactor during its life-time and it will allow the establishment of
Supercell burnup model for the physics design of BWR fuel assemblies
A code called SUPERB has been developed for the BWR fuel assembly burnup analyses using supercell model. Each of the characteristic heterogeneities of a BWR fuel assembly like water gap, poisoned pins, control blade etc., is treated by invoking appropriate supercell concept. The burnup model of SUPERB is so devised as to strike a balance between accuracy and speed. This is achieved by building isotopic densities in each fuel pin separately while the depletion equations are solved only in a few groups of pins or burnup zones and the multigroup neutron spectra are differentiated in fewer group of pincell types. Multiple fuel ring burnup is considered only for Gd isotopes. A special empirical formula allows the microscopic cross section of Gd isotopes to be varied even during burnup integration. The supercell model has been tested against Monte Carlo results for the fresh cold clean Tarapur fuel assembly with two Gd fuel pins. The burnup model of SUPERB has been validated against one of the most sophisticated codes LWR-WIMS for a benchmark problem involving all the complexities of a BWR fuel assembly. The agreement of SUPERB results with both Monte Carlo and LWR-WIMS results is found to be excellent. (auth.)
Development of the SUHAM-U code for burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel is described. Developed SUHAM-U code has capacity to calculate burnup in each fuel or poison zone of each cell of VVER-1000 fuel assembly. In so doing Surface Harmonics Method is used for calculation of the detail neutron spectra in fuel assembly at separated burnup values. Verification of SUHAM-U code by burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel has been carried out. Comparisons were carried out with calculations by UNK and RECOL codes. UNK code uses the first collisions probabilities method for solution of the neutron transport equation and RECOL code uses Monte-Carlo method with point-wise continues energy presentation of cross-sections. The main conclusion of all comparisons is the SUHAM-U code calculates the fuel burnup of VVER-1000 fuel assemblies with uranium and MOX fuel with enough high accuracy. Time expenditures are adduced. (authors)
Boyarinov, V. F.; Davidenko, V. D.; Polismakov, A. A.; Tsibulsky, V. F. [Russian Research Center Kurchatov Inst., Nuclear Reactor Inst., 123182, Moscow (Russian Federation)
2006-07-01
Development of the SUHAM-U code for burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel is described. Developed SUHAM-U code has capacity to calculate burnup in each fuel or poison zone of each cell of VVER-1000 fuel assembly. In so doing Surface Harmonics Method is used for calculation of the detail neutron spectra in fuel assembly at separated burnup values. Verification of SUHAM-U code by burnup calculations of VVER-1000 fuel assemblies with uranium and MOX fuel has been carried out. Comparisons were carried out with calculations by UNK and RECOL codes. UNK code uses the first collisions probabilities method for solution of the neutron transport equation and RECOL code uses Monte-Carlo method with point-wise continues energy presentation of cross-sections. The main conclusion of all comparisons is the SUHAM-U code calculates the fuel burnup of VVER-1000 fuel assemblies with uranium and MOX fuel with enough high accuracy. Time expenditures are adduced. (authors)
On the inner workings of Monte Carlo codes
Dubbeldam, D.; Torres Knoop, A.; Walton, K.S.
2013-01-01
We review state-of-the-art Monte Carlo (MC) techniques for computing fluid coexistence properties (Gibbs simulations) and adsorption simulations in nanoporous materials such as zeolites and metal-organic frameworks. Conventional MC is discussed and compared to advanced techniques such as reactive MC, configurational-bias Monte Carlo and continuous fractional MC. The latter technique overcomes the problem of low insertion probabilities in open systems. Other modern methods are (hyper-)parallel...
Proton therapy Monte Carlo SRNA-VOX code
Ilić Radovan D.
2012-01-01
The most powerful feature of the Monte Carlo method is the possibility of simulating all individual particle interactions in three dimensions and performing numerical experiments with a preset error. These facts were the motivation behind the development of a general-purpose Monte Carlo SRNA program for proton transport simulation in technical systems described by standard geometrical forms (plane, sphere, cone, cylinder, cube). Some of the possible applications of the SRNA program are:...
The influence of pitch, burnup and absorber rods on the spent fuel pool criticality
It has been shown that supercriticality might occur for some postulated accidents for the TRIGA spent fuel pool at ''Josef Stefan'' Institute in Ljubljana, Slovenia. However, in the previous studies, the effect of burnup was not accounted for. In this work the dependence of criticality on fuel burnup, the pitch among the elements and the number of uniformly mixed absorber rods for a square arrangement is presented. The Monte Carlo computer code MCNP4B with ENDF-B/VI library and detailed three dimensional geometry was used. WIMS-D code was used to model the isotopic composition of the fuel for 5, 10, 20 and 30 % burnup without cooling time. The results show, that out of the three studied effects: pitch from contact (3.75cm) up to rack design pitch (8cm), number of absorbers from 0 to 8 and burnup up to 30 %, the pitch has the greatest influence on the multiplication factor keff. In the interval in which the pitch was changed, keff decreased for up to 0.45. The number of absorber rods affects the multiplication factor much less. This effect is bigger for more compact arrangements, e.g. for contact of fuel elements with 8 absorber rods among them, keff values are smaller for almost 0.20 than for arrangement without absorber rods regardless of the burnup. The effect of burnup is the smallest since in no case keff decreases for more than 0.10, even for high burnups of 30 %. (author)
Recent developments of JAEA’s Monte Carlo code MVP for reactor physics applications
Highlights: • This paper describes the recent development status of the Monte Carlo code MVP. • The basic features and capabilities of MVP are briefly described. • New capabilities useful for reactor analysis are also described. - Abstract: This paper describes the recent development status of a Monte Carlo code MVP developed at Japan Atomic Energy Agency. The basic features and capabilities of MVP are overviewed. In addition, new capabilities useful for reactor analysis are also described
In the report, research results discussed in 1999 fiscal year at Nuclear Code Evaluation Committee of Nuclear Code Research Committee were summarized. Present status of Monte Carlo simulation on nuclear energy study was described. Especially, besides of criticality, shielding and core analyses, present status of applications to risk and radiation damage analyses, high energy transport and nuclear theory calculations of Monte Carlo Method was described. The 18 papers are indexed individually. (J.P.N.)
Development of a Monte-Carlo Radiative Transfer Code for the Juno/JIRAM Limb Measurements
Sindoni, G.; Adriani, A.; Mayorov, B.; Aoki, S.; Grassi, D.; Moriconi, M.; Oliva, F.
2013-09-01
The Juno/JIRAM instrument will acquire limb spectra of the Jupiter atmosphere in the infrared spectral range. The analysis of these spectra requires a radiative transfer code that takes into account the multiple scattering by particles in a spherical-shell atmosphere. Therefore, we are developing a code based on the Monte-Carlo approach to simulate the JIRAM observations. The validation of the code was performed by comparison with DISORT-based codes.
Monte Carlo capabilities of the SCALE code system
Highlights: • Foundational Monte Carlo capabilities of SCALE are described. • Improvements in continuous-energy treatments are detailed. • New methods for problem-dependent temperature corrections are described. • New methods for sensitivity analysis and depletion are described. • Nuclear data, users interfaces, and quality assurance activities are summarized. - Abstract: SCALE is a widely used suite of tools for nuclear systems modeling and simulation that provides comprehensive, verified and validated, user-friendly capabilities for criticality safety, reactor physics, radiation shielding, and sensitivity and uncertainty analysis. For more than 30 years, regulators, licensees, and research institutions around the world have used SCALE for nuclear safety analysis and design. SCALE provides a “plug-and-play” framework that includes three deterministic and three Monte Carlo radiation transport solvers that can be selected based on the desired solution, including hybrid deterministic/Monte Carlo simulations. SCALE includes the latest nuclear data libraries for continuous-energy and multigroup radiation transport as well as activation, depletion, and decay calculations. SCALE’s graphical user interfaces assist with accurate system modeling, visualization, and convenient access to desired results. SCALE 6.2 will provide several new capabilities and significant improvements in many existing features, especially with expanded continuous-energy Monte Carlo capabilities for criticality safety, shielding, depletion, and sensitivity and uncertainty analysis. An overview of the Monte Carlo capabilities of SCALE is provided here, with emphasis on new features for SCALE 6.2
Application of Monte Carlo code EGS4 to calculate gamma exposure buildup factors
Exposure buildup factors up to 40 mean free paths ranging from 0.015 MeV to 15 MeV photon energy were calculated by using the Monte Carlo simulation code EGS4 for ordinary concrete. The calculation involves PHOTX cross section library, a point isotropic source, infinite uniform medium model and a particle splitting method and considers the Bremsstrahlung, fluorescent effect, correlative (Rayleigh) scatter. The results were compared with the relevant data. Results show that the data of the buildup factors calculated by the Monte Carlo code EGS4 was reliable. The Monte Carlo method can be used widely to calculate gamma-ray exposure buildup factors. (authors)
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)
The analog linear interpolation approach for Monte Carlo simulation of PGNAA: The CEARPGA code
Zhang, Wenchao; Gardner, Robin P.
2004-01-01
The analog linear interpolation approach (ALI) has been developed and implemented to eliminate the big weight problem in the Monte Carlo simulation code CEARPGA. The CEARPGA code was previously developed to generate elemental library spectra for using the Monte Carlo - library least-squares (MCLLS) approach in prompt gamma-ray neutron activation analysis (PGNAA). In addition, some other improvements to this code have been introduced, including (1) adopting the latest photon cross-section data, (2) using an improved detector response function, (3) adding the neutron activation backgrounds, (4) generating the individual natural background libraries, (5) adding the tracking of annihilation photons from pair production interactions outside of the detector and (6) adopting a general geometry package. The simulated result from the new CEARPGA code is compared with those calculated from the previous CEARPGA code and the MCNP code and experimental data. The new CEARPGA code is found to give the best result.
On the inner workings of Monte Carlo codes
D. Dubbeldam; A. Torres Knoop; K.S. Walton
2013-01-01
We review state-of-the-art Monte Carlo (MC) techniques for computing fluid coexistence properties (Gibbs simulations) and adsorption simulations in nanoporous materials such as zeolites and metal-organic frameworks. Conventional MC is discussed and compared to advanced techniques such as reactive MC
Data libraries as a collaborative tool across Monte Carlo codes
Augelli, Mauro; Han, Mincheol; Hauf, Steffen; Kim, Chan-Hyeung; Kuster, Markus; Pia, Maria Grazia; Quintieri, Lina; Saracco, Paolo; Seo, Hee; Sudhakar, Manju; Eidenspointner, Georg; Zoglauer, Andreas
2010-01-01
The role of data libraries in Monte Carlo simulation is discussed. A number of data libraries currently in preparation are reviewed; their data are critically examined with respect to the state-of-the-art in the respective fields. Extensive tests with respect to experimental data have been performed for the validation of their content.
Monte Carlo Capabilities of the SCALE Code System
Rearden, B. T.; Petrie, L. M.; Peplow, D. E.; Bekar, K. B.; Wiarda, D.; Celik, C.; Perfetti, C. M.; Ibrahim, A. M.; Hart, S. W. D.; Dunn, M. E.
2014-06-01
SCALE is a widely used suite of tools for nuclear systems modeling and simulation that provides comprehensive, verified and validated, user-friendly capabilities for criticality safety, reactor physics, radiation shielding, and sensitivity and uncertainty analysis. For more than 30 years, regulators, licensees, and research institutions around the world have used SCALE for nuclear safety analysis and design. SCALE provides a "plug-and-play" framework that includes three deterministic and three Monte Carlo radiation transport solvers that can be selected based on the desired solution, including hybrid deterministic/Monte Carlo simulations. SCALE includes the latest nuclear data libraries for continuous-energy and multigroup radiation transport as well as activation, depletion, and decay calculations. SCALE's graphical user interfaces assist with accurate system modeling, visualization, and convenient access to desired results. SCALE 6.2, to be released in 2014, will provide several new capabilities and significant improvements in many existing features, especially with expanded continuous-energy Monte Carlo capabilities for criticality safety, shielding, depletion, and sensitivity and uncertainty analysis. An overview of the Monte Carlo capabilities of SCALE is provided here, with emphasis on new features for SCALE 6.2.
Aurora T: a Monte Carlo code for transportation of neutral atoms in a toroidal plasma
This paper contains a short description of Aurora code. This code have been developed at Princeton with Monte Carlo method for calculating neutral gas in cylindrical plasma. In this work subroutines such one can take in account toroidal geometry are developed
MCNP, a general Monte Carlo code for neutron and photon transport: a summary
The general-purpose Monte Carlo code MCNP can be used for neutron, photon, or coupled neutron-photon transport, including the capability to calculate eigenvalues for critical systems. The code treats an arbitrary three-dimensional configuration of materials in geometric cells bounded by first- and second-degree surfaces and some special fourth-degree surfaces
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
MCNP: a general Monte Carlo code for neutron and photon transport. Version 3A. Revision 2
This manual is a practical guide for the use of our general-purpose Monte Carlo code MCNP. The first chapter is a primer for the novice user. The second chapter describes the mathematics, data, physics, and Monte Carlo simulation found in MCNP. This discussion is not meant to be exhaustive - details of the particular techniques and of the Monte Carlo method itself will have to be found elsewhere. The third chapter shows the user how to prepare input for the code. The fourth chapter contains several examples, and the fifth chapter explains the output. The appendices show how to use MCNP on particular computer systems at the Los Alamos National Laboratory and also give details about some of the code internals that those who wish to modify the code may find useful. 57 refs
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.
Using deterministic codes to accelerate continuous energy Monte-Carlo standards calculations
Deterministic codes are usually used for critical parameters or one dimension geometry calculations. Advantages of the use of deterministic codes are speed of the calculation and the absence of standard deviation on the keff results. Nevertheless, the deterministic results are affected by several intrinsic uncertainties as energetic condensation or self-shielding. So the way to proceed at CEA expert criticality group (CEA/SERMA/CP2C) is to always check the main results (minimum critical or maximal permissible values and un-moderated values) with a punctual Monte Carlo calculation. These last years, in particular cases (pure actinide fissile media, exotic reflectors), large discrepancies have been observed between the keff calculated by the CRISTAL V1 route reference (continuous energy Monte Carlo code TRIPOLI-4) and the keff target (by the standard route APOLLO2-Sn). The problematic for these cases was how to transpose the keff discrepancies observed between standard and reference routes to the dimensions (mass, thickness...) or how to reduce the keff discrepancies using optimized options of the deterministic code. One solution to transpose discrepancies is to iterate on dimensions using a punctual Monte Carlo code to achieve the desired keff eigenvalue. But, the amount of time for obtaining a good standard deviation and also the desired keff eigenvalue inside the Monte Carlo calculation uncertainty can quickly increase. The principle of the method presented in this paper is that the discrepancy between deterministic code and Monte-Carlo code, calculated at the same dimension, is low variable with the dimension. Therefore, correcting the keff eigenvalue on which the deterministic code converge with the discrepancy observed, leads to a dimension nearer to the true dimension (i.e. the dimension where Monte-Carlo code keff calculation is close to the keff eigenvalue). If the keff eigenvalue is outside the Monte Carlo uncertainty, the discrepancy is recalculated and
Calculations of neutron penetration through graphite medium with Monte Carlo code MCNP
Experiments for fast neutron penetration through graphite are analysed with the continuous energy Monte Carlo code MCNP. Reaction rates and energy spectra obtained with the MCNP are compared with measured values and calculated ones with McBEND code. And validity of penetration calculation with the MCNP is comfirmed. In addition, it is revealed that the MCNP code using Weight-Window method is well applicable to calculations of neutron penetration through graphite up to 70 cm in depth. (author)
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
Extension of the TRANSURANUS burn-up model
The validation range of the model in the TRANSURANUS fuel performance code for calculating the radial power density and burn-up in UO2 fuel has been extended from 64 MWd/kgHM up to 102 MWd/kgHM, thereby improving also its precision. In addition, the first verification of calculations with post-irradiation examination data is reported for LWR-MOX fuel with a rod average burn-up up to 45 MWd/kgHM. The extension covers the inclusion of new isotopes in order to account for the production of 238Pu. The corresponding one-group cross-sections used in the equations rely on results obtained with ALEPH, a new Monte Carlo burn-up code. The experimental verification is based on electron probe microanalysis (EPMA) and on secondary ion mass spectrometry (SIMS) as well as radiochemical data of fuel irradiated in commercial power plants. The deviations are quantified in terms of frequency distributions of the relative errors. The relative errors on the burn-up distributions in both fuel types remain below 12%, corresponding to the experimental scatter
MCMG: a 3-D multigroup P3 Monte Carlo code and its benchmarks
In this paper a 3-D Monte Carlo multigroup neutron transport code MCMG has been developed from a coupled neutron and photon transport Monte Carlo code MCNP. The continuous-energy cross section library of the MCNP code is replaced by the multigroup cross section data generated by the transport lattice code, such as the WIMS code. It maintains the strong abilities of MCNP for geometry treatment, counting, variance reduction techniques and plotting. The multigroup neutron scattering cross sections adopt the Pn (n ≤ 3) approximation. The test results are in good agreement with the results of other methods and experiments. The number of energy groups can be varied from few groups to multigroup, and either macroscopic or microscopic cross section can be used. (author)
A vectorized Monte Carlo code for modeling photon transport in SPECT
A vectorized Monte Carlo computer code has been developed for modeling photon transport in single photon emission computed tomography (SPECT). The code models photon transport in a uniform attenuating region and photon detection by a gamma camera. It is adapted from a history-based Monte Carlo code in which photon history data are stored in scalar variables and photon histories are computed sequentially. The vectorized code is written in FORTRAN77 and uses an event-based algorithm in which photon history data are stored in arrays and photon history computations are performed within DO loops. The indices of the DO loops range over the number of photon histories, and these loops may take advantage of the vector processing unit of our Stellar GS1000 computer for pipelined computations. Without the use of the vector processor the event-based code is faster than the history-based code because of numerical optimization performed during conversion to the event-based algorithm. When only the detection of unscattered photons is modeled, the event-based code executes 5.1 times faster with the use of the vector processor than without; when the detection of scattered and unscattered photons is modeled the speed increase is a factor of 2.9. Vectorization is a valuable way to increase the performance of Monte Carlo code for modeling photon transport in SPECT
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.
Accelerating execution of the integrated TIGER series Monte Carlo radiation transport codes
Execution of the integrated TIGER series (ITS) of coupled electron/photon Monte Carlo radiation transport codes has been accelerated by modifying the FORTRAN source code for more efficient computation. Each member code of ITS was benchmarked and profiled with a specific test case that directed the acceleration effort toward the most computationally intensive subroutines. Techniques for accelerating these subroutines included replacing linear search algorithms with binary versions, replacing the pseudo-random number generator, reducing program memory allocation, and proofing the input files for geometrical redundancies. All techniques produced identical or statistically similar results to the original code. Final benchmark timing of the accelerated code resulted in speed-up factors of 2.00 for TIGER (the one-dimensional slab geometry code), 1.74 for CYLTRAN (the two-dimensional cylindrical geometry code), and 1.90 for ACCEPT (the arbitrary three-dimensional geometry code)
Tokyo Metropolitan University of Health Sciences has done The Information Education using EGS4 Monte Carlo code since the 1998 fiscal year. Two items under practical training item were done. 1. The interaction between photon of 0.1 ∼ 10 MeV (Mega Electron Volt: MeV) and Aluminum (Al), Iron (Fe) and Lead (Pb). 2. The simulation of gamma ray energy measurement of the radiation detector. As the result, the student was possible the understanding of the radiation physics for the easiness at Practical training of EGS4 Monte Carlo code. (author)
The Monte Carlo code MCBEND - where it is and where it's going
The Monte Carlo method forms a corner stone to the calculational procedures established in the UK for shielding design and assessment. The emphasis of the work in the shielding area is centred on the Monte Carlo code MCBEND. The work programme in support of the code is broadly directed towards utilisation of new hardware, the development of improved modelling algorithms, the development of new acceleration methods for specific applications and enhancements to user image. This paper summarises the current status of MCBEND and reviews developments carried out over the past two years and planned for the future. (author)
Parallelization of MCATNP MONTE CARLO particle transport code by using MPI
A Monte Carlo code for simulating Atmospheric Transport of Neutrons and Photons (MCATNP) is used to simulate the ionization effects caused by high altitude nuclear detonation (HAND) and it was parallelized in MPI by adopting the leap random number producer and modifying the original serial code. The parallel results and serial results are identical. The speedup increases almost linearly with the number of processors used. The parallel efficiency is up to to 97% while 16 processors are used, and 94% while 32 are used. The experimental results show that parallelization can obviously reduce the calculation time of Monte Carlo simulation of HAND ionization effects. (authors)
Generalized Albedo option on the Morse Monte Carlo code
The advisability of using the albedo procedure for solving deep penetration shielding problems which have ducts and other penetrations is investigated. It is generally accepted that the use of albedo data can dramatically improve the computational efficiency of certain Monte Carlo calculations - however the accuracy of these results may be unacceptable because of lost information during the albedo event and serious errors in the available differential albedo data. This study has been done to evaluate and appropriately modify the MORSE/BREESE package, to develop new methods for generating the required albedo data, and to extend the adjoint capability to the albedo modified calculations. The major modifications include the tracking of special particles inside albedo media, an option to displace the point-of-emergence during an albedo event, and an option to read, process, and use spatially-dependent albedo data for both forward and adjoint calculations. (author)
The three-dimensional Monte-Carlo code TRIPOLI-02
TRIPOLI-2 solves the transport equation for neutrons or gamma rays in tridimensional geometrical configurations. TRIPOLI uses the Monte Carlo method. This method allows to treat exactly the geometrical configurations, the energy losses and the scattering laws. TRIPOLI 2 allows to treat the following problems: gamma transport problems, neutrons transport problems with fixed source (the problems can be time dependent or not), critical problems without fixed source and research of multiplication factor due to fissions, subcritical problems with fixed source and with multiplication by fission. These problems can be separate in two types. First type: shielding problems essentially with deep penetration and streaming through voids. Biasing technics are used to reduce the computing time. Second type: core problems for cell calculations or for small core calculations. In this case, it is necessary to have a fine representation of the cross sections. The thermalization is also treated exactly
Importance function by collision probabilities for Monte Carlo code Tripoli
We present a completely automatic biasing technique where the parameters of the biased simulation are deduced from the solution of the adjoint transport equation calculated by collision probabilities. In this study we shall estimate the importance function through collision probabilities method and we shall evaluate its possibilities thanks to a Monte Carlo calculation. We have run simulations with this new biasing method for one-group transport problems with isotropic shocks (one dimension geometry and X-Y geometry) and for multigroup problems with anisotropic shocks (one dimension geometry). For the anisotropic problems we solve the adjoint equation with anisotropic collision probabilities. The results show that for the one-group and homogeneous geometry transport problems the method is quite optimal without Splitting and Russian Roulette technique but for the multigroup and heterogeneous X-Y geometry ones the figures of merit are higher if we add Splitting and Russian Roulette technique
Longitudinal development of extensive air showers: hybrid code SENECA and full Monte Carlo
Ortiz, J A; De Souza, V; Ortiz, Jeferson A.; Tanco, Gustavo Medina
2004-01-01
New experiments, exploring the ultra-high energy tail of the cosmic ray spectrum with unprecedented detail, are exerting a severe pressure on extensive air hower modeling. Detailed fast codes are in need in order to extract and understand the richness of information now available. Some hybrid simulation codes have been proposed recently to this effect (e.g., the combination of the traditional Monte Carlo scheme and system of cascade equations or pre-simulated air showers). In this context, we explore the potential of SENECA, an efficient hybrid tridimensional simulation code, as a valid practical alternative to full Monte Carlo simulations of extensive air showers generated by ultra-high energy cosmic rays. We extensively compare hybrid method with the traditional, but time consuming, full Monte Carlo code CORSIKA which is the de facto standard in the field. The hybrid scheme of the SENECA code is based on the simulation of each particle with the traditional Monte Carlo method at two steps of the shower devel...
Depletion of a BWR lattice using the racer continuous energy Monte Carlo code
In the past several years there has been a renewed interest in the accuracy of a new generation of lattice physics codes. Most of the time these codes are benchmarked against Monte Carlo codes only at beginning of cycle. In this paper a highly heterogeneous BWR lattice depletion benchmark problem is presented. Results of a 40% void depletion using the RACER continuous energy Monte Carlo code are also presented. Complete problem specifications are given so that comparisons with lattice physics codes or other Monte Carlo codes is possible. The RACER calculations were performed with the ENDF/B-V cross section set. Each flux calculation utilized 2.7 million histories resulting in 95% confidence intervals of ∼1 milli-k on the eigenvalue and ∼1% uncertainties on pin-wise power fractions. Timing statistics for the calculation using the vectorized RACER code averaged ∼ 24,000 neutrons/minute on a single processor of a CRAY-C90 computer
MKENO-DAR: a direct angular representation Monte Carlo code for criticality safety analysis
Improving the Monte Carlo code MULTI-KENO, the MKENO-DAR (Direct Angular Representation) code has been developed for criticality safety analysis in detail. A function was added to MULTI-KENO for representing anisotropic scattering strictly. With this function, the scattering angle of neutron is determined not by the average scattering angle μ-bar of the Pl Legendre polynomial but by the random work operation using probability distribution function produced with the higher order Legendre polynomials. This code is avilable for the FACOM-M380 computer. This report is a computer code manual for MKENO-DAR. (author)
Full-core pin-power calculations using Monte Carlo codes
Pin wise calculations of core power distribution have been performed for a criticality mock up installation that models a WWER-1000 reactor. Two Monte Carlo codes have been applied for solving of this problem: the MCNP4B code and the KENO-VI code from the SCALE 4.4 system. The codes use different kinds of neutron cross section data: pointwise continuous-energy ENDF/B-VI data and multigroup ENDF/B-V data. Comparisons of calculated results show that the MCNP4B and KENO-VI results are in good agreement. (authors)
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)
Parallel processing of Monte Carlo code MCNP for particle transport problem
Higuchi, Kenji; Kawasaki, Takuji
1996-06-01
It is possible to vectorize or parallelize Monte Carlo codes (MC code) for photon and neutron transport problem, making use of independency of the calculation for each particle. Applicability of existing MC code to parallel processing is mentioned. As for parallel computer, we have used both vector-parallel processor and scalar-parallel processor in performance evaluation. We have made (i) vector-parallel processing of MCNP code on Monte Carlo machine Monte-4 with four vector processors, (ii) parallel processing on Paragon XP/S with 256 processors. In this report we describe the methodology and results for parallel processing on two types of parallel or distributed memory computers. In addition, we mention the evaluation of parallel programming environments for parallel computers used in the present work as a part of the work developing STA (Seamless Thinking Aid) Basic Software. (author)
MCAM 5: an advanced interface program for multiple Monte Carlo Codes
The Automatic Modeling Program for Neutronics and Radiation Transport Simulation (MCAM) developed in China, is an advanced interface program between CAD (Computer Aided Design) systems and Monte Carlo (MC) codes. It can significantly reduce the manpower and enhance reliability for constructing MC codes input of complex systems. The latest version MCAM 4.8 was a mature and efficient version which was benchmarked with ITER benchmark model and has been used by hundreds of institutes in more than 40 countries all over the world. It can deal with MCNP and TRIPOLI models. The main function of MCAM is to convert geometries in CAD systems to geometries in MC codes input files. The MCAM version 5.2 is going to be released with added capabilities to support SuperMC, Geant4 and FLUKA Monte Carlo codes
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
Modelling photon transport in non-uniform media for SPECT with a vectorized Monte Carlo code.
Smith, M F
1993-10-01
A vectorized Monte Carlo code has been developed for modelling photon transport in non-uniform media for single-photon-emission computed tomography (SPECT). The code is designed to compute photon detection kernels, which are used to build system matrices for simulating SPECT projection data acquisition and for use in matrix-based image reconstruction. Non-uniform attenuating and scattering regions are constructed from simple three-dimensional geometric shapes, in which the density and mass attenuation coefficients are individually specified. On a Stellar GS1000 computer, Monte Carlo simulations are performed between 1.6 and 2.0 times faster when the vector processor is utilized than when computations are performed in scalar mode. Projection data acquired with a clinical SPECT gamma camera for a line source in a non-uniform thorax phantom are well modelled by Monte Carlo simulations. The vectorized Monte Carlo code was used to stimulate a 99Tcm SPECT myocardial perfusion study, and compensations for non-uniform attenuation and the detection of scattered photons improve activity estimation. The speed increase due to vectorization makes Monte Carlo simulation more attractive as a tool for modelling photon transport in non-uniform media for SPECT. PMID:8248288
Systems guide to MCNP (Monte Carlo Neutron and Photon Transport Code)
The subject of this report is the implementation of the Los Alamos National Laboratory Monte Carlo Neutron and Photon Transport Code - Version 3 (MCNP) on the different types of computer systems, especially the IBM MVS system. The report supplements the documentation of the RSIC computer code package CCC-200/MCNP. Details of the procedure to follow in executing MCNP on the IBM computers, either in batch mode or interactive mode, are provided
Perfetti, Christopher M [ORNL; Martin, William R [University of Michigan; Rearden, Bradley T [ORNL; Williams, Mark L [ORNL
2012-01-01
Three methods for calculating continuous-energy eigenvalue sensitivity coefficients were developed and implemented into the SHIFT Monte Carlo code within the Scale code package. The methods were used for several simple test problems and were evaluated in terms of speed, accuracy, efficiency, and memory requirements. A promising new method for calculating eigenvalue sensitivity coefficients, known as the CLUTCH method, was developed and produced accurate sensitivity coefficients with figures of merit that were several orders of magnitude larger than those from existing methods.
QCDMPI - pure QCD Monte Carlo simulation code with MPI
QCDMPI is a pure QCD simulation code with MPI calls. QCDMPI is very portable because; - you can simulate any-dimensional QCD, - on any-dimensional partitioning, - on any number of processors, - with rather small working area. Also by this program, you can get two performances, - calculation (link update time) - communication (MB/sec). In this paper, outline of QCDMPI is reported. Comparison of the performances on several parallel machines; AP1000, AP1000+, AP3000, Cenju-3, Paragon, SR2201 and Workstation Cluster, is also reported. (orig.)
Verification of the shift Monte Carlo code with the C5G7 reactor benchmark
Sly, N. C.; Mervin, B. T. [Dept. of Nuclear Engineering, Univ. of Tennessee, 311 Pasqua Engineering Building, Knoxville, TN 37996-2300 (United States); Mosher, S. W.; Evans, T. M.; Wagner, J. C. [Oak Ridge National Laboratory, PO Box 2008, Oak Ridge, TN 37831 (United States); Maldonado, G. I. [Dept. of Nuclear Engineering, Univ. of Tennessee, 311 Pasqua Engineering Building, Knoxville, TN 37996-2300 (United States)
2012-07-01
Shift is a new hybrid Monte Carlo/deterministic radiation transport code being developed at Oak Ridge National Laboratory. At its current stage of development, Shift includes a parallel Monte Carlo capability for simulating eigenvalue and fixed-source multigroup transport problems. This paper focuses on recent efforts to verify Shift's Monte Carlo component using the two-dimensional and three-dimensional C5G7 NEA benchmark problems. Comparisons were made between the benchmark eigenvalues and those output by the Shift code. In addition, mesh-based scalar flux tally results generated by Shift were compared to those obtained using MCNP5 on an identical model and tally grid. The Shift-generated eigenvalues were within three standard deviations of the benchmark and MCNP5-1.60 values in all cases. The flux tallies generated by Shift were found to be in very good agreement with those from MCNP. (authors)
The general purpose code BALTORO was written for coupling the three-dimensional Monte-Carlo /MC/ with the one-dimensional Discrete Ordinates /DO/ radiation transport calculations. The quantity of a radiation-induced /neutrons or gamma-rays/ nuclear effect or the score from a radiation-yielding nuclear effect can be analysed in this way. (author)
The present of shielding analysis with nuclear data for continuous energy Monte Carlo code MCNP
Following three problems are analyzed by continuous energy Monte Carlo code MCNP with JENDL-3.2, 3.3, and ENDF/B-VI. 1. Shielding analysis of WINFRITH-Aspins iron deep penetration experiment. 2. Shielding analysis of TN-12A spent fuel transport cask experiment. 3. Shielding analysis of modular shielding house keeping spent fuel transportable casks. (author)
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)
Subroutines to Simulate Fission Neutrons for Monte Carlo Transport Codes
Lestone, J P
2014-01-01
Fortran subroutines have been written to simulate the production of fission neutrons from the spontaneous fission of 252Cf and 240Pu, and from the thermal neutron induced fission of 239Pu and 235U. The names of these four subroutines are getnv252, getnv240, getnv239, and getnv235, respectively. These subroutines reproduce measured first, second, and third moments of the neutron multiplicity distributions, measured neutron-fission correlation data for the spontaneous fission of 252Cf, and measured neutron-neutron correlation data for both the spontaneous fission of 252Cf and the thermal neutron induced fission of 235U. The codes presented here can be used to study the possible uses of neutron-neutron correlations in the area of transparency measurements and the uses of neutron-neutron correlations in coincidence neutron imaging.
Homma, Yuto; Moriwaki, Hiroyuki; Ohki, Shigeo; Ikeda, Kazumi
2014-06-01
This paper deals with verification of three dimensional triangular prismatic discrete ordinates transport calculation code ENSEMBLE-TRIZ by comparison with multi-group Monte Carlo calculation code GMVP in a large fast breeder reactor. The reactor is a 750 MWe electric power sodium cooled reactor. Nuclear characteristics are calculated at beginning of cycle of an initial core and at beginning and end of cycle of equilibrium core. According to the calculations, the differences between the two methodologies are smaller than 0.0002 Δk in the multi-plication factor, relatively about 1% in the control rod reactivity, and 1% in the sodium void reactivity.
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.
Computational Monte Carlo (MC) codes have been used for simulation of nuclear installations mainly for internal monitoring of workers, the well known as Whole Body Counters (WBC). The main goal of this project was the modeling and simulation of the counting efficiency (CE) of a WBC system using three different MC codes: MCNPX, EGSnrc and VMC in-vivo. The simulations were performed for three different groups of analysts. The results shown differences between the three codes, as well as in the results obtained by the same code and modeled by different analysts. Moreover, all the results were also compared to the experimental results obtained in laboratory for meaning of validation and final comparison. In conclusion, it was possible to detect the influence on the results when the system is modeled by different analysts using the same MC code and in which MC code the results were best suited, when comparing to the experimental data result. (author)
MCNP: a general Monte Carlo code for neutron and photon transport
MCNP is a very general Monte Carlo neutron photon transport code system with approximately 250 person years of Group X-6 code development invested. It is extremely portable, user-oriented, and a true production code as it is used about 60 Cray hours per month by about 150 Los Alamos users. It has as its data base the best cross-section evaluations available. MCNP contains state-of-the-art traditional and adaptive Monte Carlo techniques to be applied to the solution of an ever-increasing number of problems. Excellent user-oriented documentation is available for all facets of the MCNP code system. Many useful and important variants of MCNP exist for special applications. The Radiation Shielding Information Center (RSIC) in Oak Ridge, Tennessee is the contact point for worldwide MCNP code and documentation distribution. A much improved MCNP Version 3A will be available in the fall of 1985, along with new and improved documentation. Future directions in MCNP development will change the meaning of MCNP to Monte Carlo N Particle where N particle varieties will be transported
Applications of FLUKA Monte Carlo code for nuclear and accelerator physics
Battistoni, Giuseppe; Brugger, Markus; Campanella, Mauro; Carboni, Massimo; Empl, Anton; Fasso, Alberto; Gadioli, Ettore; Cerutti, Francesco; Ferrari, Alfredo; Ferrari, Anna; Lantz, Matthias; Mairani, Andrea; Margiotta, M; Morone, Christina; Muraro, Silvia; Parodi, Katerina; Patera, Vincenzo; Pelliccioni, Maurizio; Pinsky, Lawrence; Ranft, Johannes; Roesler, Stefan; Rollet, Sofia; Sala, Paola R; Santana, Mario; Sarchiapone, Lucia; Sioli, Maximiliano; Smirnov, George; Sommerer, Florian; Theis, Christian; Trovati, Stefania; Villari, R; Vincke, Heinz; Vincke, Helmut; Vlachoudis, Vasilis; Vollaire, Joachim; Zapp, Neil
2011-01-01
FLUKA is a general purpose Monte Carlo code capable of handling all radiation components from thermal energies (for neutrons) or 1keV (for all other particles) to cosmic ray energies and can be applied in many different fields. Presently the code is maintained on Linux. The validity of the physical models implemented in FLUKA has been benchmarked against a variety of experimental data over a wide energy range, from accelerator data to cosmic ray showers in the Earth atmosphere. FLUKA is widely used for studies related both to basic research and to applications in particle accelerators, radiation protection and dosimetry, including the specific issue of radiation damage in space missions, radiobiology (including radiotherapy) and cosmic ray calculations. After a short description of the main features that make FLUKA valuable for these topics, the present paper summarizes some of the recent applications of the FLUKA Monte Carlo code in the nuclear as well high energy physics. In particular it addresses such top...
Simulating fast transients with fuel behavior feedback using the Serpent 2 Monte Carlo code
Simulating transients with reactivity feedback effects using Monte Carlo neutron transport codes can be used for validating deterministic transient codes or estimating for example the total deposited energy in a fuel rod following a known reactivity insertion in the system. Recent increases in computational power as well as developments in calculation methodology makes it possible to obtain a coupled solution for several aspects of the multi-physics problem in a single calculation. This paper describes the different methods implemented in Serpent 2 Monte Carlo code that enable it to model fast transients with fuel behavior feedback. The capability is demonstrated in a prompt critical pin-cell case, where the transient is shut down by the negative reactivity from rising fuel temperature. (author)
Dose conversion coefficients for ICRP110 voxel phantom in the Geant4 Monte Carlo code
Martins, M. C.; Cordeiro, T. P. V.; Silva, A. X.; Souza-Santos, D.; Queiroz-Filho, P. P.; Hunt, J. G.
2014-02-01
The reference adult male voxel phantom recommended by International Commission on Radiological Protection no. 110 was implemented in the Geant4 Monte Carlo code. Geant4 was used to calculate Dose Conversion Coefficients (DCCs) expressed as dose deposited in organs per air kerma for photons, electrons and neutrons in the Annals of the ICRP. In this work the AP and PA irradiation geometries of the ICRP male phantom were simulated for the purpose of benchmarking the Geant4 code. Monoenergetic photons were simulated between 15 keV and 10 MeV and the results were compared with ICRP 110, the VMC Monte Carlo code and the literature data available, presenting a good agreement.
ERSN-OpenMC, a Java-based GUI for OpenMC Monte Carlo code
Jaafar EL Bakkali
2016-07-01
Full Text Available OpenMC is a new Monte Carlo transport particle simulation code focused on solving two types of neutronic problems mainly the k-eigenvalue criticality fission source problems and external fixed fission source problems. OpenMC does not have any Graphical User Interface and the creation of one is provided by our java-based application named ERSN-OpenMC. The main feature of this application is to provide to the users an easy-to-use and flexible graphical interface to build better and faster simulations, with less effort and great reliability. Additionally, this graphical tool was developed with several features, as the ability to automate the building process of OpenMC code and related libraries as well as the users are given the freedom to customize their installation of this Monte Carlo code. A full description of the ERSN-OpenMC application is presented in this paper.
A new Monte Carlo code for absorption simulation of laser-skin tissue interaction
Afshan Shirkavand; Saeed Sarkar; Marjaneh Hejazi; Leila Ataie-Fashtami; Mohammad Reza Alinaghizadeh
2007-01-01
In laser clinical applications, the process of photon absorption and thermal energy diffusion in the target tissue and its surrounding tissue during laser irradiation are crucial. Such information allows the selection of proper operating parameters such as laser power, and exposure time for optimal therapeutic. The Monte Carlo method is a useful tool for studying laser-tissue interaction and simulation of energy absorption in tissue during laser irradiation. We use the principles of this technique and write a new code with MATLAB 6.5, and then validate it against Monte Carlo multi layer (MCML) code. The new code is proved to be with good accuracy. It can be used to calculate the total power bsorbed in the region of interest. This can be combined for heat modelling with other computerized programs.
PyMercury: Interactive Python for the Mercury Monte Carlo Particle Transport Code
Iandola, F N; O' Brien, M J; Procassini, R J
2010-11-29
Monte Carlo particle transport applications are often written in low-level languages (C/C++) for optimal performance on clusters and supercomputers. However, this development approach often sacrifices straightforward usability and testing in the interest of fast application performance. To improve usability, some high-performance computing applications employ mixed-language programming with high-level and low-level languages. In this study, we consider the benefits of incorporating an interactive Python interface into a Monte Carlo application. With PyMercury, a new Python extension to the Mercury general-purpose Monte Carlo particle transport code, we improve application usability without diminishing performance. In two case studies, we illustrate how PyMercury improves usability and simplifies testing and validation in a Monte Carlo application. In short, PyMercury demonstrates the value of interactive Python for Monte Carlo particle transport applications. In the future, we expect interactive Python to play an increasingly significant role in Monte Carlo usage and testing.
PyMercury: Interactive Python for the Mercury Monte Carlo Particle Transport Code
Monte Carlo particle transport applications are often written in low-level languages (C/C++) for optimal performance on clusters and supercomputers. However, this development approach often sacrifices straightforward usability and testing in the interest of fast application performance. To improve usability, some high-performance computing applications employ mixed-language programming with high-level and low-level languages. In this study, we consider the benefits of incorporating an interactive Python interface into a Monte Carlo application. With PyMercury, a new Python extension to the Mercury general-purpose Monte Carlo particle transport code, we improve application usability without diminishing performance. In two case studies, we illustrate how PyMercury improves usability and simplifies testing and validation in a Monte Carlo application. In short, PyMercury demonstrates the value of interactive Python for Monte Carlo particle transport applications. In the future, we expect interactive Python to play an increasingly significant role in Monte Carlo usage and testing.
The use of an inbuilt importance generator for acceleration of the Monte Carlo code MCBEND
Monte Carlo is currently the most accurate method for the analysis of neutron and gamma-ray transport. However its application, especially to deep penetration studies, is costly in terms of the man-days to set up the calculation and in terms of computer usage. The MAGIC module, developed at the Winfrith Technology Centre, addresses both these problems. It employs an automated procedure based upon the established technique of splitting/roulette with an importance function derived from the solution of the adjoint diffusion equation. Examples are given of the application of the module with Monte Carlo code MCBEND
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)
The development of depletion program coupled with Monte Carlo computer code
The paper presents the development of depletion code for light water reactor coupled with MCNP5 code called the MCDL code (Monte Carlo Depletion for Light Water Reactor). The first order differential depletion system equations of 21 actinide isotopes and 50 fission product isotopes are solved by the Radau IIA Implicit Runge Kutta (IRK) method after receiving neutron flux, reaction rates in one group energy and multiplication factors for fuel pin, fuel assembly or whole reactor core from the calculation results of the MCNP5 code. The calculation for beryllium poisoning and cooling time is also integrated in the code. To verify and validate the MCDL code, high enriched uranium (HEU) and low enriched uranium (LEU) fuel assemblies VVR-M2 types and 89 fresh HEU fuel assemblies, 92 LEU fresh fuel assemblies cores of the Dalat Nuclear Research Reactor (DNRR) have been investigated and compared with the results calculated by the SRAC code and the MCNPREBUS linkage system code. The results show good agreement between calculated data of the MCDL code and reference codes. (author)
A new assembly-level Monte Carlo neutron transport code for reactor physics calculations
This paper presents a new assembly-level Monte Carlo neutron transport code, specifically intended for diffusion code group-constant generation and other reactor physics calculations. The code is being developed at the Technical Research Centre of Finland (VTT), under the working title 'Probabilistic Scattering Game', or PSG. The PSG code uses a method known as Woodcock tracking to simulate neutron histories. The advantages of the method include fast simulation in complex geometries and relatively simple handling of complicated geometrical objects. The main drawback is the inability to calculate reaction rates in optically thin volumes. This narrows the field of application to calculations involving parameters integrated over large volumes. The main features of the PSG code and the Woodcock tracking method are introduced. The code is applied in three example cases, involving infinite lattices of two-dimensional LWR fuel assemblies. Comparison calculations are carried out using MCNP4C and CASMO-4E. The results reveal that the code performs quite well in the calculation cases of this study, especially when compared to MCNP. The PSG code is still under extensive development and there are both flaws in the simulation of the interaction physics and programming errors in the source code. The results presented here, however, seem very encouraging, especially considering the early development stage of the code. (author)
Development of a space radiation Monte Carlo computer simulation based on the FLUKA and ROOT codes
Pinsky, L; Ferrari, A; Sala, P; Carminati, F; Brun, R
2001-01-01
This NASA funded project is proceeding to develop a Monte Carlo-based computer simulation of the radiation environment in space. With actual funding only initially in place at the end of May 2000, the study is still in the early stage of development. The general tasks have been identified and personnel have been selected. The code to be assembled will be based upon two major existing software packages. The radiation transport simulation will be accomplished by updating the FLUKA Monte Carlo program, and the user interface will employ the ROOT software being developed at CERN. The end-product will be a Monte Carlo-based code which will complement the existing analytic codes such as BRYNTRN/HZETRN presently used by NASA to evaluate the effects of radiation shielding in space. The planned code will possess the ability to evaluate the radiation environment for spacecraft and habitats in Earth orbit, in interplanetary space, on the lunar surface, or on a planetary surface such as Mars. Furthermore, it will be usef...
Srna - Monte Carlo codes for proton transport simulation in combined and voxelized geometries
Ilić Radovan D.
2002-01-01
Full Text Available This paper describes new Monte Carlo codes for proton transport simulations in complex geometrical forms and in materials of different composition. The SRNA codes were developed for three dimensional (3D dose distribution calculation in proton therapy and dosimetry. The model of these codes is based on the theory of proton multiple scattering and a simple model of compound nucleus decay. The developed package consists of two codes: SRNA-2KG and SRNA-VOX. The first code simulates proton transport in combined geometry that can be described by planes and second order surfaces. The second one uses the voxelized geometry of material zones and is specifically adopted for the application of patient computer tomography data. Transition probabilities for both codes are given by the SRNADAT program. In this paper, we will present the models and algorithms of our programs, as well as the results of the numerical experiments we have carried out applying them, along with the results of proton transport simulation obtained through the PETRA and GEANT programs. The simulation of the proton beam characterization by means of the Multi-Layer Faraday Cup and spatial distribution of positron emitters obtained by our program indicate the imminent application of Monte Carlo techniques in clinical practice.