Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum entropy”

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 361 records · Page 20

Pixon deconvolution of far-infrared images from the UT multichannel photometer

Experiments with Pixon-based image deconvolution demonstrate that this technique can significantly enhance the spatial resolution of two-dimensional images acquired by the UT multichannel far-infrared photometer. A deconvolved 100 micron image of the giant HII region NGC 3603 reveals a bright, asymmetrical central peak some 40 x 55 arcsec in size, possibly surrounded by several secondary peaks. The results using the pixon technique are compared with those obtained using the standard Maximum Entropy and Richardson-Lucy methods.

Koresko, Chris D.↗

Microwave imaging of metal objects

The procedure of microwave imaging by maximum entropy method is discussed. First, the relationship between the induced current on the metal object surface and the scattered field is introduced. Our imaging concept is to reconstruct the induced current on the object surface from the measured scattered field. The object configuration will be provided by the induced current which is zero everywhere except on the object surface. Future work is also included with focus on the application of microwave imaging to both NASA and industry.

Krishen, Kumar↗

The source and distribution of Galactic positrons

The oriented scintillation spectrometer experiment (OSSE) observations of the Galactic plane and the Galactic center region were combined with observations acquired with other instruments in order to produce a map of the Galactic 511 keV annihilation radiation. Two mapping techniques were applied to the data: the maximum entropy method, and the basis pursuit inversion method. The resulting maps are qualitatively similar and show evidence for a central bulge and a weak galactic disk component. The weak disk is consistent with that expected from positrons produced by the decay of radioactive Al-26 in the interstellar medium. Both maps suggest an enhanced region of emission near l = -4 deg, b = 7 deg, with a flux of approximately 50 percent of that of the bulge. The existence of this emission appears significant, although the location is not well determined. The source of this enhanced emission is presently unknown.

Purcell, W. R.↗

Microwave Imaging on Metal Objects

This final report for the project discusses the attempts to model, using different methods, microwave image reconstruction. Maximum Entropy Method was not successful. Attempts to use Singular Value Decomposition (SVD) got some good results after initial failure. SVD is based upon a theory of linear algebra, to the effect that any M X N Matrix A whose number of rows M is greater than or equal to its number of columns, N can be written as the product of an M X N column-orthogonal matrix U, an N X N diagonal Matrix, W, with m positive or zero elements (the singular values) and the transposition of an N X N orthogonal matrix V. In microwave imaging, the scattered fields can be expressed by the induced current distribution. The SVD method required more contiguous computer memory than was available. Work was also done on the Conjugate Gradient Method (CGM), which didn't work well when tried earlier. It was found that separation of the imaginary part and the real part during calculation may work. This work was considered incomplete as of the end of the grant period.

Tolliver, C. L.↗

On the Spatial Distribution of High Velocity Al-26 Near the Galactic Center

We present results of simulations of the distribution of 1809 keV radiation from the decay of Al-26 in the Galaxy. Recent observations of this emission line using the Gamma Ray Imaging Spectrometer (GRIS) have indicated that the bulk of the AL-26 must have a velocity of approx. 500 km/ s. We have previously shown that a velocity this large could be maintained over the 10(exp 6) year lifetime of the Al-26 if it is trapped in dust grains that are reaccelerated periodically in the ISM. Here we investigate whether a dust grain velocity of approx. 500 km/ s will produce a distribution of 1809 keV emission in latitude that is consistent with the narrow distribution seen by COMPTEL. We find that dust grain velocities in the range 275 - 1000 km/ s are able to reproduce the COMPTEL 1809 keV emission maps reconstructed using the Richardson-Lucy and Maximum Entropy image reconstruction methods while the emission map reconstructed using the Multiresolution Regularized Expectation Maximization algorithm is not well fit by any of our models. The Al-26 production rate that is needed to reproduce the observed 1809 keV intensity yields in a Galactic mass of Al-26 of approx. 1.5 - 2 solar mass which is in good agreement with both other observations and theoretical production rates.

Sturner, Steven J.↗

Determining Biosignatures by Complexity Analysis in Antarctic Cryptoendolithic Communities

