Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error propagation”

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 541 records · Page 30

Filter Tuning Using the Chi-Squared Statistic

The Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF) performs orbit determination (OD) for the Aqua and Aura satellites. Both satellites are located in low Earth orbit (LEO), and are part of what is considered the A-Train satellite constellation. Both spacecraft are currently in the science phase of their respective missions. The FDF has recently been tasked with delivering definitive covariance for each satellite.The main source of orbit determination used for these missions is the Orbit Determination Toolkit developed by Analytical Graphics Inc. (AGI). This software uses an Extended Kalman Filter (EKF) to estimate the states of both spacecraft. The filter incorporates force modelling, ground station and space network measurements to determine spacecraft states. It also generates a covariance at each measurement. This covariance can be useful for evaluating the overall performance of the tracking data measurements and the filter itself. An accurate covariance is also useful for covariance propagation which is utilized in collision avoidance operations. It is also valuable when attempting to determine if the current orbital solution will meet mission requirements in the future.This paper examines the use of the Chi-square statistic as a means of evaluating filter performance. The Chi-square statistic is calculated to determine the realism of a covariance based on the prediction accuracy and the covariance values at a given point in time. Once calculated, it is the distribution of this statistic that provides insight on the accuracy of the covariance.For the EKF to correctly calculate the covariance, error models associated with tracking data measurements must be accurately tuned. Over estimating or under estimating these error values can have detrimental effects on the overall filter performance. The filter incorporates ground station measurements, which can be tuned based on the accuracy of the individual ground stations. It also includes measurements from the NASA space network (SN), which can be affected by the assumed accuracy of the TDRS satellite state at the time of the measurement.The force modelling in the EKF is also an important factor that affects the propagation accuracy and covariance sizing. The dominant force in the LEO orbit regime is the drag force caused by atmospheric drag. Accurate accounting of the drag force is especially important for the accuracy of the propagated state. The implementation of a box and wing model to improve drag estimation accuracy, and its overall effect on the covariance state is explored.The process of tuning the EKF for Aqua and Aura support is described, including examination of the measurement errors of available observation types (Doppler and range), and methods of dealing with potentially volatile atmospheric drag modeling. Predictive accuracy and the distribution of the Chi-square statistic, calculated based of the ODTK EKF solutions, are assessed versus accepted norms for the orbit regime.

Covariance Analysis↗

Power series evaluation of transition and covariance matrices.

Reexamination power series solutions to the matrix covariance differential equation and the transition differential equation. Truncation error bounds are derived which are computationally attractive and which extend previous results. Polynomial approximations are obtained by exploiting the functional equations satisfied by the transition and covariance matrices. The series-functional equation propagation technique represents a fast and accurate alternative to the numerical integration of the time-invariant transition and covariance equations.

Bierman, G. J.↗

Radio propagation through solar and other extraterrestrial ionized media

The present S- and X-band communications needs in deep space are addressed to illustrate the aspects which are affected by propagation through extraterrestrial plasmas. The magnitude, critical threshold, and frequency dependence of some eight propagation effects for an S-band propagation path passing within 4 solar radii of the Sun are described. The theory and observation of propagation in extraterrestrial plasmas are discussed and the various plasma states along a near solar propagation path are illustrated. Classical magnetoionic theory (cold anisotropic plasma) is examined for its applicability to the path in question. The characteristics of the plasma states found along the path are summarized and the errors in some of the standard approximations are indicated. Models of extraterrestrial plasmas are included. Modeling the electron density in the solar corona and solar wind, is emphasized but some cursory information on the terrestrial planets plus Jupiters is included.

Smith, E. K.↗

Solution of the linear shallow water equations by the fourth-order leapfrog scheme

Numerical schemes of the first and second order of approximation introduce numerical distortion when the wave propagation over a long distance is investigated. To alleviate this problem, the fourth-order leapfrog scheme is constructed. The standard leapfrog method is based on the truncated Taylor series expansion which depicts an error proportional to the second-order terms. In the proposed method the numerical solution is corrected for these terms. The space and time corrections work well in diminishing numerical dispersion and dissipation.

