Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative regularization”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗
NASA NTRS2009

Differential Absorption Lidar to Measure Sub-Hourly Variation of Tropospheric Ozone Profiles

A tropospheric ozone Differential Absorption Lidar (DIAL) system, developed jointly by the University of Alabama at Huntsville and NASA, is making regular observations of ozone vertical distributions between 1 and 8 km with two receivers under both daytime and nighttime conditions using lasers at 285 and 291 nm. This paper describes the lidar system and analysis technique with some measurement examples. An iterative aerosol correction procedure reduces the retrieval error arising from differential aerosol backscatter in the lower troposphere. Lidar observations with coincident ozonesonde flights demonstrate that the retrieval accuracy ranges from better than 10% below 4 km to better than 20% below 8 km with 750-m vertical resolution and 10-min temporal integration

Kuang, Shi↗
NASA NTRS2011

Differential Absorption Lidar to Measure Subhourly Variation of Tropospheric Ozone Profiles

A tropospheric ozone Differential Absorption Lidar system, developed jointly by The University of Alabama in Huntsville and the National Aeronautics and Space Administration, is making regular observations of ozone vertical distributions between 1 and 8 km with two receivers under both daytime and nighttime conditions using lasers at 285 and 291 nm. This paper describes the lidar system and analysis technique with some measurement examples. An iterative aerosol correction procedure reduces the retrieval error arising from differential aerosol backscatter in the lower troposphere. Lidar observations with coincident ozonesonde flights demonstrate that the retrieval accuracy ranges from better than 10% below 4 km to better than 20% below 8 km with 750-m vertical resolution and 10-min 17 temporal integration.

Kuang, Shi↗
NASA NTRS1984

Dynamic characteristics of a hydrostatic gas bearing driven by oscillating exhaust pressure

Vibration of a statically loaded, inherently compensated hydrostatic journal bearing due to oscillating exhaust pressure is investigated. Both angular and radial vibration modes are analyzed. The time-dependent Reynolds equation governing the pressure distribution between the oscillating journal and sleeve is solved together with the journal equation of motion to obtain the response characteristics of the bearing. The Reynolds equation and the equation of motion are simplified by applying regular perturbation theory for small displacements. The numerical solutions of the perturbation equations are obtained by discretizing the pressure field using finite-difference aproximations with a discrete, nonuniform line-source model which excludes effects due to feeding hole volume. An iterative scheme is used to simultaneously satisfy the equations of motion for the journal. The results presented include Bode plots of bearing-oscillation gain and phase for a particular bearing configuration for various combinations of parameters over a range of frequencies, including the resonant frequency.

Watkins, C. B.↗
NASA NTRS2004

Accumulate repeat accumulate codes

In this paper we propose an innovative channel coding scheme called Accumulate Repeat Accumulate codes. This class of codes can be viewed as trubo-like codes, namely a double serial concatenation of a rate-1 accumulator as an outer code, a regular or irregular repetition as a middle code, and a punctured accumulator as an inner code.

LDPC codes turbo-like codes iterative decoding↗
NASA NTRS2018

Global High Resolution Crustal Magnetic Field at the Surface of the Moon from Low-Altitude Lunar Prospector Magnetic Gradient Data

We derived new vector gradients based models of crustal magnetic field at the lunar surface with data from the Lunar Prospector (LP) satellite using two model parameterization approaches: a global set of 35820 1° spaced (~30 km) equal area monopoles at 20 km below the surface (O’Brien and Parker, 1994; Olsen et al., 2017) and combined results of subsets of 100000 0.66° spaced monopoles at the same depth. We use the scheme of iteratively reweighted least-squares inversion to compute the initial model. Then the amplitudes of these monopoles are determined by minimizing the misfit to the components together with the global average of |Br| at the ellipsoid surface (i.e. applying a L1 model regularization of Br). In previous approaches using vector fields for modeling, we found that external field contamination leads to spurious anomalies in the downward continued field models even with stringent data selection criteria and ad-hoc noise removal techniques (e.g., satellite’s position in the Moon’s wake w.r.t. the solar wind and in the Earth’s magnetotail, internal/external dipoles fields removal, low-order polynomial removal, joint equivalent source cross-validation technique and visually removing remaining anomalous segments). On the other hand, with the use of gradients-only data (along-track first finite differences), we were able to completely bypass the ad-hoc techniques. Similar processing of Kaguya magnetic data, which have only higher altitude coverage over the most of the Moon except in the region of the South Pole–Aitken (SPA) basin, completely misses some of the anomalies seen in the Lunar Prospector data. The combined Lunar Prospector and Kaguya gradient-based models also severely degrade the derivation of the anomaly fields in many regions. Euler analysis of isolated anomaly features from the Reiner Gamma swirl suggests top depths of about 0.3 to 1.5 km and center depths of 10-14 km; in the region between Stein and Vallier craters north of the SPA basin our analysis suggests top depths of around 1.5 km and center depths of 13-15 km. With the spectral depth determination techniques, the SPA basin region yields depths to the base of magnetization ranging between 15 and 40 km. Three-dimensional modeling and the bulk magnetization determinations of the sources constrained by the Euler and spectral methods is underway.

