Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameter estimation using data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 343 records · Page 19

Multi-Modal Bayesian Neural Network Surrogates with Conjugate Last-Layer Estimation

As data collection and simulation capabilities advance, multi-modal learning, the task of learning from multiple modalities and sources of data, is becoming an increasingly important area of research. Surrogate models that learn from data of multiple auxiliary modalities to support the modeling of a highly expensive quantity of interest have the potential to aid outer loop applications such as optimization, inverse problems, or sensitivity analyses when multi-modal data are available. We develop two multi-modal Bayesian neural network surrogate models and leverage conditionally conjugate distributions in the last layer to estimate model parameters using stochastic variational inference (SVI). We provide a method to perform this conjugate SVI estimation in the presence of partially missing observations. Here, we demonstrate improved prediction accuracy and uncertainty quantification compared to unimodal surrogate models for both scalar and time series data.

97 MATHEMATICS AND COMPUTING↗

Machine learning modeling and model predictive control of a closed-circuit reverse osmosis system

Closed-circuit reverse osmosis (CCRO) offers a flexible and energy-efficient alternative to conventional reverse osmosis by operating in a semi-batch mode that recycles brine, enabling higher recovery rates and reduced specific energy consumption (SEC). However, developing accurate, system-level dynamic models for CCRO remains challenging due to its nonlinear, multi-phase operation and sensitivity to variable feed water conditions. Traditional modeling approaches, such as NARMAX (nonlinear autoregressive moving average with exogenous inputs), often struggle to generalize across varying inlet feed concentrations, necessitating frequent parameter re-estimation and limiting their utility for real-time control applications. To address these limitations, we developed a long short-term memory (LSTM) neural network model trained on an extensive experimental data set from a CCRO pilot plant. The model accepts three inputs, feed flow rate, recirculation flow rate, and initial feed conductivity, and predicts three key outputs: reject conductivity, feed pump power draw, and recirculation pump power draw. We validated the LSTM model against experimental data, demonstrating its ability to distinguish between different feed conductivities and adapt to variable flow rates. Subsequently, we incorporated the LSTM model within a nonlinear model predictive control (MPC) scheme and conducted closed-loop simulations to optimize the integrated SEC (iSEC). In conclusion, the results project up to a 6% reduction in iSEC by using MPC to optimize performance over the entire experiment duration, without requiring any random excitation for data collection or parameter re-estimation.

Desalination↗

Soil Water Characteristics of Cores from Low- and High-Centered Polygons, Utqiagvik (Barrow), Alaska, 2012

This dataset includes soil water characteristic curves for soil and permafrost in two representative frozen cores collected from a high-center polygon (HCP) and a low-center polygon (LCP) from the Barrow Environmental Observatory. Data include soil water content and soil water potential measured using the simple evaporation method for hydrological and biogeochemical simulations and experimental data analysis. Data can be used to generate a soil moisture characteristic curve, which can be fit to a variety of hydrological functions to infer critical parameters for soil physics. Considering the measured the soil water properties, the van Genuchten model predicted well the HCP, in contrast, the Kosugi model well fitted LCP which had more saturated condition. This dataset includes 8 zipped files, one txt files, and one pdf file.Data were collected to estimate soil water potential and water availability at various temperatures and conditions for hydrological and biogeochemical simulations and experimental data analysis. Data can be used to generate a soil moisture characteristic curve, which can be fit to a variety of hydrological functions to infer critical parameters for soil physics.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Landscape-scale consequences of differential tree mortality from catastrophic wind disturbance in the Amazon

