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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗
DOE OSTI2025

A semi-analytic estimate for the effective sound speed counterterm in the EFTofLSS

The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. Here, we apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value.

Cosmological perturbation theory in GR and beyond↗
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Elliptic flow of charged particles at midrapidity relative to the spectator plane in Pb–Pb and Xe–Xe collisions

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons v 2 {$Ψ_{SP}$} in collisions of Pb ions at center-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}=$ 2.76 TeV and Xe ions at $\sqrt{s_{NN}}=$ 5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum for the 5–70% and 0.2–6 GeV/c ranges, respectively. The ratio between v 2 {$Ψ_{SP}$} and the elliptic flow coefficient relative to the participant plane v 2 {4}, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the T R ENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models describing the initial state fluctuations. A significant transverse momentum dependence of the ratio v 2 {$Ψ_{SP}$} /v 2 {4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of v 2 {$Ψ_{SP}$} and v 2 {4} to the corresponding initial state eccentricities for Xe–Xe and Pb–Pb collisions at similar initial entropy density show a difference of (7.0 ± 0.9)% with an additional variation of +1.8% when including RHIC data in the T R ENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark–gluon plasma produced in heavy-ion collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Structure-aware Initialization via Numerical Continuation and Informed Priors

Scientific machine learning (SciML) often operates in ill-conditioned, weakly identifiable regimes due to limited data or indirect observations. In such settings, optimization and inference are highly sensitive to the starting point, making initialization--often under-reported--a consequential degree of freedom. Random initialization is not a neutral default as it induces an implicit prior over candidate solutions and can systematically bias the result, producing large run-to-run variability. Here, we formalize this view by treating initialization as a hidden confounder in SciML and develop a unifying theory for structure-aware initialization via numerical continuation, constructing warm starts from related problem instances. Across representative tasks, including physics-informed neural networks, maximum likelihood estimation, and variational inference, warm starts have been shown to consistently reduce optimization effort and improve reliability.

Data integrity↗
DOE OSTI2026

Conditional Latent Diffusion for High-Resolution Prediction of Electrochemical Surface Morphology

A conditionally guided generative latent diffusion process that is trained on a set of experimental processing parameters and their associated resulting electron microscope images of the electrodeposition process is able to interpolate between processing parameters in a physically consistent way. Electrodeposition of rhenium with pulse and pulse-reverse waveforms is used as a model system, and the process is adaptable to other electrodeposition, electropolishing, or corrosion processes. The method is able to extrapolate, predicting estimates of material morphologies for experimental setups unseen in the training data. The results are demonstrated with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Predictive Model for Internal Magazine Temperatures

This report describes data and methodology used to predict temperatures inside magazines used to store explosive devices, based on outside temperatures from a local weather station and data logged by sensors inside a single magazine over a limited period. It follows on earlier work and predictive models. In conjunction with aging models, predicted temperatures can be used to estimate degradation in the performance of explosives over long periods of storage.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1977

The spectrum of the geoid from altimeter data

A variety of sources of detailed information has been analyzed to arrive at a geoid power spectrum from global altimeter data. Using the equivalent of only two revolutions of data (mostly from GEOS-3) from all the major oceans, the high frequency geoid power (rms) is estimated (most simply) to be 80.7 n to the minus 1.47th power meters, where n is in cycles/global revolutions. This law is valid for all frequencies above 19 cycles but includes sea state. The (simple) law has more power than predicted by Kaula's rule for the geopotential. However, the data shows significantly less power for frequencies below 100 cycles. A closer approximation to the altimetry accumulates 2.18m (rss) for all frequencies higher than 19 cycles/rev. (including sea state), somewhat less power than predicted by the rule. The data permits up to 1.25 (rms) non-gravitational departures from the high frequency marine geoid.

Wagner, C. A.↗
NASA NTRS1989

Cratering of the terrestrial planets by Apollo objects

An asteroidal collision model and Monte Carlo program used for studies of the terrestrial meteorite flux, the steady-state number of Apollo-Amor objects, and the orbital distribution of both meteorites and Apollo-Amor objects is used to calculate absolute and relative cratering rates on the terrestrial planets. It is found that the 'best' estimates of the predicted asteroidal cratering rate are three times lower than estimates of the observed terrestrial cratering rate. If this is due to errors in the asteroidal production rate of Apollo-Amor objects, the predicted present-day cratering rate per unit area on Mars is four times that on earth, whereas that on Mercury is twice that on earth.

Wetherill, G. W.↗
NASA NTRS1991

Approximate effect of parameter pseudonoise intensity on rate of convergence for EKF parameter estimators

When using parameter estimation methods based on extended Kalman filter (EKF) theory, it is common practice to assume that the unknown parameter values behave like a random process, such as a random walk, in order to guarantee their identifiability by the filter. The present work is the result of an ongoing effort to quantitatively describe the effect that the assumption of a fictitious noise (called pseudonoise) driving the unknown parameter values has on the parameter estimate convergence rate in filter-based parameter estimators. The initial approach is to examine a first-order system described by one state variable with one parameter to be estimated. The intent is to derive analytical results for this simple system that might offer insight into the effect of the pseudonoise assumption for more complex systems. Such results would make it possible to predict the estimator error convergence behavior as a function of the assumed pseudonoise intensity, and this leads to the natural application of the results to the design of filter-based parameter estimators. The results obtained show that the analytical description of the convergence behavior is very difficult.

Hill, Bryon K.↗
NASA NTRS2008

Assessments of Hollow Cathode Wear in the Xenon Ion Propulsion System (XIPs(c)) by Numerical Analyses and Wear Tests

The standard approach presently followed by NASA to qualify electric propulsion for the required mission throughput has been based largely on life tests, which can be costly and time consuming. Revised electric propulsion lifequalification approaches are being formulated that combine analytical and/or computational methods with (shorter-duration) wear tests. As a model case, a wear test is being performed at JPL to assess the lifetime of the discharge hollow cathode in the Xenon Ion Propulsion System (XIPS(c)), a 25-cm ion engine developed by L-3 Communications Electron Technologies, Inc. for commercial applications. Wear and plasma data accumulated throughout this life-assessment program are being used to validate the existing 2-D hollow cathode code OrCa2D. We find that the OrCa2D steady-state solution predicts very well the time-averaged plasma data and the keeper voltage after 5500 hrs of operation in high-power mode. When the wave motion that occurs naturally in these devices is accounted for, based on an estimate of the maximum wave amplitude, the molybdenum-keeper erosion profile observed in the XIPS(c) discharge cathode is also reproduced within a factor of two of the observation. When the same model is applied to predict the erosion of a tantalum keeper we find that erosion is reduced by more than two orders of magnitude compared to the molybdenum keeper due the significantly lower sputtering yield of tantalum. A tantalum keeper would therefore allow keeper lifetimes that greatly exceed the present requirements for deep-space robotic missions considered by NASA. Moreover, such large reduction of the erosion renders the largest uncertainties in the models, which are associated with the wave amplitude estimates and the electron transport model, negligible.

plasma simulations↗
NASA NTRS2018

The SMAP Level-4 ECO Product - Phase 1: Improving Vegetation Simulations Through Observation-Driven Parameter Estimation

Simulations of hydrological fields as well as water, energy and carbon fluxes from the land surface to the atmosphere are crucial for a wide range of applications, including agricultural advisories, forecasts of (short-term) atmospheric behavior and seasonal weather predictions including forecasts of extreme events, such as heatwaves or droughts. The NASA Soil Moisture Active Passive (SMAP) mission Level-4 (L4) Eco-Hydrology (ECO) project aims to improve modeled estimates of the terrestrial water, energy and carbon fluxes and states by developing a fully-coupled hydrology-vegetation data assimilation (DA) algorithm. The DA system is developed for the NASA Goddard Earth Observing System version 5 (GEOS-5) Catchment-CN land surface model, which combines land hydrology components of the GEOS-5 Catchment model with dynamic vegetation components of the Community Land Model version 4. Catchment-CN fully couples the terrestrial water, energy and carbon cycles, allowing feedbacks from the land hydrology to the biosphere and vice versa. For SMAP L4 ECO a calibration of the Catchment-CN vegetation parameterization against observations of the fraction of absorbed photosynthetically active radiation (FPAR) from the Moderate Resolution Imaging Spectroradiometer (MODIS) is implemented to improve the model's standalone skill. Next, the DA algorithm used to produce the SMAP L4 soil moisture product is adapted to Catchment-CN to assimilate SMAP brightness temperatures and inform the model's land hydrology component. The DA system is further extended to assimilate MODIS FPAR observations in order to constrain the model's dynamic vegetation component. In this presentation, we demonstrate that the Catchment-CN parameter calibration leads to more realistic vegetation simulations and reduces the root mean squared error between modeled and observed vegetation states across the model's various plant functional types. We also show that the assimilation of SMAP observations is able to improve the average correlation, bias and unbiased RMSE between the modeled surface and root zone soil moisture estimates, and ground observations from the SMAP core validation sites.

Kolassa, J.↗
DOE OSTI2022

Phoebe: a high-performance framework for solving phonon and electron Boltzmann transport equations

Understanding the electrical and thermal transport properties of materials is critical to the design of electronics, sensors, and energy conversion devices. Computational modeling can accurately predict material properties but, in order to be reliable, requires accurate descriptions of electron and phonon states and their interactions. While first-principles methods are capable of describing the energy spectrum of each carrier, using them to compute transport properties is still a formidable task, both computationally demanding and memory intensive, requiring integration of fine microscopic scattering details for estimation of macroscopic transport properties. To address this challenge, we present Phoebe—a newly developed software package that includes the effects of electron–phonon, phonon–phonon, boundary, and isotope scattering in computations of electrical and thermal transport properties of materials with a variety of available methods and approximations. This open source C++ code combines MPI-OpenMP hybrid parallelization with GPU acceleration and distributed memory structures to manage computational cost, allowing Phoebe to effectively take advantage of contemporary computing infrastructures. We demonstrate that Phoebe accurately and efficiently predicts a wide range of transport properties, opening avenues for accelerated computational analysis of complex crystals.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Resolution of Orbit Determination Prediction Instabilities at Titan During Cassini's Solstice Mission

The Cassini spacecraft has been in orbit about Saturn since 2004. Exploration of the Saturn system is driven by gravitational flybys of the moon Titan which alter the spacecraft trajectory. The Cassini Navigation Team receives regular updates to the Saturn satellites ephemeris from JPL's Solar System Dynamics group. The difference between subsequent ephemeris deliveries can be hundreds of meters in the position of Titan at the time of a flyby. Errors in Titan's position propagate downstream to the next flyby through the estimated spacecraft trajectory. Prior to 2013, the Cassini Orbit Determination Team estimated the Saturn satellite ephemeris parameters and used the a posteriori states and covariance of an operations arc as a priori inputs to subsequent estimation arcs. Since 2013, the OD Team has only been considering errors in the ephemeris and not estimating a correction to the satellite positions. The T119 Titan flyby exhibited a 3D miss distance of 2.44 km and the following T120 flyby yielded a smaller miss of 1.06 km at the 2.9s error level. These discrepancies between pre-flyby prediction and post-flyby trajectory reconstruction were due to errors in the Titan ephemeris. In order to improve the targeting of Titan in future flybys, the team restarted the satellite ephemeris estimation process for orbit determination solutions. Subsequent flybys had target misses of less than 1 km at the sub-3 error level. This paper describes the method of scaling the a priori satellite ephemeris covariance in the orbit determination process to allow larger corrections to the satellite system and improve the prediction of the spacecraft’s Titan-relative position at the time of encounters.

Boone, Dylan↗
DOE OSTI2024

Elastic deformation as a tool to investigate watershed storage connectivity

Storage-discharge relationships and dynamic changes in storage connectivity remain key unknowns in understanding and predicting watershed behavior. In this study, we use Global Positioning System measurements of load-induced Earth surface displacement as a proxy for total water storage change in four climatologically diverse mountain watersheds in the western United States. Comparing total water storage estimates with stream-connected storage derived from hydrograph analysis, we find that each of the investigated watersheds exhibits a characteristic seasonal pattern of connection and disconnection between total and stream-connected storage. We investigate how the degree and timing of watershed-scale connectivity is related to the timing of precipitation and seasonal changes in dominant hydrologic processes. Our results show that elastic deformation of the Earth due to water loading is a powerful new tool for elucidating dynamic storage connectivity and watershed discharge response across scales in space and time.

Environmental Sciences & Ecology↗
DOE OSTI2023

Constraining the $\overline{K}$N coupled channel dynamics using femtoscopic correlations at the LHC

The interaction of K – with protons is characterised by the presence of several coupled channels, systems like $\overline{K}$ 0 n and πΣ with a similar mass and the same quantum numbers as the K – p state. The strengths of these couplings to the K – p system are of crucial importance for the understanding of the nature of the Λ(1405) resonance and of the attractive K – p strong interaction. In this article, we present measurements of the K – p correlation functions in relative momentum space obtained in pp collisions at $\sqrt{s}$=13 Te, in p–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te, and (semi)peripheral Pb–Pb collisions at $\sqrt{s_{NN}}$= 5.02 Te. The emitting source size, composed of a core radius anchored to the K + p correlation and of a resonance halo specific to each particle pair, varies between 1 and 2 fm in these collision systems. The strength and the effects of the $\overline{K}$ 0 n and πΣ inelastic channels on the measured K – p correlation function are investigated in the different colliding systems by comparing the data with state-of-the-art models of chiral potentials. A novel approach to determine the conversion weights ω, necessary to quantify the amount of produced inelastic channels in the correlation function, is presented. In this method, particle yields are estimated from thermal model predictions, and their kinematic distribution from blast-wave fits to measured data. The comparison of chiral potentials to the measured K – p interaction indicates that, while the πΣ –K – p dynamics is well reproduced by the model, the coupling to the $\overline{K}$ 0 n channel in the model is currently underestimated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