Journal articles
A free surface capturing discretization for the staggered grid finite difference scheme
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.
Cited literature [29 references]
https://hal-insu.archives-ouvertes.fr/insu-01271289
Contributor : Isabelle Dubigeon <>
Submitted on : Thursday, February 11, 2016 - 8:16:20 AM
Last modification on : Thursday, January 14, 2021 - 11:41:19 AM
Long-term archiving on: : Thursday, May 12, 2016 - 5:29:37 PM
Contributor : Isabelle Dubigeon <>
Submitted on : Thursday, February 11, 2016 - 8:16:20 AM
Last modification on : Thursday, January 14, 2021 - 11:41:19 AM
Long-term archiving on: : Thursday, May 12, 2016 - 5:29:37 PM
File
PL03542.pdf
Publisher files allowed on an open archive