Multiple right-hand side techniques for the numerical simulation of quasistatic electric and magnetic fields

Authors:
Markus Clemens;Moritz Helias;Thorsten Steinmetz;Georg Wimmer
Affiliations:
Helmut-Schmidt-Universität, Universität der Bundeswehr Hamburg, Theoretische Elektrotechnik und Numerische Feldberechnung, Holstenhofweg 85, D-22043 Hamburg, Germany;Helmut-Schmidt-Universität, Universität der Bundeswehr Hamburg, Theoretische Elektrotechnik und Numerische Feldberechnung, Holstenhofweg 85, D-22043 Hamburg, Germany;Helmut-Schmidt-Universität, Universität der Bundeswehr Hamburg, Theoretische Elektrotechnik und Numerische Feldberechnung, Holstenhofweg 85, D-22043 Hamburg, Germany;Helmut-Schmidt-Universität, Universität der Bundeswehr Hamburg, Theoretische Elektrotechnik und Numerische Feldberechnung, Holstenhofweg 85, D-22043 Hamburg, Germany
Venue:
Journal of Computational and Applied Mathematics
Year:
2008

Citing 4
Cited 0

An efficient preconditioned CG method for the solution of a class of layered problems with extreme contrasts in the coefficients

Journal of Computational Physics
A Deflated Version of the Conjugate Gradient Algorithm

SIAM Journal on Scientific Computing
An Augmented Conjugate Gradient Method for Solving Consecutive Symmetric Positive Definite Linear Systems

SIAM Journal on Matrix Analysis and Applications
Iterative Methods for Sparse Linear Systems

Iterative Methods for Sparse Linear Systems

Quantified Score

Hi-index	7.29

Visualization

Abstract

The simulation of slowly varying transient electric high-voltage fields and magnetic fields requires the repeated and successive solution of high-dimensional linear algebraic systems of equations with identical or near-identical system matrices and different right-hand side vectors. For these solution processes which are required within implicit time integration schemes and nonlinear (quasi-)Newton-Raphson methods an iterative multiple right-hand side (mrhs) scheme is used which recycles vector subspaces resulting from previous preconditioned conjugate gradient iteration runs. The combination of this scheme with a subspace projection extrapolation start value generation scheme is discussed. Numerical results for three-dimensional electric and magnetic field simulations are presented and the efficiency of the new schemes re-using eigenvector information from previous iteration processes with different tolerance criteria are compared to those of standard conjugate gradient iterations.