HAL will be down for maintenance from Friday, June 10 at 4pm through Monday, June 13 at 9am. More information
Skip to Main content Skip to Navigation
Journal articles

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

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.
Complete list of metadata

https://hal-insu.archives-ouvertes.fr/insu-03641599
Contributor : Nathalie Pothier Connect in order to contact the contributor
Submitted on : Thursday, April 14, 2022 - 2:04:41 PM
Last modification on : Saturday, April 16, 2022 - 3:39:34 AM

File

Geochem Geophys Geosyst - 2013...
Publisher files allowed on an open archive

Licence


Copyright

Identifiers

Citation

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, AGU and the Geochemical Society, 2013, 14, pp.1200-1225. ⟨10.1002/ggge.20078⟩. ⟨insu-03641599⟩

Share

Metrics

Record views

27

Files downloads

1