Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tabulation”

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 37 records · Page 2

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗

Correct Interpretations of ENDF-102 Definitions for Resonance Effects

My Uncle Willie circa 1600 wrote “What’s in a name; a rose by any other name would smell as sweet.” I fear in this case we have a somewhat similar problem in that we may be using the same word but are not using the same definition; specifically, the word Unresolved. The simplest physics definition as it applies to neutron resonances, is the energy point where we can no longer see/measure ALL – let me repeat that – ALL - of the individual resonances. That seems simple and clear, but the question is: how to represent resonances beyond this point in order to accurately reproduce the effects we have seen in measurements and expect/need to reproduce in our applications. We know there are more, unseen resonances, otherwise we wouldn’t say Unresolved. The ENDF approach is well defined in ENDF-102 and simple: for ENDF data the only way to represent Unresolved data is by using a theoretical model to define the distribution of resonances, including those that are too narrow to measure (i.e., are unresolved). It is important to note that in ENDF this is the one and only Unresolved model, e.g., there is no provision in ENDF to accurately define individually ALL resonances above the Resolved energy range – by ALL here I mean both those that we can measure and those that we cannot individually measure, but that theory and integral measurements tells us are present. An alternative approach, which would appear to be equally valid, would be to include the latest measured data as tabulated energy expendent data extending upwards in energy above the Resolved energy range. In this approach the evaluation would not include an ENDF style Unresolved energy range; it would only include a Resolved resonance region, followed by tabulated higher energy points, representing the resonances that could be measured beyond the Resolved range. But an important point to note: By listing these resonances above the resolved energy one admits that at least some resonances in this energy range are missing as Unresolved; i.e., they are too narrow or overlapping to measure. The purpose of this paper is to illustrate that the later approach, while done with good intentions, and appearing to be valid/adequate in plots, does not meet the need of our engineering applications. Why? As we will see below, of these two possible approaches, only the ENDF use of a model to statistically include the missing, i.e., unresolved, resonances, can meet our engineering needs to reproduce the integral effects we have measured and understand. Only with this statistical model can we predict and include in our calculated results the important effects of temperature (Doppler broadening), and energy integrals (self-shielding). Below I will first present results using two ENDF/B-VIII.1 evaluations, U235 and U238, that use the correct ENDF-102 definition of an Unresolved resonance region, using a statistical model to include the effects of resonances that theory predicts are present, but are too narrow to measure. These two evaluations reproduce the expected temperature (Doppler) and energy integral (self-shielding) effects that we expect. Next I will present results using one ENDF/B-VIII.1 evaluation, 26-Fe-56, that does not use an ENDF-102 Unresolved resonance region; instead above its Resolved energy range it lists many tabulated energy points, that look like measured data, but by definition, since they are included above the ENDF Resolved energy range there are missing Unresolved resonances, i.e., there are missing the resonances that are too narrow to resolve, i.e., are unresolved. My conclusion, and I hope yours, is that the below figures illustrate that this approach does not reproduce the temperature and energy integrals that we expect and need to accurately calculate results for our fission reactor calculations. As such this approach should not be used in ENDF formatted evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Robust second-order approximation of the compressible Euler equations with an arbitrary equation of state

Here, this paper is concerned with the approximation of the compressible Euler equations supplemented with an arbitrary or tabulated equation of state. The proposed approximation technique is robust, formally second-order accurate in space, invariant-domain preserving, and works for every equation of state, tabulated or analytic, provided the pressure is nonnegative. An entropy surrogate functional that grows across shocks is proposed. The numerical method is verified with novel analytical solutions and then validated with several computational benchmarks seen in the literature including problems with composite waves.

97 MATHEMATICS AND COMPUTING↗

Electron Bifurcating Hydrogenases

