Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive state 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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Knowledge-informed deep learning for hydrological model calibration: an application to Coal Creek Watershed in Colorado

Abstract. Deep learning (DL)-assisted inverse mapping has shown promise in hydrological model calibration by directly estimating parameters from observations. However, the increasing computational demand for running the state-of-the-art hydrological model limits sufficient ensemble runs for its calibration. In this work, we present a novel knowledge-informed deep learning method that can efficiently conduct the calibration using a few hundred realizations. The method involves two steps. First, we determine decisive model parameters from a complete parameter set based on the mutual information (MI) between model responses and each parameter computed by a limited number of realizations (∼50). Second, we perform more ensemble runs (e.g., several hundred) to generate the training sets for the inverse mapping, which selects informative model responses for estimating each parameter using MI-based parameter sensitivity. We applied this new DL-based method to calibrate a process-based integrated hydrological model, the Advanced Terrestrial Simulator (ATS), at Coal Creek Watershed, CO. The calibration is performed against observed stream discharge (Q) and remotely sensed evapotranspiration (ET) from the water year 2017 to 2019. Preliminary MI analysis on 50 realizations resulted in a down-selection of 7 out of 14 ATS model parameters. Then, we performed a complete MI analysis on 396 realizations and constructed the inverse mapping from informative responses to each of the selected parameters using a deep neural network. Compared with calibration using observations covering all time steps, the new inverse mapping improves parameter estimations, thus enhancing the performance of ATS forward model runs. The Nash–Sutcliffe efficiency (NSE) of streamflow predictions increases from 0.53 to 0.8 when calibrating against Q alone. Using ET observations, on the other hand, does not show much improvement on the performance of ATS modeling mainly due to both the uncertainty of the remotely sensed product and the insufficient coverage of the model ET ensemble in capturing the observation. By using observed Q only, we further performed a multiyear analysis and show that Q is best simulated (NSE > 0.8) by including in the calibration the dry-year flow dynamics that show more sensitivity to subsurface characteristics than the other wet years. Moreover, when continuing the forward runs till the end of 2021, the calibrated models show similar simulation performances during this evaluation period as the calibration period, demonstrating the ability of the estimated parameters in capturing climate sensitivity. Our success highlights the importance of leveraging data-driven knowledge in DL-assisted hydrological model calibration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Technoeconomic and Life Cycle Analysis of an Integrated Fermentation and Microbial Electrochemical Process for Volatile Fatty Acid Production from Food Waste

Techno-economic analysis (TEA) and life cycle assessment (LCA) were conducted for an integrated system designed for the production of volatile fatty acid (VFA) from food waste. The TEA estimated a production cost of $\$$3.12/kg VFA, and the LCA predicted negative greenhouse gas (GHG) emissions of -0.4 kg CO 2 e/kg VFA, driven primarily by diverting organic waste from landfills and avoiding methane emissions while producing valuable chemical products. Hotspot analysis showed arrested methanogenesis (AM) fermentation as the largest contributor to costs (37%) and environmental burden (47%), driven by high sodium hydroxide (NaOH) consumption. Distillation and microbial electrosynthesis (MES) units were the next-largest environmental contributors (28% and 18%). Major cost drivers also included residuals management (biosolids and wastewater) and the equipment and operating costs for AM, MES, and sonication pretreatment units. Although the new integrated system is environmentally benign, its costs and environmental impacts can be further reduced by integrating alternative energy sources, minimizing chemical and energy inputs through process optimization, and improving efficiency. In conclusion, this work highlighted the viability of waste-derived VFA production and provided a clear, data-driven strategy to accelerate the commercialization of waste valorization technology.

Carboxylic Acid Production↗
DOE OSTI2024

Aboveground woody biomass estimation of young bioenergy plantations of Populus and its hybrids using mobile (backpack) LiDAR remote sensing

