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

How to solve the boundary layer without wall functions?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 17, 2016, 08:31
Default How to solve the boundary layer without wall functions?
  #1
Member
 
Join Date: Dec 2015
Posts: 40
Rep Power: 3
WhiteW is on a distinguished road
Dear Foamers,
I'm using OpenFOAM 2.3.0. What conditions have I to set on k and omega parameters if I want to solve the boundary layer with k-omega SST model?
I have a mesh with a y+ less than 4 in the wall patches; however, I think the conditions with wall functions are not correct:

omega:
Code:
    "(wall_1|wall_2)"
    {
      type            omegaWallFunction;
      value           $internalField;
     }

k
Code:
   "(wall_1|wall_2)"
    {
		type            kqRWallFunction;
		value   	$internalField;
	}
     }
If I want to solve the boundary layer with a fine mesh, without wall functions, which conditions have I to set?
I tried to use:

omega:
Code:
    "(wall_1|wall_2)"
    {
 	type zeroGradient;
     }

k
Code:
   "(wall_1|wall_2)"
    {
	type zeroGradient;
 }
     }
However when I run the simulation, OF told me:

Code:
Selecting RAS turbulence model kOmegaSST
--> Upgrading k to employ run-time selectable wall functions
    Backup original k to k.old
    Writing updated k
--> Upgrading omega to employ run-time selectable wall functions
    Backup original omega to omega.old
    Writing updated omega
--> Creating mut to employ run-time selectable wall functions
    Writing new mut
--> Creating alphat to employ run-time selectable wall functions
    Writing new alphat
So OF set the wallfunctions again.
Witch is the correct setting?
Thanks,
WhiteW
WhiteW is offline   Reply With Quote

Old   March 20, 2016, 06:31
Default
  #2
Member
 
Join Date: Dec 2015
Posts: 40
Rep Power: 3
WhiteW is on a distinguished road
Any idea?
WhiteW
WhiteW is offline   Reply With Quote

Old   March 21, 2016, 04:34
Default
  #3
Member
 
Alex
Join Date: Jan 2014
Posts: 94
Rep Power: 4
Kina is on a distinguished road
Well you would not set k and omega to zeroGradient, but to a very small, non zero fixedValue. For example:

k: fixedValue, value uniform 1e-8

However, the omega wall function is known to work well with lowRe meshes as it only applies the correct omega value to the wall without much of damping or blending. So i'd suggest just switching k to a fixedValue.
Kina is offline   Reply With Quote

Old   March 21, 2016, 06:41
Default
  #4
Member
 
Join Date: Dec 2015
Posts: 40
Rep Power: 3
WhiteW is on a distinguished road
Thanks for the reply. Now I left omega with the wallFunction and I set k as:

Code:
    "(wall_1|wall_2)"
    {
        type           fixedValue;
	value          uniform 1e-8;
    }
However, when the solver (rhosimplecFoam) run, it still changes the conditions using the wallFunctions with the value I have set:

Code:
Selecting RAS turbulence model kOmegaSST
--> Upgrading k to employ run-time selectable wall functions
    Backup original k to k.old
    Writing updated k
--> Upgrading omega to employ run-time selectable wall functions
    Backup original omega to omega.old
    Writing updated omega
--> Creating mut to employ run-time selectable wall functions
    Writing new mut
--> Creating alphat to employ run-time selectable wall functions
    Writing new alphat
the new k is:

Code:
    "(wall_1|wall_2)"
    {
        type            compressible::kqRWallFunction;
        value           uniform 1e-08;
    }
is it correct?
WhiteW 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
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 09:44
Wind turbine simulation Saturn CFX 45 February 8, 2016 05:42
Prismatic boundary layer KateEisenhower enGrid 5 September 15, 2015 07:48
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32


All times are GMT -4. The time now is 19:41.