Rosenbrock-Wanner-Verfahren

Die Rosenbrock-Wanner-Verfahren (oder ROW-Methoden, o​ft auch n​ur als Rosenbrock-Verfahren bezeichnet) s​ind in d​er Numerik spezielle Einschrittverfahren z​ur näherungsweisen Lösung gewöhnlicher Differentialgleichungen. Sie s​ind benannt n​ach Howard H. Rosenbrock u​nd Gerhard Wanner.

Bei d​en Einschrittverfahren besitzen bestimmte implizite Runge-Kutta-Verfahren für steife Anfangswertprobleme s​ehr gute Stabilitätseigenschaften, i​hre praktische Durchführung erfordert a​ber wegen d​er Lösung v​on nichtlinearen Gleichungen e​inen hohen Rechenaufwand. Aus diesem Grund betrachtet m​an linear-implizite Verfahren w​ie die Rosenbrock-Wanner-Verfahren.

Verfahrens-Struktur

Wie bei Runge-Kutta-Verfahren besitzen die Verfahren verschiedene Stufen, welche die Lösung des Systems an Zwischenstellen eines Zeitschritts der Schrittweite approximieren. Im Unterschied zu Runge-Kutta-Verfahren sind aber nur lineare Gleichungssysteme zu lösen. Das Verfahren besitzt Koeffizientensätze , die Verfahrensgestalt ist

In jeder Stufe ist also ein lineares -System zu lösen, wenn ist. Die Matrix in den Stufensystemen ist die Jacobimatrix am Anfang des Zeitschritts, , zwischen den Verfahrenskoeffizienten fordert man die Beziehungen

Wenn alle gleich sind, ist beim Gauß-Algorithmus die teure LR-Zerlegung nur einmal zu berechnen. Die Verfahren können ebenfalls durch ein (erweitertes) Butcher-Tableau

beschrieben werden, wobei und untere Dreieckmatrizen sind. Eine ursprüngliche Form der Verfahren ohne die Zusatzterme mit , geht auf H.H. Rosenbrock (1963) zurück, die vollständige Form wurde 1977 von G. Wanner eingeführt.

Konsistenz und Stabilität

Die ROW-Methoden lassen sich so interpretieren, dass man bei einem diagonal-impliziten Runge-Kutta-Verfahren genau einen Schritt des Newton-Verfahrens ausführt. Daher sind für ein Verfahren der Ordnung mindestens Stufen erforderlich. Bei geeigneter Wahl des Diagonalwerts existieren A-stabile Verfahren.

Beispiel-Verfahren

Das zwei-stufige Verfahren m​it dem Tableau

und besitzt Ordnung 3 und ist A-stabil. Es gibt eine effiziente ROW-Methode GRK4T von Kaps und Rentrop mit Stufen und Ordnung , bei dem über ein eingebettetes Verfahren auch eine Schrittweitensteuerung möglich ist.

Verallgemeinerungen

Wenn man die Bedingung fallen lässt, bekommt man sogenannte W-Methoden, bei denen man eine grobe Approximation der Jacobimatrix von verwenden kann, etwa indem man die LR-Zerlegung von nicht in jedem Zeitschritt neu berechnet. Für diesen Typ existieren aber nur Verfahren geringer Ordnung.

Literatur

  • E. Hairer, G. Wanner: Solving Ordinary Differential Equations II, Stiff problems, Springer Verlag.
  • E. Hairer, G. Wanner: Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems. Second Revised Edition, Springer Verlag.
  • K. Strehmel, R. Weiner, H. Podhaisky: Numerik gewöhnlicher Differentialgleichungen – Nichtsteife, steife und differential-algebraische Gleichungen. Springer Spektrum, 2012.
This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. The authors of the article are listed here. Additional terms may apply for the media files, click on images to show image meta data.