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 595 records · Page 33

United States Forest Disturbance Trends Observed Using Landsat Time Series

Disturbance events strongly affect the composition, structure, and function of forest ecosystems; however, existing U.S. land management inventories were not designed to monitor disturbance. To begin addressing this gap, the North American Forest Dynamics (NAFD) project has examined a geographic sample of 50 Landsat satellite image time series to assess trends in forest disturbance across the conterminous United States for 1985-2005. The geographic sample design used a probability-based scheme to encompass major forest types and maximize geographic dispersion. For each sample location disturbance was identified in the Landsat series using the Vegetation Change Tracker (VCT) algorithm. The NAFD analysis indicates that, on average, 2.77 Mha/yr of forests were disturbed annually, representing 1.09%/yr of US forestland. These satellite-based national disturbance rates estimates tend to be lower than those derived from land management inventories, reflecting both methodological and definitional differences. In particular the VCT approach used with a biennial time step has limited sensitivity to low-intensity disturbances. Unlike prior satellite studies, our biennial forest disturbance rates vary by nearly a factor of two between high and low years. High western US disturbance rates were associated with active fire years and insect activity, while variability in the east is more strongly related to harvest rates in managed forests. We note that generating a geographic sample based on representing forest type and variability may be problematic since the spatial pattern of disturbance does not necessarily correlate with forest type. We also find that the prevalence of diffuse, non-stand clearing disturbance in US forests makes the application of a biennial geographic sample problematic. Future satellite-based studies of disturbance at regional and national scales should focus on wall-to-wall analyses with annual time step for improved accuracy.

Masek, Jeffrey G.↗

Combined Satellite and Surface-Based Estimation of the Intracloud/ Cloud-to-Ground Lightning Ratio Over the Continental United States

Four years of observations from the NASA Optical Transient Detector (OTD) and Global Atmospherics National Lightning Detection Network (NLDN) are combined to determine the geographic distribution of the intracloud/cloud to ground lightning ratio over the continental United States. The mean ratio over this region is 2.64-2.94, with a standard deviation of 1.1-1.3 and anomalies as low as 1.0 or less over the Rocky and Appalachian mountains and as high as 8-10 in the central-upper midwest. There is some indication that the ratio covaries with ground elevation, although the relationship is nonunique. Little evidence is found to support a latitudinal covariance, despite significant variation in the climatological mean tropopause pressure over the latitudes considered. The dynamic range of local variability is comparable to the range of values cited by previous studies for latitudinal variation from the deep tropics to midlatitudes. Local high anomalies of this ratio in the midwest are coincident with anomalies in the climatological percentage of positive CG occurrence, as well as in the occurrence of large positive CGs characteristic of MCS convective and trailing stratiform regions. This suggests that storm type, morphology and level of organization may dominate over environmental cofactors in the local determination of this ratio.

Boccippio, D. J.↗

Digital adaptive model following control

Development of an implementable digital adaptive control system requires a control law that can be efficiently adjusted during system operation and an identification algorithm that is capable of supplying parameter estimates at an accuracy and rate specified by the controller characteristics. To this effect a control law has been derived based upon the minimization of a weighted combination of the instantaneous squared error between the states of a linear plant and model at the next stage and the present control energy. The result is a control signal that is a linear combination of plant states, model states, and any model inputs. Explicit formula are available for computing the pertinent gains. To illustrate the utility of such a control law, its application to typical linearized lateral aircraft dynamics is currently being considered.

Alag, G.↗

A methodology for designing robust multivariable nonlinear feedback control systems

This paper presents an outline of a methodology for the design of nonlinear dynamic compensators for nonlinear multivariable systems with guarantees of closed-loop stability, robustness, and performance. The method is an extension of the Linear-Quadratic-Gaussian with Loop-Transfer-Recovery (LQG/LTR) methodology for linear systems, thus hinging upon the idea of constructing an approximate nonlinear inverse operator for the plant. A major feature of the method is a unification of both the state-space and Input-Output formulations. The major results of the paper are: (1) Properties of the solution of certain optimal control problems; (2) New results on nonlinear estimation, including a new nonlinear observer with certain guaranteed properties; (3) The Loop Operator Recovery (LOR) procedure, and (4) Inclusion of command-following and disturbance rejection with a goal of a closed-loop stable, robust system.

Grunberg, D. B.↗

Adjoints and Low-rank Covariance Representation

Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.

Tippett, Michael K.↗

HLPW-4/GMGW-3: Mesh Adaptation for RANS Technology Focus Group Workshop Summary

