dorsal/arxiv
View SchemaMultigrid iterative algorithm using pseudo-compressibility for three-dimensional mantle convection with strongly variable viscosity
| Authors | Masanori Kameyama, Akira Kageyama, Tetsuya Sato |
|---|---|
| Categories | |
| ArXiv ID | physics/0410249 |
| URL | https://arxiv.org/abs/physics/0410249 |
| DOI | 10.1016/j.jcp.2004.11.030 |
Abstract
A numerical algorithm for solving mantle convection problems with strongly variable viscosity is presented. Equations for conservation of mass and momentum for highly viscous and incompressible fluids are solved iteratively by a multigrid method in combination with pseudo-compressibility and local time stepping techniques. This algorithm is suitable for large-scale three-dimensional numerical simulations, because (i) memory storage for any additional matrix is not required and (ii) vectorization and parallelization are straightforward. The present algorithm has been incorporated into a mantle convection simulation program based on the finite-volume discretization in a three-dimensional rectangular domain. Benchmark comparisons with previous two- and three-dimensional calculations including the temperature- and/or depth-dependent viscosity revealed that accurate results are successfully reproduced even for the cases with viscosity variations of several orders of magnitude. The robustness of the numerical method against viscosity variation can be significantly improved by increasing the pre- and post-smoothing calculations during the multigrid operations, and the convergence can be achieved for the global viscosity variations up to 1e10.
{
"annotation_id": "4a9eac57-6b96-40cf-82ed-00118e1b4e42",
"date_created": "2026-03-02T18:00:53.450000Z",
"date_modified": "2026-03-02T18:00:53.450000Z",
"file_hash": "8045ef1df15eeb4dd70aaf3ec7bb82f8b403ee5a4eb60cf10f010453580a47ed",
"private": false,
"record": {
"abstract": "A numerical algorithm for solving mantle convection problems with strongly\nvariable viscosity is presented. Equations for conservation of mass and\nmomentum for highly viscous and incompressible fluids are solved iteratively by\na multigrid method in combination with pseudo-compressibility and local time\nstepping techniques. This algorithm is suitable for large-scale\nthree-dimensional numerical simulations, because (i) memory storage for any\nadditional matrix is not required and (ii) vectorization and parallelization\nare straightforward. The present algorithm has been incorporated into a mantle\nconvection simulation program based on the finite-volume discretization in a\nthree-dimensional rectangular domain. Benchmark comparisons with previous two-\nand three-dimensional calculations including the temperature- and/or\ndepth-dependent viscosity revealed that accurate results are successfully\nreproduced even for the cases with viscosity variations of several orders of\nmagnitude. The robustness of the numerical method against viscosity variation\ncan be significantly improved by increasing the pre- and post-smoothing\ncalculations during the multigrid operations, and the convergence can be\nachieved for the global viscosity variations up to 1e10.",
"arxiv_id": "physics/0410249",
"authors": [
"Masanori Kameyama",
"Akira Kageyama",
"Tetsuya Sato"
],
"categories": [
"physics.geo-ph"
],
"doi": "10.1016/j.jcp.2004.11.030",
"title": "Multigrid iterative algorithm using pseudo-compressibility for three-dimensional mantle convection with strongly variable viscosity",
"url": "https://arxiv.org/abs/physics/0410249"
},
"schema_id": "dorsal/arxiv",
"source": {
"execution_id": "61081e91-5d95-48f4-a804-e1e36f3e5e0b",
"id": "arXiv Dataset IDs",
"type": "Model",
"variant": "snapshot-2026-03-01",
"version": "0.1.0"
},
"user_id": 1000002
}