Kostadin Ivanov

Improved cross-section modeling methodology for coupled three-dimensional transient simulations

Abstract The purpose of this paper is to introduce a more accurate and sophisticated methodology for use in three-dimensional coupled neutronic/thermal hydraulic analysis. The approach described in this paper is an original method of modeling cross-section variations for off-nominal core conditions, which is becoming an important issue with the increased use of coupled three-dimensional neutronic/thermal hydraulic simulations. This proposed method improves the accuracy of the cross-section modeling for transient applications and it is called the adaptive high-order table lookup method (AHTLM). During nuclear power plant (NPP) transient and accident simulations AHTLM interpolates into multi-dimensional cross-sections tables, which form a box envelope bounding the expected range of change of both nominal and off-nominal NPP conditions. This paper further addresses the methodologies for the development of the cross-section libraries and issues that affect the proper formulation of accurate data. The automated generation procedure outlined in this paper gives the user the tools and the ability of generating accurate cross-sections that cover a large range of thermal hydraulic parameters. Further improvements and expansions for future applications are also discussed.

Nodal kinetics model upgrade in the Penn State coupled TRAC/NEM codes

Abstract The Pennsylvania State University currently maintains and does development and verification work for its own versions of the coupled three-dimensional kinetics/thermal-hydraulics codes TRAC-PF1/NEM and TRAC-BF1/NEM. The subject of this paper is nodal model enhancements in the above mentioned codes. Because of the numerous validation studies that have been performed on almost every aspect of these codes, this upgrade is done without a major code rewrite. The upgrade consists of four steps. The first two steps are designed to improve the accuracy of the kinetics model, based on the nodal expansion method. The polynomial expansion solution of 1D transverse integrated diffusion equation is replaced with a solution, which uses a semi-analytic expansion. Further the standard parabolic polynomial representation of the transverse leakage in the above 1D equations is replaced with an improved approximation. The last two steps of the upgrade address the code efficiency by improving the solution of the time-dependent NEM equations and implementing a multi-grid solver. These four improvements are implemented into the standalone NEM kinetics code. Verification of this code was accomplished based on the original verification studies. The results show that the new methods improve the accuracy and efficiency of the code. The verification of the upgraded NEM model in the TRAC-PF1/NEM and TRAC-BF1/NEM coupled codes is underway.

Methodology of internal assessment of uncertainty and extension to neutron kinetics/thermal-hydraulics coupled codes

Abstract The best-estimate calculation results from complex system codes are affected by approximations that are unpredictable without the use of computational tools that account for the various sources of uncertainty. The code with (the capability of) internal assessment of uncertainty (CIAU) has been previously proposed by the University of Pisa to realize the integration between a qualified system code and an uncertainty methodology and to supply proper uncertainty bands each time a nuclear power plant (NPP) transient scenario is calculated. The derivation of the methodology and the results achieved by the use of CIAU are discussed to demonstrate the main features and capabilities of the method. In a joint effort between the University of Pisa and The Pennsylvania State University, the CIAU method has been recently extended to evaluate the uncertainty of coupled three-dimensional neutronics/thermal-hydraulics calculations. The result is CIAU-TN. The feasibility of the approach has been demonstrated, and sample results related to the turbine trip transient in the Peach Bottom NPP are shown. Notwithstanding that the full implementation and use of the procedure requires a database of errors not available at the moment, the results give an idea of the errors expected from the present computational tools.

Analysis of the Peach Bottom turbine trip 2 experiment by coupled RELAP5-PARCS three-dimensional codes

Abstract Thanks to continuous progress in computer technology, it is now possible to perform best-estimate simulations of complex scenarios in nuclear power plants. This method is carried out through the coupling of three-dimensional (3-D) neutron modeling of a reactor core into system codes. It is particularly appropriate for transients that involve strong interactions between core neutronics and reactor loop thermal hydraulics. For this purpose, the Peach Bottom boiling water reactor turbine trip test was selected to challenge the capability of such coupled codes. The test is characterized by a power excursion induced by rapid core pressurization and a self-limiting course behavior. In order to perform the closest simulation, the coupled thermal-hydraulic system code RELAP5 and 3-D neutron kinetic code PARCS were used. The obtained results are compared to those available from experimental data. Overall, the coupled code calculations globally predict the most significant observed aspects of the transient, such as the pressure wave amplitude across the core and the power course, with an acceptable agreement. However, sensitivity studies revealed that more-accurate code models should be considered in order to better match the void dynamic and the cross-section variations during transient conditions.

