Search NASA⌕ Search

SEARCH · Search NASA

Results for “CORE FLOW”

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 379 records · Page 21

Nuclear Thermal Rocket Emulator for a Hardware-in-the-Loop Test Bed

To support NASA’s mission to use nuclear thermal rockets for future Mars missions, an instrumentation and control test bed has been built at Oak Ridge National Laboratory. The system is designed as a hardware-in-the-loop test bed for testing control elements and autonomous control algorithms for nuclear thermal propulsion rockets. The mock reactor system consists of a modular and scalable framework, using inexpensive components and open-source software. The hardware system consists of a two-phase flow loop and a mock reactor with six control drums. A single-board computer (NVIDIA Jetson) handles reactor core emulation and hosts a message queuing telemetry transport broker that allows user-deployed control algorithms to interact with the system hardware. The reactor emulator receives sensor data from the hardware and provides the simulated performance of the reactor under steady-state, transient, and fault conditions. The emulator uses a reactivity lookup table and the point kinetics equations to solve for the reactor dynamics in real time. Emulated reactor dynamics and sensor input inform the autonomous control algorithm’s decision-making in a closed-loop manner. The current system is capable of operating at 10 Hz, but faster cycle rates are an area of ongoing research. This test bed will enable NASA and other space vendors to rigorously test their autonomous control systems for NTP rockets under transient (reactor startup and shutdown), steady-state, and fault conditions to reduce development time and risk for autonomous control systems in future missions.

autonomous control↗

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗

Nanoscale Hybrid Electrolytes with Viscosity Controlled Using Ionic Stimulus for Electrochemical Energy Conversion and Storage

As renewable energy is rapidly integrated into the grid, the challenge has become storing intermittent renewable electricity. Technologies including flow batteries and CO 2 conversion to dense energy carriers are promising storage options for renewable electricity. To achieve this technological advancement, the development of next generation electrolyte materials that can increase the energy density of flow batteries and combine CO 2 capture and conversion is desired. Liquidlike nanoparticle organic hybrid materials (NOHMs) composed of an inorganic core with a tethered polymeric canopy (e.g., polyetheramine (HPE)) have a capability to bind chemical species of interest including CO 2 and redox-active species. In this study, the unique response of NOHM-I-HPE-based electrolytes to salt addition was investigated, including the effects on solution viscosity and structural configurations of the polymeric canopy, impacting transport behaviors. The addition of 0.1 M NaCl drastically lowered the viscosity of NOHM-based electrolytes by up to 90%, reduced the hydrodynamic diameter of NOHM-I-HPE, and increased its self-diffusion coefficient, while the ionic strength did not alter the behaviors of untethered HPE. This study is the first to fundamentally discern the changes in polymer configurations of NOHMs induced by salt addition and provides a comprehensive understanding of the effect of ionic stimulus on their bulk transport properties and local dynamics. These insights could be ultimately employed to tailor transport properties for a range of electrochemical applications.

25 ENERGY STORAGE↗

Design of an Additively Manufactured Recuperator With 800 °C Inlet Temperature for sCO2 Power Cycle Application

Abstract Recuperators with design temperatures at and above 800 °C can further increase the thermal efficiency of supercritical CO2 power cycles by enabling higher turbine exhaust temperatures. Mar-M247 is a well-suited nickel-based superalloy for high temperature service due to its high creep strength that prevents excessive material thickness being required for pressure containment. Additive manufacturing using a high-speed laser-directed energy deposition (L-DED) process presents a promising solution, with build trials demonstrating the ability to produce nonconventional flow channels for enhanced heat transfer. A design process is presented that includes aerothermal and mechanical evaluation to maximize performance within the constraints of the manufacturing process. A 2-D heat transfer network and pressure drop code allows prediction of flow distribution and its effect on overall thermal performance. Established literature correlations, along with CFD simulation, inform the prediction of heat transfer coefficients and friction factors for the flowpaths and enhancement features in the heat exchanger core. Mechanical evaluation using finite element analysis (FEA) modeling with the intent of the ASME Boiler & Pressure Vessel Code (BPVC) Section VIII, Div. 2 assesses the operational safety of the design. The detailed design features annular finned passages that take advantage of helical flow paths to distribute the flow from separated headers to shared heat transfer surfaces. Performance predictions for the recuperator at a 50 kW scale provide insights into the feasibility of the additively manufactured (AM) process to produce recuperators on a commercial scale that extend existing operating envelopes.

