Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mixed precision”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Feasibility study: Monodisperse polymer particles containing laser-excitable dyes

The optical properties associated with small particles, which include aerosols, hydrosols and solid microspheres have an impact on several areas of science and engineering. Since the advent of high-speed computers and lasers, the interaction of light with matter in the form of small particles with a discontinuous optical boundary relative to the surroundings has been much better understood. Various nonlinear optical effects have been observed involving interaction of a laser beam with both solid microspheres and liquid microdroplets. These include observation of second and third harmonic generation, four wave mixing, optical visibility, two photon absorption, observation of stimulated emission and lasing, and Stimulated Raman Scattering. Many of these effects are observed with laser intensities which are orders of magnitude less than that required by threshold condition for interactions in macroscopic bulk medium. The primary reason for this is twofold. The front surface of the microsphere acts as a thick lens to enhance the internal intensity of the input laser radiation, and the spherical shape of the droplet acts as an optical cavity to provide feedback at specific wavelengths corresponding to the whispering gallery modes or the morphology dependent resonances (MDR's). The most interesting and significant recent finding in this field is undoubtedly the existence of resonance peaks in linear and nonlinear optical spectra. Such resonance peaks are only dependent on the particle morphology, which means the size, shape and refractive index of the particle. Because of the simultaneous presence of these resonances, they have been referred to by many names, including structural resonances, whispering modes or whispering gallery modes, creeping waves, circumferential waves, surfaces modes, and virtual modes. All of these names refer to the same phenomena, i.e. morphology dependent resonances (MDR's) which has already been described and predicted precisely by electromagnetic theory and Loentz-Mie theory since 1908. MDR's can become important when the particle size (radius a) approaches and exceeds the wavelength of the electromagnetic wave (lambda) and the refractive index of the particle is greater than that of the surrounding medium. Such resonances are easiest to observe from a single particle with high symmetry, such as a sphere, spheroid, or cylinder. MDR's correspond to solutions of the characteristic equations of the electromagnetic fields in the presence of a sphere.

Venkateswarlu, Putcha↗
NASA NTRS2009

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗
NASA NTRS2012

Twenty Years of Progress on Global Ocean Tides: The Impact of Satellite Altimetry

At the dawn of the era of high-precision altimetry, before the launch of TOPEX/Poseidon, ocean tides were properly viewed as a source of noise--tidal variations in ocean height would represent a very substantial fraction of what the altimeter measures, and would have to be accurately predicted and subtracted if altimetry were to achieve its potential for ocean and climate studies. But to the extent that the altimetry could be severely contaminated by tides, it also represented an unprecedented global-scale tidal data set. These new data, together with research stimulated by the need for accurate tidal corrections, led to a renaissance in tidal studies in the oceanographic community. In this paper we review contributions of altimetry to tidal science over the past 20 years, emphasizing recent progress. Mapping of tides has now been extended from the early focus on major constituents in the open ocean to include minor constituents, (e.g., long-period tides; non-linear tides in shelf waters, and in the open ocean), and into shallow and coastal waters. Global and spatially local estimates of tidal energy balance have been refined, and the role of internal tide conversion in dissipating barotropic tidal energy is now well established through modeling, altimetry, and in situ observations. However, energy budgets for internal tides, and the role of tidal dissipation in vertical ocean mixing remain controversial topics. Altimetry may contribute to resolving some of these important questions through improved mapping of low-mode internal tides. This area has advanced significantly in recent years, with several global maps now available, and progress on constraining temporally incoherent components. For the future, new applications of altimetry (e.g., in the coastal ocean, where barotropic tidal models remain inadequate), and new mission concepts (studies of the submesoscale with SWOT, which will require correction for internal tides) may bring us full circle, again pushing further development of tidal models as corrections.

Egbert, Gary↗
NASA NTRS1999

Global Geochemical Variation on the Lunar Surface: A Three-Element Approach

We present a method for displaying the relative abundances of three important elements (Th, Fe, and Ti) on the same map projection of the lunar surface. Using Th-, Fe-, and Ti-elemental abundances from orbital geochemical data and assigning each element a primary color, a false-color map of the lunar surface was created. This approach is similar to the ternary diagram approach presented by Davis and Spudis with some important differences, discussed later. For the present maps, Th abundances were measured by the Lunar Prospector (LP) Gamma-Ray Spectrometer(GRS).The new LPGRS low-altitude dataset was used in this analysis. Iron and Ti weight percentages were based on Clementine spectral reflectance data smoothed to the LP low altitude footprint. This method of presentation was designed to aid in the location and recognition of three principal lunar compositions: ferroan anorthosite (FAN), mare basalts (MB), and the Mg suite/ KREEP-rich rocks on the lunar surface, with special emphasis on the highlands and specific impact basins. In addition to the recognition of these endmember rock compositions, this method is an attempt to examine the relationship between elemental compositions that do not conform readily to previously accepted or observed endmember rocks in various specific regions of interest, including eastern highlands regions centered on 150 deg longitude, and a northern highlands Th-rich region observed. The LP low-altitude data has full width at half-maximum spatial resolution of about 40 km. The Clementine spectral reflectance datasets were adapted using an equal-area, gaussian smoothing routine to this footprint. In addition, these datasets, reported in weight percent of FeO and of Ti02, were adjusted to Fe and Ti weight percentages. Each dataset was then assigned one of the three primary colors: blue for Th, red for Fe, and green for Ti. For each element, the data range was normalized to represent the ratio of each point to the maximum in the dataset. (To view the color table, go to http://cass.jsc.nasa.gov/meetings/moon99/pdf/8033.pdf.) The full range of lunar longitudes is represented, but due to the lack of coverage of the Clementine data for latitudes > 70 deg and <-70 deg, the data for these regions is excluded. The differences between this approach and the ternary diagram approach of Davis and Spudis eliminate some of the uncertainty and ambiguity of the ternary diagram approach. Rather than using a ratio of Th to Ti normalized to CI chondritic ratios, and a ternary diagram with ternary apexes located at specific endmember compositional values, elemental compositions were used independently, eliminating the errors resulting from dividing numbers that can have high uncertainties, especially at low concentration. The three elements used in this method of presentation were chosen for several reasons. One reason for the inclusion of Th in this study is that it is an accurate indicator of KREEP. Iron and Ti concentrations are both low in highland regolith, causing any small fluctuations in Th to stand out very well. In addition, Fe and Ti are good compositional indicators of different mare basalts. Mixed with red for Fe, the green for Ti produces a yellow signal in high-Ti basalts. While universally high in Fe relative to the surrounding highlands, mare basalts have a diverse range of Ti values, making Ti concentration a valuable asset to the classification and identification of different basalt types. Finally, an important constraint in element selection is the availability of the global data, both from LP and Clementine results. Data for Th, Fe, and Ti are among the highest quality of existing lunar remote-sensing data. In addition, LP data for Fe and Ti will become available, enabling these data to be incorporated into the analysis. Using upper-limit values for end member rock compositions calculated from Korotev et al., attempts were made to locate the different endmember compositions of terranes on this diagram. Most strikingly, ferroan anorthosite (Th < and = 0.37 micro g/g; Fe (wt%)< and =2.29; Ti (wt%) < and = 0.22), which should appear as an almost black, reddish color, does not appear on the diagram at any noticeable frequency. Based on this analysis, the suggestion of extensive FAN regions on the lunar surface is not strong, especially at the presently accepted values for Fe and Th. However, to make sure this effect is not due to systematic errors, a thorough investigation of the precision, accuracy, and uncertainties of the Fe, Ti, and Th abundances needs to be carried out, especially at low concentrations. A particular region of interest is an area of high Th concentrations relative to Fe and Ti content north and east of Humboldtianum Crater. First observed by Lawrence et al., this region does not coincide with any visible impact structure and comprises one of the closest approximations to pure blue (high Th, very low Ti and Fe) on the lunar surface. Such an elemental composition does not lend itself readily to classification, and presents something of an anomaly. More detailed analysis of this region is needed to understand its structure and origin. There seems to be a longitudinal asymmetry in the Th concentrations of the highlands regolith. High-Th, low-Ti, and Fe regions are located between 135 deg and 180 deg longitude and between -30 deg and +30 deg latitude. While the Th levels are not high enough to attract attention in a single elemental display, the variation in the abundance of Th relative to Fe and Ti abundances can be clearly seen. The composition that these data suggest is not well represented in the sample return suite. In addition, these regions were largely missed by the Apollo orbital ground tracks, which only covered the outer edge of the areas of interest. The LP orbital Th data represent the first information about the Th concentrations in these regions of the highlands. Additional information contained in original.

Thomsen, D. R.↗
NASA NTRS2001

Our Sun V: A Bright Young Sun Consistent with Helioseismology and Warm Temperatures on Ancient Earth and Mars

The relatively warm temperatures required on early Earth and Mars have been difficult to account for with warming from greenhouse gases. A slightly more massive young Sun would be brighter than predicted by the standard solar model, simultaneously resolving this problem for both Earth and Mars. We computed high-precision solar models with seven initial masses, from Mi = 1.01 to 1.07 solar mass - the latter being the maximum permitted if the early Earth is not to lose its water via a moist greenhouse effect. The relatively modest early mass loss that is required remains consistent with observational limits on mass loss from young stars and with estimates of the past solar wind obtained from lunar rocks. We considered three types of mass loss rates: (1) a reasonable choice of a simple exponential decline, (2) an extreme step-function case that gives the maximum effect consistent with observations, and (3) the radical case of a linear decline which is inconsistent with the solar wind mass loss estimates from lunar rocks. Our computations demonstrated that mass loss leaves a fingerprint oil the Sun's internal structure large enough to be detectable with helioseismic observations. All of our mass-losing solar models were consistent with the helioseismic observations; in fact, our preferred mass-losing cases were in marginally better agreement with the helioseismology than the standard solar model was, although this difference was smaller than the effects of other uncertainties in the input physics and in the solar composition. Mass loss has only a relatively minor effect on the predicted lithium depletion; the major portion of the solar lithium depletion must still be due to rotational mixing. Thus the modest mass loss cases considered here cannot be ruled out by observed lithium depletions. For the three mass loss types considered, the preferred initial masses were 1.07 solar mass for the exponential case and 1.04 solar mass for the step-function and linear cases; all of these provided high enough solar fluxes at Mars 3.8 Gyr ago to be consistent with the existence of liquid water. For a more massive early Sun, the planets would have had to be closer to the young Sun in order to end up in their present orbits; the orbital radii of the planets would vary inversely with the solar mass. Both of these effects contribute to the fact that the early solar flux at the planets would have been considerably higher than that of the standard solar model at that time. In fact, the 1.07 solar mass exponential case has a flux at birth 5% higher than the present solar flux, while the radical 1.04 solar mass linear case has a nearly constant flux over the first 3 Gyr only about 10% lower than at present. The early solar evolution would be in the opposite direction in the H-R diagram to that of the standard Sun.

Sackmann, I.-Juliana↗
NASA NTRS2003

A Physical Validation Program for the GPM Mission

The GPM mission is currently planned for start in the late 2007 - early 2008 time frame. Its main scientific goal is to help answer pressing scientific problems arising within the context of global and regional water cycling. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like core satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe -- continuously. The constellation's orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. A major requirement before the retrieved rainfall information generated by the GPM mission can be used effectively by prognostic models to improve weather forecasts, hydrometeorological forecasts, and climate model reanalysis simulations is a capability to quantify the error characteristics of the retrievals. A solution for this problem has been upheld in past precipitation missions because of the lack of suitable error modeling systems incorporated into the validation programs and data distribution systems. An overview of how NASA intends to overcome this problem for the GPM mission using a physically-based error modeling approach within a multi-faceted validation program is described. The solution is to first identify specific user requirements and then determine the most stringent of these requirements that embodies all essential error characterization information needed by the entire user community. In the context of NASA s scientific agenda for the GPM mission, the most stringent user requirement is found within the data assimilation community. The fundamental theory of data assimilation vis-a-vis ingesting satellite precipitation information into the pre-forecast initializations is based on quantifying the conditional bias and precision errors of individual rain retrievals, and the space-time structure of the precision error (i.e., the spatial-temporal error covariance). By generating the hardware and software capability to produce this information in a near real-time fashion, and to couple the derived quantitative error properties to the actual retrieved rainrates, all key validation users can be satisfied. The talk will describe the essential components of the hardware and software systems needed to generate such near real-time error properties, as well as the various paradigm shifts needed within the validation community to produce a validation program relevant to the precipitation user community.

Smith, Eric A.↗
NASA NTRS2000

Analysis of the Radar Reflectivity of Aircraft Vortex Wakes

Radar has been proposed as a way to track wake vortices to reduce aircraft spacing and tests have revealed radar echoes from aircraft wakes in clear air. The results are always interpreted qualitatively using Tatarski's theory of weak scattering by isotropic atmospheric turbulence. The goal of the present work was to predict the value of the radar cross-section (RCS) using simpler models. This is accomplished in two steps. First, the refractive index is obtained. Since the structure of the aircraft wakes is different from atmospheric turbulence, three simple mechanisms specific to vortex wakes are considered: (1) Radial density gradient in a two-dimensional vortex, (2) three-dimensional fluctuations in the vortex cores, and (3) Adiabatic transport of the atmospheric fluid in a two-dimensional oval surrounding the pair of vortices. The index of refraction is obtained more precisely for the two-dimensional mechanisms than for the three-dimensional ones. In the second step, knowing the index of refraction, a scattering analysis is performed. Tatarski's weak scattering approximation is kept but the usual assumptions of a far-field and a uniform incident wave are dropped. Neither assumption is generally valid for a wake that is coherent across the radar beam. For analytical insight, a simpler approximation that invokes, in addition to weak scattering, the far-field and wide cylindrical beam assumptions, is also developed and compared with the more general analysis. The predicted RCS values for the oval surround the vortices (mechanism C) agree with the experiments of Bilson conducted over a wide range of frequencies. However, the predictions have a cut-off away from normal incidence which is not present in the measurements. Estimates suggest that this is due to turbulence in the baroclinic vorticity generated at the boundary of the oval. The reflectivity of a vortex itself (mechanism A) is comparable to that of the oval (mechanism C) but cuts-off at frequencies lower than those considered in all the experiments to date. The RCS of a vortex happens to peak at the frequency (about 49 MHz) where atmospheric radars (known as ST radars) operate and so the present prediction could be verified in the future. Finally , we suggest that hot engine exhaust could increase RCE by 40 db and reveal vortex circulation, provided its mixing with the surroundings is prevented in the laminarising flow of the vortices.

Shariff, Karim↗
NASA NTRS2003

Archetypal TRMM Radar Profiles Identified Through Cluster Analysis

It is widely held that identifiable 'convective regimes' exist in nature, although precise definitions of these are elusive. Examples include land / Ocean distinctions, break / monsoon beahvior, seasonal differences in the Amazon (SON vs DJF), etc. These regimes are often described by differences in the realized local convective spectra, and measured by various metrics of convective intensity, depth, areal coverage and rainfall amount. Objective regime identification may be valuable in several ways: regimes may serve as natural 'branch points' in satellite retrieval algorithms or data assimilation efforts; one example might be objective identification of regions that 'should' share a similar 2-R relationship. Similarly, objectively defined regimes may provide guidance on optimal siting of ground validation efforts. Objectively defined regimes could also serve as natural (rather than arbitrary geographic) domain 'controls' in studies of convective response to environmental forcing. Quantification of convective vertical structure has traditionally involved parametric study of prescribed quantities thought to be important to convective dynamics: maximum radar reflectivity, cloud top height, 30-35 dBZ echo top height, rain rate, etc. Individually, these parameters are somewhat deficient as their interpretation is often nonunique (the same metric value may signify different physics in different storm realizations). Individual metrics also fail to capture the coherence and interrelationships between vertical levels available in full 3-D radar datasets. An alternative approach is discovery of natural partitions of vertical structure in a globally representative dataset, or 'archetypal' reflectivity profiles. In this study, this is accomplished through cluster analysis of a very large sample (0[107) of TRMM-PR reflectivity columns. Once achieved, the rainconditional and unconditional 'mix' of archetypal profile types in a given location and/or season provides a description of the local convective spectrum which retains vertical structure information. A further cluster analysis of these 'mixes' can identify recurrent convective spectra. These are a first step towards objective identification of convective regimes, and towards answering the question: 'What are the most convectively similar locations in the world?'

Boccippio, Dennis J.↗
NASA NTRS2011

An Empirical State Error Covariance Matrix for Batch State Estimation

State estimation techniques serve effectively to provide mean state estimates. However, the state error covariance matrices provided as part of these techniques suffer from some degree of lack of confidence in their ability to adequately describe the uncertainty in the estimated states. A specific problem with the traditional form of state error covariance matrices is that they represent only a mapping of the assumed observation error characteristics into the state space. Any errors that arise from other sources (environment modeling, precision, etc.) are not directly represented in a traditional, theoretical state error covariance matrix. Consider that an actual observation contains only measurement error and that an estimated observation contains all other errors, known and unknown. It then follows that a measurement residual (the difference between expected and observed measurements) contains all errors for that measurement. Therefore, a direct and appropriate inclusion of the actual measurement residuals in the state error covariance matrix will result in an empirical state error covariance matrix. This empirical state error covariance matrix will fully account for the error in the state estimate. By way of a literal reinterpretation of the equations involved in the weighted least squares estimation algorithm, it is possible to arrive at an appropriate, and formally correct, empirical state error covariance matrix. The first specific step of the method is to use the average form of the weighted measurement residual variance performance index rather than its usual total weighted residual form. Next it is helpful to interpret the solution to the normal equations as the average of a collection of sample vectors drawn from a hypothetical parent population. From here, using a standard statistical analysis approach, it directly follows as to how to determine the standard empirical state error covariance matrix. This matrix will contain the total uncertainty in the state estimate, regardless as to the source of the uncertainty. Also, in its most straight forward form, the technique only requires supplemental calculations to be added to existing batch algorithms. The generation of this direct, empirical form of the state error covariance matrix is independent of the dimensionality of the observations. Mixed degrees of freedom for an observation set are allowed. As is the case with any simple, empirical sample variance problems, the presented approach offers an opportunity (at least in the case of weighted least squares) to investigate confidence interval estimates for the error covariance matrix elements. The diagonal or variance terms of the error covariance matrix have a particularly simple form to associate with either a multiple degree of freedom chi-square distribution (more approximate) or with a gamma distribution (less approximate). The off diagonal or covariance terms of the matrix are less clear in their statistical behavior. However, the off diagonal covariance matrix elements still lend themselves to standard confidence interval error analysis. The distributional forms associated with the off diagonal terms are more varied and, perhaps, more approximate than those associated with the diagonal terms. Using a simple weighted least squares sample problem, results obtained through use of the proposed technique are presented. The example consists of a simple, two observer, triangulation problem with range only measurements. Variations of this problem reflect an ideal case (perfect knowledge of the range errors) and a mismodeled case (incorrect knowledge of the range errors).

Frisbee, Joseph H., Jr.↗
NASA NTRS1997

The Complex Stratigraphy of the Highland Crust in the Serenitatis Region of the Moon Inferred from Mineral Fragment Chemistry

Large impact basins are natural drill holes into the Moon, and their ejecta carries unique information about the rock types and stratigraphy of the lunar crust. We have conducted an electron microprobe study of mineral fragments in the poikilitic melt breccias collected from the Taurus Mountains at the Apollo 17 landing site. These breccias are virtually unanimously agreed to be impact melt produced in the Serenitatis impact event. They contain lithic fragments and much more abundant mineral fragments of crustal origin. We have made precise microprobe analyses of minor element abundances in fragments of olivine, pyroxene, and plagioclase to provide new information on the possible source rocks and the crustal stratigraphy in the Serenitatis region. These data were also intended to elucidate the nature of the cryptic geochemical component in breccias such as these with low-K Fra Mauro basalt compositions. We chose the finest-grained (i.e., most rapidly quenched) breccias for study, to avoid reacted and partly assimilated fragments as much as possible. Most of the mineral fragments appear to have been derived from rocks that would fall into the pristine igneous Mg-suite as represented by lithic fragments in the Apollo collection, or reasonable extensions of it. Gabbroic rocks were more abundant in the target stratigraphy than is apparent from the Apollo sample collection. Some pyroxene and plagiociase, but probably not much olivine, could be derived from feldspathic granulites, which are metamorphosed polymict breccias. Some mineral fragments are from previously unknown rocks. These include highly magnesian olivines (up to Fo(sub 94)), possibly volcanic in origin, that exacerbate the difficulty in explaining highly magnesian rocks in the lunar crust. It appears that some part of the lunar interior has an mg*(= 100 x Mg/(Mg/Fe) atomic) greater than the conventional bulk Moon value of 80-84. Other volcanic rocks, including mare basalts, and rapidly- cooled impact melt rocks do not contribute significantly to the fragment population. Nor do ferroan anorthosites contribute more than a tiny part of even the plagiociase fragment population. A few mineral fragments that are consistent with the cryptic low-K Fra Mauro chemical component were found, and these appear to be from gabbroic sources. The mineral fragment populations cannot be mixed in their observed proportions to produce the whole rock composition, because the fragments are more refractory and deficient in Ti, P, and alkalis. A preferential contribution to the melt from a rock similar to sodic ferrogabbro can partly resolve the discrepancy. The population of mineral fragments requires a very diverse population of igenous rocks that are not all related to each other, demonstrating the existence of a complex crust built of numerous separate igneous plutons. Many of these plutons may have crystallized at shallow depths. The chemical composition of the melt breccias, in combination with the mineral fragment data and an understanding of the cratering process, suggests that the deepest crust sampled by the Serenitatis impace (not necessarily the deepest crust) was basaltic in composition, including KREEP and gabbroic rocks like sodic ferrogabbro, and lacking abundant olivine-rich material. These were overlain by Mg-suite rocks of varied types, including norites and troctolites that supplied most of the olivine mineral fragments. Granulities, which are metamorphosed and more feldspathic breccias, were abundant near the surface. Remote sensing indicates that the entire Serenitatis region lacks ferroan anorthosite, consistent with the results of our study.

Ryder, Graham↗
NASA NTRS2019

Bounding Methods for Heterogeneous Lidar-derived Navigational Geofences

Safe Unmanned Aerial Vehicle (UAV) operations near the ground require navigation methods that avoid fixed obstacles such as buildings, power lines and trees. Aerial lidar surveys of ground structures are available with the precision and accuracy to geolocate obstacles, but the high volume of raw survey data can exceed the compute power of onboard processors and the rendering ability of ground-based flight planning maps. Representing ground structures with bounding polyhedra instead of point clouds greatly reduces the data size and can enable effective obstacle avoidance, as long as the bounding geometry envelopes the structures with high spatial fidelity. This report describes in detail four methods to compute bounding geometries of ground obstacles from lidar point clouds. The four methods are: 1) 2.5D Maximum Elevation Box, 2) 2.5D Ground Map Extrusion, 3) 3D Bounding Cylinder, and 4) 3D Bounding Box. The methods are applied to five point cloud datasets from lidar surveys of UAV flight research sites in Georgia and Virginia with an average point spacing that ranges from 0.1m to 0.6m. The methods are assessed using survey areas with geometrically heterogeneous ground structures: buildings, vegetation, power lines, and sub-meter structures such as road signs and guy wires. The 2.5D Maximum Elevation Box method is useful for simple structures. The 2.5D Ground Map Extrusion method efficiently encloses vegetation, but requires handdrawn ground footprints. The 3D Bounding Cylinder method excels at enclosing linear structures such as power lines and fences. The 3D Bounding Box method excels at enclosing planar structures such as buildings. The methods are compared on the basis of data compression and boundary fidelity on selected areas. The 2.5D methods yield the highest data compression but the polyhedra produced by them enclose significant amounts of empty space. Boundary fidelity is superior for the 3D methods, though this fidelity comes at the cost of a roughly thirtyfold lower data compression ratio than the 2.5D Maximum Elevation Box method. A mix of these output geometries is proposed for autonomous UAV navigation with limited on-board computing. Both the accuracy and spatial detail of emerging satellite-based survey technology lower than that of aerial lidar scanning survey technology. Sub-meter structures and thin linear structures are not reliably mapped at present by satellite-based surveys.

