Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function mapping”

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 631 records · Page 35

Unoccupied electronic resonances of Sc adsorbed on W(001) by k-resolved inverse photoemission

Scandium adsorbed on the (001) face of tungsten has been studied using Auger-electron spectroscopy, low-energy electron diffraction, k-resolved inverse photoelectron spectroscopy, work-function measurements, and relativistic-electronic-structure calculations. We find that the work function of W(001) does not decrease monotonically as a function of scandium adsorption but reaches a minimum value of Phi = 3.25 eV at 50% coverage. For a complete monolayer, an order 1 X 1 Sc overlayer is formed. Inverse photoemission in the isochromat mode was used to map the unoccupied energy levels of this overlayer along the bar-Gamma(bar-Mu) symmetry direction of the surface Brillouin zone. Local-density-functional calculations using the muffin-tin orbital method were performed for a 1 X 1 Sc overlayer on W and are compared with the experimental two-dimensional band structure.

Lamouri, A.↗

Continental Spatio-Temporal Data Analysis with Linear Spectral Mixture Model Using FOSS

This work demonstrates the development and implementation of a Fully Constrained Least Squares (FCLS) unmixing model developed in C++ programming language with OpenCV package and boost C++ libraries in the NASA Earth Exchange (NEX). Visualization of the results is supported by GRASS GIS and statistical analysis is carried in R in a Linux system environment. FCLS was first tested on computer simulated data with Gaussian noise of various signal-to-noise ratio, and Landsat data of an agricultural scenario and an urban environment using a set of global end members of substrate (soils, sediments, rocks, and non-photosynthetic vegetation), vegetation that includes green photosynthetic plants and dark objects which encompasses absorptive substrate materials, clear water, deep shadows, etc. For the agricultural scenario, a spectrally diverse collection of 11 scenes of Level 1 terrain corrected, cloud free Landsat-5 TM data of Fresno, California, USA were unmixed and the results were validated with the corresponding ground data. To study an urbanized landscape, a clear sky Landsat-5 TM data were unmixed and validated with coincident World View-2 abundance maps (of 2 m spatial resolution) for an area of San Francisco, California, USA. The results were evaluated using descriptive statistics, correlation coefficient, RMSE, probability of success, boxplot and bivariate distribution function. Finally, FCLS was used for sub-pixel land cover analysis of the monthly WELD (Wen-enabled Landsat data) repository from 2008 to 2011 of North America. The abundance maps in conjunction with DMSP-OLS nighttime lights data were used to extract the urban land cover features and analyze their spatial-temporal growth.

Landsat Satellites↗

Geophysical models of Western Aphrodite-Niobe region: Venus

The new topography and gravitational field data for Venus expressed in spherical harmonics of degree and order up to 50 allow us to analyze the crust-mantle boundary relief and stress state of the Venusian lithosphere. In these models, we consider models in which convection is confined beneath a thick, buoyant lithosphere. We divide the convection regime into an upper mantle and lower mantle component. The lateral scales are smaller than on Earth. In these models, relative to Earth, convection is reflected in higher order terms of the gravitational field. On Venus geoid height and topography are highly correlated, although the topography appears to be largely compensated. We hypothesize that Venus topography for those wavelengths that correlate well with the geoid is partly compensated at the crust-mantle boundary, while for the others compensation may be distributed over the whole mantle. In turn the strong sensitivity of the stresses to parameters of the models of the external layers of Venus together with geological mapping allows us to begin investigations of the tectonics and geodynamics of the planet. For stress calculations we use a new technique of space- and time-dependent Green's response functions using Venus models with rheologically stratified lithosphere and mantle and a ductile lower crust. In the basic model of Venus the mean crust is 50-70 km thick, the density contrast across the crust-mantle boundary is in the range from 0.3 to 0.4 g/cm(exp -3). The thickness of a weak mantle zone may be from 350 to 1000 km. Strong sensitivity of calculated stress to various parameters of the layered model of Venus together with geological mapping and analysis of surface tectonic patterns allow us to investigate the tectonics and geodynamics of the planet. The results are presented in the form of maps of compression-extension and maximum shear stresses in the lithosphere and maps of crust-mantle boundary relief, which can be presented as a function of time. We have modeled the region of Western Aphrodite and the Niobe plains to get reasonable depths of compensation. Crust mantle boundary relief is calculated for Western Aphrodite-Niobe relative to a mean crustal thickness of 50 km. The calculations include the consequences of simple crust models and more complicated models with a weak, ductile lower crust, a strong upper mantle and a weak lower mantle layer.

