Tecplot shock detection
 

Hi everyone,
I'm looking for an algorithm to detect a shock during post processing. I found a tool in Tecplot 360 which does this work quite well: does anyone know which algorithm is employed in Tecplot?
Thank you very much

Andrea Ferrero

Andrea,

On page 384 the Tecplot 360 user manual states ( http://download.tecplot.com/360/current/tpum.pdf ) :


http://raphael.mit.edu/fx/. These features are accessed through the Extract Flow Features dialog.

This paragraph and reference may be an answer to your question.

Scott

Thank you very much, Scott.
I find the article "Shock Detection From Computational Fluid Dynamics Results" cited in your link very useful, in particular the algorithm for moving shocks.

Andrea

Has anyone any idea that why the value of "shock" variable in Tecplot can be greater than 1 or less than -1? I tried to capture a shock interface in a 2D convergent-divergent nozzle. Although the position is deriable, but the lengend given in Tecplot has a value which can be as high as 30 and as low as -30. According to the definition,
it should be (U/a) dot (grad(P))/|grad(P)|, how could it be larger than 1?
I have read the reference listed above, it says for stationary shock the the location is the place where the value of Ma vector dot pressure gradient should be one. Can anybody tell me what's going on there? Thank you in advance.

Due to numerical error of the interpolation scheme, small and undirected pressure gradient may be calculated in uniform flow region. Then (U/a) dot (grad(P))/|grad(P)| will vary from -Ma to Ma, where Ma is the local Mach number.

The following link http://authors.elsevier.com/sd/artic...00936113001143 is our paper related to shock detection and post-processing, which would be in your interest.