Moore, Andrew J.↗
NASA NTRS2020

Using Two-Stream Theory to Capture Fluctuations of Satellite-Perceived TOA SW Radiances Reflected from Clouds Over Ocean

Shortwave (SW) fluxes estimated from broadband radiometry rely on empirically gathered and hemispherically resolved fields of outgoing top-of-atmosphere (TOA) radiances. This study aims to provide more accurate and precise fields of TOA SW radiances reflected from clouds over ocean by introducing a novel semiphysical model predicting radiances per narrow sun-observer geometry. This model was statistically trained using CERES-measured radiances paired with MODIS-retrieved cloud parameters as well as reanalysis-based geophysical parameters. By using radiative transfer approximations as a framework to ingest the above parameters, the new approach incorporates cloud-top effective radius and above-cloud water vapor in addition to traditionally used cloud optical depth, cloud fraction, cloud phase, and surface wind speed. A two-stream cloud albedo – serving to statistically incorporate cloud optical thickness and cloud-top effective radius – and Cox–Munk ocean reflectance were used to describe an albedo over each CERES footprint. Effective-radius-dependent asymmetry parameters were obtained empirically and separately for each viewing-illumination geometry. A simple equation of radiative transfer, with this albedo and attenuating above-cloud water vapor as inputs, was used in its log-linear form to allow for statistical optimization. We identified the two-stream functional form that minimized radiance residuals calculated against CERES observations and outperformed the state-of-the-art approach for most observer geometries outside the sun-glint and solar zenith angles between 20 and 70∘, reducing the median SD of radiance residuals per solar geometry by up to 13.2 % for liquid clouds, 1.9 % for ice clouds, and 35.8 % for footprints containing both cloud phases. Geometries affected by sun glint (constituting between 10 % and 1 % of the discretized upward hemisphere for solar zenith angles of 20 and 70∘, respectively), however, often showed weaker performance when handled with the new approach and had increased residuals by as much as 60 % compared to the state-of-the-art approach. Overall, uncertainties were reduced for liquid-phase and mixed-phase footprints by 5.76 % and 10.81 %, respectively, while uncertainties for ice-phase footprints increased by 0.34 %. Tested for a variety of scenes, we further demonstrated the plausibility of scene-wise predicted radiance fields. This new approach may prove useful when employed in angular distribution models and may result in improved flux estimates, in particular dealing with clouds characterized by small or large droplet/crystal sizes.