Marchenkov, K. I.↗

MicroShell Minimalist Shell for Xilinx Microprocessors

MicroShell is a lightweight shell environment for engineers and software developers working with embedded microprocessors in Xilinx FPGAs. (MicroShell has also been successfully ported to run on ARM Cortex-M1 microprocessors in Actel ProASIC3 FPGAs, but without project-integration support.) Micro Shell decreases the time spent performing initial tests of field-programmable gate array (FPGA) designs, simplifies running customizable one-time-only experiments, and provides a familiar-feeling command-line interface. The program comes with a collection of useful functions and enables the designer to add an unlimited number of custom commands, which are callable from the command-line. The commands are parameterizable (using the C-based command-line parameter idiom), so the designer can use one function to exercise hardware with different values. Also, since many hardware peripherals instantiated in FPGAs have reasonably simple register-mapped I/O interfaces, the engineer can edit and view hardware parameter settings at any time without stopping the processor. MicroShell comes with a set of support scripts that interface seamlessly with Xilinx's EDK tool. Adding an instance of MicroShell to a project is as simple as marking a check box in a library configuration dialog box and specifying a software project directory. The support scripts then examine the hardware design, build design-specific functions, conditionally include processor-specific functions, and complete the compilation process. For code-size constrained designs, most of the stock functionality can be excluded from the compiled library. When all of the configurable options are removed from the binary, MicroShell has an unoptimized memory footprint of about 4.8 kB and a size-optimized footprint of about 2.3 kB. Since MicroShell allows unfettered access to all processor-accessible memory locations, it is possible to perform live patching on a running system. This can be useful, for instance, if a bug is discovered in a routine but the system cannot be rebooted: Shell allows a skilled operator to directly edit the binary executable in memory. With some forethought, MicroShell code can be located in a different memory location from custom code, permitting the custom functionality to be overwritten at any time without stopping the controlling shell.

Werne, Thomas A.↗

Compact Augmented Spark Igniters for Liquid Rocket Engines

The objective of this dissertation is to characterize augmented spark impinging (torch) igniters. Torch igniters inject propellants into a small prechamber where a spark igniter imparts the activation energy to initiate combustion. The flame stabilizes within the prechamber and exits into the main combustion chamber where it ignites the main flow directly or by mixing with a coaxial flow to create a secondary flame. The underlying phenomena that promote reliable ignition are multifaceted and required a broad scope of investigations. Initial efforts focused on spark discharges across an annular geometry. Gaseous propellent projects the electrical arc toward the prechamber, where the arc and the thermal energy is imparted to the surrounding fluid then meets with a combustible mixture. Discharges are parametrically examined against pressure, spark-gap size, and exciter types to emulate transient engine startup. The energy imparted to the flow, approximate velocity, and distance traveled by the exhaust plume are determined. Spark quenching (no energy deposition) and its effects on perceived versus actual measurements are shown. Results outline the prechamber target zone where a combustible mixture is necessary for ignition. Next, three-dimensional, time-accurate, and nonreacting computational fluid dynamics simulations assess effects of geometric and mass flow differences on the size, location, and composition of a combustible mixture. Fluid dynamic phenomena are elucidated. Correlations between nondimensional variables are found. Results are coupled with a conservation-law analysis to produce novel objective functions. Lastly, a full-scale, modular torch igniter is tested. Injector sizes, mixture ratios, and injector momentum ratios are systematically varied to map ignition probability and ignition delay. Results show the importance of local mixture ratios, provide injector momentum ratios guidelines, and demonstrate reliable ignition of core mixture ratios typically outside flammability limits. The deliverable of this project is not a torch igniter or derivative thereof. Instead, this dissertation produces novel objective-functions for igniter design. These equations are functions of variables that may be constrained by engine-level requirements, physics, and experimental results from this work. They provide a foundation for future igniter design and will minimize iterative processes from the design cycle for next-generation ignition hardware.

