Analysis of VIP-BWR reactor core physics experiments on UO2 and MOX mockup fuel assemblies - part 2: analysis with SIMULATE5

Analysis of the three test cores in the VIP-BWR program was performed in a two-dimensional geometrical model with CASMO5 coupled with the JENDL-4.0-based neutron data library, and reported in the previous paper. Following the study, interpretation of the experiments were carried out in a three-dimensional geometrical model with SIMULATE5 for the code validation study. The nuclear libraries for the SIMULATE5 calculations were generated with CASMO5 with the JENDL-4.0-based neutron data library. The effective multiplication factors of the critical cores ranged from 0.9983 to 1.0023 with measurement uncertainties of 0.0003 to 0.0004 (one σ). The root mean squares of (the calculated/the measured-1) for the fission rates at the core-mid plain of all the measured fuel rods were about 3% for the three cores. It was noticed that the calculations underestimated the fission rates of the UO₂ fuel rods and overestimated those of the MOX fuel rods for the test cores loaded with MOX fuel rods, which was consistent with trends in the preceding analysis studies of the VIP-BWR program and other MOX core experiments, and the biases were confirmed in the calculation results of power distributions in MOX-fueled light water reactor cores.     

Development of Spectral History and Pin Power Correction Methods for Pin-by-Pin Core Analysis Method Using Three-Dimensional Direct Response Matrix

Spectral history and pin power correction methods have been developed for the pin-by-pin core analysis method using the three-dimensional direct response matrix (3D-DRM). The direct response matrix is formalized using four subresponse matrices in order to respond to a core eigenvalue k and thus it can be recomposed at each outer iteration in the core analysis. For core analysis, it is necessary to take into account the historical effect, which is related to spectral heterogeneity. The spectral history method is used to evaluate the nodal burn-up spectrum obtained by using the outgoing neutron current instead of the nodal flux because the 3D-DRM method does not use the nodal flux. The pin power correction method corrects the fuel rod neutron production rates obtained in the pin-by-pin calculation. These two methods were tested in a heterogeneous system. The test results show that the neutron multiplication factor error and nodal neutron production rate errors can be reduced by half during burn-up. The root-mean-square differences between the relative fuel rod neutron production rate distributions and the maximum error of relative fuel rod production rate can also be reduced by half. This means that the developed methods can reflect the effects of intra- and interassembly heterogeneities during burn-up and can be used for core analysis.     

Analysis of core physics experiments on irradiated BWR MOX fuel in the REBUS program

Critical experiments were performed in the REBUS program on a core loaded with a test bundle including 16 irradiated BWR-type MOX rods of average burnup of 61 GWd/t. The experimental data were analyzed using diffusion, transport, and continuous-energy Monte Carlo calculation codes coupled with nuclear data libraries based on JENDL-3.2 or JENDL-3.3. Biases in effective multiplication factors of the critical cores were ?1.0%Δk for the diffusion calculations (JENDL-3.2), ?0.3%Δk for the transport calculations (JENDL-3.3), and 0.2%Δk for the Monte Carlo calculations (JENDL-3.2). The measured core fission rate and co-activation rate distributions were generally well reproduced using the three types of calculations. The burnup reactivity determined using the measured water level reactivity coefficients was ?2.41 ± 0.08%Δk/kk’, which also agreed with the results of the three type of calculations within the measurement and calculation errors. The most probable isotopic inventories in the irradiated MOX rods was tentatively obtained by using the ratios of the calculation to chemical assay data on a pellet sample, and the burnup reactivity was reanalyzed to split the calculation error into those due to the inventory and reactivity calculations. This approach showed that the inventory calculation error compensated the reactivity calculation error.     

Hafnium clad fuels for fast spectrum BWRs