The AIAA workshop series helps guide the Computational Fluid Dynamics (CFD) community through publicly available test cases and impartial evaluation to establish the state-of-the-art in applied CFD. A summary of the Mesh Adaptation Technical Focus Group submissions for the joint Fourth CFD High Lift Prediction and Third Geometry & Mesh Generation Workshop is provided. The Common Research Model High-Lift variant (CRM-HL) is the target of this investigation. A 2D high-lift airfoil solution verification case was extracted from the 3D CRM-HL. Low variation is shown between 2D solutions with mesh adaptation to control estimated solution interpolation and output errors. As expected, variation for these mesh adaptation techniques is larger for the 3D CRM-HL than the 2D airfoil. More importantly, the 3D mesh adapted results have a lower variation than expert-crafted meshes of equal or larger size. Suspected multiple solutions to Reynolds-averaged Navier-Stokes with the Spalart-Allmaras turbulence model appear to resolve into a common solution with mesh resolution.

HLPW-4↗

Supersonic Stall Flutter of High Speed Fans

An analytical model is developed for predicting the onset of supersonic stall bending flutter in axial flow compressors. The analysis is based on a modified two dimensional, compressible, unsteady actuator disk theory. It is applied to a rotor blade row by considering a cascade of airfoils whose geometry and dynamic response coincide with those of a rotor blade element at 85 percent of the span height (measured from the hub). The rotor blades are assumed to be unshrouded (i.e., free standing) and to vibrate in their first flexural mode. The effects of shock waves and flow separation are included in the model through quasi-steady, empirical, rotor total-pressure-loss and deviation-angle correlations. The actuator disk model predicts the unsteady aerodynamic force acting on the cascade blading as a function of the steady flow field entering the cascade and the geometry and dynamic response of the cascade. Calculations show that the present model predicts the existence of a bending flutter mode at supersonic inlet Mach numbers. This flutter mode is suppressed by increasing the reduced frequency of the system or by reducing the steady state aerodynamic loading on the cascade. The validity of the model for predicting flutter is demonstrated by correlating the measured flutter boundary of a high speed fan stage with its predicted boundary. This correlation uses a level of damping for the blade row (i.e., the log decrement of the rotor system) that is estimated from the experimental flutter data. The predicted flutter boundary is shown to be in good agreement with the measured boundary.

Adamczyk, J. J.↗

Identification of helicopter rotor dynamic models

A recursive, extended Kalman-filter approach is applied to the identifiction of rotor damping levels of representative helicopter dynamic systems. The general formulation of the approach is presented in the context of a typically posed stochastic estimation problem, and the method is analytically applied to determining the damping levels of a coupled rotor-body system. The identified damping covergence characteristics are studied for sensitivity to both constant-coefficient and periodic-coefficient measurement models, process-noise covariance levels, and specified initial estimates of the rotor-system damping. A second application of the method to identifying the plant model for a highly damped, isolated flapping blade with a constant-coefficient state model (hover) and a periodic-coefficient state model (forward flight) is also investigated. The parameter-identification capability is evaluated for the effect of periodicity on the plant model coefficients and the influence of different measurement noise levels.

Molusis, J. A.↗

James Webb Space Telescope Navigation Optimization Challenges

