Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration method”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Multigrid Acceleration of Time-Accurate DNS of Compressible Turbulent Flow

An efficient scheme for the direct numerical simulation of 3D transitional and developed turbulent flow is presented. Explicit and implicit time integration schemes for the compressible Navier-Stokes equations are compared. The nonlinear system resulting from the implicit time discretization is solved with an iterative method and accelerated by the application of a multigrid technique. Since we use central spatial discretizations and no artificial dissipation is added to the equations, the smoothing method is less effective than in the more traditional use of multigrid in steady-state calculations. Therefore, a special prolongation method is needed in order to obtain an effective multigrid method. This simulation scheme was studied in detail for compressible flow over a flat plate. In the laminar regime and in the first stages of turbulent flow the implicit method provides a speed-up of a factor 2 relative to the explicit method on a relatively coarse grid. At increased resolution this speed-up is enhanced correspondingly.

Broeze, Jan↗
NASA NTRS1996

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Kreider, Kevin L.↗
NASA NTRS1996

Wavelet Sparse Approximate Inverse Preconditioners

There is an increasing interest in using sparse approximate inverses as preconditioners for Krylov subspace iterative methods. Recent studies of Grote and Huckle and Chow and Saad also show that sparse approximate inverse preconditioner can be effective for a variety of matrices, e.g. Harwell-Boeing collections. Nonetheless a drawback is that it requires rapid decay of the inverse entries so that sparse approximate inverse is possible. However, for the class of matrices that, come from elliptic PDE problems, this assumption may not necessarily hold. Our main idea is to look for a basis, other than the standard one, such that a sparse representation of the inverse is feasible. A crucial observation is that the kind of matrices we are interested in typically have a piecewise smooth inverse. We exploit this fact, by applying wavelet techniques to construct a better sparse approximate inverse in the wavelet basis. We shall justify theoretically and numerically that our approach is effective for matrices with smooth inverse. We emphasize that in this paper we have only presented the idea of wavelet approximate inverses and demonstrated its potential but have not yet developed a highly refined and efficient algorithm.

Chan, Tony F.↗
NASA NTRS1996

Research in Computational Aeroscience Applications Implemented on Advanced Parallel Computing Systems

Improving the numerical linear algebra routines for use in new Navier-Stokes codes, specifically Tim Barth's unstructured grid code, with spin-offs to TRANAIR is reported. A fast distance calculation routine for Navier-Stokes codes using the new one-equation turbulence models is written. The primary focus of this work was devoted to improving matrix-iterative methods. New algorithms have been developed which activate the full potential of classical Cray-class computers as well as distributed-memory parallel computers.

Wigton, Larry↗
NASA NTRS1997

Correction for the STIS echelle blaze function

Using the Early Release Observation of 9 Comae, we demonstrate an iterative method for correcting Space Telescope Imaging Spectrograph (STIS) echelle spectra for the effects of the echelle ripple. This analytic approach allows the actual spectrum of interest to be used in the determination of its calibration. The late F star 9 Comae is not an ideal candidate for this method, due to the many absorption lines present in its E230M spectrum, yet, given this difficulty, the method still works quite well.

Heap, Sara R.↗
NASA NTRS1960

Calculation of Aerodynamic Loading and Twist Characteristics of a Flexible Wing at Mach Numbers Approaching 1.0 and Comparison with Experiment

An iteration method is presented by which the detailed aerodynamic loading and twist characteristics of a flexible wing with known elastic properties may be calculated. The method is applicable at Mach numbers approaching 1.0 as well as at subsonic Mach numbers. Calculations were made for a wing-body combination; the wing was swept back 45 deg and had an aspect ratio of 4. Comparisons were made with experimental results at Mach numbers from.0.80 to 0.98.

Mugler, John P., Jr.↗
NASA NTRS2000

Nonlinear Transient Problems Using Structure Compatible Heat Transfer Code

