CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Too much turbulence

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

Like Tree1Likes
  • 1 Post By andimb

Reply
 
LinkBack Thread Tools Display Modes
Old   March 20, 2006, 09:39
Default Hi! I'm still calculating a
  #1
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

I'm still calculating a flow around a NACA airfoil and compare it to a CFX solution.
I use the LaunderSharmaKE turbulence model with an yplus near 1.
If I start a calculation with the same boundary condiations as in CFX, where epsilon is calculated by epsilon = density*C_mu*k^2/eta_T and eta_T=10*eta where eta is the dynamic viscoity, the solver reaches no steady solution, seen in the first image.

If I use the following formulationfor epsilon: epsilon = k^(3/2)/L with L as an charactersitic length as the chordlength, I get a much lower epsilon and the solution is steady but with a too big boundary layer ass seen in the picture.


I also try other turbulence models and guessed epsilons and finer or coarser meshes, with no dramatic change.

I don't see the mistake I make.

Thanks for help!
Andreas
andimb is offline   Reply With Quote

Old   March 20, 2006, 10:31
Default Are the numerical set-up ident
  #2
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 8
pierre is on a distinguished road
Are the numerical set-up identical in both CFX and Foam?

Pierre
pierre is offline   Reply With Quote

Old   March 20, 2006, 10:45
Default Hi! I think completly ident
  #3
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

I think completly identical is not possible. I used the Hiresolution scheme in CFX which, as I understood, should be equal to the gamma scheme in OpenFOAM.
So it should be equal.
You just have not as much possibilities to choose in CFX as in OpenFOAM. And you need something better than upwind to get right results.

Andreas
andimb is offline   Reply With Quote

Old   March 21, 2006, 12:02
Default Hi! I tried the easiest cas
  #4
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

I tried the easiest case I can think of with simpleFoam. Coarse mesh with no turbulence model (exactly laminar dummy model) and upwind scheme where possible.
There are still seperation bubbles. Perhaps the boundary consditions are wrong but I don't know which?

Perhaps someone can help me.




Thanks
Andreas
andimb is offline   Reply With Quote

Old   March 21, 2006, 12:07
Default Hi, again! I use an own wri
  #5
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi, again!

I use an own written mesh generator. I just give the profil coordinates and the mesh is autogenerated. Perhaps there is something wrong. But checkMesh doesn't throw any error.

Andreas
andimb is offline   Reply With Quote

