Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “input 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.

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

At least 1,045 records · Page 58

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↗

3D orthorhombic earth model effects on seismic source characterization

Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).

58 GEOSCIENCES↗

235-F GoldSim Fate and Transport Model: Uncertainty Quantification

Building 235-F was configured with two missions in mind: Actinide Billet Line (ABL) and the fabrication of Pu-238 oxide for space program applications. ABL produced Np-237 billets for use in SRS reactors, whereas the design process, fabrication, and examination of Pu-238 oxide powder occurred within the following areas, respectively: Plutonium Experimental Facility (PEF), Plutonium Fuel Form (PuFF), and Old Metallography Lab (OML). By 1990 production ceased and by 2006 de-inventory occurred; however, assays have shown significant holdup remains within ABL and PuFF. As a result, 235-F is a Category 2 nuclear facility, with plans to undergo deactivation and decommission (D and D) via In-Situ Disposal (ISD). The purpose of this project is to ensure United States Environmental Protection Agency (USEPA) groundwater radiation maximum contaminant level (MCL) and dosage standards are met during the D and D of 235-F by quantifying uncertainty through probabilistic modeling and evaluation of various ISD alternatives. GoldSim is a dynamic modeling software package with a graphical, object-oriented interface capable of capturing the influence of complex system input variability on probabilistic system outcomes. A GoldSim stochastic fate and transport model for 235-F was developed and matched with a PORFLOW deterministic model to simulate probabilistic release and flow of radionuclides from ABL and PuFF into the vadose zone, the Upper Three Runs (UTR) Aquifer, and UTR Creek. The GoldSim model was used to probabilistically evaluate four ISD scenarios against USEPA groundwater radiation MCLs and dosage standards. The deterministic 235-F GoldSim fate and transport model continues to be refined to match the results of the PORFLOW deterministic model to ensure the model accurately represents radionuclide movement through the groundwater system. The stochastic variables that are utilized within the GoldSim model are founded on the most current data; a conservative perspective is taken where needed. Alignment with the PORFLOW deterministic model, coupled with input stochastic variability, allows the probabilistic 235-F GoldSim model to capture the conservative breadth of possible outcomes for radionuclide fate within this particular system. This ensures that the USEPA MCLs and dosage limits hold even in the worst case scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Mesoscale modeling of microstructure-dependent thermal conductivity in U-Zr fuels

In uranium-zirconium (U-Zr) based metallic fuels, different phases can form at different compositions and temperatures. Typically, lamellar δ-UZr 2 and α-U phases are the dominant microstructures in U-rich U-Zr alloys at temperatures below 880 K. In this work, a finite element method based mesoscale modeling technique is used to calculate the effective thermal conductivities of such heterogeneous microstructures, using the thermal conductivities of two individual phases and their interphase thermal resistance (Kapitza resistance) as input parameters. The Kapitza resistance between δ-UZr 2 and α-U is determined at different temperatures, which shows an approximately T 3 dependence in the temperature range between 300 and 800 K. In addition, the Kapitza resistance exhibits a strong dependence on the aspect ratio of the δ-UZr 2 phase. Further, an analytical model is therefore developed to quantify the effects of both temperature and δ-UZr 2 aspect ratio on the Kapitza resistance. Using this newly developed Kapitza resistance model, the effective thermal conductivities of a number of δ-UZr 2 + α-U heterogeneous microstructures in U-Zr alloys, including non-lamellar microstructures, can be estimated accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel’dovich measurements from BOSS CMASS and LOWZ halos

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel’dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope DR5 and Planck in combination with the CMASS (mean redshift $\langle z \rangle$ = 0.55 and host halo mass $\langle M_{vir} \rangle$ = 3 × 10 13 M ⊙ ) and LOWZ ( $\langle z \rangle$ = 0.31 , $\langle M_{vir} \rangle$ = 5 × 10 13 M ⊙ ) galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothesis at 6.5σ, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally > 90 σ for the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e., the electron thermal pressure profile of the same CMASS objects, and reject the no-tSZ hypothesis at 10σ. We combine tSZ and kSZ measurements to estimate the electron temperature to 20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ) signal at 2.9 (13.9) σ , and leave their interpretation to future work. This paper and the companion paper demonstrate that current CMB experiments can detect and resolve gas profiles in low mass halos and at high redshifts, which are the most sensitive to feedback in galaxy formation and the most difficult to measure any other way. They will be a crucial input to cosmological hydrodynamical simulations, thus improving our understanding of galaxy formation. These precise gas profiles are already sufficient to reduce the main limiting theoretical systematic in galaxy-galaxy lensing: baryonic uncertainties. Future such measurements will thus unleash the statistical power of weak lensing from the Rubin, Euclid and Roman observatories. Finally, our stacking software ThumbStackis publicly available and directly applicable to future Simons Observatory and CMB-S4 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Soil phosphorus drawdown by perennial bioenergy cropping systems in the Midwestern US