Kowalik, Z.↗

Star Tracker Based ATP System Conceptual Design and Pointing Accuracy Estimation

A star tracker based beaconless (a.k.a. non-cooperative beacon) acquisition, tracking and pointing concept for precisely pointing an optical communication beam is presented as an innovative approach to extend the range of high bandwidth (> 100 Mbps) deep space optical communication links throughout the solar system and to remove the need for a ground based high power laser as a beacon source. The basic approach for executing the ATP functions involves the use of stars as the reference sources from which the attitude knowledge is obtained and combined with high bandwidth gyroscopes for propagating the pointing knowledge to the beam pointing mechanism. Details of the conceptual design are presented including selection of an orthogonal telescope configuration and the introduction of an optical metering scheme to reduce misalignment error. Also, estimates are presented that demonstrate that aiming of the communications beam to the Earth based receive terminal can be achieved with a total system pointing accuracy of better than 850 nanoradians (3 sigma) from anywhere in the solar system.

star tracker↗

The Next Generation of Lightning Mapping

With the removal of the Lightning Imaging Sensor from the International Space Station, a gap has opened in lightning observation from low-Earth orbit. The CubeSpark mission concept aims to fill this role using a constellation of satellites with radio frequency (RF) sensors and bi-spectral optical imagers to observe lightning flashes more completely and with better resolution than is currently possible from space. In this study, we assess the feasibility of multiple methods of not only mapping lightning locations, but also inferring 3D flash and charge structures. This is done primarily by simulating lightning emissions in the very high frequency (VHF) band, modeling their propagation to orbital sensors, and reconstructing their locations using time-of-arrival (TOA) minimization algorithms. Constellation shape, number, and atmospheric conditions are varied in order to quantify the three main sources of geolocation error: geometric, ionospheric, and instrumental effects. The promising results presented here demonstrate 3D resolution better than 1-2 km in many cases, enabling new applications in meteorology and climate sciences.

CubeSpark↗

An improved empirical model for diversity gain on Earth-space propagation paths

An empirical model was generated to estimate diversity gain on Earth-space propagation paths as a function of Earth terminal separation distance, link frequency, elevation angle, and angle between the baseline and the path azimuth. The resulting model reproduces the entire experimental data set with an RMS error of 0.73 dB.

Hodge, D. B.↗

Forecasts of the 500 mb height using a dynamically oriented statistical model

The forecast skill of a simple dynamically inspired statistical model of the Northern Hemisphere 500 mb height field is evaluated in spectral and physical space for a variety of forecast lead times (1-32 days) and predictand averaging times (1-32 days). The model includes viscous damping, wave propagation, climatology and implicit stochastic forcing. The largest model skill was found for forecasts of the zonal flow and the largest waves. In general, the largest forecast skills were also associated with the largest forecast error, there being a slight geographic phase shift of the skill with respect to the error. Model skills for climate (time-averaged) forecasts are greater when using instantaneous rather than time averages to forecast time averages. Analysis of model errors suggests areas for improvement in representing forcing terms and model physics. However, the model error fields are largely 'white noise' which suggests that global forecast skills substantially larger than those obtained here are unlikely to be achieved by more sophisticated models.

Roads, J. O.↗

African Easterly Waves in 30-day High-Resolution Global Simulations: A Case Study During the 2006 NAMMA Period

In this study, extended -range (30 -day) high-resolution simulations with the NASA global mesoscale model are conducted to simulate the initiation and propagation of six consecutive African easterly waves (AEWs) from late August to September 2006 and their association with hurricane formation. It is shown that the statistical characteristics of individual AEWs are realistically simulated with larger errors in the 5th and 6th AEWs. Remarkable simulations of a mean African easterly jet (AEJ) are also obtained. Nine additional 30 -day experiments suggest that although land surface processes might contribute to the predictability of the AEJ and AEWs, the initiation and detailed evolution of AEWs still depend on the accurate representation of dynamic and land surface initial conditions and their time -varying nonlinear interactions. Of interest is the potential to extend the lead time for predicting hurricane formation (e.g., a lead time of up to 22 days) as the 4th AEW is realistically simulated.

