Search NASA⌕ Search

SEARCH · Search NASA

Results for “code comparison study”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Quantitative measurement of transverse injector and free stream interaction in a nonreacting SCRAMJET combustor using laser-induced iodine fluorescence

A preliminary quantitative study of the compressible flowfield in a steady, nonreacting model SCRAMJET combustor using laser-induced iodine fluorescence (LIIF) is reported. Measurements of density, temperature, and velocity were conducted with the calibrated, nonintrusive, optical technique for two different combustor operating conditions. First, measurements were made in the supersonic flow over a rearward-facing step without transverse injection for comparison with calculated pressure profiles. The second configuration was staged injection behind the rearward-facing step at an injection dynamic pressure ratio of 1.06. These experimental results will be used to validate computational fluid dynamic (CFD) codes being developed to model supersonic combustor flowfields.

Fletcher, D. G.↗
NASA NTRS1993

Alternate vehicles for engine/vehicle optimization

The space shuttle main engine, full flow staged combustion, integrated modular engine concept, split expander, and expander bleed rocket engine cycles were studied for use on Single-Stage-to-Orbit fully reusable vehicles. The vehicle uses vertical take-off and horizontal landing and liquid hydrogen rockets for propulsion. The PROPSIZE computer code was modified for use on the available computers. Vehicle sizing was done after trajectory optimization for various engine designs. The various rocket engine cycles were compared and it was found that the full flow staged combustion cycle had the lowest vehicle dry mass. The split expander cycle could have advantages in safety and cooling which could make it the preferred choice in some cases. Comparison to previous work with two-stage heavy-lift vehicles showed that the same ranking of the cycles was found for both vehicle concepts.

Venkatasubramanyam, G.↗
NASA NTRS2009

Investigating Response from Turbulent Boundary Layer Excitations on a Real Launch Vehicle using SEA

Statistical Energy Analysis (SEA) response has been fairly well anchored to test observations for Diffuse Acoustic Field (DAF) loading by others. Meanwhile, not many examples can be found in the literature anchoring the SEA vehicle panel response results to Turbulent Boundary Layer (TBL) fluctuating pressure excitations. This deficiency is especially true for supersonic trajectories such as those required by this nation s launch vehicles. Space Shuttle response and excitation data recorded from vehicle flight measurements during the development flights were used in a trial to assess the capability of the SEA tool to predict similar responses. Various known/measured inputs were used. These were supplemented with a range of assumed values in order to cover unknown parameters of the flight. This comparison is presented as "Part A" of the study. A secondary, but perhaps more important, objective is to provide more clarity concerning the accuracy and conservatism that can be expected from response estimates of TBL-excited vehicle models in SEA (Part B). What range of parameters must be included in such an analysis in order to land on the conservative side in response predictions? What is the sensitivity of changes in these input parameters on the results? The TBL fluid structure loading model used for this study is provided by the SEA module of the commercial code VA One.

Harrison, Phillip↗
NASA NTRS2000

Parametric Weight Comparison of Advanced Metallic, Ceramic Tile, and Ceramic Blanket Thermal Protection Systems

A parametric weight assessment of advanced metallic panel, ceramic blanket, and ceramic tile thermal protection systems (TPS) was conducted using an implicit, one-dimensional (I-D) finite element sizing code. This sizing code contained models to account for coatings fasteners, adhesives, and strain isolation pads. Atmospheric entry heating profiles for two vehicles, the Access to Space (ATS) vehicle and a proposed Reusable Launch Vehicle (RLV), were used to ensure that the trends were not unique to a certain trajectory. Ten TPS concepts were compared for a range of applied heat loads and substructural heat capacities to identify general trends. This study found the blanket TPS concepts have the lightest weights over the majority of their applicable ranges, and current technology ceramic tiles and metallic TPS concepts have similar weights. A proposed, state-of-the-art metallic system which uses a higher temperature alloy and efficient multilayer insulation was predicted to be significantly lighter than the ceramic tile stems and approaches blanket TPS weights for higher integrated heat loads.

Myers, David E.↗
DOE OSTI2022

Misalignment of magnetic field in DIII-D assessed by post-mortem analysis of divertor targets

