Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel performance codes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

Advancing the Fork detector for quantitative spent nuclear fuel verification

The Fork detector is widely used by the safeguards inspectorate of the European Atomic Energy Community (EURATOM) and the International Atomic Energy Agency (IAEA) to verify spent nuclear fuel. Fork measurements are routinely performed for safeguards prior to dry storage cask loading. Additionally, spent fuel verification will be required at the facilities where encapsulation is performed for acceptance in the final repositories planned in Sweden and Finland. The use of the Fork detector as a quantitative instrument has not been prevalent due to the complexity of correlating the measured neutron and gamma ray signals with fuel inventories and operator declarations. A spent fuel data analysis module based on the ORIGEN burnup code was recently implemented to provide automated real-time analysis of Fork detector data. This module allows quantitative predictions of expected neutron count rates and gamma units as measured by the Fork detectors using safeguards declarations and available reactor operating data. This study describes field testing of the Fork data analysis module using data acquired from 339 assemblies measured during routine dry cask loading inspection campaigns in Europe. Assemblies include both uranium oxide and mixed-oxide fuel assemblies. More recent measurements of 50 spent fuel assemblies at the Swedish Central Interim Storage Facility for Spent Nuclear Fuel are also analyzed. An evaluation of uncertainties in the Fork measurement data is performed to quantify the ability of the data analysis module to verify operator declarations and to develop quantitative go/no-go criteria for safeguards verification measurements during cask loading or encapsulation operations. The goal of this approach is to provide safeguards inspectors with reliable real-time data analysis tools to rapidly identify discrepancies in operator declarations and to detect potential partial defects in spent fuel assemblies with improved reliability and minimal false positive alarms. Finally, the results are summarized, and sources and magnitudes of uncertainties are identified, and the impact of analysis uncertainties on the ability to confirm operator declarations is quantified.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Refinements of a treat model for transient analysis

In order to prepare for the TREAT reactor being taken off standby, the Gateway for Advanced Modeling and Simulation Lab has used a set of KENO models to simulate certain TREAT temperature- limited transients using the Transient Reactor excursion simulator (T-ReX) code, a hybrid deterministic-stochastic code. These transients are from the M8CAL experiment series, a set of experiments done to calibrate the reactor for the planned M8 series of fuel tests. Recently, the model has been updated in media definitions and geometry. Models were examined and calculations were performed to low uncertainty values to find the k{sub eff} and reactivity insertions. Alterations to modeled values of shutdown and initial transient rod insertion depth were made to create models which had experimental values of reactivity insertion, pre-transient k{sub eff}, or both. It was found that two of the altered media inputs, fuel and Zircaloy-3 cladding, had significant effect on the k{sub eff}. Surprisingly, the Zircaloy-3 composition had a larger effect on k{sub eff}. The reactivity insertions in this model's 2855, 2856 and 2857 temperature-limited transients were found to be above experimental values by a factor of 5.5% to 7%, and that pre-transient k{sub eff} was too high, at 1.002. It was found that increasing shutdown rod insertion by 2 to 2.5 cm decreased pre-transient k{sub eff} to 1.0, and that increasing shutdown rod insertion by 3 to 5 cm and decreasing initial transient rod insertion by 1 to 2 cm gave perfect pre-transient k{sub eff} and total reactivity insertion. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Advancing Understanding of Fission Gas Behavior in Nuclear Fuel through Leadership Class Computing

This report describes the outcome of the pilot SciDAC project jointly funded by the DOE Office of Nuclear Energy and the DOE Office of Science, Advanced Scientific Computing Research to advance understanding of fission gas behavior in nuclear fuel. The computational approach involved utilizing both a bottom up and top down multiscale materials modeling to evaluate the processes of xenon fission gas diffusion, xenon clustering and precipitation and the mechanisms of fission gas resolution by energetic fission fragments in uranium dioxide nuclear fuel. The computer science and applied mathematics research activities were focused on optimizing the uranium – oxygen – xenon interatomic potential implementation in the LAMMPS molecular dynamics code for large-scale simulations of fission gas resolution, in addition to improvements in visualization of accelerated molecular dynamics simulations and performance optimization of the Xolotl continuum cluster dynamics with a new grouping scheme to reduce the computational cost of modeling very large fission gas bubbles containing more than 10,000 xenon atoms. The most impactful result from this project stems from a systematic series of atomistic molecular dynamics to evaluate fission gas resolution from xenon bubbles that, combined with the fission product yields and consideration of electron-phonon coupling, informed a first principles model that predicts the magnitude of fission gas resolution in good agreement with empirical estimates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Blended Wing Body (BWB) Boundary Layer Ingestion (BLI) Inlet Configuration and System Studies