The James Webb Space Telescope (JWST) is a NASA flagship mission launched on December 25, 2021. The early orbit phase was highlighted by three midcourse correction burns that were performed to maneuver the vehicle into a libration orbit at the second Sun-Earth-Moon (SEM) libration point (L2). During the coast to L2, the vehicle’s primary observatory and sunshield were deployed. Once the tennis-court-sized sunshield was unfurled and tensioned into place, the large area exposed to solar radiation pressure (SRP) dictated that the SRP force model would need to account for the vehicle’s geometry, reflective properties, and orientation relative to the Sun. Following the insertion of JWST into its L2 libration orbit on January 24, 2022, the commissioning phase for the vehicle’s observatory commenced along with the first cycle of routine station-keeping maneuvers. The NASA Goddard Space Flight Center’s Flight Dynamics Facility (FDF) provides navigation services for the JWST mission. The FDF serves as the prime or backup navigation operations center for more than 30 active missions spanning a wide array of flight regimes performing orbit determination, maneuver planning, trajectory optimization, and tracking data evaluation. Definitive orbit determination for JWST is performed by the FDF using an Extended Kalman Filter (EKF) to estimate the vehicle’s trajectory and reflective properties using Deep Space Network (DSN) TRK-2-34 tracking data and spacecraft attitude telemetry. For the purposes of orbit prediction and station-keeping maneuver targeting, the sensitivity of the SRP force to the vehicle’s orientation requires that predictive attitude information must be incorporated to generate accurate orbit predictions and maneuver plans. Predictive attitude plans are provided to the FDF by the JWST Spacecraft Operations Center (SOC). Short-term (28 day) orbit predictions are propagated using a Short Range Attitude Plan (SRAP) to model the vehicle’s future attitude states for up to a week. SRAP attitude data is highly reliable, as it reflects finalized planned attitude states that the vehicle will be commanded to attain during the ensuing week. A Long Range Attitude Plan (LRAP) can be implemented to model predicted attitude states beyond one week, but these predictions are tentative and subject to revision due to changes in science operations plans. Instead, a conservative approach of applying a Sun-Pointing Neutral (SPN) attitude configuration is utilized for long-term (2 year) orbit predictions where reliable planned attitude data is not available. SPN attitude mode aligns the net SRP force along the JWST-to-Sun vector where it becomes independent of the vehicle’s angle of rotation about this vector, defined as the Sun yaw. Given the dynamic instability of the libration orbit, station-keeping thrust must be applied in either the sunward or anti-sunward direction to balance the resulting orbit. The attitude constraints of the vehicle also impose limits on the available pointing directions for station-keeping thrusters. The thrusters cannot be aligned with the optimal pointing direction for sunward maneuvers, rendering sunward maneuvers to be significantly less fuel efficient than anti-sunward maneuvers. Consequently, station-keeping maneuvers must be targeted to balance the orbit while ensuring that the next maneuver will also be performed in the anti-sunward direction to optimize propellant usage. When targeting a station-keeping maneuver, the attitude and SRP modeling configuration which is applied to the predicted post-maneuver trajectory will dictate the direction of the subsequent station-keeping maneuver, assuming nominal propulsion system performance. If the predicted post-maneuver attitude states result in an SRP model which under-predicts the cumulative SRP impact on the post-maneuver orbit, a targeted anti-sunward maneuver will be larger than necessary and the next maneuver will need to be executed sunward in order to compensate. In contrast, an anti-sunward maneuver which is targeted using an SRP model which over-predicts the cumulative SRP impact will achieve station-keeping while ensuring that the next maneuver will likewise be performed anti-sunward. For this reason, station-keeping maneuvers are targeted while applying the SPN attitude mode, as this mode entails the largest possible SRP area cross-section and therefore a larger modeled cumulative SRP impact. This paper documents the NASA Goddard Space Flight Center’s FDF support for JWST on-orbit operations and the analysis projects undertaken utilizing the experiences and data accumulated throughout the first full year of routine science operations. The results of these analysis efforts have been used to implement improvements to orbit prediction accuracy and maneuver efficiency which have the potential to prolong the lifespan of JWST to continue to conduct ground-breaking infrared astronomy for decades to come.

Flight Dynamics↗

Geometrical evidence for dark matter: X-ray constraints on the mass of the elliptical galaxy NGC 720

