Search NASA⌕ Search

SEARCH · Search NASA

Results for “function approximation”

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

Mock data sets for the Eboss and DESI Lyman-α forest surveys

We present a publicly-available code to generate sets of mock Lyman-α (Lyα) forest data that have realistic large-scale correlations including those due to the Baryonic Acoustic Oscillations (BAO). The primary purpose of these mocks is to test the analysis procedures of the Extended Baryon Oscillation Survey (eBOSS) and the Dark Energy Spectroscopy Instrument (DESI) surveys. The transmitted flux fraction, F(λ), of background quasars due to Lyα absorption in the intergalactic medium (IGM) is simulated using the Fluctuating Gunn-Petterson Approximation (FGPA) applied to Gaussian random fields produced through the use of fast Fourier transforms (FFT). The output includes the IGM-Lyα transmitted flux fraction along quasar lines of sight and a catalog of high-column-density systems appropriately placed at high-density regions of the IGM. This output serves as input to additional code that superimposes the IGM tranmission on realistic quasar spectra, adds absorption by high-column-density systems and metals, and simulates instrumental transmission and noise. Redshift space distortions (RSD) of the flux correlations are implemented by including the large-scale velocity-gradient field in the FGPA resulting in a correlation function of F(λ) that can be accurately predicted. One hundred realizations have been produced over the 14,000 deg 2 DESI survey footprint with 100 quasars per deg 2 . The analysis of these realizations shows that the correlations of F(λ) follows the prediction within the accuracy of eBOSS survey. Here, the most time-consuming part of the mock production occurs before application of the FGPA, and the existing pre-FGPA forests can be used to easily produce new mock sets with modified redshift-dependent bias parameters or observational conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tunneling magnetoresistance in MgO tunnel junctions with Fe-based leads in empirically corrected density functional theory

The minority-spin Fe/MgO interface states are at the Fermi level in density functional theory (DFT), but experimental evidence and GW calculations place them slightly higher in energy. This small shift can strongly influence tunneling magnetoresistance (TMR) in junctions with a thin MgO barrier and its dependence on the concentration of Co in the electrodes. Here, in this study, an empirical potential correction to DFT is introduced to shift the interface states up to match the tunnel spectroscopy data. With this shift, TMR in Fe/MgO/Fe junctions exceeds 800% and 3000% at 3 and 4 monolayers (ML) of MgO, respectively. We further consider the effect of alloying of the Fe electrodes with up to 30% Co or 10% V, treating them in the coherent potential approximation (CPA). Alloying with Co broadens the interface states and brings a large incoherent minority-spin spectral weight to the Fermi level. Alloying with V brings the minority-spin resonant states close to the Fermi level. However, in both cases the minority-spin spectral weight at the Fermi level resides primarily at the periphery of the Brillouin zone, which is favorable for spin filtering. Using convolutions of k II -resolved barrier densities of states calculated in CPA, it is found that TMR is strongly reduced by alloying with Co or V but still remains above 500% at 4 ML of MgO up to 30% of Co or 5% V. At 5 ML, the TMR increases above 1000% in all systems considered. However, while TMR declines sharply with increasing bias up to 0.2 eV in the tunnel junctions with pure Fe leads, it remains almost constant up to 0.5 eV if leads are alloyed with Co.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Heterostructural Alloy Phase Diagram for (Cd 1-x Zn x ) 3 As 2

Alloying the topological semimetal Cd 3 As 2 with Zn 3 As 2 provides a potential route for controlling the electronic properties. We predict the alloy phase diagram from first-principles calculations, considering that both end members have a crystal structure derived from the antifluorite lattice, but with different arrangements of the unoccupied cation sites. To overcome the limitations of the regular solution approximation and to include short-range order effects, we perform Monte Carlo simulations, parameterize the temperature dependence of the mixing enthalpy ΔH m , and perform thermodynamic integration of the free energy. The resulting phase diagram exhibits features that are unique to heterostructural alloy systems and provides computational predictions of solubility limits and composition ranges that are stable against spinodal decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Sensitivity of Variational Bayesian Neural Network Performance to Hyperparameters

In scientific applications, predictive modeling is often of limited use without accurate uncertainty quantification (UQ) to indicate when a model may be extrapolating or when more data needs to be collected. Bayesian Neural Networks (BNNs) produce predictive uncertainty by propagating uncertainty in neural network (NN) weights and offer the promise of obtaining not only an accurate predictive model but also accurate UQ. However, in practice, obtaining accurate UQ with BNNs is difficult due in part to the approximations used for model training (such as those made in variational inference) and in part to the need to choose a suitable set of hyperparameters; these hyperparameters outnumber those needed for traditional NNs and often have opaque effects on the results. We aim to shed light on the effects of hyperparameter choices for variational BNNs by performing a global sensitivity analysis of variational BNN performance under varying hyperparameter settings. Our results indicate that many of the hyperparameters interact with each other to affect both predictive accuracy and UQ. For improved usage of variational BNNs in real-world applications, we suggest that thorough hyperparameter tuning, including tuning of prior hyperparameters and loss function parameters, is essential for accurate UQ in variational BNNs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A novel transient infrared imaging method for non-intrusive, low-cost, fast, and accurate air leakage detection in building envelopes

