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 811 records · Page 45

An actuator extension transformation for a motion simulator and an inverse transformation applying Newton-Raphson's method

A set of equations which transform position and angular orientation of the centroid of the payload platform of a six-degree-of-freedom motion simulator into extensions of the simulator's actuators has been derived and is based on a geometrical representation of the system. An iterative scheme, Newton-Raphson's method, has been successfully used in a real time environment in the calculation of the position and angular orientation of the centroid of the payload platform when the magnitude of the actuator extensions is known. Sufficient accuracy is obtained by using only one Newton-Raphson iteration per integration step of the real time environment.

Dieudonne, J. E.↗

Two-dimensional analytical and experimental performance comparison for a compressor stator section with D-factor of 0.47

Analytically computed flow parameters were compared to measured values for the midspan double-circular-arc section of a stator in subsonic flow. Analytical procedures included calculations for inviscid flow, blade surface boundary layers, and loss coefficients. Comparisons were made at three incidence angles. Methods for prescribing the exit fluid angle in inviscid flow calculations were investigated. Measured loss coefficients were compared to calculated values. Two methods for predicting performance involving iterative use of ideal flow and boundary layer calculations were investigated.

Sanger, N. L.↗

Resizing procedure for structures under combined mechanical and thermal loading

The fully-stressed design (FSD) appears to be the most widely used approach for sizing of flight structures under strength and minimum-gage constraints. Almost all of the experience with FSD has been with structures primarily under mechanical loading as opposed to thermal loading. In this method the structural sizes are iterated with the step size, depending on the ratio of the total stress to the allowable stress. In this paper, the thermal fully-stressed design (TFSD) procedure developed for problems involving substantial thermal stress is extended to biaxial stress members using a Von Mises failure criterion. The TFSD resizing procedure for uniaxial stress is restated and the new procedure for biaxial stress members is developed. Results are presented for an application of the two procedures to size a simplified wing structure.

Adelman, H. M.↗

Dipole and quadrupole synthesis of electric potential fields

A general technique for expanding an unknown potential field in terms of a linear summation of weighted dipole or quadrupole fields is described. Computational methods were developed for the iterative addition of dipole fields. Various solution potentials were compared inside the boundary with a more precise calculation of the potential to derive optimal schemes for locating the singularities of the dipole fields. Then, the problem of determining solutions to Laplace's equation on an unbounded domain as constrained by pertinent electron trajectory data was considered.

Tilley, D. G.↗

Transfer of line radiation in differentially expanding atmospheres. VI The plane parallel atmosphere with expanding and contracting regions

The non-LTE radiative transfer problem for a two level atom with complete redistribution over a Doppler profile is solved for a plane parallel slab (overlying a radiating photosphere) that has a velocity field which rises symmetrically from zero at either face to a central maximum. Since the velocity gradient reverses, distant layers of the slab become coupled by radiation that jumps intervening layers. The Feautrier method is used, but an iterative variant is also employed as a check in cases where poorly conditioned matrices are encountered. Approximations are developed to explain some of the principal features. It is found that the source function S tends to have two plateaus with values near 2/3 I sub 0 and 1/3 I sub 0, where I sub 0 is the photospheric continuum incident from below; the larger value lies nearer the photosphere. The upper layers sometimes exhibit a rise in S owing to interconnection by radiation to the base. It is noted that the radiation force is largest at the two faces and the midplane. Some line profiles are found to have unusually steep absorptions at rest frequency because of the low excitation in the uppermost, stationary layers.

Noerdlinger, P. D.↗

Laser remote sensing of atmospheric temperature by observing resonant absorption of oxygen

A dual-frequency system is used to demonstrate the measurement of atmospheric temperature, through laser energy absorption monitoring at the center of an O2 resonant absorption line near 770 nm. It is shown that the average temperature of a 1 km path can be determined within 1.0 C, with a noise level of 0.3 C. An iterative algebraic expression for this method was developed, and is shown to be applicable in the troposphere; the effects of pressure and humidity on temperature determination were made clear by the algorithm and found to be small near the earth's surface.

Kalshoven, J. E., Jr.↗

Atmospheric transfer of radiation above an inhomogeneous non-Lambertian reflective ground. I - Theory

A method for solving the three-dimensional equation of transfer for a vertically inhomogeneous atmosphere bounded by a reflecting surface of non-uniform bidirectional reflectance is presented. The technique incorporates a two-dimensional spatial Fourier transform of the transfer equation and solution of the resulting expressions for each Fourier component of the radiation field using the method of Gauss-Seidel iteration. The intensity field in the spatial domain, which is calculated for a variety of altitudes, zenith angles, and azimuths, is reconstructed using the inverse Fourier transform. An empirical surface bidirectional reflectance function is employed, permitting the consideration of non-Lambertian reflective properties.

Diner, D. J.↗

Aerothermal modeling