A study was conducted to determine the potential reduction in fuel burned for BLI (boundary layer ingestion) inlets on a BWB (blended wing body) airplane employing AFC (active flow control). The BWB is a revolutionary type airplane configuration with engines on the aft upper surface where thick boundary layer offers the greatest opportunity for ram drag reduction. AFC is an emerging technology for boundary layer control. Several BLI inlet configurations were analyzed in the NASA-developed RANS Overflow CFD code. The study determined that, while large reductions in ram drag result from BLI, lower inlet pressure recovery produces engine performance penalties that largely offset this ram drag reduction. AFC could, however, enable a short BLI inlet that allows surface mounting of the engine which, when coupled with a short diffuser, would significantly reduce drag and weight for a potential 10% reduction in fuel burned. Continuing studies are therefore recommended to achieve this reduction in fuel burned considering the use of more modest amounts of BLI coupled with both AFC and PFC (Passive Flow Control) to produce a fail-operational system.

Kawai, Ronald T.↗

Development of a Smart Signal Detection Method for Cyclic Voltammetry via Artificial Neural Intelligence

The electrorefiner (ER) is the heart of pyroprocessing technology which is a high-temperature method for separating uranium from Used Nuclear Fuel (UNF). It is important to improve this technology with respect to nuclear materials detection and accountability. Artificial Neural Intelligence (ANI) is a novel data analysis and simulation method that can be applied to electrochemical data sets. This computational code, which has been performed using the commercial software MATLAB, can be trained to generalize adequate electrical current and potential simulated data sets for the unseen data with a high accuracy of prediction. For this purpose, a massive collection of cyclic voltammetry (CV) data sets by Hoover (2014), for 0.5, 1, 2.5, and 5 wt% of zirconium chloride in LiCl-KCl molten salt with different scan rates at 773 K has been considered. The computer is trained via ANI to predict the unseen data after providing suitable hidden layers and validation numbers. In addition, this work can trace the CV plot for a blind condition by interpolating between two simulated data set. The different hidden layers with various neurons (from 5 to 30) at several validation numbers (from 5 to 30) has been studied and the average percent error between experimental and theoretical data for 0.5 wt% with 200 and 450 mV/s has been calculated. Preliminary results demonstrate that if the number of hidden layers increases from one to three, the average error falls down from 44% to around 8%. The best condition which gives a minimum average percent error has been discovered and the simulated CV graph has been compared with the experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Study of reference burnup steps optimization in fuel segment data file generation for NEXUS/ANC9 code system

For any two-step core design code system, the cross-section files are the primary factor to determine the accuracy of the system prediction. For a once-through cross-section system to cover all potential applicable conditions, the system may need to perform tens of thousands state points of lattice calculations. These calculations take a significant amount of CPU times and computing power, especially when a fine or ultra-fine energy group library is used. As computer power increases, so does the calculational complexity. Therefore, it is important to make the lattice calculations effective and efficient. Based on the Westinghouse core design code system NEXUS/ANC9, this work focuses on the reference burnup steps used by the lattice code when performing the calculations. Following the fundamental cross-section methodology and analyzing the contribution of each individual terms, this study provides an applicable solution of the optimized reference burnup steps. The test results show that, with the existing cross-section methodology, it is possible to significantly reduce the lattice calculation cases and cross-section file generation time without sacrificing the accuracy of system prediction. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Benchmark for evaluation and validation of reactor simulations (BEAVRS)

