SEARCH · Engineering Papers

Results for “approximation algorithm”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Millimeter Wave Synthetic Aperture Imaging System with a Unique Rotary Scanning System

In recent years, millimeter wave imaging techniques, using synthetic aperture focusing and holographical approaches, have shown tremendous potential for nondestructive testing applications, involving materials and structures used in space vehicles, including the space shuttle external fuel tank spray on foam insulation and its acreage heat tiles. The ability of signals at millimeter wave frequencies (30 - 300 GHz) to easily penetrate inside of low loss dielectric materials, their relatively small wavelengths, and the possibility of detecting coherent (magnitude and phase) reflections make them suitable for high resolution synthetic aperture focused imaging the interior of such materials and structures. To accommodate imaging requirements, commonly a scanning system is employed that provides for a raster scan of the desired structure. However, most such scanners, although simple in design and construction, are inherently slow primarily due to the need to stop and start at the beginning and end of each scan line. To this end, a millimeter wave synthetic aperture focusing system including a custom-designed transceiver operating at 35 - 45 GHz (Q-band) and unique and complex rotary scanner was designed and developed. The rotary scanner is capable of scanning an area with approximately 80 cm in diameter in less than 10 minutes at step sizes of 3 mm and smaller. The transceiver is capable of producing accurate magnitude and phase of reflected signal from the structure under test. Finally, a synthetic aperture focusing algorithm was developed that translates this rotary-obtained magnitude and phase into a synthetic aperture focusing image of inspected structures. This paper presents the design of the transceiver and the rotary scanning system along with showing several images obtained with this system from various complicated structures.

Ghasr, M. T.↗

NASA NTRS2005

Kurtosis Approach Nonlinear Blind Source Separation

In this paper, we introduce a new algorithm for blind source signal separation for post-nonlinear mixtures. The mixtures are assumed to be linearly mixed from unknown sources first and then distorted by memoryless nonlinear functions. The nonlinear functions are assumed to be smooth and can be approximated by polynomials. Both the coefficients of the unknown mixing matrix and the coefficients of the approximated polynomials are estimated by the gradient descent method conditional on the higher order statistical requirements. The results of simulation experiments presented in this paper demonstrate the validity and usefulness of our approach for nonlinear blind source signal separation Keywords: Independent Component Analysis, Kurtosis, Higher order statistics.

Duong, Vu A.↗

NASA NTRS2005

Kurtosis Approach for Nonlinear Blind Source Separation

kurtosis↗

NASA NTRSDate not supplied

Probabilistic Calibration of Expensive Models using Efficiently Trained Surrogates

Calibration of computational models in the presence of uncertainty is often cast as a Bayesian inference problem and solved via sampling methods, e.g., Markov chain Monte Carlo. When the computational model is expensive, this task becomes intractable due to the large number of samples required to accurately estimate the posterior distribution of the calibration parameters. A popular solution to this problem is to use machine learning to develop a faster-to-evaluate, lower-fidelity substitute for the original model to serve as a surrogate while solving the inference problem. Although considered an offline cost, generating training data to construct this surrogate model can still be an expensive task in practice. An active learning algorithm is presented that focuses training on improving surrogate accuracy specifically in and around the bulk of the posterior distribution, as this is where the model is exercised during calibration. Candidate samples are drawn from families of distributions related to an approximation of the posterior. The sample maximizing predictive variance is then selected for evaluation by the original computational model, yielding a label for the training point. Iterating this approach increases efficiency relative to space filling designs (e.g., Latin hypercube sampling) by avoiding low probability points. Practical considerations are discussed, including the benefits of using a sequential Monte Carlo sampling approach, convergence heuristics, and the importance of both exploration and exploitation given that the true posterior is unknown a priori.

uncertainty quantification↗

NASA NTRS1972

Compression of transmission bandwidth requirements for a certain class of band-limited functions.

A study of source-encoding techniques that afford a reduction of data-transmission rates is made with particular emphasis on the compression of transmission bandwidth requirements of band-limited functions. The feasibility of bandwidth compression through analog signal rooting is investigated. It is found that the N-th roots of elements of a certain class of entire functions of exponential type possess contour integrals resembling Fourier transforms, the Cauchy principal values of which are compactly supported on an interval one N-th the size of that of the original function. Exploring this theoretical result, it is found that synthetic roots can be generated, which closely approximate the N-th roots of a certain class of band-limited signals and possess spectra that are essentially confined to a bandwidth one N-th that of the signal subjected to the rooting operation. A source-encoding algorithm based on this principle is developed that allows the compression of data-transmission requirements for a certain class of band-limited signals.

