Search NASA⌕ Search

SEARCH · Search NASA

Results for “truncation uncertainty”

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.

33 records · Page 2

Uncertainty estimates in geomagnetic field modeling

This paper presents an extension of the conventional uncertainty analysis which characterizes the sources of uncertainty in the coefficients of the geomagnetic field models. The new formalism accounts for the systematic errors introduced by the omission of such sources as the presence of crustal fields, the external fields, and the field from the truncated terms. The usefulness of this formalism depends on two critical conditions. The first of these is the knowledge of the statistical properties of the fields whose parameters are not solved in the analysis, i.e., the crustal field and the external field. The second critical point in the practical use of the method is the approximation used for the weight matrix.

Langel, R. A.↗

Space Station dynamic modeling, control and simulation

Dynamic disturbances of many orders of magnitude greater than those of conventional spacecraft will be routine for the Space Station. Accurate knowledge of inflight structural dynamics and disturbances will be lacking. System identification will reduce uncertainties. To deal with the remaining model errors and time-varying elements, adaptive control may be required. Dynamic models for two Space Station configurations are derived. A direct model reference adaptive control algorithm is synthesized and evaluated with respect to the station model parameter errors and plant dynamics truncations. Both the rigid body and the flexible modes are treated. Simulation results show that convergence of the adaptive algorithm can be achieved in 100 to about 300 seconds with reasonable performance even during Shuttle hard docking operations in which station mass and inertia are instantaneously changed by more than 100 percent.

Ih, C.-H. C.↗

Scramjet Engine Flowpath Development for the Hyper-X Mach 10 Flight Test

The third flight of the Hyper-X Research Vehicle at nearly Mach 10 in November 2004 demonstrated the potential for airbreathing propulsion at hypersonic speeds. The flowpath of the scramjet engine used in the Hyper-X Research Vehicle was developed using computational fluid dynamic based design method in combination with a systematic series of ground tests in the NASA HyPulse Shock Tunnel at conditions duplicating Mach 10 flight. The test model was a true scale height and length of the internal flight engine flowpath that was partial width and truncated fore and aft. Tests were conducted using a mixture of hydrogen augmented with two-percent silane. Silane gas served as an ignition aid during the short duration of the pulse facility tests. The engine inflow values of pressure, temperature, and Mach number were parametrically varied during the ground test entries to broaden the database over the expected uncertainty in the flight conditions. A review of the ground test technique and comparisons of the ground test pressures along with selected data from the third flight are presented.

Rogers, R. C.↗

Details of Side Load Test Data and Analysis for a Truncated Ideal Contour Nozzle and a Parabolic Contour Nozzle

Two cold flow subscale nozzles were tested for side load characteristics during simulated nozzle start transients. The two test article contours were a truncated ideal and a parabolic. The current paper is an extension of a 2009 AIAA JPC paper on the test results for the same two nozzle test articles. The side load moments were measured with the strain tube approach in MSFC s Nozzle Test Facility. The processing techniques implemented to convert the strain gage signals into side load moment data are explained. Nozzle wall pressure profiles for separated nozzle flow at many NPRs are presented and discussed in detail. The effect of the test cell diffuser inlet on the parabolic nozzle s wall pressure profiles for separated flow is shown. The maximum measured side load moments for the two contours are compared. The truncated ideal contour s peak side load moment was 45% of that of the parabolic contour. The calculated side load moments, via mean-plus-three-standard-deviations at each nozzle pressure ratio, reproduced the characteristics and absolute values of measured maximums for both contours. The effect of facility vibration on the measured side load moments is quantified and the effect on uncertainty is calculated. The nozzle contour designs are discussed and the impact of a minor fabrication flaw in the nozzle contours is explained.

Ruf, Joseph H.↗

Orthogonal series generalized likelihood ratio test for failure detection and isolation

A new failure detection and isolation algorithm for linear dynamic systems is presented. This algorithm, the Orthogonal Series Generalized Likelihood Ratio (OSGLR) test, is based on the assumption that the failure modes of interest can be represented by truncated series expansions. This assumption leads to a failure detection algorithm with several desirable properties. Computer simulation results are presented for the detection of the failures of actuators and sensors of a C-130 aircraft. The results show that the OSGLR test generally performs as well as the GLR test in terms of time to detect a failure and is more robust to failure mode uncertainty. However, the OSGLR test is also somewhat more sensitive to modeling errors than the GLR test.

Hall, Steven R.↗

Estimating the Concentration of Large Raindrops from Polarimetric Radar and Disdrometer Observations

