<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Dec 4, 2013 at 11:23 PM, Eric Lindsey <span dir="ltr"><<a href="mailto:elindsey@ucsd.edu" target="_blank">elindsey@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Hi guys,<div><br></div><div>I'm having some trouble getting the model with a fault in an elastoplastic rheology (plane strain) to converge. I'm pretty sure I've fixed all the boundary conditions to behave correctly now. The LU method fails entirely once yielding occurs, and the approximate method (multiplicative fieldsplit using ml and jacobi for the two domains) does not converge. More details:</div>

<div><br></div><div>The first thing I tried was using (solver04.cfg) from the tutorial, which is a schur fieldsplit with LU and the custom Pylith preconditioner. In the case of a fault in a purely elastic material, this procedure is working fine. The nonlinear solve does take over 100 iterations, I'm not sure if this is just a function of my mesh; anyway I'm not too worried about optimization yet.</div>
</div></blockquote><div><br></div><div>Is this a through going fault? The linear solve is exact, so we would not expect this.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div dir="ltr"><div>In the plastic case, whenever stresses near the bends in the fault begin to exceed the yield criterion (I increase them slowly to this value over several time steps), I get a Zero Pivot error. But maybe this is not unexpected for a plastic rheology?</div>
</div></blockquote><div><br></div><div>Using this option</div><div><br></div><div>fs_fieldsplit_0_pc_factor_shift_type = nonzero<br></div><div><br></div><div>will prevent a 0-pviot, but the preconditioner becomes weaker.</div>
<div><br></div><div>   Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">
<div>I also tried the solver options for an approximate solution with multiplicative fieldsplit suggested by Brad (solver08.cfg): The elastic case still works fine; it uses way more KSP iterates but still finishes in about the same total time. But the plastic case still throws a Zero Pivot. When I ran it earlier without increasing snes_max_it, I got a string of these messages instead:</div>

<div>
<br></div><div>"Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 0"</div><div><br></div><div>In either case, I think I must need a better solver. I've gone back through the tutorial videos for the solvers, but I'm not very clear on how to apply the nonlinear options from the driven cavity problem to my situation, or really where to start. My input files are attached, but I'm happy to send any additional files or output as needed. Any insight, suggestions, or wild guesses?</div>

<div><br></div><div>Thanks</div><span class=""><font color="#888888"><div>Eric</div>
<div><br></div><div><br></div></font></span></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>