<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hi all,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">I've been doing some runs of supernova remnant evolution with thermal conduction. I was surprised to find that the shock propagates faster when I include thermal conduction, even with flux limitation. I attribute this to conduction in the shock front, which is necessarily spread over a few grid points by numerical effects. </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">This is mostly, I believe, not a physically realistic result, so I'd like to turn off conduction in the shock fronts. The approach that I'm thinking of using is related to the one used in the Burn unit wherein Hydro_detectShock is called which fills an array shock with 0 (no shock) or 1 (shock). Then in the Diffuse unit, diff_advanceTherm to be precise, I would call Hydro_detectShock and set cond_zone to 0 if shock is 1. Does this sound like a reasonable approach? Has anyone else done something similar? Are there pitfalls I need to watch out for?</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Thanks in advance for any advice.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Regards,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Jon</div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<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: (617) 496-7981 Cambridge, MA 02138-1516<br>cell: (781) 363-0035 USA<br>________________________________________________________<br><br></div></div></div></div>
</div>