Research Article |
Corresponding author: Andrei A. Andrianov ( andreyandrianov@yandex.ru ) Academic editor: Yury Kazansky
© 2022 Andrei A. Andrianov, Ilya S. Kuptsov, Tatiana A. Osipova, Anastasia A. Spiridonova, Olga N. Andrianova, Tatiana V. Utianskaya.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Andrianov AA, Kuptsov IS, Osipova TA, Spiridonova AA, Andrianova ON, Utianskaya TV (2022) Nuclear energy system modelling application package: functional overview and examples. Nuclear Energy and Technology 8(4): 289-295. https://doi.org/10.3897/nucet.8.97651
|
The paper provides a brief description of the functionality of the nuclear energy system modelling application package (NESAPP), including a list of modelling objects, key assumptions and areas of application. NESAPP consists of the following main modules: NUDAPS (a module for calculating thermal neutron cross-sections, resonance integrals and one-group neutron cross-sections, and associated uncertainties), NUCLEX (a module for calculating the evolution of the nuclide composition and characteristics of nuclear fuel in reactors and at the nuclear fuel cycle front-end and back-end steps), NUCAB (a module for adjusting isotopic composition and blending), FANES (a module for analysing material flows and integrating data in nuclear energy evolution scenarios), ECNES (a module for assessing economic performance metrics for the nuclear energy evolution scenarios). Each of the modules is a calculation tool that can be used as independent or integrated into the software for technical and economic modelling of nuclear energy systems. Various calculation models are implemented in the modules, allowing users to evaluate the methodological component of the calculation uncertainty in scenario modelling studies and the functionality for assessing the impact of initial data uncertainties on the resulting indicators. The authors also provide some examples of applying NESAPP.
Nuclear energy system, nuclear fuel cycle, technical and economic modelling, scenario analysis, material flow analysis, economic analysis
Software tools for conducting technical and economic modelling of nuclear energy systems (NES) and scenario (system) studies in the field of the nuclear fuel cycle (NFC) make it possible to estimate the change in time of the volume of material flows along the NFC stages, the need for various NFC services and goods, indicators of efficiency and competitiveness of the NES option, taking into account various external and internal factors (the NES refers to the NPP system and the NFC enterprises associated with this system). The corresponding computational models can be constructed at the level of an energy company, a country, a region, or even the whole world and assume a significant representation of the key processing stages of the associated NFC and physical processes that occur with nuclear materials during their circulation in the NES fuel cycle (
Such software tools are used for a comprehensive analysis and feasibility study of possible NES deployment options, optimisation of NES options in terms of a set of performance indicators (taking into account resource and infrastructure constraints), comparison of alternative NES options, NFC organization approaches, approaches to radioactive waste management and selection of the most acceptable ones under given conditions, studies of the impact of technological improvements and new developments on the NES performance indicators.
The corresponding calculation tools integrate submodels of various levels: physical models that make it possible to assess changes in the nuclide composition and characteristics of nuclear fuel in reactors and in the different NFC steps during the NES evolution; models of the dynamics of the movement of nuclear materials in the NFC; models for assessing NES economic efficiency indicators; and market models of NPPs, NFC goods and services. The software tools are used in conjunction with specialised databases containing information on the current state and prehistory of the NES development at the appropriate level (global, regional, country, corporate), as well as the technical and economic parameters of operating and planned reactor plants and NFC facilities.
Despite the existing arsenal of models, approaches, and software tools for systematic studies of the prospects for the development of nuclear energy, systematic work is currently underway on further upgrading and improving them in order to provide the possibility, on their basis, of developing technical and economic NES models with increased accuracy (correct modelling of material flows in the NFC, taking into account nuclide evolution, refinement of estimates of economic indicators, etc.), as well as on translating into computer codes the optimisation and prioritisation mechanisms for overcoming technological forks, the functionality of sensitivity and uncertainty analysis, exploratory data analysis, aggregation of system-wide indicators and multi-criteria comparative evaluation (
The NESAPP software modules are designed for calculation support of studies in the field of nuclear energy planning, scenario modelling of the nuclear energy development, system-level assessment of reactor and NFC technologies both as independent calculation tools and as part of applied software packages for technical and economic modelling of nuclear energy systems. On the basis of the software modules, it is possible to construct complex calculation models for assessing material flows and needs for NFC services (needs for natural uranium, services for uranium conversion and enrichment, fresh fuel fabrication, SNF storage and reprocessing, etc.) and economic performance metrics of NES development scenarios. The main software modules include:
NUDAPS (Nuclear Data Processing Spreadsheets): preparation of one-group neutron cross-sections averaged over standardised and (or) user-specified neutron spectra;
NUCLEX (Nuclide Evolution Explorer): calculation of the evolution of the nuclide composition and characteristics of nuclear fuel in reactors and at the main NFC front-end and back-end steps;
NUCAB (Nuclide Composition Adjustment and Blending Tool): correction of the nuclide composition of the fuel and neutron flux; if necessary, the modelling of the reactor operation on the new isotopic composition of fresh nuclear fuel;
FANES (Material Flow Analysis Data Integrator for Nuclear Energy System): assessment of material flows and needs for NFC goods and services in nuclear energy development scenarios;
ECNES (Economic Assessment Tool for Nuclear Energy System): assessment of economic performance metrics for nuclear energy development scenarios.
The general flow chart of NESAPP is shown in Fig.
Characteristic features of NESAPP are the modular principle of software organization, the autonomy of the modules, the multiplicity of implemented calculation models in each module, the assessment of the uncertainties of the calculation functionals in accordance with the uncertainty of the initial data. In addition, NESAPP is accompanied by a database of typical reactor characteristics, an atlas of one-group neutron cross-sections for various neutron spectra (for nuclide evolution modelling), and a database with data on neutron spectra. The user can form sets of initial data in accordance with the task under consideration. The user is responsible for generating a set of initial data for conducting economic assessments according to the study context and model assumptions implemented in the ECNES module. The assumptions of the calculation models implemented in the NESAPP software modules are specified in the user instructions (
Module | Purpose | Initial data | Assumptions / functional |
---|---|---|---|
NUDAPS | Calculation of one-group neutron cross-sections, various functionals, which are linear combinations, products and ratios of one-group neutron cross-sections and their uncertainties | Evaluated neutron data and covariance matrices, neutron spectrum (these data are placed in a specialised database of evaluated neutron data of actinide nuclei) | The module contains a set of evaluated neutron data and covariance matrices for the main actinide nuclei in the GENDF and BOXER formats (for covariance matrices), a set of standardized neutron spectra (fission, Maxwell, Fermi, etc.) and neutron spectra characteristic of various types of power reactors (the spectra can be specified in both analytical and tabular forms). |
NUCLEX | Calculation of nuclide composition evolution (mass and concentration), activity, radiotoxicity, heat release, neutron and gamma sources of nuclear fuel | Initial isotopic composition of the fuel, one-group neutron fluxes and neutron cross-sections, both determined at various stages of fuel irradiation in the reactor and one-group flux and cross-sections averaged over the entire irradiation cycle and constant in time, specific power density, nuclear fuel burnup | Basic set of actinide nuclei: 234U, 235U, 236U, 238U, 237Np, 238Pu, 239Pu, 240Pu, 241Pu, 242Pu, 241Am, 242mAm, 243Am, 242Cm, 243Cm, 244Cm, 245Cm, 246Cm, 247Cm, 248Cm. Basic set of fission and activation products: 3H, 14C, 36Cl, 54Mn, 55Fe, 60Co, 59Ni, 63Ni, 79Se, 85Kr, 87Rb, 90Sr, 93Mo, 94Nb, 93Zr, 98Tc,99Tc, 102Rh, 106Ru, 107Pd, 110mAg, 113mCd, 121mSn, 125Sb, 126Sn, 125mTe, 127mTe, 126Sn, 129I, 134Cs, 135Cs, 137Cs, 144Ce, 146Pm, 147Pm, 146Sm,147Sm, 151Sm, 152Eu, 154Eu, 155Eu, 166mHo. Output data (masses and concentrations of the main actinides, fission products, radioactivity, decay heat, radiotoxicity, neutron and gamma radiation) are given for the end of the irradiation cycle, after cooling, for the 100th, 1000th and 10000th years. The functionality for estimating the nuclear data and technological uncertainty of the design characteristics is implemented. |
NUCAB | Calculation of the values and uncertainties of the corrections for the nuclear fuel enrichment and the total neutron flux, taking into account the conservation of the fraction of heavy nuclei in the fuel, reactivity and power of the reactor plant | Fuel isotope composition (design and actual), one-group neutron cross-sections, covariance matrices and neutron flux | Various methods for estimating the enrichment and total neutron flux corrections and their uncertainties are implemented in order to select the most appropriate method for its subsequent implementation. The uncertainties of the corrections associated with the initial neutron data are estimated. |
FANES | Calculation of material flows and needs for NFC goods and services for the main stages of the NFC front-end and back-end, assessment of uncertainties in material indicators | Dynamics of commissioning and decommissioning of power units, parameters of fuel use and SNF/RW production, isotope composition of fresh and spent fuel, characteristics of the strategy for handling SNF/RW, fissile materials, principles of prioritisation/management of stocks and reserves of nuclear materials | It is possible to take into account the radioactive decay of nuclides in external part of the NFC when they appear in the NFC in different time periods, simulate delays at various stages of the NFC and in SNF storage facilities, set priorities for the consumption of limited resources and the utilisation of limited capacities of NFC facilities, evaluate the uncertainty in the material indicators caused by the uncertainty in the initial technical data. |
ECNES | Calculation for the NES development scenario of the total discounted costs, net present value, internal rate of return, profitability index, payback period, levelised unit energy cost (LUEC), assessment of the uncertainties in economic indicators due to the uncertainty of the initial cost data | Dynamics of commissioning and decommissioning of power units, technical and economic parameters of power units, material flows and needs for NFC goods and services (needs for natural uranium, uranium conversion and enrichment services, fresh fuel fabrication, SNF storage and processing, etc.), cost data on construction, operation and decommissioning of power units, NFC goods and services, discount rate, electricity tariff | It is possible to consider different distributions of capital costs (uniform, sinusoidal distribution, distribution of a second degree polynomial, arbitrary distribution), adjust capital costs for units that remain in operation beyond the forecast horizon, and account for depreciation deductions during the operation of historical capacities. It is also possible to evaluate the uncertainty of economic indicators caused by the uncertainty of the initial cost data (uniform and triangular distribution of costs, correlated and uncorrelated distributions of input values). |
The NESAPP modules are created using the built-in spreadsheet language – MS Excel Visual Basic for Applications (VBA) (basic option) (Software for Working with MS Excel Spreadsheets).
There are also versions of NESAPP modules developed using the Python programming language, the MathCAD computer algebra system (Engineering Mathematical Software Mathcad), and the Stella Architect system dynamics modelling environment (Stella Architect Modelling Tool) using built-in programming languages. Accordingly, to work with NESAPP, the Microsoft Office application package is required. If the versions of the software modules for MathCAD and Stella Architect are used, the appropriate modelling environments must be available.
Work with NESAPP is carried out through the standard MS Excel interface. If the versions of the software modules for MathCAD and Stella Architect are used, the appropriate engineering graphical interfaces are available. All the programme modules are alienable; the need for author’s support is determined by the user’s qualifications.
The potential user of NESAPP is assumed to have knowledge of the basics of nuclear engineering, the principles of nuclear energy planning and the economics of nuclear energy. The computation time using the models developed with NESAPP is determined by their complexity and can range from a few seconds to several hours.
The software modules were verified using software tools similar in functionality (Modelling Nuclear Energy Systems with MESSAGE, INPRO Methodology, Nuclear Fuel Cycle Simulation System). Identical results were obtained for similar model assumptions and initial data.
NESAPP was comprehensively tested using test problems developed for the purpose of cross-verification of Russian software tools for technical and economic modelling of nuclear energy systems. The cross-verification was carried out in the format of voluntary (initiative) activity in the period of 2020–2021 (
The NESAPP software modules have been used since 2008 for calculation support of studies in the field of nuclear energy planning, scenario modelling of nuclear energy development, assessment of nuclear energy technologies and NFC technologies at the Obninsk Institute for Nuclear Power Engineering of the National Research Nuclear University MEPhI as part of various initiative and contractual research projects (
Using NESAPP, a simulation of scenarios for the partitioning and transmutation of RW generated during the operation of a hypothetical NES with thermal reactors for 1000 years was carried out under two assumptions: the simultaneous appearance of the entire volume of RW (Option 1) and the gradual production of RW in accordance with the scenario of the NES operation (Option 2) (
It is shown that, under comparable conditions, Option 1 (without taking into account the RW production over time) gives an underestimated time for reaching the given reference level by the selected indicator compared to Option 2 (taking into account the RW accumulation over time) (Fig.
NESAPP was used in modelling scenarios for the development of the Russian two-component NES with thermal and fast sodium reactors (
VVER(100%): share of VVER-TOI = 100%;
VVERmox(10%): shares of VVER-TOI, VVER-TOI MOX = 90 and 10%;
VVERmox(30%): shares of VVER-TOI, VVER-TOI MOX = 70 and 30%;
VVERmox(50%): shares of VVER-TOI, VVER-TOI MOX = 50 and 50%;
BN(20%): shares of VVER-TOI, BN = 80 and 20%;
BN(50%): shares of VVER-TOI, BN = 50 and 50%;
BN(90%): shares of VVER-TOI, BN = 10 and 90%;
VVERmox(10%) BN(20%): shares of VVER-TOI, VVER-TOI MOX, BN = 70, 10 and 20%;
VVERmox(50%) BN(20%): shares of VVER-TOI, VVER-TOI MOX, BN = 30, 50 and 20%;
VVERmox(10%) BN(50%): shares of VVER-TOI, VVER-TOI MOX, BN = 40, 10 and 50%.
It is shown that, when the assessment of NES deployment scenarios is based on exclusively of material and economic indicators, the involvement of MOX-fuelled thermal reactors in the two-component NES leads to a decrease in the attractiveness of the corresponding scenario compared to scenarios without MOX fuel in thermal reactors (Table
Scenario (technology share in 2100 is in parentheses) / indicator | Cumulative uranium consumption, kt | Cumulative needs for uranium enrichment services, 106 SWU | Cumulative needs for spent fuel reprocessing services, kt h.m. | Amount of SNF in 2100, kt h.m. | Amount of RW in 2100, kt h.m. | Amount of plutonium in the NFC in 2100, kt | Amount of depleted uranium in 2100, kt | LUEC, mills/kWh |
---|---|---|---|---|---|---|---|---|
VVER(100%) | 787.65 | 666.75 | 0 | 126.98 | 0 | 1.09 | 1669.03 | 29.48 |
VVERmox(10%) | 782.92 | 662.72 | 27.36 | 106.12 | 26.92 | 0.84 | 1658.34 | 29.89 |
VVERmox(30%) | 776.18 | 656.97 | 71.81 | 69.11 | 70.76 | 0.54 | 1644.80 | 30.40 |
VVERmox(50%) | 772.07 | 653.46 | 92.71 | 51.06 | 91.42 | 0.41 | 1638.11 | 30.57 |
BN(20%): | 658.10 | 556.25 | 12.15 | 111.54 | 10.96 | 1.02 | 1544.24 | 29.53 |
BN(50%): | 492.27 | 414.81 | 41.77 | 77.70 | 38.93 | 0.91 | 1384.26 | 30.22 |
BN(90%): | 284.97 | 237.99 | 100.28 | 13.96 | 95.37 | 0.76 | 1184.07 | 31.09 |
VVERmox(10%) BN(20%) | 651.51 | 550.63 | 346.08 | 81.58 | 46.76 | 0.77 | 1531.93 | 30.29 |
VVERmox(50%) BN(20%) | 640.97 | 541.64 | 126.27 | 14.46 | 123.68 | 0.34 | 1511.78 | 31.57 |
VVERmox(10%) BN(50%) | 470.99 | 396.66 | 79.63 | 44.59 | 76.26 | 0.65 | 1359.36 | 31.70 |
NESAPP was used to estimate the levelised unit energy cost (LUEC) at NPPs with different types of nuclear reactors, with account taken of the uncertainty in the cost data for the construction and operation of reactor plants as well as for the NFC goods and services (Table
Type of NES | LUEC and its components, mills/kW∙h | |||
---|---|---|---|---|
LUAC | LUOM | LUFC | LUEC | |
LWR* | 37.6±9.3 | 12.2±1.2 | 8.9±2.3 | 58.7±9.6 |
LWR** | 37.6±9.5 | 12.2±1.2 | 9.9±2.3 | 59.8±9.8 |
HWR* | 36.2±9.1 | 12.2±1.3 | 10.6±2.0 | 59.2±9.3 |
HWR** | 36.2±9.1 | 12.2±1.2 | 17.6±3.2 | 65.8±9.8 |
ALWR* | 35.4±8.8 | 11.6±1.2 | 6.7±1.5 | 53.8±9.0 |
ALWR** | 35.3±8.7 | 11.6±1.2 | 7.5±1.5 | 54.5±9.1 |
FR-1 | 40.2±12.0 | 12.5±1.2 | 13.5±2.4 | 66.2±12.4 |
FR-2 | 40.0±12.1 | 12.5±1.2 | 14.2±2.4 | 67.3±12.4 |
AFR | 40.4±12.1 | 12.6±1.2 | 8.9±1.5 | 61.4±12.1 |
FR-U | 40.4±12.0 | 12.6±1.2 | 12.9±2.0 | 65.7±12.2 |
Three types of thermal reactors (LWR, HWR and ALWR) operating in an once-through NFC and four variants of fast reactors (FR-1, FR-2, AFR and FR-U) operating in a closed uranium-plutonium NFC were considered (data for the relevant technologies were taken from the international IAEA INPRO GAINS project, cost data were from a report prepared as part of one of the studies in the interests of the US Department of Energy (
LWR is a PWR type reactor with a burnup of 45 GW∙day/t h.m., fuel enrichment of 4%, a specific core energy density of 38.5 MW/t, and a capacity factor of 85%.
ALWR is an advanced LWR with higher fuel burnup; compared to LWR, equilibrium fuel loading in ALWR is 30% less than in LWR, initial fuel enrichment is 3.4%, equilibrium fuel enrichment is 4.95%.
HWR is a pressurised heavy water nuclear reactor with a burnup of 7 GW∙day/t h.m., fuel based on natural uranium, a specific power density of 24.0 MW/t, and a capacity factor of 85%.
FR-1 is a sodium-cooled fast reactor with a breeding factor close to unity and an average fuel burnup (core and blankets) of about 38 GW∙day/t h.m.
FR-2 is a prototype sodium-cooled fast breeder reactor with an average breeding factor of 1.16 and an average fuel burnup (core and blankets) of about 31 GW∙day/t h.m.
AFR is a commercial sodium-cooled fast reactor with an average breeding factor of 1.2 and an average (core and blankets) fuel burnup of 54 GW∙day/t h.m. Unlike FR-1 and FR-2, fresh AFR fuel contains about 1% minor actinides (MA).
FR-U is a lead-cooled fast reactor that uses enriched uranium fuel (enrichment of about 15% in 235U) for initial core loading and first refuelling; subsequent SNF reprocessing and the use of secondary nuclear fuel (Pu + U + MA) are assumed. The blanket is not provided. The breeding factor is 1.05, fuel burnup is about 72.8 GW∙day/t h.m.
Based on the results of assessing the levelised unit energy cost and with account taken of the uncertainty in its values due to the spread in cost data, it can be concluded that it is impossible to make an unambiguous judgment about the greatest attractiveness of a particular concept of fast reactors, relying only on the analysis of the levelised unit energy cost, and it is also incorrect to make categorical statements about the lower economic efficiency and competitiveness of fast reactors compared to thermal reactors.
The nuclear energy system modelling application package (NESAPP) is intended for calculation support of scenario studies and technical and economic modelling of nuclear energy systems. The performed verification and cross-verification calculations using similar in functionality software tools and analogous model assumptions have demonstrated the correctness of the methods and algorithms implemented in the NESAPP. The experience of using NESAPP allows us to conclude that the developed toolkit is an effective means of supporting studies in the field of nuclear energy planning and scenario modelling, especially in the cases where it is necessary to assess the impact of uncertainties in the initial data and model assumptions on the calculation results.