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 1,189 records · Page 66

A numerical and experimental investigation of electrochemical aircraft deicing

This study was composed of three parts. The first part involved the extension of an existing transient two dimensional numerical code for an electrothermal deicer so that it would simulate the situation where a variable thickness ice layer existed at the outer surface. The Enthalpy Method was used to simulate the phase change, and Gauss-Seidel iteration was used to solve the resulting system of finite difference equations. A set of criteria were developed for determining when a variable thickness ice layer had an effect on deicer performance. The second part was the acquisition and analysis of experimental data. The test model was a section of a Bell UH-1H helicopter blade equipped with an electrothermal deicer. A total of fifty-two thermocouples were utilized to document the thermal response of the blade and decier assembly. In the deicing runs, the experimental temperature response data clearly showed when melting, shedding or refreezing occurred. The tests illustrated that the criterion for shedding in the three cases where it did occur was that the abrasion shield interface temperature was 32 to 34 F. The third part concerned the validation of a one dimensional transient thermal model of an electrothermal deicer by comparison of the predictions with the experimental data. The Enthalpy Method was found to effectively model the phase change which occurred, and the ice shedding algorithm employed in the simulation was also evaluated.

Leffel, K. L.↗

First MEM task IRMEO in IRAF

The first task for image restoration using the Maximum Entropy Method (MEM) in IRAF, called IRMEO, is described. The underlining algorithm is the approximate Newton method for optimization. The basic input images and parameters for deconvolution are described in some detail. Results of preliminary tests, including the number of iterations, required CPU time on a variety of computers, and deconvolved images are reported and compared with those from other deconvolution methods. The merits and limitations of this task are pointed out. The possible development of better MEM tasks on the basis of IRMEO is also discussed.

Wu, Nailong↗

Efficient probabilistic fracture mechanics analysis

A systematic and efficient method for probabilistic fracture mechanics analysis is proposed. The method is based on a most-probable-point-locus concept. The locus is obtained iteratively where the initial locus is determined using the linear approximation of the service life N(X) function about the mean values of the random variables X. Linear and quadratic approximations of N(X) are established locally at the most probable points, and the reliability analysis methods are used to compute the cumulative probabilities. By using two examples, the proposed method is demonstrated to be efficient and accurate. One example involved a random loading and N(X) was computed using cycle-by-cycle integration. The method is general and can be applied to other performance functions. It is particularly suitable when the computation of the performance function is time consuming such that Monte Carlo simulation is prohibitively costly.

Wu, Y.-T.↗

Navier-Stokes solutions using Stetter's method

Compressible Navier-Stokes equations for quasi-one-dimensional flow in a converging-diverging nozzle have been solved using Stetter's three-step predictor-corrector technique. Particular emphasis is given to the minimum iterative step feature for steady-state solutions. It is found that for the nonviscous-dominated case, Stetter's method attains the steady-state solution in the fewest steps when compared to four other currently used techniques.

Graves, R. A., Jr.↗

A look-ahead variant of the Lanczos algorithm and its application to the quasi-minimal residual method for non-Hermitian linear systems

The Lanczos algorithm can be used both for eigenvalue problems and to solve linear systems. However, when applied to non-Hermitian matrices, the classical Lanczos algorithm is susceptible to breakdowns and potential instabilities. In addition, the biconjugate gradient (BCG) algorithm, which is the natural generalization of the conjugate gradient algorithm to non-Hermitian linear systems, has a second source of breakdowns, independent of the Lanczos breakdowns. Here, we present two new results. We propose an implementation of a look-ahead variant of the Lanczos algorithm which overcomes the breakdowns by skipping over those steps where a breakdown or a near-breakdown would occur. The new algorithm can handle look-ahead steps of any length and requires the same number of matrix-vector products and inner products per step as the classical Lanczos algorithm without look-ahead. Based on the proposed look-ahead Lanczos algorithm, we then present a novel BCG-like approach, the quasi-minimal residual (QMR) method, which avoids the second source of breakdowns in the BCG algorithm. We present details of the new method and discuss some of its properties. In particular, we discuss the relationship between QMR and BCG, showing how one can recover the BCG iterates, when they exist, from the QMR iterates. We also present convergence results for QMR, showing the connection between QMR and the generalized minimal residual (GMRES) algorithm, the optimal method in this class of methods. Finally, we give some numerical examples, both for eigenvalue computations and for non-Hermitian linear systems.

