Abstract
We outline a mathematical formulation for mantle convection which can deal with the viscoelastic-plastic rheology of the cool parts of the lithosphere. This formulation is then analyzed to expose the numerical challenges inherent in the equations and a suitable solution strategy is outlined. With this strategy in place, we discuss a parallel implementation, explaining how we maintain computational efficiency, in a tiered and modular software environment. We show two examples from geodynamics which demonstrate the capability of the numerical method and its implementation.
Original language | English |
---|---|
Pages (from-to) | 69-82 |
Number of pages | 14 |
Journal | Physics of the Earth and Planetary Interiors |
Volume | 163 |
Issue number | 1-4 |
DOIs | |
Publication status | Published - 15 Aug 2007 |
Externally published | Yes |