discontinuity detection
hi,
i am looking for literatures in this area with little luck. if you know any, please recommend. thank you very much. linhbao 
Re: discontinuity detection
Which discontinuities? Shock waves, contact discontinuities, or? And detection is in a numerical simulation or in a experiment? With best wishes

Re: discontinuity detection
both forms of discontinuities, shock and contact, in a numerical simulation. Regards.

Re: discontinuity detection
There are two approaches, namely interfacetracking and interfacecapturing, and a lot of variants. You can find almost all in J. Comput. Phys. and in Int. J. Numer. Methods in Fluids for the last fourfive years. With best wishes

Re: discontinuity detection
maybe you misunderstood my post. given a flow field with a captured shock or contact surface, i'd like to know how to detect the location of these discontinuities, e.g. by fourier analysis, etc. thanks for responses and best regards.

Re: discontinuity detection
(1).pressure contour plot and velocity contour plot should give you the pressure jump or shear layer locations. By changing the contour levels, you can zoom in a particular location for detail picture. (2). If you use banded contour plot to paint area between two contour lines, you can get a bitmap representation of the flow field ( instead of line contour plot). (3). I think, you can use tools in pattern recognition to study the bitmap picture just like any digitized photo to remove the low gradient region. (4). Regardless of the approach used, the result can only be determined by your eyes and brain. I think the first method of using selective contour level plot of pressure and velocity is fun and easy to do . This is also an effective way to perform adaptive mesh refinement ( either based on the gradient or on the location).

Re: discontinuity detection
not visually but numerically.

Re: discontinuity detection
The method for discontinuity capturing in VolumeofFluid (VOF) method will work in any case, but some expensively, especially for highorder representation. Best regards

Re: discontinuity detection
As I understand the problem is that of detecting shock and interface discontinuties during the computation. If you are solving in Lagragian (reference) description with shock fitting then no problem you always know where these discontinuties are. In Eulerian description they are a bit difficult to track. In this case for shock wave: look for the pressure as well as pressure gradient maxima. The shock is located where both coincide. For material interfaces you can look for maximum density gradients. A small routine can do it. But interfaces are always difficult to track in Eulerian description.

Re: discontinuity detection
at last...someone understood. you mentioned the shock is located where both presure and pressure gradient maxima coincide, yes, but in a computation, there's errors and it's not so easy...i think looking into maximum entropy is probably on way. i've read somewhere..have to go back and check.

Re: discontinuity detection
May be you are right. I told you my experience. Now I recall something more straight forward and easy. The shock front is located in your computational domain where the local acceleration a_i+1/2 = (p_i+1  p_i)/((rho_i+1/2)(x_i+1 x_i)) is maximum. Please look into the book "Numerical Simulation of reactive flow" by Oran and Boris, Elsevier, 1987.

Re: discontinuity detection
... and about material interface look for maximum density gradient in the domain where pressure gradient is zero (pressure is continous across interfaces).

Re: discontinuity detection
Hi
If you are looking at "entropy criterion" for designing schemes which are used for capturing shocks and contact discontinuities, you can refer to papers by Harten in Journal of Computational Physics on Total Variation Diminishing schemes. The essential idea is at the discontinuity (shock), the entropy maximises. Some numerical schemes exploit this concept to design the numerical algorithm (TVD schemes). The contact discontinuity, for such schemes would appear as an expansion shock, (which is obviously unphysical) and so the thermodynamic entorpy criterion is enforced to obtain the physical solution. Several other papers notably by H. C. Yee have been published in J. Comp. Phy. Vols (58 80 range) You can also refer to Ch Hirsch (Text book) for a discussion on the entropy criterion. Bye Sundar 
Re: discontinuity detection
Hi,
the following paper describes a shock indicator: U. Göhner, G. Warnecke: A Shock Indicator for Adaptive Transonic Flow Computations. Numer. Math. 66 (1994) 423448. the indicator is used to move nodes exactly on the shockline to reduce blurring. This article uses a potential flow (instead of velocity and pressure). You can find more information about the second author at G. Warnecke's Curriculum vitae. greetings ThoLi 
Re: discontinuity detection
>at last...someone understood. you mentioned the shock is >located where both presure and pressure gradient maxima >coincide, yes, but in
>a computation, there's errors and it's not so easy...i >think looking into maximum entropy is probably on way. i've >read somewhere..have >to go back and check. You can get very sophisticated and use things like VOF to get contact discontinuities, and shock tracking for shocks. However, if you just have a regular old Euler code, gradient detection algorithms work fine. I used this approach in my own work, for tracking detonations, and had no real problems. The main drawback is that checking for large gradients is a little crude, only works well for simple shock structures. More complex flow structures require a lot of tuning and fiddling around before the tracking works well. Deep down, I never really liked it all that much, but still, it worked fine for me at the time. Dan W. 
All times are GMT 4. The time now is 03:15. 