Gradient↗
NASA NTRS1983

Computation of viscous flows over airfoils, including separation, with a coupling approach

Viscous incompressible flows over single or multiple airfoils, with or without separation, were computed using an inviscid flow calculation, with modified boundary conditions, and by a method providing calculation and coupling for boundary layers and wakes, within conditions of strong viscous interaction. The inviscid flow is calculated with a method of singularities, the numerics of which were improved by using both source and vortex distributions over profiles, associated with regularity conditions for the fictitious flows inside of the airfoils. The viscous calculation estimates the difference between viscous flow and inviscid interacting flow, with a direct or inverse integral method, laminar or turbulent, with or without reverse flow. The numerical method for coupling determines iteratively the boundary conditions for the inviscid flow. For attached viscous layers regions, an underrelaxation is locally calculated to insure stability. For separated or separating regions, a special semi-inverse algorithm is used. Comparisons with experiments are presented.

Leballeur, J. C.↗
NASA NTRS1983

A time-dependent approach for calculating steady inverse boundary-layer flows with separation

An unsteady inverse boundary-layer method is developed which can be used to calculate steady flows with separation. Two versions of Keller's box method with the Mechul function formulation developed by Cebeci (1976) are employed, depending on the complexity of the flow. The regular box is employed in regions of positive streamwise velocity component u, whereas the zigzag box is employed in regions where u becomes negative (t greater than 0). The regular box with the FLARE approximation is employed when t = 0 and u becomes negative in some region across the layer. Results of calculations show that the use of a time-dependent inverse boundary-layer method in which time is used as an iteration parameter provides a good approach in improving the accuracy of the solutions obtained from the FLARE approximation.

Cebeci, T.↗
NASA NTRS2003

Universal Behavior of a Cyclic Oxidation Model

A mathematical model has been generated to represent the iterative, discrete growth and spallation processes associated with cyclic oxidation. Parabolic growth kinetics (k(sub p)) over and a constant spall area (F(sub A)) were assumed, with spalling occurring interfacially at the thickest regions of the scale. Although most models require numerical techniques, the regularity and simplicity of this progression permitted an approximation by algebraic expressions. Normalization could now be performed to reflect all parametric effects, and a universal cyclic oxidation response was generated: W(sub u) = 1/2 {3J(sub u)(sup 1/2)+ J(sub u)(sup 3/2)} where W, is weight change normalized by the maximum and J(sub u) is the cycle number normalized by the number to reach maximum. Similarly, the total amount of metal consumed was represented by a single normalized curve. The factor [(S(sub c)-l)(raised dot)sqrt(F(sub A)k(sub p)DELTAt)] was identified as a general figure of merit, where S(sub c) is the mass ratio of oxide to oxygen and DELTAt is the cycle duration. A cyclic oxidation failure map was constructed, in normalized k(sub p)-F(sub A) space, as defined by the locus of points corresponding to a critical amount of metal consumption in a given time. All three constructions describe behavior for every value of growth rate, spall fraction, and cycle duration by means of single curves, but with two branches corresponding to the times before and after steady state is achieved.

Smialek, James L.↗
NASA NTRS1990

Implementation of two projection methods on a shared memory multiprocessor - DEC VAX 6240

The relative performance of two iterative schemes, based on projection techniques, is compared on a shared memory multiprocessor - VAX 6240. The CG accelerated Block-SSOR method and the CG accelerated Symmetric-Kaczmarz method are considered for the solution of large sparse nonsymmetric systems of linear equations. It is shown that the regular structure of many matrices can be exploited by the CG-accelerated Block-SSOR method to provide good speedup in a multiprocessing environment. However, the CG accelerated Symmetric-Kaczmarz method, while being a viable alternative on a scalar machine, is unable to benefit from multiprocessing.

Kamath, C.↗
NASA NTRS2011

A New Global Mascon Solution Tuned for High-Latitude Ice Studies

