SEARCH · Engineering Papers

Results for “generator parameter estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Correlation approach for quality assurance of additive manufactured parts based on optical metrology

Surface topography and surface finish are two significant factors for evaluating the quality of products in additive manufacturing (AM). AM parts are fabricated layer by layer, which is quite different from traditional formative or subtractive methods. Despite rapid progress in additive manufacturing and associated optical metrology for quality control and in-situ monitoring, limited research has been conducted to investigate the reliability of 3D surface measurement data. The surface topologies scanned by multiple optical systems demonstrated significant differences due to varying sampling mechanisms, resolutions, system noises, etc. The 3D datasets should be trustworthy in order to extract parameters for quality assurance or feedback control from 3D surface measurements. In this paper, we set up new standards to evaluate the reliability of 3D surface measurement data and analyze the variation in the topographical profile. In this study, two non-contact optical methods based on Focus Variation Microscopy (FVM) and Structured Light System (SLS) were adopted to measure the surface topography of the target components. The two optical metrology systems generated two entirely different point cloud datasets. Statistical methods were applied to test the difference between the data obtained from the two systems. By using a data analytics approach for comparison, it was found that the surface roughness estimated from the point cloud data sets of FVM and SLS has no significant difference, though the point cloud data sets were completely different. This paper provides standard validation approach to evaluate the plausibility of metrology data from in-situ real-time surface analysis for process planning of AM.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Synchrotron X-ray diffraction analysis of constituent phases in transition joint between nickel alloy 738LC and a MnFeCoNiCu alloy

Synchrotron x-ray diffraction (XRD) analysis was performed on transition joints between a single-phase MnFeCoNiCu alloy and Ni-base Alloy 738LC to efficiently identify the constituent phases across the interface, with different levels of material intermixing generated by laser-welding with variable power. Here, the intermixing extent was quantified by postmortem energy dispersive x-ray spectroscopy mapping. Diffraction-based analyses on complex composition spaces with limited prior data present significant challenges because elemental substitution in both disordered and ordered phases is more extensive than in conventional alloy spaces, which may effectuate relatively large shifts in the observed lattice parameters that convolute the analysis. Therefore, thermodynamic simulations and crystallographic literature data were employed to construct a system-specific diffraction library of twelve prospective phases for the composition space investigated. Subsequently, for predicted disordered cubic phases, statistical hard-sphere models were established to estimate the lattice parameters and predict diffraction peak positions for inclusion in the library. The library was then employed to analyze diffraction profiles measured from the variably intermixed transition joints, with focus on accounting for both high and low-intensity peaks. 99.0 % of diffraction peaks with relative intensity greater than 0.001 were assigned to phases from the system-specific library, exemplifying rigorous peak accounting and indicating that no unexpected phases were present. Up to six of the twelve library phases were experimentally found in the transition joints. The lattice parameters predicted by the statistical hard-sphere model based upon thermodynamic simulations agree reasonably well with the measured values for the disordered FCC matrix phase.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2008

Defining Constellation Suit Helmet Field of View Requirements Employing a Mission Segment Based Reduction Process

Field of view has always been a design feature paramount to helmet design, and in particular space suit design, where the helmet must provide an adequate field of view for a large range of activities, environments, and body positions. For Project Constellation, a slightly different approach to helmet requirement maturation was utilized; one that was less a direct function of body position and suit pressure and more a function of the mission segment in which the field of view is required. Through taxonimization of various parameters that affect suited FOV, as well as consideration for possible nominal and contingency operations during that mission segment, a reduction process was able to condense the large number of possible outcomes to only six unique field of view angle requirements that still captured all necessary variables without sacrificing fidelity. The specific field of view angles were defined by considering mission segment activities, historical performance of other suits, comparison between similar requirements (pressure visor up versus down, etc.), estimated requirements from other teams for field of view (Orion, Altair, EVA), previous field of view tests, medical data for shirtsleeve field of view performance, and mapping of visual field data to generate 45degree off-axis field of view requirements. Full resolution of several specific field of view angle requirements warranted further work, which consisted of low and medium fidelity field of view testing in the rear entry ISuit and DO27 helmet prototype. This paper serves to document this reduction progress and followup testing employed to write the Constellation requirements for helmet field of view.

McFarland, Shane M.↗

