# ENO scheme asymmetry

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 September 27, 2004, 18:00 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.

 September 27, 2004, 18:12 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

 September 28, 2004, 20:11 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

 September 30, 2004, 14:14 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)

 September 30, 2004, 23:56 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

 October 7, 2004, 01:01 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];

 October 7, 2004, 11:05 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

 October 7, 2004, 11:21 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post J. K. Main CFD Forum 8 March 28, 2013 09:09 Shiv Main CFD Forum 4 January 29, 2007 16:31 Vinod Venugopalan Main CFD Forum 5 July 22, 2004 13:06 Question Man. Main CFD Forum 5 January 17, 2003 08:10 J. K Main CFD Forum 3 July 18, 2002 12:51

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

 Contact Us - CFD Online - Top