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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Quantum utility-scale error mitigation for quantum quench dynamics in Heisenberg spin chains

Here, we implement a quantum error mitigation method termed self-mitigation, which is comparable to zero-noise extrapolation, at large scales to achieve quantum utility on near-term, noisy quantum computers. We investigate the effectiveness of several quantum error mitigation strategies, including self-mitigation, by simulating quantum quench dynamics for Heisenberg spin chains with system sizes up to 104 qubits using IBM quantum processors. In particular, we discuss the limitations of zero-noise extrapolation and the advantages offered by self-mitigation at large scales. The self-mitigation method demonstrates stable accuracy with large systems of 104 qubits comprising more than 3,000 CNOT gates. Also, we combine the discussed quantum error mitigation methods with practical entanglement entropy measuring methods, and it shows a good agreement with the theoretical estimation. Our study illustrates the usefulness of near-term noisy quantum hardware in examining the quantum quench dynamics of many-body systems at large scales and lays the groundwork for surpassing classical simulations with quantum methods prior to the development of fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sensitive detection of structural dynamics using a statistical framework for comparative crystallography

Chemical and conformational changes are crucial to protein function and its pharmacological control. X-ray crystallography can reveal these changes in atomic detail, but standard analysis methods, which refine separate datasets, often overlook differences that are subtle or arise in only a subset of molecules. Direct comparison of crystallographic datasets is, in principle, more powerful, but systematic errors (“scales”) often mask changes in the crystallographic observables (“structure factors”). Machine learning algorithms that jointly estimate scales and structure factors can address this limitation. Here, we augment this approach with multivariate, structured priors derived from crystallographic theory, implemented in the variational deep learning framework Careless. Doing so strongly improves the detection of protein dynamics, element-specific anomalous signals, and the binding of drug candidates, offering a robust approach to comparative crystallography and, potentially, to detection of protein dynamics by other structure determination methods.

