Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulated”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Initial Demonstration of New Griffin Technologies for Simulating the Running-In Phase of Pebble Bed Reactors

Griffin is a reactor multiphysics modeling application based on MOOSE (Multiphysics Object-Oriented Simulation Environment) and specifically targeting transient modeling of advanced reactors. Griffin has been used recently to model pebble-bed reactors for the Nuclear Regulatory Commission (NRC) Office of Nuclear Regulatory Research and the Advanced Reactor Technology program. This modeling work has focused thus far on the direct calculation of equilibrium cores. This report documents an initial demonstration of a new running-in simulation capability. The new running-in capability is verified using the existing direct equilibrium core calculation capability. A simplified pebble-bed reactor model is then used to demonstrate the running-in simulation capability. This demonstration shows that Griffin is able to simulate years of operation during the running-in phase efficiently with each depletion step taking only several seconds. Two new technologies are also presented in this report which have been developed in Griffin that will be essential for improved accuracy both of the direct equilibrium core computation and the new running-in simulation capability. The first technology is an online cross section generation capability specifically targeted for pebble-bed reactors. This will improve the accuracy of the depletion calculation as the cross sections are generated at the exact core status. This also avoids the difficult step of pre-generating a separate standalone multigroup cross section set. Secondly, a newly implemented discretization for discontinuous finite element method (DFEM) SN transport in cylindrical (RZ) coordinates, which can be solved efficiently using the existing SN sweep solver, is discussed and some results are shown demonstrating the usefulness of the additional accuracy transport provides over a diffusion approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ensemble Simulation Techniques and Fast Randomized Algorithms

The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Report on the "Symposium on Advancements in Simulating Neutron Scattering Instruments and Experiments"

An international group of researchers gathered for the “Symposium on Advancements in Simulating Neutron Scattering Instruments and Experiments” at Oak Ridge National Laboratory (ORNL) on October 29- 30, 2024. The symposium highlighted state-of-the-art advancements in neutron scattering instrument simulations and the development of digital twins for scattering instruments. This document summarizes the recurring themes from the 19 talks and discussions. A key focus was on the role of neutronics simulations in the design, analysis, and performance optimization of neutron instruments. Presentations featured the most widely used software packages, such as McStas, McVine, Vitess, and MCNP. Another major theme addressed the analysis and interpretation of experimental neutron data, where materials models and simulations play a critical role. Talks on relevant software tools, including NCrystal, OCLimax, and Paradyse, were delivered by their respective developers. During the symposium, both talks and breakout sessions explored strategies to enhance the interplay between neutronics simulations and materials models. A key question emerged: How can we develop a RealTime Explainable Digital Twin of a neutron scattering experiment? Discussions include the idea that the Digital Twin concept has implications for instrument control and testing. Through this symposium, we gained valuable insights into the challenges and requirements for developing a Real-Time Explainable Digital Twin of a neutron scattering experiment. The discussions and presentations paved the way for future innovations in this rapidly advancing field.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Simulation of Deuterium Oxidation in a Room with a Fire

A series of experiments where deuterium was released in trace amounts into a room with a fire were performed at Sandia National Laboratories. This report describes the corresponding effort to model the test series using SIERRA/FUEGO. The objective of this modeling effort was to produce a simulation test matrix that can be utilized to help interpret the corresponding experiments and be used to assess the credibility of using the SIERRA/FUEGO simulations as a surrogate for real tritium reaction data for mock fire and release studies. The simulations assessed multiple configurations of fire size and location and overall found that there was low conversion of deuterium to oxidized deuterium. There are some notable differences between the simulation predictions and the experimental results, such as the simulations generally predicting higher amounts of deuterium conversion compared to the experiment. A strategy to address modeling uncertainties and unknowns is also given in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Report on Lightning Simulation Standards and Facilities

Lightning strikes pose a significant threat to various forms of infrastructure and electrical equipment. Aircraft and high explosives are of particular interest due to the massive ramifications of the device being compromised due to exposure to lightning damage. However, due to lightning being a varying weather phenomenon, testing a device’s vulnerability to lightning is generally not as simple as simply exposing it to real world lightning. To consistently test a device’s ability to survive the worst-case scenario of lightning exposure, the effects of the most severe lightning strikes must be able to be reproduced consistently and measurably. One such facility for this purpose is the Sandia Lightning Simulator (SLS), which is designed to simulate the current waveform of a severe lightning flash. However, the SLS’s design and age lead to difficulties in its operation and inconsistencies in its output. The methodology used in other lightning simulation facilities may offer information on how to improve the SLS’s operation, or at the very least, expose the SLS’s hidden strengths. This report explores lightning simulation, lightning testing standards, and lightning simulation facilities.