Shen, Bo-Wen↗

Statistical Analysis of Instantaneous Frequency Scaling Factor as Derived From Optical Disdrometer Measurements At KQ Bands

The rain rate data and statistics of a location are often used in conjunction with models to predict rain attenuation. However, the true attenuation is a function not only of rain rate, but also of the drop size distribution (DSD). Generally, models utilize an average drop size distribution (Laws and Parsons or Marshall and Palmer. However, individual rain events may deviate from these models significantly if their DSD is not well approximated by the average. Therefore, characterizing the relationship between the DSD and attenuation is valuable in improving modeled predictions of rain attenuation statistics. The DSD may also be used to derive the instantaneous frequency scaling factor and thus validate frequency scaling models. Since June of 2014, NASA Glenn Research Center (GRC) and the Politecnico di Milano (POLIMI) have jointly conducted a propagation study in Milan, Italy utilizing the 20 and 40 GHz beacon signals of the Alphasat TDP#5 Aldo Paraboni payload. The Ka- and Q-band beacon receivers provide a direct measurement of the signal attenuation while concurrent weather instrumentation provides measurements of the atmospheric conditions at the receiver. Among these instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which yields droplet size distributions (DSD); this DSD information can be used to derive a scaling factor that scales the measured 20 GHz data to expected 40 GHz attenuation. Given the capability to both predict and directly observe 40 GHz attenuation, this site is uniquely situated to assess and characterize such predictions. Previous work using this data has examined the relationship between the measured drop-size distribution and the measured attenuation of the link]. The focus of this paper now turns to a deeper analysis of the scaling factor, including the prediction error as a function of attenuation level, correlation between the scaling factor and the rain rate, and the temporal variability of the drop size distribution both within a given rain event and across different varieties of rain events. Index Terms-drop size distribution, frequency scaling, propagation losses, radiowave propagation.

Statistics↗

Statistical Analysis of Instantaneous Frequency Scaling Factor as Derived From Optical Disdrometer Measurements At KQ Bands

The rain rate data and statistics of a location are often used in conjunction with models to predict rain attenuation. However, the true attenuation is a function not only of rain rate, but also of the drop size distribution (DSD). Generally, models utilize an average drop size distribution (Laws and Parsons or Marshall and Palmer [1]). However, individual rain events may deviate from these models significantly if their DSD is not well approximated by the average. Therefore, characterizing the relationship between the DSD and attenuation is valuable in improving modeled predictions of rain attenuation statistics. The DSD may also be used to derive the instantaneous frequency scaling factor and thus validate frequency scaling models. Since June of 2014, NASA Glenn Research Center (GRC) and the Politecnico di Milano (POLIMI) have jointly conducted a propagation study in Milan, Italy utilizing the 20 and 40 GHz beacon signals of the Alphasat TDP#5 Aldo Paraboni payload. The Ka- and Q-band beacon receivers provide a direct measurement of the signal attenuation while concurrent weather instrumentation provides measurements of the atmospheric conditions at the receiver. Among these instruments is a Thies Clima Laser Precipitation Monitor (optical disdrometer) which yields droplet size distributions (DSD); this DSD information can be used to derive a scaling factor that scales the measured 20 GHz data to expected 40 GHz attenuation. Given the capability to both predict and directly observe 40 GHz attenuation, this site is uniquely situated to assess and characterize such predictions. Previous work using this data has examined the relationship between the measured drop-size distribution and the measured attenuation of the link [2]. The focus of this paper now turns to a deeper analysis of the scaling factor, including the prediction error as a function of attenuation level, correlation between the scaling factor and the rain rate, and the temporal variability of the drop size distribution both within a given rain event and across different varieties of rain events. Index Terms-drop size distribution, frequency scaling, propagation losses, radiowave propagation.

Radio Attenuation↗

Investigation of a Parabolic Iterative Solver for Three-dimensional Configurations

A parabolic iterative solution procedure is investigated that seeks to extend the parabolic approximation used within the internal propagation module of the duct noise propagation and radiation code CDUCT-LaRC. The governing convected Helmholtz equation is split into a set of coupled equations governing propagation in the positive and negative directions. The proposed method utilizes an iterative procedure to solve the coupled equations in an attempt to account for possible reflections from internal bifurcations, impedance discontinuities, and duct terminations. A geometry consistent with the NASA Langley Curved Duct Test Rig is considered and the effects of acoustic treatment and non-anechoic termination are included. Two numerical implementations are studied and preliminary results indicate that improved accuracy in predicted amplitude and phase can be obtained for modes at a cut-off ratio of 1.7. Further predictions for modes at a cut-off ratio of 1.1 show improvement in predicted phase at the expense of increased amplitude error. Possible methods of improvement are suggested based on analytic and numerical analysis. It is hoped that coupling the parabolic iterative approach with less efficient, high fidelity finite element approaches will ultimately provide the capability to perform efficient, higher fidelity acoustic calculations within complex 3-D geometries for impedance eduction and noise propagation and radiation predictions.

Nark, Douglas M.↗

Improving Computational Efficiency of Prognostics Algorithms in Resource-Constrained Settings

In engineering and aerospace applications, it is vital to operational success to have insight into the expected performance and health of physical systems. The field of prognostics and health management provides quantitative methods for monitoring, predicting, and managing system health. Prognostics algorithms can be employed to assess the current state of a system, propagate the system throughout time, and predict potential anomalies or failures that may occur. While they can provide accurate prediction results, effective prognostics algorithms can be challenging to use in resource-constrained settings due to computational limitations and high computational latency, leading to obsolete predictions. Thus, computationally efficient and accurate algorithms are necessary for future remaining useful life predictions. In this work, we implement new algorithmic approaches for prediction, quantitatively compare them via a battery degradation use-case, and provide recommendations of potential improvements to a prognostics framework. One approach to prediction is through sampling, whereby the current state of a physical system is sampled many times and each sample is propagated forward until failure is reached, resulting in a distribution of failure values. To improve the efficiency of this process, we implemented five new algorithmic approaches to prediction, including three distinct sampling methods (standard Monte Carlo, Quasi-Monte Carlo, and Latin Hypercube Sampling), a variable time step algorithm, and a variable sample size algorithm. To compare the algorithms, we employ a variety of metrics designed specifically to analyze both computational efficiency and model accuracy. Our metrics include accuracy to compare the average predicted value to ground truth, mean absolute deviation to illustrate dispersion, specific percentile error to describe accuracy within a user-defined risk tolerance, and code run-time. To quantitatively analyze our results, we employ a use-case of degradation of a Lithium-ion battery. We use an electrochemistry-based model to describe the current health state of the battery, and implement our prediction algorithms to propagate forward in time until end-of-discharge (EOD) is reached. Notably, through this work it was found that none of our sampling approaches had a significant impact on computational efficiency or model accuracy in predicting EOD of the battery. We find that while the sampling methods are unique, the distributions they generate are similar, ultimately producing final predictions that are nearly identical. In exploring the effect of the time step within the prediction algorithm, we found that prediction accuracy was highly dependent on the time step used, and that implementing a variable time step within a particular prediction may provide an increase in computational efficiency while also maintaining prediction accuracy. Finally, implementing a variable sample size also affected prediction, and our results show that tuning both the magnitude and timing of the sample size adjustment can result in improved computation speed and maintained prediction accuracy. Taken together, our findings highlight the challenge of performing prognostics in resource-constrained settings, and illustrate the potential of developing new prediction algorithms to improve computational efficiency.

prognostics↗

Theoretical Modeling for Numerical Weather Prediction

The goal is to utilize predictability theory and numerical experimentation to identify and understand some of the dynamical processes which must be modeled more realistically if large-scale numerical weather predictions are to be improved. The emphasis is on the use of relatively simple models to exlore the properties of physically comprehensive general circulation models (GCM's). A global linear quasi-geostrophic model and the Goddard Laboratory for Atmospheric Sciences (GLAS) GCM were used to investigate several mechanisms which are responsible for the decay of large-scale forecast skill in mid-latitude numerical weather predictions. Five-day forecasts for an ensemble of cases were made using First GARP Global Experiment data. It was found that forecast skill depends crucially on the specification of the stationary forcing. A lack of stationary forcing leads to spurious westwad propagation of the ultralong waves. Forecasts made with stationary forcings derived from climatological data are superior to those using forcings inferred from observations immediately preceding the forecast period. Interhemispheric forecast differences were analyzed, and the model errors were compared to errors of a simple persistence-damped-to-climatology scheme and to errors of the GLAS GCM.

Somerville, R. C. J.↗

Improved Coplanar Waveguides

Three new types of coplanar waveguide transmission lines proposed: suspended, striplinelike suspended, and inverted. Numerical analyses of propagation of electromagnetic fields of new and of more-conventional plain and conductor-backed coplanar waveguides showed new waveguides offered several advantages over conventional ones. Configurations promise lower attenuation and reduce effects of manufacturing errors.

Simons, Rainee N.↗

Comparison of several finite-difference methods

The present numerical study establishes that, with the MacCormack and 'two-four' methods, the solution accuracy obtainable in problems involving wave propagation, shock-wave and contact discontinuities, and viscous effects, will be strongly dependent on Courant number. The application of flux correction to the MacCormack and two-four methods is noted to significantly attenuate dispersion errors; the ensuing solutions capture the discontinuities in the shock-tube problem with improved accuracy and resolution, and are free of dispersion errors for the viscous Burgers' equation.

Biringen, S.↗

Horizontal fields generated by return strokes

Horizontal fields generated by return strokes play an important role in the interaction of lightning generated electric fields with power lines. In many of the recent investigations on the interaction of lightning electromagnetic fields with power lines, the horizontal field was calculated by employing the expression for the tilt of the electric field of a plane wave propagating over finitely conducting earth. The method is suitable for calculating horizontal fields generated by return strokes at distances as close as 200m. At these close ranges, the use of the wavetilt expression can cause large errors.

Cooray, Vernon↗

Estimate Low and High Order Wavefront Using P1640 Calibrator Measurements

P1640 high contrast imaging system on the Palomar 200 inch Telescope consists of an apodized-pupil Lyot coronagraph, the PALM-3000 adaptive optics (P3K-AO), and P1640 Calibrator (CAL). Science images are recorded by an integral field spectrograph covering J-H bands for detecting and characterizing stellar companions. With aberrations from atmosphere corrected by the P3K-AO, instrument performance is limited mainly by the quasi-static speckles due to noncommon path wavefront aberrations for the light to propagate to the P3K-AO wavefront sensor and to the coronagraph mask. The non-common path wavefront aberrations are sensed by CAL, which measures the post-coronagraph E-field using interferometry, and can be effectively corrected by offsetting the P3K-AO deformable mirror target position accordingly. Previously, we have demonstrated using CAL measurements to correct high order wavefront aberrations, which is directly connected to the static speckles in the image plane. Low order wavefront, on the other hand, usually of larger amplitudes, causes light to leak through the coronagraph making the whole image plane brighter. Knowledge error in low order wavefront aberrations can also affect the estimation of the high order wavefront. Even though, CAL is designed to sense efficiently high order wavefront aberrations, the low order wavefront front can be inferred with less sensitivity. Here, we describe our method for estimating both low and high order wavefront aberrations using CAL measurements by propagating the post-coronagraph E-field to a pupil before the coronagraph. We present the results from applying this method to both simulated and experiment data.

aberrations↗