scientific article; zbMATH DE number 1061627

From MaRDI portal
Publication:4354473





Die Arbeit befaßt sich mit der Herleitung von a posteriori Fehlerabschätzungen mit Hilfe hierarchischer Unterraumzerlegungen für die Randelementmethode (BEM) und für die Kopplung der Finite-Element-Methode (FEM) mit der Randelementmethode. Mit den in den Fehlerabschätungen auftretenden lokalen Fehlerindikatoren werden adaptive Algorithmen formuliert und anhand numerischer Beispiele getestet. Mit Hilfe der abstrakten Theorie der Schwarz-Methoden wird eine a posteriori Abschätzung für den bei der Galerkin-Diskretisierung auftretenden Fehler bewiesen. Im Fall linearer Probleme wird auch der Fehler, der durch inexaktes Lösen der Gleichungssysteme entsteht, in der Abschätzung berücksichtigt. Voraussetzung für die Abschätzung ist neben einer Saturationsannahme die Konstruktion stabiler Zwei-Level-Zerlegungen der beteiligten diskreten Räume. Zum Nachweis der Stabilität wird gezeigt, daß die Konditionszahlen der entsprechenden additiven Schwarz-Operatoren unabhängig von den Diskretisierungsparametern sind oder nur langsam (d.h. logarithmisch bzw. polylogarithmisch) anwachsen. Im Idealfall sollten die Resultate auch für den Fall adaptiver, d.h. lokal verfeinerter, Gitter gültig sein. In der vorliegenden Arbeit wird die Stabilität von Zwei-Level-Zerlegungen für lineare Finite Elemente im \(\mathbb{R}^2\) und für trilineare Finite Elemente im \(\mathbb{R}^3\) gezeigt. Der zweidimensionale Fall ergibt sich als Folgerung aus dem Resultat von \textit{H. Yserentant} [Numer. Math. 49, 379-412 (1986; Zbl 0608.65065)] über hierarchische Basen (HB) im \(\mathbb{R}^2\). Ferner wird polylogarithmisches Wachstum der Konditionszahlen von HB-Operatoren für schwach singuläre und für hypersinguläre Integralgleichungen erster Art im \(\mathbb{R}^2\) und für uniforme Gitter nachgewiesen. Daraus ergibt sich ein Zwei-Level-Resultat für den Fall lokal verfeinerter Gitter. Für Integralgleichungen erster Art im \(\mathbb{R}^3\) wird die Stabilität von Zwei-Level-Zerlegungen für uniforme Gitter nachgewiesen und durch numerische Ergebnisse für allgemeinere Zerlegungen ergänzt. Für den schwach singulären Fall wird ein Ergebnis von \textit{J. Weiße} [Schwarz-Methoden für indefinite und nicht symmetrische Probleme, Dissertation, Inst. Angew. Math., Univ. Hannover (in Vorbereitung)] über Gebietszerlegungen verwendet. Als konkrete Anwendung der abstrakten a posteriori Abschätzung werden die adaptive \(h\)-Version und die adaptive \(p\)-Version des Galerkin-Verfahrens für schwach singuläre Integralgleichungen erster Art auf ebenen und gekrümmten Flächen im \(\mathbb{R}^3\) untersucht und ein adaptiver Algorithmus, der eine Richtungssteuerung der Verfeinerungen erlaubt, formuliert und implementiert. Durch Kombination der Zwei-Level-Resultate für die FEM und die BEM mit den abstrakten Abschätzungen von \textit{R. E. Bank} und \textit{R. K. Smith} [SIAM J. Numer. Anal. 30, No. 4, 921-935 (1993; Zbl 0787.65078)] für indefinite und nichtlineare Probleme wird eine a posteriori Fehlerabschätzung für nichtlineare Übergangsprobleme im \(\mathbb{R}^2\) und im \(\mathbb{R}^3\) hergeleitet. Für die Diskretisierung wird die symmetrische FEM-BEM-Kopplung, die auf \textit{M. Costabel} [The mathematics of finite elements and applications VI, MAFELAP 1987, Proc. 6th Conf., Uxbridge/UK 1987, 281-288 (1988; Zbl 0682.65069)] und \textit{E. P. Stephan} [Comput. Methods Appl. Mech. Eng. 101, No. 1-3, 61-72 (1992; Zbl 0778.73076)] zurückgeht, verwendet. Durch eine geeignete Gewichtung der hierarchischen Fehlerindikatoren für die FEM und die BEM und eine unabhängige Verfeinerung der FEM-Gitter und der BEM-Gitter wird in den numerischen Beispielen eine optimale Verteilung der Freiheitsgrade erreicht.












This page was built for publication:

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4354473)