We show that by use of hafnium cladding, a fast neutron spectrum is achievable in the top of uprated BWRs. Monte Carlo calculations have been made for Hf clad inert matrix nitride and low fertile MOX fuels, with fuel segments located in the upper part of an uprated BWR, where the coolant void fraction exceeds 70%. The nitride fuel results in the hardest neutron spectrum, but the low fertile MOX fuel still yields fission probabilities for even neutron number nuclides similar to those of sodium cooled reactors. The inert matrix nitride fuel configuration yields high burning rates, permitting to stabilise TRU inventories with less than 50% BWR cores of the here suggested type in the power park. The core with low fertile MOX fuel is less efficient, but still a zero net producer of TRU. Fuel and coolant temperature feedbacks are affected by introduction of absorbing elements in the fuel, but remain within acceptable ranges for the low fertile MOX fuel. Although control rod worths are reduced, shutdown margins are sufficient to ensure sub-criticality in cold conditions. From a materials point of view, the behaviour of hafnium clad MOX fuel would be similar to zircalloy clad MOX fuel already used extensively in nuclear industry. Thus, if dynamic stability of the core can be ensured, the here proposed fuel may be considered as a low cost solution for transmutation of minor actinides on industrial scale.     

Plutonium recycling in hydride fueled PWR cores

The classic approach to the recycling of Pu in PWR is to use mixed U-oxide Pu-oxide (MOX) fuel. The mono-recycling of plutonium in PWR transmutes less than 30% of the loaded plutonium, providing only a limited reduction in the long-term radiotoxicity and in the inventory of TRU to be stored in the repository. The primary objective of this study is to assess the feasibility of plutonium recycling in PWR in the form of plutonium hydride, PuH₂, mixed with uranium and zirconium hydride, ZrH_1.6, referred to as PUZH, that is loaded uniformly in each fuel rod. The assessment is performed by comparing the performance of the PUZH fueled core to that of the MOX fueled core. Performance characteristics examined are transmutation effectiveness, proliferation resistance of the discharged fuel and fuel cycle economics. The PUZH loaded core is found superior to the MOX fueled core in terms of the transmutation effectiveness and proliferation resistance. For the reference cycle duration and reference fuel rod diameter and pitch, the percentage of the plutonium loaded that is transmuted in one recycle is 53% for PUZH versus 29% for MOX fuel. That is, the net amount of plutonium transmuted in the first recycle is 55% higher in cores using PUZH than in cores using MOX fuel. Relative to the discharged MOX, the discharged PUZH fuel has smaller fissile plutonium fraction - 45% versus 60%, 15% smaller minor actinides (MA) inventory and more than double spontaneous fission neutron source intensity and decay heat per gram of discharged TRU. Relative to the MOX fuel assembly, the radioactivity of the PUZH fuel assembly is 26% smaller and the decay heat and the neutron yield are only 3% larger. The net effect is that the handling of the discharged PUZH fuel assembly will be comparable in difficulty to that of the discharged MOX assembly while the proliferation resistance of the TRU of the discharged PUZH fuel is enhanced.     

A macroscopic cross-section model for BWR pin-by-pin core analysis

A macroscopic cross-section model used in boiling water reactor (BWR) pin-by-pin core analysis is studied. In the pin-by-pin core calculation method, pin-cell averaged cross sections are calculated for many combinations of core state and depletion history variables and are tabulated prior to core calculations. Variations of cross sections in a core simulator are caused by two different phenomena (i.e. instantaneous and history effects). We treat them through the core state variables and the exposure-averaged core state variables, respectively. Furthermore, the cross-term effect among the core state and the depletion history variables is considered. In order to confirm the calculation accuracy and discuss the treatment of the cross-term effect, the k-infinity and the pin-by-pin fission rate distributions in a single fuel assembly geometry are compared. Some cross-term effects could be negligible since the impacts of them are sufficiently small. However, the cross-term effects among the control rod history (or the void history) and other variables have large impacts; thus, the consideration of them is crucial. The present macroscopic cross-section model, which considers such dominant cross-term effects, well reproduces the reference results and can be a candidate in practical applications for BWR pin-by-pin core analysis on the normal operations.     

Monte Carlo Modelling of Increasing Void Fraction in 100% MOX ABWR: Lessons Drawn from the FUBILA Program

The French Atomic Energy Commission CEA and the Japanese Incorporated Administrative Agency JNES (Japan Nuclear Energy Safety Organisation) have undertaken first-of-a-kind full MOX core physics experiments, FUBILA, in the EOLE critical facility of the CEA Cadarache Centre. The experiments have been designed to obtain core physics data under high-burn-up 9 × 9 and 10 × 10 BWR MOX assemblies operating conditions. The experimental program, consisting of eight different core configurations, started in January 2005 and ended on September 1, 2006. The analysis of the void increase part of the experimental data between 0 and 70% void has been carried out using the French TRIPOLI-4.5 continuous-energy Monte Carlo calculation code with the newly released JEFF3.1.1 nuclear data library. The average C/E discrepancies obtained on critical masses, reactivity worth, and pin-by-pin power distributions enable us to estimate all the integral and local parameters with uncertainties largely within the target uncertainties, demonstrating the capability of the code to treat complex geometries with a high degree of accuracy. Additional k_eff calculations performed with the latest ENDF/B-VII evaluation exhibit a clear tendency to overestimate the k_eff by about 500 to 650 pcm and the void worth by more than 4%, showing that the JEFF3.1.1 library is more precise for MOX lattices.     

A methodology to predict a fission gas release ratio of MOX fuel with heterogeneous microstructure

The objective of this study is to formulate a methodology to predict a fission gas release ratio of MIMAS MOX. An irradiated MIMAS MOX fuel with plutonium rich agglomerates was subjected to elemental analyses by electron probe micro analysis and secondary ion mass spectrometry in order to investigate xenon distribution. The results of the elemental analyses showed that the plutonium rich agglomerates at the periphery of the fuel pellet sample retained a high concentration of xenon as gas bubbles. Then, the results were used as reference data for modification of models in a fuel rod analysis code, FEMAXI-7. Using the modified FEMAXI-7, we applied an approach to prediction of fission gas release ratio of MOX fuel with plutonium rich agglomerates. In the approach, two separated analyses using FEMAXI-7 were performed for the plutonium rich agglomerates and the matrix. Fission gas release ratios obtained from the two analyses were processed through weighted-average with burnup ratios of the plutonium rich agglomerates and the matrix. Finally, the fission gas release ratios were compared with results of rod puncture tests. As a result of the comparison, it was confirmed that the proposed approach could well predict fission gas release ratio of MOX fuel with plutonium rich agglomerates.     

Analysis of MOX fuel behavior in reduced-moderation water reactor by fuel performance code FEMAXI-RM

To assess the feasibility of the 31% Pu-MOX fuel rod design of reduced-moderation water reactor (RMWR) in terms of thermal and mechanical behaviors, a single rod assumed to be irradiated in the core of RMWR up to 106 GWd/tHM has been analyzed by a fuel performance code FEMAXI-RM which is an extended version of FEMAXI-6 code. In the analysis, design specifications of fuel rod and irradiation conditions have been input, and available models of both MOX fuel and UO₂ fuel have been used as appropriate. The results are: fission gas release is several tens of percent, rod internal pressure does not exceed the coolant pressure, and the highest fuel center temperature is 2400 K, while cladding diameter increase caused by pellet swelling is within 1% strain. These predictions suggest that the MOX fuel rod integrity will be held during irradiation in RMWR, though actual behavior of MOX pellet swelling and cladding oxidation require to be investigated in detail.     

基于蒙特卡罗方法的三维燃耗计算研究

采用通过编写连接MCNP程序和ORIGEN2程序的接口处理程序的方法进行快中子系统的燃耗计算。由MCNP、ORIGEN2、接口处理程序和截面文件组成的软件系统可用于燃料或堆芯非均匀布置快中子系统的燃料同位素成分和燃耗反应性损失计算,在燃耗反应性损失计算中采用了伪裂变产物的方法。介绍程序系统的研制情况,并给出用该软件系统计算中国实验快堆首炉堆芯和OECD/NEAMOX燃料快堆基准题的燃耗计算结果。     

MOX Co-deposition Tests at RIAR for SF Reprocessing Optimization

The Oxide Electrowinning method has been studied as one of the candidate dry reprocessing concepts of the future fuel cycle systems. On the MOX co-deposition process, the main process of that method, some fundamental experiments have been performed to confirm its feasibility. In the experiments, several parameters were set to study the suitable electrolysis condition to obtain desired granule of MOX. The concentrations of uranium, plutonium, fission products(FP) simulators, and corrosion products(CP) simulators were adopted as the parameters. The blowing gas composition (O₂, Cl₂, Ar) during the electrolysis was also set as the variable condition. Through these experiments, it was clarified that the partial pressure of chlorine gas during electrolysis was important to obtain MOX granule with high Pu concentration (about 30%) without generating bottom precipitation in melt. Finally, adequacy of the process control method for MOX co-electrolysis was confirmed through the test using spent fast reactor(FR) fuel.     

