Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic 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 775 records · Page 43

Lunar gravity - A long-term Keplerian rate method

Nearside tracking data from Apollo subsatellites and Lunar Orbiter 5 are used to construct a model of the farside lunar gravity field. An indirect method of calculating the average effect of all the gravity anomalies which affect spacecraft orbits is reduced to a linear estimation problem by eliminating the integration of spacecraft dynamics. Doppler frequency shifts between the lunar satellite and earth-based tracking station are used to measure the line of sight velocity component between the orbiter and the tracker; gravity determination is then performed in a weighted least squares estimator which employs the orbital rates as data and formulates the normal equations using the associated Keplerian state. Simulations performed demonstrate that farside gravity features can be economically determined by use of gravity estimates obtained from long-period models of nearside orbiting tracks.

Ferrari, A. J.↗

State-space self-tuner for on-line adaptive control

Dynamic systems, such as flight vehicles, satellites and space stations, operating in real environments, constantly face parameter and/or structural variations owing to nonlinear behavior of actuators, failure of sensors, changes in operating conditions, disturbances acting on the system, etc. In the past three decades, adaptive control has been shown to be effective in dealing with dynamic systems in the presence of parameter uncertainties, structural perturbations, random disturbances and environmental variations. Among the existing adaptive control methodologies, the state-space self-tuning control methods, initially proposed by us, are shown to be effective in designing advanced adaptive controllers for multivariable systems. In our approaches, we have embedded the standard Kalman state-estimation algorithm into an online parameter estimation algorithm. Thus, the advanced state-feedback controllers can be easily established for digital adaptive control of continuous-time stochastic multivariable systems. A state-space self-tuner for a general multivariable stochastic system has been developed and successfully applied to the space station for on-line adaptive control. Also, a technique for multistage design of an optimal momentum management controller for the space station has been developed and reported in. Moreover, we have successfully developed various digital redesign techniques which can convert a continuous-time controller to an equivalent digital controller. As a result, the expensive and unreliable continuous-time controller can be implemented using low-cost and high performance microprocessors. Recently, we have developed a new hybrid state-space self tuner using a new dual-rate sampling scheme for on-line adaptive control of continuous-time uncertain systems.

Shieh, L. S.↗

Ensemble Forecasting/Assimilation/Emissions Estimation with WRF-Chem/DART: Accomplishments, Lessons Learned, and Future Plans

Over the past ten years, we have been conducting research on regional ensemble atmospheric composition forecasting/data assimilation/emissions estimation with WRF-Chem/DART. WRF- Chem/DART integrates the Weather Research and Forecasting model (WRF) with online chemistry (WRF-Chem) into the Data Assimilation Research Testbed (DART). DART is an ensemble data assimilation system based on the ensemble adjustment Kalman filter (EAKF) with adaptive inflation, localization (physical and state space), and an optional non-Gaussian formulation of the EAKF. DART includes assimilation of meteorological and limited chemical observations. WRF-Chem/DART extends DART to include assimilation of: MOPITT CO; IASI CO and O3; MODIS AOD; OMI O3, NO2, and SO2; TROPOMI CO, O3, NO2, and SO2, TES CO, CO2 (research mode), O3, NH3, and CH4 (research mode); CrIS CO, O3, NH3, CH4 (research mode), and PAN; SCIAMACHY NO2; GOME2a NO2; MLS O3 and HNO3; and proxy TEMPO O3, and NO2 satellite retrievals as raw retrievals or as ‘compact phase space retrievals’ (CPSRs) for profile retrievals. WRF-Chem/DART also assimilates in situ atmospheric composition measurements and uses the ‘state augmentation method’ for emissions estimation. In our presentation, we will provide an overview of WRF-Chem/DART: • applications and results; • lessons learned from: (i) independent versus joint assimilation; (ii) total/partial column versus profile retrieval assimilation; (iii) joint in situ and retrieval assimilation; (iv) assimilation at grid resolutions ranging from 100 km to 4 km; (v) dynamic emissions estimation; and • future work related to wildfire emissions estimation and intercomparison of CMAQ (online)/JEDI and CMAQ (online and offline)/DART.

WRF-Chem/DART↗

Strategy for Risk Quantification of Spaceflight Crew Health and Performance Using Dynamic Probabilistic Risk Assessment