Torch Igniter↗

Iterative-Transform Phase Retrieval Using Adaptive Diversity

A phase-diverse iterative-transform phase-retrieval algorithm enables high spatial-frequency, high-dynamic-range, image-based wavefront sensing. [The terms phase-diverse, phase retrieval, image-based, and wavefront sensing are defined in the first of the two immediately preceding articles, Broadband Phase Retrieval for Image-Based Wavefront Sensing (GSC-14899-1).] As described below, no prior phase-retrieval algorithm has offered both high dynamic range and the capability to recover high spatial-frequency components. Each of the previously developed image-based phase-retrieval techniques can be classified into one of two categories: iterative transform or parametric. Among the modifications of the original iterative-transform approach has been the introduction of a defocus diversity function (also defined in the cited companion article). Modifications of the original parametric approach have included minimizing alternative objective functions as well as implementing a variety of nonlinear optimization methods. The iterative-transform approach offers the advantage of ability to recover low, middle, and high spatial frequencies, but has disadvantage of having a limited dynamic range to one wavelength or less. In contrast, parametric phase retrieval offers the advantage of high dynamic range, but is poorly suited for recovering higher spatial frequency aberrations. The present phase-diverse iterative transform phase-retrieval algorithm offers both the high-spatial-frequency capability of the iterative-transform approach and the high dynamic range of parametric phase-recovery techniques. In implementation, this is a focus-diverse iterative-transform phaseretrieval algorithm that incorporates an adaptive diversity function, which makes it possible to avoid phase unwrapping while preserving high-spatial-frequency recovery. The algorithm includes an inner and an outer loop (see figure). An initial estimate of phase is used to start the algorithm on the inner loop, wherein multiple intensity images are processed, each using a different defocus value. The processing is done by an iterative-transform method, yielding individual phase estimates corresponding to each image of the defocus-diversity data set. These individual phase estimates are combined in a weighted average to form a new phase estimate, which serves as the initial phase estimate for either the next iteration of the iterative-transform method or, if the maximum number of iterations has been reached, for the next several steps, which constitute the outerloop portion of the algorithm. The details of the next several steps must be omitted here for the sake of brevity. The overall effect of these steps is to adaptively update the diversity defocus values according to recovery of global defocus in the phase estimate. Aberration recovery varies with differing amounts as the amount of diversity defocus is updated in each image; thus, feedback is incorporated into the recovery process. This process is iterated until the global defocus error is driven to zero during the recovery process. The amplitude of aberration may far exceed one wavelength after completion of the inner-loop portion of the algorithm, and the classical iterative transform method does not, by itself, enable recovery of multi-wavelength aberrations. Hence, in the absence of a means of off-loading the multi-wavelength portion of the aberration, the algorithm would produce a wrapped phase map. However, a special aberration-fitting procedure can be applied to the wrapped phase data to transfer at least some portion of the multi-wavelength aberration to the diversity function, wherein the data are treated as known phase values. In this way, a multiwavelength aberration can be recovered incrementally by successively applying the aberration-fitting procedure to intermediate wrapped phase maps. During recovery, as more of the aberration is transferred to the diversity function following successive iterations around the ter loop, the estimated phase ceases to wrap in places where the aberration values become incorporated as part of the diversity function. As a result, as the aberration content is transferred to the diversity function, the phase estimate resembles that of a reference flat.

Dean, Bruce H.↗

Passive microwave mapping of ice thickness

Basic calculations are presented for evaluating the feasibility of a scanning microwave radiometer system for mapping the thickness of lake ice. An analytical model for the apparent brightness temperature as a function of ice thickness has been developed, and elaborated to include such variables as galactic and atmospheric noise, aspect angle, polarization, temperature gradient in the ice, the presence of transition layers such as snow, slush, and water, increased loss due to air inclusions in the ice layer, and the presence of multiple ice thicknesses within the antenna footprint. It was found that brightness temperature measurements at six or seven frequencies in the range of 0.4 to 0.7 GHz were required to obtain unambiquous thickness estimates. A number of data processing methods were examined. The effects of antenna beamwidth, scanning rate, receiver bandwidth, noise figure, and integration time were studied.