Application of Quick Subchannel Analysis Method for Three-Dimensional Pin-by-Pin BWR Core Calculations

Three-dimensional pin-by-pin core analysis is considered to be a candidate for the next-generation BWR core calculation method. In our previous study, the applicability of the transport and burnup calculations for a three-dimensional pin-by-pin BWR core analysis was investigated. However, the thermal-hydraulics calculation has not yet been studied in this framework. In the conventional core analysis code, the bundlewise thermal-hydraulics calculation is adopted. In the actual core analysis, the power distribution inside a fuel assembly is tilted at the region adjacent to a control blade or the core peripheral region. In these regions, the consideration of the subchannel-wise void distribution has an impact on the fission rate distribution. Therefore, an evaluation of the detailed void distribution inside an assembly, i.e., the incorporation of the subchannel wise void distribution, is desirable for the pin-by-pin BWR core analysis. Although several subchannel analysis codes have been developed, these subchannel analysis codes generally require a large computational effort to estimate the subchannel-wise void distribution in a whole BWR core. Therefore, to analyze a whole BWR core within a reasonable computation time, it was necessary to apply a fast subchannel analysis code. In this paper, a quick subchannel analysis code dedicated to pin-by-pin BWR core analysis is newly developed, and the void distribution of the present subchannel analysis code is compared with the prevailing subchannel analysis code NASCA using three-dimensional single-assembly geometries. Since the present subchannel analysis code is used for a coupled neutronics/thermal-hydraulics analysis, the results of the coupling calculation are also compared with those of NASCA. The calculation result indicates that the void distribution difference between NASCA and the present subchannel analysis code is slightly less than 10%. This result indicates that the prediction accuracy of the present subchannel analysis code will be reasonably appropriate for a pin-by-pin BWR core analysis. Furthermore, the results show that the calculation time of the present subchannel analysis code is only 10 min for a hypothetical three-dimensional ABWR quarter-core geometry using a single CPU. This calculation time is sufficient for a pin-by-pin BWR core analysis.     

基于非均匀变分节块法的pin-by-pin计算加速算法研究

计算效率是制约pin-by-pin计算工程应用的主要因素之一。本文利用三维扩散的非均匀变分节块法的非均匀节块的描述能力,在不改变原问题栅元均匀化材料分布的前提下,将传统pin-by-pin计算中使用的均匀材料细网剖分方式替代为非均匀材料粗网剖分方式（粗网加速方法）,既能保证pin-by-pin的计算分辨率,又能显著降低红 黑迭代所需的浮点数操作数目,减小内迭代的计算代价。针对外迭代,运用广义矩阵分离加速（GPM）算法和粗网有限差分（CMFD）算法提高源迭代的收敛速度,降低计算时间。数值结果表明,提出的加速算法能在保证计算精度的前提下,有效提高pin-by-pin计算的效率。     

Pu recycling in a full Th-MOX PWR core. Part I: Steady state analysis

Current practice of Pu recycling in existing Light Water Reactors (LWRs) in the form of U-Pu mixed oxide fuel (MOX) is not efficient due to continuous Pu production from U-238. The use of Th-Pu mixed oxide (TOX) fuel will considerably improve Pu consumption rates because virtually no new Pu is generated from thorium. In this study, the feasibility of Pu recycling in a typical pressurized water reactor (PWR) fully loaded with TOX fuel is investigated.Detailed 3-dimensional 100% TOX and 100% MOX PWR core designs are developed. The full MOX core is considered for comparison purposes. The design stages included determination of Pu loading required to achieve 18-month fuel cycle assuming three-batch fuel management scheme, selection of poison materials, development of the core loading pattern, optimization of burnable poison loadings, evaluation of critical boron concentration requirements, estimation of reactivity coefficients, core kinetic parameters, and shutdown margin.The performance of the MOX and TOX cores under steady-state condition and during selected reactivity initiated accidents (RIAs) is compared with that of the actual uranium oxide (UOX) PWR core.Part I of this paper describes the full TOX and MOX PWR core designs and reports the results of steady state analysis. The TOX core requires a slightly higher initial Pu loading than the MOX core to achieve the target fuel cycle length. However, the TOX core exhibits superior Pu incineration capabilities.The significantly degraded worth of control materials in Pu cores is partially addressed by the use of enriched soluble boron and B₄C as a control rod absorbing material. Wet annular burnable absorber (WABA) rods are used to flatten radial power distribution. The temperature reactivity coefficients of the TOX core were found to be always negative. The TOX core has a slightly reduced, as compared to UOX core, but still sufficient shutdown margin.In the TOX core β_eff is smaller by about a factor of two in comparison to the UOX core and even lower than that of the MOX core. The combination of small β_eff and reduced control materials worth may potentially deteriorate the performance under RIA conditions and requires an additional examination. The behavior of the considered cores during the most limiting RIAs, such as rod ejection, main steam line break, and boron dilution, is further investigated and reported in Part II of the paper.     

