<div dir="ltr">that should have read: "to separate shocks from contacts"</div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 23, 2015 at 12:04 PM, Michael Zingale <span dir="ltr"><<a href="mailto:michael.zingale@stonybrook.edu" target="_blank">michael.zingale@stonybrook.edu</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">I don't have yt code, but I agree that this would be a useful standard derived field to add.  As for the methodology, the shock detection algorithm that is used in a lot of codes looks for regions where there is compression and a pressure jump (to separate shocks from rarefactions).  Most PPM or similar hydro codes (even those that are not directionally split) will do this check in 1-d when determining whether to apply flattening to shocks (to prevent their self-steepening from making them too steep), and then look in all directions to figure out if a zone should be flattened.  You want to do a single multi-d shock detection, which will mean computing the velocity divergence and the pressure jump in the normal direction.  Here's the code we use in Castro which inherits some ideas from other codes (Flash, and the ppm papers).  Look for the subroutine shock().<div><br></div><div><a href="https://github.com/BoxLib-Codes/Castro/blob/master/Source/Src_3d/riemann_3d.f90" target="_blank">https://github.com/BoxLib-Codes/Castro/blob/master/Source/Src_3d/riemann_3d.f90</a><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Tue, Jun 23, 2015 at 11:51 AM, Yuan Li <span dir="ltr"><<a href="mailto:yuan@astro.columbia.edu" target="_blank">yuan@astro.columbia.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Dear all,<div><br></div><div>I am trying to figure out a way to mask out cells that contain shocks. One way I can think of is to create a derived field by identifying shocks. It seems that there are different ways of finding shocks: Enzo has a refinement criterion (refine by shocks) and there is also shock finding written by Sam. </div><div><br></div><div>Has anybody tried to find shocks using yt? I would appreciate it if you would like to share your code/experience/advice. Thank you!</div><span><font color="#888888"><div><br></div><div>Yuan</div></font></span></div>
<br></div></div>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div>Michael Zingale</div><div>Associate Professor</div><div><br></div><div>Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY 11794-3800</div><div><i>phone</i>:  <a href="tel:631-632-8225" value="+16316328225" target="_blank">631-632-8225</a></div><div><i>e-mail</i>: <a href="mailto:Michael.Zingale@stonybrook.edu" target="_blank">Michael.Zingale@stonybrook.edu</a></div><div><i>web</i>: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div>Michael Zingale</div><div>Associate Professor</div><div><br></div><div>Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY 11794-3800</div><div><i>phone</i>:  631-632-8225</div><div><i>e-mail</i>: <a href="mailto:Michael.Zingale@stonybrook.edu" target="_blank">Michael.Zingale@stonybrook.edu</a></div><div><i>web</i>: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a></div></div>
</div>