[CIG-SHORT] Pylith question

[Charles Williams]

Dear Xiaoxi,

Looking at your mesh, it is very coarse for this type of problem, particularly since you are using linear tets.  I suspect that your results will improve with a finer mesh.  As Matt mentions, getting stresses from an initial elastic problem may also help.  Keep in mind that the FE solution in general is not exact, particularly for coarse meshes.  It might be helpful to look at the difference in stresses between the elastic solution and the analytical solution.  You can either do this in Python or using the calculator in Paraview.

Cheers,
Charles


> On 29/03/2017, at 5:17 PM, Xiaoxi Zhao <xiaoxiz at usc.edu> wrote:
> 
> To whom it will concern,
> 
> My name is Xiaoxi. I'm learning Pylith these days and trying to build a model with multiple layers and faults. The model is assigned gravity and initial stress equilibrium to gravity. Five surfaces are fixed with free ground surface. I expect the result to be zero displacement in the domain, however, there is some subsidence above the layers. I'm not sure where the problem is. Can you please help me check the files and let me know if anything need to be changed to solve the problem
> 
> Thanks for your time and help,
> Best regards,
> Xiaoxi
> 
> -- 
> Sincerely,  
> 
> Xiaoxi Zhao
> 
> Mork Family Dept, Petroleum Engineering
> USC Viterbi School of Engineering
> University of Southern California
> 
> Los Angeles, CA 90089
> Tel: 213-447-6366 <>
> 
> xiaoxiz at usc.edu <mailto:xiaoxiz at usc.edu><test.zip>_______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short

Charles Williams I Geodynamic Modeler
GNS Science I Te Pῡ Ao
1 Fairway Drive, Avalon 5010, PO Box 30368, Lower Hutt 5040, New Zealand
Ph 0064-4-570-4566 I Mob 0064-22-350-7326 I Fax 0064-4-570-4600
http://www.gns.cri.nz/ <http://www.gns.cri.nz/> I Email: C.Williams at gns.cri.nz <mailto:your.email at gns.cri.nz>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170330/4e889c1e/attachment.html>