Applicability of the Diffusion and Simplified P3 Theories for Pin-by-Pin Geometry of BWR

The pin-by-pin fine-mesh core calculation method is considered as a candidate next-generation core calculation method for BWR. In this study, the diffusion and simplified P3 (SP3) theories are applied to the BWR pin-by-pin fine-mesh calculation. The performances of the diffusion and SP3 theories for cell-homogeneous pin-by-pin fine-mesh calculation for BWR are evaluated through comparison with a cell-heterogeneous detailed transport calculation by the method of characteristics (MOC). Two-dimensional, 2 × 2 multi-assemblies geometry is used to compare the prediction accuracies of the diffusion and SP3 theories. The 2 × 2 multi-assemblies geometry consists of 9 × 9 UO₂ fuel assemblies that have two different enrichment splittings. To minimize the cell-homogenization error, the SPH method is applied for the pin-by-pin fine-mesh calculation. The SPH method is a technique that reproduces a result of heterogeneous calculation using that of homogeneous calculation. The calculation results indicated that the diffusion theory shows a discrepancy larger than that of the SP3 theory on the pin-wise fission rate distribution. In contrast to the diffusion theory, the SP3 theory shows a much better accuracy on the pin-wise fission rate distribution. The computation time using the SP3 theory is about 1.5 times longer than that using the diffusion theory. The BWR core analysis consists of various calculations, e.g., the cross section interpolation, neutron flux calculation, thermal hydraulic calculation, and burn-up calculation. The function of the calculation time for the neutron flux calculation is usually less than half in the typical BWR core analysis. Therefore, the difference in the calculation time between the diffusion and SP3 theories would have no significant impact on the calculation time of the BWR core analysis. For these reasons, the SP3 theory is more suitable than the diffusion theory and is expected to have sufficient accuracy for the 2 × 2 multi-assemblies geometry used in this study, which simulates a typical situation of the actual BWR core.     

Study of the PWR REA pulse width for realistic UO2 and MOX core designs using 3D kinetics methods

For the development of new acceptance criteria for the analysis of Rod-Ejection-Accidents (REA) in Pressurized Water Reactor UO2 cores, full-width-at-half-maximum power pulse widths in the range 25–40 ms were employed in the analytical transient fuel behaviour studies. Due to the operation of MOX cores in Switzerland, a study using the CORETRAN 3D kinetics code was carried out at PSI to investigate and compare on the basis of real operated cycles MOX and UO2 cores, the pulse width magnitude as function of rod reactivity for postulated super-prompt critical REAs. This study is presented in a first part of this paper. Thereafter, to quantify the variation in pulse width when applying a different 3D kinetic solver, a benchmarking with the SIMULATE-3K code is carried out. This code-to-code comparison aims also to serve as an assessment of SIMULATE-3K for REA analyses of the Swiss cores since only CORETRAN had so far been benchmarked at PSI for this type of applications. Finally, to provide confidence in the results, particularly for pulse widths below 25 ms, sensitivity studies are performed to assess the effects of specific modelling options and assumptions applied in the REA 3D kinetic analysis and to compare the sensitivities between UO2 and MOX cores. On the basis of all these studies, it is found that the MOX core pulse width is usually around 5 to 10 ms lower than for UO2 cores. For rod worths close to 2.0$, corresponding to around 500 pcm above prompt criticality and which was obtained for many of the investigated MOX cycles particularly at end-of-cycle conditions, a representative MOX pulse width magnitude is found to be around 15 ms. To shortly address the applicability of these results, i.e. obtained for very specific MOX cores, additional calculations were performed for hypothetical cases with increased MOX core fractions and with higher MOX assembly Pu enrichments. And it was found that even for these types of cores, the pulse width would, although being slightly reduced, remain in the range 10–15 ms for rod reactivities up to 500 pcm above prompt criticality.     