Engineering↗

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗

A review on the state of thermal hydraulics research on air ingress scenarios in High-Temperature Gas-cooled Reactors following a D-LOFC

With the expectation of near-immediate carbon neutrality, widespread implementation of proven High-Temperature Gas-cooled Reactors (HTGRs) embodies a viable solution pathway given their inherent, passive safety features and high thermal efficiency. This study provides an overview of the current state of research involving the thermal hydraulics associated with air ingress from a depressurized loss of forced cooling (D-LOFC) in HTGRs. Accurately characterizing and predicting the physical phenomena underlying air ingress is of paramount concern, as the integrity of the fuel and core graphite support structures are threatened by the presence of oxygen. Broadly speaking, the air ingress scenario can be delineated into three main stages: (1) Depressurization, (2) Density-Driven Flow, and (3) Natural Convection. In tandem with the underlying fundamental theory, this review collates and synthesizes the existing body of contemporary research concerning the air ingress scenario following a D-LOFC. As evinced by this review, our current understanding and predictive abilities have benefited from extensive research, predominantly concentrated on the rate of air ingestion into the core. Here, additional research is necessary to holistically capture the phenomenology of an air ingress scenario following a D-LOFC by considering an additional variable: the oxygen content of the ingressing air. The latter variable requires investigation into the complex interactions of the fully integrated system. Additionally, while numerical tools are evolving domestically through the Nuclear Energy Advanced Modeling and Simulation program, a sufficiently validated code remains absent.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the effects of neutron irradiation on first-generation corrosion mitigation coatings on SiC for accident-tolerant fuel cladding

In this work, high purity SiC and SiC/SiC composites coated with commercial TiN, Cr, CrN, or CrN/Cr multilayer coatings were irradiated in Ar or flowing PWR water in the Massachusetts Institute of Technology Nuclear Reactor Laboratory (MITR). Irradiation in Ar was performed in the core. In the water environment, identical samples were placed in one of three different locations: in-core, providing exposure to neutron damage and radiolysis-affected water; above-core, where samples were exposed to radiolysis-affected water but not neutron damage, or outside of the core, where samples were exposed to the coolant water without the effects of radiation. Radiation in Ar revealed significant cracking of all but the TiN coatings, attributed to differential swelling between the coating and substrate. Lattice swelling was not observed in any of the coatings, but 0.2% void swelling was observed in the Cr coating. All of the coatings failed during water exposures in the core. Likewise, CrN/Cr spalled in each condition. Cr was protective, except under radiation damage as a result of cracking, and TiN severely degraded in the core with no coating was found following exposure. A SiC/coating ATF cladding system is anticipated to perform adequately following improvements in coating ductility and purity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Observation of Double Excitations in the Resonant Inelastic X-ray Scattering of Nitric Oxide

The nitrogen K-edge resonant inelastic X-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft X-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s → π* and 1s → Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s → π* and π → π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis

Persister cells are genetically identical variants in a bacterial population that have phenotypically modified their physiology to survive environmental stress. In bacterial pathogens, persisters are able to survive antibiotic treatment and reinfect patients in a frustrating cycle of chronic infection. To better define core persistence mechanisms for therapeutics development, we performed transcriptomics analyses of Burkholderia thailandensis populations enriched for persisters via three methods: flow sorting for low proton motive force, meropenem treatment, and culture aging. Although the three persister-enriched populations generally displayed divergent gene expression profiles that reflect the multimechanistic nature of stress adaptations, there were several common gene pathways activated in two or all three populations. These include polyketide and nonribosomal peptide synthesis, Clp proteases, mobile elements, enzymes involved in lipid metabolism, and ATP binding cassette (ABC) transporter systems. In particular, identification of genes that encode polyketide synthases (PKSs) and fatty acid catabolism factors indicates that generation of secondary metabolites, natural products, and complex lipids could be part of the metabolic program that governs the persistence state. We also found that loss-of-function mutations in the PKS-encoding gene locus BTH_I2366, which plays a role in biosynthesis of histone deacetylase (HDAC) inhibitors, resulted in increased sensitivity to antibiotics targeting DNA replication. Furthermore, treatment of multiple bacterial pathogens with a fatty acid synthesis inhibitor, CP-640186, potentiated the efficacy of meropenem against the persister populations. Altogether, our results suggest that bacterial persisters may exhibit an outwardly dormant physiology but maintain active metabolic processes that are required to maintain persistence.

59 BASIC BIOLOGICAL SCIENCES↗

GaN Core-shell Nanofin Vertical Transistor (CoNVerT): A New Direction for Power Electronics (Final Scientific/Technical Report)

A novel power transistor architecture, the GaN c ore-shell n anofin ver tical transistor (CoNVerT) to address fundamental challenges in realizing the ultimate limit of GaN power transistor performance was explored experimentally. This technology promises ultra-high-efficiency high voltage/high power applications (e.g. DC/DC converters, motor control, fast charging, actuation), as well as to operate in harsh environments. The device exploits a vertical superjunction structure based on an experimentally-validated core-shell nanofin growth process in which lateral p-n heterojunctions are formed in a single growth step, while still maintaining vertical current flow for compact die size and low cost. The concept leverages the best properties of GaN for mid-range voltage applications: high mobility, high breakdown voltage, and native heterojunction enhancement-mode operation. Due to the crystallographic nature of the nanofin growth by molecular beam epitaxy, the sidewall heterojunctions occur on non-polar planes, resulting in ultra-smooth interfaces for high mobility, no sidewall etch damage and related surface/interface states, and elimination of piezoelectric effects that can limit reliability in conventional structures. This also facilitates superjunction formation for maximum device performance, and the selective-area growth of the nanofin results in dislocation-free growth, even on low-cost Si (111) substrate. In this program, core-shell nanofins were grown by molecular beam epitaxy, test structures to evaluate the doping, resistivity, and other electrical properties were fabricated, and the material and test structures were characterized in detail. The work identified clear potential (e.g., the doping was well controlled as required for superjunction concepts), but also additional areas that require additional effort to resolve (some unexpected crystal defects were encountered that require additional engineering to overcome). Simulation studies of the proposed concept validate that the fundamental approach is very promising, but additional effort in experimental realization is needed.

42 ENGINEERING↗

A multi-scale experimental investigation of flow properties in coarse-grained hydrate reservoirs during production (Final Scientific/Technical Report)

In this project, “A multi-scale experimental investigation of flow properties in coarse-grained hydrate reservoirs during production” we explored the relative permeability and dissipation behavior of hydrate-bearing coarse-grained sediment at the core scale and we observed the formation and dissociation of these materials at the pore scale with micro-CT and Raman experiments. The 3-year project results inform reservoir simulation efforts, critical to determining the viability of the coarse-grained hydrate reservoir as an energy resource. We performed the first 3 phase relative permeability experiments in a hydrate-bearing medium using the steady state method (gas and water flow in the presence of hydrate). Results support a model where water is the most wetting phase, hydrate is the intermediate wetting phase, and gas is the least wetting phase. This model can be used to predict relative permeability of gas and water in the presence of hydrate. We showed that when hydrate dissociates, it does so at a pressure and temperature predicted for fresh water (no salinity) conditions. At the scale of a grid block, reservoir simulation models should assume that hydrate dissociates at the freshwater phase boundary. We observed that the chemistry and the pore habit of methane hydrates change dramatically over the timescales of experiments (hours to weeks). It takes weeks to months for hydrates to evolve from a non-stoichiometric to stoichiometric compound with 3:1 large cage vs small cage occupancy. We documented hydrates forming initially in small clayey silt surfaces but gradually concentrated in large pores in sand-sized sediment.