We describe (1) a new test for dark matter and alternate theories of gravitation based on the relative geometries of the X-ray and optical surface brightness distributions and an assumed form for the potential, of the optical light, (2) a technique to measure the shapes of the total gravitating matter and dark matter of an ellipsoidal system which is insensitive to the precise value of the temperature of the gas and to modest temperature gradients, and (3) a new method to determine the ratio of dark mass to stellar mass that is dependent on the functional forms for the visible star, gas and dark matter distributions, but independent of the distance to the galaxy or the gas temperature. We apply these techniques to X-ray data from the ROSAT Position Sensitive Proportional Counter (PSPC) of the optically flattened elliptical galaxy NGC 720; the optical isophotes have ellipticity epsilon approximately 0.40 extending out to approximately 120 sec. The X-ray isophotes are significantly elongated, epsilon = 0.20-0.30 for semimajor axis a approximately 100 sec. The major axes of the optical and X-ray isophotes are misaligned by approximately 30 deg +/- 15 deg. Spectral analysis of the X-ray data reveals no evidence of temperature gradients or anisotropies and demonstrates that a single-temperature plasma (T approximately 0.6 keV) having subsolar heavy element abundances and a two-temperature model having solar abundances describe the spectrum equally well. Considering only the relative geometries of the X-ray and optical surface brightness distributions and an assumed functional form for the potential of the optical light, we conclude that matter distributed like the optical light cannot produce the observed ellipticities of the X-ray isophotes, independent of the gas pressure, the gas temperature, and the value of the stellar mass; this comparison assumes a state of quasi-hydrostatic equilibrium so that the three-dimensional surfaces of the gas emissivity trace the three-dimensional isopotential surfaces -- we discuss the viability of this assumption for NGC 720. Milgrom's Modification of Newtonian Dynamics (MOND) cannot dispel this manifestation of dark matter. Hence, geometrical considerations require, without mention of pressure or temperature, the presence of an extended, massive dark matter halo in NGC 720. Employing essentially the technique of Buote & Canizares (1992; Buote 1992) we use the shape of the X-ray surface brightness to constrain the shape of the total gravitating matter. The total matter is modeled as either an oblate or prolate spheriod of constant shape and orientation having either a Ferrers (rho approximately r(exp -n)) or Hernquist density. Assuming the X-ray gas is in hydrostatic equilibrium, we construct a model X-ray gas distribution for various temperature profiles. We determine the ellipticity of the total gravitating matter to be epsilon approximately 0.50-0.70. Using the single-temperature model we estimate a total mass approximately (0.41-1.4) x 10(exp 12) h(sub 80) solar mass interior to the ellipsoid of semimajor axis 43.6 h(sub 80) kpc. Ferrers densities as steep as r(exp -3) do not fit the data, but the r(exp -2) and Hernquist models yield excellent fits. We estimate the mass distributions of the stars and the gas and fit the dark matter directly. For a given gas equation of state and functional forms for the visible stars, gas, and dark matter, these models yield a distance-independent and temperature-independent measurement of the ratio of dark mass to stellar mass M(sub DM)/M(sub stars). We estimate a minimum M(sub DM)/M(sub stars) greater than or equal to 4 which corresponds to a total mass slightly greater than that derived from the single-temperature models for distance D = 20h(sub 80) Mpc.

Buote, David A.↗

Effects of surface baroclinicity on frontal overrunning along the central Gulf coast

Along the United States Gulf coast and over the northern Gulf of Mexico, frontal overrunning occurs frequently. Cyclogenesis over the Gulf is often associated with this type of weather system. Effects of baroclinic fields on frontal overrunning are investigated from synoptic and climatological points of view. It is found that, from October through April, the orientation of the shelf break is a very important baroclinic characteristic because fronts tend to stall there rather than at the physical coastline. To further substantiate this deduction dynamically, the local geostrophic vorticity field over the western Louisiana-upper Texas shelf region is estimated monthly. The correlation coefficient between the vorticity field and the frequency of frontal overrunning along the central Gulf coast was 0.86. For forecasting applications, a simple formula is provided to estimate this local vorticity from the temperature difference between Lake Charles, Louisiana, and buoy station 42002 in the deep Gulf.

Hsu, S. A.↗

Initial Evaluations of LoC Prediction Algorithms Using the NASA Vertical Motion Simulator

Flying near the edge of the safe operating envelope is an inherently unsafe proposition. Edge of the envelope here implies that small changes or disturbances in system state or system dynamics can take the system out of the safe envelope in a short time and could result in loss-of-control events. This study evaluated approaches to predicting loss-of-control safety margins as the aircraft gets closer to the edge of the safe operating envelope. The goal of the approach is to provide the pilot aural, visual, and tactile cues focused on maintaining the pilot's control action within predicted loss-of-control boundaries. Our predictive architecture combines quantitative loss-of-control boundaries, an adaptive prediction method to estimate in real-time Markov model parameters and associated stability margins, and a real-time data-based predictive control margins estimation algorithm. The combined architecture is applied to a nonlinear transport class aircraft. Evaluations of various feedback cues using both test and commercial pilots in the NASA Ames Vertical Motion-base Simulator (VMS) were conducted in the summer of 2013. The paper presents results of this evaluation focused on effectiveness of these approaches and the cues in preventing the pilots from entering a loss-of-control event.

Pilot Cue↗

Dynamical response of the dayside ionosphere of Venus to the solar wind

Dayside ion composition measurements made by the orbiter ion mass spectrometer and the orbiter electron temperature probe on the Pioneer Venus orbiter are used to infer the dominant processes involved in the dynamic response of the Venus ionosphere to the solar wind. The analysis is confined to the topside ionosphere in the vicinity of the subsolar point, where the ionosphere-solar wind interaction is expected to be maximized. Height profiles of the ion composition and plasma temperatures in the main body of the topside ionosphere, lying between the ionopause and chemical equilibrium regions, reveal that the ionosphere is in a compressed state. This region of the ionosphere is interpreted in terms of a stationary equilibrium where the compression is derived from the ponderomotive force j x B. The estimated magnitude of this force is confirmed by the magnetic field measurements made by the orbiter magnetometer.