Smith, I. R.↗

NASA NTRS1972

A curve fitting method for solving the flutter equation

A curve fitting approach was developed to solve the flutter equation for the critical flutter velocity. The psi versus nu curves are approximated by cubic and quadratic equations. The curve fitting technique utilized the first and second derivatives of psi with respect to nu. The method was tested for two structures, one structure being six times the total mass of the other structure. The algorithm never showed any tendency to diverge from the solution. The average time for the computation of a flutter velocity was 3.91 seconds on an IBM Model 50 computer for an accuracy of five per cent. For values of nu close to the critical root of the flutter equation the algorithm converged on the first attempt. The maximum number of iterations for convergence to the critical flutter velocity was five with an assumed value of nu relatively distant from the actual crossover.

Cooper, J. L.↗

NASA NTRS1980

A software-based system which produces SEASAT SAR imagery

A digital processing algorithm and its associated system design for producing images from SEASAT synthetic aperture radar (SAR) data is described. The proposed system uses the fast Fourier transform approach to perform the two dimensional correlation process. The range migration problem can be alleviated by approximating the locus of echoes from a point target by several linear segments. Data corresponding to each segment are correlated separately and the results are coherently summed to produce full resolution images. This processing approach exhibits high computation efficiency and simple processing control functions. It is particularly attractive for software implementation based on general purpose computers. Results of this implementation and examples of digitally correlated SEASAT SAR imagery are discussed.

Wu, C.↗

NASA NTRS1989

Frequency estimation techniques for high dynamic trajectories

A comparison is presented of four different estimation techniques applied to the problem of continuously estimating the rapidly varying parameters of a sinusoidal signal, observed in the presence of additive noise. Frequency estimates are emphasized, although phase and/or frequency rate are also estimated by some of the algorithms. These parameters are related to the velocity, position, and acceleration of the maneuvering receiver or transmitter. Estimated performance at low carrier-to-noise ratios and high dynamics is investigated for the purpose of determining the useful operating range of an approximate maximum-likelihood estimator, an extended Kalman filter, a cross-product automatic frequency loop, and a phase-locked loop. Numerical simulations are used to evaluate performance while tracking a common trajectory exhibiting high dynamics.

Vilnrotter, Victor A.↗

NASA NTRS1989

Acoustical effects of a large ridge on low-frequency sound propagation in stationary and moving atmospheres

The effects of a ridge on a low-frequency acoustic propagation in quiescent and windy atmospheres are investigated using a parabolic approximation. A logarithmic wind-speed profile, commonly employed to model atmospheric wind currents, is modified and used to model two-dimensional atmospheric flow over a triangularly-shaped hill. The parabolic equation is solved using an implicit finite-difference algorithm. Several examples are examined to determine the combined effects of source-ridge distance, ridge dimensions, wind-speed profile, and CW source frequency on the received acoustic field.

Robertson, J. S.↗

NASA NTRS1990

Gamma guidance of trajectories for coplanar, aeroassisted orbital transfer

The optimization and guidance of trajectories for coplaner, aeroassisted orbital transfer (AOT) from high Earth orbit (HEO) to low Earth orbit (LEO) are examined. In particular, HEO can be a geosynchronous Earth orbit (GEO). It is assumed that the initial and final orbits are circular, that the gravitational field is central and is governed by the inverse square law, and that at most three impulses are employed: one at HEO exit, one at atmospheric exit, and one at LEO entry. It is also assumed that, during the atmospheric pass, the trajectory is controlled via the lift coefficient. The presence of upper and lower bounds on the lift coefficient is considered. First, optimal trajectories are computed by minimizing the total velocity impulse (hence, the propellant consumption) required for AOT transfer. The sequential gradient-restoration algorithm (SGRA) is used for optimal control problems. The optimal trajectory is shown to include two branches: a relatively short descending flight branch (branch 1) and a long ascending flight branch (branch 2). Next, attention is focused on guidance trajectories capable of approximating the optimal trajectories in real time, while retaining the essential characteristics of simplicity, ease of implementation, and reliability. For the atmospheric pass, a feedback control scheme is employed and the lift coefficient is adjusted according to a two-stage gamma guidance law. Further improvements are possible via a modified gamma guidance which is more stable with respect to dispersion effects arising from navigation errors, variations of the atmospheric density, and uncertainties in the aerodynamic coefficients than gamma guidance trajectory. A byproduct of the studies on dispersion effects is the following design concept. For coplaner aeroassisted orbital transfer, the lift-range-to-weight ratio appears to play a more important role than the lift-to-drag ratio. This is because the lift-range-to-weight ratio controls mainly the minimum altitude (hence, the peak heating rate) of the guidance trajectory; on the other hand, the lift-to-drag ratio controls mainly the duration of the atmospheric pass of the guidance trajectory.

