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

BC's for k-omega SST (low / hi Re setups)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 25, 2012, 18:22
Default BC's for k-omega SST (low / hi Re setups)
  #1
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 150
Rep Power: 8
Jonathan is on a distinguished road
hi everyone,

i have been testing my k-omega SST setup for OF using a simple flat plate geometry, and BC's i have worked out from looking at a number of other posts on the issue. although there is a *lot* out there on the setups you need to use, i dont seem to get the impression that there has been widespread agreement on how the model is supposed to be implemented

anyway, so i have a low Re and a hi Re mesh. the BC's for the two meshes are:

For the hi Re mesh:
Code:
'k'
WALL
{
  type   kqRWallFunction;
  value  uniform 1;
}

'omega'
WALL
{
  type   omegaWallFunction;
  value  uniform 1;
}

'nut'
WALL
{
  type   nutkWallFunction;
  value   uniform 0;
}
and for the low Re mesh:

Code:
'k'
WALL
{
  type   fixedValue;
  value  uniform 1e-12;
}

'omega'
WALL
{
  type   omegaWallFunction;
  value  uniform 1e9;
}

'nut'
WALL
{
  type   nutUSpaldingWallFunction;
  value   uniform 0;
}
Also attached is a plot of the non-dimensional parameters for the wall, using the standard definitions (i.e. not Cmu etc!). As you can see, the low Re result looks ok, but the hi Re BL is too think.

Using yPlusRAS to check the meshes, i see that the first node points for each mesh are in the correct region, so i cant quite work out what the problem is. BC's perhaps??

anyway, i would appreciate any feedback on anything obvious i have gotten wrong + am quite keen to 'finalise' the 'proper' approach for the current implementation of k-omega SST (OF v2.1.1)

thanks a lot in advance!
Jonathan
Attached Images
File Type: jpg bl_comparisons.jpg (32.7 KB, 129 views)
Jonathan is offline   Reply With Quote

Old   September 26, 2012, 08:33
Default
  #2
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 150
Rep Power: 8
Jonathan is on a distinguished road
never mind ... solved it ...
Jonathan is offline   Reply With Quote

Old   September 28, 2012, 16:37
Default
  #3
mrt
New Member
 
Rick Tejeda
Join Date: Sep 2012
Location: Madison, WI
Posts: 7
Rep Power: 4
mrt is on a distinguished road
Hi Jonathan,

How did you resolve this?
mrt is offline   Reply With Quote

Old   October 4, 2012, 06:39
Default
  #4
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 300
Rep Power: 9
openfoam_user is on a distinguished road
Hi Jonathan,

what formula do you use to set 'value uniform 1e9'
for k at the wall for your low Reynolds mesh (y+~1) ?

Best regards,
Stephane.
openfoam_user is offline   Reply With Quote

Old   October 4, 2012, 07:06
Default
  #5
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 150
Rep Power: 8
Jonathan is on a distinguished road
Hi Stephane,

I didnt use a formula - looking through the forums and theory, i realise that omega tends to infinity at the wall - therefore, i just used a large number as the omega BC ...

regards
jonathan
Jonathan is offline   Reply With Quote

Old   February 2, 2014, 19:15
Default
  #6
Member
 
Ye Zhang
Join Date: Dec 2009
Location: Delft,Netherland
Posts: 89
Rep Power: 7
kiddmax is on a distinguished road
Hi Jonathan

I also did this validation with komega SST turbulence model. But with Low Re and high Re model both failed. There is mismatch in the log-law region. Can you please give a hint how you solve it?

Best regards,
Ye
kiddmax is offline   Reply With Quote

Old   February 3, 2014, 04:07
Default
  #7
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 150
Rep Power: 8
Jonathan is on a distinguished road
hi Kid,

I havent looked at OF for a little while, but in the last version i was using (2.1.1), i finally used the following for a low Re mesh:

'k' -> set to SMALL no at wall
'omega' -> use omegaBC built into OF. This calculates omega at the wall using the correct formula from Menter

'nut' -> I used 'calculated' condition and not nut wall function as in my first post

PS I didnt do anymore with the hi-Re meshes, as my project needed low-Re only.

hope this helps a little.

cheers
jonathan

PPS I dont know what the latest versions of OF are using for turbulence wall BC's these days. Might pay to download the latest UG and look through the release notes and code to see what is actually going on ... cheers jon

Quote:
Originally Posted by kiddmax View Post
Hi Jonathan

I also did this validation with komega SST turbulence model. But with Low Re and high Re model both failed. There is mismatch in the log-law region. Can you please give a hint how you solve it?

Best regards,
Ye
Jonathan is offline   Reply With Quote

Old   February 10, 2014, 16:36
Default
  #8
Member
 
Ye Zhang
Join Date: Dec 2009
Location: Delft,Netherland
Posts: 89
Rep Power: 7
kiddmax is on a distinguished road
Hi Jonathan,

Thanks a lot for your reply. I finally got the correct result after several trying.

Best regards,
Ye
kiddmax 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
Low Reynolds Number SST Model Josh CFX 39 February 19, 2014 09:47
overshooting of Omega in SST komega using simpleFoam cm_jubayer OpenFOAM 1 August 26, 2011 13:18
Boundary Conditions for k omega SST dancfd OpenFOAM Pre-Processing 0 June 9, 2011 23:25
Very Low Re: SST vs. v2-f? CAVT Main CFD Forum 0 September 25, 2010 04:34
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


All times are GMT -4. The time now is 16:58.