<div dir="ltr">On Wed, Oct 30, 2013 at 5:43 PM, Eric Lindsey <span dir="ltr"><<a href="mailto:elindsey@ucsd.edu" target="_blank">elindsey@ucsd.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div><div>I'd like to apply background stresses to a simple elastic 2D plane strain model in pylith, then add fault slip or other quasistatic deformation; later I'll modify the material properties etc. Of course the stresses also cause some deformation in the volume that I'd like to ignore. So I'm imposing the slip at the second time step, with the hopes of subtracting out the deformation from the first time step. (Is this the right way to do this?)</div>

<div><br></div><div>However, I noticed that the effects from the boundary conditions continue to increase, and then oscillate, and don't stabilize until I let it run to the 8th step or so (ideally I'd only use 2 time steps, I think). I think the problem is with the solver's convergence; I see the message </div>

<div><br></div><div>...</div><div>499 KSP Residual norm 9.952574507238e-04</div><div>500 KSP Residual norm 9.952555275971e-04</div><div>Linear solve did not converge due to DIVERGED_ITS iterations 500</div><div><br></div>

<div>I've taken most of the configuration straight from examples/2d/greensfns/strikeslip, and just modified the boundary conditions from Dirichlet to Neumann. Config file is attached; any advice would be great.</div>
</div></blockquote><div><br></div><div>This is a really bad solver. I am sure Brad and Charles can comment on the model setup, but this will not solve</div><div>anything fast. I see that other examples have this basic setup. We definitely did not use it in the online tutorials.</div>
<div>My solver configuration files are in share/solver in the repository. Did you checkout from Git?</div><div><br></div><div>  Thanks,</div><div><br></div><div>      Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>Thanks,</div><div>Eric</div><div><br></div></div>
<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></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
</div></div>