An adaptive staggered grid finite difference method for modeling geodynamic Stokes flows with strongly variable viscosity - Archive ouverte HAL Access content directly
Journal Articles Geochemistry, Geophysics, Geosystems Year : 2013

An adaptive staggered grid finite difference method for modeling geodynamic Stokes flows with strongly variable viscosity

, , (1)
1
T. V. Gerya
  • Function : Author
D. A. May
  • Function : Author

Abstract

Here we describe a new staggered grid formulation for discretizing incompressible Stokes flow which has been specifically designed for use on adaptive quadtree-type meshes. The key to our new adaptive staggered grid (ASG) stencil is in the form of the stress-conservative finite difference constraints which are enforced at the "hanging" velocity nodes between resolution transitions within the mesh. The new ASG discretization maintains a compact stencil, thus preserving the sparsity within the matrix which both minimizes the computational cost and enables the discrete system to be efficiently solved via sparse direct factorizations or iterative methods. We demonstrate numerically that the ASG stencil (1) is stable and does not produce spurious pressure oscillations across regions of grid refinement, which intersect discontinuous viscosity structures, and (2) possesses the same order of accuracy as the classical nonadaptive staggered grid discretization. Several pragmatic error indicators that are used to drive adaptivity are introduced in order to demonstrate the superior performance of the ASG stencil over traditional nonadaptive grid approaches. Furthermore, to demonstrate the potential of this new methodology, we present geodynamic examples of both lithospheric and planetary scales models.
Fichier principal
Vignette du fichier
Geochem Geophys Geosyst - 2013 - Gerya - An adaptive staggered grid finite difference method for modeling geodynamic Stokes.pdf (3.52 Mo) Télécharger le fichier
Origin : Publisher files allowed on an open archive

Dates and versions

insu-03641599 , version 1 (14-04-2022)

Licence

Copyright

Identifiers

Cite

T. V. Gerya, D. A. May, T. Duretz. An adaptive staggered grid finite difference method for modeling geodynamic Stokes flows with strongly variable viscosity. Geochemistry, Geophysics, Geosystems, 2013, 14, pp.1200-1225. ⟨10.1002/ggge.20078⟩. ⟨insu-03641599⟩
41 View
13 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More