03 NATURAL GAS↗

TChem v2.0 - A Software Toolkit for the Analysis of Complex Kinetic Models

TChem is an open source software library for solving complex computational chemistry problems and analyzing detailed chemical kinetic models. The software provides support for: complex kinetic models for gas-phase and surface chemistry; thermodynamic properties based on NASA polynomials; species production/consumption rates; stable time integrator for solving stiff time ordinary differential equations; and, reactor models such as homogenous gas-phase ignition (with analytical Jacobian matrices), continuously stirred tank reactor, plug-flow reactor. This toolkit builds upon earlier versions that were written in C and featured tools for gas-phase chemistry only. The current version of the software was completely refactored in C++, uses an object-oriented programming model, and adopts Kokkos as its portability layer to make it ready for the next generation computing architectures i.e., multi/many core computing platforms with GPU accelerators. We have expanded the range of kinetic models to include surface chemistry and have added examples pertaining to Continuously Stirred Tank Reactors (CSTR) and Plug Flow Reactor (PFR) models to complement the homogenous ignition examples present in the earlier versions. To exploit the massive parallelism available from modern computing platforms, the current software interface is designed to evaluate samples in parallel, which enables large scale parametric studies, e.g. for sensitivity analysis and model calibration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Update on Cardinal: Toward Full Core Pebble Simulations

The new DOE-NE Center of Excellence for Thermal-fluids applications in Nuclear Energy inaugurated in April 2018 considers and researches novel new solution strategies for historically challenging flow issues that still plague the current fleet of deployed Light Water nuclear Reactors (LWRs) as well as predicting various fluid flow and fluid related issues with advanced reactor technologies. Our advanced thermal-fluids research and development approach synergistically combines three natural, though overlapping, length and time scales in a hierarchal multi-scale approach to avoid the temptation and pitfalls of attempting to develop a single solve-all algorithm for physical fluid flow problems that will span nine orders of magnitude in spatial and temporal scales. To demonstrate the multi-scale philosophy of the center we focus on Fluoride Cooled High Temperature Reactors (FHRs), and in particular on the Berkley’s PB-FHR Mark-I design. The Fluoride salt cooled High temperature Reactor (FHR) is a class of advanced nuclear reactors that combine the robust coated particle fuel form from high temperature gas cooled reactors, direct reactor auxiliary cooling system (DRACS) passive decay removal of liquid metal fast reactors, and the transparent, high volumetric heat capacitance liquid Fluoride salt working fluids - such as Flibe - from molten salt reactors. This combination of fuel and coolant enables FHRs to operate in a high-temperature low-pressure design space that has beneficial safety and economic implications. The PB-FHR reactor relies on a pebble bed approach and pebble bed reactors are in a sense the poster child for multi-scale analysis. The lower length-scale simulator for pebble reactor cores comprises three physics: neutronics (OpenMC), thermal-fluids (Nek5000/NekRS) and fuel performance (BISON). As part of Center of Excellence ongoing research efforts we have developed Cardinal, a new tool platform for lower length-scale simulation. Cardinal tightly couples all three physics and leverages advances in MOOSE such as the MultiApp system and the concept of MOOSE-wrapped Apps. The present report provides an update on the development of Cardinal with an extension of Cardinal to GPUs. We also perform a first-of-a-kind demonstration simulation on Summit representing a 10x capability increase in terms of pebble count for Cardinal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