[aspect-devel] Van Keken et al. 2008 subduction benchmark

[Max Rudolph]

I was wondering if anyone has been successful in setting up the Van Keken
et al. 2008 subduction benchmark using ASPECT? I have been using this as a
starting point to set up a more complicated corner flow model and have not
yet been successful.

This benchmark solution requires a rigid overriding plate and kinematic
subducting plate. There is a discontinuity along the subduction interface.
I tried imposing kinematic boundary conditions on the left and bottom of
the box and using Jonathan Perry-Houts' prescribed velocity plugin
(cookbooks/prescribed_velocity), but the Stokes solver fails to converge.
Looking at the pressure field, it's not surprising to see very large
oscillations in this region where the velocity field is discontinuous. It
occurred to me that using the locally conservative discretization might
help. However, when the locally conservative discretization is enabled, the
internal prescribed velocities are no longer enforced anywhere. Is it
obvious why this might be the case?

Regards,
Max
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161027/b4c8b750/attachment.html>