Issue
EPJ Nuclear Sci. Technol.
Volume 8, 2022
Fuel Cycle Simulation TWoFCS 2021
Article Number 1
Number of page(s) 8
DOI https://doi.org/10.1051/epjn/2021025
Published online 11 January 2022

© H. Barale et al., Published by EDP Sciences, 2022

Licence Creative CommonsThis 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

Scenario studies can be used to explore future options of the nuclear energy development, investigating the advantages and disadvantages, technological constraints, costs and overall attractiveness of different fuel cycles and reactor fleet developments. To further improve nuclear fuel use and nuclear waste management, advanced nuclear fuel cycles involving innovative reprocessing plants and reactors are being evaluated. As highlighted by the French PPE (Programmation Pluriannuelle de l’Énergie) [1], managing the spent fuel arising from nuclear power plants until its reuse or disposal is an important step of the nuclear fuel cycle and constitutes the so-called back-end. It is essential to address spent fuel management properly because this issue directly involves the sustainable development of nuclear energy, the nuclear non-proliferation, the nuclear energy costs, the natural resource utilization etc.

In order to assess the performances of these nuclear fuel cycles, the equilibrium simulation tool SEPAR (Equilibrium Simulator of Advanced Reactor Parks) has been recently developed at the Department of Reactor Studies (DER) at CEA Cadarache. It relies on a neural network approach [2] to estimate the equilibrium material composition and fluxes for each isotope to and from reactors, enrichment and reprocessing facilities. It models advanced nuclear fuel cycles in equilibrium conditions, i.e. in conditions which stabilize selected nuclear inventories such as spent nuclear fuel constituents, plutonium or some minor actinides. One of the major innovations brought by SEPAR is the calculation of the equilibrium of the material flows in a reactor fleet, this aspect is absent from the usual dynamic scenario simulation codes (such as COSI [3]), which on the other hand allows to simulate the transition from one fleet to another, year after year.

France is a major player in the recycling of nuclear materials in the world since it recycles plutonium from its spent Low Enriched Uranium (LEU) fuel in the form of MOX fuel. Recycling of the fissile uranium remaining in spent LEU fuel has also been undertaken in the past in the form of ERU fuel, and is expected to resume over the next decade [4]. However, the recycling of spent fuel remains incomplete today. To advance the recovery of materials and the reduction of nuclear waste, it is necessary to study advanced reactor fleets that use innovative recycling plants and/or new reactor concepts. For example, in the medium term, the renewal of the French fleet by a new generation of EPRs offers the opportunity to deploy new fuel management systems, capable of recycling all the spent fuel in the fleet [5]. SEPAR has the advantage of being extremely fast in calculation time, this feature allows the study and the future complete validation of the code to different reactors and related fuel cycle strategies at steady state. The fuel types took in consideration during this study are the Uranium OXide (UOX), the Mixed OXide (MOX), Enriched Reprocessed Uranium (ERU) and the MOX Enriched Uranium Support also known as MIX. This paper focuses on the application of SEPAR to different fuel options four different strategies: the UOX - MOX mono-recycling of Pu (the actual French strategy [6]), the UOX - ERU mono-recycling of U, the UOX - MIX full multi-recycling of Pu and the UOX - MOX - AMR partial multi-recycling. In Section 2 the first description of the SEPAR tool is presented. In Section 3 an analysis of the SEPAR results concerning the four aforementioned cases is proposed, as well as a primary comparison/validation with COSI6.

2 Methodologies

The modeling of a complex system such as a nuclear energy system requires physical modeling of all major processes occurring in the system. A nuclear energy scenario is generally defined by the components of the system, i.e. enrichment, fuel fabrication, reprocessingfacilities and reactors. The scenario aims to forecast the evolution or equilibrium state of the system depending on input parameters. The evolution is mostly captured in the material fluxes between the different components of the system.

Different types of scenario simulation tools have been developed, this study focuses on dynamic and equilibrium codes. Dynamic fuel cycle simulation codes model evolving nuclear fuel cycles, and calculate nuclides inventories and material flows in each unit of the cycle over a given period of time [7]. In equilibrium simulation the study is made at steady state, defined by the fact that all the material fluxes, at every step of the recycling processes are constant and that the isotopic composition of the fuel is the same at each cycle [8].

