<div dir="ltr"><div><div>Jeevan,<br><br></div>Oscillations in temperature like that around shocks seem to be fairly common in my experience. Trying to understand how to get a better temperature field, especially for real astrophysical equations of state, was one of the motivations for <a href="http://adsabs.harvard.edu/abs/2015ApJS..216...31Z">Zingale and Katz (2015)</a>. We found that the temperature field is sensitive to both the method you use for predicting interface states to the Riemann solver, and the Riemann solver itself. While our technique was mainly implemented in the code CASTRO, we do show comparisons to FLASH as well. One of the things we were able to concretely say is that turning off the characteristic variable limiting in FLASH improved our results on these particular test problems. However it does not make the oscillations near the interface go away, it just makes them smaller.<br><br></div>However, while we were able to make the temperature field look better for these 1D test problems we used, I don't think we came up with a general understanding of how to improve the temperature field in all cases; possibly there is none. Even with these new interface prediction methods and an accurate Riemann solver, I still see temperature fields just like that in multi-D in some cases (though note that CASTRO is using an unsplit solver and maybe you are not). Still, if you are interested, it might be worth looking into trying these other methods we investigated -- maybe they will work for you.<br><div><div><div><div><br></div></div></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Max Katz<br>Ph.D. Candidate, Department of Physics and Astronomy<br>Stony Brook University<br><a href="http://astro.sunysb.edu/mkatz/" target="_blank">http://astro.sunysb.edu/mkatz/</a><br></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On Wed, Mar 9, 2016 at 2:45 PM, jeevan dahal <span dir="ltr"><<a href="mailto:jeevan.dahal1@gmail.com" target="_blank">jeevan.dahal1@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Dear all, <br><br></div>I am trying to run simulation of  a sine wave interface by shock strength of M 1.5. When I run my simulation using the PPM in FLASH code, I see an alternating temperature profile on the interface of sine wave. I have attached a picture of temperature profile to make it clear. In addition, I see a lot of diffusion in my simulation which is around 5 times the width of the zone even though I have turned the diffusion model to be off. <br><br></div>I would appreciated any insights on solving the temperature problem. In addition, how can diffusion be reduced in a simulation?<br><br></div>Thanks<br></div>Regards,<br></div>Jeevan <br></div>
</blockquote></div><br></div>