Issue 
EPJ Nuclear Sci. Technol.
Volume 8, 2022
Fuel Cycle Simulation TWoFCS 2021



Article Number  4  
Number of page(s)  14  
DOI  https://doi.org/10.1051/epjn/2021030  
Published online  31 January 2022 
https://doi.org/10.1051/epjn/2021030
Regular Article
DONJON5/CLASS coupled simulations of MOX/UO_{2} heterogeneous PWR core
^{1}
École Polytechnique de Montréal, Institut de Génie Nucléaire,
P.O. Box 6079, Station “CentreVille”,
Montréal,
QC,
Canada
^{2}
IJCLab, IN2P3CNRS, Université ParisSaclay,
91406
Orsay Cedex,
France
^{3}
Laboratoire Subatech, IMTAIN2P3/CNRSUniversité,
Nantes
44307,
France
^{*} email: maxime2.paradis@polymtl.ca
Received:
15
July
2021
Received in final form:
9
December
2021
Accepted:
9
December
2021
Published online: 31 January 2022
Most fuel cycle simulation tools are based either on fixed recipes or assembly calculations for reactor modeling. Due to the high number of calculations and extensive computational power requirements, fullcore computations are often seen as not viable for this purpose. However, this leads to additional hypotheses and modeling biases, thus limiting the realism of the resulting fuel cycle. For several applications, the current modeling method is sufficient, but precise calculations of discharged fuel composition may require further refinements. CLASS (Core Library for Advanced Simulation Scenarios) is a dynamic fuel cycle simulation code developed since 2012 with reactor models based on neural networks to produce nuclear data and physical quantities. Past work has shown a first coupling between CLASS and DONJON5 to quantify neural networks approach biases. This work assesses the applicability of 3D fullcore diffusion calculations using the DONJON5 code coupled with nuclear scenario simulations involving a realistic PWR core at equilibrium cycle conditions. DONJON5 interpolates burnup dependent diffusion coefficients and cross sections generated beforehand by DRAGON5, a deterministic lattice calculation tool. Whereas previous studies considered only homogeneous reactors (i.e. homogeneous assembly in terms of composition and enrichment as well as homogeneous core), the present contribution focuses on the integration of fullcore calculations in CLASS for fuel cycles involving a MOX/UO_{2} PWR core (i.e. 1/3 MOx–2/3 UOx). The DONJON5 model considered in this work describes a core with critical boron concentration at each time step partially loaded with MOx heterogeneous assemblies composed of three enrichments. In fuel cycle calculations, the main issue is to adapt, in the fabrication stage, the fresh fuel composition for the reactor with regards to the isotopic composition of the available stocks. This work presents a fuel loading model based on power peaking factors minimization that respects irradiation cycle length, ^{235}U enrichment as well as Pu concentration and fissile quality, hence, ensuring a more uniform power distribution in the core.
© M. Paradis et al., Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
In fuel cycle simulators, reactor models are usually based on infinite assembly calculations both for the depletion simulation and for the fuel loading model that calculate each fresh fuel composition in fuel cycle simulations [1]. In the CLASS package, reactor models rely mostly on ANN (Artificial Neural Networks) to predict infinite neutron multiplication coefficients and average crosssections ((n, f), (n, γ) and (n, 2n)) to solve Bateman equations. Those neural networks are trained on a dedicated 1group crosssections database built with any fuel depletion software. Although these models are pretty convenient since they are easy to build and very effective in terms of computational costs, they do not consider core scale physics and, thus, lead to modeling biases.
In order to increase the reliability of fuel cycle simulations, reference [2] proposed a first coupling between DONJON5 [3] and CLASS [4] allowing fullcore calculations for PWRs (Pressurized Water Reactors) in a fuel cycle simulator. This work quantified biases induced by the use of infinite calculations in a fullscale scenario exploring MOx fuels deployment in PWRs after two decades of UOx irradiation. The scenario studied involved homogeneous UOx and MOx reactors (i.e. homogeneous assembly in terms of fresh composition as well as homogeneous core) whereas it is well known that PWRs are loaded with only 30% of MOx for safety reasons.
This work offers a coupling between CLASS and DONJON5 for heterogeneous reactors: PWR loaded with both MOx and UOx fuels. The impact of fullcore physics is then investigated on both the fresh fuel loading model and the irradiation model.
The purpose of loading models is to calculate simultaneously the enrichment of the UOx assemblies with the plutonium content (as a function of its isotopic quality) in the MOx assemblies in order to match two reactor criteria, namely the targeted burnup and the power peaking factor. For reactor operations, MOx assemblies are supposed to release as much energy as UOx assemblies. A minimal power peaking factor will then ensure an acceptable power distribution in the core but also an equivalence between MOx and UOx fuels in terms of energy production.
The irradiation model relies on fullcore calculations, with critical boron concentration at all time steps. All lattice calculations have been performed with DRAGON5 [5]. These are used to build the reactor database used by DONJON5 for fullcore simulations and to build artificial neural networks for models based on assembly considerations.
The contents of this paper is as follows. First, in Section 2, DRAGON5 lattice calculations are detailed. Section 3 describes how diffusion and ANN databases are built from lattice calculations while Section 4 describes calculation options for fullcore calculations with DONJON5. The fuel loading model (also referred to as the equivalence model) that estimates the Pu content as a function of the uranium enrichment in the UOx fuels and the plutonium isotopic quality are presented in Section 5 followed respectively, in Sections 6 and 7, by CLASS analyses of an elementary scenario involving a single MOx/UO_{2} fuelled PWR and complex scenario that involves both UOx and MOx/UO_{2} heterogeneously fuelled PWRs.
2 Assembly calculations
The computer code DRAGON5 is a lattice code developed at École Polytechnique de Montréal that enables the resolution of the neutron transport equation [5]. Specifically, this code allows the characterization of the internal reactor lattice, namely unit cells or fuel assemblies, via physics models articulated to achieve the following processes: interpolation of microscopic crosssections supplied by means of standard libraries; resonance selfshielding calculation; multigroup and multidimensional neutron flux calculation; transporttransport or transportdiffusion equi valence calculation; editing of condensed and homogenized nuclear properties for reactor calculations; and finally isotopic depletion calculations.
2.1 Fuel assembly
The heterogeneous PWR core studied here is composed of different assembly designs for UOx and MOx fuels (cf. Fig. 1). The MOx assembly is divided in three zones characterized by different Pu contents to attenuate power variations at the MOxUOx interface. Pu contents in the three zones were chosen to agree with the loading patterns introduced by the French operator EdF [6].
While guide tube cells are simulated, control rods are not considered in this fuel assembly modeling. Hence, core reactivity management is entirely ensured via boron diluted in the moderator.
To simplify the present study, all DRAGON5 infinite assembly calculations were performed for thermal power densities as well as fuel, clad and moderator temperatures fixed to a specific and theoretical mean value. A 4rings discretization was used in fuel rods in order to reduce calculation costs while keeping sufficient precision in flux simulation [7]. Finally, a brief study of the guide tube cell discretization pointed out that 3rings in the internal moderator suffice to properly take into account neutron moderation in water holes [8]. Table 1 summarize the dimension of geometrical elements previously discussed and Table 2 gives additional information on assembly material properties. Note that the dimension of geometrical elements are considered at cold state (20 °C), and thus, no thermal expansion was considered.
The tracking of the fuel assembly geometry is generated via the NXT: module (this tracking is being used to compute collision probabilitiesrequired for solving the transport equation) with white boundary conditions due to greater efficiency in terms of computational costs, when compared to mirrorlike boundary conditions. Note that B1 homogeneous models are used here [10].
It is considered preferable to use a simplified tracking for selfshielding calculations which can lead to substantial computing costs. The approach adopted here involves tracking individual cells in the assembly coupled by interface current coupling (SYBILT: module). The USS: module is then used to create physical probability tables required by the subgroup selfshielding approach.
We considered a SHEM multigroup library composed of 295 energy groups detailed in [11]. Simply, the SHEM295 defers from SHEM361 in its coarser meshing of the energetic domain from 4.0 to 22.5 eV [12]. This library accuracy was evaluated by measuring discrepancies generated on unloaded isotopic inventories resulting from deterministic calculations via DRAGON5 when compared to a continuous energy calculation performed via SERPENT(i.e. Monte Carlo code [13]). It was shown that the 295groups library induced negligible discrepancies (< 1% on ^{235}U and < 0.5% on ^{239}Pu inventories) for reasonable computing time scale [8].
It should be noted that, although MOx assembly calculations should be carried out in a UOx environment to properly characterize neutronic behaviors expected at UOx/MOx interfaces [14], we considered infinite lattice calculations for UOx and MOx fuels. Previous work showed that this simplification can lead to discrepancies reaching 6.5% on Pu^{239} inventories at EOC for a typical MOx assembly [15]. However, this work focuses on the quantification of bias generated by the current ANN approach used in CLASS and, thus, it was considered preferable to use ANN and diffusion databases built from identical lattice calculation schemes.
Fig. 1 UOx (a) and MOx (b) assembly geometries. Materials are identified as follows: moderator in green and fuel clad in purple. For the UOx assembly, fuel pellets are in yellow. For the MOx assembly, yellow, orange and red fuel pellets are characterized by high (9.77%_{HM}), medium (6.49%_{HM}) and low (3.66%_{HM}) Pu contents. 
Other assembly properties.
2.2 Reflector
Crosssection calculations for the core reflector require a more complex geometry due to the following modeling challenges:
 1.