One of the most difficult problems of life detection is that of identifying biosignatures across a wide range of scales using multiple co-registered probes. The technique should be of equal utility across a wide range of search spaces from remote sensors probing volumes of space or planetary surfaces, visual eye or camera searches across the surface of a rock in Antarctica, low resolution microscopic scanning of a rock or a space craft in situ, or high resolution electron microscope and computerized tomography scanning of geobiological samples. We describe here an approach to this problem which derives in large part from past work done in the area of astrophysics - namely the analysis of complexity in galactic signals by data compression methods. This approach is a radically new one for geobiology and astrobiology, and allows us to assess the complexity (and thus potential biogenicity) of an object being examined. This is done by considering the information within pixels of an image (regardless the sensor used to gather the information) as an energetic system capable of description in terms of classical thermodynamics. The image data space is searched by an algorithm that judges complexity via data compression (e.g., the more compressible it is, the less complex, and vice versa) and maximum entropy as originally outlined by Shannon. At present we are implementing methods to utilize images from multiple sensors gathering different kinds of information (e.g., visible gray-scale data, color analyses, UV fluorescence, chemical information, etc). We present here preliminary data from deep UV fluorescence and ESEM (Environmental Scanning Electron Microscope) images from a layered cryptoendolithic community of an Antarctic rock.

Storrie-Lombardi, M. C.↗

Long-Slit Spectra of Mars in the Thermal Infrared

We have obtained long-slit spectra of Mars in the 7-14 micrometer regime. Our two integrations, taken 1993 March 2, are spaced about 2 hours apart and are centered at longitudes of roughly 210 and 240 degrees. The data cover the Elysium Planitia and volcanic regions, the classical dark albedo regions of Cerburus and central and eastern Hesperia, and the classical bright albedo region Utopia Planitia. We have deconvolved the point spread function from the images using a maximum entropy reconstruction algorithm; the resulting spatial resolution of the images is better than 1". Mars occupies approximately 9" of the slit, and we have detected variations in the spectral emission from different areas of the planet. We are in the process of interpreting these spectral variations in terms of the temperature and composition of the surface and the transmission of the Martian atmospheric dust and gases.

Rowland, C. M.↗

Long-slit Spectra of the Red Rectangle at 3 Microns

We have observed the Red Rectangle with the Cryogenic Spectrometer (CRSP) at the 2.1 m telescope of Kitt Peak National Observatory. The Red Rectangle is a well-known source of the infrared emission features associated with polycyclic aromatic hydrocarbons (PAHs). Our long-slit spectra cover the 3.2-3.6 micrometer regime, which includes the strong PAH band at 3.29 micrometers, a weaker component at 3.40 micrometers, and an emission plateau which continues out to 3.6 micrometers. Emission in the PAH bands extends several arcseconds to the north and south of the central source (HD 44179). The seeing during our observations was approximately 1.5 arcseconds, but we have deconvolved the point spread function from our spectral images at each wavelength using maximum entropy reconstruction, allowing us to examine the relative behavior of these emission components at finer spatial scales. With the new 256x256 array installed on CRSP, each pixel covers 0.6 arcseconds of sky. Our reconstructions reveal that the strength of the small 3.4 micrometer emission feature increases with respect to the 3.3 micrometer feature as the distance from the central source increases. In particular, there appears to be a region of enhanced 3.4 micrometer emission roughly 3" to the north and south of the central source. We are analyzing our results in terms of the competing models for the origin of the 3.4 micrometer emission feature, which could be a vibrational overtone of the 3.3 micrometer feature, or could result from the attachment of molecular sidegroups to the PAH molecule.

Sloan, G. C.↗

First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Foreground Emission

The WMAP mission has mapped the full sky to determine the geometry, content, and evolution of the universe. Full sky maps are made in five microwave frequency bands to separate the temperature anisotropy of the cosmic microwave background (CMB) from foreground emission, including diffuse Galactic emission and Galactic and extragalactic point sources. We define masks that excise regions of high foreground emission, so CMB analyses can became out with minimal foreground contamination. We also present maps and spectra of the individual emission components, leading to an improved understanding of Galactic astrophysical processes. The effectiveness of template fits to remove foreground emission from the WMAP data is also examined. These efforts result in a CMB map with minimal contamination and a demonstration that the WMAP CMB power spectrum is insensitive to residual foreground emission. We use a Maximum Entropy Method to construct a model of the Galactic emission components. The observed total Galactic emission matches the model to less than 1% and the individual model components are accurate to a few percent. We find that the Milky Way resembles other normal spiral galaxies between 408 MHz and 23 GHz, with a synchrotron spectral index that is flattest (beta(sub s) approx. -2.5) near star-forming regions, especially in the plane, and steepest (beta(sub s) approx. -3) in the halo. This is consistent with a picture of relativistic cosmic ray electron generation in star-forming regions and diffusion and convection within the plane. The significant synchrotron index steepening out of the plane suggests a diffusion process in which the halo electrons are trapped in the Galactic potential long enough to suffer synchrotron and inverse Compton energy losses and hence a spectral steepening. The synchrotron index is steeper in the WMAP bands than in lower frequency radio surveys, with a spectral break near 20 GHz to beta(sub s) less than -3. The modeled thermal dust spectral index is also steep in the WMAP bands, with beta(sub d) approx. = 2.2. Our model is driven to these conclusions by the low level of total foreground contamination at approx. 60 GHz. Microwave and Ha measurements of the ionized gas agree well with one another at about the expected levels. Spinning dust emission is limited to less than 5% of the Ka-band foreground emission. A catalog of 208 point sources is presented. The reliability of the catalog is 98%, i.e., we expect five of the 208 sources to be statistically spurious. The mean spectral index of the point sources is alpha approx. 0(beta approx. -2). Derived source counts suggest a contribution to the anisotropy power from unresolved sources of (15.0 +/- 1.4) x 10(exp -3)micro sq K sr at Q-band and negligible levels at V-band and W-band. The Sunyaev-Zeldovich effect is shown to be a negligible "contamination" to the maps.

Bennett, C. L.↗

Lightning and Other Influences On Tropical Tropospheric Ozone: Empirical Studies of Covariation

Tropical and subtropical tropospheric ozone are important radiatively active species, with particularly large effects in the upper third of the troposphere. Temporal variability of O3 has proved difficult to simulate day by day in process models. Thus, individual roles of lightning, biomass burning, and other pollution in providing precursor NO(x), radicals, and chain carriers (CO, hydrocarbons) remain unquantified by simulation, and it is theoretically reasonable that individual roles are magnified by a joint synergy. We use wavelet analysis and Burg-algorithm maximum entropy spectral analyses to describe time-scales and correlation of ozone with proxies for processes controlling its concentration. Our empirical studies link time variations apparent in several datasets: the SHADOZ (Southern Hemisphere Additional Ozonesondes) network stations (Nairobi, Fiji), and auxiliary series with power to explain ozone-determining processes, with some interpretation based on the TTO (Tropical Tropospheric Ozone) product derived from TOMS (the Total Ozone Mapping Spectrometer). The auxiliary series are The OTD/LIS(Optical Transient Detector/Lightning Imaging Sensor) measurements of the lightning NO(x) source, the OLR (Outgoing Longwave Radiation)measurement of high-topped clouds, and standard meteorological variables from the United States NCEP (National Centers for Environmental Prediction) and Data Assimilation Office analyses. Concentrating on equatorial ozone, we compare the statistical evidence on the variability of tropospheric ozone. Important variations occur on approximately two-week, two-month (Madden-Julian Oscillation) and annual scales, and relations with OLR suggest controls associated with continental clouds. Hence we are now using the Lightning Imaging Sensor data set to indicate NO(x) sources. We report initial results defining relative roles of the process mentioned affecting O3 using their covariance properties.

Chatfield, Robert B.↗

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

Spray Characteristics of a Hybrid Twin-Fluid Pressure-Swirl Atomizer

The spray performance of a fuel injection system applicable for use in main combustion chamber of an oxidizer-rich staged combustion (ORSC) cycles is presented. The experimental data reported here include mean drop size and drop size distribution, spray cone half-angle, and momentum rate (directly related to spray penetration). The maximum entropy formalism, MEF, method to predict drop size distribution is applied and compared to the experimental data. Geometric variables considered include the radius of the injector inlet orifice plate through which oxidizer flows (&) and the exposed length from the fuel inlet to the injector exit plane (L2). Operating conditions that were varied include the liquid mass flow rate and air mass flow rate. For orifices B and C there is a significant dependence of D3Z on both the air and liquid mass flow rates, as well as on L2. For the A orifice, the momentum rate of the air flow appears to exceed a threshold value above which a constant D32 is obtained. Using the MEF method, a semi-analytical process was developed to model the spray distribution using two input parameters (q = 0.4 and Dso). The momentum rate of the spray is directly related to the air and liquid mass flow rates. The cone half angle of the spray ranges from 25 to 17 degrees. The data resulting from this project will eventually be used to develop advanced rocket systems.

Durham, M. J.↗

Irregular Wavelike Structure in Saturn's Rings

We have searched Saturn's A, B, and C rings for irregular wavelike structure using Voyager Photopolarimeter (PPS), Ultraviolet Spectrometer (UVS), and Radio Science (RSS) occultation datasets, as well as ring reflectivity profiles derived from Voyager images. A maximum entropy technique for conducting spectral analysis was used to estimate wave frequency power in relation to radial location for each dataset. Using this method we have found irregular structure in the PPS and UVS inner B Ring occultation datasets previously identified in Voyager imaging data. Both finer structure, with a wavelength of around 20 km, and large structure with wavelengths of 200 to 1000 km, are visible in the occultation data and appear similar to that seen in the imaging data. After removing ringlets from the C-Ring data, we have identified what appears to be a 1000-km wave sustained throughout the ring. The large dominant wavelength appears in all datasets; however, tests are currently being conducted in an attempt to verify its existence. Irregular structure with a wavelength of approximately 20 km has been observed in the C Ring reflectivity profiles, but not within the occultation datasets. This leads us to doubt it is caused by ring surface mass density fluctuations detectable by the occultation experiments.