Estimation of rainfall integral parameters, including radar observables, and empirical relations between them are sensitive to the truncation of the drop size distribution (DSD), particularly at the large drop end. The sensitivity of rainfall integral parameters to the maximum drop diameter (D(sub max)) is exacerbated at C‐band since resonance effects are pronounced for large drops in excess of 5 mm diameter (D). Due to sampling limitations, it is often difficult to reliably estimate D(sub max) with disdrometers. The resulting uncertainties in D(sub max0 potentially increase errors in radar retrieval methods, particularly at C‐band, that rely on disdrometer observations for DSD input to radar models. In fact, D(sub max) is typically an assumed DSD parameter in the development of radar retrieval methods. Because of these very uncertainties, it is difficult to independently confirm disdrometer estimates of D(sub max) with polarimetric radar observations. A couple of approaches can be taken to reduce uncertainty in large drop measurement. Longer integration times can be used for the collection of larger disdrometer samples. However, integration periods must be consistent with a radar resolution volume (RRV) and the temporal and spatial scales of the physical processes affecting the DSD therein. Multiple co‐located disdrometers can be combined into a network to increase the sample size within a RRV. However, over a reasonable integration period, a single disdrometer sample volume is many orders of magnitudes less than a RRV so it is not practical to devise a network of disdrometers that has an equivalent volume to a typical RRV. Since knowledge of DSD heterogeneity and large drop occurrence in time and space is lacking, the specific accuracy or even general representativeness of disdrometer based D(sub max) and large drop concentration estimates within a RRV are currently unknown. To address this complex issue, we begin with a simpler question. Is the frequency of occurrence of large rain drops (D > 5 mm) in disdrometer observations, either stand alone or networked, generally representative and consistent with polarimetric radar observations? We first show from simulations that the concentration of large (D > 5 mm) rain drops (N(sub T5)) can be estimated from polarimetric observations of specific differential phase (K(sub dp)) and differential reflectivity (Z(sub dr)), N(sub T5)=F(K(sub dp),Z(sub dr)), or horizontal reflectivity (Z(sub h)) and Z(sub dr), N(sub T5)=(Z(sub h),Z(sub dr)). We assess the error associated with polarimetric retrieval of N(sub T5), including sensitivity to D(sub max) parameterization assumptions and measurement error in the radar simulations. Polarimetric measurements at S‐band and C‐band will then be used to retrieve estimates of N(sub T5) and compared to disdrometer estimates of N(sub T5). After careful consideration of retrieval error, we will check consistency between disdrometer and polarimetric radar estimates of N(sub T5) and the frequency of occurrence of large rain drops in a variety of precipitating regimes using data from NASA's Global Precipitation Measurement (GPM) Ground Validation (GV) program, including field campaigns such as MC3E (Oklahoma) and IFloodS (Iowa) and extended measurements over Huntsville, Alabama and NASA Wallops Flight Facility in coastal Virginia.

Carey, Lawrence D.↗

Spatial distribution of the X-ray-emitting plasma of U Geminorum in quiescence and outburst

We present our analysis of the Suzaku data of U Geminorum (U Gem) from 2012 both in quiescence and outburst. Unlike SS Cygni (SS Cyg), the hard X-ray flux of U Gem is known to increase at times of optical outburst. A sophisticated spectral model and reliable distance estimate now reveal that this can be attributed to the fact that the mass accretion rate onto the white dwarf (WD) does not exceed the critical rate that causes the optically thin to thick transition of the boundary layer. From comparison of the X-ray and optical light curves, the X-ray outburst peak seems to be retarded by 2.1 ± 0.5 d, although there remains uncertainty in the X-ray peak identification, due to short data coverage. The larger delay than SS Cyg (0.9–1.4 d) also supports the lower accretion rate in U Gem. A fluorescent iron 6.4 keV emission line bears significant information about the geometry of the X-ray-emitting hot plasma and the accretion disk (AD) that reflects the hard X-ray emission. Our reflection simulation has shown that the optically thick AD is truncated at a distance of 1.20–1.25 times the white dwarf radius (RWD) in quiescence, and the accreting matter in the disk turns into the optically thin hard-X-ray-emitting plasma at this radius. In outburst, on the other hand, our spectral analysis favors the picture that the optically thick disk reaches the WD surface, although disk truncation can take place in the region of <1.012 RWD. From the profile of the 6.4 keV line, we have also discovered that the accreting matter is heated up close to the maximum temperature immediately after the matter enters the boundary layer at the disk truncation radius. This is consistent with the fact that the hard X-ray spectra of dwarf novae, in general, can be well represented with the cooling flow model.

Mai Takeo↗

Numerical Experimentation with Maximum Likelihood Identification in Static Distributed Systems

Many important issues in the control of large space structures are intimately related to the fundamental problem of parameter identification. One might also ask how well this identification process can be carried out in the presence of noisy data since no sensor system is perfect. With these considerations in mind the algorithms herein are designed to treat both the case of uncertainties in the modeling and uncertainties in the data. The analytical aspects of maximum likelihood identification are considered in some detail in another paper. The questions relevant to the implementation of these schemes are dealt with, particularly as they apply to models of large space structures. The emphasis is on the influence of the infinite dimensional character of the problem on finite dimensional implementations of the algorithms. Those areas of current and future analysis are highlighted which indicate the interplay between error analysis and possible truncations of the state and parameter spaces.

Scheid, R. E., Jr.↗

Dynamic modeling and adaptive control for space stations

Of all large space structural systems, space stations present a unique challenge and requirement to advanced control technology. Their operations require control system stability over an extremely broad range of parameter changes and high level of disturbances. During shuttle docking the system mass may suddenly increase by more than 100% and during station assembly the mass may vary even more drastically. These coupled with the inherent dynamic model uncertainties associated with large space structural systems require highly sophisticated control systems that can grow as the stations evolve and cope with the uncertainties and time-varying elements to maintain the stability and pointing of the space stations. The aspects of space station operational properties are first examined, including configurations, dynamic models, shuttle docking contact dynamics, solar panel interaction, and load reduction to yield a set of system models and conditions. A model reference adaptive control algorithm along with the inner-loop plant augmentation design for controlling the space stations under severe operational conditions of shuttle docking, excessive model parameter errors, and model truncation are then investigated. The instability problem caused by the zero-frequency rigid body modes and a proposed solution using plant augmentation are addressed. Two sets of sufficient conditions which guarantee the globablly asymptotic stability for the space station systems are obtained.

Ih, C. H. C.↗

Lunar Prospector Orbit Determination Uncertainties Using the High Resolution Lunar Gravity Models

The Lunar Prospector (LP) mission began on January 6, 1998, when the LP spacecraft was launched from Cape Canaveral, Florida. The objectives of the mission were to determine whether water ice exists at the lunar poles, generate a global compositional map of the lunar surface, detect lunar outgassing, and improve knowledge of the lunar magnetic and gravity fields. Orbit determination of LP performed at the Jet Propulsion Laboratory (JPL) is conducted as part of the principal science investigation of the lunar gravity field. This paper will describe the JPL effort in support of the LP Gravity Investigation. This support includes high precision orbit determination, gravity model validation, and data editing. A description of the mission and its trajectory will be provided first, followed by a discussion of the orbit determination estimation procedure and models. Accuracies will be examined in terms of orbit-to-orbit solution differences, as a function of oblateness model truncation, and inclination in the plane-of-sky. Long term predictions for several gravity fields will be compared to the reconstructed orbits to demonstrate the accuracy of the orbit determination and oblateness fields developed by the Principal Gravity Investigator.

Carranza, Eric↗

Spectro-Timing Study of GX 339-4 in a Hard Intermediate State

We present an analysis of Nuclear Spectroscopic Telescope Array observations of a hard intermediate state of the transient black hole GX 339-4 taken in 2015 January. With the source softening significantly over the course of the 1.3 day long observation we split the data into 21 sub-sets and find that the spectrum of all of them can be well described by a power-law continuum with an additional relativistically blurred reflection component. The photon index increases from approx. 1.69 to approx. 1.77 over the course of the observation. The accretion disk is truncated at around nine gravitational radii in all spectra. We also perform timing analysis on the same 21 individual data sets, and find a strong type-C quasi-periodic oscillation (QPO), which increases in frequency from approx. 0.68 to approx. 1.05 Hz with time. The frequency change is well correlated with the softening of the spectrum. We discuss possible scenarios for the production of the QPO and calculate predicted inner radii in the relativistic precession model as well as the global disk mode oscillations model. We find discrepancies with respect to the observed values in both models unless we allow for a black hole mass of approx. 100 Mass compared to the Sun, which is highly unlikely. We discuss possible systematic uncertainties, in particular with the measurement of the inner accretion disk radius in the relativistic reflection model. We conclude that the combination of observed QPO frequencies and inner accretion disk radii, as obtained from spectral fitting, is difficult to reconcile with current models.

accretion disks↗

Three CA-Related Covariance Issues and Their Solutions

This presentation addresses three issues that arise in the use of DoD-produced satellite state estimate covariances in the conjunction assessment process: the realism of the provided covariances, how to address correlated error between two satellites' covariance matrices, and how to proceed when a furnished covariance is non-positive-definite. To address the first, DoD has implemented a set of two "consider parameters" with which to alter/expand the covariance to account for atmospheric density forecast error and satellite frontal area uncertainty, the two largest sources of unmodeled position prediction error for LEO orbits; these values are governed by a satellite's orbital parameters and ballistic coefficient histories, as well as the current and predicted space weather situation. This approach substantially improves the realism of the covariance by accounting for expected prediction errors that are not part of the fit process that generates the covariance. The second issue of covariance correlation is a complicated one, but the main shared error source—global atmospheric density error that is common to both satellites can be characterized and, through sensitivity vectors, quantified at the conjunction's time of closest approach and removed from the joint covariance, which is used in the CA probability of collision (Pc) calculation. Finally, while according to the orbit determination theory a non-positive-definite (NPD) covariance is not possible, numerical truncation and covariance interpolation can conspire to produce NPD results, which render the covariance unusable for certain CA risk assessment calculations, such as Monte Carlo Pc determination. Three similar techniques are profiled and the simplest of them recommended as a reasonable remediation technique when NPD covariances are encountered by CA practitioners.

Dynamic Consider Parameter↗

Photoionization of disk galaxies: An explanation of the sharp edges in the H I distribution

We have reproduced the observed radial truncation of the H I distribution in isolated spiral galaxies with a model in which extragalactic radiation photoionizes the gaseous disk. For a galactic mass distribution model that reproduces the observed rotation curves, including dark matter in the disk and halo, the vertical structure of the gas is determined self-consistently. The ionization structure and column densities of H and He ions are computed by solving the radiation transfer equation for both continuum and lines. Our model is similar to that of Maloney, and the H I structure differs by less than 10%. The radial structure of the column density of H I is found to be more sensitive to the extragalactic radiation field than to the distribution of mass. For this reason, considerable progress can be made in determining the extragalactic flux of ionizing photons, phi(sub ex), with more 21 cm observations of isolated galaxies. However, owing to the uncertainty of the radial distribution of total hydrogen at large radii, inferring the extragalactic flux by comparing the observed edges to photoionization models is somewhat subjective. We find 1 x 10(exp 4)/sq cm/s is less than or approximately phi(sub ex) is less than or approximately 5 x 10(exp 4)/sq cm/s, corresponding to 2.1 is less than or approximately iota(sub 0) is less than or approximately 10.5 x 10(exp -23) ergs/sq cm/s/Hz/sr for a 1/nu spectrum. Although somewhat higher, our inferred range of iota(sub 0) is consistent with the large range of values obtained by Kulkarni & Fall from the 'proximity effect' toward Quasi-Stellar Objects (QSOs) at approximately 0.5.

Dove, James B.↗

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗

Time Changes of the European Gravity Field from GRACE: A Comparison with Ground Measurements from Superconducting Gravimeters and with Hydrology Model Predictions

We investigate the time-variable gravity changes in Europe retrieved from the initial GRACE monthly solutions spanning a 18 month duration from April 2002 to October 2003. Gravity anomaly maps are retrieved in Central Europe from the monthly satellite solutions we compare the fields according to various truncation levels (typically between degree 10 and 20) of the initial fields (expressed in spherical harmonics to degree 120). For these different degrees, an empirical orthogonal function (EOF) decomposition of the time-variable gravity field leads us to its main spatial and temporal characteristics. We show that the dominant signal is found to be annual with an amplitude and a phase both in agreement with predictions in Europe modeled using snow and soil-moisture variations from recent hydrology models. We compare these GRACE gravity field changes to surface gravity observations from 6 superconducting gravimeters of the GGP (Global Geodynamics Project) European sub-network, with a special attention to loading corrections. Initial results suggest that all 3 data sets (GRACE, hydrology and GGP) are responding to annual changes in near-surface water in Europe of a few microGal (at length scales of approx.1000 km) that show a high value in winter and a summer minimum. We also point out that the GRACE gravity field evolution seems to indicate that there is a trend in gravity between summer 2002 and summer 2003 which can be related to the 2003 heatwave in Europe and its hydrological consequences (drought). Despite the limited time span of our analysis and the uncertainties in retrieving a regional solution from the network of gravimeters, the calibration and validation aspects of the GRACE data processing based on the annual hydrology cycle in Europe are in progress.

Hinderer, J.↗