Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative 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 343 records · Page 19

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing/acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1) the elimination of large matrix storage requirements, and (2) shorter solution times under most conditions.

Baumeister, K. J.↗

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗

Influence of exit impedance on finite difference solutions of transient acoustic mode propagation in ducts

The time-dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard-wall duct without flow and with spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the steady solutions of the acoustic field. The transient method did not converge to the steady-state solution for cutoff acoustic duct modes. This has implications as to its use in a variable-area duct, where modes may become cutoff in the smal-area portion of the duct. For single cutoff mode propagation the steady-state impedance boundary condition produced acoustic reflections during the initial transient that caused finite instabilities in the numerical calculations. The stability problem is resolved by reformulating the exit boundary condition. Example calculations show good agreement with exact analytical and numerical results for forcing frequencies above, below, and nearly at the cutoff frequency.

Baumeister, K. J.↗

Sensitivity of constrained linear inversions to the selection of the Lagrange multiplier

The influence of the choice of the Lagrange multiplier on constrained linear inversions is explored, with reference made to applications in inferring the columnar aerosol size distributions from spectral aerosol optical depth measurements. A range of the Lagrange multiplier is examined to find all positive solutions for the solution vector, which represents modifying factors to the assumed form of the size distribution. An iterative method is devised to constrain the calculations to consideration of only positive quantities and a requirement that the regression fit to data be consistent with measurement errors. The determination of the variances and covariances is formulated and applied to existing data sets for optical depth. Variances in the solution are found to be large for particle radii when the information content of the data is small.

King, M. D.↗

Mathematical enhancement of data from scientific measuring instruments

The accuracy of any physical measurement is limited by the instruments performing it. The proposed activities of this grant are related to the study of and application of mathematical techniques of deconvolution. Two techniques are being investigated: an iterative method and a function continuation Fourier method. This final status report describes the work performed during the period July 1 to December 31, 1982.

Ioup, J. W.↗

Accelerating an iterative process by explicit annihilation

A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.

Jespersen, D. C.↗

Experimental and finite element study of a standing torus under normal and tangential loads

The effect of a combined external load consisting of a normal component and a tangential braking force applied to an inflated torus is addressed. Experimental results obtained by photographic study of the contact area between the torus and load plate are presented as well as measurements of the vertical and horizontal displacement of the torus under load. A numerical procedure for displacement analysis is developed in which the finite element program STAGS (Structural Analysis of General Shells) is used in an iterative manner to produce a flat, horizontal footprint surface under force loading. The redistributed force distribution obtained by the iterative process is displayed along with computed meridional and circumferential stresses. Finally, an extension of the iterative method is introduced which eliminates the need to experimentally determine the footprint area.

Flugrad, D. R.↗

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.

Phillips, T. N.↗

Adapting iterative algorithms for solving large sparse linear systems for efficient use on the CDC CYBER 205

Adapting and designing mathematical software to achieve optimum performance on the CYBER 205 is discussed. Comments and observations are made in light of recent work done on modifying the ITPACK software package and on writing new software for vector supercomputers. The goal was to develop very efficient vector algorithms and software for solving large sparse linear systems using iterative methods.

Kincaid, D. R.↗

The use of satellite data in understanding and predicting convective and large-scale dynamical processes

Mesoscale convective processes and how they affect and interact with mid-latitude cyclones were studied. The ageostrophic and associated vertical motion field was calculated using a highly accurate iterative method of solving the semigeostrophic omega equation. The tendencies for convective destabilization in the 850-750 mb layer due to differential geostrophic and ageostrophic advection and differential moist adiabatic ascent, were found. The spectral models of the index oscillation, one barotropic and the other baroclinic, were developed. Theoretical and observational studies of cloud streets were conducted.

Dutton, J. A.↗

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem. Previously announced in STAR as N83-30109

Phillips, T. N.↗

Fidelity of Satellite Soundings

Two pattern recognition procedures to provide improvements to first guess fields for satellite temperature retrievals were developed. The first technique where a radiometer measurement is used to select one or more historical radiosonde temperature profiles as analog estimates of ambient thermal structure. Radiosondes are the vertical scales of the analogs. The vertical resolving power of the satellite radiometer is relevant only to a decision process. The analog selection process is much more effective if implemented in an orthogonalized space of measurement information. The second procedure partitions a priori dependent data into shape coherent pattern libraries using structure information inherent in the data itself. It is shown that the pattern recognition techniques are able to reduce first guess profile errors by nearly 50%. The impact of pattern recognition on temperature retrieval error is assessed by regression and physical iterative retrieval algorithms. The influence of improved first guess fields is markedly different on these two types of algorithm. Pattern recognition was a strong positive impact on the physical iterative method.

Thompson, O.↗

Accelerating an iterative process by explicit annihilation

A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.

Jespersen, D. C.↗

Boundary integral techniques for multi-connected domains

Attention is given to iterative methods applicable to such changing domain cases of elliptic problems in multiconnected domains as those of incompressible inviscid flow with free surfaces. One such formulation is developed and tested on circular and elliptic annuli, using interpolated quadrature points to maintain accuracy when annuli regions are thin. The technique is noted to be applicable to the study of accelerating thin fluid shells.

Baker, G. R.↗

Progress and supercomputing in computational fluid dynamics; Proceedings of U.S.-Israel Workshop, Jerusalem, Israel, December 1984

Current developments and future trends in the application of supercomputers to computational fluid dynamics are discussed in reviews and reports. Topics examined include algorithm development for personal-size supercomputers, a multiblock three-dimensional Euler code for out-of-core and multiprocessor calculations, simulation of compressible inviscid and viscous flow, high-resolution solutions of the Euler equations for vortex flows, algorithms for the Navier-Stokes equations, and viscous-flow simulation by FEM and related techniques. Consideration is given to marching iterative methods for the parabolized and thin-layer Navier-Stokes equations, multigrid solutions to quasi-elliptic schemes, secondary instability of free shear flows, simulation of turbulent flow, and problems connected with weather prediction.

Murman, E. M.↗

Torsional vibrations and dynamic loads in a basic planetary gear system

An iterative method has been developed for analyzing dynamic loads in a light weight basic planetary gear system. The effects of fixed, semi-floating, and fully-floating sun gear conditions have been emphasized. The load dependent variable gear mesh stiffness were incorporated into a practical torsional dynamic model of a planetary gear system. The dynamic model consists of input and output units, shafts, and a planetary train. In this model, the sun gear has three degrees of freedom; two transverse and one rotational. The planets, ring gear, and the input and output units have one degree of freedom, (rotation) thus giving a total of nine degrees of freedoms for the basic system. The ring gear has a continuous radial support. The results indicate that the fixed sun gear arrangement with accurate or errorless gearing offers in general better performance than the floating sun gear system.

August, R.↗

A three dimensional calculation of elastic equilibrium for composite materials

A compact scheme is applied to three-dimensional elasticity problems for composite materials, involving simple geometries. The mathematical aspects of this approach are discussed, in particular the iteration method. A vector processor code implementing the compact scheme is presented, and several numerical experiments are summarized.

Lustman, Liviu R.↗

Rapid computation of chemical equilibrium composition - An application to hydrocarbon combustion

A scheme for rapidly computing the chemical equilibrium composition of hydrocarbon combustion products is derived. A set of ten governing equations is reduced to a single equation that is solved by the Newton iteration method. Computation speeds are approximately 80 times faster than the often used free-energy minimization method. The general approach also has application to many other chemical systems.

Erickson, W. D.↗