The report documents the recent effort to enhance a transient linear heat transfer code so as to solve nonlinear problems. The linear heat transfer code was originally developed by Dr. Kim Bey of NASA Largely and called the Structure-Compatible Heat Transfer (SCHT) code. The report includes four parts. The first part outlines the formulation of the heat transfer problem of concern. The second and the third parts give detailed procedures to construct the nonlinear finite element equations and the required Jacobian matrices for the nonlinear iterative method, Newton-Raphson method. The final part summarizes the results of the numerical experiments on the newly enhanced SCHT code.

Hou, Gene↗
NASA NTRS1997

Correction of Rayleigh Scattering Effects in Cloud Optical Thickness Retrievals

We present results that demonstrate the effects of Rayleigh scattering on the 9 retrieval of cloud optical thickness at a visible wavelength (0.66 Am). The sensor-measured radiance at a visible wavelength (0.66 Am) is usually used to infer remotely the cloud optical thickness from aircraft or satellite instruments. For example, we find that without removing Rayleigh scattering effects, errors in the retrieved cloud optical thickness for a thin water cloud layer (T = 2.0) range from 15 to 60%, depending on solar zenith angle and viewing geometry. For an optically thick cloud (T = 10), on the other hand, errors can range from 10 to 60% for large solar zenith angles (0-60 deg) because of enhanced Rayleigh scattering. It is therefore particularly important to correct for Rayleigh scattering contributions to the reflected signal from a cloud layer both (1) for the case of thin clouds and (2) for large solar zenith angles and all clouds. On the basis of the single scattering approximation, we propose an iterative method for effectively removing Rayleigh scattering contributions from the measured radiance signal in cloud optical thickness retrievals. The proposed correction algorithm works very well and can easily be incorporated into any cloud retrieval algorithm. The Rayleigh correction method is applicable to cloud at any pressure, providing that the cloud top pressure is known to within +/- 100 bPa. With the Rayleigh correction the errors in retrieved cloud optical thickness are usually reduced to within 3%. In cases of both thin cloud layers and thick ,clouds with large solar zenith angles, the errors are usually reduced by a factor of about 2 to over 10. The Rayleigh correction algorithm has been tested with simulations for realistic cloud optical and microphysical properties with different solar and viewing geometries. We apply the Rayleigh correction algorithm to the cloud optical thickness retrievals from experimental data obtained during the Atlantic Stratocumulus Transition Experiment (ASTEX) conducted near the Azores in June 1992 and compare these results to corresponding retrievals obtained using 0.88 Am. These results provide an example of the Rayleigh scattering effects on thin clouds and further test the Rayleigh correction scheme. Using a nonabsorbing near-infrared wavelength lambda (0.88 Am) in retrieving cloud optical thickness is only applicable over oceans, however, since most land surfaces are highly reflective at 0.88 Am. Hence successful global retrievals of cloud optical thickness should remove Rayleigh scattering effects when using reflectance measurements at 0.66 Am.

Wang, Meng-Hua↗
NASA NTRS2000

A Two Colorable Fourth Order Compact Difference Scheme and Parallel Iterative Solution of the 3D Convection Diffusion Equation

A new fourth order compact difference scheme for the three dimensional convection diffusion equation with variable coefficients is presented. The novelty of this new difference scheme is that it Only requires 15 grid points and that it can be decoupled with two colors. The entire computational grid can be updated in two parallel subsweeps with the Gauss-Seidel type iterative method. This is compared with the known 19 point fourth order compact differenCe scheme which requires four colors to decouple the computational grid. Numerical results, with multigrid methods implemented on a shared memory parallel computer, are presented to compare the 15 point and the 19 point fourth order compact schemes.

Zhang, Jun↗
NASA NTRS1995

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗
NASA NTRS2002

Enhancement of Optical Nonlinearities Via Whispering Gallery Mode Splitting

An iterative method is applied to the analysis of N coupled ring-resonators, and the results are compared with multilayered spheres. Splitting of the whispering gallery modes into N higher-Q modes occurs when the round-trip phase shifts in each ring (or layer) are equal, in agreement with results for planar resonators. This mode-splitting is, therefore, a universal phenomenon for resonant structures, and can lead to reduced thresholds for nonlinear optical effects.

Chang, Hongrok↗
NASA NTRS2003