Two- stream theory↗
NASA NTRS2013

Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays

Single-slope analog-to-digital converters (ADCs) are particularly useful for onchip digitization in focal plane arrays (FPAs) because of their inherent monotonicity, relative simplicity, and efficiency for column-parallel applications, but they are comparatively slow. Squareroot encoding can allow the number of code values to be reduced without loss of signal-to-noise ratio (SNR) by keeping the quantization noise just below the signal shot noise. This encoding can be implemented directly by using a quadratic ramp. The reduction in the number of code values can substantially increase the quantization speed. However, in an FPA, the fixed pattern noise (FPN) limits the use of small quantization steps at low signal levels. If the zero-point is adjusted so that the lowest column is onscale, the other columns, including those at the center of the distribution, will be pushed up the ramp where the quantization noise is higher. Additionally, the finite frequency response of the ramp buffer amplifier and the comparator distort the shape of the ramp, so that the effective ramp value at the time the comparator trips differs from the intended value, resulting in errors. Allowing increased settling time decreases the quantization speed, while increasing the bandwidth increases the noise. The FPN problem is solved by breaking the ramp into two portions, with some fraction of the available code values allocated to a linear ramp and the remainder to a quadratic ramp. To avoid large transients, both the value and the slope of the linear and quadratic portions should be equal where they join. The span of the linear portion must cover the minimum offset, but not necessarily the maximum, since the fraction of the pixels above the upper limit will still be correctly quantized, albeit with increased quantization noise. The required linear span, maximum signal and ratio of quantization noise to shot noise at high signal, along with the continuity requirement, determines the number of code values that must be allocated to each portion. The distortion problem is solved by using a lookup table to convert captured code values back to signal levels. The values in this table will be similar to the intended ramp value, but with a correction for the finite bandwidth effects. Continuous-time comparators are used, and their bandwidth is set below the step rate, which smoothes the ramp and reduces the noise. No settling time is needed, as would be the case for clocked comparators, but the low bandwidth enhances the distortion of the non-linear portion. This is corrected by use of a return lookup table, which differs from the one used to generate the ramp. The return lookup table is obtained by calibrating against a stepped precision DC reference. This results in a residual non-linearity well below the quantization noise. This method can also compensate for differential non-linearity (DNL) in the DAC used to generate the ramp. The use of a ramp with a combination of linear and quadratic portions for a single-slope ADC is novel. The number of steps is minimized by keeping the step size just below the photon shot noise. This in turn maximizes the speed of the conversion. High resolution is maintained by keeping small quantization steps at low signals, and noise is minimized by allowing the lowest analog bandwidth, all without increasing the quantization noise. A calibrated return lookup table allows the system to maintain excellent linearity.