2.1 Equilibrium searching tool: SEPAR

SEPAR is a calculation module that simulates the operation of a nuclear reactor fleet and its associated fuel cycle facilities in a steady-stateregime. It is developed by SPRC (Reactor Physics and fuel Cycle Service) at the IRESNE Institute of the CEA/Cadarache. SEPAR is used to size nuclear fleets that operate a total or partial reprocessing of irradiated fuels. It characterizes recycling fleets with respect to the services they may provide in terms of isotopic streams at cycle level. This may concern the management of recyclable materials or waste, such as the reduction of a plutonium inventory or quantities of minor actinides to store in-depth, etc.

SEPAR development follows a first study that showed the interest in this method on EPR - SFR symbiotic fleet [9]. Currently, SEPAR-based studies aim to quickly provide detailed characteristics of advanced fleets under various operating conditions. In particular, it can highlight and explain the effects of prospective fuel batches that may be operated in any concept of existing or prospective nuclear reactors.

The SEPAR module is integrated into a code environment comprising:

  • upstream level: any irradiation code able to simulate the fuel composition after irradiation. These results feed SEPAR, which then builds neural network meta-models for irradiation as well as equivalence model for fresh fuel composition.

  • Equilibrium level: SEPAR is the module for material balance calculations, the system of non-linear equations is solved applying the GNU Scientific Library (GSL) using the Broyden solver [10]. Additional calculation modules allowing access to fleet performance beyond material balances were added at this stage.

  • Downstream level: once steady-state fleets are characterized, they can be introduced into dynamic fleet simulations, for instance using the CEA scenario code COSI6. At this stage, economic or environmental indicators in SEPAR are planned to be investigated.

Figure 1 illustrates how SEPAR schematizes an elementary nuclear system. X represents the fuel batch fraction, K is the isotopic composition treatment function at the reprocessing stage, t defines the transfer times between different facilities, while ϵ characterizes the transferred material fuel fraction between them.