Miele, A.↗

NASA NTRS1990

Real-time simulation of supersonic inlets

A previously published real-time simulation algorithm, the matrix stability region placement (MSRP) method, is used to simulate a small perturbation model of the NASA Lewis Mach 2.5 40-60 mixed compression inlet. The model is representative of high-speed internal flow propulsion systems which can be approximated as quasi-one-dimensional flows. The resulting system of equations, which is stiff, is also simulated by the second-order Adam-Bashforth method. It is shown that the MSRP method can be used to simulate small perturbation models of high-speed internal flow propulsion systems in real time.

Mossayebi, F.↗

NASA NTRS1989

Totally parallel multilevel algorithms for sparse elliptic systems

The fastest known algorithms for the solution of a large elliptic boundary value problem on a massively parallel hypercube all require O(log(n)) floating point operations and O(log(n)) distance-1 communications, if massively parallel is defined to mean a number of processors proportional to the size n of the problem. The Totally Parallel Multilevel Algorithm (TPMA) that has, as special cases, four of these fast algorithms is described. These four algorithms are Parallel Superconvergent Multigrid (PSMG), Robust Multigrid, the Fast Fourier Transformation (FFT) based Spectral Algorithm, and Parallel Cyclic Reduction. The algorithm TPMA, when described recursively, has four steps: (1) project to a collection of interlaced, coarser problems at the next lower level; (2) apply TPMA, recursively, to each of these lower level problems, solving directly at the lowest level; (3) interpolate these approximate solutions to the finer grid, and to verage them to form an approximate solution on this grid; and (4) refine this approximate solution with a defect-correction step, using a local approximate inverse. Choice of the projection operator (P), the interpolation operator (Q), and the smoother (S) determines the class of problems on which TPMA is most effective. There are special cases in which the first three steps produce an exact solution, and the smoother is not needed (e.g., constant coefficient operators).

Frederickson, Paul O.↗

NASA NTRS2007

The Plane-parallel Albedo Bias of Liquid Clouds from MODIS Observations

In our most advanced modeling tools for climate change prediction, namely General Circulation Models (GCMs), the schemes used to calculate the budget of solar and thermal radiation commonly assume that clouds are horizontally homogeneous at scales as large as a few hundred kilometers. However, this assumption, used for convenience, computational speed, and lack of knowledge on cloud small scale variability, leads to erroneous estimates of the radiation budget. This paper provides a global picture of the solar radiation errors at scales of approximately 100 km due to warm (liquid phase) clouds only. To achieve this, we use cloud retrievals from the instrument MODIS on the Terra and Aqua satellites, along with atmospheric and surface information, as input into a GCM-style radiative transfer algorithm. Since the MODIS product contains information on cloud variability below 100 km we can run the radiation algorithm both for the variable and the (assumed) homogeneous clouds. The difference between these calculations for reflected or transmitted solar radiation constitutes the bias that GCMs would commit if they were able to perfectly predict the properties of warm clouds, but then assumed they were homogeneous for radiation calculations. We find that the global average of this bias is approx.2-3 times larger in terms of energy than the additional amount of thermal energy that would be trapped if we were to double carbon dioxide from current concentrations. We should therefore make a greater effort to predict horizontal cloud variability in GCMs and account for its effects in radiation calculations.

Oreopoulos, Lazaros↗

NASA NTRS2004

Software for Generating Strip Maps from SAR Data

Jurassicprok is a computer program that generates strip-map digital elevation models and other data products from raw data acquired by an airborne synthetic-aperture radar (SAR) system. This software can process data from a variety of airborne SAR systems but is designed especially for the GeoSAR system, which is a dual-frequency (P- and X-band), single-pass interferometric SAR system for measuring elevation both at the bare ground surface and top of the vegetation canopy. Jurassicprok is a modified version of software developed previously for airborne-interferometric- SAR applications. The modifications were made to accommodate P-band interferometric processing, remove approximations that are not generally valid, and reduce processor-induced mapping errors to the centimeter level. Major additions and other improvements over the prior software include the following: a) A new, highly efficient multi-stage-modified wave-domain processing algorithm for accurately motion compensating ultra-wideband data; b) Adaptive regridding algorithms based on estimated noise and actual measured topography to reduce noise while maintaining spatial resolution; c) Exact expressions for height determination from interferogram data; d) Fully calibrated volumetric correlation data based on rigorous removal of geometric and signal-to-noise decorrelation terms; e) Strip range-Doppler image output in user-specified Doppler coordinates; f) An improved phase-unwrapping and absolute-phase-determination algorithm; g) A more flexible user interface with many additional processing options; h) Increased interferogram filtering options; and i) Ability to use disk space instead of random- access memory for some processing steps.

