Beschreibung der Forschungsarbeit:
Forschungsgegenstand sind parallelisierbare numerische Verfahren zur Lösung von Anfangswertproblemen großer Systeme von Algebro-Differentialgleichungen (DAE)
, , gegeben, gesucht, wie sie etwa bei der Prozeßsimulation chemischer Anlagen oder in der Schaltkreissimulation auftreten.
Solche Systeme sind oftmals schon durch die Modellierung entsprechend ihrer funktionellen Bestandteile in Teilsysteme strukturiert. Die von uns betrachteten Block-Waveform-Iterationsverfahren [1] weisen günstige Voraussetzungen für eine Parallelisierung auf der Ebene der DAE-Systeme auf. Sie sind jedoch nur für solche Aufgabenklassen anwendbar, für die eine für diese Verfahren geeignete Blockzerlegung bestimmt werden kann. Für DAE-Systeme, bei denen aufgrund starker Kopplungen zwischen den Teilsystemen keine solche Blockzerlegungen existieren, betrachten wir parallelisierbare strukturierte Newton-Verfahren. Hierbei werden die DAE-Systeme nach ihrer Zeitdiskretisierung in geeignet erweiterte nichtlineare Gleichungssysteme mit spezieller Blockstruktur überführt [2] und nach Linearisierung über die parallele Behandlung linearer Teilgleichungssysteme gelöst.
Bei einstufig strukturierten Verfahren werden die inneren Variablen und die äußeren Variablen der Teilsysteme bestimmt und das um Identifikationsgleichungen für äußere Variablen erweiterte System
erzeugt. Wir partitionieren dann so, daß regulär ist. Dabei erfolgt die Auswahl der Gleichungen mit Hilfe eines Algorithmus zur Pivotwahl in den überbestimmten Teilmatrizen . Die Bestimmung der Inkremente für einen Schritt des Newtonverfahrens erfolgt in 3 Teilschritten:
Einstufig strukturierte Newtonverfahren wurden von uns implementiert und an praxisrelevanten Beispielen getestet. Sie haben sich dabei als moderat parallele Verfahren erwiesen.
Projektliteratur: