Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

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
DOE OSTI2015

Anderson acceleration of the Jacobi iterative method: An efficient alternative to Krylov methods for large, sparse linear systems

We employ Anderson extrapolation to accelerate the classical Jacobi iterative method for large, sparse linear systems. Specifically, we utilize extrapolation at periodic intervals within the Jacobi iteration to develop the Alternating Anderson–Jacobi (AAJ) method. We verify the accuracy and efficacy of AAJ in a range of test cases, including nonsymmetric systems of equations. We demonstrate that AAJ possesses a favorable scaling with system size that is accompanied by a small prefactor, even in the absence of a preconditioner. In particular, we show that AAJ is able to accelerate the classical Jacobi iteration by over four orders of magnitude, with speed-ups that increase as the system gets larger. Moreover, we find that AAJ significantly outperforms the Generalized Minimal Residual (GMRES) method in the range of problems considered here, with the relative performance again improving with size of the system. As a result, the proposed method represents a simple yet efficient technique that is particularly attractive for large-scale parallel solutions of linear systems of equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

Subspace Iteration Method for Complex Eigenvalue Problems with Nonsymmetric Matrices in Aeroelastic System

Modern airplane design is a multidisciplinary task which combines several disciplines such as structures, aerodynamics, flight controls, and sometimes heat transfer. Historically, analytical and experimental investigations concerning the interaction of the elastic airframe with aerodynamic and in retia loads have been conducted during the design phase to determine the existence of aeroelastic instabilities, so called flutter .With the advent and increased usage of flight control systems, there is also a likelihood of instabilities caused by the interaction of the flight control system and the aeroelastic response of the airplane, known as aeroservoelastic instabilities. An in -house code MPASES (Ref. 1), modified from PASES (Ref. 2), is a general purpose digital computer program for the analysis of the closed-loop stability problem. This program used subroutines given in the International Mathematical and Statistical Library (IMSL) (Ref. 3) to compute all of the real and/or complex conjugate pairs of eigenvalues of the Hessenberg matrix. For high fidelity configuration, these aeroelastic system matrices are large and compute all eigenvalues will be time consuming. A subspace iteration method (Ref. 4) for complex eigenvalues problems with nonsymmetric matrices has been formulated and incorporated into the modified program for aeroservoelastic stability (MPASES code). Subspace iteration method only solve for the lowest p eigenvalues and corresponding eigenvectors for aeroelastic and aeroservoelastic analysis. In general, the selection of p is ranging from 10 for wing flutter analysis to 50 for an entire aircraft flutter analysis. The application of this newly incorporated code is an experiment known as the Aerostructures Test Wing (ATW) which was designed by the National Aeronautic and Space Administration (NASA) Dryden Flight Research Center, Edwards, California to research aeroelastic instabilities. Specifically, this experiment was used to study an instability known as flutter. ATW was a small-scale airplane wing comprised of an airfoil and wing tip boom. This wing was formulated based on a NACA-65A004 airfoil shape with a 3.28 aspect ratio. The wing had a span of 18 inch with root chord length of 13.2 inch and tip chord length of 8.7 inch. The total area of this wing was 197 square inch. The wing tip boom was a 1 inch diameter hollow tube of length 21.5 inch. The total weight of the wing was 2.66 lbs.

Pak, Chan-gi↗
DOE OSTI2018

Parallel iterative methods for a finite family of sequences of nearly nonexpansive mappings in Hilbert spaces

In this paper, we propose some new parallel iterative methods by the hybrid and shrinking projection to find a common fixed point of a finite family of sequences of nearly nonexpansive mappings in Hilbert spaces when the domain C satisfies diam(C) < ∞. We also give some applications of our main results for the problem of finding a common fixed point of nonexpansive mappings, nonexpansive semi-groups, the problem of finding a common zero point of monotone operators, the system of generalized mixed equilibrium problems and the system of variational inequalities without the condition diam(C) < ∞. Three numerical examples also are given to illustrate the effectiveness of the proposed algorithms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Applications of partial orderings to the study of positive definiteness, monotonicity, and convergence of iterative methods for linear systems.

Most of the known results concerning convergence of iterative methods for solving linear systems involve either positive definiteness or monotonicity. In this paper a new concept, called K-semipositivity, is introduced, which provides a link between convergence theory, monotonicity, and positive definiteness. By using this concept, together with partial orderings in Euclidean n-space, several new convergence theorems are proved. Application to Jacobi's methods and the theory of regular splittings shows the usefulness of these new results.

Vandergraft, J. S.↗
NASA NTRS1996

Practical Aspects of Krylov Subspace Iterative Methods in CFD