NASA NTRS2019

Constraining Stellar Parameters and Atmospheric Dynamics of the Carbon AGB star V Oph

Molecules and dust produced by the atmospheres of cool evolved stars contribute to a significant amount of the total material found in the interstellar medium. To understand the mechanism behind the mass loss of these stars, it is of pivotal importance to investigate the structure and dynamics of their atmospheres. Our goal is to verify if the extended molecular and dust layers of the carbon-rich asymptotic giant branch (AGB) star V Oph, and their time variations, can be explained by dust-driven winds triggered by stellar pulsation alone, or if other mechanisms are in play. We model V Oph mid-infrared interferometric VLTI-MIDI data (8–13 μm), at phases 0.18, 0.49, and 0.65, together with literature photometric data, using the latest-generation self-consistent dynamic atmosphere models for carbon-rich stars: DARWIN. We determine the fundamental stellar parameters: T(eff) = 2600 K, L(bol) = 3585 L⊙, M = 1.5 M⊙, C/O = 1.35, Mass flow rate= 2.50 x 10(exp -6) M⊙/yr. We calculate the stellar photospheric radii at the three phases: 479, 494, 448 R⊙; and the dust radii: 780, 853, 787 R⊙. The dynamic models can fairly explain the observed N-band visibility and spectra, although there is some discrepancy between the data and the models, which is discussed in the text. We discuss the possible causes of the temporal variations of the outer atmosphere, deriving an estimate of the magnetic field strength, and computing upper limits for the Alfvén waves velocity. In addition, using period–luminosity sequences, and interferometric modeling, we suggest V Oph as a candidate to be reclassified as a semi-regular star.

Gioia Rau↗

NASA NTRS2021

Synthesizing Disparate LiDAR and Satellite Datasets through Deep Learning to Generate Wall-to-Wall Regional Inventories for the Complex, Mixed-Species Forests of the Eastern United States

Light detection and ranging (LiDAR) has become a commonly-used tool for generating remotely-sensed forest inventories. However, LiDAR-derived forest inventories have remained uncommon at a regional scale due to varying parameters among LiDAR data acquisitions and the availability of sufficient calibration data. Here, we present a model using a 3-D convolutional neural network (CNN), a form of deep learning capable of scanning a LiDAR point cloud, combined with coincident satellite data (spectral, phenology, and disturbance history). We compared this approach to traditional modeling used for making forest predictions from LiDAR data (height metrics and random forest) and found that the CNN had consistently lower uncertainty. We then applied the CNN to public data over six New England states in the USA, generating maps of 14 forest attributes at a 10 m resolution over 85% of the region. Aboveground biomass estimates produced a root mean square error of 36 Mg ha−1 (44%) and were within the 97.5% confidence of independent county-level estimates for 33 of 38 or 86.8% of the counties examined. CNN predictions for stem density and percentage of conifer attributes were moderately successful, while predictions for detailed species groupings were less successful. The approach shows promise for improving the prediction of forest attributes from regional LiDAR data and for combining disparate LiDAR datasets into a common framework for large-scale estimation.

Elias Ayrey↗

DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗

NASA NTRS2003

First Year Wilkinson Microwave Anisotropy Probe(WMAP) Observations: Data Processing Methods and Systematic Errors Limits

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of l/f noise. We infer and correct for a small (approx. 1 %) transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

Hinshaw, G.↗

NASA NTRS2018

Improving the SMAP Level-4 Soil Moisture Product

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level 4 Soil Moisture (L4 SM) product. The L4 SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter. The L4 SM data are published with a mean latency of approx. 2.5 days from the time of observation and provide global, three-hourly, 9 km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. In 2018, the product was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high quality space-borne remote sensing observations. Land cover inputs were updated to the GlobCover2009 product, which is based on satellite observations from the Medium Resolution Imaging Spectrometer. Topographic statistics now rely on observations from the Shuttle Radar Topography Mission. Finally, vegetation height inputs are derived from space-borne Lidar measurements. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of soil moisture from the models root-zone excess reservoir into the surface excess reservoir. Specifically, the replenishment of soil moisture near the surface from below under non-equilibrium conditions was substantially reduced, which brings the models surface soil moisture more in line with the SMAP Level 2 and in situ soil moisture. Finally, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4. Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not, and on balance, the overall improvement is modest at best. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation minus forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4 SM product and the lessons learned from the effort.