The importance of electron-bifurcating enzymes is manifest by their ability to maximize energy efficiency. Specifically, they couple a downhill oxidation-reduction (redox) reaction with an uphill redox reaction. Since the rapid increase in the discovery of bifurcating enzymes starting in 2008, there has been interest in incorporating their mechanistic principles into artificial/semiartificial systems to drive chemically challenging reactions. This has yet to be achieved, partly because the details of electron bifurcation, i.e. mechanisms, are largely elusive. Nevertheless, much progress has been made in understanding reactivities, structures, and some mechanistic aspects of these enzymes. Notable examples are electron-bifurcating hydrogenases, which are the focus of this chapter. The chapter is organized as follows. Section 11.1 provides an overview of hydrogenases and electron bifurcation. In Section 11.2, some physiological roles of electron-bifurcating hydrogenases are highlighted. Additionally, electron-bifurcating subunit compositions and biochemical reactivities are comprehensively tabulated, and some key points/considerations about these are noted. In Section 11.3, we discuss the known structures of these enzymes, which provide insight into their complex arrangements of redox cofactors, such as iron-sulfur clusters. Also provided are tabulations and discussions of some biophysical properties of the cofactors. In Section 11.4, we discuss the mechanistic proposals of these enzymes, which are primarily based on structural information. Areas of research that are much needed are outlined in Section 11.5. We conclude on the note that what is learned from electron-bifurcating hydrogenases has applicability to other bifurcating enzymes, nonbifurcating analogs, and mechanistic enzymology at large.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Free Energies of Proton-Coupled Electron Transfer Reagents and Their Applications

Here, we present an update and revision to our 2010 review on the topic of proton- coupled electron transfer (PCET) reagent thermochemistry. Over the past decade, the data and thermochemical formalisms presented in that review have been of value to multiple fields. Concurrently, there have been advances in the thermochemical cycles and experimental methods used to measure these values. This Review (i) summarizes those advancements, (ii) corrects systematic errors in our prior review that shifted many of the absolute values in the tabulated data, (iii) provides updated tables of thermochemical values, and (iv) discusses new conclusions and opportunities from the assembled data and associated techniques. We advocate for updated thermochemical cycles that provide greater clarity and reduce experimental barriers to the calculation and measurement of Gibbs free energies for the conversion of X to XHn in PCET reactions. In particular, we demonstrate the utility and generality of reporting potentials of hydrogenation, E°(V vs H 2 ), in almost any solvent and how these values are connected to more widely reported bond dissociation free energies (BDFEs). The tabulated data demonstrate that E°(V vs H 2 ) and BDFEs are generally insensitive to the nature of the solvent and, in some cases, even to the phase (gas versus solution). This Review also presents introductions to several emerging fields in PCET thermochemistry to give readers windows into the diversity of research being performed. Some of the next frontiers in this rapidly growing field are coordination-induced bond weakening, PCET in novel solvent environments, and reactions at material interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimized spatial information for 1990, 2000, and 2010 U.S. census microdata

Abstract We report on the successful completion of a project to upgrade the positional accuracy of every response to the 1990, 2000, and 2010 U.S. decennial censuses. The resulting data set, called Optimized Spatial Census Information Linked Across Time (OSCILAT), resides within the restricted-access data warehouse of the Federal Statistical Research Data Center (FSRDC) system where it is available for use with approval from the U.S. Census Bureau. OSCILAT greatly improves the accuracy and completeness of spatial information for older censuses conducted prior to major quality improvements undertaken by the Bureau. Our work enables more precise spatial and longitudinal analysis of census data and supports exact tabulations of census responses for arbitrary spatial units, including tabulating responses from 1990, 2000, and 2010 within 2020 block boundaries for precise measures of change over time for small geographic areas.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

A compendium and meta-analysis of flatband potentials for TiO2, ZnO, and SnO2 semiconductors in aqueous media