Hartle, R. E.↗

Resonant enhanced multiphoton ionization studies of atomic oxygen

In resonant enhanced multiphoton ionization (REMPI), an atom absorbs several photons making a transition to a resonant intermediate state and subsequently ionizing out of it. With currently available tunable narrow-band lasers, the extreme sensitivity of REMPI to the specific arrangement of levels can be used to selectively probe minute amounts of a single species (atom) in a host of background material. Determination of the number density of atoms from the observed REMPI signal requires a knowledge of the multiphoton ionization cross sections. The REMPI of atomic oxygen was investigated through various excitation schemes that are feasible with available light sources. Using quantum defect theory (QDT) to estimate the various atomic parameters, the REMPI dynamics in atomic oxygen were studied incorporating the effects of saturation and a.c. Stark shifts. Results are presented for REMPI probabilities for excitation through various 2p(3) (4S sup o) np(3)P and 2p(3) (4S sup o) nf(3)F levels.

Dixit, S. N.↗

Observability and Estimation of Distributed Space Systems via Local Information-Exchange Networks

In this work, we develop an approach to formation estimation by explicitly characterizing formation's system-theoretic attributes in terms of the underlying inter-spacecraft information-exchange network. In particular, we approach the formation observer/estimator design by relaxing the accessibility to the global state information by a centralized observer/estimator- and in turn- providing an analysis and synthesis framework for formation observers/estimators that rely on local measurements. The noveltyof our approach hinges upon the explicit examination of the underlying distributed spacecraft network in the realm of guidance, navigation, and control algorithmic analysis and design. The overarching goal of our general research program, some of whose results are reported in this paper, is the development of distributed spacecraft estimation algorithms that are scalable, modular, and robust to variations inthe topology and link characteristics of the formation information exchange network. In this work, we consider the observability of a spacecraft formation from a single observation node and utilize the agreement protocol as a mechanism for observing formation states from local measurements. Specifically, we show how the symmetry structure of the network, characterized in terms of its automorphism group, directly relates to the observability of the corresponding multi-agent system The ramification of this notion of observability over networks is then explored in the context of distributed formation estimation.

observability↗

A novel multistage estimation of the signal parameters of a possibly data-modulated sinusoid under very high dynamics

A novel multistage estimation scheme is presented for estimating the parameters of a received carrier signal possibly phase-modulated by unknown data, and experiencing very high Doppler, Doppler rate, etc. Such a situation arises, for example, in the case of Global Positioning Systems (GPS) where the signal parameters are directly related to the position, velocity, acceleration and jerk of the GPS receiver. In the proposed multistage scheme, the first stage estimator operates as a coarse estimator resulting in higher rms estimation errors but with a relatively small probability of the frequency estimation error exceeding one-half of the sampling frequency (an event termed cycle slip). The second state of the estimator operates on the error signal available from the first stage, refining the overall estimates, and in the process also reduces the number of cycle slips. The first stage algorithm is selected to be a modified least squares algorithm operating upon the differential signal model and referred to as differential least squares. This estimation stage provides relatively coarse estimates of the frequency and its derivatives. The second algorithm is simply an extended Kalman filter which also yields the estimate of the phase along with a more refined estimate of frequency as well. A major advantage of the proposed algorithm is a reduction in the threshold on received carrier power-to-noise power spectral density ratio as compared to the threshold achievable by either of these algorithms alone. In fact, it appears from the simulations that for the case of an unmodulated carrier, the proposed scheme achieves the same threshold as for an almost exact and computationally intensive implementation of the maximum likelihood estimator.

Kumar, Rajendra↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

Second coefficient of viscosity in air

Acoustic attenuation measurements in air were analyzed in order to estimate the second coefficient of viscosity. Data over a temperature range of 11 C to 50 C and at relative humidities between 6 percent and 91 percent were used. This analysis showed that the second coefficient of viscosity varied between 1900 and 20,000 times larger than the dynamic or first coefficient of viscosity over the temperature and humidity range of the data. In addition, the data showed that the molecular relaxation effects, which are responsible for the magnitude of the second coefficient of viscosity, place severe limits on the use of time-independent, thermodynamic equations of state. Compressible flows containing large streamwise velocity gradients, like shock waves, which cause significant changes in particle properties to occur during time intervals shorter than hundredths of seconds, must be modeled using dynamic equations of state. The dynamic model approach is described briefly.

Ash, Robert L.↗