CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Why Menter's SST model low-Re issue has not been seriously investigated?

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

Like Tree66Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 30, 2013, 11:16
Default
  #41
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: India
Posts: 205
Rep Power: 18
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Sure, I am doing it.

I have a small doubt, for epsilon do I need to use epsilonWallFunction??

top0 (wall)
{
type epsilonWallFunction;
value uniform 0.26295;

}
sivakumar is offline   Reply With Quote

Old   July 30, 2013, 11:18
Default
  #42
Senior Member
 
Joachim
Join Date: Mar 2012
Location: Paris, France
Posts: 145
Rep Power: 15
Joachim is on a distinguished road
If you have y+~1, you can resolve the full boundary layer, so that you don't need to use any wall function.

Simply put k=1e-9 and epsilon = 1e-9 at the wall and you should be good!
Joachim is offline   Reply With Quote

Old   July 30, 2013, 11:21
Default
  #43
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: India
Posts: 205
Rep Power: 18
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
the key word should be fixedValue right?

top0 (wall)

type fixedValue;
value uniform 1e-9;
sivakumar is offline   Reply With Quote

Old   July 30, 2013, 11:22
Default
  #44
Senior Member
 
Joachim
Join Date: Mar 2012
Location: Paris, France
Posts: 145
Rep Power: 15
Joachim is on a distinguished road
yep, that's right.
Joachim is offline   Reply With Quote

Old   July 30, 2013, 14:47
Default
  #45
Senior Member
 
Joachim
Join Date: Mar 2012
Location: Paris, France
Posts: 145
Rep Power: 15
Joachim is on a distinguished road
Out of curiosity, I tried to use the k-omega model on some basic test cases. The version implemented in OpenFOAM is the original Wilcox 1988 model. Impossible to get a solution!

Has someone already faced this problem? I checked in the .C and .H, but could not find any mistake in the implementation so far...
Joachim is offline   Reply With Quote

Old   July 30, 2013, 15:09
Default
  #46
Senior Member
 
Joachim
Join Date: Mar 2012
Location: Paris, France
Posts: 145
Rep Power: 15
Joachim is on a distinguished road
Actually, it seems to work. I just needed to be more careful with the initial and boundary conditions...
Joachim is offline   Reply With Quote

Old   August 1, 2013, 06:45
Default
  #47
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: India
Posts: 205
Rep Power: 18
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Dear All,
As Joachim and others instruction, I have tried both LaunderSharmaKE and kOmegaSST for low re calculation. but none of them working fine.

the Boundary conditions which i used are

For kEpsilon model:

0/k:

internalField uniform 0.39;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.39;
}
wall
{
type fixedValue;
value uniform 1e-12;

0/e:

internalField uniform 0.26295;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.26295;
}

wall
{
type fixedValue;
value uniform 1e-12;
}

0/nut:

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}

wall
{
type nutLowReWallFunction;
value uniform 0;
}

for kOmegaSST:

0/k:

internalField uniform 0.39;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.39;
}
wall
{
type fixedValue;
value uniform 1e-12;


0/Omega:

internalField uniform 3.7;

boundaryField
{
inlet
{
type fixedValue;
value uniform 3.7;
}

wall
{
type omegaWallFunction;
value uniform 1441;
}

0/nut:

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}

wall
{
type nutkWallFunction;
value uniform 0;
}

it seems for low re kepsilon the boundary conditions are fine but the simulation crashes due to floating point, the same error for kOmegaSST.

I tried different BC as well but none of them working fine.

can you guys spend few minutes to check my BCs to find what is wrong?

Thanks,
Sivakumar
sivakumar is offline   Reply With Quote

Old   August 1, 2013, 10:48
Default
  #48
New Member
 
Andre Weiner
Join Date: Aug 2012
Posts: 29
Rep Power: 13
andre.weiner is on a distinguished road
Hello Sivakumar!

I did some low Re simulations with k-omega-SST model and it worked fine for me. Your bc's look good to me (only refering to k-omega-SST), so there may be another problem with your simulation.
For k you could also try a zero gradient bc. For nut you could try to set the value to zero (like in the nutLowReWallFunction) or just don't correct it (calculated).
Maybe your freestream bc's are not appropriate but thats hard to say without knowing the characteristics of the flow.

Regards, Andre
andre.weiner is offline   Reply With Quote

Old   August 1, 2013, 12:28
Default
  #49
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: India
Posts: 205
Rep Power: 18
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Hi Andre,
Thanks for your reply, I believe that the problem could be mesh.
however checkMesh is not giving any warning and error message.

Thanks,
Sivakumar
sivakumar is offline   Reply With Quote

Old   October 22, 2013, 07:18
Default
  #50
Senior Member
 
MadsR's Avatar
 
Mads Reck
Join Date: Aug 2009
Location: Copenhagen, Denmark
Posts: 177
Rep Power: 17
MadsR is on a distinguished road
Sivakumar, did you resolve the problem? I have the same issue. I have two meshes which, in essence, are identical but one is fine with around 2000k cells and the other is coarse with around 500k cells. The fine mesh spends most of the extra cells near the walls and so the y+ min/max/avg are significantly lower.

OpenFOAM 2.2.2 blows up on the fine mesh but solves perfectly on the coarse one. I've been changing nut/k/omega BCs to whatever is suggested here on cfd-online without luck. I will continue to investigate on the mesh distribution close to the walls since some say that smoothness is important for OF. My mesh is made from blockMesh so it is already quite smooth, but there is certainly stretching near walls. No checkMesh issues.

Mads
__________________
Online free airfoil-mesher for OpenFOAM here
MadsR is offline   Reply With Quote

Old   October 22, 2013, 09:09
Default
  #51
Senior Member
 
sivakumar selvaraju
Join Date: Mar 2009
Location: India
Posts: 205
Rep Power: 18
sivakumar is on a distinguished road
Send a message via Skype™ to sivakumar
Dear Mads,
Yes, I have solved that problem, the cases are running fine.
I have few questions, which turbulence model you want to use? what is the yplus value for the fine mesh? , I think KOSST will be fine with all kind of y+ values, but for LaunderSharmaKE we need be little bit careful.

my problem was yplus value, it was less than 1(people usually says y+ less than or equal to 1 will be fine, but in my experience not), that is why my case was diverging, now I changed the y+ value to in between 3 to 5 its working fine.

I hope you will be fine with BC's, if you need any info regarding that please go through the below link.

http://www.cfd-online.com/Forums/ope...ase-setup.html


Thanks,
Sivakumar
sivakumar is offline   Reply With Quote

Old   October 22, 2013, 10:17
Default
  #52
Senior Member
 
MadsR's Avatar
 
Mads Reck
Join Date: Aug 2009
Location: Copenhagen, Denmark
Posts: 177
Rep Power: 17
MadsR is on a distinguished road
Dear Sivakumar,

thank you for your reply. My geometry is quite complex so y+ varies quite a lot from one area to another. I try to get as low on y+ as possible because I know that I need low-Re treatment in certain areas but in other areas I save cells and get a larger y+ (maybe even into the hundreds).
A code like CFX copes nicely with this using their automatic wall functions but I am not sure how this works out in OpenFOAM - it crashes for now It seems that it should be ok to use omegaWallFunction as it "automatically" adopts to local, sensible, y+ whereas I am more uncertain about the wall treatment for k. I am not going to use LaunderSharmaKE.

I think I will try to remesh so I don't have y+<1 anywhere which seems to have been your problem.

Mads
__________________
Online free airfoil-mesher for OpenFOAM here
MadsR is offline   Reply With Quote

Old   April 4, 2014, 06:07
Default
  #53
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Joachim,

I ran the case in Fluent and get the following:
k-omega-sst-tau395.png

And also compared to Spalart-Allmaras:
tau_w.jpg
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   May 9, 2014, 11:17
Default
  #54
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Jonathan, I don't know if you still use it, but the function you call "gammaStar0" is not correct. At least not compared to Fluent.
In Fluent they don't use the "SST" beta_i but the regular "k-Omega" beta_i for the "beta_i / 3.0" calculation. So you neet to change it to

return 0.072 / 3.0;

instead of

return betai(F1) / 3.0;
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   May 13, 2014, 04:31
Default
  #55
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Please find a working version of the incompressible sst model with low-Re wall damping functions here:
http://www.cfd-online.com/Forums/ope...ng-fluent.html
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   May 14, 2014, 07:10
Default
  #56
Senior Member
 
Join Date: Mar 2010
Posts: 172
Rep Power: 17
Jonathan is on a distinguished road
Hi Phillip,

thanks for the post + correction. I havent had time to check the correction you made, but surely it is correct - I didnt get anyone willing to help debug first time around

when i get a chance i will look and see why i wrote it like i did!

cheers and regards
jonathan
Jonathan is offline   Reply With Quote

Old   May 14, 2014, 07:14
Default
  #57
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
It's not your fault. It's the Fluent manual that is misleading.
There is also a second typo in the Fluent manual in the omega production term, that I corrected in my code.
__________________
The skeleton ran out of shampoo in the shower.

Last edited by RodriguezFatz; May 14, 2014 at 09:18.
RodriguezFatz is offline   Reply With Quote

Old   June 12, 2014, 22:21
Default Problems with validating kOmegaSST for channel flow
  #58
Member
 
Manan
Join Date: Oct 2013
Location: Göteborg
Posts: 37
Rep Power: 12
MaLa is on a distinguished road
Hi Vesselin and Joachim

While using OpenFOAM 2.2x to implement two dimensional developed, channel flow with kOmegaSST, I have encountered the same results as Joachim and have found it quite difficult to figure out the issue.

Here's a post where I talk about it.

http://www.cfd-online.com/Forums/ope...nnel-flow.html

Maybe I should try the same on OF 1.7x as suggested by Vesselin.
MaLa is offline   Reply With Quote

Old   January 8, 2018, 16:20
Default
  #59
Senior Member
 
Reviewer #2
Join Date: Jul 2015
Location: Knoxville, TN
Posts: 141
Rep Power: 10
randolph is on a distinguished road
Quote:
Originally Posted by Jonathan View Post
Hi V,

Ok, I have done all the homework on this and thought perhaps i should jot them down for anyone else wanting to know this stuff

1) Regarding the implicit "low Re" nature of the k-omega models + variants - yes, you are spot on. I knew they were able to be integrated to the wall w/out any modifications, but thought that damping functions were needed for the production terms etc and to extract the proper / full low Re behaviour. So - as you said, they are not and seems they were added specifically by Wilcox for trying to model laminar-to-turbulent BL transition (Wilcox 1994), so even the Fluent terminology is misleading on this.