the distance between peripheral fuel assemblies and the radial reflector is not constant around the core.
 2.
In the absence of fissile material in the reflector, a fissile environment must be simulated, namely surrounding fuel assemblies, in order to induce a neutron flux in the reflector itself.
 3.
The combined thickness of all elements composing the reflector amounts to a much larger geometry than the original assembly defined for fullcore calculations.
Taking into consideration the whole reflector would tend toward very costly lattice calculations. As for the varying distance between the baffle and core skin, it would lead to an important database dedicated to the reflector material in addition to the database summarizing the fuel crosssections ultimately entailing costly interpolations in fullcore calculations. Thus, the reflector characterization justifies the development of a judicious crosssections calculation scheme based on diverse approximations. The approach adopted in the previous coupling of CLASS and DONJON5 [16] and inspired by [17] goes as follows:
 1.
reflector material calculation is performed only at BOC (i.e. t = 0 s). Indeed, while the neutron flux varies in regards to the burnup, the neutron spectrum variation observed in the simulated fuel assembly has a limited impact on the reflector composition itself.
 2.
Reflector geometry is simplified in two ways: the reflector elements simulated are limited to elements included in the volume of a fuel assembly (cf. Tab. 3); the fuel assembly simulated to account for fissile environment is loaded with UOx (cf. Fig. 2).
Note that two fuel assemblies are simulated. Giving the symmetric conditions imposed at the upper and lower limit of the geometry shown in Figure 2, simulated assemblies can be reduced to their halves. It should also be noted that reflective conditions are used on the left of the fuel assemblies, while void conditions are defined beyond the reflective material.
Fig. 2 Reflector geometry. Materials are identified as follows: moderator (fuel assembly) in green, UOx fuel pellet in purple, iron (baffle) in yellow, iron (skin) in red and moderator (reflector) in orange. 
3 Database characteristics
Considering that the isotopic composition of the MOx introduced in PWRs is a priori unknown in fuel cycle calculations, the multiparameter databank built via lattice calculations and used for the interpolation of neutron characteristics requires a large sampling of significant parameters. The databank created in the present study, denoted MULTICOMPO in the following, includes 8 dimensions : boron concentration, ^{235}U enrichment of UOx fuels and isotopic composition (which amount to 5 dimensions) as well as Pu content of MOx fuels. It should be noted that, while the Pu content and the isotopic composition are varied here, the Pu content ratios between the three enrichment zones of MOx assembly designs are fixed.
3.1 Diffusion database
The sampling space must include all possible configurations that one may encounter in fuel cycle simulations. Consequently, minimal and maximal ^{235}U enrichments are chosen equal to 3.0 and 5.0%_{HM}, respectively.It was shown a posteriori that those limits enable fuel irradiation between 260 EFPD and 390 EFPD thereby embedding scenarios of interest studied via CLASS. An analog approach was performed to also ensure that the Pu content limits (i.e. 4.5 and 13.5%_{HM}) allowed a critical operation of the heterogeneous core with either very good or very bad isotopic compositions.
The sampling for ^{235}U enrichment, Pu contents and boron concentration are summarized in Table 4. Cubic interpolation method was shown to be satisfactory [16].
It is worth mentioning that crosssections saved in this hyperspace are burnup dependent. That is to say they are calculated at 25 nonregular depletion steps (first fine then coarse) as detailed in Table 5. This burnup sampling enables a precise transcription of the significant neutron spectrum evolution expected in the early stage of irradiation for thermal nuclear reactors (i.e. xenon and samarium effect) and coarse representation of the long time behavior of the irradiation. This ensures a reasonable database size while providing adequate interpolation of macroscopic crosssections in DONJON5 calculations [8]. As an example, Figure 3 presents k_{∞} evolution obtained for different fuel compositions.
Table 7 presents the plutonium and americium 241 sampling points considered for the diffusion database. The ^{239}Pu concentration acts here as a buffer for the plutonium isotopic vector. It is worth noting that, while only six isotopes are perturbed, the crosssections for 22 isotopes selected either for their preponderant impact or because of their relatively high concentration are saved in the diffusion database. Reliable isotopic inventories may be obtained with only a few isotopic perturbations. Adding sampling points (for instance additional ^{238}Pu or ^{241}Am perturbations) leads to far greater computational costs (i.e. 150% increase when compared to samplings presented in Tab. 7) for each crosssection interpolation in module NCR: in fullcore calculation. A contrario, reducing the ^{242}Pu sampling leads to discrepancies of approximately 1.2% on this isotope and 0.2% on ^{235}U after only one irradiation cycle of 280 days (equivalent full power).
Finally, note that crosssections resulting from the lattice calculation are 2groups condensed and homogenized spatially on the whole assembly before they are stored on the diffusion database.
Diffusion database fuel perturbations.
Diffusion database time sampling.
Fig. 3 Infinite multiplication factor evolution resulting from lattice calculations performed for UOx (green line) and MOx (other lines)fuel assemblies. Isotopic compositions for UOX1 as well as MOX1, MOX2 and MOX3 are presented in Table 6. 
Isotopic compositions studied for UOx and MOx assemblies characterization.
Diffusion database isotopic perturbations.
3.2 Neural network database
For the sake of clean comparisons, an independent database was built via DRAGON5 calculations in order to train ANNs employed in usual reactor models for CLASS. This enables us to eliminate sources of errors in fuel cycle calculations that are not induced by neural networks or physic models (i.e. different depletion calculation software, nuclear data or different geometrical descriptions).
For ANNs, it is well known that regular sampling of input parameters may lead to significant training biases. Hence, fresh fuel compositions were considered as parameters yielding 600 initial compositions for UOx fuels and 2000 for MOX fuels sampled with the LHS technique (La tin Hyper Square). Furthermore, it has been suggested to set boron concentration to zero at all computed burnup steps to yield correct evaluation of plutonium evolution [16]. All the other parameters remained identical.
4 Fullcore calculations
4.1 Geometry description
In fuel cycle studies, dozens of core calculations are performed for each reactor. Considering the resulting computational requirements, it appears necessary to simplify the fullcore geometry illustrated in Figure 4. The elementary cell for the core is ahomogenized assembly corresponding to a cubic element with height equal to the fuel assembly side. Note that no additional spatial discretization is considered.
As one can see, fuel assemblies are radially surrounded by cubic elements corresponding to the reflector introduced in Section 2. Furthermore, the core geometry proposed here includes additional upper and lower planes representing the reflective mixture (i.e. axial reflectors). To simplify fullcore geometry tracking, these two planes have the same height as the core cells. It is also worth mentioning that the core height is split in 9 nonregular planes. This coarse discretization led to reasonable calculation costs and was shown to be adequate given the constant moderator and fuel temperatures within the core [8]. The core dimension and geometrical elements are summarized in Table 8.
Fig. 4 Real (a) and simulated (b) PWR quarter core layout. Assemblies are displayed in white, whereas baffle (light gray), barrel (dark gray), vessel (gray) and thermal insulator (orange) are replaced by reflector elements (gray). 
Core geometrical data.
4.2 Heterogeneous modeling for CLASS
We considered a fixed cyclic fuel reloading scheme consisting of 40 new fuel assemblies (i.e. 28 identical UOx and 12 identical MOx). Notethat the number of UOx assembly decreases from 28 to 25 after the third irradiation cycle to maintain a total of 157 assemblies. The proportion of MOx assemblies in the core and the four cycles reloading scheme are defined in accordance with the present industrial practice (i.e. PARITÉ MOX introduced by EdF [18]) as well as the previous coupling considering UOx and MOx homogeneous cores [16].
The heterogeneous reactor models developed are divided in two successive steps involving 3D fullcore calculations.
The first calculation step involves a core fully fuelled with UOx in order to determine the minimal ^{235}U enrichment required to maintain criticality for a given output burnup objective. In detail, two UOx core calculations characterized respectively by enrichment equal to 3.0 and 5.0 %_{HM} are first performed. Then, the enrichment is adjusted in subsequent UOx fullcore calculations (i.e. dichotomic search) until the targeted burnup is reached. The loading pattern used in those core calculations is identical to the heterogeneous pattern displayed on Figure 5. MOx assemblies are then simply substituted by their UOx counterpart.
A second calculation step is then performed to characterize the irradiation of a heterogeneous core composed of UOx assemblies defined in the previous homogeneous DONJON5 calculations as well as energyequivalent MOx assemblies. This equivalence between UOx and MOx is discussed in Section 5.
Both UOx and heterogeneous fullcore calculations include 9 computed burnup steps as detailed in Table 9.
Fig. 5 Loading pattern used in 3D fullcore modeling of heterogeneous PWR core. The 26 fuel assemblies indicated correspond to combustion zones composing one eighth of the core. 
Time discretization for fullcore calculations.
4.3 Acceleration of calculations
As fullcore calculations lead to computational costs far greater than ANNs typically used in CLASS, several ways to improve calculation speed have been studied.
4.3.1 Critical ^{235}U enrichment predictions
Considering that UOx fullcore calculations are essentially performed to obtain critical fissile contents, we propose a direct relation between ^{235}U enrichment and cycle length, allowing us to dismiss entirely UOx calculations in CLASS. 14 UOx fullcore calculations were performed apriori with burnup targets covering evenly irradiation cycles between 260 and 390 days (equivalent full power). Note that these cycles were defined in such a way as to incorporate a multitude of fuel management schemes put forward in the literature [18–20]. Critical uranium enrichments were then extracted and used to generate a cubic polynomial regression fit as displayedin Figure 6.
One can see that it is possible to correctly predict (for this loading pattern) the critical ^{235}U enrichment evolution in relation of the cycle length with this relation: (1)
where ϵ_{crit.} and l_{cycle} correspond to the critical ^{235}U enrichment (%_{HM}) and the cycle length (equivalent full power days), respectively. Moreover, constants A to D are determined empirically and are presented in Table 10. Figure 6 shows negligible discrepancies between predicted and calculated enrichments resulting from DONJON5. The biases induced by this accelerated approach are in the range 0.023%_{HM} ; 0.014%_{HM}.
In the interest of quantifying errors generated on heterogeneous fullcore calculations by critical enrichment predictions, we compare unloaded inventories resulting from two heterogeneous PWR cores : one fuelled with the predicted enrichment characterized by the largest deviation noted in Figure 6 and another one fuelled with the corresponding calculated enrichment. We obtained relative errors smaller than 1.5% on all Pu and U isotopes with significant concentration. It should be noted that the largest mass deviation, that is to say 1.5%, is observed on ^{235}U inventory and is larger than the initial bias on the predicted enrichment. Considering that it yields calculations cost 12.5% less important, this accelerated approach is adopted albeit the minor bias on initial fissile inventory.
Fig. 6 Critical ^{235}U enrichment evaluation for UOx fullcore calculations performed with DRAGON5DONJON5. 
Details of the critical enrichment prediction fit.
4.3.2 Reloading burnup predictions
Due to the a priori unknown used assemblies burnup at each reloading step, the usual practice relies on an initial calculation with fuel assemblies reloaded with burnups in direct proportion to the mean irradiation rhythm. In order to calculate real incore burnups for all assemblies composing the heterogeneous core at equilibrium, preliminary irradiation cycles are performed until convergence is reached in all combustion zones. To avoid these preliminary cycles altogether (from 4 to 8cycles in our case) and greatly accelerate 3D fullcore calculations, we propose to predict reloading burnup values using ANNs. Based on the core symmetry visible on Figure 5, we limit burnup predictions to 26 combustion zones (i.e. 20 used assemblies and 6 new assemblies) located in one eighth of the core.
To generate these neural networks, a multidimensional training databank covering campaign length included between 260 and 390 irradiation days (equivalent full power) was constructed. This hyperspace also take into consideration the ^{235}U enrichment in UOx assemblies and the Pu content as well as the isotopic composition in MOx assemblies. Again, an LHS method was used to randomly generate core characteristics within the hyperspace of interest. A total of 3000 fullcore calculations were performed and resulting unloaded burnups were extracted for all 20 used assemblies and saved in 20 distinct training databases.
For accuracy validation, 300 additional and independent fullcore calculations were performed in order to validate the 20 trained neural networks. Overall, 95.52% of burnups predicted via the 20 ANNs tested in this experiment are characterized by a discrepancy under 2%, when compared to the correspondent calculated burnup value.
Furthermore, relative discrepancies between calculated and predicted burnups obtained for a UOX1MOX1 heterogeneous core (cf. Tab. 6) are indicated in Figure 7. In this figure, fresh fuels are not represented as they have a zero initial burnup.
It is worth recalling that, although calculated burnups are determined for all 157 assemblies, predicted values are obtained only for one eighth of the core and redistributed by symmetry. Power asymmetries generated in preliminary cycles, which in itself leads to varying burnups within symmetric combustion zones, are not reproduced by the accelerated approach proposed here. Consequently, slight deviations will be observed in symmetric combustion zones (p. ex. combustion zones A8, H1, O8, and H15). However, these variations are less then 0.2% and can be decreased significantly by reducing the convergence criteria imposed to the diffusion equation solutions obtained with the FLUD: module in fullcore calculations [8].
Fig. 7 Burnup discrepancies between predicted and calculated values at reloading. 
4.3.3 Critical boron calculations
3D fullcore calculations are carried out with critical boron concentration at every computed burnup steps. Critical boron is determined in an iterative process similar to critical enrichment calculations as introduced in Section 4.3.1. Based on the stationary irradiation (i.e. the reloading scheme is cyclic and uniform), the core’s reactivity evolution can be assumed identical for every cycle. Thus, the critical boron concentration is calculated only for the first cycle and is reproduced in the subsequent ones. Previous works have shown that this approach yields negligible boron concentration discrepancies while decreasing substantially computational costs [8,16].
4.3.4 Resulting accelerated calculations
In order to estimate errors generated by the three simplifications (^{235}U enrichment fit, burnup neural networks and boron simplified calculation), we examined unloaded isotopic inventories resulting from one irradiation cycle of the UOX1MOX1 heterogeneous core. Figure 8 compares isotopic inventories obtained with typical and accelerated fullcore calculations. As one can see, deviations between isotopic inventories resulting from these approaches are very small. This is also the case for ^{238}U which is not introduced in this figure in order to simplify presentation. On the other hand, these two approaches are characterized by costs amounting to 4014 seconds (i.e. 480 seconds for the preliminary UOx fullcore calculation and 3534 seconds for the heterogeneous fullcore calculation) and 1185 seconds, respectively, leading to almost four times less computational costs for each core calculation.
Fig. 8 Unloaded actinides inventories obtained with full and accelerated calculations. 
5 Fuel fabrication: matching MOx composition to UO_{2} enrichment
Previous sections were dedicated to the elaboration of an optimized fullcore calculation scheme that may be integrated in the CLASS package for fuel cycle simulations. Although the resulting calculation scheme enables us to reliably simulate fuel irradiation in the PWR considered here, further development is required to insure adequate fuel loading in regards to the isotopic composition of the available stocks.
Whereas previous works (either based on infinite calculations [1] or on fullcore calculation [16]) depended entirely on core criticality as a criterion to calculate optimal Pu content, we propose an alternative approach based essentially on power peaking factor minimization for fresh MOx fabrication.
5.1 Power equivalence for fresh fuel determination
The approach consists of adjusting the plutonium content in fresh MOx fuels according to its isotopic composition and critical ^{235}U enrichment to obtain similar unloaded burnups for UOx and MOx assemblies. Considering that similar power measured on UOx and MOx assemblies lead to small burnup differences in unloaded assemblies for uniform UOx and MOx irradiation lengths (i.e. four cycles reloading scheme), we adopted the peaking factor at EOC (End Of Cycle) as a measure of equivalence between these two fuels [8]. Literature usually defines the peaking factor as the maximum ratio between a fuel pin power over pin mean power. Given that the elementary cell defined for 3D fullcore calculations is the assembly, the peaking factor definition is measured on fuel assemblies. As an example, Figure 9 recapitulates peaking factors resulting from DONJON5 calculations carried out with a fixed standard Pu quality (i.e. MOX1 in Tab. 6) and for a wide array of ^{235}U enrichment in UOx fuels and Pu content in MOx fuels. One can clearly see that for a given isotopic composition and ^{235}U enrichment, there is one Pu content that leads to a minimal peaking factor and thus, according to the relation previously mentioned, a small burnup differences in unloaded assemblies.
Fig. 9 Power peaking factors (PF) calculated for different Pu content and ^{235}U enrichment combinations with a fixed standard isotopic composition. 
5.2 ANN calculations for peaking factor estimation
The layout of the peaking factor surface displayed in Figure 9, which can be reproduced for any isotopic composition, gave rise to the methodology adopted for fuel loading: predicting peaking factors via an ANN and identifying which Pu content leads to the minimal factor. Consequently, another ANN is built to calculate peaking factors with respect to ^{235}U enrichment in UOx fuel and Pu content as well as isotopic composition in MOx fuel. Pseudorandom core characteristics were defined via an LHS method and lead to a training database built on 3000 fullcore calculations. To validate the neural network precision, predicted peaking factors were compared to calculated factors resulting from thousands of independent heterogeneous fullcore calculations. From this comparison, we found that 94.26% of predicted factors are characterized by relative deviations under 5%. On the other hand, it should be mentioned that 1.28% of predictions display absolute discrepancies over 0.2. Albeit their small proportion, large negative discrepancies significantly affect equivalent Pu content determination and thus limit the precision of this equivalence. Let us recall that the approach adopted here consists of predicting peaking factors at a fixed enrichment and for a fine and regular Pu content vector, and thus enables us to identify which content leads to the minimal factor.
5.3 Fresh fuel loading model for heterogeneous reactors
Finally, our reactor modeling with DONJON5 for dynamic fuel cycle calculations relies on three steps:
 1.