42 ENGINEERING↗
DOE OSTI2026

High-Resolution Mesoscale Simulations of Historical U.S. Gulf Coast Hurricanes

This dataset provides coupled mesoscale and microscale (large-eddy simulation; LES) atmospheric simulations of five historical U.S. Gulf Coast hurricanes. Meso-microscale coupled simulations of Hurricanes Harvey (2017), Nate (2017), Michael (2018), Laura (2020), and Delta (2020) were performed with the Weather Research and Forecasting (WRF) model v4.1.5. The mesoscale component captures the multi-day evolution of each storm's track, intensity, and large-scale structure across the Gulf of Mexico, while the LES component explicitly resolves the dominant turbulent eddies in the hurricane boundary layer at a horizontal grid spacing of 33.33 m. Together, these simulations characterize tropical cyclone wind fields across atmospheric scales ranging from ~100 km down to ~100 m, spanning storm intensities from Category 2 through Category 4 on the Saffir–Simpson scale. This dataset contains the mesoscale simulations, separated by each storm.

17 WIND ENERGY↗
DOE OSTI2026

High-Resolution Regional Atmosphere–Ocean–Wave Coupled Simulations of Hurricane Henri (2021)

To explore the integrated effects of ocean and ocean surface wave related physical processes on tropical cyclone simulations, a set of three model simulations is performed. * In experiment 'A', the event is modeled using the Weather Research Forecasting (WRF) model alone with prescribed Sea Surface Temperature (SST) at 6-hour intervals. * In experiment ‘AO,’ WRF is coupled with the Finite Volume Community Ocean Model (FVCOM), enabling variable exchange between atmosphere and ocean, but without considering ocean surface wave-related physical processes. * In experiment ‘AOW’, WRF, FVCOM, and Simulating WAves Nearshore (SWAN) exchange variables with each other every hour through the OASIS3-MCT Coupler to allow direct and indirect atmosphere-ocean-wave interactions. * Observational data are also included in this dataset (Dropsonde, HRD-Radar, NDBC_wave). All simulations are initialized at 18:00 UTC on August 19, 2021, within a domain encompassing the western North Atlantic Ocean. The atmospheric domain features a horizontal resolution of 3 km. The ocean domain, which covers a substantial portion of the WRF ocean domain, employs an unstructured triangular grid with resolutions ranging from 3 km near the coast to 9 km in the open ocean, effectively resolving the complex coastline of the U.S. Northeast Coast. Initial and boundary conditions for the atmosphere model are obtained from the 6-hourly 0.25° NCEP (National Centers for Environmental Prediction) Global Forecast System (GFS; NCEP, 2015) data. These CSV files are derived from the NetCDF files in the c0 dataset. Unlike the original format, where geographic coordinates were stored in a separate file, each CSV now embeds the corresponding latitude and longitude values alongside the measured or simulated variables.

17 WIND ENERGY↗
DOE OSTI2025

Damping Wing-Like Features in the Spectra of High Redshift Quasars: a Challenge for Fully-Coupled Simulations

Recently, several observational detections of damping-wing-like features at the edges of “dark gaps” in the spectra of distant quasars (the “Malloy-Lidz effect”) have been reported, rendering strong support for the existence of “neutral islands” in the universe at redshifts as low as z<5.5. We apply the procedure from one of these works, Zhu et al 2024, to the outputs of fully coupled cosmological simulations from two recent large projects, “Cosmic Reionization On Computers” (CROC) and “Thesan”. Synthetic spectra in both simulations have statistics of dark gaps similar to observations, but do not exhibit the damping wing features. Moreover, a toy model with neutral islands added “by hand” only reproduces the observational results when the fraction of neutral islands among all dark gaps approaches 90%. I.e., simulations and observations appear to produce two distinct “populations” of dark gaps. In addition, in the simulations, the neutral islands at z=5.9 should be short-lived and should not extend to z<5.5. A plausible explanation for this discrepancy is that both simulations underestimate the fluctuations in the photoionization rate and, hence, are missing a population of long-lived neutral islands, located in the large downward fluctuations of the photoionization rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.↗
DOE OSTI2024

Evaluation of Global Fire Simulations in CMIP6 Earth System Models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of the El Niño–Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in CMIP5: (1) the simulated global burned area being less than half of that of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past 2 decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere midlatitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet–dry conditions, indicating the need to improve fuel wet-ness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest a postprocessing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products.

Wildfire, Earth system models↗
DOE OSTI2025

Assessing modifications to the Abdul-Razzak and Ghan aerosol activation parameterization (version ARG2000) to improve simulated aerosol–cloud radiative effects in the UK Met Office Unified Model (UM version 13.0)

