CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

discontinuity detection

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 15, 1999, 18:30
Default discontinuity detection
  #1
linhbao tran
Guest
 
Posts: n/a
hi,

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

linhbao
  Reply With Quote

Old   March 16, 1999, 06:15
Default Re: discontinuity detection
  #2
Shigunov
Guest
 
Posts: n/a
Which discontinuities? Shock waves, contact discontinuities, or? And detection is in a numerical simulation or in a experiment? With best wishes

  Reply With Quote

Old   March 16, 1999, 06:23
Default Re: discontinuity detection
  #3
linhbao tran
Guest
 
Posts: n/a
both forms of discontinuities, shock and contact, in a numerical simulation. Regards.
  Reply With Quote

Old   March 17, 1999, 05:42
Default Re: discontinuity detection
  #4
Shigunov
Guest
 
Posts: n/a
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
  Reply With Quote

Old   March 17, 1999, 06:37
Default Re: discontinuity detection
  #5
linhbao tran
Guest
 
Posts: n/a
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.
  Reply With Quote

Old   March 17, 1999, 10:24
Default Re: discontinuity detection
  #6
John C. Chien
Guest
 
Posts: n/a
(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).
  Reply With Quote

Old   March 18, 1999, 05:11
Default Re: discontinuity detection
  #7
linhbao tran
Guest
 
Posts: n/a
not visually but numerically.
  Reply With Quote

Old   March 18, 1999, 07:17
Default Re: discontinuity detection
  #8
Shigunov
Guest
 
Posts: n/a
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
  Reply With Quote

Old   March 20, 1999, 14:41
Default Re: discontinuity detection
  #9
Muhammad Akram
Guest
 
Posts: n/a
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.
  Reply With Quote

Old   March 21, 1999, 15:43
Default Re: discontinuity detection
  #10
linhbao tran
Guest
 
Posts: n/a
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.
  Reply With Quote

Old   March 22, 1999, 13:49
Default Re: discontinuity detection
  #11
Muhammad Akram
Guest
 
Posts: n/a
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.
  Reply With Quote

Old   March 23, 1999, 13:59
Default Re: discontinuity detection
  #12
Muhammad Akram
Guest
 
Posts: n/a
... and about material interface look for maximum density gradient in the domain where pressure gradient is zero (pressure is continous across interfaces).
  Reply With Quote

Old   March 23, 1999, 21:15
Default Re: discontinuity detection
  #13
K. Sundar
Guest
 
Posts: n/a
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
  Reply With Quote

Old   March 24, 1999, 13:06
Default Re: discontinuity detection
  #14
ThoLi
Guest
 
Posts: n/a
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
  Reply With Quote

Old   March 24, 1999, 21:46
Default Re: discontinuity detection
  #15
Dan Williams
Guest
 
Posts: n/a
>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.
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Autogrid: discontinuity of 90 degrees in shroud Eric Joubarne NUMECA 1 July 21, 2009 04:34
Density discontinuity! Ilias FLUENT 0 September 17, 2004 08:04
How to deal with BC with discontinuity? HT Main CFD Forum 0 April 26, 2004 10:37
Stepwise Discontinuity problem Alan Kang Main CFD Forum 1 February 23, 2000 00:46
Solid-Solid Heat Transfer with Contact Discontinuity Jonas Larsson Main CFD Forum 9 September 10, 1999 15:46


All times are GMT -4. The time now is 17:12.