skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An Implementation and Evaluation of the AMLS Method for SparseEigenvalue Problems

Abstract

We describe an efficient implementation and present aperformance study of an algebraic multilevel sub-structuring (AMLS)method for sparse eigenvalue problems. We assess the time and memoryrequirements associated with the key steps of the algorithm, and compareitwith the shift-and-invert Lanczos algorithm in computational cost. Oureigenvalue problems come from two very different application areas: theaccelerator cavity design and the normal mode vibrational analysis of thepolyethylene particles. We show that the AMLS method, when implementedcarefully, is very competitive with the traditional method in broadapplication areas, especially when large numbers of eigenvalues aresought.

Authors:
; ; ;
Publication Date:
Research Org.:
Ernest Orlando Lawrence Berkeley NationalLaboratory, Berkeley, CA (US)
Sponsoring Org.:
USDOE Director. Office of Science. Advanced ScientificComputing Research
OSTI Identifier:
929006
Report Number(s):
LBNL-57438
R&D Project: KS1210; BnR: KJ0101010; TRN: US0804069
DOE Contract Number:
DE-AC02-05CH11231
Resource Type:
Journal Article
Resource Relation:
Journal Name: ACM Transactions on Mathematical Software; Journal Volume: 34; Journal Issue: 4; Related Information: Journal Publication Date: 09/19/2007
Country of Publication:
United States
Language:
English
Subject:
99; ACCELERATORS; ALGORITHMS; DESIGN; EIGENVALUES; EVALUATION; IMPLEMENTATION; PERFORMANCE; POLYETHYLENES; multilevel substructuring method sparse eigenvalueproblems

Citation Formats

Gao, Weiguo, Li, Xiaoye S., Yang, Chao, and Bai, Zhaojun. An Implementation and Evaluation of the AMLS Method for SparseEigenvalue Problems. United States: N. p., 2006. Web.
Gao, Weiguo, Li, Xiaoye S., Yang, Chao, & Bai, Zhaojun. An Implementation and Evaluation of the AMLS Method for SparseEigenvalue Problems. United States.
Gao, Weiguo, Li, Xiaoye S., Yang, Chao, and Bai, Zhaojun. Tue . "An Implementation and Evaluation of the AMLS Method for SparseEigenvalue Problems". United States. doi:. https://www.osti.gov/servlets/purl/929006.
@article{osti_929006,
title = {An Implementation and Evaluation of the AMLS Method for SparseEigenvalue Problems},
author = {Gao, Weiguo and Li, Xiaoye S. and Yang, Chao and Bai, Zhaojun},
abstractNote = {We describe an efficient implementation and present aperformance study of an algebraic multilevel sub-structuring (AMLS)method for sparse eigenvalue problems. We assess the time and memoryrequirements associated with the key steps of the algorithm, and compareitwith the shift-and-invert Lanczos algorithm in computational cost. Oureigenvalue problems come from two very different application areas: theaccelerator cavity design and the normal mode vibrational analysis of thepolyethylene particles. We show that the AMLS method, when implementedcarefully, is very competitive with the traditional method in broadapplication areas, especially when large numbers of eigenvalues aresought.},
doi = {},
journal = {ACM Transactions on Mathematical Software},
number = 4,
volume = 34,
place = {United States},
year = {Tue Feb 14 00:00:00 EST 2006},
month = {Tue Feb 14 00:00:00 EST 2006}
}
  • The standard implementation of the differential operator (Taylor series) perturbation method for Monte Carlo criticality problems has previously been shown to have a wide range of applicability. In this method, the unperturbed fission distribution is used as a fixed source to estimate the change in the k{sub eff} eigenvalue of a system due to a perturbation. A new method, based on the deterministic perturbation theory assumption that the flux distribution (rather than the fission source distribution) is unchanged after a perturbation, is proposed in this paper. Dubbed the F-A method, the new method is implemented within the framework of themore » standard differential operator method by making tallies only in perturbed fissionable regions and combining the standard differential operator estimate of their perturbations according to the deterministic first-order perturbation formula. The F-A method, developed to extend the range of applicability of the differential operator method rather than as a replacement, was more accurate than the standard implementation for positive and negative density perturbations in a thin shell at the exterior of a computational Godiva model. The F-A method was also more accurate than the standard implementation at estimating reactivity worth profiles of samples with a very small positive reactivity worth (compared to actual measurements) in the Zeus critical assembly, but it was less accurate for a sample with a small negative reactivity worth.« less
  • This work represents a first-of-its-kind successful application to employ advanced numerical methods in solving realistic two-phase flow problems with two-fluid six-equation two-phase flow model. These advanced numerical methods include high-resolution spatial discretization scheme with staggered grids (high-order) fully implicit time integration schemes, and Jacobian-free Newton–Krylov (JFNK) method as the nonlinear solver. The computer code developed in this work has been extensively validated with existing experimental flow boiling data in vertical pipes and rod bundles, which cover wide ranges of experimental conditions, such as pressure, inlet mass flux, wall heat flux and exit void fraction. Additional code-to-code benchmark with the RELAP5-3Dmore » code further verifies the correct code implementation. The combined methods employed in this work exhibit strong robustness in solving two-phase flow problems even when phase appearance (boiling) and realistic discrete flow regimes are considered. Transitional flow regimes used in existing system analysis codes, normally introduced to overcome numerical difficulty, were completely removed in this work. As a result, this in turn provides the possibility to utilize more sophisticated flow regime maps in the future to further improve simulation accuracy.« less
  • Generalizations drawn from implementation research can offer useful guidance in formulating new implementation plants, as the development of a revenue-recycling program illustrates. According to the conventional wisdom drawn from implementation research, the program design should be simple and all the steps necessary to produce the desired results should be identified at the outset. The program that was developed in this case was an extensive one, intended to provide financial relief to all US citizens in an energy emergency. The program would be under the control of one existing agency where equipment, personnel, space, and standard operating procedures were most suitablemore » to the goals of the plan. Although the admonitions derived from past research are straightforward, it turns out that adhering to them in the practical process of resolving specific implementation problems requires imagination and tenacity. Bacause implementation takes place in settings that are extraordinarily diverse, new generalizations regarding implementation will not be easy. 20 references.« less
  • Cited by 1