Search NASASearch

SEARCH · Search NASA

Results for “Uniform Methods Project”

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.

70 records · Page 4

Predicting the Seawater Chemistry of an Ocean World Using Machine Learning on Isotopic Measurements of Volatile CO2

Introduction: Given the long time intervals required for data transmission to and from ocean worlds targets, low bandwidth for data transmission, time required for data processing and analysis, and potentially extreme radiation environments (e.g., Europa), it is clear that ocean worlds missions will need more autonomous flight instruments and software in order to achieve established science goals. Protracted time intervals for data analysis (e.g., Europa Lander) strongly motivates the development of rapid, consistent and streamlined methods for interpreting data from flight mass spectrometers to e.g., determine how mass spectra from a plume or surface liquid/ice relates to the surface/subsurface. Since mass spectrometry also has the potential to correctly identify biosignatures[1], it is imperative that such methods for interpreting data are consistent and accurate. We used 848 isotope ratio mass spectra from laboratory analyses of CO2 that interacted with ocean worlds-relevant seawaters as a ‘training’ dataset for ‘unsupervised’ machine learning. In unsupervised learning, characteristics of the data are not labeled or linked, and any similarities found only result from the neural network. CO2 isotopologues analyzed for this dataset mimic the remote measurements of CO2 by a flight mass spectrometer, and are detailed in Theiling [2]. From this dataset, we used measured features of the spectra, such as retention time, intensity, and (isotopologue) mass ratios as inputs for our autoencoder neural network. Our neural network was trained to find similarities in these and other spectral features for seawaters of a particular composition and amount of initial CO2. Successful training then created an output of these similarities for various seawaters, which included MgSO4, Na2SO4, NaCl, MgCl2, KCl, and NaHCO3, and combinations of these salts. We then applied dimensionality reduction techniques such as Principal Component Analysis (PCA), T-Distributed Stochastic Neighbor Embedding (TSNE), and Uniform Manifold Approximation and Projection (UMAP) to demonstrate latent data features as a two-dimensional projection in a unitless, high-dimensional space. In this projection, a data point represents the combined effect of spectral features such as intensity, retention time, and isotope ratio. Our initial UMAP demonstrates data clustering (organization of the data by the neural network) based on the amount of CO2 that had initially interacted with each seawater. Further training using more ‘supervised’ learning techniques demonstrate strong clustering of preliminary data based on initial CO2 concentration, seawater chemical composition, and ionic strength (salinity). Our preliminary work therefore suggests that machine learning has the potential to identify compositional variants of an ocean world seawater based on mass spectra from volatile CO2 measurements. Acknowledgments: This work was funded through a Strategic Task Group at NASA Goddard Space Flight Center. The training dataset was collected through funding from the Oklahoma Space Grant Consortium. References: [1] Pappalardo, R. et al. (2013) Astrobiology, 13, 740–773. [2] Theiling (2020) Icarus, 114216.

Europa

SO2 on Venus: A final cross-calibration with Pioneer Venus

In our previous report, we noted that we had successfully completed all of the IUE observations which were proposed. A total of 7 Venus spectra were obtained. In the present reporting period, essentially all of the associated data reduction and analysis activities were completed, and Drs. Stern, Barker, and Na met to collaborate on data interpretation. The model atmosphere used in this analysis is similar to the one used in the analysis of the Pioneer Venus observations (Esposito, et al. 1979, 1988) and the previous IUE observations (Na, et al. 1990). In this model, sulfuric acid aerosols with a radius of about 1 micron are mixed uniformly with Rayleigh scattering gases throughout the cloud layer. There are two wavelength-independent absorbers. First, the aerosols are assumed to have a single scattering albedo of 0.98, and second a pure absorbing layer with an optical depth of 0.2 is added at an altitude of 75 mbar. The variable parameters in this model are the mixing ratios of SO2 at the cloud top (40 mbar level), and the scale height of SO2 at the same altitude. Model spectra were calculated using a radiative transfer code utilizing a Markov Chain method. Each model calculation takes into account the multiple scattering and the vertical inhomogeneity of the atmosphere. We reported those initial results at the 1992 meeting of the AAS Division of Planetary Sciences (Barker, Stern, and Na 1992) The important finding we have reached is that the 1991 IUE observations indicate that the long-term decline in Venus SO2 abundance has either halted or reversed. Our simultaneous groundbased measurements made on the McDonald Observatory 2.7 m corroborate this result. We are now preparing a paper describing these IUE results for publication in Icarus. With the IUE results reduced, we are also beginning a project to compare the 1991 IUE results to Pioneer Venus UV spectra made on the same dates.