Wrigley, Chris J.↗
NASA NTRS2010

Airborne and Ground-Based Measurements Using a High-Performance Raman Lidar

A high-performance Raman lidar operating in the UV portion of the spectrum has been used to acquire, for the first time using a single lidar, simultaneous airborne profiles of the water vapor mixing ratio, aerosol backscatter, aerosol extinction, aerosol depolarization and research mode measurements of cloud liquid water, cloud droplet radius, and number density. The Raman Airborne Spectroscopic Lidar (RASL) system was installed in a Beechcraft King Air B200 aircraft and was flown over the mid-Atlantic United States during July August 2007 at altitudes ranging between 5 and 8 km. During these flights, despite suboptimal laser performance and subaperture use of the telescope, all RASL measurement expectations were met, except that of aerosol extinction. Following the Water Vapor Validation Experiment Satellite/Sondes (WAVES_2007) field campaign in the summer of 2007, RASL was installed in a mobile trailer for groundbased use during the Measurements of Humidity and Validation Experiment (MOHAVE-II) field campaign held during October 2007 at the Jet Propulsion Laboratory s Table Mountain Facility in southern California. This ground-based configuration of the lidar hardware is called Atmospheric Lidar for Validation, Interagency Collaboration and Education (ALVICE). During theMOHAVE-II field campaign, during which only nighttime measurements were made, ALVICE demonstrated significant sensitivity to lower-stratospheric water vapor. Numerical simulation and comparisons with a cryogenic frost-point hygrometer are used to demonstrate that a system with the performance characteristics of RASL ALVICE should indeed be able to quantify water vapor well into the lower stratosphere with extended averaging from an elevated location like Table Mountain. The same design considerations that optimize Raman lidar for airborne use on a small research aircraft are, therefore, shown to yield significant dividends in the quantification of lower-stratospheric water vapor. The MOHAVE-II measurements, along with numerical simulation, were used to determine that the likely reason for the suboptimal airborne aerosol extinction performance during theWAVES_2007 campaign was amisaligned interference filter. With full laser power and a properly tuned interference filter,RASL is shown to be capable ofmeasuring themain water vapor and aerosol parameters with temporal resolutions of between 2 and 45 s and spatial resolutions ranging from 30 to 330 m from a flight altitude of 8 km with precision of generally less than 10%, providing performance that is competitive with some airborne Differential Absorption Lidar (DIAL) water vapor and High Spectral Resolution Lidar (HSRL) aerosol instruments. The use of diode-pumped laser technology would improve the performance of an airborne Raman lidar and permit additional instrumentation to be carried on board a small research aircraft. The combined airborne and ground-based measurements presented here demonstrate a level of versatility in Raman lidar that may be impossible to duplicate with any other single lidar technique.