Apinis, J. J.↗

Non-Gaussian statistical models of surface wave fields for remote sensing applications

Based on the complete Stokes wave model with the bias term and using a simple mapping approach and an iteration solution method, we established a formula for the joint probability density function of the surface slope elevation of a nonlinear random wave field. The formula requires three parameters to define the whole density function: the rms surface elevation and slope values and the significant slope. This model represents the dynamics of the wave in a more direct way than the Gram-Charlier approximation. Based on this new statistical model and laboratory experiments, formula and numerical values of EM bias and dynamics bias are derived. The results indicate that various biases should be considered seriously if accuracy of the altimeter measurement is required in centimeter range.

Huang, N. E.↗

A parametric study of various synthetic aperture telescope configurations for coherent imaging applications

The comparative advantages of synthetic aperture telescopes (SATs) of segmented primary mirror and common secondary mirror type, on the one hand, and on the other those employing an array of independent telescopes, are discussed. The diffraction-limited optical performance of both redundant and nonredundant subaperture configurations are compared in terms of point spread function characteristics and encircled energy plots. Coherent imaging with afocal telescope SATs involves a pupil-mapping operation followed by a Fourier transform one. A quantitative analysis of the off-axis optical performance degradation due to pupil-mapping errors is presented, together with the field-dependent effects of residual design aberrations of independent telescopes.

Harvey, James E.↗

Disparity coding - An approach for stereo reconstruction

As the possibility of stereo-controlled robots becomes a reality, the need to transmit the stereo pair of images to a ground station or space station for man-in-the-loop supervision will be a necessity. The complexity of transmitting stereo images by coding the preprocessed disparity is presently discussed. The approach demonstrates the quantization, modulation, and reconstruction of the stereo images. Results indicate the accuracy of reconstruction in terms of mean-square-error criterion as a function of the signal-to-noise ratio. Key research issues of interpolation from sparse disparity maps and reconstruction of the stereo pairs in the presence of spatial noise are presented. It is concluded that stereo reconstruction is possible, and the noise constraints are given.

Griswold, N. C.↗

Head-Strictness is Not a Monotonic Abstract Property

A property P of a language is said to be definable by abstract interpretation if there is a monotonic map abs from the domain of standard semantics to an abstract domain A of finite height, and a partition of the abstract domain into two parts A(sub p) and A(sub non p), such that any value has property P if and only if abs maps it to an element of A(sub p). Head-strictness is a property of functions over lists which asserts, roughly speaking, that whenever the function looks at the tail of a list, it looks at the head of the tail. We prove that head-strictness is not definable by abstract interpretation. We then present a non-monotonic abstract interpretation for head-strictness and prove its safety.

Kamin, Samuel↗

Precision Interval Estimation of the Response Surface by Means of an Integrated Algorithm of Neural Network and Linear Regression

The integration of Radial Basis Function Networks and Back Propagation Neural Networks with the Multiple Linear Regression has been accomplished to map nonlinear response surfaces over a wide range of independent variables in the process of the Modem Design of Experiments. The integrated method is capable to estimate the precision intervals including confidence and predicted intervals. The power of the innovative method has been demonstrated by applying to a set of wind tunnel test data in construction of response surface and estimation of precision interval.

Lo, Ching F.↗

Multispectral Microimager for Astrobiology

A primary goal of the astrobiology program is the search for fossil records. The astrobiology exploration strategy calls for the location and return of samples indicative of environments conducive to life, and that best capture and preserve biomarkers. Successfully returning samples from environments conducive to life requires two primary capabilities: (1) in situ mapping of the mineralogy in order to determine whether the desired minerals are present; and (2) nondestructive screening of samples for additional in-situ testing and/or selection for return to laboratories for more in-depth examination. Two of the most powerful identification techniques are micro-imaging and visible/infrared spectroscopy. The design and test results are presented from a compact rugged instrument that combines micro-imaging and spectroscopic capability to provide in-situ analysis, mapping, and sample screening capabilities. Accurate reflectance spectra should be a measure of reflectance as a function of wavelength only. Other compact multispectral microimagers use separate LEDs (light-emitting diodes) for each wavelength and therefore vary the angles of illumination when changing wavelengths. When observing a specularly-reflecting sample, this produces grossly inaccurate spectra due to the variation in the angle of illumination. An advanced design and test results are presented for a multispectral microimager which demonstrates two key advances relative to previous LED-based microimagers: (i) acquisition of actual reflectance spectra in which the flux is a function of wavelength only, rather than a function of both wavelength and illumination geometry; and (ii) increase in the number of spectral bands to eight bands covering a spectral range of 468 to 975 nm.