Implementation issues associated with the application of Krylov subspace iterative methods, such as Newton-GMRES, are presented within the framework of practical computational fluid dynamic (CFD) applications. This paper categorizes, evaluates, and contrasts the major ingredients (function evaluations, matrix-vector products, and preconditioners) of Newton-GMRES Krylov subspace methods in terms of their effect on the local linear and global nonlinear convergence, memory requirements, and accuracy. The discussion focuses on Newton-GMRES in both a structured multi-zone incompressible Navier-Stokes solver and an unstructured mesh finite-volume Navier-Stokes solver. Approximate versus exact matrix-vector products, effective preconditioners, and other pertinent issues are addressed.

Pulliam, Thomas H.↗
NASA NTRS1990

A two-step iterative method for evolving nonlinear acoustic systems to a steady-state

A new approach for evolving two-dimensional nonlinear acoustic systems with flow to a steady state is presented. The approach is a two-step iterative method which is tested on a benchmark acoustic problem for which an exact analytical solution is available. Results are also calculated for a nonlinear acoustic problem for which an exact analytical solution is not known. Results indicate that the two-step method represents a powerful, efficient, and stable method for evolving two-dimensional acoustic systems to a steady state, and that the method is applicable to any number of spatial dimensions and to other hyperbolic systems. It is noted that for the benchmark problem only a single iteration on the method is required when the transient and steady-state field are of the same order of magnitude; however, four iterations are required when the steady-state field is several orders of magnitude smaller than the transient field. This method requires six iterations before achieving a steady state for the nonlinear test problem.

Watson, Willie R.↗
NASA NTRS1989

Iterative methods for design sensitivity analysis

A numerical method is presented for design sensitivity analysis, using an iterative-method reanalysis of the structure generated by a small perturbation in the design variable; a forward-difference scheme is then employed to obtain the approximate sensitivity. Algorithms are developed for displacement and stress sensitivity, as well as for eignevalues and eigenvector sensitivity, and the iterative schemes are modified so that the coefficient matrices are constant and therefore decomposed only once.

Belegundu, A. D.↗
DOE OSTI2018

A modified iterative method for split problem of variational inclusions and fixed point problems

In this paper, we study a new split problem of variational inclusions and fixed point problems. With the help of averaged mappings, we introduce a new iterative method to solve the stated split problem. Strong convergence results are obtained under mild conditions. Also, we present some preliminary numerical examples to demonstrate the convergence and efficiency of our proposed algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

A fast-converging iterative method for X-ray in-line phase contrast tomography

X-ray in-line phase contrast tomography holds great promise for the quantitative analysis of soft materials. However, its applications have been limited, so far, by the fact that direct methods based on the transport-of-intensity equation and the contrast transfer function are sensitive to noise and applicable only to limited types of samples. Here, we propose an iterative method based on the Gerchberg-Saxton algorithm (R. W. Gerchberg and W. O. Saxton, Optik 35, 237 (1972)), but overcoming its slow convergence by an acceleration technique, named random signed feedback, which shows an excellent performance, both in numerical simulation and tomographic experiment, of discriminating various polymers even when using 53 keV synchrotron X-rays.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Efficient iterative method for solving the Dirac-Kohn-Sham density functional theory

We present for the first time an efficient iterative method to directly solve the four-component Dirac-Kohn-Sham (DKS) density functional theory. Due to the existence of the negative energy continuum in the DKS operator, the existing iterative techniques for solving the Kohn-Sham systems cannot be efficiently applied to solve the DKS systems. The key component of our method is a novel filtering step (F) which acts as a preconditioner in the framework of the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The resulting method, dubbed the LOBPCG-F method, is able to compute the desired eigenvalues and eigenvectors in the positive energy band without computing any state in the negative energy band. The LOBPCG-F method introduces mild extra cost compared to the standard LOBPCG method and can be easily implemented. We demonstrate our method in the pseudopotential framework with a planewave basis set which naturally satisfies the kinetic balance prescription. Numerical results for Pt$_{2}$, Au$_{2}$, TlF, and Bi$_{2}$Se$_{3}$ indicate that the LOBPCG-F method is a robust and efficient method for investigating the relativistic effect in systems containing heavy elements.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Iterative Method to Obtain Raman Spectra of Organic Dyes by Their Chaotic Generation in Multiple-Scattering Media

An iterative method is developed to obtain vibrational spectra of laser organic dyes from their secondary radiation spectra in multiple-scattering media, which manifests as a result of coupled process of chaotic generation (CG) and stimulated Raman scattering (SRS). The developed method increases significantly the capability to detect weak Raman lines by using batch processing of a set of experimental SRS-CG spectra, which makes it possible to optimize the deconvolution of these spectra and to average vibrational spectra obtained from the set of experimental SRS-CG spectra. The method is tested on spectra of the SRS-CG laser dye pyrromethene 597 in a vesicular film. The obtained results agree well with Raman spectra calculated by the Hartree–Fock method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