Stern, S. Alan

Study of the Radiative Properties of Inhomogeneous Stratocumulus Clouds

Clouds play an important role in the radiation budget of the atmosphere. A good understanding of how clouds interact with solar radiation is necessary when considering their effects in both general circulation models and climate models. This study examined the radiative properties of clouds in both an inhomogeneous cloud system, and a simplified cloud system through the use of a Monte Carlo model. The purpose was to become more familiar with the radiative properties of clouds, especially absorption, and to investigate the excess absorption of solar radiation from observations over that calculated from theory. The first cloud system indicated that the absorptance actually decreased as the cloud's inhomogeneity increased, and that cloud forcing does not indicate any changes. The simplified cloud system looked at two different cases of absorption of solar radiation in the cloud. The absorptances calculated from the Monte Carlo is compared to a correction method for calculating absorptances and found that the method can over or underestimate absorptances at cloud edges. Also the cloud edge effects due to solar radiation points to a possibility of overestimating the retrieved optical depth at the edge, and indicates a possible way to correct for it. The effective cloud fraction (Ne) for a long time has been calculated from a cloud's reflectance. From the reflectance it has been observed that the N, for most cloud geometries is greater than the actual cloud fraction (Nc) making a cloud appear wider than it is optically. Recent studies we have performed used a Monte Carlo model to calculate the N, of a cloud using not only the reflectance but also the absorptance. The derived Ne's from the absorptance in some of the Monte Carlo runs did not give the same results as derived from the reflectance. This study also examined the inhomogeneity of clouds to find a relationship between larger and smaller scales, or wavelengths, of the cloud. Both Fourier transforms and wavelet transforms were used to analyze the liquid water content of marine stratocumulus clouds taken during the ASTEX project. From the analysis it was found that the energy in the cloud is not uniformly distributed but is greater at the larger scales than at the smaller scales. This was determined by examining the slope of the power spectrum, and by comparing the variability at two scales from a wavelet analysis.

Batey, Michael

Embracing Open Source for NASA's Earth Science Data Systems

The overarching purpose of NASAs Earth Science program is to develop a scientific understanding of Earth as a system. Scientific knowledge is most robust and actionable when resulting from transparent, traceable, and reproducible methods. Reproducibility includes open access to the data as well as the software used to arrive at results. Additionally, software that is custom-developed for NASA should be open to the greatest degree possible, to enable re-use across Federal agencies, reduce overall costs to the government, remove barriers to innovation, and promote consistency through the use of uniform standards. Finally, Open Source Software (OSS) practices facilitate collaboration between agencies and the private sector. To best meet these ends, NASAs Earth Science Division promotes the full and open sharing of not only all data, metadata, products, information, documentation, models, images, and research results but also the source code used to generate, manipulate and analyze them. This talk focuses on the challenges to open sourcing NASA developed software within ESD and the growing pains associated with establishing policies running the gamut of tracking issues, properly documenting build processes, engaging the open source community, maintaining internal compliance, and accepting contributions from external sources. This talk also covers the adoption of existing open source technologies and standards to enhance our custom solutions and our contributions back to the community. Finally, we will be introducing the most recent OSS contributions from NASA Earth Science program and promoting these projects for wider community review and adoption.

Earth Science

A Regional Analysis of Cloudy Mean Spherical Albedo over the Marine Stratocumulus Region and the Tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. This study examines the relationships between cloud properties and cloud fraction in order to supplement grid scale parameterizations. The satellite data used in this study is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.50 x 2.50 latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakleys (1991) theory of reflection for uniform and broken layered cloud and to Kedem, et al.(1990) findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 in) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al. theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.