Whiteman, David N.↗
NASA NTRS2013

The W-W02 Oxygen Fugacity Buffer at High Pressures and Temperatures: Implications for f02 Buffering and Metal-silicate Partitioning

Oxygen fugacity (fO2) controls multivalent phase equilibria and partitioning of redox-sensitive elements, and it is important to understand this thermodynamic parameter in experimental and natural systems. The coexistence of a metal and its oxide at equilibrium constitutes an oxygen buffer which can be used to control or calculate fO2 in high pressure experiments. Application of 1-bar buffers to high pressure conditions can lead to inaccuracies in fO2 calculations because of unconstrained pressure dependencies. Extending fO2 buffers to pressures and temperatures corresponding to the Earth's deep interior requires precise determinations of the difference in volume (Delta) V) between the buffer phases. Synchrotron x-ray diffraction data were obtained using diamond anvil cells (DAC) and a multi anvil press (MAP) to measure unit cell volumes of W and WO2 at pressures and temperatures up to 70 GPa and 2300 K. These data were fitted to Birch-Murnaghan 3rd-order thermal equations of state using a thermal pressure approach; parameters for W are KT = 306 GPa, KT' = 4.06, and αKT = 0.00417 GPa K-1. Two structural phase transitions were observed for WO2 at 4 and 32 GPa with structures in P21/c, Pnma and C2/c space groups. Equations of state were fitted for these phases over their respective pressure ranges yielding the parameters KT = 190, 213, 300 GPa, KT' = 4.24, 5.17, 4 (fixed), and αKT = 0.00506, 0.00419, 0.00467 GPa K-1 for the P21/c, Pnma and C2/c phases, respectively. The W-WO2 buffer (WWO) was extended to high pressure by inverting the W and WO2 equations of state to obtain phase volumes at discrete pressures (1-bar to 100 GPa, 1 GPa increments) along isotherms (300 to 3000K, 100 K increments). The slope of the absolute fO2 of the WWO buffer is positive with increasing temperature up to approximately 70 GPa and is negative above this pressure. The slope is positive along isotherms from 1000 to 3000K with increasing pressure up to at least 100 GPa. The WWO buffer is at a higher fO2 than the IW buffer at pressures lower than 40 GPa, and the magnitude of this difference decreases at higher pressures. This qualitatively indicates an increasingly lithophile character for W at higher pressures. The WWO buffer was quantitatively applied to W metal-silicate partitioning by using the WWO-IW buffer difference in combination with literature data on W metal-silicate partitioning to model the exchange coefficient (KD) for the Fe-W exchange reaction. This approach captures the pressure dependence of W metal-silicate partitioning using the WWO-IW buffer difference and models the activities of the components in the silicate and metallic phases using an expression of the Gibbs excess energy of mixing. Calculation of KD along a peridotite liquidus predicts a decrease in W siderophility at higher pressures that supports the qualitative behavior predicted by the WWO-IW buffer difference, and agrees with findings of others. Comparing the competing effects of temperature and pressure on W metal-silicate partitioning, our results indicate that pressure exerts a greater effect.