The objectives, approach, and status of a program to develop the computational fluid dynamics tools needed to improve combustor design and analysis are outlined. The calculation procedure selected consists of a finite difference solution of the time averaged, steady state, primitive variable, elliptic form of the Reynolds equations. Standard TEACH type numerics are used to solve the resulting equations. These include hybrid differencing, SIMPLE algorithm for the pressure field, line by line iterative solution using the ADI method and the tridiagonal matrix algorithm (TDMA). Convergence is facilitated by using under relaxation. The physical processes are modeled by a two equation eddy viscosity model for turbulence; combustion is represented by a simple, irreversible, one step chemical reaction whose rate is influenced only by the time scale of the turbulence. The model evaluation procedure is also described.

Sturgess, G. J.↗

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗

The three-dimensional compressible flow in a radial inflow turbine scroll

This work presents the results of an analytical study and an experimental investigation of the three-dimensional flow in a turbine scroll. The finite element method is used in the iterative numerical solution of the locally linearized governing equations for the three-dimensional velocity potential field. The results of the numerical computations are compared with the experimental measurements in the scroll cross sections, which were obtained using laser Doppler velocimetry and hot wire techniques. The results of the computations show a variation in the flow conditions around the rotor periphery which was found to depend on the scroll geometry.

Hamed, A.↗

Solving Finite-Element Problems on a Concurrent Processor

By use of "conjugate gradients" technique, concurrent efficiency greater than 90 percent. Algorithm applies method of conjugate gradients to iterative solution of finite-element problems on concurrent processor. With algorithm, iteration rates nearly proportional to number of processors. For sufficiently large problems, fraction of proportional speedup achieved, called concurrent efficiency, exceeds 90 percent. Results indicate future application of this and related algorithms to large finite-element problems depend primarily upon applicability of iteratiave techniques, not upon issues of concurrency or efficiency.

Lyzenga, G. A.↗

Star counts and galactic structure

Star counts are discussed as an aid in studies of galactic structure. Iterative computations using a method-of-color-bins to match numerical models with observational data are described. The star counts provide bounds for the luminosity function and the color-magnitude diagram and allow spheroid normalization in a copycat Galaxy model which posits a galaxy with a Population I disk and a Population II spheroid. Good agreement has been obtained between star counts predicted by the model and observed distributions, including color and magnitude ranges and direction. Modifications being considered for the model include consideration of the halo, a galactic bulge and a thick disk. Several galactic structural features which have yet to be modeled and which may be tractable to the copycat model corrected with star counts are identified.

Bahcall, John N.↗

A hybrid numerical technique for predicting the aerodynamic and acoustic fields of advanced turboprops

A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.

Homicz, G. F.↗

An experimental study of an adaptive-wall wind tunnel

A series of adaptive wall ventilated wind tunnel experiments was carried out to demonstrate the feasibility of using the side wall pressure distribution as the flow variable for the assessment of compatibility with free air conditions. Iterative and one step convergence methods were applied using the streamwise velocity component, the side wall pressure distribution and the normal velocity component in order to investigate their relative merits. The advantage of using the side wall pressure as the flow variable is to reduce the data taking time which is one the major contributors to the total testing time. In ventilated adaptive wall wind tunnel testing, side wall pressure measurements require simple instrumentation as opposed to the Laser Doppler Velocimetry used to measure the velocity components. In ventilated adaptive wall tunnel testing, influence coefficients are required to determine the pressure corrections in the plenum compartment. Experiments were carried out to evaluate the influence coefficients from side wall pressure distributions, and from streamwise and normal velocity distributions at two control levels. Velocity measurements were made using a two component Laser Doppler Velocimeter system.

Celik, Zeki↗

Retrieval of constituent mixing ratios from limb thermal emission spectra

An onion-peeling iterative, least-squares relaxation method to retrieve mixing ratio profiles from limb thermal emission spectra is presented. The method has been tested on synthetic data, containing various amounts of added random noise for O3, HNO3, and N2O. The retrieval method is used to obtain O3 and HNO3 mixing ratio profiles from high-resolution thermal emission spectra. Results of the retrievals compare favorably with those obtained previously.

Shaffer, William A.↗

An approach to probabilistic finite element analysis using a mixed-iterative formulation

An efficient algorithm for computing the response sensitivity of finite element problems based on a mixed-iterative formulation is proposed. This method does not involve explicit differentiation of the tangent stiffness array and can be used with formulations for which a consistent tangent stiffness is not readily available. The method has been successfully applied to probabilistic finite element analysis of problems using the proposed mixed formulation, and this exercise has provided valuable insights regarding the extension of the method to a more general class of problems to include material and geometric nonlinearities.

Dias, J. B.↗

Run-time parallelization and scheduling of loops

Run time methods are studied to automatically parallelize and schedule iterations of a do loop in certain cases, where compile-time information is inadequate. The methods presented involve execution time preprocessing of the loop. At compile-time, these methods set up the framework for performing a loop dependency analysis. At run time, wave fronts of concurrently executable loop iterations are identified. Using this wavefront information, loop iterations are reordered for increased parallelism. Symbolic transformation rules are used to produce: inspector procedures that perform execution time preprocessing and executors or transformed versions of source code loop structures. These transformed loop structures carry out the calculations planned in the inspector procedures. Performance results are presented from experiments conducted on the Encore Multimax. These results illustrate that run time reordering of loop indices can have a significant impact on performance. Furthermore, the overheads associated with this type of reordering are amortized when the loop is executed several times with the same dependency structure.

Saltz, Joel H.↗