[CIG-SHORT] PyLith: Stresses not transmitted throughout the domain from Drichlet BCs on CUBIT Mesh

[Mamun Miah]

Thanks, Charles. Your comments are very helpful. I have been going through
each of your comments and discovered that both KSPSolve and SNESSolve do
not converge when I uncommented corresponding flag options. So I figured
the problem is with the mesh creation. Also set the the normalizer
according to the smallest mesh dimension. Attached here is the log file in
case you want to see.

Thanks again,
Mamun

On Tue, Oct 13, 2015 at 12:37 PM, Charles Williams <willic3 at gmail.com>
wrote:

> Dear Mamun,
>
> It is a bit difficult to diagnose your problem without more information.
> However, here are several potential problems:
>
>
> 1.  First of all, the mesh you’re using is not of very high quality.  The
> method you’ve used for refinement is pretty crude, and will never result in
> high quality elements. Trelis/Cubit has a number of more powerful methods
> for creating a high quality mesh, and I would take advantage of those.
> 2.  It is unclear why you have set the normalization length scale to 1 m,
> since I believe this is much smaller than your smallest element.
> 3.  There is no indication of whether the linear and nonlinear solutions
> are converging, since you have turned off all the solution monitoring
> options, including the journal options.  You should never do this, as there
> is then no way to tell what has happened with your solution.
>
> You should try running your problem again with all the solution monitoring
> options turned back on, and save the results to a log file.  This will help
> us diagnose your problem further.  I would also try improving your mesh
> quality and selecting a normalization length scale that is appropriate for
> your mesh (~equal to the smallest cell dimension).
>
> Cheers,
> Charles
>
>
> On 14/10/2015, at 6:26 am, Mamun Miah <mmiah at go.olemiss.edu> wrote:
>
> Hello Developers,
>
> I have created the attached geometry (5kmx5kmx4.1km) in CUBIT. It has
> three layer of materials (2000m+100m+200m) and a vertical fault (along Z
> axis) 1000m to the left from the positive X-plane boundary. The exported
> Exodus mesh is then used to run a PyLith simulation with Dirichlet
> displacement and displacement rate BCs (e.g. PyLith Tutorial 'step12.cfg')
> on +x and -x sides. The simulation runs fine but even after running for 100
> years I don't see any stress being transmitted to the interior domain. All
> I can see is the stress changes in the boundary elements where the
> Dirichlet BCs are applied. I am not sure if it is caused by improper
> meshing or boundary conditions.
>
> Attached also are the .cfg files, geometry journal file, and stress output
> file from simulation. I would highly appreciate any help.
>
> --
> Thanks & Best Regards,
>
> Mamun Miah
> Graduate Research Assistant
> Department of Civil Engineering
> The University of Mississippi
> Oxford, MS 38655
> Phone: (662)-801-1904
> Website: http://www.linkedin.com/pub/mamun-miah-eit/21/53/339
> <geometry.jpg><application.cfg><pylithapp.cfg><visit0000.png>
> <application.jou>_______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
>
> Charles A. Williams
> Scientist
> GNS Science
> 1 Fairway Drive, Avalon
> PO Box 30368
> Lower Hutt  5040
> New Zealand
> ph (office): 0064-4570-4566
> fax (office): 0064-4570-4600
> C.Williams at gns.cri.nz
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>



-- 
Thanks & Best Regards,

Mamun Miah
Graduate Research Assistant
Department of Civil Engineering
The University of Mississippi
Oxford, MS 38655
Phone: (662)-801-1904
Website: http://www.linkedin.com/pub/mamun-miah-eit/21/53/339
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20151014/26cb4377/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log
Type: application/octet-stream
Size: 18698 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20151014/26cb4377/attachment-0001.obj>