Search NASA⌕ Search

Engineering topics

Carpenter, J. Russell

Publications and source records attributed to Carpenter, J. Russell.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Covariance Realism Is Not Enough

A great deal of effort has been put into improving the practice of space situational awareness such that covariance data associated with predicted close approaches is more “realistic.” However, “realistic” usually has meant “larger” and this presents a problem. In many cases, there exist multiple sources for predictive ephemerides, which may be fused to produce predictive states with smaller associated covariances. Ancillary to the fusion computation is the capability to assess consistency of the estimates. If actionable covariance information becomes available, interval estimates for the miss distance provide a more informative alternative to collision probability for risk assessment.

Carpenter, J. Russell↗
NASA NTRS2018

Satellite Formation-Flying and Rendezvous

GNSS has come to play an increasingly important role in satellite formation-flying and rendezvous applications. In the last decades, the use of GNSS measurements has provided the primary technique for determining the relative position of cooperative co-orbiting satellites in low Earth orbit.

D'Amico, Simone↗
NASA NTRS2017

Relevance of the American Statistical Society's Warning on P-Values for Conjunction Assessment

On March 7, 2016, the American Statistical Association (ASA) issued an online editorial paper on the "context, process, and purpose of p-values." According to the paper, "the statement articulates in non-technical terms a few select principles that could improve the conductor interpretation of quantitative science, according to widespread consensus in the statistical community." This "consensus" statement was accompanied by 21 "commentaries" expressing a diversity of opinions among the panel of experts ASA convened. In the present work, we express our view that the ASA p-value warning has relevance to the space object Conjunction Assessment (CA) community.

p-values↗
NASA NTRS2017

Precomputing Process Noise Covariance for Onboard Sequential Filters

Process noise is often used in estimation filters to account for unmodeled and mismodeled accelerations in the dynamics. The process noise covariance acts to inflate the state covariance over propagation intervals, increasing the uncertainty in the state. In scenarios where the acceleration errors change significantly over time, the standard process noise covariance approach can fail to provide effective representation of the state and its uncertainty. Consider covariance analysis techniques provide a method to precompute a process noise covariance profile along a reference trajectory using known model parameter uncertainties. The process noise covariance profile allows significantly improved state estimation and uncertainty representation over the traditional formulation. As a result, estimation performance on par with the consider filter is achieved for trajectories near the reference trajectory without the additional computational cost of the consider filter. The new formulation also has the potential to significantly reduce the trial-and-error tuning currently required of navigation analysts. A linear estimation problem as described in several previous consider covariance analysis studies is used to demonstrate the effectiveness of the precomputed process noise covariance, as well as a nonlinear descent scenario at the asteroid Bennu with optical navigation.

onboard↗
NASA NTRS2016

A Collision Avoidance Strategy for a Potential Natural Satellite around the Asteroid Bennu for the OSIRIS-REx Mission

The cadence of proximity operations for the OSIRIS-REx mission may have an extra induced challenge given the potential of the detection of a natural satellite orbiting the asteroid Bennu. Current ground radar observations for object detection orbiting Bennu show no found objects within bounds of specific size and rotation rates. If a natural satellite is detected during approach, a different proximity operation cadence will need to be implemented as well as a collision avoidance strategy for mission success. A collision avoidance strategy will be analyzed using the Wald Sequential Probability Ratio Test.

satellite↗
NASA NTRS2016

Global Positioning System Navigation Above 76,000 km for NASA's Magnetospheric Multiscale Mission

NASA's Magnetospheric Multiscale (MMS) mission, launched in March of 2015, consists of a controlled formation of four spin-stabilized spacecraft in similar highly elliptic orbits reaching apogee at radial distances of 12 and 25 Earth radii (RE) in the first and second phases of the mission. Navigation for MMS is achieved independently on-board each spacecraft by processing Global Positioning System (GPS) observables using NASA Goddard Space Flight Center (GSFC)'s Navigator GPS receiver and the Goddard Enhanced Onboard Navigation System (GEONS) extended Kalman filter software. To our knowledge, MMS constitutes, by far, the highest-altitude operational use of GPS to date and represents a high point of over a decade of high-altitude GPS navigation research and development at GSFC. In this paper we will briefly describe past and ongoing high-altitude GPS research efforts at NASA GSFC and elsewhere, provide details on the design of the MMS GPS navigation system, and present on-orbit performance data from the first phase. We extrapolate these results to predict performance in the second phase orbit, and conclude with a discussion of the implications of the MMS results for future high-altitude GPS navigation, which we believe to be broad and far-reaching.