Woody aboveground biomass (AGB) including short-rotation Populus is used as a feedstock for renewable and carbon-neutral bioenergy. While woody AGB can be estimated with allometric equations requiring labor-intensive field data, remote sensing technologies like mobile terrestrial light detection and ranging (LiDAR) can estimate woody AGB quickly and accurately. Therefore, the goals of this study were to develop a model to predict woody AGB of 2-year-old Populus spp. from three taxa (P. deltoides, P. deltoides × P. maximowiczii and P. deltoides × P. trichocarpa) using allometric (height and diameter at breast height (DBH)) or LiDAR-derived metrics from a mobile terrestrial (backpack) system. Likewise, we sought to compare LiDAR-estimated tree height and DBH with field-measured values. We found that a taxa-specific model containing LiDAR-measured tree height, crown volume, and taxa interactions with the height of the 10 th percentile, and the density of the lowest interval (density metric 0) explained 84 % of the variation in woody AGB with a root mean square error (RMSE) of 28.7 % and performed slightly better than the allometric model. The best model excluding taxa had a slightly higher RMSE but lower bias than the allometric model. LiDAR-derived tree heights were highly correlated with field-measured heights, but DBH could not be estimated accurately. Therefore, terrestrial mobile LiDAR systems can accurately estimate woody AGB and tree height of Populus in short rotation systems to aid in the fast and efficient quantification of woody bioenergy production and renewable energy resources.

AGB↗
DOE OSTI2026

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2021

Characterization of shale using Helium and Argon at high pressures

In order to estimate the shale gas in place and the eventual recovery during shale gas operations, one of the key requirements is to accurately characterize the shale’s petrophysical and transport properties such as porosity, permeability, diffusivity, and storage capacity. Despite the many efforts reported in the technical literature aiming to provide an improved understanding of the complex pore structures and the associated fluid flow in gas shales, complete characterization of organic-rich shale samples still poses a big challenge. Here, we have characterized mass transfer and sorption in shale at different length scales using Helium (He) and Argon (Ar) as probe gases. Thermogravimetric analysis (TGA) with a shale cube of ~1 cm 3 in size and gas expansion experiments with a full-diameter core (3.5” in diameter) were used to measure sorption kinetics/isotherms and mass transfer, respectively. Both samples are from the same depth/location in the Marcellus shale formation. The TGA steady-state technique was initially used to generate excess sorption isotherms for Ar, while dynamic TGA experiments were used to study its sorption kinetics. The TGA experiments demonstrate that Ar, which has a similar sorption potential as Methane, but is generally assumed to be inert, adsorbs onto the surfaces of the mesoporous and microporous regions of the shale samples according to a Langmuir-type behavior. Helium expansion experiments, on the full-diameter core, were used to measure the overall porosity, on the basis that He is a non-sorbing and inert gas as compared to Ar. The He expansion experiments, furthermore, allow us to delineate the mass transfer of gas across the inherent hierarchy of pore sizes, including macropores (macro- and microcracks), mesopores and micropores. Similar expansion experiments were also performed with Ar to study the combined impact of mass transfer and sorption. A triple-porosity model (TPM) was utilized to interpret the He expansion experiments with the shale core and to extract (estimate) relevant transport parameters. We report and compare here the diffusivities and permeabilities of the whole core for both He and Ar, as calculated from the modeling and fitting of the experimental data. On the premise that the shale cube is representative of the matrix region of the core, the Ar sorption kinetics from the cube experiments were subsequently combined with the transport parameters extracted from the He experiments to predict the behavior of the Ar expansion test with the full-diameter core. An excellent agreement is observed between the model predictions and the experimental data. The experimental observations and their interpretation indicate that one must be cautious when using Ar to estimate the true porosity and permeability of shales. In addition, we demonstrate that He and Ar probe gases, when used in tandem, can be employed effectively as a tool to characterize shales in terms of mass transfer and sorption dynamics across scales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

Verification of Reactor Activation Modelling with Physical Characterization Data and the Impact on Long Term Assessments - 20303