Air leakage through the building envelope in the U.S. accounts for about four quads of energy annually, costing approximately $40 billion per year. However, a high-fidelity and non-intrusive method to detect air leakage has not been demonstrated to date. Here, in this paper, we propose a novel non-intrusive and low-cost method called Transient Infrared (IR) Imaging (TIRI) that can rapidly and accurately identify air leakage locations and relative rates on building envelopes. When the interior and exterior temperatures are different, and a small internal pressure pulse is created by HVAC, the temperature at locations with air leakages will change rapidly, while the areas without a leakage do not change. Based on a heat transfer model, we have derived the temperature change as a function of time after the HVAC is turned on. By tracking the temperature change, which depends on leakage rate and size, we have obtained the air leakage map in the case studies. Using an exterior door as an example, we took transient IR images in different seasons and different times of the day, and successfully obtained the leakage map in all the scenarios. Successfully obtained the air leakage map even when the indoor-outdoor air temperature difference is as small as 2 °C. We have also realized a detection speed of 10s and demonstrated that this method also worked for windows, which have mirror-like IR reflections. Our TIRI method will accelerate the improvement of airtightness in buildings, save building energy, and help reduce greenhouse gas emissions.

42 ENGINEERING↗
DOE OSTI2026

Observed Temperature and Moisture Advection Characteristics in the Southern Great Plains

Until recently, advection profiles in the planetary boundary layer have been difficult to quantify from observations. As a result, not much is known about its basic characteristics, like magnitude and diurnal cycles, or how these differ as a function of large-scale environments. Here, to provide insight into advective characteristics, a Green’s Theorem–based method for calculating profiles of advection in the lowest 3 km from an array of ground-based thermodynamic and kinematic profiling instruments has been applied to 2 years of observations from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observing site. Since advection largely depends on synoptic-scale forcing, a self-organizing map (SOM) was used to classify days based on mean sea level pressure analyses; from that, the diurnal evolution of vertical profiles of advection for different synoptic environments was quantified. The overall mean magnitude of the potential temperature advection is approximately ±3 K h −1 , while the moisture advection is approximately ±1.5 g kg −1 h −1 , with substantial variability in the sign and magnitude of the advective tendencies at different heights and throughout the diurnal cycle. Advection magnitude is strongly connected to the strength of synoptic-scale forcing, which varies depending on the time of year. Variability in advection is larger for potential temperature than it is for moisture, and strongly forced environments, like midlatitude cyclones, feature greater variability in the advection magnitude than weakly forced or quiescent environments.

advection↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Near-half-metallic state in the half-Heusler PtMnSb film on a III-V substrate

The interplay between half-metallic ferromagnetism and spin-orbit coupling within the inversion-symmetry-broken structure of half Heuslers provides an ideal platform for various spintronics functionalities. Taking advantage of good lattice matching, it is highly desired to epitaxially integrate promising Heuslers into III-V semiconductor-based devices. PtMnSb is one of the first half Heuslers predicted to be an above-room-temperature half-metal with large spin orbit coupling, however, its half-metallicity and potential as a spintronics material has remained elusive due to lack of high quality samples. Here, in this study, we demonstrate epitaxial growth of single crystal PtMnSb(001) film on GaSb(001) substrates using molecular beam epitaxy. Direct observation of the band structure via angle-resolved photoemission spectroscopy and many-body perturbation theory within the quasi-particle self-consistent GW approximation (QPGW) reveal that PtMnSb hosts rather a near-half-metallic state with both spin bands crossing the Fermi level and with high spin polarization over 90%. Temperature dependence of magnetization also shows an anomalous enhancement below 60 K, which can be associated with the development of such a near-half-metallic state at low temperatures. Epitaxial growth of high crystalline PtMnSb on a III-V paves the way for systematic clarification of its spin transport properties with fine-tuning of strain in heterostructure devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assembly, comparative analysis, and utilization of a single haplotype reference genome for soybean

