SEARCH · Engineering Papers
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.
Our understanding of snow distributions in the polar regions is severly restricted due to the heterogeneity, both in space and time, of this solid precipitate. Processes such as vapor and mass fluxes across the interface are, to a large extent, controlled by the presence and geophysical state of the snow cover on sea ice.
Beginning in 1997, the Merged Precipitation Group at NASA Goddard has distributed gridded global precipitation products built by combining satellite and surface gauge data. This started with the Global Precipitation Climatology Project (GPCP), then the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA), and recently the Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (GPM) mission (IMERG). This 20+-year (and on-going) activity has yielded an important set of insights and lessons learned for making state-of-the-art precipitation data accessible to the diverse communities of users. Merged-data products critically depend on the input sensors and the retrieval algorithms providing accurate, reliable estimates, but it is also important to provide ancillary information that helps users determine suitability for their application. We typically provide fields of estimated random error, and recently reintroduced the quality index concept at user request. Also at user request we have added a (diagnostic) field of estimated precipitation phase. Over time, increasingly more ancillary fields have been introduced for intermediate products that give expert users insight into the detailed performance of the combination algorithm, such as individual merged microwave and microwave-calibrated infrared estimates, the contributing microwave sensor types, and the relative influence of the infrared estimate.
Atmospheric carbon dioxide (CO2) inversions for estimating natural carbon fluxes typically do not allow for adjustment of fossil fuel CO2 emissions, despite significant uncertainties in emission inventories and inadequacies in the specification of international bunker emissions in inversions. Also, most inversions place CO2 release from fossil fuel combustion and biospheric sources entirely at the surface. However, a non-negligible portion of the emissions actually occurs in the form of reduced carbon species, which are eventually oxidized to CO2 downwind. Omission of this ‘chemical pump’ can result in a significant redistribution of the inferred total carbon fluxes among regions. We assess the impacts of different prescriptions of fossil fuel emissions and accounting for the chemical pump on flux estimation, with a novel aspect of conducting both satellite CO2 observation-based and surface in situ-based inversions. We apply 3-D carbon monoxide (CO) loss rates archived from a state-of-the-art GEOS chemistry and climate model simulation in a forward transport model run to simulate the distribution of CO2 originating from oxidation of carbon species. We also subtract amounts from the prior surface CO2 fluxes that are actually emitted in the form of fossil and biospheric CO, methane, and nonmethane volatile organic compounds (VOCs). We find that the posterior large-scale fluxes are generally insensitive to the finer-scale spatial differences between the ODIAC and CDIAC fossil fuel CO2 gridded datasets and assumptions about international bunker emissions. However, accounting for 3-D chemical CO2 production and the surface correction shifts the global carbon sink, e.g., from land to ocean and from the tropics to the north, with a magnitude and even direction that depend on assumptions about the surface correction. A GOSAT satellite-based inversion is more sensitive to the chemical pump than one using in situ observations, exhibiting substantial flux impacts of 0.28, 0.53, and -0.47 Pg C yr(exp -1) over tropical land, global land, and oceans, due to differences in the horizontal and vertical sampling of the wo observation types. Overall, the biases from neglecting the chemical pump appear to be minor relative to the flux estimate uncertainties and the differences between the in situ and GOSAT inversions, but their relative importance will grow in the future as observational coverage further increases and satellite retrieval biases decrease.
The difference equations of Kalman filtering and smoothing recursively factor and invert the covariance of the output of a linear state-space system driven by a white-noise process. Here it is shown that similar recursive techniques factor and invert the inertia matrix of a multibody robot system. The random field models are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. They are easier to describe than the models based on classical mechanics, which typically require extensive derivation and manipulation of equations of motion for complex mechanical systems. With the spatially random models, more primitive locally specified computations result in a global collective system behavior equivalent to that obtained with deterministic models. The primary goal of applying random field estimation is to provide a concise analytical foundation for solving robot control and motion planning problems.
S-band time delay measurements were collected from the spacecraft Helios A and B during three solar occultations in 1975/76 within heliocentric distances of about 3 and 215 earth radius in terms of range, Doppler frequency shift, and electron content. Characteristic features of measurement and data processing are described. Typical data sets are discussed to probe the electron density distribution near the sun (west and east limb as well) including the outer and extended corona. Steady-state and dynamical aspects of the solar corona are presented and compared with earth-bound-K-coronagraph measurements. Using a weighted least squares estimation, parameters of an average coronal electron density profile are derived in a preliminary analysis to yield electron densities at r = 3, 65, 215 earth radius. Transient phenomena are discussed and a velocity of propagation v is nearly equal to 900 km/s is determined for plasma ejecta from a solar flare observed during an extraordinary set of Helios B electron content measurements.
Reference is made to a study by Casasent et al. (1983), which gave a description of a frequency-multiplexed acoustooptic processor and showed how it was capable of performing all the individual operations required in Kalman filtering. The data flow and organization of all required operations however, were not detailed in that study. Consideration is given here to a simpler Kalman filter state estimation problem. Equally spaced time-sampled intervals (k times T sub s, with k the iterative time index) are assumed. It is further assumed that the system noise vector w and the measurement noise vector v are uncorrelated and Gaussian distributed and that the noise statistics (Q and R) and the system model (Phi, Gamma, H) are known. The error covariance matrix P and the extrapolated error covariance matrix M can thus be precomputed and the Kalman gain matrix K sub k can be precomputed and stored for each input time sample.
The VLBI data set spanning Aug. 1919 to the present is a major constituent of ITRF2000 because of the time span, the geometric strength and the accuracy of the modeling. Formal velocity errors for the best stations are less than 0.3 mm/yr while the probable frame definition error for station velocities is approximately 0.4 mm/yr. The primary systematic errors are the troposphere and tidal loading (solid Earth and ocean). However, the scale of the VLBI TRF is probably only weakly affected because of continued improvements in the tropospheric mapping function and gradient estimation. There are some weaknesses in the spatial distribution of data, particularly in the southern hemisphere, and the quality of the mobile VLBI results is now suboptimal. The deployment of the TIGO system to Chile will improve the geometric coverage. The current state of the art of VLBI analysis will be discussed.
This paper presents a decentralized observer with a consensus filter for the state observation of a discrete-time linear distributed systems. In this setup, each agent in the distributed system has an observer with a model of the plant that utilizes the set of locally available measurements, which may not make the full plant state detectable. This lack of detectability is overcome by utilizing a consensus filter that blends the state estimate of each agent with its neighbors' estimates. We assume that the communication graph is connected for all times as well as the sensing graph. It is proven that the state estimates of the proposed observer asymptotically converge to the actual plant states under arbitrarily changing, but connected, communication and sensing topologies. As a byproduct of this research, we also obtained a result on the location of eigenvalues, the spectrum, of the Laplacian for a family of graphs with self-loops.
Model-based diagnostics and prognostics rely on state estimation and uncertainty management algorithms to produce useful information for system operators and maintainers. This information enables more informed operational decisions, condition-based maintenance, and overall mission safety assurance. Typically, uncertainty is associated with vehicle state-of-health estimation and prediction results because of modeling errors, internal or external sources of noise, and sensor inaccuracy. Probabilistic uncertainty management methods including Sequential Monte Carlo simulation are commonly used to reason about state-of-health estimates and predictions in the presence of these sources of uncertainty. However, such algorithms can be computationally expensive as they require a very large number of samples to obtain a sufficiently accurate quantification of the end of life probability distribution. As a result, highly mobile autonomous systems that leverage the prognostic results for mission-level replanning are often constrained in their processing capability because of these computationally expensive simulation approaches. Therefore, in this paper, we investigate algorithmic methods for dynamically adjusting simulation time step as well as number of samples to achieve highly efficient prognostic results while maintaining results accuracy. Results obtained from simulated flight experiments of an electric unmanned aerial vehicle are presented to verify the efficacy of such algorithms.
This work explores on-board planning for the single spacecraft, multiple ground station Earth-observing satellite scheduling problem through artificial neural network function approximation of state–action value estimates generated by Monte Carlo tree search (MCTS). An extensive hyperparameter search is conducted for MCTS on the basis of performance, safety, and downlink opportunity utilization to determine the best hyperparameter combination for data generation. A hyperparameter search is also conducted on neural network architectures. The learned behavior of each network is explored, and each network architecture’s robustness to orbits and epochs outside of the training distributions is investigated. Furthermore, each algorithm is compared with a genetic algorithm, which serves to provide a baseline for optimality. MCTS is shown to compute near-optimal solutions in comparison to the genetic algorithm. The state–action value networks are shown to match or exceed the performance of MCTS in six orders of magnitude less execution time, showing promise for execution on board spacecraft.
Many a time an objective of a statistical analysis is to estimate a limit value like 3-sigma 95% confidence upper limit from a data sample. The generalized Extreme Value Distribution method can be profitably employed in many situations for such an estimate. . .. It is well known that according to the Central Limit theorem the mean value of a large data set is normally distributed irrespective of the distribution of the data from which the mean value is derived. In a somewhat similar fashion it is observed that many times the extreme value of a data set has a distribution that can be formulated with a Generalized Distribution. In space shuttle entry with 3-string GPS navigation the Figure Of Merit (FOM) value gives a measure of GPS navigated state accuracy. A GPS navigated state with FOM of 6 or higher is deemed unacceptable and is said to form a FOM 6 or higher chimney. A FOM chimney is a period of time during which the FOM value stays higher than 5. A longer period of FOM of value 6 or higher causes navigated state to accumulate more error for a lack of state update. For an acceptable landing it is imperative that the state error remains low and hence at low altitude during entry GPS data of FOM greater than 5 must not last more than 138 seconds. I To test the GPS performAnce many entry test cases were simulated at the Avionics Development Laboratory. Only high value FoM chimneys are consequential. The extreme value statistical technique is applied to analyze high value FOM chimneys. The Maximum likelihood method is used to determine parameters that characterize the GEV distribution, and then the limit value statistics are estimated.
Transonic buffet forces are a major source of unsteady loading on launch vehicles, thus requiring accurate estimation for efficient vehicle design. The state of the art in modeling these unsteady loads utilizes wind-tunnel tests where the fluctuating pressures are measured by pressure transducers (PTs) at discrete locations on a rigid buffet model. These pressures are then integrated over the surface of the vehicle to yield a series of orthogonal centerline loads called buffet forcing functions (BFFs). Typically, the PT layout aims at resolving the pressure correlation along the longitudinal axis of the vehicle. As a result, the distribution of azimuthal correlation and its impact on the estimated BFFs are not well known. To fill these gaps, extremely high-spatial-resolution uPSP data were collected for three different configurations of the Space Launch System in the NASA Ames Research Center 11-Foot Transonic Unitary Plan Wind Tunnel. The spatio-temporal behavior of the pressure correlation on these vehicles is analyzed and flow features of interest are investigated. It is shown that terminal shocks interacting with turbulence are a source of increased azimuthal coherence, especially when the shock develops at a junction. Vortex shedding off the forward attachment hardware that connects the core stage to the solid rocket boosters (SRBs) is the most severe buffet environment on the vehicle. The associated fluctuating pressures are shown to be highly coherent as far as the vehicle tail and up to 40 degrees away from the boosters. For selected areas of the vehicle, factoring the azimuthal coherence into the attenuation of discrete-measurements-based BFFs results in under prediction relative to the BFFs obtained from full integration of the uPSP data.
S-band time delay measurements were collected from the spacecrafts Helios A and B during three solar occultations in 1975/76 within heliocentric distances of about 3 and 215 solar radii in terms of range, Doppler frequency shift, and electron content. A description is given concerning some characteristic features of the methods of measurement and data processing. Typical data sets are discussed to probe the electron density distribution near the sun (west and east limb as well) including the outer and extended corona. Steady-state and dynamical aspects of the solar corona are presented and compared with earth-bound K-coronagraph measurements. Using a weighted least squares estimation 3 parameters of an average coronal electron density profile are derived in a preliminary analysis to yield electron densities of about 130 billion; 100 million; 7 million/cu m at r?3; 65; 215 solar radii. Transient phenomena are discussed and a velocity of propagation v approximately 900 km/s for plasma ejecta from a solar flare is determined from an extraordinary set of Helios B electron content measurements on April 30/May 1, 1976.
The economic benefits of NASA programs are discussed. Emphasis is given to an analysis of indirect economic benefits which estimates the effect of NASA programs on employment, personal income, corporate sales and profits, and government tax revenues in the U.S. and in each state. Data are presented that show that NASA programs have widely varying multipliers by industry and that illustrate the distribution of jobs by industry as well as the distribution of sales.
An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.
An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).
Theoretical calculations are presented which estimate the possible magnitude of the O3/H2O derived OH interference signal resulting from the use of the laser-induced fluorescence technique in measuring natural levels of tropospheric OH. Critical to this new assessment has been the measurement of the nascent OH quantum state distribution resulting from the reaction O(1D) + H2O yields 2OH, and an assessment of the subsequent rotational relaxation of the OH species when formed in high k levels.
The ceramics analysis and reliability evaluation of structures (CARES) computer program is described. The primary function of the code is to calculate the fast-fracture reliability or failure probability of macro-scopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. CARES uses results from MSC/NASTRAN or ANSYS finite-element analysis programs to evaluate how inherent surface and/or volume type flaws component reliability. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effects of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or uniform uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for a single or multiple failure modes by using a least-squares analysis or a maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-to-fit-tests, 90 percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan 90 percent confidence band values are also provided. Examples are provided to illustrate the various features of CARES.