Electron-H Elastic Scattering

Precision calculations for e^{-}-H and e^{-}-He^{+} for S-wave scattering in the elastic region have been carried out using the optical potential approach. This formalism is now extended to e^{-}-H P-wave scattering in the elastic region. The scattering equations are solved by the non-iterative method. Phase shifts are calculated using Hylleraas-type correlation functions up to 84 terms. Results are rigorous lower bounds to the exact phase shifts and they are compared to those obtained in previous calculations.

Bhatia, A. K.↗
NASA NTRS2003

Retrieval of Ozone Column Content from Airborne Sun Photometer Measurements During SOLVE II: Comparison with SAGE III, POAM III,THOMAS and GOME Measurements

During the Second SAGE 111 Ozone Loss and Validation Experiment (SOLVE II), the 14- channel NASA Ames Airborne Trackmg Sunphotometer (AATS-14) was mounted on the NASA DC-8 and successfully measured spectra of total and aerosol optical depth (TOD and AOD) during the sunlit portions of eight science flights. Values of ozone column content above the aircraft have been derived from the AATS-14 data by using a linear least squares method. For each AATS-14 measured TOD spectrum, this method iteratively finds the ozone column content that yields the best match between measured and calculated TOD. The calculations assume the known Chappuis ozone band shape and a three-parameter AOD shape (quadratic in log-log space). Seven of the AATS-14 channels (each employing an interference filter with a nominal full-width at half maximum bandpass of -5 nm) are within the Chappuis band, with center wavelengths between 452.9 nm and 864.5 nm. One channel (604.4 nm) is near the peak, and three channels (499.4, 519.4 and 675.1 nm) have ozone absorption within 30-40% of that at the peak. For the typical DC-8 SOLVE II cruising altitudes of approx. 8-12 km and the background stratospheric aerosol conditions that prevailed during SOLVE 11, absorption of incoming solar radiation by ozone comprised a significant fraction of the aerosol-plus-ozone optical depth measured in the four AATS-14 channels centered between 499.4 and 675.1 nm. Typical AODs above the DC-8 ranged from 0.003-0.008 in these channels. For comparison, an ozone overburden of 0.3 atm-cm (300 DU) translates to ozone optical depths of 0.009,0.014, 0.041, and 0.012, respectively, at these same wavelengths. In this paper, we compare AATS-14 values of ozone column content with temporally and spatially near-coincident values derived from measurements acquired by the Stratospheric Aerosol and Gas Experiment III (SAGE III) and the Polar Ozone and Aerosol Measurement 111 (POAM III) satellite sensors. We also compare AATS-14 ozone retrievals during selected DC-8 latitudinal and longitudinal transects with total column ozone data acquired by the Total Ozone Mapping Spectrometer (TOMS) and the Global Ozone Monitoring Experiment (GOME) satellite sensors. To enable this comparison, the amount of ozone in the column below the aircraft is estimated by combining SAGE and/or POAM data with high resolution, fast response in-situ ozone measurements acquired during the DC-8 ascent at the start of each science flight.

Livingston, J.↗
NASA NTRS1997

Implicit Kalman filtering

For an implicitly defined discrete system, a new algorithm for Kalman filtering is developed and an efficient numerical implementation scheme is proposed. Unlike the traditional explicit approach, the implicit filter can be readily applied to ill-conditioned systems and allows for generalization to descriptor systems. The implementation of the implicit filter depends on the solution of the congruence matrix equation (A1)(Px)(AT1) = Py. We develop a general iterative method for the solution of this equation, and prove necessary and sufficient conditions for convergence. It is shown that when the system matrices of an implicit system are sparse, the implicit Kalman filter requires significantly less computer time and storage to implement as compared to the traditional explicit Kalman filter. Simulation results are presented to illustrate and substantiate the theoretical developments.

Non-NASA Center↗
NASA NTRS2004

Frequency Domain Calculations Of Acoustic Propagation