Shofner, G. A.↗
NASA NTRS1995

Concept Designed and Developed for Distortion- Tolerant, High-Stability Engine Control

Engine Control Future aircraft turbine engines, both commercial and military, must be able to successfully accommodate expected increased levels of steady-state and dynamic engine-face distortion. Advanced tactical aircraft are likely to use thrust vectoring to enhance their maneuverability. As a result, the engines will see more extreme aircraft angles-of-attack and sideslip levels than are currently encountered with present-day aircraft. Also, the mixed-compression inlets needed for the High Speed Civil Transport will likely encounter disturbances similar to those seen by tactical aircraft, in addition to planar pulse, inlet buzz, and high distortion levels at low flight speed and off-design operation. The current approach of incorporating a sufficient component design stall margin to tolerate these increased levels of distortion would significantly reduce performance. The objective of the High Stability Engine Control (HISTEC) program is to design, develop, and flight demonstrate an advanced, high-stability, integrated engine-control system that uses measurement-based, real-time estimates of distortion to enhance engine stability. The resulting distortion-tolerant control reduces the required design stall margin, with a corresponding increase in performance and decrease in fuel burn. The HISTEC concept has been designed and developed, and the software implementing the concept has successfully accommodated time-varying distortion. The NASA Lewis Research Center is currently overseeing the development and validation of the hardware and software necessary to flight test the HISTEC concept. HISTEC is a contracted effort with Pratt & Whitney of West Palm Beach, Florida. The HISTEC approach includes two major systems: A Distortion Estimation System (DES) and Stability Management Control (SMC). DES is an aircraft-mounted, high-speed processor that estimates the amount and type of distortion present and its effect on the engine. It uses high-response pressure measurements at the engine face to calculate indicators of the type and extent of distortion in real time. From these indicators, DES determines the effects of distortion on the propulsion systems and the corresponding engine match point necessary to accommodate it. DES output consists of fan and compressor pressure ratio trim commands that are passed to the SMC. In addition, DES uses maneuver information, consisting of angle-of-attack and sideslip from the flight control, to anticipate high inlet distortion conditions. The SMC, which is contained in the engine-mounted, Improved Digital Electronic Engine Control (IDEEC), includes advanced control laws to directly control the fan and compressor transient operating line (pressure ratio). These advanced control laws, with a multivariable design, have the potential for higher bandwidth and the resulting more precise control of engine match. The ability to measure and assess the distortion effects in real time coupled with a high-response controller improves engine stability at high levels of distortion. The software algorithms implementing DES have been designed, developed, and demonstrated, and integration testing of the DES and SMC software has been completed. The results show that the HISTEC system will be able to sense inlet distortion, determine the effect on engine stability, and accommodate distortion by maintaining an adequate margin for engine surge. The Pratt &Whitney Comprehensive Engine Diagnostic Unit was chosen as the DES processor. An instrumented inlet case for sensing distortion was designed and fabricated. HISTEC is scheduled for flight test on the ACTIVE F-15 aircraft at the NASA Dryden Flight Research Center in Edwards, California, in late 1996.