Old   March 21, 2006, 14:33
Default I recall that I had something
  #6
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 8
pierre is on a distinguished road
I recall that I had something similar early on in my PhD, i.e. my calculation would always show behaviour that looked physical to a point but not what was wanted. eventually it was due to the fact that my generated mesh had inconsistent coordinates (I can't recall if i had used checkMesh to verify ny mesh). In the blockMesh I had messed up some of the points ordering. I don't if that helps. I have you tried to look at the mass continuity, not just locally but what come in and out of the domain.
I am not sure whether that helps.

Pierre
pierre is offline   Reply With Quote

Old   March 22, 2006, 09:14
Default Hi! I checked the mass cont
  #7
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

I checked the mass continuity, but everthing is like it should be. And I checked the mesh by hand and found no mistake.

Andreas
andimb is offline   Reply With Quote

Old   March 22, 2006, 09:42
Default What is your boundary conditio
  #8
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
What is your boundary condition for k and epsilon?
If it is zero gradient, that might be your problem. For low Re models, the BC should be "fixedValue" with the value set to some small non-zero number e.g. 1e-10;
eugene is offline   Reply With Quote

Old   March 22, 2006, 09:53
Default Hi! First I want to thank y
  #9
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

First I want to thank you all for your help!!

Yes, the boundary conditions were zero gradient, so I'll try fixed value. But why could this matter with turbulence model laminar? k and epsilon won't be calculated, right?
This was the reason I did the run laminar, to check if the turbulence boundary conditions are a problem.

Andreas
andimb is offline   Reply With Quote

Old   March 22, 2006, 10:08
Default Laminar flow seperates easier
  #10
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
Laminar flow seperates easier than turbulent flow. Have you run laminar using a different solver?
eugene is offline   Reply With Quote

Old   March 22, 2006, 14:42
Default Hi! I tried a run with the
  #11
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

I tried a run with the "fixed value" boundary conditions. But there was no difference.

So I did a calculation with 0° angle of attack. Even there the flow "seperates", means there is an subboundary layer with goes the other way. I made a picture to show.



I did a laminar run with simpleFoam and icoFoam.

Andreas
andimb is offline   Reply With Quote

Old   March 23, 2006, 11:31
Default Hey Andreas, Are you sure,
  #12
Senior Member
 
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 8
hartinger is on a distinguished road
Hey Andreas,

Are you sure, that the flow shouldn't seperate.
I mean, why not?

markus
hartinger is offline   Reply With Quote

Old   March 23, 2006, 12:59
Default By different solver, I meant o
  #13
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
By different solver, I meant other than OpenFOAM.

Laminar flow seperates very easily so the laminar seperation results might be correct. For instance, on a cylinder you get seperation quite early before 90 deg if the flow is laminar.

Turbulent seperation is probably wrong. Have you rerun the turbulent case with small fixedValue BCs for k and epsilon?

Also, give the SpalartAllmaras model with nuTilda fixedValue 0 on the wall a try. As far as I know it is one of the better turbulence models for wings.
eugene is offline   Reply With Quote

Old   March 24, 2006, 06:33
Default Hi! I think, I've got the p
  #14
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi!

I think, I've got the problem now.

As I wrote in the start post, I used two definitions of epsilon at the inlet. The first one gave a small and the second (CFX) a great epsilon. As I understood a large epsilon means nearly laminar calculation because all turbulent energy dissipates to velocity. In the OpenFOAM formulation this is just a wrong value of epsilon and that's why the calculation shows a laminar behaviour.

So I have to use the "small" epsilon. There is the problem of the "stagnation point anomaly", which is an overpredicition of turbulent energy at stagnation points in k-epsilon models. The only High-Re models that doesn't have this problem are the "LRR" and the "NonlinearKEShih" turbulence models. They have another formulation of the production term. But as I said and you know, they are High-Re models. There is no low-Re model without this problem.
This overprediction causes the huge boundary layer over the profile and as a consequence of this to less velocity.

But there is still the question why CFX didn't have these problem. They use a scalable wallfunction for some of there turbulence models. So as a consequence they are not as yPlus sensitive as a normal formulations. When I used the standard wallfunction CFX has the same Problems.
And this was the point I wasn't able to see. There boundary contitions are not portable to OpenFOAM, because of the differnce in the wall function theory.

So there is the posibility of the Kato Launder modification of k-epsilon models, that should solve the "stagnation point anomaly". It's an reformulation of the turbulent production.

A laminar run with CFX gave the same results as OpenFOAM. So it's ok.

The fixedValue run gave the same results as the zeroGradient with the right epsilon at the inlet.

So if I don't want to change stuff in the code, I have to live with a too low lift and too high drag. But this is ok as long as I know why this is happening.

As far as I know the Spalart Allmaras model give a good lift prediction but also too high drag (drag is a near wall problem, so you need a good describtion what is happening there). But I'll try, because it is easier and perhaps give better results.

I really want to thank you for your help. This is a great forum with great guys, who want to help.

Please tell me, if I understood something wrong or desripted wrong.

Andreas
owayz likes this.
andimb is offline   Reply With Quote

Old   March 24, 2006, 06:48
Default Well, if it is the stagnation
  #15
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Well, if it is the stagnation point anomaly that's giving you trouble, try using the high-Re RNG k-epsilon: that should do a better job.

If you are set on Kato-Launder, its implementation in OpenFOAM is (almost) trivial: copy k-epsilon, rename the files and class names and change the form of the generation term to use grad U - grad U^T and you're done (it would be worth checking my words in the original reference though).

As for the results by CFX, I would consider it a distinct possibility that they have messed about with the k-epsilon model (to deal with the kind of problem you are seeing) without giving you chapter and verse of what they have done. As you know, CFX have got a very competent guy leading the turbulence work (Dr. Florian Menter) and it would be by no means unusual to do this kind of thing. A correction could be done in several places and you will probably find a brief comment about it somewhere in the documentation if you look closely.

Finally, I can tell you that the k-epsilon implementation in OpenFOAM is strictly "by the book" in all the details including the wall functions, which you can easily check by inspecting the source.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 2, 2006, 04:43
Default Hi Eugene & Andreas, i am m
  #16
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hi Eugene & Andreas,

i am modeling flow over a NACA12 hydrofoil at 5 deg at Re about 10^5 . when i used a k-epsilon i got results similar to klaus wittig ( please search with naca23012 + klaus , for his results - basically the stagnation point cp is larger than 1 and the suction peak is underpredicted ) . when i chose the launderSharma low Reynolds number i got results similar to Andreas.

Prediction of lift is very important in my case since the free surface charecteristics is affected by the suction . so i would like to use the spalartAllmaras turb.model .

Eugene has mentioned that the boundary condition for nuTilda is fixedValue 0 on the wall . isn't the turbulence model used as a low Re model in this case ? is the yplus < approx 5 in this case ? is the physicalType of wall wall & ( not wallFunctions )? If i use the model as a high Re case will the nuTilda be zeroGradient on the wall and the physicalType of wall wallFunctions ? if i know the inlet velocity , what is a quick way of guessing nuTilda at inlet ?

Andreas , did you try out with the spalartAllmaras model and what was your experiance ?

Thanks a lot
Regards
Kumar
kumar2 is offline   Reply With Quote

Old   May 2, 2006, 04:57
Default Hi Eugene & Andreas, i am m
  #17
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hi Eugene & Andreas,

i am modeling flow over a NACA12 hydrofoil at 5 deg at Re about 10^5 . when i used a k-epsilon i got results similar to klaus wittig ( please search with naca23012 + klaus , for his results - basically the stagnation point cp is larger than 1 and the suction peak is underpredicted ) . when i chose the launderSharma low Reynolds number i got results similar to Andreas.

Prediction of lift is very important in my case since the free surface charecteristics is affected by the suction . so i would like to use the spalartAllmaras turb.model .

Eugene has mentioned that the boundary condition for nuTilda is fixedValue 0 on the wall . isn't the turbulence model used as a low Re model in this case ? is the yplus < approx 5 in this case ? is the physicalType of wall wall & ( not wallFunctions )? If i use the model as a high Re case will the nuTilda be zeroGradient on the wall and the physicalType of wall wallFunctions ? if i know the inlet velocity , what is a quick way of guessing nuTilda at inlet ?

Andreas , did you try out with the spalartAllmaras model and what was your experiance ?

Thanks a lot
Regards
Kumar
kumar2 is offline   Reply With Quote

Old   May 3, 2006, 12:21
Default Hi Kumar, the SpalartAllmar
  #18
Member
 
Andreas Hauffe
Join Date: Mar 2009
Location: Dresden, Germany
Posts: 36
Rep Power: 8
andimb is on a distinguished road
Hi Kumar,

the SpalartAllmaras model did a good job. The lift prediction was ok.

As far as I understood the SpalartAllmaras model, it is like a low-Re formulation. The formulation has no real Re-limits.

The boundary conditions are fixedValue 0 at walls and nuTilda=0.1*nu at inlets.

Andreas
andimb is offline   Reply With Quote

Old   May 3, 2006, 16:45
Default Hi Andreas Thanks a lot for
  #19
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hi Andreas

Thanks a lot for the information . i am also going to try out the spalartAllmaras turb.model on my naca12 wing. ( i also started with the standard k-epsilon , then RNG kepsilon , then a low re launderSharma ... and my experiances were very similar to yours )

thanks again for sharing your experiances

regards

kumar
kumar2 is offline   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
turbulence Amir Main CFD Forum 1 January 2, 2008 20:50
turbulence:turbulence viscosity vs kinetic energy sam FLUENT 0 June 25, 2007 23:51
Turbulence AFP Main CFD Forum 2 January 14, 2004 11:05
HELP! TURBULENCE k-e OR k-omega TURBULENCE MODEL? Mirek Kabacinski FLUENT 5 August 24, 2003 22:31
turbulence Burkhard Gölling Main CFD Forum 1 April 7, 2000 15:38


All times are GMT -4. The time now is 22:26.