A new global mascon solution has been developed with I-arc-degree spatial and IO-day temporal sampling. The global mas cons are estimated from the reduction of nearly 8 years of GRACE K-band range-rate data. Temporal and anisotropic spatial constraints have been applied for land, ocean and ice regions. The solution construction and tuning is focused towards the Greenland and Antarctic ice sheets (GIS and AIS) as well as the Gulf of Alaska mountain glaciers (GoA). Details of the solution development will be discussed, including the mascon parameter definitions, constraints, and the tuning of the constraint damping factor. Results will be presented, exploring the spatial and temporal variability of the ice sheets and GoA regions. A detailed error analysis will be discussed, including solution dependence on iteration, damping factor, forward modeling, and multitechnique comparisons. We also investigate the fundamental resolution of the solution and the spatial correlation of ice sheet inter-annual change. Finally, we discuss future improvements, including specific constraint application for the rest of the major land ice regions and improvements in solution regularization.

Luthcke, S. B.↗
NASA NTRS1988

Stability sensitivity analysis for the aeroelastic optimization of a helicopter rotor

A sensitivity study of blade stability in forward flight for a hingeless rotor with respect to structural design variables is carried out using a direct analytical method. Structural design variables include nonstructural mass distribution (spanwise and chordwise), chordwise offset of center of gravity, and blade bending stiffnesses (flap, lag and torsion). The formulation for blade steady response is based on a finite element method in space and time. The vehicle trim and blade steady response are calculated iteratively as one coupled solution using a modified Newton method. Eigenvalues corresponding to different blade modes are calculated using Floquet transition matrix theory. The formulation for derivatives of the eigenvalues with respect to design variables is implemented using a direct analytical approach (chain rule differentiation), and constitutes an integral part of the regular stability analysis. The stability sensitivity derivatives were obtained at a fraction of computation time compared to the frequently adopted finite difference method. A parametric study showed that nonstructural mass and chordwise cg offset of outboarad elements, and lag bending stiffness of inboard elements, have powerful influence on blade stability.

Lim, Joon W.↗
NASA NTRS1987

Simulated annealing in networks for computing possible arrangements for red and green cones

Attention is given to network models in which each of the cones of the retina is given a provisional color at random, and then the cones are allowed to determine the colors of their neighbors through an iterative process. A symmetric-structure spin-glass model has allowed arrays to be generated from completely random arrangements of red and green to arrays with approximately as much disorder as the parafoveal cones. Simulated annealing has also been added to the process in an attempt to generate color arrangements with greater regularity and hence more revealing moirepatterns than than the arrangements yielded by quenched spin-glass processes. Attention is given to the perceptual implications of these results.

Ahumada, Albert J., Jr.↗
NASA NTRS1992

A VLSI implementation of DCT using pass transistor technology

A VLSI design for performing the Discrete Cosine Transform (DCT) operation on image blocks of size 16 x 16 in a real time fashion operating at 34 MHz (worst case) is presented. The process used was Hewlett-Packard's CMOS26--A 3 metal CMOS process with a minimum feature size of 0.75 micron. The design is based on Multiply-Accumulate (MAC) cells which make use of a modified Booth recoding algorithm for performing multiplication. The design of these cells is straight forward, and the layouts are regular with no complex routing. Two versions of these MAC cells were designed and their layouts completed. Both versions were simulated using SPICE to estimate their performance. One version is slightly faster at the cost of larger silicon area and higher power consumption. An improvement in speed of almost 20 percent is achieved after several iterations of simulation and re-sizing.

Kamath, S.↗
NASA NTRS1981

An empirical investigation of methods for nonsymmetric linear systems

The present investigation is concerned with a comparison of methods for solving linear algebraic systems which arise from finite difference discretizations of the elliptic convection-diffusion equation in a planar region Omega with Dirichlet boundary conditions. Such linear systems are typically of the form Ax = b where A is an N x N sparse nonsymmetric matrix. In a discussion of discretizations, it is assumed that a regular rectilinear mesh of width h has been imposed on Omega. The discretizations considered include central differences, upstream differences, and modified upstream differences. Six methods for solving Ax = b are considered. Three variants of Gaussian elimination have been chosen as representatives of state-of-the-art software for direct methods under different assumptions about pivoting. Three iterative methods are also included.

Sherman, A. H.↗
NASA NTRS1988

Analysis of spectra using correlation functions

A novel method is presented for the quantitative analysis of spectra based on the properties of the cross correlation between a real spectrum and either a numerical synthesis or laboratory simulation. A new goodness-of-fit criterion called the heteromorphic coefficient H is proposed that has the property of being zero when a fit is achieved and varying smoothly through zero as the iteration proceeds, providing a powerful tool for automatic or near-automatic analysis. It is also shown that H can be rendered substantially noise-immune, permitting the analysis of very weak spectra well below the apparent noise level and, as a byproduct, providing Doppler shift and radial velocity information with excellent precision. The technique is in regular use in the Atmospheric Trace Molecule Spectroscopy (ATMOS) project and operates in an interactive, realtime computing environment with turn-around times of a few seconds or less.

Beer, Reinhard↗
NASA NTRS2000

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absorption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗
NASA NTRS1999

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: ##EQU1## wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absoption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