Abstract We assess the toroidal magnetic field B t asymmetry in DIII-D due to a misalignment of the toroidal field coils with respect to the poloidal magnetic field coils and vacuum vessel. The peak-to-peak variation of the divertor strike point (SP) radius is measured to be 1 cm, with an n = 1 toroidal pattern. We use the centre of a narrow carbon deposition band on tungsten-coated divertor tiles just inside the outer strike point (OSP) as a proxy for the divertor SP location. The band occurred in a series of reverse B t discharges with the OSP positioned on the divertor inserts due to strong E × B drift transport of C from the inner to the outer SP through the private flux region. The variation in band radius (and hence the magnetic SP) is a (4.89 ± 0.31) mm shift toward (310 ± 4)° toroidal direction. These measurements agree well with previous measurements of the 3D magnetic field distribution (Luxon 2003 Nucl. Fusion 43 1813), simulations performed by the mafot field line integration code, and recent Langmuir probe measurements in the small-angle-slot (SAS) divertor (Watkins et al 2019 Nucl. Mater. Energy 18 46). Comparison of these measurements in the SAS divertor also indicates that there is the possibility of a tilt (in conjunction with the shift) of the B t coil field of (0.04 ± 0.07)° towards the toroidal angle of (215 ± 25)°. Previous measurements suggested a field misalignment of (4.6 ± 0.3) mm in the 270° toroidal direction, and a tilt of (0.06 ± 0.02)° toward the 114° toroidal direction, which is similar to the results reported here. These studies will be important for better understanding the radial variation of the toroidal strike line in DIII-D, for designing the new generation of SAS divertor, and for developing an understanding of the impact of error fields on tokamaks with tightly baffled slot divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

ChatMPI: LLM-Driven MPI Code Generation for HPC Workloads

The Message Passing Interface (MPI) standard plays a crucial role in enabling scientific applications for parallel computing and is an essential component in high-performance computing (HPC). However, implementing MPI code manually—especially applying a proper domain decomposition and communication pattern—is a challenging and error-prone task. We present ChatMPI, an AI assistant for MPI parallelization of sequential C codes. In our analysis, we focus on testing six essential HPC workloads, which are based on Basic Linear Algebra Subprograms levels 1, 2, and 3 as well as sparse, stencil, and iterative operations. We analyze the process of creating ChatMPI by using the ChatHPC library. This lightweight large language model (LLM)–based infrastructure enables HPC experts to efficiently create and supervise trustworthy AI capabilities for critical HPC software tasks. We study the data required for training (fine-tuning) ChatMPI to generate parallel codes that not only use MPI syntax correctly but also apply HPC techniques to reduce memory communication and maximize performance by using proper work decomposition. With a relatively small training dataset composed of a few dozen prompts and fewer than 15 minutes of fine-tuning on one node equipped with two NVIDIA H100 GPUs, ChatMPI elevates trustworthiness for MPI code generation of current LLMs (e.g., Code Llama, ChatGPT-4o and ChatGPT 5). Additionally, we evaluate the performance of the MPI codes generated by ChatMPI in comparison with the ones generated by ChatGPT-4o and ChatGPT-5. The codes generated by ChatMPI provide up to a 4 × boost in performance by using better problem decomposition, communication patterns, and HPC techniques (e.g., communication avoiding).

Valero Lara, Pedro [ORNL] (ORCID:0000000214794310)↗
DOE OSTI2023

OC6 Phase III Definition Document

The objective of Phase III of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was to evaluate the accuracy of aerodynamic load predictions by offshore wind modeling tools for a floating offshore wind turbine as it experiences large surge-translational and pitch-rotational motion, as would occur during normal operation. A variety of models were examined in the project, with the goal of performing a three-way validation between engineering-level tools, higher-fidelity tools, and measurement data from two wind tunnel experimental campaigns. This document provides the necessary information to model and simulate the system studied in the OC6 Phase III project. This information is not only useful for project participants, but also future modelers who want to replicate this work for verification and validation purposes.

17 WIND ENERGY↗
NASA NTRS1993

High Reynolds number and turbulence effects on aerodynamics and heat transfer in a turbine cascade