Behavior of Irradiated ATR/MOX Fuel under Reactivity Initiated Accident Conditions

Pulse irradiation experiments with irradiated ATR/MOX fuel rods of 20MWd/kgHM were conducted at the NSRR in Japan Atomic Energy Research Institute to study the transient behavior of MOX fuel rod under reactivity initiated accident conditions. Four pulse irradiation experiments were performed with peak fuel enthalpy ranging from 335 J/g to 586 J/g, resulted in no failure of fuel rods. Relatively large radial deformation of the fuel rods due to pellet-cladding mechanical interaction occurred in the experiments with peak fuel enthalpy above 500 J/g. Significant fission gas release up to 20% was measured by rod puncture measurement. The generation of fine radial cracks in pellet periphery, micro-cracks and boundary separation over the entire region of pellet were observed. These microstructure changes might contribute to the swelling of fuel pellets during the pulse irradiation. This could cause the large radial deformation of fuel rod and high fission gas release when the pulse irradiation conducted at relatively high peak fuel enthalpy. In addition, fine grain structures around the plutonium spot and cauliflower structure in cavity of the plutonium spot were observed in the outer region of the fuel pellet.     

Analysis of Core Physics Experiments on Fresh and Irradiated PWR UO2 Fuels in the REBUS Program

Critical experiments of two cores each loaded with fresh 5 × 5 test PWR-type fuel rods of ²³⁵U enrichment of 3.8 wt% or irradiated 5 × 5 test rods of rod average burnup of 55 GWd/t in the REBUS program were analyzed using diffusion, transport, and continuous-energy Monte Carlo calculation codes coupled with nuclear data libraries based on JENDL-3.2 and JENDL-3.3. Biases in effective multiplication factors k _eff's of the critical cores were about ?1:2%Δk for the diffusion calculations (JENDL-3.2), ?0:5%Δk for the transport calculations (JENDL-3.3), and ?0:5 and 0.1%Δk for the Monte Carlo calculations (JENDL-3.3 and JENDL-3.2, respectively). The measured core fission rate and Sc- or Co-activation rate distributions were generally well reproduced using the three types of calculation. The burnup reactivity determined using the measured water level reactivity coefficients was ?2:35 ± 0:07Δk/kk′. The calculated result of the Monte Carlo calculations agreed with it; however, the diffusion and transport calculations overestimated the absolute value by about 7%, which would be mainly attributed to the errors in the calculation of the reactivity caused by changing the fuel compositions from fresh fuel to irradiated fuel.     

Analysis of Mixed Oxide Fuel Critical Experiments with Neutronics ANalysis Codes for Boiling Water Reactors

Critical experiments of UO₂ and full mixed oxide (MOX) fuel cores conducted at the Tank-type Critical Assembly (TCA) were aNalyzed using BWR design-purpose codes HINES and CERES with ENDF/B files and Monte Carlo fine analysis codes VMONT and MVP with the JENDL-3.2 library.

The averaged values of the multiplication factors calculated with HINES/CERES, VMONT and MVP agreed with those of experiments within 0.3%ΔAk. The values by the design-purpose codes showed a small difference of 0.1%Δk between UO₂ and MOX cores. Monte Carlo code results showed that the JENDL-3.2 library had a tendency to overestimate the multiplication factors of UO₂ cores by about 0.3%Δk compared with those values of MOX cores. The root mean square errors of calculated power distributions were less than 1% for HINES/CERES and VMONT.

These results showed that (1) the accuracy of these codes when applied to full MOX cores was almost the same as their accuracy for UO₂ cores, which confirmed the accuracy of present core design codes for full MOX cores; and (2) the accuracy of the 190-energy-group Monte Carlo calculation code VMONT was almost the same as that of the continuous-energy Monte Carlo calculation code MVP.