Wind disturbance can create large forest blowdowns, which greatly reduces live biomass and adds uncertainty to the strength of the Amazon carbon sink. Observational studies from within the central Amazon have quantified blowdown size and estimated total mortality but have not determined which trees are most likely to die from a catastrophic wind disturbance. Also, the impact of spatial dependence upon tree mortality from wind disturbance has seldom been quantified, which is important because wind disturbance often kills clusters of trees due to large treefalls killing surrounding neighbors. We examine (1) the causes of differential mortality between adult trees from a 300-ha blowdown event in the Peruvian region of the northwestern Amazon, (2) how accounting for spatial dependence affects mortality predictions, and (3) how incorporating both differential mortality and spatial dependence affect the landscape level estimation of necromass produced from the blowdown. Standard regression and spatial regression models were used to estimate how stem diameter, wood density, elevation, and a satellite-derived disturbance metric influenced the probability of tree death from the blowdown event. The model parameters regarding tree characteristics, topography, and spatial autocorrelation of the field data were then used to determine the consequences of non-random mortality for landscape production of necromass through a simulation model. Tree mortality was highly non-random within the blowdown, where tree mortality rates were highest for trees that were large, had low wood density, and were located at high elevation. Of the differential mortality models, the non-spatial models overpredicted necromass, whereas the spatial model slightly underpredicted necromass. When parameterized from the same field data, the spatial regression model with differential mortality estimated only 7.5% more dead trees across the entire blowdown than the random mortality model, yet it estimated 51% greater necromass. As a result, we suggest that predictions of forest carbon loss from wind disturbance are sensitive to not only the underlying spatial dependence of observations, but also the biological differences between individuals that promote differential levels of mortality.

54 ENVIRONMENTAL SCIENCES↗

Cs-137 Body-burdens in Wild Pigs on the Department of Energy's Savannah River Site: Determination of an Appropriate Framework for Spatiotemporal Analysis of Hunter Dose Tracking System Data

The Savannah River Site (SRS) is a U.S. Department of Energy facility with an area of 310 square miles. Since its opening in 1951, SRS has produced and processed nuclear materials for national defense and industrial purposes. These activities have resulted in the release of radiocesium (Cs-137) which has led to contamination of associated streams, floodplains, man-made reservoirs, soil and vegetation. Contamination is spread to wild pigs as they move about SRS, consume vegetation, drink water and wallow in soil. SRS provides annual deer hunts open to the public in which wild pigs may be taken as well. The Cs-137 concentration of each animal is determined before it is released to the hunter for consumption and a cumulative dose is retained for each individual hunter. These tasks are accomplished using the Hunter Dose Tracking System (HDTS). The data collected on SRS by the HDTS for White-Tailed deer has been extensively analyzed. A repeatable method for calculating the background radiocesium body burden for the SRS deer population was established by Gaines and Novak (2009,2016). This model was framed using Cs-137 contamination obtained via gamma over-flight data, the decay rate of radiocesium and the historical behavior observed in the SRS deer population. These parameters were used in a spatiotemporal analysis to estimate background Cs-137 body burdens (from global fallout, etc.), as well as the effect SRS activities may have. The purpose of this project is to modify the framework to represent the behavior of the wild pig population on SRS such that HDTS data collected may be analyzed similarly. The 2018 data were collected via the HDTS in the field to screen wild pigs dispatched by Forest Service - Savannah River. These data will provide insight for an analysis of seasonal variance in the body burdens of the SRS wild pig population. The HDTS quickly provided Cs-137 concentrations in the field, as well as automatically updated the electronic database. Space: Historically, the average home range of wild pigs on the Savannah River Site has been reported at approximately six square miles. This range was confirmed accurate as of July 2018, with the observed minimum and maximum at one and ten square miles respectively. This range was geoprocessed using ArcMap to project buffers based on 1991 gamma over-flight data. The buffers were constructed to originate at all Cs-137 contamination point coordinates and protrude outward in all directions for 1.5, 2.45 and 6 square miles. While the 2.45 mile buffer is theoretically representative of potential hog exposure based on the average home range, it does not consider the possibility of home ranges not uniform in shape. Therefore, the 1.5 and 6 mile buffers were also constructed as visual aids for determining the appropriate hunt compartments for background level considerations. The SRS is divided into numbered hunt compartments of varying size and shape. These compartment numbers are used during sponsored hunts to record the location in which a wild pig was dispatched. Time: The Gaines report determined deer analysis should be carried out at the half-decade level for the best representation of true body burdens. This is due to the impact natural factors, such as drought or flood, can have on the bioaccumulation of Cs-137. To determine the appropriate time scale for hog analysis, all available data was totaled and averaged at 1, 2, 3 and 7 year increments. The annual averages and running average are provided. The appropriate compartments for the determination of the SRS hog population background level should be: 2, 3, 4, 9, 10, 14 and 15. Hogs dispatched within other compartments are most likely exposed to Cs-137 contamination and should not be included in the background level. The current dataset is small, but was still able to demonstrate the relationship between where hogs on SRS and the Cs-137 body burdens they incur. Geoprocessing the data using the average home range of wild pigs highlighted the importance of collecting additional samples from areas surrounding SRS if an accurate background is to be determined. When the average concentration of Cs-137 in pigs from background compartments is compared against that of pigs from contaminated compartments, a significant difference can be seen. A Table provides a summary of the analyzed data from the 15 years processed thus far.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimation of Uncertainties of Wall Temperature Accomplished by Model Calibration Technique Using Small to Large Scale Reflooding Test Data