Experimental data on pressure distribution and heat transfer on a turbine airfoil were obtained over a range of Reynolds numbers from 0.75 to 7.5 x 10 exp 6 and a range of turbulence intensities from 1.8 to about 15 percent. The purpose of this study was to obtain fundamental heat transfer and pressure distribution data over a wide range of high Reynolds numbers and to extend the heat transfer data base to include the range of Reynolds numbers encountered in the Space Shuttle main engine (SSME) turbopump turbines. Specifically, the study aimed to determine (1) the effect of Reynolds number on heat transfer, (2) the effect of upstream turbulence on heat transfer and pressure distribution, and (3) the relationship between heat transfer at high Reynolds numbers and the current data base. The results of this study indicated that Reynolds number and turbulence intensity have a large effect on both the transition from laminar to turbulent flow and the resulting heat transfer. For a given turbulence intensity, heat transfer for all Reynolds numbers at the leading edge can be correlated with the Frossling number developed for lower Reynolds numbers. For a given turbulence intensity, heat transfer for the airfoil surfaces downstream of the leading edge can be approximately correlated with a dimensionless parameter. Comparison of the experimental results were also made with a numerical solution from a two-dimensional Navier-Stokes code.

Yeh, Frederick C.↗
NASA NTRSDate not supplied

Identification and Study of Validation Level Test Cases for Computational Modeling of Non-Charring Ablators

Computational modeling of Thermal Protection System (TPS) materials, used for aerospace applications, provides numerous advantages in preliminary selection and design of a heatshield material and shape for atmospheric entry vehicles. However, to serve as a reliable tool for prediction of material thermal and ablative behavior, the modeling approach needs to be validated against real experimental and flight data, preferably at a range of applied conditions. The validation study is typically very complex as it requires reliable measured data not only for the material thermal response and surface recession, but also well characterized environmental conditions. The validation problem becomes even more complex when the material thermal response is dictated by multi-physics effects such as solid conduction, in-depth thermal decomposition, pyrolysis gas flow and chemical reactions. The multi-physics effects complicate not only the modeling effort, but also the experimental measurement for validation of various aspects of the highly coupled problem. In this study, an attempt is made to identify suitable experimental data that could serve as a source for validation of material thermal response modeling tools. To reduce the computational complexity, this study focuses only on non-charring ablators, where the material thermal response could be modeled with a single governing equation for solid conduction and the ablation is limited only to the surface of the material. With a well characterized and publicly available experimental data being sparse, the study is limited in presenting test cases for only three materials: camphor, graphite and FiberForm® in the sequence of increased modeling complexity. Graphite is a commonly used TPS material for aerospace applications, both for leading edges of high-speed vehicles and internal insulation of solid rocket motors. FiberForm® is a porous carbon pre-form used in preparation of the well known PICA material Tran et al. [1996]. Inclusion of camphor into the list is conditioned with the relative simplicity in modeling the material thermal and chemical response and the low-enthalpy flow environment. In addition, camphor has been used as a simple test material for study of flow transition behavior by Stock and Ginoux [1973] and assessment of a heatshield shape change at flight relevant conditions by Rotondi et al. [2022]. In this work, the identified experimental data was extracted from the public literature and test cases that yet have been published. As it was found from the review, not a single test case contains an exhaustive set of data that would validate every aspect of the material physics. However, in the data collected, various aspects of the material behavior can be still validated, such as surface and in-depth temperature, amount of recession and a shape change. The identified experimental data for each case is accompanied with a characterized flow environment and simulated boundary conditions predicted by a Data-Parallel Line Relaxation (DPLR) code Wright et al. [1998]. In addition, material thermal response numerical simulations in each test case are performed with Kentucky Aerothermodynamics and Thermal Response System (KATS-MR) Zibitsker et al. [2022] providing a comparative study and a sanity check for the proposed validation data. Sample results from the performed numerical study are shown below. Figure 1 shows distribution of surface heat flux and pressure values on a hemi-cylinder model made of FiberForm® and tested in HyMETS arc-jet facility. The results are shown for the high pressure condition among the two tests. Flow simulation was performed with DPLR code on a quarter of original geometry. In the figure, the quarter shape was mirrored across zx and xy planes to show the complete distribution. Figure 2 shows the material response results for the high pressure case (7500 Pa), simulated with KATS-MR and a comparison to the experimental data for the surface temperature and shape shape. The simulation was performed on a 2-D slice, extracted in the xy plane at the middle of the sample. Figure 3 shows the material response simulation for the low pressure case (3500 Pa) and a comparison to the experimental data for surface temperature and shape change.

