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

ENO scheme asymmetry

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 27, 2004, 18:00
Default ENO scheme asymmetry
  #1
vinod venugopalan
Guest
 
Posts: n/a
hi,

I am trying to code an ENO 3rd order scheme for Euler equations (conservation form). I have managed to get good results for the shock tube problem. I extended the code to 2D and then found assymetry in the results.

So I went back and tried running a 1D problem, trying to solve a problem with energy input at the center of the domain. The results ought to be a set of symmetrical left and right running pressure waves, however I get asymmetrical results. The right running wave is stronger (about 5-7%) and steeper, i.e with more gradient than the left one.

I am not sure why this is occuring. I am using ENO 3 FD component-wise Lax Friedrich flux splitting,

fluxp[i] = 0.5*(f(q1[i],q2[i],q3[i]) + alpha*us[i]); fluxn[i] = 0.5*(f(q1[i],q2[i],q3[i]) - alpha*us[i]);

After finding the hp[i] & hn[i] value from ENO reconstruction starting my stencil from i for both the above split fluxes (is this correct?), I add them together

for(i=0;i<=n;i++) h[i] = hp[i] + hn[i+1];

Are the above steps correct? They work well for the shock tube problem. No probs there. The asymmetry is really bugging me and I can't find a reason why.. is it the scheme? or something else?

I would be obliged if someone could take the time to reply.
  Reply With Quote

Old   September 27, 2004, 18:12
Default Re: ENO scheme asymmetry
  #2
vinod venugopalan
Guest
 
Posts: n/a
Follow up Q..

If this is indeed a problem with the scheme.. then is there any way I can reverse this.. as in make the left running wave stronger and right one weaker? That way I can may be alternate between the both and get a mean solution..

Thanks

Vinod
  Reply With Quote

Old   September 28, 2004, 20:11
Default Re: ENO scheme asymmetry
  #3
Pete
Guest
 
Posts: n/a
Hi, Vinod,

ENO scheme can't cause this much assymmetry. There may be something wrong with your code. Or you have some misunderstanding on the scheme itself.

The mistake is probably in the evaluation of the variable values at half-grid points. Hope you can find out the bug soon.

Pete
  Reply With Quote

Old   September 30, 2004, 14:14
Default Re: ENO scheme asymmetry
  #4
Vinod Venugopalan
Guest
 
Posts: n/a
thanks,

Is there any way you could mail me a copy or a snippet of your ENO code on how to achieve flux splitting? I am not so sure if the evaluation of variable values at half-grid points are causing the error... the code works just fine for the shock tube problem. I think the problem is with flux splitting and I just can't put my finger on it. You may look at the results @

http://f1vinod.tripod.com/results/results.emf http://f1vinod.tripod.com/results/shocktube.emf http://f1vinod.tripod.com/results/shocktube2.emf

I would appreciate it if you could mail me the code or the code snippet.

thanks Vinod (f1vinod@yahoo.com)
  Reply With Quote

Old   September 30, 2004, 23:56
Default Re: ENO scheme asymmetry
  #5
Pete
Guest
 
Posts: n/a
Funny, your email doesn't work. My email was bounced back.

And I can't view your links either.

Pete
  Reply With Quote

Old   October 7, 2004, 01:01
Default Re: ENO scheme asymmetry
  #6
Matthew
Guest
 
Posts: n/a
h[i] = hp[i] + hn[i+1]; Maybe the fault is in here. Try to change the above equation into h[i] = hp[i] + hn[i];
  Reply With Quote

Old   October 7, 2004, 11:05
Default Re: ENO scheme asymmetry
  #7
Vinod Venugopalan
Guest
 
Posts: n/a
I am sorry pete.. the email Id had a typo

f1vinod@yahoo.com or try venugov@email.uah.edu

The links do work when I click them in windows.. I checked the links on my solaris machine and it doesn't work for some reason.

http://f1vinod.tripod.com/results/results.emf

http://f1vinod.tripod.com/results/shocktube.emf

http://f1vinod.tripod.com/results/shocktube2.emf

I tried something else, I split my flux such that negative flux interpolation stencil search starts with (i,i+1) instead of (i,i-1) as in the positive flux. This gives me symmetric results but it also ends up giving me small oscillations in parts of the domain where there's supposed to be a smooth solution (none at the shock or contact discontinuity)..

Please do email me.Appreciate your help.

Thanks

Vinod
  Reply With Quote

Old   October 7, 2004, 11:21
Default Re: ENO scheme asymmetry
  #8
Vinod Venugopalan
Guest
 
Posts: n/a
hi,

I did try that, it didn't work... the code gives me NaNs, it blows up.

fluxpositive[i] --->||node i||<--- fluxnegative[i+1]

Anyway the reason I add h[i] = hp[i] + hn[i+1]; is that hn[i] would be the negative flux for node i-1, since I am performing ENO interpolation starting from (i,i-1) for both negative and postive interpolations. Therefore, hn[i+1] would be the correct negative flux for node i.

I tried something else,I split my flux such that negative flux interpolation stencil search starts with (i,i+1) instead of (i,i-1) earlier. This gives me symmetric results but it also ends up giving me small oscillations in parts of the domain where there's supposed to be a smooth solution (none at the shock or contact discontinuity)..

I am not sure what's going on. If you can send me a code or code snippet that shows how this works, please do.

Thanks Vinod
  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
ENO Scheme in a Shock tube problem Shiv Main CFD Forum 4 January 29, 2007 16:31
ENO scheme asymmetry Vinod Venugopalan Main CFD Forum 5 July 22, 2004 13:06
Entropy condition and ENO scheme Question Man. Main CFD Forum 5 January 17, 2003 08:10
ENO scheme: (Shu and Osher) J. K Main CFD Forum 3 July 18, 2002 12:51


All times are GMT -4. The time now is 03:36.