Hensley, Scott↗

NASA NTRSDate not supplied

Active Thermography Based on Tensor Rank Decomposition

Principal Component Thermography applies Singular Value Decomposition (SVD) to post-process data that are derived from active thermographic inspections. SVD provides useful compression of the data and allows for better understanding of substructure and indications of potential damage. In the standard approach, SVD is applied to a certain reshaping of a three-dimensional data stack into a two-dimensional array. This work applies the CANDECOMP-PARAFAC (CP) tensor rank decomposition directly to the three-dimensional data to avoid the initial reshaping step in order to begin to develop an inspection method that can more accurately detect defects in non-homogeneous and anisotropic materials. Tests against simulated data that compare the CP decomposition method with traditional Principal Component Thermography based on SVD are described. Finally, the method of Proper Generalized Decomposition (PGD) is used to derive the CP decomposition, and its performance against other algorithms is also discussed.

Thermography↗

NASA NTRS2022

Inferring iron-oxide species content in atmospheric mineral dust from DSCOVR EPIC observations

The iron-oxide content of dust in the atmosphere and most notably its apportionment between hematite (α-Fe2O3) and goethite (α-FeOOH) are key determinants in quantifying dust’s light absorption, its top of atmosphere UV radiances used for dust monitoring, and ultimately shortwave dust direct radiative effects (DRE). Hematite and goethite column mass concentrations and iron-oxide mass fractions of total dust mass concentration were retrieved from the DeepSpace Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) measurements in the ultraviolet–visible (UV–Vis) channels. The retrievals were performed for dust-identified aerosol plumes over land using aerosol optical depth (AOD) and spectral imaginary refractive index provided by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm over six continental regions (North America, North Africa, West Asia, Central Asia, EastAsia, and Australia). The dust particles are represented as an internal mixture of non-absorbing host and absorbing hematite and goethite. We use the Maxwell–Garnett effective medium approximation with carefully selected complex refractive indices of hematite and goethite that produce mass fractions of iron oxides species consistent with in situ values found in the literature to derive the hematite and goethite volumetric/mass concentrations from MAIAC EPIC products. We compared the retrieved hematite and goethite concentrations with in situ dust aerosol mineralogical content measurements, as well as with published data. Our data display variations within the published range of hematite, goethite, and iron-oxide mass fractions for pure mineral dust cases. A specific analysis is presented for 15 sites over the main dust source regions. Sites in the central Sahara, Sahel, and Middle East exhibit a greater temporal variability of iron oxides relative to other sites. Niger site(13.52°N, 2.63°E) is dominated by goethite over Harmattan season with median of ~2 weight percentage (wt.%) of iron-oxide. Saudi Arabia site (27.49°N, 41.98°E) over Middle East also exhibited surge of goethite content with the beginning of Shamal season. The Sahel dust is richer in iron-oxide than Saharan and northern China dust except in Summer. The Bodélé Depression area shows a distinctively lower iron-oxide concentration (~1 wt.%) throughout the year. Finally, we show that EPIC data allow to constrain the hematite refractive index. Specifically, we select 5 out of 13 different number of hematite refractive indices widely variable in published laboratory studies by constraining the iron-oxide mass ratio to the known measured values. Provided climatology of hematite and goethite mass fractions across main dust regions of the Earth will be useful for dust shortwave DRE studies and climate modeling.

Sujung Go↗

NASA NTRS2012

NASA's Robotic Lunar Lander Development Project