ablation↗
DOE OSTI2021

Impurity transport studies at the HSX stellarator using active and passive CVI spectroscopy

The transport of carbon impurities has been studied in the helically symmetric stellarator experiment (HSX) using active and passive charge exchange recombination spectroscopy (CHERS). For the analysis of the CHERS signals, the STRAHL impurity transport code has been re-written in the python programming language and optimized for the application in stellarators. In addition, neutral hydrogen densities both along the NBI line of sight as well as for the background plasma have been calculated using the FIDASIM code. By using the basinhopping algorithm to minimize the difference between experimental and predicted active and passive signals, significant levels of impurity diffusion are observed. In this work, comparisons with neoclassical calculations from DKES/PENTA show that the inferred levels exceed the neoclassical transport by about a factor of four in the core and more than 100 times towards the plasma edge, thus indicating a high level of anomalous transport. This observation is in agreement with experimental heat diffusivites determined from a power balance analysis which exhibits strong anomalous transport as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Budding yeasts in the subphylum Saccharomycotina Genome sequencing and assembly

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available. By extracting mitochondrial reads from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades. Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements. As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

diversity↗
NASA NTRS1994

High Reynolds Number and Turbulence Effects on Turbine Heat Transfer

Experimental data on pressure distribution and heat transfer on a turbine airfoil were obtained over a range of Reynolds numbers from 0.75 to 7.0 x 10(exp 6) and a range of turbulence intensities from 1.8 to about 15%. The purpose of this study was to obtain fundamental heat transfer and pressure distribution data over a wide range of high Reynolds numbers and to extend the heat transfer data base to include the range or Reynolds numbers encountered in the Space Shuttle main engine turbopump turbines. The results of this study indicated that Reynolds number and turbulence intensity have a large effect on both the transition from laminar to turbulent flow and the resulting heat transfer. For a given turbulence intensity, heat transfer for all Reynolds numbers at the leading edge can generally be correlated with the Frossling number developed for lower Reynolds numbers. For a given turbulence intensity, heat transfer for the airfoil surfaces downstream of the leading edge can be approximately correlated with a dimensionless parameter. Comparisons of the experimental results were also made with a numerical solution from a two-dimensional Navier-Stokes code.

Yeh, Frederick C.↗
NASA NTRS2013

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

DeBonis, James R.↗
NASA NTRS2017

Getting to the Heart of Cardiovascular Risk Assessment in Astronauts for Exploration Class Missions

Since the beginning of manned spaceflight, NASA has recognized the potential risk of cardiovascular decrements due to stressors in the space environment. Of particular concern is the effect of space radiation on cardiovascular disease since astronauts will be exposed to higher levels of galactic cosmic rays outside the Earth's protective magnetosphere. To date, only a few studies have examined the effects of heavy ion radiation on cardiovascular disease, and at lower, space-relevant doses, the association between radiation exposure and cardiovascular pathology is more varied and unclear. Furthermore, other spaceflight conditions such as microgravity, circadian shifts, and confinement stress pose unique challenges in estimating the health risks that can be attributed to exposure to ionizing radiations. In this work, we review age, cause of mortality, and radiation exposure amongst early NASA astronauts in selection groups and discuss the limitations of assessing such a cohort when attempting to characterize the risk of space flight, including stressors such as space radiation and microgravity exposure, on cardiovascular health. METHODS: NASA astronauts in selection groups 1-7 were chosen and the comparison population was white men of the same birth cohort as drawn from data from the CDC Wonder Database and CDC National Center for Health Statistics Life Tables. Cause of death information was obtained from the Lifetime Surveillance of Astronaut Health program and deceased astronauts were classified based on ICD-10 codes: ischemic heart disease (IHD), stroke, cancer, acute occupational events, non-NASA accidents, and other/unknown. Expected years of life left and expected age at death were calculated for the cohort. RESULTS AND CONCLUSIONS: There were 32 deaths in this early astronaut population, 12 of which were due to accidents or acute occupational events that impacted lifespan considerably. The average age at death from these causes is 30 years lower than the average expected ~70 years of age in the general population. Remarkably, all 41 living early astronauts outlived our calculated expected age at death for members of their birth cohort; furthermore, 13 of the 20 deceased astronauts who did not die in NASA/non-NASA accidents exceeded this age. There was no difference in IHD between the astronaut cohort and the comparison population; therefore, it is not possible to associate IHD mortality with radiation in that astronaut cohort. As NASA looks toward future exploration-class missions, early astronaut cohorts provide a convenient option for assessing these risks and for developing mitigation strategies. However, many challenges still exist when assessing such limited evidence, including small cohort size, health and lifestyle confounders (such as smoking and drinking), the high accident mortality rate, and the fact that many of these astronauts are still alive, outliving many of their birth-cohort peers. Future analysis should include a longitudinal study, monitoring cases as they occur in the cohort. As this cohort is currently followed-up over time, and as more IHD cases are anticipated in a population of this age, this type of study is not as resource-intensive as would normally be the case.