At NASA, the Crew Health and Performance (CHP) system represents the span of countermeasures, capabilities, interventions, and tested processes and procedures that in combination work to mitigate the human component of spaceflight mission risk. Across the varying NASA mental models of the CHP system, the different functionalities needed to meet human flight systems standards can be broken down into specific categories (i.e. medical capability, environmental health, behavioral health). These categories can be further broken down into specific subgroups generally associated with the CHP functionalities meant to mitigate or buy down individual human system risks. Taking a similar development approach we seek to leverage dynamic probabilistic risk assessment as a means to quantify and relatively assess the human risk state within the crew health and performance domain. By utilizing existing tools as integrators, we propose a rapid development strategy for incorporating research and operational data that represent the influence of the CHP system functionalities, in order to provide order of magnitudes estimates of the influence on most human system risks outcomes. The model system utilizes a modest cumulative risk approach and that limits the scope to primary paths of influence between the CHP functionalities and human system risks, thus enabling quick prototypes of the integrative effects of CHP functional combinations to solicit valuable feedback from stakeholders and customers on the data, relationship, and structure of the integration.

Drayton Munster↗

Anticipating impacts of introducing aluminum-cerium alloys into the United States automotive market

Increasing demand of abundant co-product rare earths such as cerium and lanthanum has become a focus for improving mine economics. The recent development of an aluminum-cerium (Al-Ce) alloy targeted toward broad market adoption in the automotive industry is one recent example. This alloy has shown exceptional mechanical properties and cost-saving potential. However, since market adoption is unknown, impacts of this alloy on Ce demand remains uncertain. We estimate Al-Ce impacts by simulating adoption in different automotive sectors using a dynamic material flow analysis model. We then project future Ce content in the recycled Al stream given current recycling practices. We also compute the required segregation given different Ce impurity tolerances of existing automotive Al alloys. Finally, once secondary supply is derived from segregation, we calculate primary Ce demand. Results showed that segregation of Al-Ce components from non-AlCe automotive Al components is needed for most adoption scenarios and impurity thresholds. Ce metal demand was found to be lowest at 1000 mt/yr if the target market is light duty trucks. If Al-Ce could penetrate the High Efficiency Spark Ignition market (i.e., turbocharger market), Ce demand could reach 23,200 mt Ce/year in 2045, without accounting for Chinese and Australian Ce primary production capacity. Furthermore, when accounting for primary Ce production capacity, demand could reach up to 16,200 mt/yr.

99 GENERAL AND MISCELLANEOUS↗

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

Changes in soil water content and lateral flow exert large effects on soil thermal dynamics across Alaskan landscapes

Both lateral surface and subsurface water flow affect soil moisture dynamics, yet most land surface models only solve subsurface water movement vertically. Here, we use a 3D ecosystem model that considers both land surface and subsurface hydrologic processes to simulate soil moisture, which is then used to drive a 1-D vertical soil thermal model to simulate the soil moisture effects on soil thermal dynamics in central Alaska. Our coupled model improves soil temperature (ST) estimates by 43.5% in comparison with observational data. Soil moisture has little effect on ST during the wet season (-1.5%) and a substantial influence during the dry season (60%). Spatially, water lateral flow has significant impacts on both soil moisture and ST, causing model estimates for thawed areas in the transition season to increase by ~10% in the study area. Our results highlight the importance of considering dynamical soil moisture, as well as lateral flow effects, on soil thermal dynamics in permafrost regions.

54 ENVIRONMENTAL SCIENCES↗

Biermann battery magnetic fields in ICF capsules: Total magnetic flux generation

This paper focuses on the process of magnetic flux generation in inertial confinement fusion implosions. Hot-spots are shown to be dominated by fields generated during stagnation when the temperature and density gradients are largest. A scaling of hot-spot magnetic flux is derived and compared with simulations, revealing that perturbations with both larger amplitudes and higher mode numbers generate more magnetic flux. Furthermore, the model allows for greater understanding of which target designs will be susceptible to magnetohydrodynamic effects. For example, the model can be used to ascertain the time when most magnetic flux is generated. If generation is weighted more toward early times, then more high-mode magnetic field loops will be present. A hot-spot with no high-mode perturbations at time of peak neutron production can still contain significant magnetic flux on those scales. By assuming that magnetic flux is deposited at the hot-spot edge by Nernst advection, the model can be used to post-process radiation-hydrodynamics data to estimate magnetic field strengths and magnetizations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗

Variation in hyphal production rather than turnover regulates standing fungal biomass in temperate hardwood forests