Advances in parallel computing have made possible the development of high-fidelity tools for the design and analysis of nuclear reactor cores, and such tools require extensive verification and validation. This paper introduces BEAVRS, a new multi-cycle full-core Pressurized Water Reactor (PWR) depletion benchmark based on two operational cycles of a commercial nuclear power plant that provides a detailed description of fuel assemblies, burnable absorbers, in-core fission detectors, core loading patterns, and numerous in-vessel components. This benchmark enables analysts to develop extremely detailed reactor core models that can be used for testing and validation of coupled neutron transport, thermal-hydraulics, and fuel isotopic depletion. The benchmark also provides measured reactor data for Hot Zero Power (HZP) physics tests, boron letdown curves, and three-dimensional in-core flux maps from fifty-eight instrumented assemblies. Initial comparisons between calculations performed with MIT's OpenMC Monte Carlo neutron transport code and measured cycle 1 HZP test data are presented, and these results display an average deviation of approximately 100 pcm for the various critical configurations and control rod worth measurements. Computed HZP radial fission detector flux maps also agree reasonably well with the available measured data. All results indicate that this benchmark will be extremely useful in validation of coupled-physics codes and uncertainty quantification of in-core physics computational predictions. The detailed BEAVRS specification and its associated data package is hosted online at the MIT Computational Reactor Physics Group web site (http://crpg.mit.edu/), where future revisions and refinements to the benchmark specification will be made publicly available. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Comparison of silver release predictions using PARFUME with results from the AGR-2 irradiation experiment

The PARFUME (PARticle FUel ModEl) code was used to predict silver release from tristructural isotropic (TRISO) coated fuel particles and compacts during the second irradiation experiment (AGR-2) of the Advanced Gas Reactor Fuel Development and Qualification program. The PARFUME model for the AGR-2 experiment used the fuel compact volume average temperature for each of the 559 days of irradiation to calculate the release of fission product silver from a representative particle for a select number of AGR-2 compacts and individual fuel particles containing either mixed uranium carbide/oxide (UCO) or 100% uranium dioxide (UO2) kernels. Post-irradiation examination (PIE) measurements were performed to provide data on release of silver from these compacts and individual fuel particles. The available experimental fractional releases of silver were compared to their corresponding PARFUME predictions. Preliminary comparisons show that PARFUME under-predicts the PIE results in UCO compacts and is in reasonable agreement with experimental data for UO2 compacts. The accuracy of PARFUME predictions is impacted by the code limitations in the modeling of the temporal and spatial distributions of the temperature across the compacts. Nevertheless, the comparisons on silver release lie within the same order of magnitude.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

CELCAP: A Computer Model for Cogeneration System Analysis

A description of the CELCAP cogeneration analysis program is presented. A detailed description of the methodology used by the Naval Civil Engineering Laboratory in developing the CELCAP code and the procedures for analyzing cogeneration systems for a given user are given. The four engines modeled in CELCAP are: gas turbine with exhaust heat boiler, diesel engine with waste heat boiler, single automatic-extraction steam turbine, and back-pressure steam turbine. Both the design point and part-load performances are taken into account in the engine models. The load model describes how the hourly electric and steam demand of the user is represented by 24 hourly profiles. The economic model describes how the annual and life-cycle operating costs that include the costs of fuel, purchased electricity, and operation and maintenance of engines and boilers are calculated. The CELCAP code structure and principal functions of the code are described to how the various components of the code are related to each other. Three examples of the application of the CELCAP code are given to illustrate the versatility of the code. The examples shown represent cases of system selection, system modification, and system optimization.

Source record↗

Startup Test Plan and Predictions for Highly Enriched Uranium to Low-Enriched Uranium Fuel Conversion at the University of Missouri Research Reactor

Nonpower reactors licensed by the U.S. Nuclear Regulatory Commission require a startup test plan as part of any facility modification to verify operability prior to resumption of operations. In order to support conversion of the University of Missouri Research Reactor from the use of highly enriched uranium to low-enriched uranium (LEU) fuel, a startup test plan has been devised to measure certain reactor physics parameters for the initial all-fresh LEU core licensing documentation that will be submitted. These parameters include the approach to critical, primary coolant void coefficient of reactivity, flux trap void coefficient of reactivity, determination of flux trap sample reactivity worth, radial and axial thermal neutron flux mapping, control blade worth calibration, primary and pool coolant temperature coefficient of reactivity, and flux mapping of experimental positions. Here, predictions for these parameters made using the Monte Carlo N-Particle Version 5 (MCNP5) radiation transport code are reported. These predictions will support the startup tests by providing a baseline set of expectations and additional insight into the performance of the LEU core.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis with the ORIGEN Module of PNAR Spent Fuel Measurements for Nuclear Safeguards Applications

