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 19 records

A comparison between progressive extension method (PEM) and iterative method (IM) for magnetic field extrapolations in the solar atmosphere

This paper presents a comparison between two numerical methods for the extrapolation of nonlinear force-free magnetic fields, viz the Iterative Method (IM) and the Progressive Extension Method (PEM). The advantages and disadvantages of these two methods are summarized, and the accuracy and numerical instability are discussed. On the basis of this investigation, it is claimed that the two methods do resemble each other qualitatively.

Wu, S. T.↗

Comparative Analysis of Soft-Error Detection Strategies: A Case Study with Iterative Methods

Undetected soft errors caused by transient bit flips can lead to silent data corruption (SDC), an undesirable outcome where invalid results pass for valid ones. This has motivated the design of soft error detectors to minimize SDCs. However, the detectors have been studied under different contexts, making comparative evaluation difficult. In this paper, we present the first comprehensive evaluation of four online soft error detection techniques in detecting the adverse impact of soft errors on iterative methods. We observe that, across five iterative methods, the detectors studied achieve high but not perfect detection rates. To understand the potential for improved detection, we evaluate a machine-learning based detector that takes as features that are the runtime features observed by the individual detectors to arrive at their conclusions. Our evaluation demonstrates improved but still far from perfect detection accuracy for the machine learning based detectors. This extensive evaluation demonstrates the need for designing error detectors to handle the evolutionary behavior exhibited by iterative solvers.

Silent data corruption, Soft error detection, Iter↗

Analysis of Monte Carlo accelerated iterative methods for sparse linear systems: Analysis of Monte Carlo accelerated iterative methods for sparse linear systems

Here, we consider hybrid deterministic-stochastic iterative algorithms for the solution of large, sparse linear systems. Starting from a convergent splitting of the coefficient matrix, we analyze various types of Monte Carlo acceleration schemes applied to the original preconditioned Richardson (stationary) iteration. We expect that these methods will have considerable potential for resiliency to faults when implemented on massively parallel machines. We also establish sufficient conditions for the convergence of the hybrid schemes, and we investigate different types of preconditioners including sparse approximate inverses. Numerical experiments on linear systems arising from the discretization of partial differential equations are presented.

97 MATHEMATICS AND COMPUTING↗

The solution of radiative transfer problems in molecular bands without the LTE assumption by accelerated lambda iteration methods