Nachtigal, Noel M.↗

Determining the Mass Density Along Magnetic Field Lines from Toroida Eigenfrequencies

Toroidal eigenfrequencies can be used to remotely sense the equatorial mass density rho(sub eq) and the density dependence along a magnetic field line. Here we present improvements to the method of Schulz [1996], which allows rho(sub eq) and the power law index alpha (for mass density along a field line proportional to R(sup -alpha), where R is the radial distance from the center of the Earth) to be determined from the y intercept and slope of a plot of toroidal frequency versus toroidal harmonic number n. Our modifications include a model form for eigenfrequencies with a fractional precision of 0.0005 for -6 less than or = alpha less than or = 6 and 2 less than or = L less than or = 8 (accuracy is doubtful beyond L = 5) and an iterative procedure for getting more accurate results than those found using Schulz's method. In addition, we do an analysis of the effect of random measurement errors. Observed frequencies need to be accurate to approx. 6% (3%) of the fundamental frequency in order to determine rho(sub eq) (alpha) to a precision of 30% (unity). We then apply our method to data generated using the Global Core Plasma Model for plasmaspheric mass density; our analysis demonstrates clearly bow the alpha index represents the mass density dependence on the outer part of the field line (R/(LR(sub E)) greater than or approx. 2/3).

Denton, R.↗

The core-saturation method with partial redistribution

The core-saturation technique for calculating the radiation field in spectral lines (Rybicki, 1971) is extended to cases where partial redistribution is important. The frequency space is discretized and a general vector equation relating the source function S to the mean intensities is derived. The components of S are then used in a lambda iteration to determine the specific intensities. Numerical tests show that the efficiency and accuracy of this method are maintained.

Stenholm, L. G.↗

Propeller-wing interaction using a frequency domain panel method

The unsteady aerodynamic coupling between a propeller and a wing is analyzed using linear compressible aerodynamic theory. The periodic loads are decomposed into harmonics, and the harmonic amplitudes are found iteratively. Each stage of the iteration involves the solution of an isolated propeller or wing problem; the interaction is done through the Fourier transform of the induced velocity field. The method was validated by comparing the predicted velocity field about an isolated propeller with detailed laser Doppler velocimeter measurements and by comparison with mean loads measured in a wing-propeller experiment. Comparisons have also been made between the fluctuating loads predicted by the present method and a quasisteady vortex lattice scheme.

Williams, Marc H.↗

Counter-rotating propellant analysis using a frequency domain panel method

The unsteady aerodynamic coupling between the front and rear rotors in a counter-rotating propeller system is analyzed using a frequency-domain panel method based on linear compressible aerodynamic theory. The periodic loads are decomposed into harmonics, and the harmonic amplitudes are found iteratively. Each stage of the interation involves the solution of an isolated propeller problem, the interaction being done through the Fourier transform of the induced velocity field. The method was validated by comparing mean performance parameters with measured data and by comparing the predicted velocity field with detailed laser Doppler velocimeter measurements. Comparisons have also been made between the fluctuating loads predicted by the present method and a time-domain panel method.

Williams, Marc H.↗

Optimization methods, flux conserving methods for steady state Navier-Stokes equation

Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Adeyeye, John↗

A parabolized Navier-Stokes algorithm for separated supersonic internal flows

A stable global-iteration procedure is applied to a fully implicit noniterative parabolized Navier-Stokes algorithm to improve the solutions of the single-sweep method as well as to resolve small separation bubbles due to shock-wave/boundary-layer interactions. A modified form of the FLARE approximation is employed in the separated regions and some requirements for stable iterations are discussed. The results of the global-iteration scheme for two test cases of supersonic internal flows are compared against the results of the single-sweep method as well as the results of the full Navier-Stokes solver.

Chitsomboon, T.↗

Simplified method for nonlinear structural analysis

