<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi all,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I've been trying to do supernova remnant evolution runs that include electron thermal conduction and MHD using the unsplit staggered mesh unit.  I can do runs that include either electron thermal conduction or MHD but when I try to combine them I get a rapid (after 43 steps) failure with:</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"> [gr_hypreSolve]: Ok ierr=0, component=0, converged=T, |initial RHS|^2=3.6549+240, num_iterations=1, final_res_norm=NaN</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">so a nonconvergence related to thermal conduction, I believe. </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Of course ideally one would use anisotropic thermal conduction (which has been discussed previously on this list), but as an approximation I was going to use a reduced conductivity to see what sort of results I got. I was not expecting this sort of failure, especially so quickly.  </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I tried reducing the conductivity further and found that only by reducing it by many orders of magnitude (~8) was I able to get runs that didn't crash.  In those cases then there is essentially no conduction.  I also tried a slow turn on of thermal conduction and that didn't work either.  I'm baffled about why I should get such failures.  It made me wonder if there is some numerical way that the implementations of thermal conduction and MHD interfere with each other.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Any help on this would be appreciated.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Regards,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Jon</div><div><br></div>-- <br><div class="gmail-m_6096473260365936312gmail_signature"><div dir="ltr"><div><div dir="ltr">______________________________<wbr>__________________________<br>Jonathan D. Slavin                 Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>phone: <a href="tel:(617)%20496-7981" value="+16174967981" target="_blank">(617) 496-7981</a>       Cambridge, MA 02138-1516<br>cell: <a href="tel:(781)%20363-0035" value="+17813630035" target="_blank">(781) 363-0035</a>             USA<br>______________________________<wbr>__________________________<br><br></div></div></div></div>
</div>