An iterative method based on the use of approximate transfer operators, which was designed initially to solve multilevel NLTE line formation problems in stellar atmospheres, is adapted and applied to the solution of the NLTE molecular band radiative transfer in planetary atmospheres. The matrices to be constructed and inverted are much smaller than those used in the traditional Curtis matrix technique, which makes possible the treatment of more realistic problems using relatively small computers. This technique converges much more rapidly than straightforward iteration between the transfer equation and the equations of statistical equilibrium. A test application of this new technique to the solution of NLTE radiative transfer problems for optically thick and thin bands (the 4.3 micron CO2 band in the Venusian atmosphere and the 4.7 and 2.3 micron CO bands in the earth's atmosphere) is described.

Kutepov, A. A.↗

Annual Copper Mountain Conferences on Multigrid and Iterative Methods, Copper Mountain, Colorado

This project supported the Copper Mountain Conference on Multigrid and Iterative Methods, held from 2007 to 2015, at Copper Mountain, Colorado. The subject of the Copper Mountain Conference Series alternated between Multigrid Methods in odd-numbered years and Iterative Methods in even-numbered years. Begun in 1983, the Series represents an important forum for the exchange of ideas in these two closely related fields. This report describes the Copper Mountain Conference on Multigrid and Iterative Methods, 2007-2015. Information on the conference series is available at http://grandmaster.colorado.edu/~copper/.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An Exploratory Study of a Subspace Iteration Method as an Alternative to the QR Method for Floquet Eigenanalysis

Floquet eigenanalysis requires a few dominant eigenvalues of the Floquet transition matrix (FTM). Although the QR method is used almost exclusively, it is expensive for such partial eigenanalysis; the operation counts and, thereby, the approximate machine-time grow cubically with the matrix order. Accordingly, for Floquet eigenanalysis, the Arnold-Saad method, a subspace iteration method, is investigated as an alternative to the QR method. The two methods are compared for machine-time efficiency and the residual errors of the corresponding eigenpairs. The Arnolds-Saad method takes much less machine-time than the QR method with comparable computational reliability and offers promise fpr large-scale Floquet eigenanalysis.

Achar, N. S.↗

On semi-convergence and inexact iteration of the GSS iteration method for nonsymmetric singular saddle point problems

In this paper, the generalized shift-splitting (GSS) iteration method is used to solve nonsymmetric singular saddle point systems with the positive real (1, 1)-block sub-matrix. By a careful analysis of the spectral properties of the corresponding iteration matrix, we investigate the semi-convergence property of the GSS iteration method, give the sharper bounds of the eigenvalues for the GSS iteration method and propose the inexact GSS iteration method, also. Numerical experiments are given to show the correctness of the theoretical analysis and the feasibility of the GSS preconditioner with appropriate parameters.

97 MATHEMATICS AND COMPUTING↗

Nonlocal transport in bounded two-dimensional systems: An iterative method

Here, the concept of transport mediated through the dynamics of 'jumping' particles is used to develop an iterative method for obtaining steady-state solutions to the nonlocal transport equation in two dimensions (2D). The technique is self-adjoint and capable of correctly treating spatially non-uniform, asymmetric systems. An appropriate reduced version of the iteration method is used to compare with results obtained with a self-adjoint 1D transport matrix approach. The transport 'jump' probability distribution functions (pdfs) are based on Lévy α-stable distributions. The technique can handle the entire Lévy α-parameter range from one (Lorentz distributions) to two (Gaussian distributions). Cases with α = 2 (standard diffusion) are used to establish the validity of the iterative method. The capabilities of the iterative method are demonstrated by presenting examples from systems with various source configurations, boundary shapes, boundary conditions and spatial variations in parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Leapfrog variants of iterative methods for linear algebra equations

Two iterative methods are considered, Richardson's method and a general second order method. For both methods, a variant of the method is derived for which only even numbered iterates are computed. The variant is called a leapfrog method. Comparisons between the conventional form of the methods and the leapfrog form are made under the assumption that the number of unknowns is large. In the case of Richardson's method, it is possible to express the final iterate in terms of only the initial approximation, a variant of the iteration called the grand-leap method. In the case of the grand-leap variant, a set of parameters is required. An algorithm is presented to compute these parameters that is related to algorithms to compute the weights and abscissas for Gaussian quadrature. General algorithms to implement the leapfrog and grand-leap methods are presented. Algorithms for the important special case of the Chebyshev method are also given.

Saylor, Paul E.↗

Iterative methods based upon residual averaging

Iterative methods for solving boundary value problems for systems of nonlinear partial differential equations are discussed. The methods involve subtracting an average of residuals from one approximation in order to arrive at a subsequent approximation. Two abstract methods in Hilbert space are given and application of these methods to quasilinear systems to give numerical schemes for such problems is demonstrated. Potential theoretic matters related to the iteration schemes are discussed.

Neuberger, J. W.↗

Factorizations of hyperpower family of iterative methods via least squares approach

A new method for finding proper factorizations of the hyperpower family of iterative methods for computing generalized inverses is proposed. The method is based on numerical optimization of an appropriate least squares problem. An iterative method in a non-normalized form for finding the generalized outer inverses with prescribed range and null space of any given complex matrix is proposed and considered. Computational complexity and convergence of the method are studied. Several numerical experiments are also presented.

97 MATHEMATICS AND COMPUTING↗

Ground-Truth Prediction to Accelerate Soft-Error Impact Analysis for Iterative Methods

Understanding the impact of soft errors on applications can be expensive. Often, it requires an extensive error injection campaign involving numerous runs of the full application in the presence of errors. In this paper, we present a novel approach to arriving at the ground truth--the true impact of an error on the final output--for iterative methods by observing a small number of iterations to learn deviations between normal and error-impacted execution. We develop a machine learning based predictor for three iterative methods to generate ground-truth results without running them to completion for every error injected. We demonstrate that this approach achieves greater accuracy than alternative prediction strategies, including three existing soft error detection strategies. We demonstrate the effectiveness of the ground truth prediction model in evaluating vulnerability and the effectiveness of soft error detection strategies in the context of iterative methods.

iterative methods, failure prediction, Soft error ↗

Numerical Comparison of Periodic MoM (Method of Moments) and BMIA (Banded Matrix Iteration Method)

The most popular numerical technique in rough surface scattering is the Method of Moments (MoM). Since the scattering patch size is finite, the edge current must be suppressed to obtain accurate scattering cross sections. Two standard ways to minimize the edge current are periodic boundary conditions and incident wave tapering. We compare the accuracy & computational requirements of these methods.

Method of Moments Banded Matrix Iteration Method M↗

An iterative method for computing a symplectic SVD-like decomposition

In this paper, we present an effective iterative method for computing symplectic SVD-like decomposition for a 2n-by-m rectangular real matrix A. The main purpose here is a block-power iterative method with the ortho-symplectic SR decomposition in the normalization step. We compute an k-block symplectic SVD-like decomposition, namely A{sub k}=S{sub k}Σ{sub k}V{sub k}{sup T} where S{sub k}∈ℝ{sup 2n×2k} is symplectic and V{sub k}∈ℝ{sup m×2k} is orthogonal. For large matrices, we usually have to compute only part of eigenvalues. The main interest of the proposed method is to efficiently compute the desired number of ordered in magnitude eigenvalues of structured matrices.

97 MATHEMATICS AND COMPUTING↗