The fidelity of the prediction can be evaluated by performing validation and uncertainty quantification with suitable test data. Nuclear power plant design and operation must accommodate the uncertainties in predicting system performance, where those uncertainties originating due to initial and boundary conditions, numerics, and modeling uncertainties. To improve the prediction of the reactor system performance, these various sources of uncertainties can be refined through model calibration technique. The objective of this work is to refine the model parameters and to estimate the uncertainties on parameters/responses during reflood phase. To accomplish this experiments should be performed, data collected, and the data utilized to complete model calibration for the thermal-hydraulic input parameters. Simulations were performed using Safety and Performance Analysis Code (SPACE) developed by multiple research institutes in South Korea to predict the thermal hydraulic system behaviors of nuclear power plants. Since existing experimental data were obtained from several tests performed at the small scale, scaling considerations must be addressed during the uncertainty quantification process. (authors)

42 ENGINEERING↗

Online Static Load Model Estimation in Distribution Systems

This paper focuses on parameter identification for time-varying voltage-dependent polynomial load models in distribution systems—namely, the constant impedance (Z), constant current (I), and constant power (P). The paper assumes the presence of measurements at the distribution nodes from various sources, each at varying temporal resolutions. The least- squares estimation method with bounded variables is used to convert the measurements into ZIP parameter estimations for loads at various locations on the feeder. The results and discussions focus on some of the estimation issues faced in a distribution system caused by varying voltage sensitivities of loads along a feeder, measurement resolution, impact of window size and sampling rate of measurement data. Finally, the paper uses dynamic simulations data (synthetic high-resolution measurements) to further improve the estimation results.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Dark Energy Survey Year 3 results: Simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps. II. cosmological results