One point of interest thought, Menter 1994 says "one criticism of the k-omega model [is that it] ... does not correctly predict the asymptotic behaviour of the turbulence as it approaches the wall. However ... even if the turbulence model is not asymptotically consistent, the mean flow profile and the wall skin friction are still predicted correctly. A second point [is that] ... the k-omega model does not accurately represent the k and epsilon distribution in agreement with DNS data ... In cases where the agreement with DNS data is considered important, the damping functions developed by Wilcox can be applied to the present [k-omega SST] model."
Menter, F 1994 "Two-equation eddy-viscosity turbulence models for engineering applications", AIAA Journal, 32(8) 1994.

However, as per Menter's comments, I think this is a relatively 'minor' detail unless you are actually wanting the turbulence profiles exactly from your simulation.

2) Correct BC's - As you say, the correct BC's for the k-omega model are:
i)k=0 or some small value
ii)omega=omegaWallFunction or Menter's omega BC
iii)nut=calculated

3) Regarding adding the damping functions - I am still deciding whether this will be worthwhile, but I’ll put the model up if i do. I think it might be a nice addition to have anyway ...

Again, thanks for the discussion,
best regards
jonathan

PS) Finally, I was just wondering whether you could explain what you think EXACTLY Henry was trying to do but suggesting the use of nutUSpaldingWallFunction for nut in the 'fix' for kOmegaSST. The nut values are calculated from the value of the turbulence properties - what exactly would using nutUSpaldingWallFunction do - override the way in which nut was calculated at the walls? or something else. I am a little lost as to exactly the point of specifying nut (for any of the t-models i guess), apart from an initial guess for the first solver iteration.
___________________________________________

Regarding the (1) question, Menter probably is not the best source of answer. Please check the original paper by Wilcox (1991).

Even for standard K-omega (wilcox), simple modifications (this is not what people usually called low-re correction) are still needed in order to integrate through viscous sub layer. By doing so, satisfactory prediction for ski friction and mean flow profile (But not the turbulence characteristic) can be obtain with standard K-omega formula. (This is not the case for Epsilon based model. Epsilon based model require extensive modification for just getting the mean flow profile right).

While, in some applications where the correct prediction of turbulence characteristic (such as TKE peak near the wall)
are needed in addition to just the mean flow profile, standard k-omega (simply modified) will require additional low-Reynolds modification (people usually call this low-Re correction for SST and K-omega) in order to get the right turbulence characteristic.

If your simulation does not require detailed information over turbulence (TKE ,Omega ...), standard k-omega SST will do it. If you need to predict TKE peak or something, you need to include the low-Reynolds corrections.
lev likes this.
randolph is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX-11 vs CFX-13 SST Model Zigainer CFX 10 December 2, 2011 05:40
Low Reynolds k-epsilon model YJZ ANSYS 1 August 20, 2010 14:57
Understanding k-omega SST model source code tmhonka OpenFOAM Programming & Development 1 September 8, 2009 08:33
multi fluid mixture model issue rystokes CFX 3 August 9, 2009 20:13
Convergence issue in SST for Porous model Raj CFX 0 May 2, 2008 03:43


All times are GMT -4. The time now is 05:24.