A regional analysis of cloudy mean spherical albedo over the marine stratocumulus region and the tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. The relationships between cloud properties and cloud fraction are studied in order to supplement grid scale parameterizations. The satellite data used is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.5 deg x 2.5 deg latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution, and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakley's theory of reflection for uniform and broken layered cloud and to Kedem, et al.'s findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 m) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al's theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.

Constructing LDPC Codes from Loop-Free Encoding Modules

A method of constructing certain low-density parity-check (LDPC) codes by use of relatively simple loop-free coding modules has been developed. The subclasses of LDPC codes to which the method applies includes accumulate-repeat-accumulate (ARA) codes, accumulate-repeat-check-accumulate codes, and the codes described in Accumulate-Repeat-Accumulate-Accumulate Codes (NPO-41305), NASA Tech Briefs, Vol. 31, No. 9 (September 2007), page 90. All of the affected codes can be characterized as serial/parallel (hybrid) concatenations of such relatively simple modules as accumulators, repetition codes, differentiators, and punctured single-parity check codes. These are error-correcting codes suitable for use in a variety of wireless data-communication systems that include noisy channels. These codes can also be characterized as hybrid turbolike codes that have projected graph or protograph representations (for example see figure); these characteristics make it possible to design high-speed iterative decoders that utilize belief-propagation algorithms. The present method comprises two related submethods for constructing LDPC codes from simple loop-free modules with circulant permutations. The first submethod is an iterative encoding method based on the erasure-decoding algorithm. The computations required by this method are well organized because they involve a parity-check matrix having a block-circulant structure. The second submethod involves the use of block-circulant generator matrices. The encoders of this method are very similar to those of recursive convolutional codes. Some encoders according to this second submethod have been implemented in a small field-programmable gate array that operates at a speed of 100 megasymbols per second. By use of density evolution (a computational- simulation technique for analyzing performances of LDPC codes), it has been shown through some examples that as the block size goes to infinity, low iterative decoding thresholds close to channel capacity limits can be achieved for the codes of the type in question having low maximum variable node degrees. The decoding thresholds in these examples are lower than those of the best-known unstructured irregular LDPC codes constrained to have the same maximum node degrees. Furthermore, the present method enables the construction of codes of any desired rate with thresholds that stay uniformly close to their respective channel capacity thresholds.

Divsalar, Dariush

Techniques Optimized for Reducing Instabilities in Advanced Nickel-Base Superalloys for Turbine Blades