Here, we present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years of the Dark Energy Survey. We implement the following: (1) second and third moments; (2) wavelet phase harmonics; (3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven 𝑤 Cold Dark Matter (𝑤⁢ CDM) cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the figure of merit (𝑆 8 , Ω m ) is improved by 70% (Λ ⁢CDM) and 90% (𝑤 ⁢CDM). When all the summary statistics are combined, we achieve a 2% constraint on the amplitude of fluctuations parameter 𝑆 8 ≡ 𝜎 8 ⁢(Ω m /0.3) 0.5 , obtaining 𝑆 8 = 0.794 ±0.017 (Λ⁢ CDM) and 𝑆 8 = 0.817 ±0.021 (𝑤 ⁢CDM), and a ∼10% constraint on Ω m , obtaining Ω m =0.259 ±0.025 (Λ ⁢CDM) and Ω m = 0.273 ±0.029 (𝑤⁢ CDM). In the context of the 𝑤⁢ CDM scenario, these statistics also strengthen the constraints on the parameter 𝑤, obtaining 𝑤 <−0.72. The constraints from different statistics are shown to be internally consistent (with a 𝑝-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the first three years of the Dark Energy Survey data, finding good consistency; we also compare with results from external datasets, such as planck constraints from the cosmic microwave background, finding statistical agreement, with discrepancies no greater than <2.2⁢𝜎.

79 ASTRONOMY AND ASTROPHYSICS↗

Covariance matrices for variance-suppressed simulations

ABSTRACT Cosmological N-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that ezmock (Effective Zel’dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the ezmock bias model.

79 ASTRONOMY AND ASTROPHYSICS↗

Repulsive Tomonaga-Luttinger liquid in the quasi-one-dimensional alternating spin-$\frac{1}{2}$ antiferromagnet NaVOPO 4

Here, we probe the magnetic field-induced Tomonaga-Luttinger liquid (TLL) state in the bond-alternating spin1/2 antiferromagnetic (AFM) chain compound NaVOPO 4 using thermodynamic as well as local μSR and 31 P NMR probes down to mK temperatures in magnetic fields up to 14 T. The μSR and NMR relaxation rates in the gapless TLL regime decay slowly following characteristic power-law behavior, enabling us to directly determine the interaction parameter K as a function of the magnetic field. These estimates are crosschecked using magnetization and specific heat data. The field-dependent K lies in the range of 0.4 < K < 1 and indicates the repulsive nature of interactions between the spinless fermions, in line with the theoretical predictions. This renders NaVOPO 4 the first experimental realization of TLL with repulsive fermionic interactions in hitherto studied S = 1/2 bond-alternating AFM-AFM chain compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parameter Identification for Battery Abuse Reaction Kinetics Models Using DSC Analysis [SWR-23-87]

This code estimates the parameters as inputs for NREL’s abuse reaction kinetics model using Differential scanning calorimetry (DSC) experimental data. The parameters for calculating the rates of abuse reactions among battery cell components at elevated temperatures include frequency factor, activation energy, and reaction orders. The least-squared method is employed to fit the reaction model to DSC data. The code is implemented in MATLAB.

Yang, Chuanbo↗

Uncertainties of the Yn Parameters of the Hage-Cifarelli Formalism

One method for determining the physical parameters of a multiplying system is summarized by Cifarelli [1]. In this methodology the single, double and triple rates are determined from what is commonly referred to as Feynman histograms. This paper will examine two methods for estimating the uncertainty in the parameters used in inferring these rates. These methods will be compared with simulated data in order to determine which one best approximates the sample uncertainty.

Instrumentation Related to Nuclear Science & Techn↗

Joint Bayesian Inference for Near-surface Explosion Yield

Geophysical observations of an explosion scale with the yield of the explosion. Guided by the principal of yield-scaling and observed relationships in the data, we defined and parameterized seismoacoustic forward models in Ford et al. [2014]. This work extends the model with additional data features of maximum velocity and positive phase duration for the seismic and acoustic data types, respectively, and a new optic data type. Parameters for the seismoacoustic models are estimated via Bayesian inference. The parameterized models can then be used to infer the yield of a new explosion.

58 GEOSCIENCES↗

Toward Development of a Stochastic Wake Model: Validation Using LES and Turbine Loads

Wind turbines within an array do not experience free-stream undisturbed flow fields. Rather, the flow fields on internal turbines are influenced by wakes generated by upwind unit and exhibit different dynamic characteristics relative to the free stream. The International Electrotechnical Commission (IEC) standard 61400-1 for the design of wind turbines only considers a deterministic wake model for the design of a wind plant. This study is focused on the development of a stochastic model for waked wind fields. First, high-fidelity physics-based waked wind velocity fields are generated using Large-Eddy Simulation (LES). Stochastic characteristics of these LES waked wind velocity field, including mean and turbulence components, are analyzed. Wake-related mean and turbulence field-related parameters are then estimated for use with a stochastic model, using Multivariate Multiple Linear Regression (MMLR) with the LES data. To validate the simulated wind fields based on the stochastic model, wind turbine tower and blade loads are generated using aeroelastic simulation for utility-scale wind turbine models and compared with those based directly on the LES inflow. The study’s overall objective is to offer efficient and validated stochastic approaches that are computationally tractable for assessing the performance and loads of turbines operating in wakes.

17 WIND ENERGY↗

Search for an Exotic S=-2, Q=-2 Baryon Resonance in Proton-Proton Interactions at $\sqrt{s_{NN}}$ = 17.3 GeV

Searches and analyses of strong gravitational lenses are challenging due to the rarity and image complexity of these astronomical objects. Next-generation surveys (both ground- and space-based) will provide more opportunities to derive science from these objects, but only if they can be analyzed on realistic time-scales. Currently, these analyses are expensive. In this work, we present a regression analysis with uncertainty estimates using deep learning models to measure four parameters of strong gravitational lenses in simulated Dark Energy Survey data. Using only $gri$-band images, we predict Einstein Radius, lens velocity dispersion, lens redshift to within $10-15\%$ of truth values and source redshift to $30\%$ of truth values, along with predictive uncertainties. This work helps to take a step along the path of faster analyses of strong lenses with deep learning frameworks.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Substrate Utilization, Metabolite Production, and Uranium Immobilization in Shewanella oneidensis Biofilms

In this study, we developed a two-dimensional mathematical model to predict substrate utilization and metabolite production rates in Shewanella oneidensis MR-1 biofilm in the presence and absence of uranium (U). In our model, lactate and fumarate are used as the electron donor and the electron acceptor, respectively. The model includes the production of extracellular polymeric substances (EPS). The EPS bound to the cell surface and the EPS distributed in the biofilm were considered bound EPS (bEPS) and loosely associated EPS (laEPS), respectively. COMSOL ® Multiphysics finite element analysis software was used to solve the model numerically (model file provided in the Supplementary Material). The input variables of the model were the lactate and fumarate concentrations, cell density, EPS density, half saturation constant for fumarate, and diffusion coefficients of the substrates and metabolites. To estimate unknown parameters and calibrate the model, we used a custom designed biofilm reactor placed inside a nuclear magnetic resonance (NMR) microimaging and spectroscopy system and measured substrate utilization and metabolite production rates. From these data we estimated the yield coefficients, maximum substrate utilization rate, half saturation constant for lactate, stoichiometric ratio of fumarate and acetate to lactate and stoichiometric ratio of succinate to fumarate. These parameters are critical to predicting the activity of biofilms and are not available in the literature. Lastly, the model was used to predict uranium immobilization in S. oneidensis MR-1 biofilms by considering reduction and adsorption processes in the cells and in the EPS. We found that the majority of immobilization was due to cells, and that EPS were less efficient at immobilizing U. Furthermore, most of the immobilization occurred within the top 10 µm of the biofilm. To the best of our knowledge, this research is one of the first biofilm immobilization mathematical models based on experimental observation. It has the ability to predict the relative contributions to U immobilization of laEPS, bEPS and cells.

59 BASIC BIOLOGICAL SCIENCES↗

A Machine Learning Based System Performance Prediction Model for Transportable FHR

A machine learning based system performance prediction model is currently created to support the development of autonomous control for small reactors, such as the Transportable Fluoride-salt-cooled High-temperature Reactor (TFHR), which has a 20 MWth compact core proposed by MIT for remote sites. The prediction model is constructed using support vector regression (SVR) with training data generated by RELAP5. A particle filtering framework is utilized to estimate and update model parameters with instrument measurements. Verifications of the prediction and filtering models have been carried out using TFHR reactivity insertion cases. Satisfactory performance in predicting the core behavior and in recognizing the inserted reactivity rate is concluded. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