Automated Design and Optimization of Pebble-bed Reactor Cores

Abstract This paper presents a conceptual design approach for high-temperature gas-cooled reactors using recirculating pebble bed cores. The method employs PEBBED, a reactor physics code specifically designed to solve for the asymptotic burnup state of pebble bed reactors in conjunction with a genetic algorithm to obtain a core with acceptable properties. The uniqueness of the asymptotic core state and the small number of independent parameters that define it suggest that core geometry and fuel cycle can be efficiently optimized toward a specified objective. A novel representation of the distribution of pebbles enables efficient coupling of the burnup and neutron diffusion solvers. Complex pebble recirculation schemes can be expressed in terms of a few parameters that are amenable to manipulation using modern optimization techniques. The user chooses the type and range of core physics parameters that represent the design space. A set of traits, each with acceptable and preferred values expressed by a simple fitness function, is used to evaluate the candidate reactor cores. The stochastic search algorithm automatically drives the generation of core parameters toward the optimal core as defined by the user. For this study, the design of two pebble bed high-temperature reactor concepts subjected to demanding physical constraints demonstrated the technique’s efficacy.

OECD/NRC BWR Turbine Trip Transient Benchmark as a Basis for Comprehensive Qualification and Studying Best Estimate Coupled Codes

Abstract An Organisation for Economic Co-operation and Development (OECD)/U.S. Nuclear Regulatory Commission (NRC)-sponsored coupled-code benchmark has been initiated for a boiling water reactor (BWR) turbine trip (TT) transient. Turbine trip transients in a BWR are pressurization events in which the coupling between core space-dependent neutronic phenomena and system dynamics plays an important role. In addition, the available real plant experimental data make this benchmark problem very valuable. Over the course of defining and coordinating the BWR TT benchmark, a systematic approach has been established to validate best-estimate coupled codes. This approach employs a multilevel methodology that not only allows for a consistent and comprehensive validation process but also contributes to the study of different numerical and computational aspects of coupled best-estimate simulations. This paper provides an overview of the OECD/NRC BWR TT benchmark activities with emphasis on the discussion of the numerical and computational aspects of the benchmark.

Features and performance of a coupled three-dimensional thermal-hydraulic/kinetics TRAC-PF1/NEM pressurized water reactor (PWR) analysis code

Abstract This paper summarizes the current status of the Pennsylvania State University (PSU) version of the coupled three-dimensional (3-D) thermal-hydraulic/kinetics TRAC-PF1/NEM code for pressurized water reactor (PWR) transient and accident analysis and describes applications to reactivity insertion accident (RIA) simulations as well as recent developments. The TRAC-PF1/NEM methodology utilizes closely coupled 3-D thermal-hydraulics and 3-D core neutronics transient models to simulate the vessel and a 1-D simulation of the primary system. An efficient and flexible cross-section generation procedure was developed and implemented into TRAC-PF1/NEM. These features make the coupled code capable of modeling PWR reactivity transients, including boron dilution transients, in a reasonable amount of computer time. Three-dimensional studies on hot zero power (HZP) rod ejection and main steam line break (MSLB) transients in a PWR, as well as a large break loss-of-coolant-accident (LBLOCA) and boron dilution transients, were accomplished using TRAC-PF1/NEM. The results obtained demonstrate that this code is appropriate for analysis of the space-dependent neutronics and thermal-hydraulic coupled phenomena related to most current safety issues.

OECD/NEA BENCHMARK FOR UNCERTAINTY ANALYSIS IN MODELING (UAM) FOR LWRS – SUMMARY AND DISCUSSION OF NEUTRONICS CASES (PHASE I)