Semiconductor/electrolyte interfaces are of great interest to numerous scientific fields including renewable energy, (photo)electrochemistry, and energy storage. The semiconductor flatband potential is a key parameter in locating the conduction band minimum or valence band maximum of the semiconductor material in electrolyte. Despite its importance for quantifying the energetic location of the semiconductor bands, literature reports for the same material demonstrate significant variability in the flatband potential. In this compendium and meta-analysis, reported flatband potentials of the common semiconductor materials TiO2, SnO2, and ZnO in aqueous electrolyte were compiled and assessed to quantify the spread in literature flatband potentials as well as determine the factors that lead to the significant spread. For TiO2, SnO2, and ZnO, literature flatband potentials referenced to the reversible hydrogen electrode span a range of nearly 2 V each. Flatband potential tabulations were separated by variables such as the solution pH, the crystalline polymorph, the crystal facet, the morphology, and the dimensions or combinations of these variables to assess the factors that contribute to the observed spread. Important and surprising findings from these categorizations are summarized: (1) Even for the narrowest categorizations, the spread in flatband potential is still large. (2) Flatband potentials of TiO2 and SnO2 follow the expected Nernstian dependence with solution pH. ZnO materials deviate from this Nernstian dependence. (3) In the aggregate, there is no statistically significant difference in the reported flatband potentials of anatase and rutile TiO2. Single crystal tabulations were the only distributions to have statistically significant differences in the flatband potential between anatase and rutile TiO2. (4) Anatase TiO2 materials with a nanotube morphology appear to have a +400 mV difference in mean flatband potential compared to all other morphologies, but we argue that this is likely due to widespread misuse of the Mott–Schottky analysis. Other interesting findings are revealed within the spread of literature flatband potentials, and possible explanations are provided to generate discussion. We also briefly review and discuss common techniques that were used to determine the flatband potential and the pitfalls/criticisms of these techniques. Last, we discuss some ways in which future research on the determination of the flatband potential can be performed to improve the reliability of reported values and the quality of the work. In total, the results from this meta-analysis suggest multiple factors can affect the measured flatband potential and that an abundance of caution should be applied when attempting to quantify the flatband potential of complex or nanostructured systems.

Patel, Milan Y. (ORCID:0000000291948141)↗

Large-eddy simulation of a multi-injection flame in a diesel engine environment using an unsteady flamelet/progress variable approach

In this work, large-eddy simulations (LESs) are conducted for a multiple-injection flame in a diesel engine environment using an unsteady flamelet/progress variable (UFPV) approach in which differential diffusion is considered. The suitability of the UFPV tabulation approach is first evaluated through a priori analyses using the state-of-the-art direct numerical simulation (DNS) dataset. Both the instantaneous data and the conditional values for the major and minor species' mass fractions are compared between the UFPV and the DNS. The comparisons show that the proposed UFPV tabulation approach can give good predictions for the multiple-injection flame at different injection phases. While the gas temperature and major species mass fractions can be accurately predicted with or without differential diffusion being considered in the UFPV flamelet library, the prediction accuracy for the highly diffusive species (e.g., hydrogen) in the main injection phase can be noticeably improved when differential diffusion is taken into account. The fully coupled LES/UFPV simulations show that the overall structure of the multiple-injection flame can be predicted, and the conditional thermo-chemical values are close to the filtered DNS dataset. The reasons for the remaining discrepancies found in the a priori analyses and the a posteriori simulations using the UFPV approach are analyzed.

Wen, Xu (ORCID:0000000331334292)↗

2D k -th nearest neighbour statistics: a highly informative probe of galaxy clustering

ABSTRACT Beyond standard summary statistics are necessary to summarize the rich information on non-linear scales in the era of precision galaxy clustering measurements. For the first time, we introduce the 2D k-th nearest neighbour (kNN) statistics as a summary statistic for discrete galaxy fields. This is a direct generalization of the standard 1D kNN by disentangling the projected galaxy distribution from the redshift-space distortion signature along the line-of-sight. We further introduce two different flavours of 2D kNNs that trace different aspects of the galaxy field: the standard flavour which tabulates the distances between galaxies and random query points, and a ‘DD’ flavour that tabulates the distances between galaxies and galaxies. We showcase the 2D kNNs’ strong constraining power both through theoretical arguments and by testing on realistic galaxy mocks. Theoretically, we show that 2D kNNs are computationally efficient and directly generate other statistics such as the popular two-point correlation function (2PCF), voids probability function, and counts-in-cell statistics. In a more practical test, we apply the 2D kNN statistics to simulated galaxy mocks that fold in a large range of observational realism and recover parameters of the underlying extended halo occupation distribution (HOD) model that includes velocity bias and galaxy assembly bias. We find unbiased and significantly tighter constraints on all aspects of the HOD model with the 2D kNNs, both compared to the standard 1D kNN, and the classical redshift-space 2PCF.

79 ASTRONOMY AND ASTROPHYSICS↗

From sensing to acclimation: The role of membrane lipid remodeling in plant responses to low temperatures