This report documents the analysis of the Passive Neutron Albedo Reactivity (PNAR) measurements for 23 boiling water reactor spent fuel assemblies that were performed in Finland under Action Sheet 65, which is an international collaboration on spent fuel safeguards verification methods in the context of the Finnish spent fuel encapsulation/repository system. PNAR measures the passive neutron and gamma emission rates from each of the spent fuel assemblies like a Fork detector, and it also measures the PNAR ratio, which is expected to correlate with the net neutron multiplication of the measured fuel assembly. The analysis was performed with the Oak Ridge Isotope Generation and Depletion (ORIGEN) Data Analysis Module, which was originally developed for predicting Fork detector spent fuel measurement signals in real time and has been integrated into the Integrated Review and Analysis Package developed by Euratom and the International Atomic Energy Agency. The module includes the ORIGEN burnup analysis code and integrates detector response functions pregenerated using the MCNP Monte Carlo transport code to predict the detector signals in several seconds per assembly. In this study, new response functions specific to the analyzed PNAR measurements were generated for the ORIGEN Module. The study also analyzes the impact of using detailed fuel design and operation information vs. standard safeguards information on the results calculated with the ORIGEN Module. Using detailed information reduced the standard deviation of the relative differences between calculated and measured neutron count rates among the 23 assemblies from ~10% to ~4%. The results obtained using standard safeguards information for these PNAR measurements were similar to those obtained for the Fork detector. A clear trend was found between the calculated net neutron multiplications and the measured PNAR ratios of the 23 assemblies. An uncertainty assessment of the calculations was also performed to estimate the potential impact of the accuracy and the completeness of declaration information on the calculated results for the count rates and net neutron multiplications. The ORIGEN Module predictions for PNAR signals and net neutron multiplication are expected to directly support the safeguards inspectors’ efforts to verify operator declarations of a spent fuel assembly in real time in safeguards practices.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

1025:1 Area Ratio Nozzle Evaluated at High Combustion Chamber Pressure

A recently completed experimental test program obtained performance data on an optimally contoured nozzle with an exit-to-throat area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 12.4 to 16.5 mPa (1800 to 2400 psia). Testing was conducted in the altitude test capsule at the NASA Lewis Research Center's Rocket Engine Test Facility (RETF), and results were compared with analytical performance predictions. This testing builds on previous work with this nozzle at Lewis, where testing was completed at nominal chamber pressure of 350 psia. High-area-ratio nozzles have long been sought as a means to increase the performance of spacebased rocket engines. However, as the area ratio increases, the physical size and weight of the nozzle also increase. As a result, engine and vehicle designers must make tradeoffs between nozzle size and performance enhancement. Until this test program, very little experimental data existed on the performance of the high-area-ratio nozzles used in rocket engine designs. The computer codes being used by rocket engine designers rely on data extrapolated from tests of low-area-ratio nozzles, and these extrapolations do not always provide the accuracy needed for a reliable design assessment. Therefore, we conducted this high-area-ratio nozzle testing program to provide performance data for use in rocket engine design and analysis computer codes. The nozzle had a nominal 2.54-cm- (1-in.-) diameter throat, an exit diameter of 81.3-cm (32.0-in.) at an exit-to-throat area ratio of 1025, and a length of 128.6 cm (50.6 in.). Testing was conducted in an altitude test capsule to simulate the static pressure at altitude by vacuum pumping. Data such as propellant mass flow, oxidizer-to-fuel mixture, and thrust were measured. These measurements were then used to calculate performance factors such as the thrust coefficient, the characteristic exhaust velocity efficiency, and the vacuum specific impulse. In addition, the nozzle temperature was measured to calculate the amount of heat transferred from the combustion gases to the nozzle.

Source record↗

Interface design of VSOP'94 computer code for safety analysis

Today, most software applications, also in the nuclear field, come with a graphical user interface. VSOP'94 (Very Superior Old Program), was designed to simplify the process of performing reactor simulation. VSOP is a integrated code system to simulate the life history of a nuclear reactor that is devoted in education and research. One advantage of VSOP program is its ability to calculate the neutron spectrum estimation, fuel cycle, 2-D diffusion, resonance integral, estimation of reactors fuel costs, and integrated thermal hydraulics. VSOP also can be used to comparative studies and simulation of reactor safety. However, existing VSOP is a conventional program, which was developed using Fortran 65 and have several problems in using it, for example, it is only operated on Dec Alpha mainframe platforms and provide text-based output, difficult to use, especially in data preparation and interpretation of results. We develop a GUI-VSOP, which is an interface program to facilitate the preparation of data, run the VSOP code and read the results in a more user friendly way and useable on the Personal 'Computer (PC). Modifications include the development of interfaces on preprocessing, processing and postprocessing. GUI-based interface for preprocessing aims to provide a convenience way in preparing data. Processing interface is intended to provide convenience in configuring input files and libraries and do compiling VSOP code. Postprocessing interface designed to visualized the VSOP output in table and graphic forms. GUI-VSOP expected to be useful to simplify and speed up the process and analysis of safety aspects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Preliminary Study of Iodine Gas Removal in Sodium Pools

