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.
From Intro: A study was initiated in the Santa Monica Mountains to investigate the use of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS)for providing improved maps of chaparral coupled with direct estimates of canopy attributes (eg. biomass, leaf area, fuel load)...Analysis focused on AVIRIS data collected on October 19, 1994.
This work demonstrates the development and implementation of a Fully Constrained Least Squares (FCLS) unmixing model developed in C++ programming language with OpenCV package and boost C++ libraries in the NASA Earth Exchange (NEX). Visualization of the results is supported by GRASS GIS and statistical analysis is carried in R in a Linux system environment. FCLS was first tested on computer simulated data with Gaussian noise of various signal-to-noise ratio, and Landsat data of an agricultural scenario and an urban environment using a set of global end members of substrate (soils, sediments, rocks, and non-photosynthetic vegetation), vegetation that includes green photosynthetic plants and dark objects which encompasses absorptive substrate materials, clear water, deep shadows, etc. For the agricultural scenario, a spectrally diverse collection of 11 scenes of Level 1 terrain corrected, cloud free Landsat-5 TM data of Fresno, California, USA were unmixed and the results were validated with the corresponding ground data. To study an urbanized landscape, a clear sky Landsat-5 TM data were unmixed and validated with coincident World View-2 abundance maps (of 2 m spatial resolution) for an area of San Francisco, California, USA. The results were evaluated using descriptive statistics, correlation coefficient, RMSE, probability of success, boxplot and bivariate distribution function. Finally, FCLS was used for sub-pixel land cover analysis of the monthly WELD (Wen-enabled Landsat data) repository from 2008 to 2011 of North America. The abundance maps in conjunction with DMSP-OLS nighttime lights data were used to extract the urban land cover features and analyze their spatial-temporal growth.
As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.
For the analysis of remotely sensed data, it is frequently necessary to design a classifier in order to locate a ground cover class of interest or to estimate the proportion of this ground cover class. Advantages of a mixture distribution formulation are discussed, and a description is presented of the results of estimating the proportion of small grains in ten Landsat data segments using the mixture model. It is found that the mixture model proportion estimates have a very low variance and coefficient of variation. The discussed investigation implies that the mixtures model is a viable method for determining the distributions of classes of interest in remote sensing problems and in estimating the proportions of these classes directly.
Radiation risks are estimated in a competing risk formalism where age or time after exposure estimates of increased risks for cancer and circulatory diseases are folded with a probability to survive to a given age. The survival function, also called the life-table, changes with calendar year, gender, smoking status and other demographic variables. An outstanding problem in risk estimation is the method of risk transfer between exposed populations and a second population where risks are to be estimated. Approaches used to transfer risks are based on: 1) Multiplicative risk transfer models -proportional to background disease rates. 2) Additive risk transfer model -risks independent of background rates. In addition, a Mixture model is often considered where the multiplicative and additive transfer assumptions are given weighted contributions. We studied the influence of the survival probability on the risk of exposure induced cancer and circulatory disease morbidity and mortality in the Multiplicative transfer model and the Mixture model. Risks for never-smokers (NS) compared to the average U.S. population are estimated to be reduced between 30% and 60% dependent on model assumptions. Lung cancer is the major contributor to the reduction for NS, with additional contributions from circulatory diseases and cancers of the stomach, liver, bladder, oral cavity, esophagus, colon, a portion of the solid cancer remainder, and leukemia. Greater improvements in risk estimates for NS s are possible, and would be dependent on improved understanding of risk transfer models, and elucidating the role of space radiation on the various stages of disease formation (e.g. initiation, promotion, and progression).
Major uncertainties exist regarding the rate and intensity of logging in tropical forests worldwide: these uncertainties severely limit economic, ecological, and biogeochemical analyses of these regions. Recent sawmill surveys in the Amazon region of Brazil show that the area logged is nearly equal to total area deforested annually, but conversion of survey data to forest area, forest structural damage, and biomass estimates requires multiple assumptions about logging practices. Remote sensing could provide an independent means to monitor logging activity and to estimate the biophysical consequences of this land use. Previous studies have demonstrated that the detection of logging in Amazon forests is difficult and no studies have developed either the quantitative physical basis or remote sensing approaches needed to estimate the effects of various logging regimes on forest structure. A major reason for these limitations has been a lack of sufficient, well-calibrated optical satellite data, which in turn, has impeded the development and use of physically-based, quantitative approaches for detection and structural characterization of forest logging regimes. We propose to use data from the EO-1 Hyperion imaging spectrometer to greatly increase our ability to estimate the presence and structural attributes of selective logging in the Amazon Basin. Our approach is based on four "biogeophysical indicators" not yet derived simultaneously from any satellite sensor: 1) green canopy leaf area index; 2) degree of shadowing; 3) presence of exposed soil and; 4) non-photosynthetic vegetation material. Airborne, field and modeling studies have shown that the optical reflectance continuum (400-2500 nm) contains sufficient information to derive estimates of each of these indicators. Our ongoing studies in the eastern Amazon basin also suggest that these four indicators are sensitive to logging intensity. Satellite-based estimates of these indicators should provide a means to quantify both the presence and degree of structural disturbance caused by various logging regimes. Our quantitative assessment of Hyperion hyperspectral and ALI multi-spectral data for the detection and structural characterization of selective logging in Amazonia will benefit from data collected through an ongoing project run by the Tropical Forest Foundation, within which we have developed a study of the canopy and landscape biophysics of conventional and reduced-impact logging. We will add to our base of forest structural information in concert with an EO-1 overpass. Using a photon transport model inversion technique that accounts for non-linear mixing of the four biogeophysical indicators, we will estimate these parameters across a gradient of selective logging intensity provided by conventional and reduced impact logging sites. We will also compare our physical ly-based approach to both conventional (e.g., NDVI) and novel (e.g., SWIR-channel) vegetation indices as well as to linear mixture modeling methods. We will cross-compare these approaches using Hyperion and ALI imagers to determine the strengths and limitations of these two sensors for applications of forest biophysics. This effort will yield the first physical ly-based, quantitative analysis of the detection and intensity of selective logging in Amazonia, comparing hyperspectral and improved multi-spectral approaches as well as inverse modeling, linear mixture modeling, and vegetation index techniques.
Thermal diffusion in systems He4-HT and He3-HD analyzed by mixture of loaded and smooth spheres as model, discussing inelastic collisions
We present new optical (0.4-0.65 micron) spectra of Mercury and lunar pure anorthosite locations, obtained quasi-simultaneously with the Nordic Optical Telescope (NOT) in 2002. A comparative study is performed with the model of Lucey et al. between iron-poor, mature, pure anorthosite (less than 90% plagioclase feldspar) Clementine spectra from the lunar farside and a combined 0.4-1.0 micron mercurian spectrum, obtained with the NOT, calculated for standard photometric geometry. Mercury is located at more extreme locations in the Lucey ratio-reflectance diagrams than any known lunar soil, specifically with respect to the extremely iron-poor mature anorthosites. Though quantitative prediction of FeO and TiO2 abundances cannot be made without a more generally applicable model, we find qualitatively that the abundances of both these oxides must be near zero for Mercury. We utilize the theory of Hapke, with realistic photometric parameters, to model laboratory spectra of matured mineral powders and lunar soils, and remotely sensed spectra of lunar anorthosites and Mercury. An important difference between fabricated and natural powders is the high value for the internal scattering parameter necessary to interpret the spectra for the former, and the requirement of rough and non-isotropically scattering surfaces in the modelling of the latter. The mature lunar anorthosite spectra were well modelled with binary mixtures of calcic feldspars and olivines, grain sizes of 25-30 micron and a concentration of submicroscopic metallic iron (SMFe) of 0.12-0.15% in grain coatings. The mercurian spectrum is not possible to interpret from terrestrial mineral powder spectra without introducing an average particle scattering function for the bulk soil that increases in backscattering efficiency with wavelength. The observed spectrum is somewhat better predicted with binary mixture models of feldspars and pyroxenes, that single-component regoliths consisting of either albite or diopside. Correct spectral reflectance values were predicted with a concentration of 0.1 wt% SMFe in coatings of 15-30 micron sized grains. Since reasonable cosmogonical formation scenarios for Mercury, or meteoritic infall, predict iron concentrations at least this high, we draw the conclusion that the average grain size of Mercury is about a factor of two smaller than for average returned lunar soil samples. The 0.6-2.5 micron spectrum of McCord and Clark is used to further limit the possible range of mineralogical composition of Mercury. It is found that an intimately mixed and matured 3 : 1 labradorite-to-enstatite regolith composition best matches both the optical and near-infrared spectra, yielding an abundance of approx. 1.2 wt% FeO and -approx. 0 wt% TiO2.
Multiyear trends in Normalized Difference Vegetation Index (NDVI) have been used as metrics of high latitude ecosystem change based on the assumption that NDVI change is associated with ecological change, generally as changes in green vegetation amount (green leaf area index [LAI] or plant cover). Further, no change in NDVI is often interpreted as no change in these variables. Three canopy reflectance models including linear mixture model, the SAIL (Scattering from Arbitrarily Inclined Leaves) model, and the GeoSail model were used to simulate scenarios representing high latitude landscape NDVI responses to changes in LAI and plant cover. The simulations showed inconsistent NDVI responses. Clear increases in NDVI are generally associated with increases in LAI and plant cover. At higher values of LAI, the change in NDVI per unit change in LAI decreases, with very little change in spruce forest NDVI where crown cover is >50% and at the tundra–taiga ecotone with transitions from shrub tundra to spruce woodland. These lower responses may bias the interpretation of greening/browning trends in boreal forests. Variations in water or snow coverage were shown to produce outsized nonbiological NDVI responses. Inconsistencies in NDVI responses exemplify the need for care in the interpretation of NDVI change as a metric of high latitude ecosystem change, and that landscape characteristics in terms of the type of cover and its characteristics, such as the initial plant cover, must be taken into account in evaluating the significance of any observed NDVI trends.
Two methods for obtaining the required spectral signatures for a particular mixture model are considered. For the model considered, the spectral signatures become signature vectors. The first method is based upon determination of the signature vectors in such a way that a measure of the inconsistency between the mixture model and the observed data is minimized. The second method is based upon determination of the signature vectors in such a way that the estimated mean percentage coverage of individual species matches apriori or ground truth estimates. The two methods proposed are applied to actual multispectral data in order to verify the concepts presented.
Thermal deformations and stresses were studied in a silicon-carbide/aluminum filamentary composite at temperatures up to 370 C (700 F). Longitudinal and transverse thermal strains were measured with strain gages and a dilatometer. An elastoplastic micromechanical analysis based on a one-dimensional rule-of-mixtures model and an axisymmetric two-material composite cylinder model was performed. It was established that beyond a critical temperature thermal strains become nonlinear with decreasing longitudinal and increasing transverse thermal-expansion coefficients. This behavior was attributed to the plastic stresses in the aluminum matrix above the critical temperature. An elastoplastic analysis of both micromechanical models was performed to determine the stress distributions and thermal deformation in the fiber and matrix of the composite. While only axial stresses can be determined by the rule-of-mixtures model, the complete triaxial state of stress is established by the composite cylinder model. Theoretical predictions for the two thermal-expansion coefficients were in satisfactory agreement with experimental results.
The approach to predicting Landsat Multispectral Scanner System (MSS) endmember signatures from fraction images derived from high resolution Landsat Thematic Mapper (TM) data is presented. The purpose, conducted in the Mogi-Guacu study site located in Sao Paulo State, Brazil, was to determine if information derived from a mixture model applied to higher resolution TM data could serve as ground truth for the lower resolution MSS sensor. The Constrained Least Squares (CLS) method was used to generate vegetation, soil, and shade fraction images from the TM data. The resulting images were then used to estimate the endmember spectral response for MSS data by regressing each MSS spectral channel against the corresponding proportion values estimated for the same resolution cells from the TM mixture model. The evaluation of the predicted signatures was performed by comparing the corresponding fraction images derived from both the MSS and TM data acquired on 14 Sep. 1986. The techniques serves as a potential tool for integrating information in global studies where remote sensors with different spectral and spatial resolutions were used.
The state of the art was assessed with regards to efforts in support of the crop production estimation problem and alternative generic proportion estimation techniques were investigated. Topics covered include modeling the greeness profile (Badhwarmos model), parameter estimation using mixture models such as CLASSY, and minimum distance estimation as an alternative to maximum likelihood estimation. Approaches to the problem of obtaining proportion estimates when the underlying distributions are asymmetric are examined including the properties of Weibull distribution.
Semi-arid regions of the island of Hawaii and of the Tucson Mountains, AZ, having less than 50 percent green vegetation cover were mapped into rock/soil types by first identifying spectral mixtures of vegetation and rock/soil, and then removing the vegetation information. Landsat MSS images were used for Hawaii, and a Landsat TM image was used for the Tucson Mountains. Following analysis of the images field transects were done in test areas to determine percent vegetation cover and vegetation community characteristics, and to observe rock and soil types. For both Hawaii and the Tucson Mountains 20 percent or more of green vegetation caused potential misidentification of the rock/soil units unless spectral mixing was taken into account.
In-situ measurements of interplanetary coronal mass ejections (ICMEs) display a wide range of properties. A distinct subset, "magnetic clouds" (MCs), are readily identifiable by a smooth rotation in an enhanced magnetic field, together with an unusually low solar wind proton temperature. In this study, we analyze Ulysses spacecraft measurements to systematically investigate five possible explanations for why some ICMEs are observed to be MCs and others are not: i) An observational selection effect; that is, all ICMEs do in fact contain MCs, but the trajectory of the spacecraft through the ICME determines whether the MC is actually encountered; ii) interactions of an erupting flux rope (PR) with itself or between neighboring FRs, which produce complex structures in which the coherent magnetic structure has been destroyed; iii) an evolutionary process, such as relaxation to a low plasma-beta state that leads to the formation of an MC; iv) the existence of two (or more) intrinsic initiation mechanisms, some of which produce MCs and some that do not; or v) MCs are just an easily identifiable limit in an otherwise corntinuous spectrum of structures. We apply quantitative statistical models to assess these ideas. In particular, we use the Akaike information criterion (AIC) to rank the candidate models and a Gaussian mixture model (GMM) to uncover any intrinsic clustering of the data. Using a logistic regression, we find that plasma-beta, CME width, and the ratio O(sup 7) / O(sup 6) are the most significant predictor variables for the presence of an MC. Moreover, the propensity for an event to be identified as an MC decreases with heliocentric distance. These results tend to refute ideas ii) and iii). GMM clustering analysis further identifies three distinct groups of ICMEs; two of which match (at the 86% level) with events independently identified as MCs, and a third that matches with non-MCs (68 % overlap), Thus, idea v) is not supported. Choosing between ideas i) and iv) is more challenging, since they may effectively be indistinguishable from one another by a single in-situ spacecraft. We offer some suggestions on how future studies may address this.
In the present work a system of two-dimensional nonsteady hydrodynamic and chemical kinetic equations was numerically integrated for an exothermic system. Assumed two-step reaction model simulates practically an oxyhydrogen mixture. The calculation starts from a plane Chapman-Jouguet detonation as an initial condition. Two-dimensional disturbances are generated by artificially placing nonuniformities ahead of the detonation front. Regardless of the difference of the given initial disturbances, a fixed number of triple shock waves were produced for a fixed combination of mixture model and geometry when the transition period was over. This shows that for a given detonation tube geometry any exothermic system has its own characteristic multidimensional structure. The obtained number of triple shock waves contained in the detonation front was in agreement with existing experimental observations under the same condition.
A computational method based on the earlier work of one of the authors has been improved and modified to provide an accurate and rapid prediction of the radiating flow about blunt conical probes entering the atmosphere of Venus. Special features include a more general approximate equation of state, a more rapidly converging iteration scheme, a 38-step absorption coefficient model for a Venus atmosphere gas mixture model, and a shock-shape formula which yields a computed body shape in excellent agreement with the desired body. The method was used to calculate distributions of radiative heating for a large spherically blunted conical probe vehicle (1.219 m base diameter) and for three small spherically blunted conical probe vehicles (0.724 m base diameter). Solutions obtained with the presently developed approximate method required short computing times of the order of 1 to 3 minutes for each complete solution and gave agreements of 10% to 15% with the more exact predictions.