A simplified inelastic analysis computer program was developed for predicting the stress-strain history of a thermomechanically cycled structure from an elastic solution. The program uses an iterative and incremental procedure to estimate the plastic strains from the material stress-strain properties and a simulated plasticity hardening model. The simplified method was exercised on a number of problems involving uniaxial and multiaxial loading, isothermal and nonisothermal conditions, and different materials and plasticity models. Good agreement was found between these analytical results and nonlinear finite element solutions for these problems. The simplified analysis program used less than 1 percent of the CPU time required for a nonlinear finite element analysis.

Kaufman, A.↗

Interactive method for computation of viscous flow with recirculation

An interactive method is proposed for the solution of two-dimensional, laminar flow fields with identifiable regions of recirculation, such as the shear-layer-driven cavity flow. The method treats the flow field as composed of two regions, with an appropriate mathematical model adopted for each region. The shear layer is computed by the compressible boundary layer equations, and the slowly recirculating flow by the incompressible Navier-Stokes equations. The flow field is solved iteratively by matching the local solutions in the two regions. For this purpose a new matching method utilizing an overlap between the two computational regions is developed, and shown to be most satisfactory. Matching of the two velocity components, as well as the change in velocity with respect to depth is amply accomplished using the present approach, and the stagnation points corresponding to separation and reattachment of the dividing streamline are computed as part of the interactive solution. The interactive method is applied to the test problem of a shear layer driven cavity. The computational results are used to show the validity and applicability of the present approach.

Brandeis, J.↗

A first-look machine analysis of ERTS-1 data

There are no author-identified significant results in this report. The preliminary analysis of a black and white image of channel 5 data together with tapes from the multispectral scanner is discussed. Four sub-projects based on the initial analysis are discussed. The sub-projects are the analysis by multispectral pattern recognition techniques of the full frame and two particular subframes and a study of the data quality. The procedures for conducting the study and the limitations in analyzing the data are examined. Methods for improving the analysis of the frames by first and second iteration are presented.

Landgrebe, D. A.↗

WIND: Computer program for calculation of three dimensional potential compressible flow about wind turbine rotor blades

A computer program is presented which numerically solves an exact, full potential equation (FPE) for three dimensional, steady, inviscid flow through an isolated wind turbine rotor. The program automatically generates a three dimensional, boundary conforming grid and iteratively solves the FPE while fully accounting for both the rotating cascade and Coriolis effects. The numerical techniques incorporated involve rotated, type dependent finite differencing, a finite volume method, artificial viscosity in conservative form, and a successive line overrelaxation combined with the sequential grid refinement procedure to accelerate the iterative convergence rate. Consequently, the WIND program is capable of accurately analyzing incompressible and compressible flows, including those that are locally transonic and terminated by weak shocks. The program can also be used to analyze the flow around isolated aircraft propellers and helicopter rotors in hover as long as the total relative Mach number of the oncoming flow is subsonic.

Dulikravich, D. S.↗

Numerical simulation of an electrothermal deicer pad

A numerical simulation is developed to investigate the removal of ice from composite aircraft blades by means of electrothermal deicing. The model considers one dimensional, unsteady state heat transfer in the composite blade-ice body. The heat conduction equations are approximated by using the Crank-Nicolson finite difference scheme, and the phase change in the ice layer is handled using the Enthalpy method. To solve the system of equations which result, Gauss-Seidel iteration is used. The simulation computes the temperature profile in the composite blade-ice body, as well as the movement of the ice-water interface, as a function of time. This information can be used to evaluate deicer performance. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Marano, J. J.↗

Predicting Wall Modifications for Adaptive Wind Tunnels

Wall shape changed iteratively until it matches streamlines. FLEXWAL predicts upper and lower wall modifications necessary to remove wall-interference effects in adaptive-wall wind tunnels. FLEXWAL aids in elimination of wall-interference effects on objects tested in typical two-dimensional wind tunnel with rigid sidewalls and flexible, solid floor and ceiling boundaries. Iterative procedure valid for subsonic and transonic test conditions, and convergence of method verified both analytically and experimentally. FLEXWAL written in FORTRAN IV for batch execution.

Everhart, J. L.↗