Search NASA⌕ Search

SEARCH · Search NASA

Results for “state space model”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reduced-order autodifferentiable ensemble Kalman filters

This paper introduces a computational framework to reconstruct and forecast a partially observed state that evolves according to an unknown or expensive-to-simulate dynamical system. Our reduced-order autodifferentiable ensemble Kalman filters (ROAD-EnKFs) learn a latent low-dimensional surrogate model for the dynamics and a decoder that maps from the latent space to the state space. The learned dynamics and decoder are then used within an EnKF to reconstruct and forecast the state. Numerical experiments show that if the state dynamics exhibit a hidden low-dimensional structure, ROAD-EnKFs achieve higher accuracy at lower computational cost compared to existing methods. If such structure is not expressed in the latent state dynamics, ROAD-EnKFs achieve similar accuracy at lower cost, making them a promising approach for surrogate state reconstruction and forecasting.

Mathematics↗
DOE OSTI2022

TAUOLA update for decay channels with e + e – pairs in the final state

With the arrival of high luminosity B-factories like the Belle II experiment, τ decay measurements have become more precise than ever, allowing rarer processes to be explored, and finer details of τ decays to be studied. These are important to understand the spectrum of intermediate particles produced in τ decays. Therefore Monte Carlo generators, like the TAUOLA program, have to facilitate precision analysis as well as confront new models that constantly emerge with the availability of high statistics experimental data. New decay channels and models may lead to large variation of matrix elements size within the available phase space, as in the case when e + e - pairs are present in final state. It requires appropriate presampler of the phase space generator and a proper documentation to help users introduce their own models. While releasing a new update, it is important for the TAUOLA Monte Carlo library to maintain the general structure of the previous versions to preserve backward compatibility. The aim is to minimize changes from the user perspective. Finally, this paper presents a demonstration of new models implementation facilitated by the current update.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Semi-supervised Learning of Dynamical Systems with Neural Ordinary Differential Equations: A Teacher-Student Model Approach

Modeling dynamical systems is crucial for a wide range of tasks, but it remains challenging due to complex nonlinear dynamics, limited observations, or lack of prior knowledge. Recently, data-driven approaches such as Neural Ordinary Differential Equations (NODE) have shown promising results by leveraging the expressive power of neural networks to model unknown dynamics. However, these approaches often suffer from limited labeled training data, leading to poor generalization and suboptimal predictions. On the other hand, semi-supervised algorithms can utilize abundant unlabeled data and have demonstrated good performance in classification and regression tasks. We propose TS-NODE, the first semi-supervised approach to modeling dynamical systems with NODE. TS-NODE explores cheaply generated synthetic pseudo rollouts to broaden exploration in the state space and to tackle the challenges brought by lack of ground-truth system data under a teacher-student model. TS-NODE employs an unified optimization framework that corrects the teacher model based on the student's feedback while mitigating the potential false system dynamics present in pseudo rollouts. TS-NODE demonstrates significant performance improvements over a baseline Neural ODE model on multiple dynamical system modeling tasks.

Wang, Yu↗
NASA NTRS1986

Grain dynamics in zero gravity

The dynamics of granular materials has proved difficult to model, primarily because of the complications arising from inelastic losses, friction, packing, and the effect of many grains being in contact simultaneously. The kinetic model of granular systems is similar to the kinetic theory of gases, except that collisional energy losses are always present in the former and must be treated explicity. Few granular materials on Earth are describable by this limiting model, since gravity tends to collapse the grains into a high density state where Coulombic friction effects are dominant. The planned Space Station offers an unusual opportunity to test the kinetic grain model and to explore its predictions. Without gravity, the regime of low interparticle velocities (where an elastic description of the collision is still valid) can be investigated. This will allow for direct interpretation by dynamical computer simulations as well as by the kinetic theory. The dynamics of spherical grains inside a clear box would be examined. Results would be compared with the predictions of the kinetic theory and computer simulations.

Werner, B. T.↗
NASA NTRS2012

Constructing an Efficient Self-Tuning Aircraft Engine Model for Control and Health Management Applications