The High-Speed Research (HSR) Airfoil Alloy program developed fourth-generation single-crystal superalloys with up to an 85 F increase in creep rupture capability over current production airfoil alloys. Recent results have been generated at the NASA Glenn Research Center on these fourth-generation alloys, but in coated form, for subsonic turbine blade applications under NASA's Ultra-Efficient Engine Technology (UEET) Program. One goal for UEET is to optimize the airfoil alloy/thermal barrier coating system for 3100 F turbine inlet temperatures. The state-of-the art turbine blade airfoil system consists of a superalloy single crystal that provides the basic mechanical performance of the airfoil. A thermal barrier coating is used to reduce the temperature of the base superalloy, and a bondcoat is deposited between the base material and the thermal barrier coating. The bondcoat improves the oxidation and corrosion resistance of the base superalloy and improves the spallation resistance of the thermal barrier coating. A commercial platinum aluminide bondcoat was applied to the HSR-developed alloys, and a diffusion zone developed as a result of interaction between the bondcoat and the superalloy. Optimized strength is obtained for superalloys when the refractory element content is high and the limits of microstructural stability are approached or exceeded slightly. For fourthgeneration alloys, instability leads to the formation of topologically close packed (TCP) phases, which form internally in the superalloy, and a secondary reaction zone (SRZ), which forms under the diffusion zone. There was a concern that excessive quantities of either TCP or SRZ might decrease the mechanical properties of the superalloy, with SRZ thought to be particularly detrimental and its formation unpredictable. Thus, an SRZreduction effort was initiated in the NASA UEET Program so that methods developed during the HSR project could be optimized further to reduce or eliminate the SRZ. An SRZ is a three-phase constituent composed of TCP and stringers of gamma phase in a matrix of gamma prime. An incoherent grain boundary separates the SRZ from the gammagamma prime microstructure of the superalloy. The SRZ is believed to form as a result of local chemistry changes in the superalloy due to the application of the diffusion aluminide bondcoat. Locally high surface stresses also appear to promote the formation of the SRZ. Thus, techniques that change the local alloy chemistry or reduce surface stresses have been examined for their effectiveness in reducing SRZ. These SRZ-reduction steps are performed on the test specimen or the turbine blade before the bondcoat is applied. Stressrelief heat treatments developed at NASA Glenn have been demonstrated to reduce significantly the amount of SRZ that develops during subsequent high-temperature exposures. Stress-relief heat treatments reduce surface stresses by recrystallizing a thin surface layer of the superalloy. However, in alloys with very high propensities to form SRZ, stress relief heat treatments alone do not eliminate SRZ entirely. Thus, techniques that modify the local chemistry under the bondcoat have been emphasized and optimized successfully at Glenn. One such technique is carburization, which changes the local chemistry by forming submicron carbides near the surface of the superalloy. Detailed characterizations have demonstrated that the depth and uniform distribution of these carbides are enhanced when a stress relief treatment and an appropriate surface preparation are employed in advance of the carburization treatment. Even in alloys that have the propensity to develop a continuous SRZ layer beneath the diffusion zone, the SRZ has been completely eliminated or reduced to low, manageable levels when this combination of techniques is utilized. Now that the techniques to mitigate SRZ have been established at Glenn, TCP phase formation is being emphasized in ongoing work under the UEET Program. The limitsf stability of the fourth-generation alloys with respect to TCP phase formation are currently being defined along with high-temperature creep rupture properties. In addition, a regression model is being developed at Glenn for the prediction of the presence of TCP phase in the microstructure and SRZ under the diffusion zone. The model is based on a design-of-experiments methodology with emphasis on the potential synergistic effects of alloying elements.

MacKay, Rebecca A.

Windowed and Wavelet Analysis of Marine Stratocumulus Cloud Inhomogeneity

To improve radiative transfer calculations for inhomogeneous clouds, a consistent means of modeling inhomogeneity is needed. One current method of modeling cloud inhomogeneity is through the use of fractal parameters. This method is based on the supposition that cloud inhomogeneity over a large range of scales is related. An analysis technique named wavelet analysis provides a means of studying the multiscale nature of cloud inhomogeneity. In this paper, the authors discuss the analysis and modeling of cloud inhomogeneity through the use of wavelet analysis. Wavelet analysis as well as other windowed analysis techniques are used to study liquid water path (LWP) measurements obtained during the marine stratocumulus phase of the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment. Statistics obtained using analysis windows, which are translated to span the LWP dataset, are used to study the local (small scale) properties of the cloud field as well as their time dependence. The LWP data are transformed onto an orthogonal wavelet basis that represents the data as a number of times series. Each of these time series lies within a frequency band and has a mean frequency that is half the frequency of the previous band. Wavelet analysis combined with translated analysis windows reveals that the local standard deviation of each frequency band is correlated with the local standard deviation of the other frequency bands. The ratio between the standard deviation of adjacent frequency bands is 0.9 and remains constant with respect to time. This ratio defined as the variance coupling parameter is applicable to all of the frequency bands studied and appears to be related to the slope of the data's power spectrum. Similar analyses are performed on two cloud inhomogeneity models, which use fractal-based concepts to introduce inhomogeneity into a uniform cloud field. The bounded cascade model does this by iteratively redistributing LWP at each scale using the value of the local mean. This model is reformulated into a wavelet multiresolution framework, thereby presenting a number of variants of the bounded cascade model. One variant introduced in this paper is the 'variance coupled model,' which redistributes LWP using the local standard deviation and the variance coupling parameter. While the bounded cascade model provides an elegant two- parameter model for generating cloud inhomogeneity, the multiresolution framework provides more flexibility at the expense of model complexity. Comparisons are made with the results from the LWP data analysis to demonstrate both the strengths and weaknesses of these models.

Gollmer, Steven M.

Aerodynamic sound of flow past an airfoil