Since early 2005, NASA's Robotic Lunar Lander Development (RLLD) office at NASA MSFC, in partnership with the Applied Physics Laboratory (APL), has developed mission concepts and preformed risk-reduction activities to address planetary science and exploration objectives uniquely met with landed missions. The RLLD team developed several concepts for lunar human-exploration precursor missions to demonstrate precision landing and in-situ resource utilization, a multi-node lunar geophysical network mission, either as a stand-alone mission, or as part of the International Lunar Network (ILN), a Lunar Polar Volatiles Explorer and a Mercury lander mission for the Planetary Science decadal survey, and an asteroid rendezvous and landing mission for the Exploration Precursor Robotics Mission (xPRM) office. The RLLD team has conducted an extensive number of risk-reduction activities in areas common to all lander concepts, including thruster testing, propulsion thermal control demonstration, composite deck design and fabrication, and landing leg stability and vibration. In parallel, the team has developed two robotic lander testbeds providing closed-loop, autonomous hover and descent activities for integration and testing of flight-like components and algorithms. A compressed-air test article had its first flight in September 2009 and completed over 150 successful flights. This small test article (107 kg dry/146 kg wet) uses a central throttleable thruster to offset gravity, plus 3 descent thrusters (~37lbf ea) and 6 attitude-control thrusters (~12lbf ea) to emulate the flight system with pulsed operation over approximately 10s of flight time. The test article uses carbon composite honeycomb decks, custom avionics (COTS components assembled in-house), and custom flight and ground software. A larger (206 kg dry/322 kg wet), hydrogen peroxide-propelled vehicle began flight tests in spring 2011 and fly over 30 successful flights to a maximum altitude of 30m. The monoprop testbed also uses a central gravity-canceling thruster and 3 descent thrusters, but has 12 attitude-control thrusters and a maximum flight time of over a minute. The testbed uses aluminum ortho-grid decks, an LN200-1 IMU, Roke Manor Radar Altimeter, Illunis optical cameras, Novatel Pro-Pak GPS truth data system, Pressure transducers & thermocouples for housekeeping, "In-Control" ground system software, and the core Flight Executive (cFE) modular software environment. The peroxide lander testbed is able to accept other sensors and algorithms for testing, both from within NASA and from other customers. Through these activities, the RLLD team has significantly reduced technical risks for all small and medium class robotic landers for the Moon and other airless planetary bodies.

Cohen, Barbara A.↗

NASA NTRS1992

On the Retrieval and Analysis of Multilevel Clouds

An accurate satellite retrieval of cloud properties depends upon the detection and analysis of multilayered, overlapping cloud systems that surface observations show to be common. Multiple cloud layers are often found, for instance, in frontal situations, where cirrus overlays boundary layer convective cloud or low-to mid-level stratus cloud. Surface observers (Hahan et al., 1982) indicate that over ocean in the Northern Hemisphere between 30 deg. N and 60 deg. N, 51 percent of observations are of multilevel clouds. A satellite analysis by Coakley (1983) over the Pacific Ocean finds that more than 50 percent of 500 (250 sq km) frames exhibit evidence of multilayered cloud systems. The questions addressed in this study are the following: What error is introduced when inferring the cloud pressure from a Field-Of-View (FOV) that contains some arbitrary amount of transparent cloud overlaying a lower-level black cloud, such as stratus, by making the assumption that there is only a single cloud layer in the FOV, and what may be done to improve the cloud retrieval? The CO2 slicing methods (e.g. McCleese and Wilson, 1976; Smith and Platt, 1978; Chahine, 1974) have been shown to provide accurate means of inferring cirrus cloud altitude from passive infrared radiance measurements. The CO2 techniques have been applied to radiometric data from several instruments, notably the High Resolution Infrared Radiometric Sounder (HIRS/2, hereafter referred to as HIRS), the VISSR Atmospheric Sounder (VAS) (e.g., Menzel et al., 1983; Wylie and Menzel, 1989), and most recently to the High Resolution Interferometer Sounder (HIS) (Smith and Frey, 1990). The methods take advantage of the fact that infrared CO2 sounding channels spaced closely in wavenumber each have varying opacity to CO2, thereby causing each channel to be sensitive to a different level in the atmosphere. The techniques have been shown to be effective for single-layered, nonblack, mid- to high-level clouds such as cirrus, but are generally applied operationally to any given cloud occurrence. The CO2 slicing algorithms are most accurate for clouds than occur in a single, well-defined layer, or for multi-layered cloud cases in which the uppermost cloud layer is nearly black. Significant cloud height retrieval errors may ensue if the HIRS Field-Of-View (FOV) is cotaminated with low cloud. McCleese and Wilson (1976) have shown that the retrieved cloud height for the case of multiple cloud layers is a weighted average of the cloud heights actually present. The weight is approximately proportional to the product of the cloud heigt and the effective cloud amount. The effect of their result is that the uppermost cloud layer dominates the cloud pressure retrieval. Beyond stating that the higher cloud dominates the cloud pressure retrieval, there is no quantitative information to provide a way of estimating the errors in cloud pressure retrieval one should expect for certain common multilevel cloud situations or any suggestions on how to reduce the errors. In this paper we estimate the magnitude of the errors and use a simple algorithm to reduce the errors in optically thin cloud height retrival.

Baum, Bryan A.↗