[CIG-SHORT] Benchmark Pylith with viscoelastic analytic solution for a pressurized sphere

[Francisco Delgado]

Dear Pylith developers,

I'm trying to simulate the ground displacement produced by a small
pressurised cavity surrounded by a viscoelastic shell (equation 7.105 from
Segall's Earthquake and Volcano Deformation book) and Pylith does a good
job within a few mm of uncertainty. The difference between the analytic and
the numerical solution is at most 0.6 mm for a signal with a maximum
amplitude of 2.49 cm, however I've noticed that this last figure increases
as time goes forward, very small but in the end negligible variations
(evident in the attached plot, the profiles are calculated for the point of
maximum uplift).

[image: Inline image 1]. Is this a typical issue of numerical viscoelastic
models or is it related to the relatively low resolution of the finite
element mesh??

Thanks

-- 
--------------------------------------------------------
Francisco Delgado
PhD student in Geological Sciences
Cornell University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20160926/9a6e2799/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sphere_visco.png
Type: image/png
Size: 107940 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20160926/9a6e2799/attachment-0001.png>