A free surface capturing discretization for the staggered grid finite difference scheme - Archive ouverte HAL Access content directly
Journal Articles Geophysical Journal International Year : 2016

A free surface capturing discretization for the staggered grid finite difference scheme

(1) , (2) , (3)
1
2
3

Abstract

The coupling that exists between surface processes and deformation within both the shallow crust and the deeper mantle-lithosphere has stimulated the development of computational geodynamic models that incorporate a free surface boundary condition. We introduce a treatment of this boundary condition that is suitable for staggered grid, finite difference schemes employing a structured Eulerian mesh. Our interface capturing treatment discretizes the free surface boundary condition via an interface that conforms with the edges of control volumes (e.g. a ‘staircase’ representation) and requires only local stencil modifications to be performed. Comparisons with analytic solutions verify that the method is first-order accurate. Additional intermodel comparisons are performed between known reference models to further validate our free surface approximation. Lastly, we demonstrate the applicability of a multigrid solver to our free surface methodology and demonstrate that the local stencil modifications do not strongly influence the convergence of the iterative solver.
Fichier principal
Vignette du fichier
PL03542.pdf (1.92 Mo) Télécharger le fichier
Origin : Publisher files allowed on an open archive
Loading...

Dates and versions

insu-01271289 , version 1 (11-02-2016)

Identifiers

Cite

Thibault Duretz, D.A. May, Philippe Yamato. A free surface capturing discretization for the staggered grid finite difference scheme. Geophysical Journal International, 2016, 204 (3), pp.1518-1530. ⟨10.1093/gji/ggv526⟩. ⟨insu-01271289⟩
173 View
429 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More