Potential iodine gas release from failed fuel pins is a critical factor in the source term analysis of oxide fuel-loaded sodium fast reactors (SFRs). The accumulated iodine-containing gas mixtures inside pin plenums are expected to be ejected during pin failures and rise through sodium pool, with potential release of gaseous iodine to the cover gas region. Due to its potential radiological impacts, a proper assessment of iodine behavior is necessary for an accurate source term assessment. Throughout the bubble rise trajectory in the sodium pool, iodine gas is continuously removed or transformed at the bubble interface by diffusion, as the combining reaction between the iodine and sodium to form sodium iodide (NaI) is a chemically preferred process. As the final amount of iodine released from the facility is strongly influenced by the removal phenomenon inside the sodium pool, experiments were previously performed by PNC (Power Reactor and Nuclear Fuel Development Corporation) to provide insight into this phenomenon. To assess the accuracy of present methods for predicting iodine gas removal within sodium pools, several candidate approaches, available in source term analysis codes, have been summarized and evaluated in this study. Spherical cap bubbles and spherical bubbles are considered in accordance with the methods adopted in each approach, and different forms of correlations for major parameters have been implemented in accordance with the original adoptions. Based on the summarized results, important aspects to be considered have been derived.

Decontamination↗

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian↗

FDNS CFD Code Benchmark for RBCC Ejector Mode Operation: Continuing Toward Dual Rocket Effects

Computational Fluid Dynamics (CFD) analysis results are compared with benchmark quality test data from the Propulsion Engineering Research Center's (PERC) Rocket Based Combined Cycle (RBCC) experiments to verify fluid dynamic code and application procedures. RBCC engine flowpath development will rely on CFD applications to capture the multi -dimensional fluid dynamic interactions and to quantify their effect on the RBCC system performance. Therefore, the accuracy of these CFD codes must be determined through detailed comparisons with test data. The PERC experiments build upon the well-known 1968 rocket-ejector experiments of Odegaard and Stroup by employing advanced optical and laser based diagnostics to evaluate mixing and secondary combustion. The Finite Difference Navier Stokes (FDNS) code [2] was used to model the fluid dynamics of the PERC RBCC ejector mode configuration. Analyses were performed for the Diffusion and Afterburning (DAB) test conditions at the 200-psia thruster operation point, Results with and without downstream fuel injection are presented.

West, Jeff↗

Development and Validation of a Conjugate Heat Transfer Model for the Two-Phase CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance CFD code NEK5000. The goal of the project is to develop a new two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) [2-3] previously developed at Argonne, but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2019. This report focuses on the FY2019 development and validation of the Conjugate Heat Transfer (CHT) model for the NEK-2P Two-Phase, Computational Fluid Dynamics (CFD) code. The NEK-2P CHT model calculates the coupled behavior of the fluid and solid domains. In the fluid domain NEK-2P calculates the vapor and liquid phase temperatures, velocities and mass fractions. In the solid domain only the temperatures are calculated. The fluid and solid domains are coupled with a heat flux boundary condition at the fluid-solid interface. The heat flux and the wall temperatures at the fluid-solid interface are updated using the NEK-2P wall heat-flux partitioning model. A wall heat flux boundary condition is also applied to the exterior surface of the solid domain. The CHT model validation was illustrated through the analysis of three of the Becker Critical Heat Flux (CHF) tests. Reasonably good agreement with measured data was obtained in predicting the CHF location and post CHF wall temperature behavior illustrating the ability of the NEK-2P CHT model to simulate the CHF phenomena for a wide range of thermal-hydraulic conditions. Apart from the CHT model development and validation, initial simulations were performed for the Virginia Tech. (VT) 3x3 bare rod bundle experiments including both cold and boiling flow simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