Self-tuning aircraft engine models can be applied for control and health management applications. The self-tuning feature of these models minimizes the mismatch between any given engine and the underlying engineering model describing an engine family. This paper provides details of the construction of a self-tuning engine model centered on a piecewise linear Kalman filter design. Starting from a nonlinear transient aerothermal model, a piecewise linear representation is first extracted. The linearization procedure creates a database of trim vectors and state-space matrices that are subsequently scheduled for interpolation based on engine operating point. A series of steady-state Kalman gains can next be constructed from a reduced-order form of the piecewise linear model. Reduction of the piecewise linear model to an observable dimension with respect to available sensed engine measurements can be achieved using either a subset or an optimal linear combination of "health" parameters, which describe engine performance. The resulting piecewise linear Kalman filter is then implemented for faster-than-real-time processing of sensed engine measurements, generating outputs appropriate for trending engine performance, estimating both measured and unmeasured parameters for control purposes, and performing on-board gas-path fault diagnostics. Computational efficiency is achieved by designing multidimensional interpolation algorithms that exploit the shared scheduling of multiple trim vectors and system matrices. An example application illustrates the accuracy of a self-tuning piecewise linear Kalman filter model when applied to a nonlinear turbofan engine simulation. Additional discussions focus on the issue of transient response accuracy and the advantages of a piecewise linear Kalman filter in the context of validation and verification. The techniques described provide a framework for constructing efficient self-tuning aircraft engine models from complex nonlinear simulations.Self-tuning aircraft engine models can be applied for control and health management applications. The self-tuning feature of these models minimizes the mismatch between any given engine and the underlying engineering model describing an engine family. This paper provides details of the construction of a self-tuning engine model centered on a piecewise linear Kalman filter design. Starting from a nonlinear transient aerothermal model, a piecewise linear representation is first extracted. The linearization procedure creates a database of trim vectors and state-space matrices that are subsequently scheduled for interpolation based on engine operating point. A series of steady-state Kalman gains can next be constructed from a reduced-order form of the piecewise linear model. Reduction of the piecewise linear model to an observable dimension with respect to available sensed engine measurements can be achieved using either a subset or an optimal linear combination of "health" parameters, which describe engine performance. The resulting piecewise linear Kalman filter is then implemented for faster-than-real-time processing of sensed engine measurements, generating outputs appropriate for trending engine performance, estimating both measured and unmeasured parameters for control purposes, and performing on-board gas-path fault diagnostics. Computational efficiency is achieved by designing multidimensional interpolation algorithms that exploit the shared scheduling of multiple trim vectors and system matrices. An example application illustrates the accuracy of a self-tuning piecewise linear Kalman filter model when applied to a nonlinear turbofan engine simulation. Additional discussions focus on the issue of transient response accuracy and the advantages of a piecewise linear Kalman filter in the context of validation and verification. The techniques described provide a framework for constructing efficient self-tuning aircraft engine models from complex nonlinear simulatns.

Armstrong, Jeffrey B.↗
DOE OSTI2023

Linearized Pair-Density Functional Theory for Vertical Excitation Energies

Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.

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

Reviving MeV-GeV indirect detection with inelastic dark matter

Thermal relic dark matter below ∼ 10 GeV is excluded by cosmic microwave background data if its annihilation to visible particles is unsuppressed near the epoch of recombination. Usual model-building measures to avoid this bound involve kinematically suppressing the annihilation rate in the low-velocity limit, thereby yielding dim prospects for indirect detection signatures at late times. In this work, we investigate a class of cosmologically viable sub-GeV thermal relics with late-time annihilation rates that are detectable with existing and proposed telescopes across a wide range of parameter space. We study a representative model of inelastic dark matter featuring a stable state χ 1 and a slightly heavier excited state χ 2 whose abundance is thermally depleted before recombination. Since the kinetic energy of dark matter in the Milky Way is much larger than it is during recombination, χ 1 χ 1 → χ 2 χ 2 upscattering can efficiently regenerate a cosmologically long-lived Galactic population of χ 2 , whose subsequent coannihilations with χ 1 give rise to observable gamma-rays in the ∼ 1 MeV − 100 MeV energy range. We find that proposed MeV gamma-ray telescopes, such as e-ASTROGAM, AMEGO, and MAST, would be sensitive to much of the thermal relic parameter space in this class of models and thereby enable both discovery and model discrimination in the event of a signal at accelerator or direct detection experiments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Control-structure interaction