microimager↗

Current-State Constrained Filter Bank for Wald Testing of Spacecraft Conjunctions

We propose a filter bank consisting of an ordinary current-state extended Kalman filter, and two similar but constrained filters: one is constrained by a null hypothesis that the miss distance between two conjuncting spacecraft is inside their combined hard body radius at the predicted time of closest approach, and one is constrained by an alternative complementary hypothesis. The unconstrained filter is the basis of an initial screening for close approaches of interest. Once the initial screening detects a possibly risky conjunction, the unconstrained filter also governs measurement editing for all three filters, and predicts the time of closest approach. The constrained filters operate only when conjunctions of interest occur. The computed likelihoods of the innovations of the two constrained filters form a ratio for a Wald sequential probability ratio test. The Wald test guides risk mitigation maneuver decisions based on explicit false alarm and missed detection criteria. Since only current-state Kalman filtering is required to compute the innovations for the likelihood ratio, the present approach does not require the mapping of probability density forward to the time of closest approach. Instead, the hard-body constraint manifold is mapped to the filter update time by applying a sigma-point transformation to a projection function. Although many projectors are available, we choose one based on Lambert-style differential correction of the current-state velocity. We have tested our method using a scenario based on the Magnetospheric Multi-Scale mission, scheduled for launch in late 2014. This mission involves formation flight in highly elliptical orbits of four spinning spacecraft equipped with antennas extending 120 meters tip-to-tip. Eccentricities range from 0.82 to 0.91, and close approaches generally occur in the vicinity of perigee, where rapid changes in geometry may occur. Testing the method using two 12,000-case Monte Carlo simulations, we found the method achieved a missed detection rate of 0.1%, and a false alarm rate of 2%.

Carpenter, J. Russell↗

Exploring Network-Related Optimization Problems Using Quantum Heuristics

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics.One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Wang, Zhihui↗

Study network-related optimization problems using quantum alternating optimization ansatz

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics. One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Zhihui Wang↗

Validations of Satellite Ozone Profiles in Austral Spring Using Ozonesonde Measurements in the Jang Bogo Station, Antarctica

Using ozonesonde measurements from 2015 to 2018 at the Jang Bogo station located in the southeastern Antarctic region, we evaluate ozone profiles retrieved from the three satellite measurements that are widely used: Ozone Monitoring Instrument (OMI), Microwave Limb Sounder (MLS), and Ozone Mapping Profiler Suite (OMPS) data. For the fair validation, ozonesonde profiles are smoothed using the weighting function of each satellite retrieval algorithm (i.e., convolution process). Compared with limb-viewing MLS and OMPS ozone profiles, the OMI ozone profiles are relatively less qualified: coarser vertical resolution and larger inter-annual variation. Nevertheless, our validation reveals that the quality of all three satellite ozone profiles looks comparable; In general, difference from ozonesonde profile is ∼1 ppm absolutely, and −20 to 30% relatively at maximum. This quantitative range well corresponds to previous work, meaning that our new validation confirms the reliability of satellite ozone profiles in the southeastern Antarctic region where the measurement data for the validation were not enough. Another interesting feature is the role of a priori ozone profile; Nadir-viewing OMI satellite can have qualified ozone profiles by a proper assumption of a priori ozone profile. Since the performance of limb-viewing ozone profiles is better, however, the careful usage of nadir-viewing ozone profile is still required. We think that the simultaneous usage of multiple satellite ozone profiles can contribute to better understanding of Antarctic ozone characteristics.

Antarctica↗

Unsteady Aerodynamic Modeling of Atmospheric Entry Vehicles in Subsonic and Incompressible flow: A Frequency Response Approach

The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.

Entry Systems Modeling↗