Iterative Solvers for Modeling Mantle Convection with Strongly Varying Viscosity
Iterative Solvers for Modeling Mantle Convection with Strongly Varying Viscosity
复制标题
用于模拟粘度剧烈变化的地幔对流的迭代求解器
DOI:
--
复制
发表时间:
2011
期刊:
影响因子:
--
通讯作者:
C. Köstler
中科院分区:
文献类型:
--
作者:
C. Köstler
Die Dissertation beschreibt Verbesserungen der FEM-Diskretisierung und des Losers der Stokes-Gleichungen im spharischen Mantelkonvektionsmodell Terra. Zunachst wurde in einem zweidimensionalen quadratischen Gitter mit jeweils stuckweise linearen Ansatzfunktionen fur Druck und Geschwindigkeit eine stabilisierte Diskretisierung nach Dohrmann & Bochev (2004) mit Projektionen auf stuckweise konstante Druckfunktionen implementiert. Deren spektrale Eigenschaften wurden systematisch untersucht. Die Stabilisierung bewirkt eine Gitterunabhangigkeit des Spektrums des Schurkomplements S. Die Viskositatsunabhangigkeit wird durch Prakonditionierung von S mit einer viskositatsabhangigen Massenmatrix Mη bzw. durch Skalierung mit deren Diagonale erreicht. Damit wurden drei Krylov-Unterraumverfahren hinsichtlich ihrer Robustheit gegenuber Viskositatsvariationen und Losertoleranzen untersucht: Druckkorrektur- (PC), Minimierte Residuen- (MINRES) und ein konjugiertes Gradientenverfahren (BPCG) mit einem von Bramble and Pasciak (1988) entwickelten Blockprakonditionierer. PC und BPCG wurden in einer auseren Schleife mit aus Eigenwertabschatzungen berechneten Abbruchkriterien mehrfach gestartet. In der Rechenzeit unterscheiden sich die Loser um weniger als Faktor 2. Bei starken Viskositatskontrasten ist PC der einfachste und schnellste Loser. In Terra kann die o.g. Stabilisierung ohne Einschrankung auf Gittern mit mindestens 85 Millionen Knoten verwendet werden. Fur grobere Gitter wurde eine adaptive Wichtung entwickelt. Das PC-Verfahren in Terra wurde gemas der o.g. Ergebnisse optimiert. Die Diagonalskalierung von S mit Mη bewirkt eine Rechenzeitreduktion um Faktor 4 bei starken lateralen Viskositatsvariationen. Bei Verwendung eines optimalen Multigrid-Losers fur den Impulsoperator ware es Faktor 30. Diese Verbesserungen sind wesentliche Schritte zur Verwendung realitatsnaherer Erdmantelmodelle.
影响因子:
--
作者:
Walzer;R. Hendel
通讯作者:
R. Hendel