A Nuclear Energy Agency (NEA), Organization for Economic Co-operation and Development (OECD) benchmark for Uncertainty Analysis in Modeling (UAM) is defined in order to facilitate the development and validation of available uncertainty analysis and sensitivity analysis methods for best-estimate Light water Reactor (LWR) design and safety calculations. The benchmark has been named the OECD/NEA UAM-LWR benchmark, and has been divided into three phases each of which focuses on a different portion of the uncertainty propagation in LWR multi-physics and multi-scale analysis. Several different reactor cases are modeled at various phases of a reactor calculation. This paper discusses Phase I, known as the “Neutronics Phase”, which is devoted mostly to the propagation of nuclear data (cross-section) uncertainty throughout steady-state stand-alone neutronics core calculations. Three reactor systems (for which design, operation and measured data are available) are rigorously studied in this benchmark: Peach Bottom Unit 2 BWR, Three Mile Island Unit 1 PWR, and VVER-1000 Kozloduy-6/Kalinin-3. Additional measured data is analyzed such as the KRITZ LEU criticality experiments and the SNEAK-7A and 7B experiments of the Karlsruhe Fast Critical Facility. Analyzed results include the top five neutronnuclide reactions, which contribute the most to the prediction uncertainty in keff, as well as the uncertainty in key parameters of neutronics analysis such as microscopic and macroscopic cross-sections, six-group decay constants, assembly discontinuity factors, and axial and radial core power distributions. Conclusions are drawn regarding where further studies should be done to reduce uncertainties in key nuclide reaction uncertainties (i.e.: 238U radiative capture and inelastic scattering (n, n’) as well as the average number of neutrons released per fission event of 239Pu).

Validation of Coupled Thermal-Hydraulic and Neutronics Codes for Safety Analysis by International Cooperations

Incorporating full three-dimensional models of the reactor core into system transient codes allows for a “best-estimate” calculation of interactions between the core behavior and plant dynamics. Considerable efforts have been made in various countries and organizations on the development of coupled thermal-hydraulic and neutronics codes. Appropriate benchmarks have been developed in international cooperations led by the Nuclear Energy Agency (NEA) of the Organization for Economic Cooperation and Development (OECD) that permit testing of the neutronics–thermal-hydraulics coupling and verification of the capability of the coupled codes to analyze complex transients with coupled core-plant interactions. Three such benchmarks are presented in this paper—the OECD/U.S. Nuclear Regulatory Commission (NRC) pressurized water reactor main steam line break benchmark, the OECD/NRC boiling water reactor turbine trip benchmark, and the OECD/U.S. Department of Energy/Commissariat à l’Energie Atomique V1000 coolant transient benchmark. To meet the objectives of the validation of best-estimate coupled codes, a systematic approach has been introduced to evaluate the analyzed transients employing a multilevel methodology. Since these benchmarks are based on both code-to-code and code-to-data comparisons, further guidance for presenting and evaluating results has been developed. During the course of the benchmark activities, a professional community has been established, which allowed our carrying out in-depth discussions of different aspects considered in the validation process of the coupled codes. This positive output has certainly advanced the state of the art in the area of coupling research.

Optimum Discharge Burnup and Cycle Length for PWRs

Abstract This paper discusses the results of a pressurized water reactor fuel management study determining the optimum discharge burnup and cycle length. A comprehensive study was performed considering 12-, 18-, and 24-month fuel cycles over a wide range of discharge burnups. A neutronic study was performed followed by an economic evaluation. The first phase of the study limited the fuel enrichments used in the study to <5 wt% 235U consistent with constraints today. The second phase extended the range of discharge burnups for 18-month cycles by using fuel enriched in excess of 5 wt%. The neutronic study used state-of-the-art reactor physics methods to accurately determine enrichment requirements. Energy requirements were consistent with today’s high capacity factors (>98%) and short (15-day) refueling outages. The economic evaluation method considers various component costs including uranium, conversion, enrichment, fabrication and spent-fuel storage costs as well as the effect of discounting of the revenue stream. The resulting fuel cycle costs as a function of cycle length and discharge burnup are presented and discussed. Fuel costs decline with increasing discharge burnup for all cycle lengths up to the maximum discharge burnup considered. The choice of optimum cycle length depends on assumptions for outage costs.