Abstract Without fertilization, harvest of perennial bioenergy cropping systems diminishes soil nutrient stocks, yet the time course of nutrient drawdown has not often been investigated. We analyzed phosphorus (P) inputs (fertilization and atmospheric deposition) and outputs (harvest and leaching losses) over 7 years in three representative biomass crops—switchgrass ( Panicum virgatum L.), miscanthus ( Miscanthus × giganteus ) and hybrid poplar trees ( Populus nigra × P. maximowiczii )—as well as in no‐till corn (maize; Zea mays L.) for comparison, all planted on former cropland in SW Michigan, USA. Only corn received P fertilizer. Corn (grain and stover), switchgrass, and miscanthus were harvested annually, while poplar was harvested after 6 years. Soil test P (STP; Bray‐1 method) was measured in the upper 25 cm of soil annually. Harvest P removal was calculated from tissue P concentration and harvest yield (or annual woody biomass accrual in poplar). Leaching was estimated as total dissolved P concentration in soil solutions sampled beneath the rooting depth (1.25 m), combined with hydrological modeling. Fertilization and harvest were by far the dominant P budget terms for corn, and harvest P removal dominated the P budgets in switchgrass, miscanthus, and poplar, while atmospheric deposition and leaching losses were comparatively insignificant. Because of significant P removal by harvest, the P balances of switchgrass, miscanthus, and poplar were negative and corresponded with decreasing STP, whereas P fertilization compensated for the harvest P removal in corn, resulting in a positive P balance. Results indicate that perennial crop harvest without P fertilization removed legacy P from soils, and continued harvest will soon draw P down to limiting levels, even in soils once heavily P‐fertilized. Widespread cultivation of bioenergy crops may, therefore, alter P balances in agricultural landscapes, eventually requiring P fertilization, which could be supplied by P recovery from harvested biomass.

09 BIOMASS FUELS↗

County Level Annual Population Projections for SSP3 and SSP5

This dataset consist of annual (2020-2100) county-level population projections for the United States (U.S.) for Shared Socioeconomic Pathways (SSPs) 3 and 5. The original decadal state-level data is used as an input to the gridded population data model to downscale the state-level projections for each SSP to a 1 km grid. The 1 km data is then aggregated to the county-level using the 2020 TIGER/Line county shapefiles from the U.S. Census Bureau. The two data files are shared as .csv files with the following structure: Rows = Data for individual counties which are identified using their Federal Information Processing Standard (FIPS) code. The FIPS code is stored in the first column. Columns = Each year from 2020-2100 is an individual column. For easier reference, the state name associated with each row is stored in the last column. Values in each cell are the downscaled estimated population for that county-year combination. Values are fractional due to the weighting scheme used in the downscaling. The paper and model source code cited in the "Related Works" section describes the initial source of the population projections.

FIPS↗

A Modified Vegetation Photosynthesis and Respiration Model (VPRM) for the Eastern USA and Canada, Evaluated With Comparison to Atmospheric Observations and Other Biospheric Models

Atmospheric CO 2 measurements from a dense surface network can help to evaluate terrestrial biosphere model (TBM) simulations of Net Ecosystem Exchange (NEE) with two key benefits. First, gridded CO 2 flux estimates can be evaluated over regional scales, not possible using flux tower observations at discrete locations for model evaluation. Second, TBM ability to explain atmospheric CO 2 fluctuations due to the biosphere can be directly tested, an important objective for anthropogenic emissions monitoring using atmospheric observations. Here, we customize the Vegetation Photosynthesis and Respiration Model (VPRM) for an eastern North American domain with strong biological activity upwind of urban areas. Parameters are optimized using flux tower observations from a historical database with sites in (and near) the domain. In addition, the respiration model (originally a linear function of temperature) is modified to account for impacts of changing foliage, non-linear temperature, and water stress. Flux estimates from VPRM, the Carnegie-Ames-Stanford Approach (CASA) model and the Simple Biosphere Model v4 (SiB4), are convolved with footprints from atmospheric transport models for evaluation with CO 2 observations at 21 towers in the domain, with roughly half of the towers used here for the first time. Results show that the new respiration model in VPRM helps to correct a growing season sink bias in the atmosphere associated with underestimated summertime respiration using the original model with annual parameters. The new VPRM also better explains fine-scale atmospheric CO 2 variability compared to other TBMs, due to higher resolution diagnostic phenology, the new respiration model, domain-specific parameters, and high-quality input data sets.