Two complex geometry problems are solved using the linearized Euler equations. The impedance mismatch method1 is used to impose the solid surfaces without the need to use a body-fitted grid. The problem is solved in the frequency domain to avoid long run times. Although the harmonic assumption eliminates all time dependence, a pseudo-time term is added to allow conventional iterative methods to be employed. A Jameson type, Runge-Kutta scheme is used to advance the solution in pseudo time. The spatial operator is based on a seven-point, sixth-order finite difference. Constant coefficient, sixth-derivative artificial dissipation is used throughout the domain. A buffer zone technique employing a complex frequency to damp all waves near the boundaries is used to minimize reflections. The results show that the method is capable of capturing the salient features of the scattering, but an excessive number of grid points are required to resolve the phenomena in the vicinity of the solid bodies because the wavelength of the acoustics is relatively short compared with the size of the bodies. Smoothly transitioning into the immersed boundary condition alleviates the difficulties, but a fine mesh is still required.

Lockard, David P.↗
NASA NTRS2006

Iterative Magnetometer Calibration

This paper presents an iterative method for three-axis magnetometer (TAM) calibration that makes use of three existing utilities recently incorporated into the attitude ground support system used at NASA's Goddard Space Flight Center. The method combines attitude-independent and attitude-dependent calibration algorithms with a new spinning spacecraft Kalman filter to solve for biases, scale factors, nonorthogonal corrections to the alignment, and the orthogonal sensor alignment. The method is particularly well-suited to spin-stabilized spacecraft, but may also be useful for three-axis stabilized missions given sufficient data to provide observability.

Sedlak, Joseph↗
NASA NTRS2007

Nonlinear Structural Analysis Methodology and Dynamics Scaling of Inflatable Parabolic Reflector Antenna Concepts

Ultra-light weight and ultra-thin membrane inflatable antenna concepts are fast evolving to become the state-of-the-art antenna concepts for deep-space applications. NASA Langley Research Center has been involved in the structural dynamics research on antenna structures. One of the goals of the research is to develop structural analysis methodology for prediction of the static and dynamic response characteristics of the inflatable antenna concepts. This research is focused on the computational studies to use nonlinear large deformation finite element analysis to characterize the ultra-thin membrane responses of the antennas. Recently, structural analyses have been performed on a few parabolic reflector antennas of varying size and shape, which are referred in the paper as 0.3 meters subscale, 2 meters half-scale, and 4 meters full-scale antenna. The various aspects studied included nonlinear analysis methodology and solution techniques, ways to speed convergence in iterative methods, the sensitivities of responses with respect to structural loads, such as inflation pressure, gravity, and pretension loads in the ground and in-space conditions, and the ultra-thin membrane wrinkling characteristics. Several such intrinsic aspects studied have provided valuable insight into evaluation of structural characteristics of such antennas. While analyzing these structural characteristics, a quick study was also made to assess the applicability of dynamics scaling of the half-scale antenna. This paper presents the details of the nonlinear structural analysis results, and discusses the insight gained from the studies on the various intrinsic aspects of the analysis methodology. The predicted reflector surface characteristics of the three inflatable ultra-thin membrane parabolic reflector antenna concepts are presented as easily observable displacement fringe patterns with associated maximum values, and normal mode shapes and associated frequencies. Wrinkling patterns are presented to show how surface wrinkle progress with increasing tension loads. Antenna reflector surface accuracies were found to be very much dependent on the type and size of the antenna, the reflector surface curvature, reflector membrane supports in terms of spacing of catenaries, as well as the amount of applied load.

Sreekantamurthy, Tham↗
NASA NTRS2008

Efficient Kriging via Fast Matrix-Vector Products

Interpolating scattered data points is a problem of wide ranging interest. Ordinary kriging is an optimal scattered data estimator, widely used in geosciences and remote sensing. A generalized version of this technique, called cokriging, can be used for image fusion of remotely sensed data. However, it is computationally very expensive for large data sets. We demonstrate the time efficiency and accuracy of approximating ordinary kriging through the use of fast matrixvector products combined with iterative methods. We used methods based on the fast Multipole methods and nearest neighbor searching techniques for implementations of the fast matrix-vector products.

Memarsadeghi, Nargess↗