An adaptive staggered grid finite difference method for modeling geodynamic Stokes flows with strongly variable viscosity - INSU - Institut national des sciences de l'Univers Accéder directement au contenu
Article Dans Une Revue Geochemistry, Geophysics, Geosystems Année : 2013

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

T. V. Gerya
  • Fonction : Auteur
D. A. May
  • Fonction : Auteur

Résumé

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
Origine : Fichiers éditeurs autorisés sur une archive ouverte

Dates et versions

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

Licence

Copyright (Tous droits réservés)

Identifiants

Citer

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⟩
49 Consultations
91 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More