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

ENO scheme: (Shu and Osher)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 16, 2002, 10:52
Default ENO scheme: (Shu and Osher)
  #1
J. K
Guest
 
Posts: n/a
Hi all,

I am implementing an ENO's scheme in a F.V. method. I am trying to write a simple 1-D code at present. However, I know that I need to use the divided differences (I want to use a 3rd order accurate polynomial) to calculate the "smoothness" and hence the relevant node to add to my stencil. However, lets say my fluid is going from left to right, and I have 100 cells. Now, if I want to implement a scheme which uses the 3rd divided differences, then, say, if my current cell is cell 98 (or 99), I don't have enough cells to the right of cell 98 (or 99) to complete the calculation of a 3rd divided difference. My question is, if this is the case, then do I have to use an ENO's scheme of lower order at the right boundary, and also, how do I apply the ooutflow B.C. The former question is the most pressing. Thank you all. John.
  Reply With Quote

Old   July 16, 2002, 12:38
Default Re: ENO scheme: (Shu and Osher)
  #2
Axel Rohde
Guest
 
Posts: n/a
Use more than one ghost cell outside the domain (say, 101, 102, 103), and fill them according to your inflow / outflow B.C. You may have to be creative there, i.e extrapolate gradients, use Riemann invariants, etc.
  Reply With Quote

Old   July 17, 2002, 06:20
Default Re: ENO scheme: (Shu and Osher)
  #3
J. K
Guest
 
Posts: n/a
Hi Alex,

Thanks for you help. I was afraid I might have to use ghost cells. I am afraid I don't know alot about riemann invariants apart than that they are the characteristic variables for the euler equations. I am using a simple linear consevation law.

Can you give me more information regarding you comments? If I create ghost cells, I will have to give them values, and this will influence what goes on within the domain. Have you any references which may be of use to me? Also, Would it be possible to use polynomials of lower order as I approach the wall, such that 2 cells from the wall, I have a polynomial which is linear etc...? I can use a polynomial which is entirely upwind at the last cellsnear the outflow....but when it comes to calculating the divided differences up to say the 3rd divided difference, I need cells outside the domain for the calculating the divided differences starting at cells n, n-1 and n-2, where n is the total number of cells. Could I simply use lower divided differences close to the exit to choose which cell to add to my stencil, or does this dramactically affect the accuracy? Any references or ideas?
  Reply With Quote

Old   July 18, 2002, 12:51
Default Re: ENO scheme: (Shu and Osher)
  #4
Axel Rohde
Guest
 
Posts: n/a
J.K.

Although I have read a few papers in the past on ENO, I am not that familiar any more with its polynomial fitting. I usually just stick to 2nd order TVD schemes, which become first order near extrema. Personally I would not bother with Riemann invariants. Other, simpler methods can be just as effective.

For a 1-D problem (pipe or shock tube) you only have two types of B.C.'s, wall or open end. For the open end, I would fill the ghost cells with ambient conditions. For the wall I would use a reflection B.C. To fill more than one wall ghost cell you can extrapolate gradients from the interior. An example of gradient extrapolation can be found on page 62 of my dissertation, which can be downloaded from my website,

www.microcfd.com/papers.htm

My dissertation also includes a big chapter on the shock tube problem, if that is what you are working on. With TVD and ENO schemes, there are a number of intermediate quantities which also have to be extrapolated somehow. These quantities are denoted in 'Fraktur' type set in my dissertation. With those intermediate quantities, one has to be 'creative' at times near B.C.'s.

I hope this helps!

Axel
  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
ENO scheme (ENO-Roe) (Shu and Osher). J. K. Main CFD Forum 8 March 28, 2013 09:09


All times are GMT -4. The time now is 00:57.