Hekstra, Doeke R. [Harvard Univ., Cambridge, MA (U↗
DOE OSTI2024

Final Report on Predictive Analyses of PRD as Function of Anomalies

This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Recent Advances toward Efficient Calculation of Higher Nuclear Derivatives in Quantum Chemistry

In this article, we provide an overview of state-of-the-art techniques that are being developed for efficient calculation of second and higher nuclear derivatives of quantum mechanical (QM) energy. Calculations of nuclear Hessians and anharmonic terms incur high costs and memory and scale poorly with system size. Three emerging classes of methods—machine learning (ML), automatic differentiation (AD), and matrix completion (MC)—have demonstrated promise in overcoming these challenges. We illustrate studies that employ unsupervised ML methods to reduce the need for multiple Hessian calculations in dynamics simulations and those that utilize supervised ML to construct approximate potential energy surfaces and estimate Hessians and anharmonic terms at reduced cost. By extension, if electronic structure operations could be written in a manner similar to functions underlying ML methods, rapid differentiation or AD routines can be employed to inexpensively calculate higher arbitrary-order derivatives. While ML approaches are typically black-box, we describe methods such as compressed sensing (CS) and MC, which explicitly leverage problem-specific mathematical properties of higher derivatives such as sparsity and low-rank, to complete higher derivative information using only a small, incomplete sample. The three classes of methods facilitate reliable predictions of observables ranging from infrared spectra to thermal conductivity and constitute a promising way forward in accurately capturing otherwise intractable higher-order responses of QM energy to nuclear perturbations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1981

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.↗
NASA NTRS1983

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.↗
NASA NTRSDate not supplied

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↗
NASA NTRS1994

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.↗
DOE OSTI2021

Chiral Effective Field Theory’s Impact on Advancing Quantum Monte Carlo Methods

Thirty years ago, Steven Weinberg published his seminal paper on “Nuclear Forces from chiral Lagrangians” which has revolutionized the field of theoretical nuclear physics. Nowadays, interactions derived from chiral effective field theory are routinely used to describe nuclear systems ranging from atomic nuclei to the dense matter explored in the core of neutron stars with theoretical uncertainty estimates. Here, in our contribution to the special issue “Celebrating 30 years of Steven Weinberg’s paper Nuclear Forces from Chiral Lagrangians”, we focus on the impact that chiral effective field theory interactions have played in advancing microscopic studies of atomic nuclei and the nuclear-matter equation of state using quantum Monte Carlo methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Long-lived period-doubled edge modes of interacting and disorder-free Floquet spin chains

Floquet spin chains have been a venue for understanding topological states of matter that are qualitatively different from their static counterparts by, for example, hosting π edge modes that show stable period-doubled dynamics. However the stability of these edge modes to interactions has traditionally required the system to be many-body localized in order to suppress heating. In contrast, here we show that even in the absence of disorder, and in the presence of bulk heating, π edge modes are long lived. Their lifetime is extracted from exact diagonalization and is found to be non-perturbative in the interaction strength. A tunneling estimate for the lifetime is obtained by mapping the stroboscopic time-evolution to dynamics of a single particle in Krylov subspace. In this subspace, the π edge mode manifests as the quasi-stable edge mode of an inhomogeneous Su-Schrieffer-Heeger model whose dimerization vanishes in the bulk of the Krylov chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

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.↗
NASA NTRS2014

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↗
DOE OSTI2023

Quantum spin ice response to a magnetic field in the dipole-octupole pyrochlore Ce 2 Zr 2 O 7

The pyrochlore magnet Ce 2 Zr 2 O 7 has attracted much attention as a quantum spin ice candidate whose novelty derives in part from the dipolar-octupolar nature of the Ce 3+ pseudospin-1/2 degrees of freedom it possesses. We report heat capacity measurements on single crystal samples of Ce 2 Zr 2 O 7 down to T~0.1K in a magnetic field along the [$1,\bar{1},0$] direction. These measurements show that the broad hump in the zero-field heat capacity moves higher in temperature with increasing field strength and is split into two separate humps by the [$1,\bar{1},0$] magnetic field at ~2T. These separate features are due to the decomposition of the pyrochlore lattice into effectively decoupled chains for fields in this direction: One set of chains (α chains) is polarized by the field while the other (β chains) remains free. This situation is similar to that observed in the classical spin ices Ho 2 Ti 2 O 7 and Dy 2 Ti 2 O 7 , but with the twist that here the strong transverse exchange interactions produce substantial quantum effects. Our theoretical modeling suggests that the β chains are close to a critical state, with nearly-gapless excitations. We also report elastic and inelastic neutron scattering measurements on single crystal Ce 2 Zr 2 O 7 in [$1,\bar{1},0$] and [0,0,1] magnetic fields at temperatures down to T=0.03K. The elastic scattering behaves consistently with the formation of independent chains for a [$1,\bar{1},0$] field, while the [0,0,1] field produces a single field-induced elastic magnetic Bragg peak at (0,2,0) and equivalent wavevectors, indicating a polarized spin ice state for fields above ~3T. For both [$1,\bar{1},0$] and [0,0,1] magnetic fields, our inelastic neutron scattering results show an approximately dispersionless continuum of scattering that increases in both energy and intensity with increasing field strength. By modeling the complete set of experimental data using numerical linked cluster and semiclassical molecular dynamics calculations, we demonstrate the dominantly multipolar nature of the exchange interactions in Ce 2 Zr 2 O 7 and the smallness of the parameter θ, which controls the mixing between dipolar and octupolar degrees of freedom. Finally, these results support previous estimates of the microscopic exchange parameters and place strong constraints on the theoretical description of this prominent spin ice candidate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

The Role of Internal Variability in Springtime Arctic Amplification from 1980 to 2022

Arctic amplification (AA) refers to the enhanced warming of the Arctic relative to the global average due to rising greenhouse gases, measured as the ratio of Arctic-mean to global-mean surface air temperature (SAT) trends. From 1980 to 2022, annual-mean AA reached 4.2 (Arctic defined as north of 70°N). Climate models simulate AA but fail to reproduce its magnitude. Sweeney et al. attributed much of this model–observation discrepancy to internal variability. AA shows seasonality and so does the discrepancy. Spring (March–May) shows the largest gap: Observed AA is 4.2, while the multimodel mean is 2.7. This raises several questions: 1) What role does internal variability play in observed spring AA? 2) How does simulated spring AA compare to observations when internal variability is removed? 3) If internal variability is significant, what mechanisms drive it? To address these, we adapted the machine learning algorithm from Sweeney et al., training on simulated multidecadal spring SAT and sea level pressure (SLP) trend maps. Our results show that internal variability enhanced spring Arctic warming by 37% and reduced global warming by 10%. Removing internal variability reconciles the spring AA discrepancy. The estimated internal contribution to Arctic spring warming is supported by an independent dynamical adjustment approach. We identify an atmospheric circulation pattern in observations associated with this internal warming. Observed internal Siberian SAT and SLP trends follow the simulated SAT–SLP relationship but lie at the distribution’s extreme, suggesting models generally underestimate internal variability unless the observed configuration reflects a rare real-world realization.

Arctic↗
NASA NTRS1980

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.↗
NASA NTRS1987

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.↗
NASA NTRS2008

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↗
NASA NTRS1989

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↗