<div dir="ltr">Hi Jon! <div><br></div><div><span style="font-family:arial,helvetica,sans-serif">>It made me wonder if there is some numerical way that the implementations of thermal conduction and MHD interfere with each other.</span></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">I have used MHD+(isotropic) Conduction without any problem of that sort. One confusing thing is that you can select Diffusion/split and Diffuse/unsplit. Be careful with that. </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Something curious is that if the simulation include unsplit/MHD + (regardless of </font><span style="font-family:arial,helvetica,sans-serif">split or </span><span style="font-family:arial,helvetica,sans-serif">unsplit) </span><span style="font-family:arial,helvetica,sans-serif">Diffusion</span><span style="font-family:arial,helvetica,sans-serif"> you have to also set the runtime parameter addThermalFlux=.false. (.true. by default)</span></div><div><font face="arial, helvetica, sans-serif">Otherwise, you are including the thermal conduction effect twice, i.e., </font><span style="font-family:arial,helvetica,sans-serif">heat </span><span style="font-family:arial,helvetica,sans-serif">diffusion via the HYPRE libraries decoupling the heat equation, and adding the thermal flux to the total energy flux in the Hydro unit.  </span></div><div><span style="font-family:arial,helvetica,sans-serif">I think the most consistent configuration would be unsplit/MHD + Diffusion/</span><span style="font-family:arial,helvetica,sans-serif">unsplit</span><span style="font-family:arial,helvetica,sans-serif"> + addThermalFlux=.false.</span></div><div><span style="font-family:arial,helvetica,sans-serif"><br></span></div><div><span style="font-family:arial,helvetica,sans-serif"><br></span></div><div><span style="font-family:arial,helvetica,sans-serif">>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 </span><span style="font-family:arial,helvetica,sans-serif"><br></span></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Maybe you can try reducing the parameter dt_diff_factor. </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Are you using super-time-stepping?  </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Best! </font></div><div><div>Ernesto. </div></div><div class="gmail_extra"><br><div class="gmail_quote"><br></div>
</div></div>