54 ENVIRONMENTAL SCIENCES↗

Statistical Downscaling of Climate Models for Solar Resource Assessment

This study presents the development of statistical models to efficiently downscale future projections of solar irradiance for solar energy applications. A climate data set simulated from a Regional Climate Model (RCM) obtained from the North American Coordinated Regional Climate Downscaling Experiment (NA-CORDEX) is selected as input to the statistical models to create high-resolution global horizontal irradiance (GHI) over the contiguous United States (CONUS). Our approach builds statistical downscaling models that (1) regrid RCM data (0.22 degree and daily spatiotemporal resolution), (2) correct bias of GHI projections, (3) downscale the future GHI project from daily-scale to hourly-scale, and (4) spatially downscale to generate GHI at 8-km resolution. To calibrate and validate the statistical models, we adapt and use the National Solar Radiation Database (NSRDB). Preliminary results show that the statistical downscaling approach downscales future projections of GHI under two climate scenarios (RCP4.5 and RCP8.5) with a nBIAS of 3%, nMAE of 34% and nRMSE of 46% estimated against NSRDB for the contiguous United State. This presentation will summarize the implemented methodology and validation results as well as future extension of this research.

climate data↗

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

Rapidly changing high-latitude seasonality: implications for the 21st century carbon cycle in Alaska

Seasonal variations in high-latitude terrestrial carbon (C) fluxes are predominantly driven by air temperature and radiation. At present, high-latitude net C uptake is largest during the summer. Recent observations and modeling studies have demonstrated that ongoing and projected climate change will increase plant productivity, microbial respiration, and growing season lengths at high-latitudes, but impacts on high-latitude C cycle seasonality (and potential feedbacks to the climate system) remain uncertain. Here we use ecosys, a well-tested and process-rich mechanistic ecosystem model that we evaluate further in this study, to explore how climate warming under an RCP8.5 scenario will shift C cycle seasonality in Alaska throughout the 21st century. The model successfully reproduced recently reported large high-latitude C losses during the fall and winter and yet still predicts a high-latitude C sink, pointing to a resolution of the current conflict between process-model and observation-based estimates of high-latitude C balance. We find that warming will result in surprisingly large changes in net ecosystem exchange (NEE; defined as negative for uptake) seasonality, with spring net C uptake overtaking summer net C uptake by year 2100. This shift is driven by a factor of 3 relaxation of spring temperature limitation to plant productivity that results in earlier C uptake and a corresponding increase in magnitude of spring NEE from –19 to –144 gC m –2 season –1 by the end of the century. Although a similar relaxation of temperature limitation will occur in the fall, radiation limitation during those months will limit increases in C fixation. Additionally, warmer soil temperatures and increased carbon inputs from plants lead to combined fall and winter C losses (163 gC m –2 ) that are larger than summer net uptake (123 gC m –2 season –1 ) by year 2100. However, this increase in microbial activity leads to more rapid N cycling and increased plant N uptake during the fall and winter months that supports large increases in spring NPP. Due to the large increases in spring net C uptake, the high-latitude atmospheric C sink is projected to sustain throughout this century. Our analysis disentangles the effects of key environmental drivers of high-latitude seasonal C balances as climate changes over the 21st century.

54 ENVIRONMENTAL SCIENCES↗

The LSST DESC data challenge 1: generation and analysis of synthetic images for next-generation surveys

Data Challenge 1 (DC1) is the first synthetic data set produced by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). DC1 is designed to develop and validate data reduction and analysis and to study the impact of systematic effects that will affect the LSST data set. DC1 is comprised of r -band observations of 40 deg 2 to 10 yr LSST depth. In this paper, we present each stage of the simulation and analysis process: (a) generation, by synthesizing sources from cosmological N -body simulations in individual sensor-visit images with different observing conditions; (b) reduction using a development version of the LSST Science Pipelines; and (c) matching to the input cosmological catalogue for validation and testing. We verify that testable LSST requirements pass within the fidelity of DC1. We establish a selection procedure that produces a sufficiently clean extragalactic sample for clustering analyses and we discuss residual sample contamination, including contributions from inefficiency in star–galaxy separation and imperfect deblending. We compute the galaxy power spectrum on the simulated field and conclude that: (i) survey properties have an impact of 50 per cent of the statistical uncertainty for the scales and models used in DC1; (ii) a selection to eliminate artefacts in the catalogues is necessary to avoid biases in the measured clustering; and (iii) the presence of bright objects has a significant impact (2σ–6σ) in the estimated power spectra at small scales (ℓ > 1200), highlighting the impact of blending in studies at small angular scales in LSST.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗

Design and analysis of a wake steering controller with wind direction variability