Pollard, Benjamin J.↗

The El Nino/Southern Oscillation and Future Soybean Prices

Recently, it was shown that the application of a method combining singular spectrum analysis (SSA) and the maximum entropy method to univariate indicators of the coupled ocean-atmosphere El Nino/Southern Oscillation (ENSO) phenomenon can be helpful in determining whether an El Nino (EN) or La Nina (LN) event will occur. SSA - a variant of principal component analysis applied in the time domain - filters out variability unrelated to ENSO and separates the quasi-biennial (QB), two-to-three year variability, from a lower-frequency (LF) four-to-six year EN-LN cycle; the total variance associated with ENSO combines the QB and LF modes. ENSO has been known to affect weather conditions over much of the globe. For example, EN events have been connected with unusually rainy weather over the Central and Western US, while the opposite phases of the oscillation (LN) have been plausibly associated with extreme dry conditions over much of the same geographical area...

Keppenne, C.↗

Space Invariant Independent Component Analysis and ENose for Detection of Selective Chemicals in an Unknown Environment

In this paper, we present a space invariant architecture to enable the Independent Component Analysis (ICA) to solve chemical detection from two unknown mixing chemical sources. The two sets of unknown paired mixture sources are collected via JPL 16-ENose sensor array in the unknown environment with, at most, 12 samples data collected. Our space invariant architecture along with the maximum entropy information technique by Bell and Sejnowski and natural gradient descent by Amari has demonstrated that it is effective to separate the two mixing unknown chemical sources with unknown mixing levels to the array of two original sources under insufficient sampled data. From separated sources, they can be identified by projecting them on the 11 known chemical sources to find the best match for detection. We also present the results of our simulations. These simulations have shown that 100% correct detection could be achieved under the two cases: a) under-completed case where the number of input (mixtures) is larger than number of original chemical sources; and b) regular case where the number of input is as the same as the number of sources while the time invariant architecture approach may face the obstacles: overcomplete case, insufficient data and cumbersome architecture.

ENose↗

Granular Contact Forces: Proof of "Self-Ergodicity" by Generalizing Boltzmann's Stosszahlansatz and H Theorem

Ergodicity is proved for granular contact forces. To obtain this proof from first principles, this paper generalizes Boltzmann's stosszahlansatz (molecular chaos) so that it maintains the necessary correlations and symmetries of granular packing ensembles. Then it formally counts granular contact force states and thereby defines the proper analog of Boltzmann's H functional. This functional is used to prove that (essentially) all static granular packings must exist at maximum entropy with respect to their contact forces. Therefore, the propagation of granular contact forces through a packing is a truly ergodic process in the Boltzmannian sense, or better, it is self-ergodic. Self-ergodicity refers to the non-dynamic, internal relationships that exist between the layer-by-layer and column-by-column subspaces contained within the phase space locus of any particular granular packing microstate. The generalized H Theorem also produces a recursion equation that may be solved numerically to obtain the density of single particle states and hence the distribution of granular contact forces corresponding to the condition of self-ergodicity. The predictions of the theory are overwhelmingly validated by comparison to empirical data from discrete element modeling.

Metzger, Philip T.↗

The Invasive Species Forecasting System (ISFS): An iRODS-Based, Cloud-Enabled Decision Support System for Invasive Species Habitat Suitability Modeling

The Invasive Species Forecasting System (ISFS) is an online decision support system that allows users to load point occurrence field sample data for a plant species of interest and quickly generate habitat suitability maps for geographic regions of interest, such as a national park, monument, forest, or refuge. Target customers for ISFS are natural resource managers and decision makers who have a need for scientifically valid, model- based predictions of the habitat suitability of plant species of management concern. In a joint project involving NASA and the Maryland Department of Natural Resources, ISFS has been used to model the potential distribution of Wavyleaf Basketgrass in Maryland's Chesapeake Bay Watershed. Maximum entropy techniques are used to generate predictive maps using predictor datasets derived from remotely sensed data and climate simulation outputs. The workflow to run a model is implemented in an iRODS microservice using a custom ISFS file driver that clips and re-projects data to geographic regions of interest, then shells out to perform MaxEnt processing on the input data. When the model completes, all output files and maps from the model run are registered in iRODS and made accessible to the user. The ISFS user interface is a web browser that uses the iRODS PHP client to interact with the ISFS/iRODS- server. ISFS is designed to reside in a VMware virtual machine running SLES 11 and iRODS 3.0. The ISFS virtual machine is hosted in a VMware vSphere private cloud infrastructure to deliver the online service.

Gill, Roger↗

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