The representation of aerosol activation is a key source of uncertainty in global composition-climate model simulations of aerosol–cloud interactions. The Abdul-Razzak and Ghan (ARG) activation parameterization is used in several global and regional models that employ modal aerosol microphysics schemes. In this study, we investigate the ability of the ARG parameterization to reproduce simulations with a cloud parcel model and find its performance is sensitive to the geometric standard deviations (widths) of the lognormal aerosol modes. We recommend adjustments to three constant parameters in the ARG equations, which improve the performance of the parameterization for small mode widths and its ability to simulate activation in polluted conditions. For the accumulation mode width of 1.4 used in the Met Office Unified Model (UM), the modifications decrease the mean bias in the activated fraction of aerosols compared to a cloud parcel model from −6.6 % to +1.2 %. We implemented the improvements in the UM and compared simulated global cloud droplet concentrations with satellite observations. The simulated cloud radiative effect changes by −1.43 W m −2 (6 %) and aerosol indirect radiative forcing over the industrial period changes by −0.10 W m −2 (10 %).

Ghosh, Pratapaditya [Carnegie Mellon University, P↗
DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis – Simulated Wind

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wind turbine. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen . While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the simulated wind energy profiles, NLR used OpenFAST to simulate a 3.4-MW International Energy Agency (IEA) reference wind turbine. The hour-long wind energy profiles varied over wind turbulence intensity (Class A or Class C) and average wind speed (5, 7, or 9 m/s). To match the power limits of the 1.25-MW electrolyzer and 3.4-MW IEA wind turbine most effectively and to maximize the efficiency of hydrogen production at a given average wind speed, the profiles were sometimes scaled by two times. This means that, in some cases, the experimental setup assumed two 1.25-MW electrolyzers were coupled with the wind turbine, representing a total maximum electrolysis load of 2.5 MW. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wind turbine electrolysis experiment and is formatted as follows: {technology}-{average wind speed}-{turbulence class}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “windIEA3.4-5ms-C_2-400.zip” represents the hour-long experiment using the IEA 3.4-MW turbine, subjected to an average wind speed of 5 m/s and Class C wind turbulence, and connected to two 1.25-MW electrolyzers with the power supply set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wind turbine power. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all simulated wind experiments combined into one dataset labeled "combined_wind_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis .

08 HYDROGEN↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗
DOE OSTI2024

Constraining Cosmology with Simulation-based inference and Optical Galaxy Cluster Abundance

We test the robustness of simulation-based inference (SBI) in the context of cosmological parameter estimation from galaxy cluster counts and masses in simulated optical datasets. We construct ``simulations'' using analytical models for the galaxy cluster halo mass function (HMF) and for the observed richness (number of observed member galaxies) to train and test the SBI method. We compare the SBI parameter posterior samples to those from an MCMC analysis that uses the same analytical models to construct predictions of the observed data vector. The two methods exhibit comparable performance, with reliable constraints derived for the primary cosmological parameters, ($\Omega_m$ and $\sigma_8$), and richness-mass relation parameters. We also perform out-of-domain tests with observables constructed from galaxy cluster-sized halos in the Quijote simulations. Again, the SBI and MCMC results have comparable posteriors, with similar uncertainties and biases. Unsurprisingly, upon evaluating the SBI method on thousands of simulated data vectors that span the parameter space, SBI exhibits worsened posterior calibration metrics in the out-of-domain application. We note that such calibration tests with MCMC is less computationally feasible and highlight the potential use of SBI to stress-test limitations of analytical models, such as in the use for constructing models for inference with MCMC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

First Estimation of Model Parameters for Neutrino-Induced Nucleon Knockout Using Simulation-Based Inference

To enable an accurate determination of oscillation parameters, accelerator-based neutrino experiments require detailed simulations of nuclear interaction physics in the GeV regime. While substantial effort from both theory and experiment is currently being invested to improve the fidelity of these simulations, their present deficiencies typically oblige experimental collaborations to resort to empirical tuning of simulation model parameters. As the precision requirements of the field continue to become more stringent, machine learning techniques may provide a powerful means of handling corresponding growth in the complexity of future neutrino interaction model tuning exercises. To study the suitability of simulation-based inference (SBI) for this physics application, in this paper we revisit a tuned configuration of the GENIE neutrino event generator that was originally developed by the MicroBooNE collaboration. Despite closely reproducing the adopted values of four physics parameters when confronted with the tuned cross-section predictions as input, we find that our trained SBI algorithm prefers modestly different values (within MicroBooNE's assigned uncertainties) and achieves slightly better goodness-of-fit when inference is run on the experimental data set originally used by MicroBooNE. We also find that our trained algorithm can create a fair approximation of an alternative neutrino scattering simulation, NuWro, that shares only a subset of its physics model parameters with GENIE.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