<FONT face="Default Sans Serif,Verdana,Arial,Helvetica,sans-serif" size=2><div>Scott,</div><div><br></div><div>I think your KSP tolerances need to be adjusted. Please send your PETSc settings used in your .cfg file.</div><div><br></div><div>Brad</div><div><font color="#990099"><br></font></div><div><font color="#990099">-----cig-short-bounces@geodynamics.org wrote: -----</font></div><font color="#990099"><br></font><blockquote style="padding-right:0px;padding-left:5px;margin-left:5px;border-left:#000000 2px solid;margin-right:0px">To: Matthew Knepley &lt;knepley@mcs.anl.gov&gt;<br>From: Scott Henderson &lt;sth54@cornell.edu&gt;<br>Sent by: cig-short-bounces@geodynamics.org<br>Date: 10/30/2012 08:15AM<br>cc: "cig-short@geodynamics.org" &lt;cig-short@geodynamics.org&gt;, Scott Henderson &lt;scottyhq@gmail.com&gt;<br>Subject: Re: [CIG-SHORT] small-scale crack opening convergence problem<br><br><!--Notes ACF
<meta http-equiv="Content-Type" content="text/html charset=utf8">--><div>Here is an output log with all the flags on. I'm thinking maybe the problem is not sufficiently constrained by the boundary conditions I've applied here...</div><div><br></div><div></div><!--Notes ACF
<meta http-equiv="Content-Type" content="text/html charset=utf8">--><div></div><div><br></div><br><div><div>On Oct 30, 2012, at 11:03 AM, Matthew Knepley &lt;<a href="mailto:knepley@mcs.anl.gov">knepley@mcs.anl.gov</a>&gt; wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">On Mon, Oct 29, 2012 at 11:25 AM, Scott Henderson <span dir="ltr">&lt;<a href="mailto:scottyhq@gmail.com" target="_blank">scottyhq@gmail.com</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello CIG-users,<br><br>
I've been trying to use Pylith for a "single edge notch in tension (SENT)" fracture mechanics simulation. In order to do this I've modified step20.cfg to be a 2D plane-stress problem. I've also followed what I could from this mailing-list post (<a href="http://www.geodynamics.org/pipermail/cig-short/2012-September/001120.html" target="_blank">http://www.geodynamics.org/pipermail/cig-short/2012-September/001120.html</a>). I'm still not able to get expected results, which makes me think I' doing something incorrectly with the pylith settings (PETSC in particular). I'm pretty new to the software, so any tips on the cubit mesh or pylith cfg files would be appreciated (zipped folder attached).<br>
<br>
Essentially, I'm not getting the expected u2 displacements from this analysis (but they are not unreasonable), and the following SNES output makes me concerned:<br><br>
500 SNES Function norm 2.483466683115e-06<br>
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 500<br><br>
Increasing snes_max_it doesn't seem to lead to convergence.<br></blockquote><div><br></div><div>Always send all output. The nonlinear solver is not converging, but the first thing to check is</div><div>how the linear solver is doing inside. Also, I want to see exactly what solver setup you have.</div>
<div><br></div><div>&nbsp; &nbsp; Matt</div><div>&nbsp;</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Any suggestions appreciated,<br>
Scott<br><br><br><br>
&nbsp;<br>_______________________________________________<br>
CIG-SHORT mailing list<br>
<a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br><br></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
_______________________________________________<br>CIG-SHORT mailing list<br><a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short<br></blockquote></div><br><font face="Courier New,Courier,monospace" size="3">_______________________________________________<br>CIG-SHORT mailing list<br>CIG-SHORT@geodynamics.org<br><a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br></font></blockquote><br><br>[attachment "log.txt" removed by Brad Aagaard/GD/USGS/DOI]<br></FONT>