EGU21-15308
https://doi.org/10.5194/egusphere-egu21-15308
EGU General Assembly 2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.

New continuity-based velocity interpolation scheme for staggered grids

Taras Gerya1, Thibault Duretz2, and Rass Ludovic3
Taras Gerya et al.
  • 1ETH-Zurich, Institute of Geophysics, Department of Earth Sciences, Zurich, Switzerland (taras.gerya@erdw.ethz.ch)
  • 2Géosciences Rennes UMR 6118, Univ Rennes, CNRS, Rennes, France
  • 3ETH-Zurich, Dept. of Civil, Environmental and Geomatic Engineering, Zurich, Switzerland

In the marker-in-cell method combined with staggered finite differences, Lagrangian markers carrying information on material properties are advected with the velocity field interpolated from the staggered Eulerian velocity grid. With existing schemes, velocity interpolation from the grid points to markers violates (to some extent) mass conservation requirement that causes excess convergence/divergence of markers and opening marker gaps after significant amount of advection. This effect is especially well visible in case of diagonal simple shear deformation along planes that are oriented at 45 degrees to the grid and marker circulation through grid corners.

Here, we present a new second order velocity interpolation scheme that guarantees exact interpolation of normal strain rate components from pressure nodes (i.e. from the locations where these components are defined by solving of the mass conservation equation). This new interpolation scheme is thus applicable to both compressible and incompressible flow and is trivially expendable to 3D and to non-regular staggered grids.

Performed tests show that, compared to other velocity interpolation approaches, the new scheme has superior performance in preserving continuity of the marker field during the long-term advection including the diagonal simple shear deformation and marker circulation through grid corners. We showcase a performance-oriented implementation of the new scheme using Julia language's shared memory parallelisation features. The Julia implementation of the new advection schemes further augments the ParallelStencil.jl  related application collection with advection routines.

How to cite: Gerya, T., Duretz, T., and Ludovic, R.: New continuity-based velocity interpolation scheme for staggered grids, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15308, https://doi.org/10.5194/egusphere-egu21-15308, 2021.

Displays

Display file