Whiteshell Laboratories, located in Manitoba, Canada, provided research facilities for the Canadian nuclear industry since the early 1960's. The facility was centered on Whiteshell Reactor 1 (WR-1), an organically cooled, heavy water moderated nuclear reactor. WR-1 was safely shut down in 1985, defueled and drained, and has been in a safe Storage with Surveillance (Deferred Decommissioning) state until present day. In support of an Environmental Assessment for proposed in-situ disposal of WR-1, detailed characterization information regarding the reactor core components (calandria, calandria tubes, and pressure tubes) was necessary to provide confidence that the long-term impacts of the facility can be assessed adequately. A simplified estimate of the total activation of core components was completed in 1992 using conservative assumptions, expected reactor configuration and computer models built with WIMS-CRNL, ONEDANT and ORIGEN-S. Due to the simplification, and age of the work, there was considerable doubt that the modelled inventory estimate of the core components was sufficient to support long-term performance modeling of the planned in-situ decommissioning. Additional data would be required to provide this confidence. Several options were considered for characterizing the WR-1 core components, but it was decided to collect only a limited number of samples in an effort to validate the 1992 simplified estimate as bounding. The limited number of samples reduced risks to workers, both radiological and non-radiological, and significantly reduced the costs required to confidently define and bound the radiological inventory of activation products in the core. The results of the sampling were compared to the simplified activation estimation. There was good agreement between the modelled and sampled results, with all major contributors to total activity showing similar relative quantities. The samples were generally lower than predicted by the models, with few exceptions. Notable exceptions included Nb-94, which were higher in the sample results than the model. This was determined to be caused by the use of Niobium containing steel alloys in the stainless steel pressure tube that were not accounted for in the models. The comparison of results provided confidence that the original inventory estimates produced through the simplified model are conservative and bounding, and therefore suitable for use in the long-term assessment of in situ disposal of the WR-1 reactor. This paper provides details of the simplified computer model results, their comparison to results of the sampling, and the changes adopted in the assessment of the in-situ decommissioning of WR-1 as a result. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Measurements of the electron energy distribution function in partially magnetized low temperature plasmas

While Langmuir probes (LPs) are relatively simple and inexpensive plasma diagnostics for the electron density, temperature, and the electron energy distribution function (EEDF), the interpretation of the measured current–voltage (I–V) characteristic is complicated considerably by the presence of a magnetic field. In regimes where the electron gyroradius is comparable to the probe radius, the electron flux to the probe surface is retarded by reduced mobility across field lines and can no longer be described by a thermal model. Predicting the current collected by the probe in these regimes requires accurate estimates of the plasma diffusion coefficients, which are usually difficult to obtain. In this work, we measure electron energy distribution functions in E×B plasmas with magnetized electrons and non-magnetized ions in argon and krypton gases, using both a LP and laser Thomson scattering (LTS) at various magnetic fields. Using the LTS measurements to provide a robust benchmark for comparison, we compare existing theories describing the flux to the probe under magnetized conditions. We find that even when the electron gyroradius associated with the effective electron temperature is small compared to the probe radius, the EEDF computed using classical probe theory is still robust at energies higher than the energy at which the gyroradius becomes larger than the probe size. For plasmas that are not strongly non-Maxwellian, we formulate a method to extract robust density and temperature measurements using physics-informed fitting techniques to analyze EEDFs computed using classical theory.