Abstract Soil fungi link above‐ and belowground carbon (C) fluxes through their interactions with plants and contribute to C and nutrient dynamics through the production, turnover, and activity of fungal hyphae. Despite their importance to ecosystem processes, estimates of hyphal production and turnover rates are relatively uncommon, especially in temperate hardwood forests. We sequentially harvested hyphal ingrowth bags to quantify the rates of Dikarya (Ascomycota and Basidiomycota) hyphal production and turnover in three hardwood forests in the Midwestern United States, where plots differed in their abundance of arbuscular (AM)‐ vs. ectomycorrhizal (ECM)‐associated trees. Hyphal production rates increased linearly with the percentage of ECM trees and annual production rates were 66% higher in ECM‐ than AM‐dominated plots. Hyphal turnover rates did not differ across the mycorrhizal gradient (plots varying in their abundance of AM vs. ECM trees), suggesting that the greater fungal biomass in ECM‐dominated plots relates to greater fungal production rather than slower fungal turnover. Differences in hyphal production across the gradient aligned with distinctly different fungal communities and activities. As ECM trees increased in dominance, fungi inside ingrowth bags produced more extracellular enzymes involved in degrading nitrogen (N)‐bearing relative to C‐bearing compounds, suggesting greater fungal (and possibly plant) N demand in ECM‐dominated soils. Collectively, our results demonstrate that shifts in temperate tree species composition that result in changes in the dominant type of mycorrhizal association may have strong impacts on Dikarya hyphal production, fungal community composition and extracellular enzyme activity, with important consequences for soil C and N cycling.

59 BASIC BIOLOGICAL SCIENCES↗

Improving estimates of sub-daily gross primary production from solar-induced chlorophyll fluorescence by accounting for light distribution within canopy

Solar-induced chlorophyll fluorescence (SIF) has long been regarded as a proxy for photosynthesis and has shown superiority in estimating gross primary production (GPP) compared to traditional vegetation indices, especially in evergreen ecosystems. However, current SIF-based GPP estimations regard the canopy as a large leaf and seldom consider the impact of interactions among light, canopy structure, and leaf physiology. In this study, we proposed GPP estimation models with different descriptions of light–structure–physiology interactions (including the layered model, the two-leaf model, and the layered two-leaf model) and compared their performances with the big-leaf model using half-hourly (or hourly) observations at evergreen needleleaf forest sites. First, we found that the big-leaf model underestimated GPP, especially at noon. All models showed higher accuracy than that of the big-leaf model. Second, we investigated the diurnal dynamics of GPP estimations in each canopy layer and found that models with a two-leaf assumption captured the diurnal variations in GPP better than that with the layered assumption. Here we also deduced that the poor performance of the big-leaf model was related to its overestimation of the overall light stress on the redox state of PSII reaction centers (qL). Finally, we noticed that the qL at the canopy scale had lower sensitivity to light change than the single-leaf qL and that the light response of canopy-scale qL was influenced by the leaf area index during seasonal cycles. Overall, this study describes methods to accurately estimate sub-daily GPP from SIF in evergreen needleleaf forests and demonstrates that the interactions among light, canopy structure, and leaf physiology regulate the SIF-GPP relationship at the canopy scale. Further, it indicates the need to consider the description of light distribution within the canopy in next-generation terrestrial biosphere models, even if they incorporate SIF to constrain their parameterization. Thus, upscaling the established leaf-scale mechanistic SIF-GPP relationship or findings to canopy-scale applications still requires much work, especially when there are significant changes in environmental conditions and their within-canopy distributions.

59 BASIC BIOLOGICAL SCIENCES↗

Use of nonlinear identification in robust attitude and attitude rate estimation for SAMPEX

A method is described for obtaining optimal attitude estimation/identification algorithms for spacecraft lacking attitude rate measurement devices (rate gyros), and then demonstrated using actual flight data from the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) spacecraft. SAMPEX does not have on-board rate sensing, and relies on sun sensors and a three-axis magnetometer for attitude determination. The absence of rate data normally reduces both the total amount of data available and the sampling density (in time) by a substantial fraction. In addition, attitude data is occasionally unavailable (for example, during sun occultation). As a result, the sensitivity of the estimates to model uncertainty and to measurement noise increases. In order to maintain accuracy in the attitude estimates, there is an increased need for accurate models of the rotational dynamics. The Minimum Model Error(MME)/Least Square Correlation(LSC) algorithm accurately identifies an improved model for SAMPEX to be used during periods of complete data loss or extreme noise. The model correction is determined by estimating only one orbit(the identification pass) just prior to the assumed data loss(the prediction pass). The MME estimator correctly predicted the states during the identification phase, but more importantly determines the necessary model correction trajectory, d(t). The LSC algorithm is then used to find this trajectory's functional form, H(x(t)). The results show significant improvement of the new corrected model's attitude estimates as compared to the original uncorrected model's estimates. The possible functional form of the correction term is limited at this point in the study to functions strictly of the estimated states. The results, however, strongly suggest that functions based on the relative position of the satellite may also be possible candidates for future consideration.