Abstract Low temperatures pose a dramatic challenge to plant viability. Chilling and freezing disrupt cellular processes, forcing metabolic adaptations reflected in alterations to membrane compositions. Understanding the mechanisms of plant cold tolerance is increasingly important due to anticipated increases in the frequency, severity, and duration of cold events. This review synthesizes current knowledge on the adaptive changes of membrane glycerolipids, sphingolipids, and phytosterols in response to cold stress. We delve into key mechanisms of low-temperature membrane remodeling, including acyl editing and headgroup exchange, lipase activity, and phytosterol abundance changes, focusing on their impact at the subcellular level. Furthermore, we tabulate and analyze current gycerolipidomic data from cold treatments of Arabidopsis, maize, and sorghum. This analysis highlights congruencies of lipid abundance changes in response to varying degrees of cold stress. Ultimately, this review aids in rationalizing observed lipid fluctuations and pinpoints key gaps in our current capacity to fully understand how plants orchestrate these membrane responses to cold stress.

Plant Sciences↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗

Reduce-Order Modeling of Multigroup Neutron Cross Sections for High-Temperature Gas-cooled Reactors

Deterministic neutronics calculations rely on multigroup neutron cross section libraries, which usually consists of a database of tabulated values, used to calculate the cross sections through multivariate linear interpolation. However, interpolation of the multidimensional cross section data becomes memory inefficient and time consuming as the number of tabulations increases, significantly slowing down the neutronics calculation, especially in the case of micro cross section libraries where every isotope (on the order of hundreds) has its own set of specific reactions and cross sections. To address this challenge, this work constructs efficient and robust reduced-order models (ROMs) of the multi-group cross sections to support the Griffin simulation of high-temperature gas-cooled reactors (HTGRs). The first part of the study investigates the linearity of the multi-group cross section data across isotopes, reaction types and energy groups on pre-generated datasets for the purpose of dimensionality reduction. Secondly, a down-selection of ROM techniques is presented on representative classical machine learning (ML) techniques, including variants of linear regression, kernel-based methods, tree-based algorithms, and artificial neural networks. The selection criteria jointly consider the memory efficiency, predictive accuracy, prediction speed, and scalability in comparison to the multidimensional interpolation. Among all the ML techniques, deep neural networks (DNNs) have proven to be the best selection with sufficient accuracy, high robustness, good memory efficiency, great scalability, and superior flexibility. DNNs for have been trained for all isotopes in this work and systematic Griffin testing is ongoing at this moment to ensure the feasibility of this ROM technique for cross section predictions.

42 - ENGINEERING↗

Reduced-Order Modeling of Multigroup Neutron Cross Sections for High-Temperature Gas-cooled Reactors

Abstract – Deterministic neutronics calculations rely on multigroup neutron cross section libraries, which consist of databases of tabulated values, used to calculate the neutron cross sections through multivariate linear interpolation. However, interpolation of the multidimensional cross section data becomes memory inefficient and time consuming as the number of tabulations increases, significantly slowing down the neutronics calculation, especially in the case of microscopic cross section libraries where every isotope (on the order of hundreds) has its own set of specific reactions and cross sections. In order to address this challenge, this work constructs efficient and robust reduced-order models (ROMs) of the multi-group cross sections to support the Griffin simulation of high-temperature gas-cooled reactors (HTGRs). The first part of the study investigates the linearity of the multigroup cross section data across isotopes, reaction types, and energy groups on pre-generated datasets for the purpose of dimensionality reduction. Secondly, a down-selection of ROM techniques is presented on representative classical machine learning (ML) techniques, including variants of linear regression, kernel-based methods, tree-based algorithms, and artificial neural networks. The selection criteria jointly consider the memory efficiency, predictive accuracy, prediction speed, and scalability in comparison to the multidimensional interpolation. Among all the ML techniques, deep neural networks (DNNs) have proven to be the best selection with sufficient accuracy, high robustness, good memory efficiency, great scalability, and superior flexibility. DNNs have been trained for all isotopes in this work and systematic Griffin testing is ongoing to ensure the feasibility of this ROM technique for predicting cross section and reducing memory requirements without a significant sacrifice in computational performance.

42 - ENGINEERING↗