Devin, E. G. [Princeton Plasma Physics Laboratory ↗
DOE OSTI2026

A dendritic strontium river isoscape for fisheries applications in the Sacramento River basin, California, USA

Objective Understanding the origins and movements of fish is fundamental to effective conservation and fisheries management. Strontium isotope ratios ( 87 Sr/ 86 Sr) in otoliths provide a powerful tracer of natal origin and migratory pathways. However, existing 87 Sr/ 86 Sr isoscapes for the Sacramento River basin, an ecosystem that supports ecologically and economically important salmon populations, rely on discrete classification approaches that overlook unsampled habitats and do not incorporate spatial uncertainty. Our objective was to develop a continuous, network-explicit 87 Sr/ 86 Sr isoscape with quantified uncertainty to fill in data gaps and enable probabilistic assignments of fish origin and movement. Methods We used river water 87 Sr/ 86 Sr data from 106 sites (1997–2021) to develop spatial stream network models that use dendritic connectivity and watershed characteristics (lithology, bedrock age, and land cover) to predict river water 87 Sr/ 86 Sr throughout the basin. Models were fitted using maximum and restricted likelihood and were evaluated via Akaike’s information criterion and leave-one-out cross validation. We produced both historical (pre-dam) and present-day (below-dam) isoscapes, delineated uncertainty-informed isotopic ranges using k -means clustering, and applied a proof-of-concept Bayesian assignment to estimate natal origins and early rearing habitats for two endangered winter-run Chinook Salmon Oncorhynchus tshawytscha. Results Cross validation indicated strong performance of the 87 Sr/ 86 Sr model (leave-one-out cross validation: R 2 = 0.91; root mean square error = 0.0005). Uncertainty-informed clustering identified 19 isotopic “suites” (reaches with indistinguishable 87 Sr/ 86 Sr values) in present-day anadromous habitats and 25 suites in the historical network. Example natal and early rearing assignments included predictions that challenged expectations for juvenile salmon migration based on predicted river 87 Sr/ 86 Sr compositions. Conclusions This study developed a continuous, network-explicit 87 Sr/ 86 Sr isoscape that integrates existing river data to predict 87 Sr/ 86 Sr in unsampled reaches and the likely achievable range and resolution of otolith-based origin and life history inference. The resulting river isoscape provides a valuable tool to predict salmon movements and identify habitats supporting their survival and growth that otherwise might remain undetected. Coupling these predictions with complementary approaches that ground-truth juvenile presence (e.g., targeted fish surveys) represents an important step toward science-informed restoration and management of critical habitats throughout the Sacramento River basin.

Environmental sciences↗
DOE OSTI2022

Machine learning-based ethylene concentration estimation, real-time optimization and feedback control of an experimental electrochemical reactor

With the increase in electricity supply from clean energy sources, electrochemical reduction of carbon dioxide (CO 2 ) has received increasing attention as an alternative source of carbon-based fuels. As CO 2 reduction is becoming a stronger alternative for the clean production of chemicals, the need to model, optimize and control the electrochemical reduction of the CO 2 process becomes inevitable. However, on one hand, a first-principles model to represent the electrochemical CO 2 reduction has not been fully developed yet because of the complexity of its reaction mechanism, which makes it challenging to define a precise state-space model for the control system. On the other hand, the unavailability of efficient concentration measurement sensors continues to challenge our ability to develop feedback control systems. Gas chromatography (GC) is the most common equipment for monitoring the gas product composition, but it requires a period of time to analyze the sample, which means that GC can provide only delayed measurements during the operation. Moreover, the electrochemical CO 2 reduction process is catalyzed by a fast-deactivating copper catalyst and undergoes a selectivity shift from the product-of-interest at the later stages of experiments, which can pose a challenge for conventional control methods. To this end, machine learning (ML) techniques provide a potential approach to overcome those difficulties due to their demonstrated ability to capture the dynamic behavior of a chemical process from data. Motivated by the above considerations, we propose a machine learning-based modeling methodology that integrates support vector regression and first-principles modeling to capture the dynamic behavior of an experimental electrochemical reactor; this model, together with limited gas chromatography measurements, is employed to predict the evolution of gas-phase ethylene concentration. The model prediction is directly used in a proportional-integral (PI) controller that manipulates the applied potential to regulate the gas-phase ethylene concentration at energy-optimal set-point values computed by a real-time process optimizer (RTO). Specifically, the RTO calculates the operation set-point by solving an optimization problem to maximize the economic benefit of the reactor. Finally, suitable compensation methods are introduced to further account for the experimental uncertainties and handle catalyst deactivation. The proposed modeling, optimization, and control approaches are the first demonstration of active control for a CO 2 electrolyzer and contribute to the automation and scale-up efforts for electrified manufacturing of fuels and chemicals starting from CO 2 .

42 ENGINEERING↗
DOE OSTI2025

Assessing the Role of Hydrodynamics in Enhancing Height-Above-the-Nearest-Drainage Derived Synthetic Rating Curves: A Comparative Study in the Wu River Basin, Taiwan

The conventional approach to generating synthetic rating curves (SRC) using the Height-Above-the-Nearest-Drainage (HAND) method typically relies on the assumption of uniform flow, such as Manning's equation, to establish stage-discharge ratings. The zero-physics application of the uniform flow equation is insufficient for capturing detailed hydraulic features (e.g., backwater effect) and neglects the hydraulic effects from adjacent channels. This lack of hydrodynamic computation can impact the accuracy and effectiveness of riverine flood risk estimation and management. To reduce this foreseeable error, we introduce the HAND-hd workflow, which integrates sophisticated hydrodynamic computations in the production of HAND-based SRC with hydrodynamic features (SRC hd ). The results indicate that SRC hd demonstrates consistent agreement with both gauge observations and benchmark solutions. Additionally, the comparative analysis suggests that SRC hd provides notable improvements in stage-discharge ratings over conventional HAND-based SRCs, particularly in channels with mild bed gradients, where it reduces water stage prediction errors and percent biases. In steeper channel segments, SRC hd maintains comparable accuracy to conventional methods. The comprehensive evaluation in this study emphasizes the potential discrepancies and inaccuracies associated with the adoption of the uniform flow assumption in the conventional HAND-SRCs and addresses the necessity of including hydrodynamic physics in the application of HAND-based SRC (e.g., inundation map) in channels with mild gradients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Progress Towards a Predictive Eagle Behavior and Risk Modeling Framework: Overview and Recent Validation Efforts

This presentation summarizes progress to date of the U.S. Department of Energy project, "Development of a computational framework for modeling golden eagles (Aquila chrysaetos) near wind farms," which focuses on stochastic behavioral modeling of soaring raptors across landscape, facility, and turbine spatiotemporal scales. This publicly available, open-source modeling framework includes behavioral models based on three different underlying principles: energy minimization at landscape scale, behavioral heuristics at landscape-facility scale, and data-driven behaviors at the facility-micro-scale. We will briefly overview the key advancements in the behavioral modeling state of the art, which leverages multiple high-resolution telemetry data sources combined with high-fidelity atmospheric flow modeling insights. We then present preliminary results from a validation study in Altamont, California. This new study involves a novel application of the Stochastic Soaring Raptor Simulator (SSRS), in a new geographic locale, to understand facility scale eagle movement patterns over time scales representative of a wind project's lifetime. For this desktop analysis (that does not depend on any high-performance computing resources), SSRS simultaneously considers a variety of wind conditions and eagle approach vectors toward a project site of interest. This work demonstrates the integration of publicly available landscape-scale atmospheric datasets, our recently improved engineering updraft models (see presentation from Thedin et al.), and our energy minimization behavioral models within the SSRS framework. While we only present results from a single behavioral model, the integration of these three modeling components forms the foundation for our more sophisticated behavioral models (see presentations from Brandes et al., Sandhu et al.) that are under active development. Results are presented in the form of presence maps, which may be applied to estimate risk to wildlife, augment ground survey data, inform wind-plant operations, or incorporated into wind-plant designs.

agent-based modeling↗
DOE OSTI2026

The Extended Baryonic Tully–Fisher Relation for SDSS MaNGA Galaxies

The baryonic Tully–Fisher relation (BTFR), a relationship between the rotational velocity and baryonic mass in spiral galaxies, probes the relative content of baryonic and total mass in galaxies and thus provides a good test of dark matter content in galaxies. Using Hα kinematics, we model the rotation curves of the Sloan Digital Sky Survey MaNGA DR17 spiral galaxies. To extend the BTFR to higher masses with elliptical galaxies, we estimate their total masses from their stellar velocity dispersions using the virial theorem and define the effective rotational velocity as the velocity a rotation-supported galaxy would exhibit given this mass. The baryonic mass of spiral galaxies is composed of stellar, H I , H 2 , and He mass, while only the stellar mass is used for the baryonic content of ellipticals. We construct joint BTFRs for 5743 MaNGA spiral and elliptical galaxies, TNG100 simulated galaxies with baryonic masses greater than 10 9 M ⊙ , and a cross-matched subsample between these two datasets (3149 spiral and 1423 elliptical galaxies). For the cross-matched subsample, we find agreement in the slopes between observed and simulated galaxies. We find a slope of $3.8{6}_{-0.62}^{+0.92}$ for the full MaNGA sample, which agrees well with the slope of 4.0 predicted by MOND and the fitted slope of $3.5{8}_{-0.38}^{+0.48}$ for the TNG100 galaxies. We find that a sample of lower-mass galaxies is necessary to differentiate between the two models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Emulator-based Bayesian calibration of a subglacial drainage model

Subglacial drainage models, often motivated by the relationship between hydrology and ice flow, sensitively depend on numerous unconstrained parameters. We explore using borehole water-pressure time series to calibrate the uncertain parameters of a popular subglacial drainage model, taking a Bayesian perspective to quantify the uncertainty in parameter estimates and in the calibrated model predictions. To reduce the computation time associated with Markov Chain Monte Carlo sampling, we construct a fast Gaussian process emulator to stand in for the subglacial drainage model. We first carry out a calibration experiment using synthetic observations consisting of model simulations with hidden parameter values as a demonstration of the method. Using real borehole water pressures measured in western Greenland, we find meaningful constraints on four of the eight model parameters and a factor-of-three reduction in uncertainty of the calibrated model predictions. These experiments illustrate Gaussian process-based Bayesian inference as a useful tool for calibration and uncertainty quantification of complex glaciological models using field data. However, significant differences between the calibrated model and the borehole data suggest that structural limitations of the model, rather than poorly constrained parameters or computational cost, remain the most important constraint on subglacial drainage modelling.

58 GEOSCIENCES↗
DOE OSTI2023

Variability in modelled reservoir greenhouse gas emissions: comparison of select US hydropower reservoirs against global estimates

Abstract Greenhouse gas (GHG) emissions from reservoirs have most often been evaluated on a global extent through areal scaling or linear-regression models. These models typically rely on a limited number of characteristics such as age, size, and average temperature to estimate per reservoir or areal flux. Such approaches may not be sufficient for describing conditions at all types of reservoirs. Emissions from hydropower reservoirs have received increasing attention as industry and policy makers seek to better understand the role of hydropower in sustainable energy solutions. In the United States (US), hydropower reservoirs span a wide range of climate regions and have diverse design and operational characteristics compared to those most heavily represented in model literature (i.e., large, tropical reservoirs). It is not clear whether estimates based on measurements and modeling of other subsets of reservoirs describe the diverse types of hydropower reservoirs in the US. We applied the Greenhouse Gas from Reservoirs (G‐res) emissions model to 28 hydropower reservoirs located in a variety of ecological, hydrological, and climate settings that represent the range of sizes and types of facilities within the US hydropower fleet. The dominant pathways for resulting GHG emissions estimates in the case-study reservoirs were diffusion of carbon dioxide, followed by methane ebullition. Among these case-study reservoirs, total post-impoundment areal flux of carbon ranges from 84 to 767 mgCm −2 d −1 , which is less variable than what has been reported through measurements at other US and global reservoirs. The net GHG reservoir footprint was less variable and towards the lower end of the range observed from modeling larger global reservoirs, with a range of 138 to 1,052 g CO 2 eq m −2 y −1 , while the global study reported a range of 115 to 145,472 g CO 2 eq m −2 y −1 . High variation in emissions normalized with respect to area and generation highlights the need to be cautious when using area or generation in predicting or communicating emissions footprints for reservoirs relative to those of other energy sources, especially given that many of the hydropower reservoirs in the US serve multiple purposes beyond power generation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Adversarial sampling of unknown and high-dimensional conditional distributions

Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Herein, two approaches, stochastic estimation, and an external neural network are compared for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