How SEPAR works can be explained with an elementary example concerning only 2 isotopes. This example is described by the system of differential equations in equation (1), in which B is the branching factor for the decay of the parent isotope p to the daughter isotope i and λ is the isotope decay constant. {dmpdt=λpmpdmidt=λimi+Bpiλpmp. \begin{equation*} \begin{cases} \frac{dm_p}{dt} = -\lambda_p m_p \\ \frac{dm_i}{dt} = -\lambda_i m_i + B_{pi} \lambda_p m_p. \end{cases}\end{equation*}(1)

The solution to this equilibrium problem is well known, equation (2) represents the evolution of the daughter nucleus: mi(t)=Bpiλpλiλpmp0(eλpteλit)+mi0eλit\begin{equation*} m_i(t) = B_{pi} \frac{\lambda_p}{\lambda_i - \lambda_p} m^0_p \bigl(e^{\lambda_p t} - e^{\lambda_i t} \bigr) + m^0_i e^{\lambda_i t}\end{equation*}(2)

here, mi0$m^0_i$ and mp0$m^0_p$ denote the initial masses of daughter and parent nuclei respectively. An analogous expression can be derived for the parent nucleus.

Since the evolution of spent fuel from the reactor to another use in the reactor is not a simple decay, but it consists of different procedures as reprocessing, the expression becomes more complicated. The process can be divided into two steps: from the reactor to the reprocessing-step (Eq. (3)) and from the reprocessing-step back to a reactor (Eq. (4)). Given the cooling and fuel fabrication times (respectively tR2U2$t_{R_2 \to U_2}$,tU2R3$t_{U_2 \to R_3}$), the evolution for both periods becomes: mi(tR2U2)=Bpiλpλiλpmp0(eλptR2U2eλitR2U2)+mi0eλitR2U2 \begin{equation*} \begin{split} m_i(t_{R_2 \to U_2}) &= B_{pi} \frac{\lambda_p}{\lambda_i - \lambda_p} m^0_p \bigl(e^{\lambda_p t_{R_2 \to U_2}} - e^{\lambda_i t_{R_2 \to U_2}} \bigr) \\&+ m^0_i e^{\lambda_i t_{R_2 \to U_2}} \end{split}\end{equation*}(3) mi(tU2R3)=Bpiλpλiλpmp0(eλptU2R3eλitU2R3)+mi0eλitU2R3 \begin{equation*} \begin{split} m_i(t_{U_2 \to R_3}) &= B_{pi} \frac{\lambda_p}{\lambda_i - \lambda_p} m^0_p \bigl(e^{\lambda_p t_{U_2 \to R_3}} - e^{\lambda_i t_{U_2 \to R_3}} \bigr) \\&+ m^0_i e^{\lambda_i t_{U_2 \to R_3}} \end{split}\end{equation*}(4)

the initial isotope masses are not the same in both equations, because during the reprocessing both parent and daughter nuclei have decayed and/or accumulated, but also have been separated completely or partially during reprocessing. For the mass of isotope i having been produced due to the decay of p during the reprocessing time, a decay factor exp(λi*tR2U2)$exp(-\lambda_i * t_{R_2 \to U_2})$ has been introduced on the produced quantity accounting for the mass that is going to decay during fuel fabrication and by the amount that will be lost in reprocessing. Similarly, during the fuel production, exp(λp*tU2R3)$exp(-\lambda_p * t_{U_2 \to R_3})$ takes into account that during reprocessing, the parent isotope has been depleted due to decay and it might also have been removed during the reprocessing process.

SEPAR is needed to explore in a systematic (on different fleets), rapid (in execution time) and precise (in results) manner, different nuclear reactor fleets at equilibrium. As seen before, the main elements of the dataset are fuel batches in reactor and reprocessing lines, linked together to define a fleet and its cycle. The actual SEPAR outputs consist in the number of reactors in the fleet, and the isotopic compositions of ingoing and outgoing fuels. Although the precision of SEPAR awaits evolutions, this new equilibrium searching tool currently allows the exploration of this domain as a new knowledge domain.

thumbnail Fig. 1

SEPAR working scheme.

2.2 Dynamic code: COSI6

Since 1985, the CEA has been developing the simulation software COSI to study different trajectories of nuclear fleet evolution and provide technical elements to decision makers [3]. The principle of COSI6, including the typical composition of the data set, is exposed in Figure 2.

COSI6 provides a detailed computation of the materials balances, including the computation of Pu content or 235U enrichment entering the fuel fabrication, based on so-called equivalence model. In COSI6, the CESAR code [11] is used for the in-pile calculations. It solves the differential equation systems that describes the fuel evolution in pile (irradiation model) from 1-group cross sections libraries.

In this study, COSI6 is applied for the comparison and validation of the UOX - MOX mono-recycling and UOX - ERU mono-recycling scenarios. The following steps of the SEPAR/COSI6 validation include the simulation of partial and fully recycling of uranium and plutonium fuel strategies. The results of COSI6 constitute a reference regarding the new developed tool for automatic research of fleet equilibrium states.

thumbnail Fig. 2

COSI6 operating diagram.

2.3 Fuel and reactor simulation

The main characteristics of the reactors involved in the following simulations are compiled in Table 1. The majority of fuels used in the simulations are assumed to reach a final burnup of 55 GWd/t, all cores are supposed to have a mass of 129 t, with 4592 MWth of thermal and 1653 MWe of electrical power. The fresh UOx isotopic composition is presented in Table 2. In spentfuel, all uranium isotopes with mass numbers ranging from 232 to 238 are present in varying proportions.

Different from the conventional MOX assembly, MIX assembly [12] is composed of fuel rods with fixed Pu content (8% in our case as mentioned in [13]), completed by enriched uranium support to counter-balance the Pu isotopic degradation. MIX also known as MOX/EUS (MOX Enriched Uranium Support), is studied at CEA within the framework of a quadripartite (Orano, Framatome and EDF) work program in which each partner is a study operator. MIX fuel is a potential method of establishing multiple recycles of plutonium in PWRs without the technical restrictions that affect conventional MOX fuel in this situation, as stated in the French PPE.

Table 1

Main characteristics of the reactors featured in the examples.

Table 2

UOx fresh isotopic composition (heavy metal weight fraction).

3 Result

An equilibrium study can be used to determine the reactor fraction which leads to break-even in regard to plutonium or stabilizationof spent fuel inventory. The reactor fraction is one essential factor which makes a nuclear fleet a net production or consumption of plutonium.

During the simulations, the required power is 61.2 GW of installed capacity, leading to a fleet comprises a total of 37 reactors. The result presented are based on the reference case in which the UOX assemblies are subjected to a cooling time of 10 years and MOX/ERU/MIX fuels are produced and introduced back into the reactor 2 years after the UOX fuel has been reprocessed.

3.1 UOX-MOX mono-recycling

A schematic representation of the UOX - MOX cycle option is shown in Figure 3. SEPAR calculates the batch fractions of the UOX and MOX reactors with the final aim of a plutonium vector stabilization of the MOX fresh fuel. These batch fractions arethen normalized to the corresponding number of nuclear reactors in the fleet, leading to the calculation of the annual natural uranium consumption and the fuel fabrication and reprocessing capacity corresponding to the UOX - MOX plutonium mono-recycling strategy analyzed.

During this exercise, the choice was to compared the main SEPAR outputs with COSI6. A UOX - MOX nuclear fleet has been simulated inCOSI6 using the same SEPAR input of Table 3. The dynamic code COSI6 requires a more detailed characterization of the fuel cycle strategy adopted, as the exact number of each reactor type, the fuel fabrication capacity, the fuel reprocessing capacity etc. Even if these data are introduced in COSI6 as inputs, they are not constant values over the time, because they’re dependent on the material disponibility and requirements. This fluctuations of e.g. annual fluxes or fuel reprocessing capacity is one of the explanations of the discrepancies found between SEPAR and COSI6 results which are related to the Pu241 and Am241 that are investigated more in detail in Section 3.3.

The actual SEPAR outputs for this exercise are presented in Table 3.

This strategy leads to the simulation of equilibrium condition in COSI6, a non-trivial exercise without knowing the SEPAR outputs. The other SEPAR output of interest is the composition of the Pu vector in the MOX fresh fuel, shown in Table 4. This result has been compared with the COSI6 ones, showing a general good agreement. As expected, discrepancies have been found on the Pu vector, and ulterior tests have been done to better understand them. Since SEPAR relies on APOLLO2 calculations and COSI6 on CESAR, a comparison was performed between these two irradiation calculations,showing a general discrepancy in the Pu vector up to 2%. The relative difference column of Table 4 has been calculated taking as reference the COSI6 values for 50 years (this condition has been taken as reference forthis kind of comparison in this work, the statistic validation will be the subject for the next studies).

For what concerns the equilibrium simulation in COSI6, the strategy was to try to have 1 year of fuel fabrication capacity in the spent fuel and therefore to decrease from 10 years to 9 years of cooling, in order to exactly simulate the 10 years cooling as in SEPAR, as can be seen in Figure 4. The same strategy has been adopted for the 2 years aging time, this value was settled to 1 year in COSI6, and the target was to have 1 year of Pu flux availability in the separated Pu stock, as can be seen in Figure 5. The 10.56 tonnes of annual Pu flux have been calculated as the 12% (MOX Pu content) of 88.04 tonnes (MOX fuel fabrication capacity). Since the stabilization of the separated Pu stock is different between SEPAR and COSI6, the 50 year samples for the realization of Table 4 have been taken from 2040 to 2090, in which the average plutonium in the stock is 10.98. The discrepancies between 10.56 and 10.98 means that the COSI6 Pu vector is affected by around 2 weeks longer cooling time per year with respect to the SEPAR one. It explains the discrepancies on 241Pu and 241Am in the Pu vector. The reason is the relatively short half life of 241Pu (t1∕2 = 14.3year), making these two isotopes very sensitive to the aging period.

thumbnail Fig. 3

Fuel cycle of a UOX - MOX fleet.

Table 3

SEPAR output for the UOX-MOX mono-recycling option.

thumbnail Fig. 4

COSI6 UOX spent fuel accumulation.

thumbnail Fig. 5

COSI6 UOX - MOX Pu stock accumulation.

Table 4

COSI6/SEPAR comparison of fresh MOX fuel Pu vector for the UOX-MOX mono-recycling option (heavy metal weight fraction).

3.2 UOX-ERU mono-recycling

A schematic representation of the UOX - ERU cycle option is shown in Figure 6. Generally, for given initial enrichment, 235U content decreases with burnup. The 236U concentration increases with initial fuel enrichment and burnup since it is the product of neutron capture of 235U. 236U is considered to be a neutron poison due to its thermal absorption cross section. The modified isotopic composition of reprocessed uranium poses two major hurdles for its re-use in a reactor: the high radiological impact caused by the presence of 232U prohibits the easy handling of the material, making more radiological protection measures necessary [14], whereas the presence of 236U, in quantities similar to 235U, reduces the macroscopic fission cross section of the fuel, thus creating the need for over-enriching in the 235U isotope. Re-enrichment can happen either directly, in a gas centrifuge plant for example, but another possibility is to mix the reprocessed uranium with already enriched uranium to obtain the desired final enrichment. Alternatively, reprocessed uranium can also be used as a support for MOx fuel, decreasing the needed concentration of plutonium in the fuel due to the higher enrichment with respect to the depleted uranium that is normally used.

The SEPAR simulation converges with 0.899 UOX reactor fraction and 0.101 ERU. For what concerns the equilibrium simulation in COSI6, the strategy on the cooling and fabrication time was the same adopted in Section 3.1. Table 5 shows the comparison between the isotopic vector of reprocessed uranium simulated in COSI6 and SEPARwith the relative differences. The U235 enrichment value for fresh ERU fuels given by SEPAR is 5.235%.

Two explanations can justify the consistent 234U discrepancy: firstly, a 20% APOLLO2 - CESAR relative difference has been found on the 234U concentration after irradiation. Secondly, SEPAR doesn’t consider decay chain longer than one parent and one daughter isotopes. Therefore the 234U SEPAR value is lower because the contribution of 242Cm (that decaysin 238Pu with t1∕2 = 160d and later in 234U with t1∕2 = 87y) is not takeninto account at the moment. Even if this contribution is relatively low, the fact that the 234U contributionof the uranium vector is also small, causes a considerable relative error between these 2 values. This longer decay chain functionality is on top of the future SEPAR development features. These two possible explanations will be investigated in the future in order to understand and reduce this discrepancy.

thumbnail Fig. 6

Fuel cycle of a UOX - ERU fleet.

Table 5

COSI6/SEPAR comparison of reprocessed uranium isotopic vector after irradiation for the UOX-ERU mono-recycling option (heavy metal weight fraction).

3.3 UOX-MIX full multi-recycling

A schematic representation of the UOX - MIX cycle option is shown in Figure 7. In this exercise the MIX fuel used has a discharged burnup of 50 GWday/t and 1404 days of EFPD for the whole irradiation time. All the plutonium from UOX/MIX spent fuel is supposed to be reprocessed and recycled to fabricate new MIX fuel, which helps to stabilize the plutonium inventory. This is a symbiotic equilibrium in the sense that all plutonium produced UOX fuels is net consumed in MIX batches [15,16]. However, there is a risk of producing significant quantities of minor actinides by capture on Pu isotopes (thermal neutron spectrum). The final objective of this exercise is to prove the SEPAR possibility to simulate a plutonium full multi-recycling.

The SEPAR simulation converges with 0.639 UOX reactor fraction and 0.361 MIX. For what concerns the equilibrium simulation in COSI6, the strategy of the cooling time was the same adopted in Section 3.1. To simulate the UOX - MIX equilibrium is a non-trivial exercise because of the complexity introduced by the full multi-recycling condition. In principle, the equilibrium state can be reached by modifying the MIX core fraction given by SEPAR in COSI6 with dichotomy.

The system is difficult to simulate at equilibrium due to the fact that before recycling the MIX fuel, the plutonium extracted from spent fuel is saved in separated Pu stock, as can be seen in Figure 8. Therefore, the separated plutonium inventory should be non-zero constant at equilibrium state, while it oscillates during dynamic simulations.

Regarding the 2 years aging time in SEPAR, the strategy in COSI6 was to reach these 2 years by the sum of the fabrication time and the time the plutonium spends in the separated stock. As can be seen in Figure 8, since the average amount of plutonium in the separated stock is 6 tonnes, to reach the 2 years of aging time, the fabrication time in COSI6 was set to 1.8 years as the result of the following formula: tFabCOSI6=tagingSEPARmPuCOSI6ΦPuCOSI6tFabCOSI6=26/30=1.8y \begin{equation*} \begin{split} t_{Fab_{COSI6}} &= t_{aging_{SEPAR}} - \frac{m_{Pu_{COSI6}}}{\Phi_{Pu_{COSI6}}} \\ t_{Fab_{COSI6}} &= 2 - 6 / 30 = 1.8 y \end{split}\end{equation*}(5)

in which mPuCOSI6$m_{Pu_{COSI6}}$ represents the plutonium average mass in the separated stock.

The other SEPAR result is the composition of the plutonium vector in the MIX fresh fuel, shown in Table 6. This result has been compared with the COSI6 ones, showing a general good agreement. The biggest different found is on the 238Pu. A discrepancy of 2.2% has been already found in Table 4 regarding the UOX - MOX mono-recycling case, therefore this difference can be the result of the error propagation due to several recycling.

thumbnail Fig. 7

Fuel cycle of a UOX - MIX fleet.

thumbnail Fig. 8

COSI6 UOX - MIX Pu stock accumulation.

Table 6

COSI6/SEPAR comparison of fresh MIX fuel Pu vector for the UOX-MIX full multi-recycling option (heavy metal weight fraction).

3.4 UOX-MOX-AMR

A schematic representation of the UOX - MOX - AMR cycle option is shown in Figure 9, where s represent the AMR fuel fraction directed to the reprocessing plant Ret1, while ϵ is the fraction of recycled MOX fuel.

There is currently no recycling strategy implemented for used MOX fuel, which is accumulating in storage. They can’t be entirely multirecycled into new MOX assemblies, because the low plutonium grade of used MOX fuel would require highplutonium content, over the generally accepted upper limit of 12%. A strategy under study is the use of sodium fast reactors, but they are reputed more expensive than their PWR counterpart. A way to reduce their cost currently consists to reduce their size. In this example, a fixed reduced fraction (5% of the installed electrical power) of Small Modular Fast Reactors (called SMR-Na in [17], here called AMR) are added to a UOX – MOX fleet. The AMR core is slightly breeder, thanks to its radial and axial fertile blankets (10.55 and 3.48 t respectively). The main characteristics of the reactors involved in this fleet are compiled in Table 7.

The retreatment plant Ret1 supplies MOX new assemblies in plutonium while Ret2 is dedicated to new AMR fuel assemblies. The fraction s of used AMR fuel, blankets included, is directed toward Ret1, the rest supplies Ret2. MOX used fuel only supplies Ret2. The fabrication and cooling times are 2 and 5 years respectively.

The goal is to maximize the fraction of recycled MOX fuel (ϵ in Fig. 9) while verifying cycle constraints. Thanks to its short calculation time, SEPAR can be used iteratively to solve this sort of optimization problems. The main results of this case are compiled in Tables 8 and 9. At best 70% of used MOX is recycled. The entirety of the used AMR fuel is directed toward MOX fabrication. The plutonium content constraint is not reached, which implies that the identified equilibrium is limited by mass, and not plutonium isotopy. SEPAR would allow for the differentiation of the retreatment of fertile blankets, for example to prioritize them for MOX retreatment. This possibility is not touched on in this paper because the plutonium grade is not the limiting factor.

This last case shows the adaptability of SEPAR to simulate complex systems. The validation of the results concerning the UOX - MOX - AMR cycle has not been tested yet with COSI6, it will be the interest for the next studies.

Table 7

Main characteristics of the reactors featured in the UOX - MOX - AMR example.

thumbnail Fig. 9

Fuel cycle of a AMR - UOX - MOX fleet.

Table 8

UOX - MOX - AMR fleet equilibrium main characteristics.

Table 9

SEPAR fresh fuel Pu vector for the UOX - MOX - AMR option (heavy metal weight fraction).

4 Conclusion

This paper aims to present for the first time the new equilibrium simulation tool SEPAR (Equilibrium Simulator of Advanced Reactor Parks) developed at the Department of Reactor Studies (DER) at CEA Cadarache. It was applied to different fuel cycle options to show how they are modeled and the type of results the SEPAR tool is actually able to produce.

Equilibrium state is the weak point of dynamic scenario codes. Industries are particularly interested in steady state conditions, especially for radioactive stock dimensioning, fabrication and reprocessing facilities, etc. SEPAR can find a unique equilibrium condition (if physically possible) that later can be simulated as a result of a fleet evolution with a dynamic code. Therefore, it can be said that equilibrium codes cannot substitute dynamic codes, while instead they can be seen as complementary tools.

The presentation of the SEPAR outputs is accompanied by the first estimate of validation. The validation work was done for the first three cases by comparing the SEPAR outputs with COSI6 simulations. It has been shown a good agreement for the main SEPAR/COSI6 results. In order to better understand and reduce the discrepancies found between the two methods, new advanced cycle options and further validations of the presented cases will be done. SEPAR is the new tool that allows to study nuclear fleets at equilibrium, precision improvements and additional modules will be added in the future.

Author contribution statement

All the authors were involved in the preparation of the manuscript. All the authors have read and approved the final manuscript.

References

  1. Ministère de la transition écologique et solidaire. Programmation pluriannuelle de l’énergie. fr. p. 400. https://www.ecologie.gouv.fr/sites/default/files/PPE-Executive%20summary.pdf [Google Scholar]
  2. G. Krivtchik, P. Blaise, C. Coquelet-Pascal, Artificial neural network surrogate development of equivalence models for nuclear data uncertainty propagation in scenario studies, EPJ Nucl. Sci. Technol. 3, 22 (2017) [CrossRef] [EDP Sciences] [Google Scholar]
  3. C. Coquelet-Pascal, M. Tiphine, G. Krivtchik, D. Freynet, C. Cany, R. Eschbach, C. Chabert, COSI6: a tool for nuclear transition scenario studies and application to SFR deployment scenarios with minor actinide transmutation, Nucl. Technol. 192, 91 (2015) [CrossRef] [Google Scholar]
  4. D. Minière, Urenco to enrich reprocessed uranium, https://www.world-nuclear-news.org/Articles/Urenco-to-enrich-reprocessed-uranium (2018) [Online; accessed 16-June-2021] [Google Scholar]
  5. G. Martin, M. Tiphine, C. Coquelet-Pascal, French transition scenarios toward a symbiotic nuclear fleet, in International Congress on Advances in Nuclear Power Plants (ICAPP - 2016) (San Francisco, United States, 2016), https://hal-cea.archives-ouvertes.fr/cea-02509781 [Google Scholar]
  6. A. Culot, MOx Strategy and the future of French nuclear plants, https://imtech.wp.imt.fr/en/2019/06/18/mox-strategy-french-nuclear-plants/ (2019), [Online; accessed 17-June-2021] [Google Scholar]
  7. F. Courtin, B. Leniau, N. Thiollière, B. Mouginot, X. Doligez, A. Somaini, A.A. Zakari-Issoufou, S. David, A. Bidaud, J.B. Clavel, Neutronic predictors for PWR fuelled with multi-recycled plutonium and applications with the fuel cycle simulation tool CLASS, Prog. Nucl. Energy 100, 33 (2017) [CrossRef] [Google Scholar]
  8. M. Ernoult, S. David, X. Doligez, O. Meplan, B. Leniau, A. Bidaud, A. Nuttin, B. Mouginot, J. Wilson, N. Capellan et al., Advanced plutonium management in PWR, complementarity of thorium and uranium, Progr. Nucl Energy 78, 330 (2015) [CrossRef] [Google Scholar]
  9. G. Martin, C. Coquelet-Pascal, Symbiotic equilibrium between sodium fast reactors and pressurized water reactors supplied with MOX fuel, Ann. Nucl. Energy 103, 356 (2017) [CrossRef] [Google Scholar]
  10. C. Broyden, A Class of Methods for Solving Nonlinear Simultaneous Equations (1965), http://pi.math.cornell.edu/~web6140/Broyden_965.pdf [Google Scholar]
  11. J.M. Vidal, R. Eschbach, A. Launay, C. Binet, J.F. Thro, in CESAR5.3: an industrial tool for nuclear fuel and waste charactherization with associated qualification (2012), p. 15 [Google Scholar]
  12. G. Youinou, A. Vasile, Plutonium multirecycling in standard PWRs loaded with evolutionary fuels, Nucl. Sci. Eng. 151, 25 (2005) [Google Scholar]
  13. G. Martin, M. Guyot, F. Laugier, G. Senentz, G. Krivtchik, B. Carlier, D. Lecarpentier, F. Descamps, C. Chabert, R. Eschbach, French scenario toward fast plutonium multi-recycling in PWR, in ICAPP, 2018, Charlotte, United States (2018), p. 11 [Google Scholar]
  14. W.B. Arthur, Uranium-232 production in current design LWRS, Tech. rep. (1997), https://www.osti.gov/servlets/purl/5963522 [Google Scholar]
  15. G. Martin, Study of a mixed fleet of breeder SFR and EPR supplied with LEU and MOX fuels to balance the plutonium inventory, in ICAPP, 2018, Charlotte, United States (2018), p. 6 [Google Scholar]
  16. OECD Nuclear Energy Agency, Plutonium management in the medium term: a review by the OECD/NEA working party on the physics of plutonium fuels and innovative fuel cycles (WPPR) (Nuclear Energy Agency, Organisation for Economic Co-operation and Development, Issy-les-Moulineaux, France, 2003), ISBN 978-92-64-02151-8, oCLC: 61895662 [Google Scholar]
  17. K. Tirel, T. Kooyman, C. Coquelet-Pascal, E. Merle, Possibility to stabilize plutonium inventory in a sodium-cooled small modular reactors – Pressurized water reactors fleet, Ann. Nucl. Energy 146, 107632 (2020) [CrossRef] [Google Scholar]

Cite this article as: Heddy Barale, Camille Laguerre, Paul Sabatini, Fanny Courtin, Kévin Tirel and Guillaume Martin, A new tool for the simulation of different nuclear fleets at equilibrium, EPJ Nuclear Sci. Technol. 8, 1 (2022)

All Tables

Table 1

Main characteristics of the reactors featured in the examples.

Table 2

UOx fresh isotopic composition (heavy metal weight fraction).

Table 3

SEPAR output for the UOX-MOX mono-recycling option.

Table 4

COSI6/SEPAR comparison of fresh MOX fuel Pu vector for the UOX-MOX mono-recycling option (heavy metal weight fraction).

Table 5

COSI6/SEPAR comparison of reprocessed uranium isotopic vector after irradiation for the UOX-ERU mono-recycling option (heavy metal weight fraction).

Table 6

COSI6/SEPAR comparison of fresh MIX fuel Pu vector for the UOX-MIX full multi-recycling option (heavy metal weight fraction).

Table 7

Main characteristics of the reactors featured in the UOX - MOX - AMR example.

Table 8

UOX - MOX - AMR fleet equilibrium main characteristics.

Table 9

SEPAR fresh fuel Pu vector for the UOX - MOX - AMR option (heavy metal weight fraction).

All Figures

thumbnail Fig. 1

SEPAR working scheme.

In the text
thumbnail Fig. 2

COSI6 operating diagram.

In the text
thumbnail Fig. 3

Fuel cycle of a UOX - MOX fleet.

In the text
thumbnail Fig. 4

COSI6 UOX spent fuel accumulation.

In the text
thumbnail Fig. 5

COSI6 UOX - MOX Pu stock accumulation.

In the text
thumbnail Fig. 6

Fuel cycle of a UOX - ERU fleet.

In the text
thumbnail Fig. 7

Fuel cycle of a UOX - MIX fleet.

In the text
thumbnail Fig. 8

COSI6 UOX - MIX Pu stock accumulation.

In the text
thumbnail Fig. 9

Fuel cycle of a AMR - UOX - MOX fleet.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.