Cultivar Williams 82 has served as the reference genome for the soybean research community since 2008, but is known to have areas of genomic heterogeneity among different sub-lines. This work provides an updated assembly (version Wm82.a6) derived from a specific sub-line known as Wm82-ISU-01 (seeds available under USDA accession PI 704477). The genome was assembled using Pacific BioSciences HiFi reads and integrated into chromosomes using HiC. The 20 soybean chromosomes assembled into a genome of 1.01Gb, consisting of 36 contigs. The genome annotation identified 48 387 gene models, named in accordance with previous assembly versions Wm82.a2 and Wm82.a4. Comparisons of Wm82.a6 with other near-gapless assemblies of Williams 82 reveal large regions of genomic heterogeneity, including regions of differential introgression from the cultivar Kingwa within approximately 30 Mb and 25 Mb segments on chromosomes 03 and 07, respectively. Additionally, our analysis revealed a previously unknown large (> 20 Mb) heterogeneous region in the pericentromeric region of chromosome 12, where Wm82.a6 matches the ‘Williams’ haplotype while the other two near-gapless assemblies do not match the haplotype of either parent of Williams 82. In addition to the Wm82.a6 assembly, we also assembled the genome of ‘Fiskeby III,’ a rich resource for abiotic stress resistance genes. A genome comparison of Wm82.a6 with Fiskeby III revealed the nucleotide and structural polymorphisms between the two genomes within a QTL region for iron deficiency chlorosis resistance. The Wm82.a6 and Fiskeby III genomes described here will enhance comparative and functional genomics capacities and applications in the soybean community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Sac1 links phosphoinositide turnover to cryptococcal virulence

Cryptococcus neoformans is an environmentally acquired fungal pathogen that causes over 140,000 deaths per year. Cryptococcal infection occurs when infectious particles are deposited into the lung, where they encounter host phagocytic cells. C. neoformans may be engulfed by these phagocytes, an important step of infection that leads to outcomes ranging from termination of infection to cryptococcal dissemination. To study this critical process, we screened approximately 4,700 cryptococcal gene deletion mutants for altered uptake, using primary mouse and human phagocytic cells. Among the hits of these two screens, we identified 93 mutants with perturbed uptake in both systems, as well as others with differences in uptake by only one cell type. We further screened the hits for changes in thickness of the capsule, a protective polysaccharide layer around the cell which is an important cryptococcal virulence factor. The combination of our three screens yielded 45 mutants, including one lacking the phosphatidylinositol-4-phosphate phosphatase Sac1. In this work, we implicate Sac1 in both host cell uptake and capsule production. We found that sac1 mutants exhibit lipid trafficking defects, reductions in secretory system function, and changes in capsule size and composition. Many of these changes occur specifically in tissue culture media, highlighting the role of Sac1 phosphatase activity in responding to the stress of host-like conditions. Overall, these findings show how genome-scale screening can identify cellular factors that contribute to our understanding of cryptococcal biology and demonstrate the role of Sac1 in determining fungal virulence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Akiachak Energy Efficiency Retrofit Project

The goal of the project is to reduce the overall energy use of the Akiachak Native Community (ANC) by implementing energy efficiency measures in five high-use Tribal buildings. This project will have the following outcomes: Projected annual energy savings of $17,369; projected annual reduction in fuel oil #1 of 1,200 gallons and electricity of 17,751 kWh; annual reduction in carbon dioxide emissions of approximately 60,340 pounds/year. ANC will install energy efficiency measures in the Laundry, Tribal Indian Reorganization Act (IRA) Office, Clinic, Daycare, and Police Station. ANC obtained energy audits on these buildings in 2018, and this project will implement high-payback recommendations such as replacing lighting with LEDs, installing setback thermostats and occupancy sensors, replacing furnaces with more efficient models, replacing the circulation pumps with variable speed ones, air tightening, and adding insulation. Buildings will see energy cost reductions from 15% to 40%. These retrofits will help build ANC’s long-term vision for sustainable energy usage and address the first goal of the Tribal IRA Council’s Energy Efficiency and Conservation Strategy, to “create and maintain functionally appropriate, sustainable, accessible, high quality tribal infrastructure and facilities.” ANC intends to replicate this project by using the resulting energy savings to address audit recommendations in other buildings as well as to demonstrate the value of energy efficiency to community members. Other outcomes will include an increase in community resiliency, reduced dependence on outside shipments of fuel oil, training for maintenance staff, and no-touch control of building appliances to reduce transmission of diseases such as COVID-19.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

A Performance and Energy Study of GPU-Resident Preconditioners for Conjugate Gradient Solvers: In the Context of Existing and Novel Approaches

