A Moment-Parity Multigrid Preconditioner for the First-Order System Least-Squares Formulation of the Boltzmann Transport Equation

This paper describes a preconditioned conjugate gradient scheme for the Pn-h finite element discretization of a first-order system least-squares (FOSLS) formulation of the Boltzmann transport equation. The preconditioner is based on the norm equivalence between the FOSLS functional and a V norm. Its...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:SIAM journal on scientific computing Jg. 25; H. 2; S. 513 - 533
Hauptverfasser: Brown, P. N., Lee, B., Manteuffel, T. A.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Philadelphia Society for Industrial and Applied Mathematics 01.01.2003
Schlagworte:
Algorithms
Applied mathematics
Atoms & subatomic particles
Laboratories
Parity
ISSN:1064-8275, 1095-7197
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:This paper describes a preconditioned conjugate gradient scheme for the Pn-h finite element discretization of a first-order system least-squares (FOSLS) formulation of the Boltzmann transport equation. The preconditioner is based on the norm equivalence between the FOSLS functional and a V norm. Its realization is an inexact inversion of the system of partial differential equations corresponding to this V norm. This preconditioner is essentially a moment-parity multigrid solver. It involves a 2 × 2 block diagonal system with each block describing only the interior like-parity coupling (even-even or odd-odd) of the spherical harmonic coefficients or moments. The interior cross-parity coupling is handled in the outer conjugate gradient iteration. Since the like-parity coupling consists of only second- and zeroth-order differential terms, whereas the full Pn system consists also of cross-parity, first-order coupling terms, the construction of a robust multigrid algorithm for each diagonal block is easier than the construction for the full Pn system. Numerical results indicate that this preconditioned conjugate gradient algorithm is more robust than a stand-alone multigrid solver for the full Pn system.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
content type line 14
ISSN:1064-8275
1095-7197
DOI:10.1137/S1064827502407172