Reichle, Rolf↗

DOE OSTI2026

FIREFLY: heat load and particle exhaust approximations for rapid evaluation of divertor designs

The divertor in a magnetic confinement fusion reactor is an essential component for power dissipation and particle removal. The FIREFLY package for rapid evaluation of divertor designs is presented as an extension of the FLARE code for field line reconstruction from a flux tube mesh. First, divertor loads are approximated with a simplified heat transport model. Neutralized particles are then sampled from the resulting load distribution, and the EIRENE code is used to track molecules and atoms in a plasma background while accounting for dissociation, charge exchange and ionization. Particles are removed on pumping surfaces in order to estimate the exhaust efficiency for a given divertor geometry. Optimization of the divertor geometry for more efficient particle exhaust is explored by using W7-X as an example, and the sensitivity to model parameters for the plasma background in the proxy calculations is evaluated.

mesh generation, magnetic field lines, scrape-off ↗

DOE OSTI2024

Misaligned Wind‐Waves Behind Atmospheric Cold Fronts

Abstract Atmospheric fronts embedded in extratropical cyclones are high‐impact weather phenomena, contributing significantly to mid‐latitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves under atmospheric cold fronts on air‐sea fluxes remain undocumented. Using the multi‐year in situ near‐surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hr. The parameterized effect of misaligned waves in a fully coupled model significantly increases the roughness length (185%), drag coefficient (19%), and air‐sea momentum flux (11%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the decrease in temperature and humidity scaling parameters being greater than the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave‐based bulk flux algorithm. Yet, considering this effect in the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind‐wave coupling in the bulk formula of the numerical models may help improve the air‐sea interaction simulations under the passing atmospheric fronts in the mid‐latitudes.

Oceanography↗

DOE OSTI2015

Limited-memory adaptive snapshot selection for proper orthogonal decomposition

Reduced order models are useful for accelerating simulations in many-query contexts, such as optimization, uncertainty quantification, and sensitivity analysis. However, offline training of reduced order models can have prohibitively expensive memory and floating-point operation costs in high-performance computing applications, where memory per core is limited. To overcome this limitation for proper orthogonal decomposition, we propose a novel adaptive selection method for snapshots in time that limits offline training costs by selecting snapshots according an error control mechanism similar to that found in adaptive time-stepping ordinary differential equation solvers. The error estimator used in this work is related to theory bounding the approximation error in time of proper orthogonal decomposition-based reduced order models, and memory usage is minimized by computing the singular value decomposition using a single-pass incremental algorithm. Results for a viscous Burgers’ test problem demonstrate convergence in the limit as the algorithm error tolerances go to zero; in this limit, the full order model is recovered to within discretization error. The resulting method can be used on supercomputers to generate proper orthogonal decomposition-based reduced order models, or as a subroutine within hyperreduction algorithms that require taking snapshots in time, or within greedy algorithms for sampling parameter space.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Field Demonstration of The Krauklis Seismic Wave In A Novel Mva Method For Geologic CO 2 Storage

The Energy & Environmental Research Center (EERC) conducted a feasibility study of a new, low-impact geophysics monitoring method designed to incrementally track the saturation front of CO 2 injected into a reservoir. The method leverages a new way of transmitting energy from the surface to the reservoir and employs a new subsurface signal called the Krauklis wave (K-wave) and other guided waves for seismic monitoring. The EERC partnered with Seismos, Inc., the company that holds the exclusive technology license for this method, to advance development of K-wave technology, and Denbury Onshore, LLC, to field-test the system. The system was field-tested at Denbury’s Bell Creek oil field, currently undergoing CO 2 -based enhanced oil recovery (EOR) in southeastern Montana. The K-wave system uses a wellhead-mounted source to propagate energy down the wellbore in the form of a tube wave and through the well perforations into the reservoir. Part of the energy transmitted through the perforations into the reservoir converts to slow moving K-waves and/or other guided waves that are expected to propagate through the formation to neighboring wells. The propagated signal would induce tube waves in the neighboring wells that return to the surface to be recorded with a wellhead-mounted receiver. The propagation time of the signal is directly affected by the presence of CO 2 between the wells, so the difference in transit time of the returning signal between pre- and postinjection data sets could be used to help define the lateral extent of a CO 2 saturation front. Objectives of the study were to demonstrate, validate, and evaluate the K-wave method’s ability to monitor the morphology of injected CO 2 and estimate saturation distribution. Results from the demonstration were to be compared to conventional 4-D surface seismic methods, dynamic reservoir simulations, and CO 2 breakthrough data to validate the technology and provide useful information that may result in improvements. Prior to mobilizing for the full-scale K-wave baseline survey, a field test was conducted to better understand wellhead hardware configurations and test the response from the Bell Creek reservoir using three wells. Useful operational information resulted but signal from the first-generation gas release/gas injection source was not recorded in neighboring wells. After the initial field test, numerical modeling of the K-wave method was initiated using parameters specific to the Bell Creek Field. Modeling results indicated that the fundamental concepts behind the K-wave method exhibited problems and revealed a significantly greater attenuation of the signal at the wellbore perforations than expected, totaling on the order of 10-8 from source well to receiver well over a 400-m well spacing. Modeling also showed that greater source frequency bandwidth was needed to produce guided waves within a 30-foot-thick reservoir layer such as that at Bell Creek and that reservoir layer seismic velocity would need to be slower than the velocities of the bounding layers in order to act as a wave guide opposite to what exists at Bell Creek. Modeling results prompted a second field test with a stronger, second-generation displacement source developed in the interim. Concurrent with the second field test, a 3-D seismic survey was acquired to serve as a baseline for project validation and to serve as a monitor data set for the portion of the project area where CO 2 injection had already commenced. The new source was tested using a four-well geometry and failed to produce a detectable signal, which appeared to confirm modeling results and indicated that the Bell Creek Field is not an ideal site for the K-wave test. It was concluded that a source with sufficient strength and bandwidth could not be built to overcome the factors discovered by the modeling within a reasonable amount of time. Several other sites undergoing CO 2 -based EOR were also evaluated; however, none of these sites was amenable because of geology, limitations of source and receiver hardware, or the advanced state of the CO 2 flood. In addition, two other configurations of the source and receiver layout were considered; however, these design modifications either changed the scope of the project or did not meet the goals of the Carbon Storage Program. In consultation with the U.S. Department of Energy (DOE) Project Manager and Technology Manager, a no-go decision was invoked to forego additional expenditures as the technology would not meet current project objectives. Project activities to date have moved the technology forward by developing modeling equations that can be applied elsewhere, advanced the state of wellhead sensors and sources, and spurred improvements to field operations. The 3-D seismic survey in the project area and the dynamic reservoir simulations done as part of this study have provided value for ongoing and future synergistic activities at Bell Creek related to DOE carbon storage goals.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The DECam MAGIC Survey: Spectroscopic Follow-up of the Most Metal-poor Stars in the Distant Milky Way Halo *

In this work, we present high-resolution spectroscopic observations for six metal-poor stars with [Fe/H] < –3 (including one with [Fe/H] < –4), selected using narrowband Ca ii HK photometry from the DECam MAGIC Survey. The spectroscopic data confirm the accuracy of the photometric metallicities and allow for the determination of chemical abundances for 16 elements, from carbon to barium. The program stars have chemical abundances consistent with the [Fe/H] < –3 range. A kinematic/dynamical analysis suggests that all program stars belong to the distant Milky Way halo population (heliocentric distances 35 < d helio /kpc ≲ 55), including three with high-energy orbits that might have been associated with the Magellanic system and one, J0026−5445, having parameters consistent with being a member of the Sagittarius stream. The remaining two stars show kinematics consistent with the Gaia-Sausage/Enceladus dwarf galaxy merger. J0433−5548, with [Fe/H] = –4.12, is a carbon-enhanced ultra metal-poor star, with [C/Fe] = +1.73. This star is believed to be a bona fide second-generation star, and its chemical abundance pattern was compared with yields from metal-free supernova models. Results suggest that J0433−5548 could have been formed from a gas cloud enriched by a single supernova explosion from an ∼11 M ⊙ star in the early Universe. The successful identification of such objects demonstrates the reliability of photometric metallicity estimates, which can be used for target selection and statistical studies of faint targets in the Milky Way and its satellite population. These discoveries illustrate the power of measuring chemical abundances of metal-poor Milky Way halo stars to learn more about early galaxy formation and evolution.

CEMP stars↗

DOE OSTI2019

On-line Generation and Error Handling for Surrogate Models within Multifidelity Uncertainty Quantification

Uncertainty quantification is recognized as a fundamental task to obtain predictive numerical simulations. However, many realistic engineering applications require complex and computationally expensive high-fidelity numerical simulations for the accurate characterization of the system responses. Moreover, complex physical models and extreme operative conditions can easily lead to hundreds of uncertain parameters that need to be propagated through high-fidelity codes. Under these circumstances, a single fidelity approach, i.e. a workflow that only uses high-fidelity simulations to perform the uncertainty quantification task, is unfeasible due to the prohibitive overall computational cost. In recent years, multifidelity strategies have been introduced to overcome this issue. The core idea of this family of methods is to combine simulations with varying levels of fidelity/accuracy in order to obtain the multifidelity estimators or surrogates with the same accuracy of their single fidelity counterparts at a much lower computational cost. This goal is usually accomplished by defining a prioria sequence of discretization levels or physical modeling assumptions that can be used to decrease the complexity of a numerical realization and thus its computational cost. However ,less attention has been dedicated to low-fidelity models that can be built directly from the small number of high-fidelity simulations available. In this work we focus our attention on Reduced-Order Models that can be considered a particular class of data-driven approaches. Our main goal is to explore the combination of multifidelity uncertainty quantification and reduced-order models to obtain an efficient framework for propagating uncertainties through expensive numerical codes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

The Next Generation Virgo Cluster Survey (NGVS). XXVI. The Issues of Photometric Age and Metallicity Estimates for Globular Clusters

Large samples of globular clusters (GC) with precise multi-wavelength photometry are becoming increasingly available and can be used to constrain the formation history of galaxies. We present the results of an analysis of Milky Way (MW) and Virgo core GCs based on 5 optical-near-infrared colors and 10 synthetic stellar population models. For the MW GCs, the models tend to agree on photometric ages and metallicities, with values similar to those obtained with previous studies. When used with Virgo core GCs, for which photometry is provided by the Next Generation Virgo cluster Survey (NGVS), the same models generically return younger ages. This is a consequence of the systematic differences observed between the locus occupied by Virgo core GCs and models in panchromatic color space. Only extreme fine-tuning of the adjustable parameters available to us can make the majority of the best-fit ages old. Although we cannot exclude that the formation history of the Virgo core may lead to more conspicuous populations of relatively young GCs than in other environments, we emphasize that the intrinsic properties of the Virgo GCs are likely to differ systematically from those assumed in the models. Thus, the large wavelength coverage and photometric quality of modern GC samples, such as those used here, is not by itself sufficient to better constrain the GC formation histories. Models matching the environment-dependent characteristics of GCs in multi-dimensional color space are needed to improve the situation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Rare Lepton Decays and Differentiable Hadronization Models - From Signatures of New Physics to Data-driven Event Generation

This dissertation is partitioned into two parts: phenomenological studies focused on rare lepton decays as probes of heavy and light new physics, and the development of differentiable, data-driven hadronization models. Part I develops the phenomenology of new physics signatures stemming from rare charged lepton flavor violating decays probed by experiments at the intensity frontier. These include interactions mediated by both high-scale effective operators and light new physics, manifesting in multi-lepton final states ($\mu \to 5e$), elastic nuclear transitions ($\mu \to e$ conversion), baryon-number-violating muon capture, and time-dependent signals from ultralight dark matter ($\mu \to e \phi, \tau \to \ell \phi$). Part II develops two distinct strategies for advancing differentiable and data-driven hadronization models. One involves comprehensive reweighting frameworks for hadronization that enable efficient uncertainty estimation, facilitate parameter tuning, and interface naturally with differentiable programming paradigms. The other introduces machine-learning-based methods for extracting microscopic fragmentation dynamics directly from macroscopic observables through the deformation of existing models -- effectively providing solutions to the inverse problem of hadronization. Altogether, these studies advance the interpretability, flexibility, and precision of theoretical predictions for both high-intensity and high-energy experiments.

Menzo, Tony [Cincinnati U.] (ORCID:000000022013457↗