hi,

i am looking for literatures in this area with little luck. if you know any, please recommend. thank you very much.

linhbao

Which discontinuities? Shock waves, contact discontinuities, or? And detection is in a numerical simulation or in a experiment? With best wishes

both forms of discontinuities, shock and contact, in a numerical simulation. Regards.

There are two approaches, namely interface-tracking and interface-capturing, 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 four-five years. With best wishes

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.

(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).

not visually but numerically.

The method for discontinuity capturing in Volume-of-Fluid (VOF) method will work in any case, but some expensively, especially for high-order representation. Best regards

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.

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.

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.

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

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

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) 423-448.

the indicator is used to move nodes exactly on the shock-line 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

>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.