On a numerical Lyapunov-Schmidt method for operator equations (Q1313248)

Es wird die Konvergenz von Diskretisierungsverfahren zur numerischen Berechnung von Lösungsmannigfaltigkeiten für Gleichungen der Form \(G(z, \lambda) = 0\) untersucht, wobei \(G:D(G) \subset E_ 0 \times \mathbb{R}^ q \to \widehat E\) abbildet und \(E_ 0\), \(\widehat E_ 0\) Hilbertsche Räume sind. Die zu berechnenden Mannigfaltigkeiten sollen von einem höheren singulären Punkt \(x_ 0: = (z_ 0, \lambda_ 0) \in E_ 0 \times \mathbb{R}^ q\) abzweigen, d.h. für den Nullraum von \(G'(x_ 0)\) soll \(\dim N(G'(x_ 0)) = m + q > q\) gelten, wobei \(m>1\) sein darf, und für das numerische Verfahren wird vorausgesetzt, daß Approximationen für \(N(G'(x_ 0))\) und \(N(G'(x_ 0)^*)\) bekannt sind. Die in der Arbeit behandelte Methode geht von der Darstellung \(x=x_ 0 = tv + t^ rw\) mit \(r \in \mathbb{N}\), \(r \geq 2\), und \(v \in N(G'(x_ 0))\), \(w \in N (G'(x_ 0))^ \perp\) für die abzweigende Mannigfaltigkeit aus. Dabei sind \(v = v(t,w)\) und \(w = w(t,v)\) durch zwei nichtlineare Gleichungen bestimmt, die sich durch Anwendung einer gewissen Version der Lyapunov-Schmidt-Methode ergeben. Es werden nun diese beiden Gleichungen mit einer gewissen Konsistenzordnung \(O(h^ p)\) diskretisiert, wodurch man Approximationen der Gestalt \((x-x_ 0)^ h = tv^ h + t^ rw^ h\) erhält, die detaillierteren Aufschluß über die Verzweigungsstruktur geben, als dies eine alleinige Approximation an \(x- x_ 0\) tun würde. Hierin liegt ein Vorteil der verwendeten Methode. Unter passenden Voraussetzungen an die diskretisierenden Operatoren, u.a. eine inverse Stabilität von \(G_ h'\), wird Konvergenz der Ordnung \(O(h^ p)\) von \(v^ h\) und \(w^ h\) gegen \(v\) bzw. \(w\) bewiesen. Als konkretes Beispiel einschließlich numerischer Berechnungen behandelt der Verf. die Aufgabe \(\Delta z + \lambda g(z) = 0\), \(\lambda \in \mathbb{R}\), \(z \in [0,1]^ 2\) unter Dirichletbedingungen, die mit dem üblichen Differenzenverfahren diskretisiert wird.