a calculation step involving a homogeneous core fuelled entirely with UOx in order to determine the critical ^{235}U enrichment.
 2.
An estimation of the Pu content according to isotopic compositions available in stocks and the critical enrichment resulting from the previous step. This content leads to a minimal peaking factor at EOC and ensures that the discharged burnups of MOx and UOx assemblies are similar.
 3.
A calculation step characterizing heterogeneous core depletion is finally performed for UOx and MOx assemblies as defined in previous steps and with critical boron concentrations follow up.
It is worth noting that the loading pattern of this reactor modeling is supposed to be constant. Hence, artificial neural networks and uranium enrichment calculations developed here are specific to this loading pattern. Further developments could include a loading pattern selection with respect to UOx and MOx fuel composition. Reference [8] shows that this methodology could be used with other typical loading patterns if required.
6 Elementary scenarios
Previous sections explained how fuel loading and fuel depletion modeling with DONJON5 were designed and implemented in a dynamic fuel cycle simulation within the CLASS package. Let us recall that this new calculation scheme differs from reactors modeling based on ANN [1] and commonly used approaches by taking into account neutrons phenomena at core level. Before investigating model biases in a full scenario calculation over several years, it appears necessary to highlight differences between the DONJON5 and ANN approach on basic park configurations involving a single fresh fuel loading followed by one irradiation cycle.
6.1 Reactor modeling with artificial neural networks
Usual CLASS reactor models rely on artificial neural networks to estimate the core criticality evolution and average crosssections of main reactions [1]. Note that these models perform heterogeneous PWRs simulation via independent UOx and MOx modeling. Mass and power of UOx and MOx fuel assemblies are then proportionally adjusted to represent reactor mass.
Specifically, each fresh composition is elaborated according to the estimation of the maximal achievable burnup for each assembly considered independently (i.e. based on their respective infinite calculations). Maximal achievable burnup is defined here as the burnup for which the average value of k_{∞} over the different batches does not compensate for neutron losses induced by structural material capture andneutron leakage. This criticality margin is represented by a single parameter denoted k_{threshold} in the following. As MOx fuels and UOx fuels have different neutron spectrum, different neutron leakage rates and neutron capture rates on structural materials are expected. While this margin depends on the fuel type in the present study, it does not depend on the fuel composition yet. Finally, the fuel loading model calculate the fresh fuel composition regarding its isotopic vector in such a way as to fulfill equation (2) during the whole fuel evolution. In other words, EOC burnups are achieved when the equality is no longer verified. (2)
UOx and MOx depletion are also calculated independently via ANNs that predict 1group crosssections of three reactions: (n, f), (n, γ) and (n, 2n) for significant isotopes. Bateman equations are then solved with a RungeKutta fourthorder method, ultimately yielding an estimation of discharged compositions.
Artificial neural networks considered in this work are MLP (Multi Layer Perceptrons) trained on a dedicated database built on neutron characteristics resulting from DRAGON5 calculations with the exact same simulation options used for the diffusion database (cf. Sect. 2). This modeling is denoted CLASSMLP in the following, whereas the notation CLASSDONJON5 refers to heterogeneous fullcore calculations.
6.2 Scenario analysis
Scenario calculations studied here (defined with the specific concerns to be as simple as possible) consider sufficient fissile materials to properly load one heterogeneous reactor. Three initial isotopic compositions have been considered for the comparison of DONJON5 and MLP approaches (cf. Tab. 6). The fuel fabrication process is set to two years. This is done to simulate chemical separation to extract plutonium and uranium inventories required for UOx and MOx fuel assemblies fabrication. Moreover, unloaded fuel assemblies are cooled for two years in order to account for radioactive decays expected in cooling pools. Further modeling options selected for this comparison are summarized in Table 11.
Figures 10(a) and 10(b) display uranium and plutonium evolution resulting from the simulated reactor, respectively. From these figures, we denote that plutonium inventories required for the reactor loading is higher for the DONJON5 approach, when compared to inventories predicted with MLP. Considering that initial heavy metals mass is fixed for all reactors in CLASS, plutonium underestimation necessarily lead to an overestimated proportion of depleted uranium in MOx fuels fabricated with MLP models. Therefore, less fissile materials are introduced in MLP cores and, given that the power is constant for all simulations, predicted ^{235}U fission rates should necessarily be higher in MLP models, when compared to DONJON5 calculations. Ultimately, this observation translates to a greater use of uranium via the CLASSMLP approach.
Regarding the plutonium inventories, note that discrepancies between CLASSDONJON5 and CLASSMLP are similar at BOC (Beginning Of Cycle) and EOC. Although some deviations may be seen for HET2, the general observation is that plutonium production from heterogeneous irradiation is reliably simulated with both approaches.
Thus, we conclude from this comparison that the main source of bias lies in the loading model. Let us recall that MLP and DONJON5 core reloading are based on distinct modeling philosophy : whereas the MLP model relies entirely on a criticality search, the DONJON5 model is based on peaking factor minimization. DONJON5 model leads to a higher plutonium inventory for fresh fuel fabrication that impacts uranium evolution and, eventually, discharged inventories. Overall these observations do not necessarily refute MLP models applicability for heterogeneous reactors, but highlight and challenge their limitations.
Reactor modeling description and denomination.
Fig. 10 Uranium (a) and plutonium (b) inventories simulated with CLASSDONJON5 and CLASSMLP approach for three different compositions denoted HET1, HET2 and HET3 (i.e. based on MOX1, MOX2 and MOX3, respectively). 
7 Complex scenarios involving UOx and heterogeneous reactors
Finally, CLASSDONJON5 coupling was evaluated on full scenarios involving MOx fuels loading in a reactor fleet composed of PWRs initially loaded with UOx assemblies. The purpose of this section is once again to quantify discrepancies between discharged isotopic inventories resulting from reactor models considering the core scale (i.e. CLASSDONJON5) and simplified models (i.e. CLASSMLP).
7.1 Scenario characteristics
As will be apparent shortly, studied scenarios were kept rather simple, yet realistic, due to their academic nature. Consequently, the reactor fleet considered here includes only 10 identical 900 MWe PWRs. Note that these PWRs are entirely loaded with UOx fuels and operated for 20 cycles. This is done in order to accumulate sufficient irradiated UOx stockpiles to constantly fabricate MOx fuels. Two UOx reactors are then partially loaded with MOx fuels. Those ten reactors (8 UOx and 2 UOx/MOx) are operated for an additional 30 cycles, amounting to approximately 50 years (standard reactor life time), before shutting down. A simplified illustration of this scenario is presented in Figure 11. Through the variation of stock management options, three different scenarios have been studied : a first scenario representing stationary fuel management as well as two scenarios describing dynamic operations characterized by distinct UOx and MOx/UOx discharged burnups. For all scenarios, UOx and MOx fuel fabrications are started at year 2.0 and year 24.5, respectively. Scenario descriptions are further detailed in Table 12.
Fig. 11 Representative schematization of the reactor fleet simulated in scenarios. 
Scenario denomination and reactors descriptions. Cycle lengths (l_{cycle}) are in full power equivalent days and discharged burnups (BU) are in GWd⋅t_{HM}.
7.2 Results and physical analysis of scenario A
Figure 12 presents the global evolution of plutonium inventories in UOx spent fuels and in MOx spent fuels for scenario A. Essentially, conclusions reached here conform to findings in Section 6. Pu content resulting from DONJON5 modeling is higher, when compared to the MLP fuel loading model, allowing for greater procurement from UOx stockpiles for MOx loading or, stated another way, less UOx spent fuel accumulation with DONJON5. It is worth clarifying that this overestimated Pu content calculated with DONJON5 does not reflect higher plutonium consumption at the fleet scale. Indeed, plutonium inventories in MOx stockpiles are also higher in CLASSDONJON5 simulations, when compared to CLASSMLP simulations.
As one can see in Figures 10 and 12, discharged plutonium inventories obtained with MLP models differed greatly from DONJON5 at the reactor scale. It is worth mentioning that these biases may not be as visible at fullcycle scale (i.e. in wastes and all cycle installations). Considering that plutonium discrepancies in MOx spent fuels predicted with MLP are primarily generated by the fuel loading model and assuming that fuel irradiation with MLP is roughly accurate, plutonium leftover is still accounted for in UOx spent fuel.
Figure 13 shows total plutonium and minor actinides inventories incycle. From this figure, we can see a remarkable compensation of biases for plutonium as simulations with MLP and DONJON5 give very similar plutonium evolution. Nonetheless, the minor actinides show strong biases, consistent with previous comparison between DONJON5 and MLP [16]. As 8 reactors are entirely fuelled with UOx fuels, differences are mainly driven by UOx behavior.
Fig. 12 Plutonium quantity accumulated UOx and heterogeneous stocks. 
Fig. 13 Global plutonium and minor actinides inventories accumulated incycle. 
7.3 Impact of reactor parameters
Scenario B differs from scenario A by its higher UOx discharged burnup and its smaller MOx discharged burnup. This leads to smaller Pu content for MOx fuels and lower isotopic plutonium quality, when compared to results from scenario A. Scenario C presents a higher MOx burnup. Also note that this scenario considers a different stock management option denoted Mix: MOx fuel fabrication is supplied in an alternating sequence between the last and first UOx spent fuel (i.e. LiFo and FiFo) resulting in a succession of degraded and improved isotopic compositions.
For the sake of brevity, only discrepancies on Pu and MA (Minor Actinides) inventories at the end of the simulation (i.e. 80 years) between CLASSMLP and CLASSDONJON5 are displayed for scenarios B and C (cf. Tab. 13). Differences between each scenario are not discussed here, but the total incycle plutonium inventory in scenario B (~100 tons) is higher than in scenarios A and C (~95 tons). Note in addition that differences between scenarios are much higher than discrepancies between MLP and DONJON5 approaches for a given scenario.
We then conclude that, although reactor models based on MLP do not represent accurately reactor physics and thus lead to possible erroneous evaluation of discharged inventories, they seem suitable for global plutonium evolution in large fuel cycle calculations.Stated another way: while calculations performed with MLP can yield acceptable values for global inventories, MLP models do not predict correctly material availability for fresh fuel fabrications as they induce a significant error on plutonium location in the fuel cycle.
Relative differences (%) for plutonium and minor actinides between simulations performed with DONJON5 and MLP models after 80 years.
8 Conclusion
In this paper, DONJON5 fullcore diffusion calculations were effectively coupled to CLASS for dynamic fuel cycle simulations. Some procedures were designed to accelerate fullcore computations while keeping the scenario code CLASS capability to dynamically adapt reactor fuel loading to respect available fissile materials specifically taking into account isotopic quality. Great efforts were invested to ensure critical depletion at all time and to adequately simulate fuel environment whilst maintaining reasonable calculation costs. Original methods were also developed to avoid the initial balancing processes and to reduce the number of iterations to determine the boron critical content in the moderator.
The coupling of DONJON5 with CLASS for heterogeneous reactors required the development of a dedicated original fuel loading model. The one proposed here relies on identifying Pu content to ensure a minimal power peaking factor. It differs from previous loading models based on infinite neutron multiplication coefficient considerations.
Fullcore calculations were used to analyze some of the parameters considered in the neural networks based approach and their impact on scenario accuracy. It was shown that the main source of discrepancy for plutonium inventory estimations is primarily due to the fuel loading model. The irradiation model was also shown unsatisfactory for ^{235}U incore inventory prediction as well as minor actinides production. The above observations are coherent with previous work like [21] and with physics [22] because ^{235}U and minor actinides are much more sensitive to neutron flux induced effects. The 2groups diffusion scheme for fullcore calculations seems to be a minimum requirement to take these differences into account.
One possible way to improve neural networks models would be to use the fullcore database built for this work to calculate fresh fuel compositions based on power peaking factors for heterogeneous reactors. This database also enables a k_{threshold} calculations that are adapted for Pu quality. In that case, it was shown that MLP processed results are both good (less than a 1% discrepancy on ore consumption and plutonium production) and very fast (few minutes compared to hours for coupled executions). This raises the question of code flexibility: although k_{threshold} models specifically developed for each reactor geometry could lead to more precise fuel fabrication with ANNs in CLASS, this would lead to considerably less flexible physic models. On the other hand, it seems that no intrinsic core effects can be simulated without considering specific core characteristics.
Funding
This work was partly funded by the Natural Sciences and Engineering Research of Council (NSERC) of Canada.
Author contribution statement
All the authors have read and approved the final manuscript.
References
 B. Leniau et al., A neural network approach for burnup calculation and its application to the dynamic fuel cycle code CLASS, Ann. Nucl. Energy 81, 125–133 (2015) [CrossRef] [Google Scholar]
 M. Guillet et al., Coupled CLASS and DONJON5 3D fullcore calculations and comparison with the neural network approach for fuel cycles involving MOX fueled PWRs, Ann. Nucl. Energy 152, 107971 (2021) [CrossRef] [Google Scholar]
 A. Hébert, D. Sekki, R. Chambon, A User Guide for DONJON Version 5, Institut de génie nucléaire, Département de génie mécanique, École Polytechnique de Montréal. Montréal QC, Canada, Tech. Rep. IGE344, 2019 [Google Scholar]
 B. Mouginot et al., Core library for advanced scenario simulation, CLASS: principle & application, in International Conference “The Role of Reactor Physics toward a Sustainable Future” (PHYSOR 2014) (2014), pp. 12 [Google Scholar]
 G. Marleau, A. Hébert, R. Roy, A User Guide for DRAGON Version 5, Institut de génie nucléaire, Département de génie physique, École Polytechnique de Montréal. Montréal QC, Canada, Tech. Rep. IGE335, 2018 [Google Scholar]
 R. Largenton, Modélisation du comportement effectif du combustible MOX par une analyse micromécanique en champsde transformation non uniformes. PhD thesis, Université de Provence AixMarseille, 2012 [Google Scholar]
 A. Santamarina, C. Collignon, C. Garat, French calculation schemes for light water reactor analysis, in The Physics of Fuel Cycles and Advanced Nuclear Systems: Global Developments (PHYSOR 2004) (2004) [Google Scholar]
 M. Paradis, Prise en compte de la chaîne de calcul DRAGON5DONJON5 dans la simulation de réacteurs à eau sous pression hybride via CLASS, Master’s thesis, École Polytechnique de Montréal, 2021 [Google Scholar]
 R. Vallerent, Développement et validation de schémas de calcul à double niveau pour les réacteurs à eau sous pression, Master’s thesis, École Polytechnique de Montréal, 2009 [Google Scholar]
 I. Petrovic, P. Benoist, G. Marleau, A quasiisotropic reflecting boundary condition for the TIBERE heterogeneous leakage model, Nucl. Sci. Eng. 122, 151–166 (1996) [CrossRef] [Google Scholar]
 A. Hébert, Development of the subgroup projection method for resonance selfshielding calculations, Nucl. Sci. Eng. 162, 56–75 (2009) [CrossRef] [Google Scholar]
 A. Santamarina, N. Hfaiedh, The SHEM energy mesh for accurate fuel depletion and BUC calculations, International Conference on Nuclear CriticalitySafety (ICNC 2007), 2007 [Google Scholar]
 J. Leppänen, Serpent – A continuousenergy Monte Carlo reactor physics burnup calculation code, VTT Technical Research Centre of Finland (2013), vol. 4 [Google Scholar]
 D. Bernard, A. Santamarina, A.M. Malvagi, L. Daudin, Experimental Validation of the APOLLO2 code for High Burnup MOx Fuel. JEF2.2 Results and JEFF3.0 Improvements, in The Physics of Fuel Cycles and Advanced Nuclear Systems: Global Developments (PHYSOR 2004) (2004) [Google Scholar]
 A. Somaini, Analyse des erreurs induites par une modélisation simplifiée sur l’évolution des combustibles REP: Impact des fuites neutroniques dans les calculs cellules. PhD thesis, Université ParisSaclay, 2017 [Google Scholar]
 M. Guillet, Prise en compte des calculs de coeur via l’utilisation de DONJON5 pour le code de scénario nucléaire CLASS et évaluation de l’impact sur les inventaires isotopiques. Master’s thesis, École Polytechnique de Montréal, 2019 [Google Scholar]
 N. Bejaoui, Simulation numérique de l’effet du réflecteur radial sur les cellules REP en utilisant les codes DRAGON et DONJON. Master’s thesis, École Polytechnique de Montréal, 2012 [Google Scholar]
 “Décision n^{o} 2017DC0608 de l’Autorité de sûreté nucléaire du 5octobre 2017 relative à la gestion de combustible ≪Parité MOX≫ de certains réacteurs électronucléaires exploités par électricité de France (EDF),” tech. rep., Autorité de sûreté nucléaire, 2017 [Google Scholar]
 A. Worrall, Effect of plutonium vector on core wide nuclear design parameters, tech. rep., IAEA, 2000 [Google Scholar]
 S. Aniel, J. Bergeron, A. Puill, Evaluation of the maximum content of a MOXfueled pressurized water reator versus isotopiccomposition with respect to the void coefficient, in IAEA TECDOC 941 (1997), pp. 181–193 [Google Scholar]
 A. Somaini et al., The impact of reactor model simplification for fuel evolution: a bias quantification for fuel cycle dynamic simulations, in 2016 International Congress on Advances in Nuclear Power Plants (ICAPP 2016) (2016), pp. 1045–1053 [Google Scholar]
 P. Reuss, Précis de neutronique (EDP Sciences, 2012) [Google Scholar]
Cite this article as: Maxime Paradis, Xavier Doligez, Guy Marleau, Marc Ernoult, Nicolas Thiollière, DONJON5/CLASS coupled simulations of MOX/UO_{2} heterogeneous PWR core, EPJ Nuclear Sci. Technol. 8, 4 (2022)
All Tables
Scenario denomination and reactors descriptions. Cycle lengths (l_{cycle}) are in full power equivalent days and discharged burnups (BU) are in GWd⋅t_{HM}.
Relative differences (%) for plutonium and minor actinides between simulations performed with DONJON5 and MLP models after 80 years.
All Figures
Fig. 1 UOx (a) and MOx (b) assembly geometries. Materials are identified as follows: moderator in green and fuel clad in purple. For the UOx assembly, fuel pellets are in yellow. For the MOx assembly, yellow, orange and red fuel pellets are characterized by high (9.77%_{HM}), medium (6.49%_{HM}) and low (3.66%_{HM}) Pu contents. 

In the text 
Fig. 2 Reflector geometry. Materials are identified as follows: moderator (fuel assembly) in green, UOx fuel pellet in purple, iron (baffle) in yellow, iron (skin) in red and moderator (reflector) in orange. 

In the text 
Fig. 3 Infinite multiplication factor evolution resulting from lattice calculations performed for UOx (green line) and MOx (other lines)fuel assemblies. Isotopic compositions for UOX1 as well as MOX1, MOX2 and MOX3 are presented in Table 6. 