Mook, D. Joseph↗

Use of microgravity sensors for quantification of space shuttle orbiter vernier reaction control system induced environments

In the modeling of spacecraft dynamics it is important to accurately characterize the environment in which the vehicle operates, including the environments induced by the vehicle itself. On the Space Shuttle these induced environmental factors include reaction control system plume. Knowledge of these environments is necessary for performance of control systems and loads analyses, estimation of disturbances due to thruster firings, and accurate state vector propagation. During the STS-71 mission, while the Orbiter was performing attitude control for the mated Orbiter/Mir stack, it was noted that the autopilot was limit cycling at a rate higher than expected from pre-flight simulations. Investigations during the mission resulted in the conjecture that an unmodelled plume impingement force was acting upon the orbiter elevons. The in-flight investigations were not successful in determining the actual magnitude of the impingement, resulting in several sequential post-flight investigations. Efforts performed to better quantify the vernier reaction control system induced plume impingement environment of the Space Shuttle orbiter are described in this paper, and background detailing circumstances which required the more detailed knowledge of the RCS self impingement forces, as well as a description of the resulting investigations and their results is presented. The investigations described in this paper applied microgravity acceleration data from two shuttle borne microgravity experiments, SAMS and OARE, to the solution of this particular problem. This solution, now used by shuttle analysts and mission planners, results in more accurate propellant consumption and attitude limit cycle estimates in preflight analyses, which are critical for pending International Space Station missions.

Friend, Robert B.↗

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↗

A Multivariate Space‐Time Dynamic Model for Characterizing the Atmospheric Impacts Following the Mt. Pinatubo Eruption

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

Dynamic Linear Model↗

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid-liquid phase transition of hydrogen and its critical point: Analysis from ab initio simulation and a machine-learned potential

We simulate high-pressure hydrogen in its liquid phase close to molecular dissociation using a machine-learned interatomic potential. The model is trained with density functional theory (DFT) forces and energies, with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. We show that an accurate NequIP model, an E(3)-equivariant neural network potential, accurately reproduces the phase transition present in PBE. Moreover, the computational efficiency of this model allows for substantially longer molecular dynamics trajectories, enabling us to perform a finite-size scaling (FSS) analysis to distinguish between a crossover and a true first-order phase transition. Here, we locate the critical point of this transition, the liquid-liquid phase transition (LLPT), at 1200-1300 K and 155-160 GPa, a temperature lower than most previous estimates and close to the melting transition.

08 HYDROGEN↗

Quantum estimation in an expanding spacetime

Highlights: • FI and QFI of estimation on Hubble parameter of de Sitter space are explored. • QFI degrades due to Gibbons–Hawking radiation. • QFI can be enhanced for proper coupling of scalar field to curvature. • Optimal QFI depends on the squeezed nature of de Sitter α-vacua. - Abstract: We investigate the quantum estimation on the Hubble parameter of an expanding de Sitter space by quantum metrological techniques. By exploring the dynamics of a freely falling Unruh–DeWitt detector, which interacts with a scalar field coupling to curvature, we calculate the Fisher information (FI) and quantum Fisher information (QFI) for the detector, which bound the highest precision of the estimation on Hubble parameter. In standard Bunch–Davies vacuum, we show that the maxima of FI/QFI are located for particular initial state of probe. Beside its dependence on the evolving time of detector and the energy spacing of atom ω, we show that the maxima of FI/QFI can be significantly enhanced once a proper coupling of scalar field to curvature is chosen. For instance, we show numerically that the estimation in the scenario with minimally/nearly minimally coupling scalar field can always outperform that with conformally coupling scalar field, corresponding to a higher FI/QFI in estimation. Moreover, we find that for general α-vacua of de Sitter space, a further improvement of estimation can be achieved, attributed to the squeezed nature of α-vacua that heavily constrains the measurement uncertainty. Some implications of our results are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