Optimizing a particular subprogram out of the set of Basic (sparse) Linear Algebra Subprograms (BLAS) for a given architecture is a common topic of research. In applications, however, these BLAS functions rarely appear in isolation; usually, many of them are used together, in various combinations and with varying inputs. As the need to solve a large, sparse linear system is ubiquitous throughout HPC applications, linear solvers constitute a realistic, sufficiently complex and well-defined representative use case for composite BLAS routines. To this end, based on a representative set of matrices drawn from a diverse set of fields, we present a framework to study, from the performance and energy perspective, the efficacy of GPU- resident parallel Conjugate Gradient (CG) linear solver with different preconditioner options, including Gauss-Seidel, Jacobi, and incomplete Cholesky. We also propose a novel GPU-based preconditioner, in which the triangular solves are approximated by an iterative process. The development of this preconditioner was motivated by solving large graph Laplacian linear systems, for which the existing preconditioners either perform slow on GPU-based platforms or are not applicable. We compare the performance of these preconditioners on different hardware accelerator architectures, i.e., AMD MI250X, MI100, Nvidia A100, V100, and Jetson. Our experiments reveal performance trade-offs and provide information on how to select the best strategy for the given linear system, dictated by its properties, and the platform of interest. We demonstrate the application of our novel preconditioner for solving CG and graph Laplacian systems. Overall, the framework can be utilized as a benchmark to guide informed decisions in choosing a specific preconditioner, i.e., whether it is better to rely on the performance of a triangular solver or on the performance of sparse matrix-vector product. Finally, by considering power consumption to solve the linear systems, we report the energy footprint for the solvers.

Preconditioned Conjugate Gradient, GPUs, iterative↗
DOE OSTI2025

Development of Li–S Battery Cells with High Energy Density and Long Cycling Life

Lithium–sulfur (Li–S) batteries are considered one of the most promising next-generation energy storage technologies for electric vehicles, owing to their high energy density (up to three times greater than current lithium-ion batteries) and the low cost of sulfur (approximately 100 times cheaper than conventional cobalt oxide). However, several challenges hinder their commercial viability, including the electronically and ionically insulating nature of sulfur, polysulfide dissolution, and lithium dendrite growth. These issues significantly limit battery energy density and cycle life and must be addressed to enable successful commercialization. Current strategies involving porous carbon structures, additives, and electrocatalyst engineering have not fully resolved the polysulfide dissolution problem. While polysulfide-free cathode materials offer a promising solution, they often suffer from low sulfur content, leading to reduced discharge capacity and sluggish reaction kinetics. To overcome these limitations, we have developed novel soluble-polysulfide-free sulfur cathode active materials with a high sulfur content (>50 wt%). By integrating these cathodes with functional binders, optimized electrolyte formulations, and refined electrode fabrication techniques, we demonstrate Li–S batteries with both high energy density and long cycle life. This advancement brings the realization of low-cost, high-energy Li–S batteries significantly closer to practical application.

25 ENERGY STORAGE↗
DOE OSTI2025

Muon Time-of-Flight studies for cosmic background rejection in the Short Baseline Near Detector

The Short-Baseline Neutrino (SBN) program at Fermilab is a cutting-edge project in experimental neutrino physics. One of its main goals is to systematically investigate the possible existence of eV-scale sterile neutrinos. This phenomenon has been hypothesized to explain some anomalies found in short-range experiments and, if confirmed, would imply a substantial extension of the Standard Model. SBN also offers an important opportunity to deepen the understanding of neutrino-nucleus interactions in the GeV energy range, through the use of Liquid Argon Time Projection Chambers (LArTPC) detectors, a fundamental technology also for the future DUNE experiment. The SBN experimental infrastructure consists of three detectors aligned along the Booster Neutrino Beamline at Fermilab. Among them, the detector located closest to the neutrino source, SBND (Short-Baseline Near Detector), positioned approximately 110 meters from the target, plays a key role in directly characterizing the initial neutrino flux. This allows for a direct comparison with the measurements from the far detector, ICARUS, located about 600 meters from the source, in order to search for potential signs of anomalous neutrino oscillations. My master's thesis focuses on the commissioning and characterization activities of the SBND detector, with particular reference to the Cosmic Ray Tagger (CRT). The CRT is a subsystem for identifying and rejecting events produced by cosmic rays, which constitute the main source of background for surface experiments like SBND. The activity began with the commissioning of the final components of the detector, as well as their validation to verify their correct functioning and signal acquisition. A central part of my work involved studying the veto efficiency of the CRT system, analyzing the rate of cosmic ray-induced events to quantify any loss of neutrino-induced events caused by cosmic background. This allowed for a more precise evaluation of the systematic impact of the CRT on the useful physics sample. A further phase of my analysis involved an in-depth study of the temporal correlation between the CRT signals and those acquired by the LArTPC's internal photodetector system, consisting of photomultiplier tubes and X-ARAPUCA devices. The objective is to explore the possibility of using combined temporal information as an additional criterion for discriminating between cosmic signals and signals genuinely due to neutrino interaction. Preliminary results indicate the presence of characteristic temporal signatures that could be exploited to improve event selection and increase the purity of the neutrino-induced sample. These methodologies will certainly contribute to the optimization of SBND analysis strategies and, more generally, to a better understanding of background mechanisms in next-generation LArTPC experiments.

Corallo, Annalea [Ferrara U.]↗
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↗