 |
| Titel |
Implementation of Newton-Rapshon iterations for parallel staggered-grid geodynamic models |
| VerfasserIn |
A. A. Popov, B. J. P. Kaus |
| Konferenz |
EGU General Assembly 2012
|
| Medientyp |
Artikel
|
| Sprache |
Englisch
|
| Digitales Dokument |
PDF |
| Erschienen |
In: GRA - Volume 14 (2012) |
| Datensatznummer |
250063926
|
|
|
|
|
|
| Zusammenfassung |
| Staggered-grid finite differences discretization has a good potential for solving highly heterogeneous geodynamic models on parallel computers (e.g. Tackey, 2008; Gerya &Yuen, 2007). They are inherently stable, computationally inexpensive and relatively easy to implement. However, currently used staggered-grid geodynamic codes employ almost exclusively the sub-optimal Picard linearization scheme to deal with nonlinearities. It was shown that Newton-Rapshon linearization can lead to substantial improvements of the solution quality in geodynamic problems, simultaneously with reduction of computer time (e.g. Popov & Sobolev, 2008).
This work is aimed at implementation of the Newton-Rapshon linearization in the parallel geodynamic code LaMEM together with staggered-grid discretization and viso-(elasto)-plastic rock rheologies. We present the expressions for the approximate Jacobian matrix, and give detailed comparisons with the currently employed Picard linearization scheme, in terms of solution quality and number of iterations.
References:
Gerya, T.V., Yuen, D.A., 2007. Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems. Phys. Earth Planet. Inter., 163: 83-105.
Tackley, P., 2008. Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yinyang grid, Phys. Earth Planet. Inter., 171: 7–18.
Popov, A.A., Sobolev, S.V., 2008. SLIM3D: A tool for three-dimensional thermomechanical modeling of lithospheric deformation with elasto-visco-plastic rheology. Phys. Earth Planet. Inter., doi:10.1016/j.pepi.2008.03.007 |
|
|
|
|
|
|