The long term objective of this project is to develop a computational method for predicting the noise of turbulence-airfoil interactions, particularly at the trailing edge. We seek to obtain the energy-containing features of the turbulent boundary layers and the near-wake using Navier-Stokes Simulation (LES or DNS), and then to calculate the far-field acoustic characteristics by means of acoustic analogy theories, using the simulation data as acoustic source functions. Two distinct types of noise can be emitted from airfoil trailing edges. The first, a tonal or narrowband sound caused by vortex shedding, is normally associated with blunt trailing edges, high angles of attack, or laminar flow airfoils. The second source is of broadband nature arising from the aeroacoustic scattering of turbulent eddies by the trailing edge. Due to its importance to airframe noise, rotor and propeller noise, etc., trailing edge noise has been the subject of extensive theoretical (e.g. Crighton & Leppington 1971; Howe 1978) as well as experimental investigations (e.g. Brooks & Hodgson 1981; Blake & Gershfeld 1988). A number of challenges exist concerning acoustic analogy based noise computations. These include the elimination of spurious sound caused by vortices crossing permeable computational boundaries in the wake, the treatment of noncompact source regions, and the accurate description of wave reflection by the solid surface and scattering near the edge. In addition, accurate turbulence statistics in the flow field are required for the evaluation of acoustic source functions. Major efforts to date have been focused on the first two challenges. To this end, a paradigm problem of laminar vortex shedding, generated by a two dimensional, uniform stream past a NACA0012 airfoil, is used to address the relevant numerical issues. Under the low Mach number approximation, the near-field flow quantities are obtained by solving the incompressible Navier-Stokes equations numerically at chord Reynolds number of 104. The far-field noise is computed using Curle's extension to the Lighthill analogy (Curle 1955). An effective method for separating the physical noise source from spurious boundary contributions is developed. This allows an accurate evaluation of the Reynolds stress volume quadrupoles, in addition to the more readily computable surface dipoles due to the unsteady lift and drag. The effect of noncompact source distribution on the far-field sound is assessed using an efficient integration scheme for the Curle integral, with full account of retarded-time variations. The numerical results confirm in quantitative terms that the far-field sound is dominated by the surface pressure dipoles at low Mach number. The techniques developed are applicable to a wide range of flows, including jets and mixing layers, where the Reynolds stress quadrupoles play a prominent or even dominant role in the overall sound generation.

Wang, Meng

MANN: A program to transfer designs for diffractive optical elements to a MANN photolithographic mask generator

There are two basic areas of interest for diffractive optics. In the first, the property of wavefront division is exploited for achieving optical fanout, analogous to the more familiar electrical fanout of electronic circuitry. The basic problem here is that when using a simple uniform diffraction grating the energy input is divided unevenly among the output beams. The other area of interest is the use of diffractive elements to replace or supplement standard refractive elements such as lenses. Again, local grating variations can be used to control the amount of bending imparted to optical rays, and the efficiency of the diffractive element will depend on how closely the element can be matched to the design requirements. In general, production restrictions limit how closely the element approaches the design, and for the common case of photolithographic production, a series of binary masks is required to achieve high efficiency. The actual design process is much more involved than in the case of elements for optical fanout, as the desired phase of the optical wavefront over some reference plane must be specified and the phase alteration to be introduced at each point by the diffraction element must be known. This generally requires the utilization of a standard optical design program. Two approaches are possible. In the first approach, the diffractive element is treated as a special type of lens and the ordinary optical design equations are used. Optical design programs tend to follow a second approach, namely, using the equations of optical interference derived from holographic theory and then allowing the introduction of phase front corrections in the form of polynomial equations. By using either of these two methods, diffractive elements can be used not only to compensate for distortions such as chromatic or spherical aberration, but also to perform the work of a variety of other optical elements such as null correctors, beam shapers, etc. The main focus of the project described in this report is how the design information from the lens design program is incorporated into the photolithographic process. It is shown that the MANN program, a photolithographic mask generator, fills the need for a link between lens design programs and mask generation controllers.The generated masks can be used to expose a resist-coated substrate which is etched and then must be re-coated, re-exposed, and re-etched for making copies, just as in the electronics industry.

Matthys, Donald R.

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.