The feasibility of using conventional proportional-integral-derivative (PID) control and an alternative optimal control to perform the pointing and tracking functions of the Space Station solar dynamic power module is investigated. A very large state model of 6 rigid body modes and 272 flexible modes is used in conjunction with classical linear-quadratic-Gaussian (LQG) optimal control to produce a full-order controller that satisfies the requirements. The results are compared with a classically designed PID controller that was implemented for a much smaller (6 rigid body, 40 flexible modes) model. The conventional control design approach is shown to be very much influenced by the order reduction of the plant model, i.e., the number of retained elastic modes from the full-order model, suggesting that for a complex, large space structure, such as the Space Station Freedom solar dynamic module, application of conventional control system design methods may not be adequate. The use of LQG control is recommended, and method for solving the large matrix. Riccati equation that arises from the optimal formulation is provided.

Cheng, Joseph K.↗
DOE OSTI2019

Enhanced correlations and superconductivity in weakly interacting partially flat-band systems: A determinantal quantum Monte Carlo study

Motivated by recent experiments realizing correlated phenomena and superconductivity in two-dimensional (2D) van der Waals devices, we consider the general problem of whether correlation effects may be enhanced by modifying band structure while keeping a fixed weak interaction strength. Using determinantal quantum Monte Carlo, we study the 2D Hubbard model for two different band structures: a regular nearest-neighbor tight-binding model and a partially flat-band structure containing a nondispersing region, with identical total noninteracting bandwidth W tot . For both repulsive and attractive weak interactions (|U| $\ll$ W$_{tot}$, correlated phenomena are significantly stronger in the partially flat model. In the repulsive case, even with $U$ being an order of magnitude smaller than W$_{tot}$, we find the presence of a Mott insulating state near half filling of the flat region in momentum space. In the attractive case, where generically the ground state is superconducting, the partially flat model exhibits significantly enhanced superconducting transition temperatures. These results suggest the possibility of engineering correlation effects in materials by tuning the noninteracting electronic dispersion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2016

Lidar-Based Estimates of Above-Ground Biomass in the Continental US and Mexico Using Ground, Airborne, and Satellite Observations

Existing national forest inventory plots, an airborne lidar scanning (ALS) system, and a space profiling lidar system (ICESat-GLAS) are used to generate circa 2005 estimates of total aboveground dry biomass (AGB) in forest strata, by state, in the continental United States (CONUS) and Mexico. The airborne lidar is used to link ground observations of AGB to space lidar measurements. Two sets of models are generated, the first relating ground estimates of AGB to airborne laser scanning (ALS) measurements and the second set relating ALS estimates of AGB (generated using the first model set) to GLAS measurements. GLAS then, is used as a sampling tool within a hybrid estimation framework to generate stratum-, state-, and national-level AGB estimates. A two-phase variance estimator is employed to quantify GLAS sampling variability and, additively, ALS-GLAS model variability in this current, three-phase (ground-ALS-space lidar) study. The model variance component characterizes the variability of the regression coefficients used to predict ALS-based estimates of biomass as a function of GLAS measurements. Three different types of predictive models are considered in CONUS to determine which produced biomass totals closest to ground-based national forest inventory estimates - (1) linear (LIN), (2) linear-no-intercept (LNI), and (3) log-linear. For CONUS at the national level, the GLAS LNI model estimate (23.95 +/- 0.45 Gt AGB), agreed most closely with the US national forest inventory ground estimate, 24.17 +/- 0.06 Gt, i.e., within 1%. The national biomass total based on linear ground-ALS and ALS-GLAS models (25.87 +/- 0.49 Gt) overestimated the national ground-based estimate by 7.5%. The comparable log-linear model result (63.29 +/-1.36 Gt) overestimated ground results by 261%. All three national biomass GLAS estimates, LIN, LNI, and log-linear, are based on 241,718 pulses collected on 230 orbits. The US national forest inventory (ground) estimates are based on 119,414 ground plots. At the US state level, the average absolute value of the deviation of LNI GLAS estimates from the comparable ground estimate of total biomass was 18.8% (range: Oregon,−40.8% to North Dakota, 128.6%). Log-linear models produced gross overestimates in the continental US, i.e., N2.6x, and the use of this model to predict regional biomass using GLAS data in temperate, western hemisphere forests is not appropriate. The best model form, LNI, is used to produce biomass estimates in Mexico. The average biomass density in Mexican forests is 53.10 +/- 0.88 t/ha, and the total biomass for the country, given a total forest area of 688,096 sq km, is 3.65 +/- 0.06 Gt. In Mexico, our GLAS biomass total underestimated a 2005 FAO estimate (4.152 Gt) by 12% and overestimated a 2007/8 radar study's figure (3.06 Gt) by 19%.

Hybrid 3-Phase sampling↗
NASA NTRS1987

Grain dynamics in zero gravity

The dynamics of granular materials has proved difficult to model, primarily because of the complications arising from inelastic losses, friction, packing, and the effect of many grains being in contact simultaneously. One interesting limit for which it was recently possible to construct a theory is that where the grain-grain interactions are dominated by binary collisions. The kinetic model of granular systems if similar to the kinetic theory of gases, except that collisional energy losses are always present in the former and must be treated explicitly. Few granular materials on Earth are describable by this limiting model, since gravity tends to collapse the grains into a high-density state where Coulombic friction effects are dominant. The planned Space Station offers an unusual opportunity to test the kinetic grain model and to explore its predictions. Without gravity, the regime of low interparticle velocities, where an elastic description of the collision is still valid, is investigated. This will allow direct interpretation by dynamical computer simulations as well as by kinetic theory.

Werner, B. T.↗
DOE OSTI2026

Structure of high-lying excited states in 47 Ca with ℓ ≥ 3 populated in fast-beam 𝛾-ray-tagged one-neutron pickup reactions

Located just one neutron removed from doubly magic, stable 48 Ca, the adjacent isotope 47 Ca has been studied extensively with neutron-removing transfer reactions and interpreted within a shell-model picture. Neutron-adding transfer reactions onto the stable 46 Ca, however, have not been performed in nearly 60 years and never with the aid of high-resolution 𝛾-ray spectroscopy. Here, we report on the 12 C⁡( 46 , Ca 47 Ca+𝛾)⁢𝑋 one-neutron pickup reaction from 12 C onto 46 Ca at more than 50 MeV/nucleon, a regime of orbital angular momentum mismatch that selectively populates 47 Ca final states based on ℓ ≥ 3 neutron configurations. Newly observed excited states and 𝛾-ray transitions are discussed in comparison to 𝑓⁡𝑝 shell-model calculations in a variety of configuration spaces. Furthermore, it is shown that the nature of the observed spectrum of final states can be well interpreted, albeit with a significant fragmentation of the 𝜈⁢0⁢𝑓 5/2 strength observed in the measurement that is not borne out in the structure calculations.

39 ≤ A ≤ 58↗
NASA NTRS1983

Interception in three dimensions - An energy formulation

The problem of minimum-time interception of a target flying in three dimensional space is analyzed with the interceptor aircraft modeled through energy-state approximation. A coordinate transformation that uncouples the interceptor's extremals from the target motion in an open-loop sense is introduced, and the necessary conditions for optimality and the optimal controls are derived. Example extremals are shown.

Rajan, N.↗
NASA NTRS1990

The role of cost analysis in manned spacecraft development

This paper describes the methods used by NASA for finding the 'best' design approach for meeting the objectives of space exploration. Models of the design process and the role that cost analysis plays in the process are described. The state-of-the-art of cost estimating is discussed, and a cost model designed for the exploration initiative is presented. The Space Exploration Initiative is discussed, and a cost estimate presented.

Cyr, Kelley J.↗
NASA NTRS1993

Linear, multivariable robust control with a mu perspective

The structured singular value is a linear algebra tool developed to study a particular class of matrix perturbation problems arising in robust feedback control of multivariable systems. These perturbations are called linear fractional, and are a natural way to model many types of uncertainty in linear systems, including state-space parameter uncertainty, multiplicative and additive unmodeled dynamics uncertainty, and coprime factor and gap metric uncertainty. The structured singular value theory provides a natural extension of classical SISO robustness measures and concepts to MIMO systems. The structured singular value analysis, coupled with approximate synthesis methods, make it possible to study the tradeoff between performance and uncertainty that occurs in all feedback systems. In MIMO systems, the complexity of the spatial interactions in the loop gains make it difficult to heuristically quantify the tradeoffs that must occur. This paper examines the role played by the structured singular value (and its computable bounds) in answering these questions, as well as its role in the general robust, multivariable control analysis and design problem.

Packard, Andy↗
NASA NTRS1993

Parametric uncertainty modeling for application to robust control

Viewgraphs and a paper on parametric uncertainty modeling for application to robust control are included. Advanced robust control system analysis and design is based on the availability of an uncertainty description which separates the uncertain system elements from the nominal system. Although this modeling structure is relatively straightforward to obtain for multiple unstructured uncertainties modeled throughout the system, it is difficult to formulate for many problems involving real parameter variations. Furthermore, it is difficult to ensure that the uncertainty model is formulated such that the dimension of the resulting model is minimal. A procedure for obtaining an uncertainty model for real uncertain parameter problems in which the uncertain parameters can be represented in a multilinear form is presented. Furthermore, the procedure is formulated such that the resulting uncertainty model is minimal (or near minimal) relative to a given state space realization of the system. The approach is demonstrated for a multivariable third-order example problem having four uncertain parameters.

Belcastro, Christine M.↗
NASA NTRS1997

Parametric State Space Structuring

Structured approaches based on Kronecker operators for the description and solution of the infinitesimal generator of a continuous-time Markov chains are receiving increasing interest. However, their main advantage, a substantial reduction in the memory requirements during the numerical solution, comes at a price. Methods based on the "potential state space" allocate a probability vector that might be much larger than actually needed. Methods based on the "actual state space", instead, have an additional logarithmic overhead. We present an approach that realizes the advantages of both methods with none of their disadvantages, by partitioning the local state spaces of each submodel. We apply our results to a model of software rendezvous, and show how they reduce memory requirements while, at the same time, improving the efficiency of the computation.

Ciardo, Gianfranco↗
NASA NTRS1992

Cosmic ray LET spectra and doses on board Cosmos-2044 biosatellite

Results of the experiments on board Cosmos-2044 (Biosatellite 9) are presented. Various nuclear track detectors (NTD) (dielectric, AgCl-based, nuclear emulsions) were used to obtain the LET spectra inside and outside the satellite. The spectra from the different NTDs have proved to be in general agreement. The results of LET spectra calculations using two different models are also presented. The resultant LET distributions are used to calculate the absorbed and equivalent doses and the orbit-averaged quality factors (QF) of the cosmic rays (CR). Absorbed dose rates inside (approximately 20 g cm-2 shielding) and outside (1 g cm-2) the spacecraft, omitting electrons, were found to be 4.8 and 8.6 mrad d-1, respectively, while the corresponding equivalent doses were 8.8 and 19.7 mrem d-1. The effects of the flight parameters on the total fluence of, and on the dose from, the CR particles are analyzed. Integral dose distributions of the detected particles are also determined. The LET values which separate absorbed and equivalent doses into 50% intervals are estimated. The CR-39 dielectric NTD is shown to detect 20-30% of the absorbed dose and 60-70% of the equivalent dose in the Cosmos-2044 orbit. The influence of solar activity phase on the magnitude of CR flux is discussed.

Non-NASA Center↗