Elgart, S. R.↗
DOE OSTI2020

Data-driven Alignment of the Tracker for the Mu2e Experiment

An alignment of the Mu2e Tracker is investigated using the Mu2e Offline Software and Millepede-II. Simulated misalignments of Tracker planes were inserted between Monte-Carlo event generation and before reconstruction using randomly generated $(x,y,z)$-translations and rotations of the rigid-body planes in the detector model. Cosmic ray muons, with no magnetic field present in the Detector Solenoid, will follow a straight-line trajectory. A maximum-likelihood fit is used to reconstruct a track and determine its parameters. Tracks traversing at least four planes, with ten registered straw hits, and a maximum time residual of $20$ ns are selected to be used in an alignment sample. The cuts were chosen to discard tracks with low information content, such as those with high levels of multiple scattering or those with low numbers of measurements contributing to a single track fit. Millepede-II produces the resulting alignment constants according to a simultaneous least-squares fit of the alignment parameters and the track parameters for all input tracks. The improved alignment constants are used to repeat the process and recover an increased number of suitable tracks, further improving the estimate of the alignment constants in subsequent iterations. A performance study of the process was conducted on a misaligned Tracker in the case of translation-only, and translation plus rotational misalignments. In both cases, the correct alignment constants were recovered to within statistical error after three global fit iterations. Additionally, a continuous integration system was developed for the \texttt{Mu2e/Offline} GitHub repository, introducing a convenient automated testing process for all proposed changes (``pull requests"). The tests developed included a full build of the software, followed by tests of the validation routines, geometry overlap checking, and multi-threading mode. Additionally, code quality checks and low-statistics physics comparison checks a re also available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Comparative verification of hydro-mechanical fracture behavior: Task G of international research project DECOVALEX–2023

Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, in this study, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.

58 GEOSCIENCES↗
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Data and code from: Selective oxidation of 3-hydroxypropionic acid to malonic acid over Pd/C: Mechanistic and kinetic study

Malonic acid (MA) is a high-value dicarboxylic acid with strong industrial demand, yet its current production heavily relies on petrochemical feedstocks. Here, we report the first systematic study for the sustainable production of MA via oxidation of 3-hydroxypropionic acid (3-HP) with a Pd/Carbon catalyst. The effects of oxidant type (O2 and H2O2), pH, and temperature on the reaction chemistry were comprehensively evaluated. Guided by experimental observation and DFT-calculated thermochemical energetics, reaction networks for 3-HP oxidation with both oxidants are proposed and validated through kinetic modeling. MA was identified as the primary oxidation product, while further oxidation yields acetic acid and oxalic acid. The kinetic model validated the network, displaying excellent agreement (R2 > 0.95). Kinetic observations also enabled a direct comparison between O2 and H2O2 and revealed their distinct behaviors. The model was further developed into a temperature-time map, providing insight into conditions that maximize MA production. Malonic acid selectivity of 56.9 % and yield of 50.5 % were achieved at 3 bar oxygen, equimolar NaOH:3-HP ratio, and 50 °C.

3-Hydroxypropionic Acid↗