Methods of Fabricating Scintillators with Radioisotopes for Beta Battery Applications

Technology has been developed for a class of self-contained, long-duration power sources called beta batteries, which harvest the energy contained in the radioactive emissions from beta decay isotopes. The new battery is a significant improvement over the conventional phosphor/solar cell concept for converting this energy in three ways. First, the thin phosphor is replaced with a thick scintillator that is transparent to its own emissions. By using a scintillator sufficiently thick to completely stop all the beta particles, efficiency is greatly improved. Second, since the energy of the beta particles is absorbed in the scintillator, the semiconductor photodetector is shielded from radiation damage that presently limits the performance and lifetime of traditional phosphor converters. Finally, instead of a thin film of beta-emitting material, the isotopes are incorporated into the entire volume of the thick scintillator crystal allowing more activity to be included in the converter without self-absorption. There is no chemical difference between radioactive and stable strontium beta emitters such as Sr-90, so the beta emitter can be uniformly distributed throughout a strontium based scintillator crystal. When beta emitter material is applied as a foil or thin film to the surface of a solar cell or even to the surface of a scintillator, much of the radiation escapes due to the geometry, and some is absorbed within the layer itself, leading to inefficient harvesting of the energy. In contrast, if the emitting atoms are incorporated within the scintillator, the geometry allows for the capture and efficient conversion of the energy of particles emitted in any direction. Any gamma rays associated with secondary decays or Bremsstrahlung photons may also be absorbed within the scintillator, and converted to lower energy photons, which will in turn be captured by the photocell or photodiode. Some energy will be lost in this two-stage conversion process (high-energy particle to low-energy photons to electric current). The geometric advantage partially offsets this as well, since the absorption depth of high-energy beta radiation is much larger than the depth of a p-n junction. Thus, in a p-n junction device, much of the radiation is absorbed far away from the junction, and the electron- hole pairs are not all effectively collected. In contrast, with a transparent scintillator the radiation can be converted to light in a larger volume, and all of the light can be collected in the active region of the photodiode. Finally, the new device is more practical because it can be used at much higher power levels without unduly shortening its lifetime. While the crystal structure of scintillators is also subject to radiation damage, their performance is far more tolerant of defects than that of semiconductor junctions. This allows the scintillator- based approach to use both higher energy isotopes and larger quantities of the isotopes. It is projected that this technology has the potential to produce a radioisotope battery with up to twice the efficiency of presently used systems.

Rensing, Noa M.

Metallic Powder Core Tubular Wire (PCTW) Development for Additive Manufacturing