GPS↗
NASA NTRS2015

Cauchy Drag Estimation For Low Earth Orbiters

Recent work on minimum variances estimators based on Cauchy distributions appear relevant to orbital drag estimation. Samples form Cauchy distributions which are part of a class of heavy-tailed distributions, are characterized by long stretches of fairly small variation, punctuated by large variations that are many times larger than could be expected from a Gaussian. Such behavior can occur when solar storms perturb the atmosphere. In this context, the present work describes an embedding of the scalar Idan-Speyer Cauchy Estimator to estimate density corrections, within an Extended Kalman Filter that estimates the state of a low Earth orbiter. In contrast to the baseline Kalman approach, the larger formal errors of the present approach fully and conservatively bound the predictive error distribution, even in the face of unanticipated density disturbances of hundreds of percent.

Low Earth↗
NASA NTRS2014

Orbit Determination Accuracy Analysis of the Magnetospheric Multiscale Mission During Perigee Raise

The Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF) will provide orbit determination and prediction support for the Magnetospheric Multiscale (MMS) mission during the mission's commissioning period. The spacecraft will launch into a highly elliptical Earth orbit in 2015. Starting approximately four days after launch, a series of five large perigee-raising maneuvers will be executed near apogee on a nearly every-other-orbit cadence. This perigee-raise operations concept requires a high-accuracy estimate of the orbital state within one orbit following the maneuver for performance evaluation and a high-accuracy orbit prediction to correctly plan and execute the next maneuver in the sequence. During early mission design, a linear covariance analysis method was used to study orbit determination and prediction accuracy for this perigee-raising campaign. This paper provides a higher fidelity Monte Carlo analysis using the operational COTS extended Kalman filter implementation that was performed to validate the linear covariance analysis estimates and to better characterize orbit determination performance for actively maneuvering spacecraft in a highly elliptical orbit. The study finds that the COTS extended Kalman filter tool converges on accurate definitive orbit solutions quickly, but prediction accuracy through orbits with very low altitude perigees is degraded by the unpredictability of atmospheric density variation.

Multiscale↗
NASA NTRS2014

Linear Covariance Analysis and Epoch State Estimators

This paper extends in two directions the results of prior work on generalized linear covariance analysis of both batch least-squares and sequential estimators. The first is an improved treatment of process noise in the batch, or epoch state, estimator with an epoch time that may be later than some or all of the measurements in the batch. The second is to account for process noise in specifying the gains in the epoch state estimator. We establish the conditions under which the latter estimator is equivalent to the Kalman filter.

Linear covariance analysis↗
NASA NTRS2014

Improving FermiI Orbit Determination and Prediction in an Uncertain Atmospheric Drag Environment

Orbit determination and prediction of the Fermi Gamma-ray Space Telescope trajectory is strongly impacted by the unpredictability and variability of atmospheric density and the spacecrafts ballistic coefficient. Operationally, Global Positioning System point solutions are processed with an extended Kalman filter for orbit determination, and predictions are generated for conjunction assessment with secondary objects. When these predictions are compared to Joint Space Operations Center radar-based solutions, the close approach distance between the two predictions can greatly differ ahead of the conjunction. This work explores strategies for improving prediction accuracy and helps to explain the prediction disparities. Namely, a tuning analysis is performed to determine atmospheric drag modeling and filter parameters that can improve orbit determination as well as prediction accuracy. A 45 improvement in three-day prediction accuracy is realized by tuning the ballistic coefficient and atmospheric density stochastic models, measurement frequency, and other modeling and filter parameters.

orbit prediction↗
NASA NTRS2014

Improving Fermi Orbit Determination and Prediction in an Uncertain Atmospheric Drag Environment

Orbit determination and prediction of the Fermi Gamma-ray Space Telescope trajectory is strongly impacted by the unpredictability and variability of atmospheric density and the spacecraft's ballistic coefficient. Operationally, Global Positioning System point solutions are processed with an extended Kalman filter for orbit determination, and predictions are generated for conjunction assessment with secondary objects. When these predictions are compared to Joint Space Operations Center radar-based solutions, the close approach distance between the two predictions can greatly differ ahead of the conjunction. This work explores strategies for improving prediction accuracy and helps to explain the prediction disparities. Namely, a tuning analysis is performed to determine atmospheric drag modeling and filter parameters that can improve orbit determination as well as prediction accuracy. A 45% improvement in three-day prediction accuracy is realized by tuning the ballistic coefficient and atmospheric density stochastic models, measurement frequency, and other modeling and filter parameters.