Wind farm control strategies are being developed to mitigate wake losses in wind farms, increasing energy production. Wake steering is a type of wind farm control in which a wind turbine's yaw position is misaligned from the wind direction, causing its wake to deflect away from downstream turbines. Current modeling tools used to optimize and estimate energy gains from wake steering are designed to represent wakes for fixed wind directions. However, wake steering controllers must operate in dynamic wind conditions and a turbine's yaw position cannot perfectly track changing wind directions. Research has been conducted on robust wake steering control optimized for variable wind directions. In this paper, the design and analysis of a wake steering controller with wind direction variability is presented for a two-turbine array using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented wake model. First, the authors propose a method for modeling the turbulent and low-frequency components of the wind direction, where the slowly varying wind direction serves as the relevant input to the wake model. Next, we explain a procedure for finding optimal yaw offsets for dynamic wind conditions considering both wind direction and yaw position uncertainty. We then performed simulations with the optimal yaw offsets applied using a realistic yaw offset controller in conjunction with a baseline yaw controller, showing good agreement with the predicted energy gain using the probabilistic model. Using the Gaussian wake model in FLORIS as an example, we compared the performance of yaw offset controllers optimized for static and dynamic wind conditions for different turbine spacings and turbulence intensity values, assuming uniformly distributed wind directions. For a spacing of five rotor diameters and a turbulence intensity of 10 %, robust yaw offsets optimized for variable wind directions yielded an energy gain equivalent to 3.24 % of wake losses recovered, compared to 1.42 % of wake losses recovered with yaw offsets optimized for static wind directions. In general, accounting for wind direction variability in the yaw offset optimization process was found to improve energy production more as the separation distance increased, whereas the relative improvement remained roughly the same for the range of turbulence intensity values considered.

17 WIND ENERGY↗

Using angular two-point correlations to self-calibrate the photometric redshift distributions of DECaLS DR9

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-z) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for χ 2 minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests, and find reasonable agreements: (1) finer photo-z bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. In conclusion, this work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Environmentally Assisted Fatigue in Light Water Reactor Environment

This report summarizes the Environmentally Assisted Fatigue (EAF) research conducted at ANL under the US DOE Light Water Reactor Sustainability (LWRS) program. Starting from a rich background in theoretical and experimental EAF, ANL previously developed an approach to evaluate fatigue performance of reactor materials in light water reactor environments with the correction factor F en . The approach was based on a large body of experimental work performed at ANL and elsewhere, and was consistent with American Society of Mechanical Engineers (ASME)’s methodology governing the design and construction of reactor components. In recent years, the program was focused on component fatigue prediction and made several major and fundamental contributions in this area. These accomplishments help meet the needs identified by the industry concerning component level fatigue predictions in complex, transient conditions. The main contribution of the ANL program involved the development of a system-level model for estimating residual strain and life of nuclear reactor coolant system components under connected-system-thermal-mechanical boundary conditions. The goal was to predict the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the finite element (FE) model results, the strain residuals, strain amplitudes and resulting fatigue lives of reactor coolant system (RCS) components were predicted. The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. In addition, the simulated component-level strain profile can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests. Building upon the system-level model, ANL developed a digital twin (DT) framework to predict the structural states and associated fatigue life of components in real-time. This framework is a comprehensive system designed to predict the structural states and fatigue lives of reactor components. It includes multiple models and integrates artificial intelligence (AI), machine learning (ML), and FE based modeling tools to evaluate the structural states and fatigue lives.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models

Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.

97 MATHEMATICS AND COMPUTING↗

Aggregation and data driven identification of building thermal dynamic model and unmeasured disturbance

An aggregate model is a single-zone equivalent of a multi-zone building, and is useful for many purposes, including model based control of large heating, ventilation and air conditioning (HVAC) equipment. This paper deals with the problem of simultaneously identifying an aggregate thermal dynamic model and unknown disturbances from input–output data of multi-zone buildings. The unknown disturbance is a key challenge since it is not measurable but non-negligible. In this paper, we first present a principled method to aggregate a multi-zone building model into a single zone model, and show the aggregation is not as trivial as it has been assumed in the prior art. We then provide a method to identify the parameters of the model and the unknown disturbance for this aggregate (single-zone) model. Finally, we test our proposed identification algorithm to data collected from a multi-zone building testbed in Oak Ridge National Laboratory. A key insight provided by the aggregation method allows us to recognize under what conditions the estimation of the disturbance signal will be necessarily poor and uncertain, even in the case of a specially designed test in which the disturbances affecting each zone are known (as the case of our experimental testbed). This insight is used to provide a heuristic that can be used to assess when the identification results are likely to have high or low accuracy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