In the text 
Fig. 4 Real (a) and simulated (b) PWR quarter core layout. Assemblies are displayed in white, whereas baffle (light gray), barrel (dark gray), vessel (gray) and thermal insulator (orange) are replaced by reflector elements (gray). 

In the text 
Fig. 5 Loading pattern used in 3D fullcore modeling of heterogeneous PWR core. The 26 fuel assemblies indicated correspond to combustion zones composing one eighth of the core. 

In the text 
Fig. 6 Critical ^{235}U enrichment evaluation for UOx fullcore calculations performed with DRAGON5DONJON5. 

In the text 
Fig. 7 Burnup discrepancies between predicted and calculated values at reloading. 

In the text 
Fig. 8 Unloaded actinides inventories obtained with full and accelerated calculations. 

In the text 
Fig. 9 Power peaking factors (PF) calculated for different Pu content and ^{235}U enrichment combinations with a fixed standard isotopic composition. 

In the text 
Fig. 10 Uranium (a) and plutonium (b) inventories simulated with CLASSDONJON5 and CLASSMLP approach for three different compositions denoted HET1, HET2 and HET3 (i.e. based on MOX1, MOX2 and MOX3, respectively). 

In the text 
Fig. 11 Representative schematization of the reactor fleet simulated in scenarios. 

In the text 
Fig. 12 Plutonium quantity accumulated UOx and heterogeneous stocks. 

In the text 
Fig. 13 Global plutonium and minor actinides inventories accumulated incycle. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.