orbit prediction↗
NASA NTRS2013

Wald Sequential Probability Ratio Test for Analysis of Orbital Conjunction Data

We propose a Wald Sequential Probability Ratio Test for analysis of commonly available predictions associated with spacecraft conjunctions. Such predictions generally consist of a relative state and relative state error covariance at the time of closest approach, under the assumption that prediction errors are Gaussian. We show that under these circumstances, the likelihood ratio of the Wald test reduces to an especially simple form, involving the current best estimate of collision probability, and a similar estimate of collision probability that is based on prior assumptions about the likelihood of collision.

Carpenter, J. Russell↗
NASA NTRS2013

Sequential Probability Ratio Test for Spacecraft Collision Avoidance Maneuver Decisions

A document discusses sequential probability ratio tests that explicitly allow decision-makers to incorporate false alarm and missed detection risks, and are potentially less sensitive to modeling errors than a procedure that relies solely on a probability of collision threshold. Recent work on constrained Kalman filtering has suggested an approach to formulating such a test for collision avoidance maneuver decisions: a filter bank with two norm-inequality-constrained epoch-state extended Kalman filters. One filter models the null hypotheses that the miss distance is inside the combined hard body radius at the predicted time of closest approach, and one filter models the alternative hypothesis. The epoch-state filter developed for this method explicitly accounts for any process noise present in the system. The method appears to work well using a realistic example based on an upcoming, highly elliptical orbit formation flying mission.

Carpenter, J. Russell↗
NASA NTRS2012

Current-State Constrained Filter Bank for Wald Testing of Spacecraft Conjunctions

We propose a filter bank consisting of an ordinary current-state extended Kalman filter, and two similar but constrained filters: one is constrained by a null hypothesis that the miss distance between two conjuncting spacecraft is inside their combined hard body radius at the predicted time of closest approach, and one is constrained by an alternative complementary hypothesis. The unconstrained filter is the basis of an initial screening for close approaches of interest. Once the initial screening detects a possibly risky conjunction, the unconstrained filter also governs measurement editing for all three filters, and predicts the time of closest approach. The constrained filters operate only when conjunctions of interest occur. The computed likelihoods of the innovations of the two constrained filters form a ratio for a Wald sequential probability ratio test. The Wald test guides risk mitigation maneuver decisions based on explicit false alarm and missed detection criteria. Since only current-state Kalman filtering is required to compute the innovations for the likelihood ratio, the present approach does not require the mapping of probability density forward to the time of closest approach. Instead, the hard-body constraint manifold is mapped to the filter update time by applying a sigma-point transformation to a projection function. Although many projectors are available, we choose one based on Lambert-style differential correction of the current-state velocity. We have tested our method using a scenario based on the Magnetospheric Multi-Scale mission, scheduled for launch in late 2014. This mission involves formation flight in highly elliptical orbits of four spinning spacecraft equipped with antennas extending 120 meters tip-to-tip. Eccentricities range from 0.82 to 0.91, and close approaches generally occur in the vicinity of perigee, where rapid changes in geometry may occur. Testing the method using two 12,000-case Monte Carlo simulations, we found the method achieved a missed detection rate of 0.1%, and a false alarm rate of 2%.

Carpenter, J. Russell↗
NASA NTRS2012

Representation of Probability Density Functions from Orbit Determination using the Particle Filter

Statistical orbit determination enables us to obtain estimates of the state and the statistical information of its region of uncertainty. In order to obtain an accurate representation of the probability density function (PDF) that incorporates higher order statistical information, we propose the use of nonlinear estimation methods such as the Particle Filter. The Particle Filter (PF) is capable of providing a PDF representation of the state estimates whose accuracy is dependent on the number of particles or samples used. For this method to be applicable to real case scenarios, we need a way of accurately representing the PDF in a compressed manner with little information loss. Hence we propose using the Independent Component Analysis (ICA) as a non-Gaussian dimensional reduction method that is capable of maintaining higher order statistical information obtained using the PF. Methods such as the Principal Component Analysis (PCA) are based on utilizing up to second order statistics, hence will not suffice in maintaining maximum information content. Both the PCA and the ICA are applied to two scenarios that involve a highly eccentric orbit with a lower apriori uncertainty covariance and a less eccentric orbit with a higher a priori uncertainty covariance, to illustrate the capability of the ICA in relation to the PCA.

Mashiku, Alinda K.↗