NASA Langley Research Center (LaRC) in Hampton, VA, and the Center for Welding, Joining, and Coating Research (CWJCR) in the Department of Metallurgical and Materials Engineering at the Colorado School of Mines (CSM) in Golden, CO, have collaborated for nearly fifteen years to develop metallic powder core tubular wire (PCTW) feedstock for additive manufacturing (AM) using electron beam deposition. The motivation for this development was to offset dealloying observed during deposition of titanium and aluminum alloy solid wire feedstock in the LaRC electron beam freeform fabrication (EBF3) system and to explore novel alloy compositions by tailoring the powder fill. EBF3 is a fusion-based AM process that is performed under vacuum and is consequently prone to vaporization loss of low vapor pressure alloying elements. Research spanning four graduate student research projects succeeded in developing the basic methodology to fabricate PCTW, the mass balance equations used to define the powder fill, microstructure control through particulate inoculation, and the feasibility of producing metal matrix composite materials. Metallic PCTW is comprised of a thin foil sheath wrapped around a powder blend. Fabrication of PCTW at CSM was accomplished using a tubular wire mill system to first bend a flat strip into a u-shaped cross section, deposit a powder fill, create a lap joint to seal the wire, and draw to final diameter. Powder management systems were developed to ensure uniform mixing of powder constituents and identify feed rates for optimum fill ratios. Processing parameters were successfully developed for titanium and aluminum alloy sheaths to account for material work hardening during drawing. Mass balance equations were developed to account for compositional contributions of the sheath and powder fill. PCTW composition can be adjusted through selection of the sheath or powder materials. PCTWs were successfully developed to offset Al loss in deposits of Ti-6-4 and Mg loss in Al 6061, resulting in deposits that were within the composition limits for each material. Ti-6-4 deposits were fabricated using PCTW comprised of commercially pure (CP) titanium sheath and a fill of pre-alloyed Ti-6-4 powder with additions of elemental Al and V powders. Additions of both elements were required to account for the CP Ti from the sheath, and Al was further enhanced to offset documented losses in EBF3 deposits. For Al 6061 deposits, composition was modulated through the choice of sheath material. The PCTW was made using an Al 5052 sheath to supplement Mg content and a fill composition of pre-alloyed 6061 powder with additions of Si and Cu.The slower solidification rate of the EBF3 process compared with laser powder bed fusion processes can lead to coarse-grained microstructures and mechanical property anisotropy. EBF3deposits of Ti-6-4 exhibited large epitaxially grown β grains, which result in degraded mechanical properties in the deposition direction. Microstructure refinement was achieved through a combination of electron beam modulation and the use of a PCTW designed for Ti-6-4 deposits and modified by the addition of powder made from Fe and B. Resulting deposits exhibited reductions in β grain size, α phase colony intercept length and lath width, and an increased volume fraction of β phase. The mechanism of microstructure refinement was related to precipitation of boride particles. Hardness testing confirmed that mechanical property anisotropy was reduced and an overall increase in hardness was achieved.Metal matrix composite (MMC) materials offer strength and stiffness benefits over unreinforced materials. Two methods were investigated using PCTW technology for producing Al 6061 basedMMC materials; the direct addition of SiC reinforcing particles and the in-situ generation of reinforcements through reaction synthesis of precursor powders. Al 6061 based MMCs reinforced with SiC exhibited particle clustering and formation of aluminum carbides during deposition. The application of a Ni coating to the SiC particles prior to incorporation into the powder fill eliminated aluminum carbide formation and aided particle distribution during deposition. Modulation of beam focus and power further reduced aluminum carbide formation and aided particle distribution. In reaction synthesis, precursor additions undergo an exothermic reaction during deposition to form new ceramic products, generally of finer scale than the initial precursor components. The product phases provide reinforcement in the MMC and act as nucleation sites during solidification, resulting in grain refinement, improved mechanical properties, and reduced susceptibility to solidification cracking. Higher tensile strength and stiffness were demonstrated in Al 6061 based MMCs with 2 vol% reinforcement. These combined research efforts demonstrate the potential of PCTW technology for customizing of alloy composition and achieving microstructure control. A manufacturing process for titanium- and aluminum-based PCTWs was developed based on iterative design, theoretical calculations, and empirical data. PCTW technology successfully mitigated Al loss from Ti-6-4 and Mg loss in Al 6061, resulted in grain refinement in Ti-6-4, and demonstrated fabrication of Al based MMCs

Marcia S Domack

Lateral and Vertical Heterogeneity of Thorium in the Procellarum KREEP Terrane: As Reflected in the Ejecta Deposits of Post-Imbrium Craters