Source record↗
NASA NTRS2002

Great (Flame) Balls of Fire! Structure of Flame Balls at Low Lewis-number-2 (SOFBALL-2)

Everyone knows that an automobile engine wastes fuel and energy when it runs with a fuel-rich mixture. 'Lean' burning, mixing in more air and less fuel, is better for the environment. But lean mixtures also lead to engine misfiring and rough operation. No one knows the ultimate limits for lean operation, for 'weak' combustion that is friendly to the environment while still moving us around. This is where the accidental verification of a decades-old prediction may have strong implications for designing and running low-emissions engines in the 21st century. In 1944, Soviet physicist Yakov Zeldovich predicted that stationary, spherical flames are possible under limited conditions in lean fuel-air mixtures. Dr. Paul Ronney of the University of Southern California accidentally discovered such 'flame balls' in experiments with lean hydrogen-air mixtures in 1984 during drop-tower experiments that provided just 2.2 seconds of near weightlessness. Experiments aboard NASA's low-g aircraft confirmed the results, but a thorough investigation was hampered by the aircraft's bumpy ride. And stable flame balls can only exist in microgravity. The potential for investigating combustion at the limits of flammability, and the implications for spacecraft fire safety, led to the Structure of Flame Balls at Low Lewis-number (SOFBALL) experiment flown twice aboard the Space Shuttle on the Microgravity Sciences Laboratory-1 (MSL-1) in 1997. Success there led to the planned reflight on STS-107. Flame balls are the weakest fires yet produced in space or on Earth. Typically each flame ball produced only 1 watt of thermal power. By comparison, a birthday candle produces 50 watts. The Lewis-number measures the rate of diffusion of fuel into the flame ball relative to the rate of diffusion of heat away from the flame ball. Lewis-number mixtures conduct heat poorly. Hydrogen and methane are the only fuels that provide low enough Lewis-numbers to produce stable flame balls, and even then only for very weak, barely flammable mixtures. Nevertheless, under these conditions flame balls give scientists the opportunity to test models in one of the simplest combustion experiments possible. SOFBALL-2 science objectives include: Improving our understanding of the flame ball phenomenon; Determining the conditions under which flame balls exist; Testing predictions of flame ball lifetimes; Acquiring more precise data for critical model comparison.

Ronney, Paul↗
NASA NTRS2013

Molecular-Based Optical Measurement Techniques for Transition and Turbulence in High-Speed Flow

High-speed laminar-to-turbulent transition and turbulence affect the control of flight vehicles, the heat transfer rate to a flight vehicle's surface, the material selected to protect such vehicles from high heating loads, the ultimate weight of a flight vehicle due to the presence of thermal protection systems, the efficiency of fuel-air mixing processes in high-speed combustion applications, etc. Gaining a fundamental understanding of the physical mechanisms involved in the transition process will lead to the development of predictive capabilities that can identify transition location and its impact on parameters like surface heating. Currently, there is no general theory that can completely describe the transition-to-turbulence process. However, transition research has led to the identification of the predominant pathways by which this process occurs. For a truly physics-based model of transition to be developed, the individual stages in the paths leading to the onset of fully turbulent flow must be well understood. This requires that each pathway be computationally modeled and experimentally characterized and validated. This may also lead to the discovery of new physical pathways. This document is intended to describe molecular based measurement techniques that have been developed, addressing the needs of the high-speed transition-to-turbulence and high-speed turbulence research fields. In particular, we focus on techniques that have either been used to study high speed transition and turbulence or techniques that show promise for studying these flows. This review is not exhaustive. In addition to the probe-based techniques described in the previous paragraph, several other classes of measurement techniques that are, or could be, used to study high speed transition and turbulence are excluded from this manuscript. For example, surface measurement techniques such as pressure and temperature paint, phosphor thermography, skin friction measurements and photogrammetry (for model attitude and deformation measurement) are excluded to limit the scope of this report. Other physical probes such as heat flux gauges, total temperature probes are also excluded. We further exclude measurement techniques that require particle seeding though particle based methods may still be useful in many high speed flow applications. This manuscript details some of the more widely used molecular-based measurement techniques for studying transition and turbulence: laser-induced fluorescence (LIF), Rayleigh and Raman Scattering and coherent anti-Stokes Raman scattering (CARS). These techniques are emphasized, in part, because of the prior experience of the authors. Additional molecular based techniques are described, albeit in less detail. Where possible, an effort is made to compare the relative advantages and disadvantages of the various measurement techniques, although these comparisons can be subjective views of the authors. Finally, the manuscript concludes by evaluating the different measurement techniques in view of the precision requirements described in this chapter. Additional requirements and considerations are discussed to assist with choosing an optical measurement technique for a given application.

Bathel, Brett F.↗