The Procellarum KREEP Terrane displays the highest concentrations of Th on the Moon. However, locations of elevated Th in this region appear to be random. As observed in the 5 deg per pixel equal-area Th data, and made more evident in the preliminary 2 deg data, Th is enhanced around the craters Aristillus, Aristarchus, Kepler, Mairan, the Apennine Bench formation, and the Fra Mauro region, while noticeably and unexpectedly lower in other locations (e.g., Archimedes, Copernicus, Eratosthenes, and Plato). We have examined the composition of the materials present in these regions with the goal of understanding the patchy nature to the distribution of Th and ultimately to decipher the geologic processes that have concentrated the Th. At present time, the published resolution of the Lunar Prospector Th gamma-ray data is low (5 deg per pixel), but this will soon be superceded by significantly higher-resolution data (2 deg per pixel). Even at this improved resolution, however, it is difficult to resolve the units that are the major source of Th. In an attempt to circumvent this problem, we employ the higher-resolution Clementine multispectral data for those regions mentioned above. We use the UV-VIS-derived compositional information and the spectral properties of craters, and their ejecta as drill holes through the mare-basalt surface to investigate the thickness and composition of underlying material. With this information we attempt to piece together the stratigraphy and geologic history of the Imbrium-Procellanim region. We processed the five-band multispectral data from the Clementine Mission (415, 750,900,950, and 1000nm) using ISIS software and calibration parameters developed by the USGS, Flagstaff, Arizona. Final image mosaics are in equal-area sinusoidal projection, and have a resolution of 250 m/pixel. Using the method of we produced maps of FeO and Ti02 composition. Here we examine the Th, FeO, and Ti02 composition and spectral properties of the craters discussed above and their ejecta, with the goal of describing the materials they excavate. One interpretation for the origin of the high-Th material is that subsurface KREEPy materials have been excavated by impact craters. The material excavated may be either volcanic KREEP (e,g., Apennine Bench Formation), KREEPy impact-melt breccia formed by the Imbrium impact (e.g., Fra Mauro Formation), or other KREEP-rich crustal material. Determining which type of material is responsible for the elevated Th and its extent is important to understanding the premare and possibly the prebasin stratigraphy of the Imbrium-Procellarum Region. Merging the 5 deg. Th data with the shaded relief map, we observe that the highest Th concentrations are not related to pre-Imbrium upper crustal materials. The Apennines, Alpes, and Caucasus Mountains represent the pre-Imbrian highlands material and do not express concentrations of Th, FeO, and TiO2 as high as the most Th-fich materials exposed within the Procellarum KREEP Terrane. We observe that, in general, these massifs contain 10-14 wt% FeO and 4-7 ppm Th. Determining whether the Th signal is from KREEP basalts or KREEPy impact-melt breccias cannot be done with the Clementine data because the two rock types are compositionally and mineralogically too similar (e.g., the Th-rich, mafic impact-melt breccias in the Apollo sample collection are dominated by a KREEP-basalt like component. Mapping-the distribution and sizes of craters and whether they display elevated Th concentrations or not, should reveal the depth and thickness of the KREEP-rich materials, and whether they are ubiquitous (i.e., impact-melt breccia) or more randomly distributed; this might be taken as an indicator of localized KREEP-basalt flows. Within the southeastern region of the Imbrium basin, there are two Th hot spots. The first is associated with the crater Aristillus, and the latter with the Apennine Bench Formation. Adjacent to these two hot spots are craters with a lower Th signature: Archimedes and Autolycus. We observe in the ejecta of Aristillus, a region of significantly lower FeO (10-14 wt%) relative to the surrounding mare basalt. The crater Autolycus, 50 km to the south, did not excavate similar low-FeO material. We suggest that the lower-FeO material in the ejecta of Aristillus corresponds to Th-rich material; the FeO content observed in Aristillus ejecta is comparable to that of KREEP basalt or mafic impact melt breccia (10-12 wt% FeO). We determine that this low FeO, Th-rich material is volcanic KREEP, as opposed to Imbrium impact melt, on the basis that the low-Fe material is exposed more prominently in ejecta in the northern portion of Aristillus. Our assumption is that if the layer underlying Aristillus was continuous, a more widespread and uniform low-Fe signature wouldbe observed in the ejecta deposit. Archimedes, 110 km southwest of Aristillus, impacted the northern portion of the Apennine Bench prior to the eruption of KREEP basalt. Archimedes rim material is not as enriched in Th as the Apennine Bench, and there are differences between the two in FeO concentration and in their continuum slope. Archimedes exhibits a much steeperor "redder" continuum slope than the Apennine Bench. This steepslope suggests the presence of glassy material. The glassy material is concentrated around an unnamed crater on the southern rim of Archimedes (4.5W, 28.2N) and along the northern rim of Archimedes. We suggest two possibilities, or a combination of the two, to explain the low-Th signal from Archimedes: (1) The Apennine Bench prior to KREEP basalt eruption was lower in Th (4-7 ppm, e.g., similar to the Apennine massifs) and KREEP basalts are absent in the rim of Archimedes; or (2) the glassy (possibly pyroclastic) material layering the rim of the Archimedes, dilutes any high-Th material present with low-Th material. (Additional information is contained in the original)

Gillis, J. J.