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

nutWallFunction

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

Like Tree24Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   March 16, 2011, 01:57
Default nutWallFunction
  #1
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear All
I am using k − ω-SST in which i am unable to understand the wall functions i.e. nutWallFunction,nutRoughWallFunction etc.Will anybody tell me what it is,How to calculate and for what it is.
Thanks In Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   March 17, 2011, 03:56
Default
  #2
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear All
Got the answer thanks to all.
__________________
Simulation Is Determination Of Imagination Towards Approximation


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   March 18, 2011, 07:52
Default
  #3
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 8
idrama is on a distinguished road
Hey kurne,

wall functions are used to model the velocity profile at the first node to the wall. Model based on magnitudes such as k, epsilon, nut, wall shear stress, friction velocity and so on. Therefore, the implemented function, e.g. epsilonWallFunction, provide the computation of the law of the wall.
Instead using them, you need to resolve the boundary layer such that the transition from laminar to turbulent can be captured.

Cheers
idrama is offline   Reply With Quote

Old   March 18, 2011, 08:05
Default
  #4
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Claus Meister

Thank you very much for the answer, in the real sense i understood the meaning of wall function but i stuck up with the nut wall function but after going through the forum and upon your reply i got the correct concept.

Well, want to ask you one silly question do you know how to calculate the nutWallFunction at the inlet if we are having the data at inlet say velocity,turbulence intensity,turbulence kinetic energy and turbulence energy dissipation rate.If you know then please tell me.But i got the information from the forum that initial value of 0 is generally used.

Thanks Once Again.
__________________
Simulation Is Determination Of Imagination Towards Approximation


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   March 18, 2011, 09:06
Default
  #5
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 8
idrama is on a distinguished road
The wall function are only used at the wall and not at a inlet or outlet patch. Since the nut value is composed of k and epsilon, you should set a "calculated" BC at the inlet and the outlet, respectively. A "calculated" BC means that the values at the respective patch are computed by other magnitudes.

Cheers
idrama is offline   Reply With Quote

Old   March 19, 2011, 00:05
Default
  #6
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Claus Yes thats right but when i go through the tutorials of OpenFOAM 1.7.1 of mixer in simpleSRFFoam solver. I found the nut directory in 0 file.When i remove the nut file from the 0 directory the solver wont run.
It means that the value of nut is required at 0 directory with the value at inlet,outlet and other patches too,Correct me if i am wrong.So i asked that how to calculate the value of nut at inlet and outlet.
Thanks In Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation


Best Regards

Mubeen K Kurne

Last edited by kurne; March 21, 2011 at 01:34.
kurne is offline   Reply With Quote

Old   March 20, 2011, 19:15
Default
  #7
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 8
idrama is on a distinguished road
You dont have to calculate the nut values for inlet and outlet. You just have to set:

inlet
{
type calculated;
value uniform 0;
}

wall
{
type nutWallFunction;
value uniform 0;
}

etc.

hope that helps
idrama is offline   Reply With Quote

Old   March 21, 2011, 01:23
Default
  #8
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Claus

Thank you very much for answering to my query.Well, is it possible to calculate the value of nutWallFunction at the patches other than inlet and outlet.If yes will you please tell me how.

I know that my queries are silly hence sorry for it.
Thanks In Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   October 28, 2011, 06:01
Default
  #9
Senior Member
 
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 152
Rep Power: 6
Anne Lincke is on a distinguished road
Quote:
Originally Posted by kurne View Post
Dear All
Got the answer thanks to all.

Hey Kurne,

I am investigating about the nutWallFunction and I do not really understand the formula used in OpenFOAM in nutWallFunction which reads

nutw[faceI] = nuw[faceI]*(yPlus*kappa_/log(E_*yPlus) - 1.0);

Can you please explain your findings? I found this question in some threads so I think me and other users will be very grateful.

Thank you in advance.
Anne Lincke is offline   Reply With Quote

Old   October 28, 2011, 07:06
Default
  #10
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 8
idrama is on a distinguished road
Loop up in Versteeg pp. 275ff; there is a description about the entire procedure. I used by my self yesterday.

http://books.google.de/books?id=RvBZ...page&q&f=false

Cheers
idrama is offline   Reply With Quote

Old   November 2, 2011, 04:00
Default
  #11
Senior Member
 
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 152
Rep Power: 6
Anne Lincke is on a distinguished road
Thank you very much. I borrowed the book and will post my findings here as soon as possible.
Anne Lincke is offline   Reply With Quote

Old   November 25, 2011, 23:26
Default nutkwallfunction
  #12
Member
 
pooyan
Join Date: Nov 2011
Posts: 62
Rep Power: 5
sam1364 is on a distinguished road
Hello everyone

I have recently started working with OpenFoam. I am now solving cavity problem in 3D. for turbulence model, I selected RANS and k-epsilon. for wall function, I selected nutkwallfunction. However, I could not find type of wall functions (how the wall function is defined) anywhere except from the c++ documentary which was not clear to me. I would be grateful if you could tell me what specific wall functions are used in nutkwallfunction.

Thanks,
sam1364 is offline   Reply With Quote

Old   November 28, 2011, 04:18
Default
  #13
Senior Member
 
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 152
Rep Power: 6
Anne Lincke is on a distinguished road
You can have a look at

~/OpenFOAM/OpenFOAM-1.7.0/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions>

There you will find all kind of wall functions.
Especially for your wall function it is written in the .H file

Class
Foam::incompressible::RASModels::nutkWallFunctionF vPatchScalarField

Description
Boundary condition for turbulent (kinematic) viscosity when using wall
functions, based on turbulence kinetic energy.
- replicates OpenFOAM v1.5 (and earlier) behaviour

At first glance I could not see any differences between nutWallFunction and nutkWallFunction....
sam1364 likes this.
Anne Lincke is offline   Reply With Quote

Old   November 28, 2011, 04:20
Default
  #14
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Quote:
Originally Posted by sam1364 View Post
Hello everyone

I have recently started working with OpenFoam. I am now solving cavity problem in 3D. for turbulence model, I selected RANS and k-epsilon. for wall function, I selected nutkwallfunction. However, I could not find type of wall functions (how the wall function is defined) anywhere except from the c++ documentary which was not clear to me. I would be grateful if you could tell me what specific wall functions are used in nutkwallfunction.

Thanks,
the boundary condition nutkWallFunction sets the turbulent viscosity in the first node point based on the logarithmic law of the wall http://en.wikipedia.org/wiki/Law_of_the_wall . the difference between the nutUWallFunction and the nutkWallFunction lies in the way yPlus is calculated (in the first case based on the velocity close to the wall and in the latter case based on the turbulent kinetic energy close to the wall)
kiddmax, fumiya, sam1364 and 4 others like this.
__________________
~roman
romant is offline   Reply With Quote

Old   November 28, 2011, 04:23
Default
  #15
Senior Member
 
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 152
Rep Power: 6
Anne Lincke is on a distinguished road
What version do you use?
I cannot find nutUWallFunction at OF-1.7.

When looking at the source code of nutWallFunction I see the same computation of yPlus like in nutkWallFunction, based on turbulent kinetic energy....
Anne Lincke is offline   Reply With Quote

Old   November 28, 2011, 04:25
Default
  #16
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Quote:
Originally Posted by Anne Lincke View Post
What version do you use?
I cannot find nutUWallFunction at OF-1.7.

When looking at the source code of nutWallFunction I see the same computation of yPlus like in nutkWallFunction, based on turbulent kinetic energy....
In 2.0 more nutWallFunction boundary conditions were introduced, one of them nutUWallFunction, where basically the way of how yPlus is calculated is changed.
__________________
~roman
romant is offline   Reply With Quote

Old   November 28, 2011, 04:28
Default
  #17
Senior Member
 
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 152
Rep Power: 6
Anne Lincke is on a distinguished road
Okay, thank you for this information!

Can you also explain the computation of nut at the wall function, especially this formula

nutw[faceI] = nuw[faceI]*(yPlus*kappa_/log(E_*yPlus) - 1.0);

Is it derived from the logLaw? How?
Anne Lincke is offline   Reply With Quote

Old   November 28, 2011, 05:36
Default log law in compact version
  #18
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Hej,

I just found how this is derived and calculated it myself, http://www.cfd-online.com/Wiki/Two_e...wall_functions

if you see how the rho is eliminated, maybe you can post it here.

solution to eliminate rho:
Code:
(mut+mu)Up/yP = (rho y uTau kappa )/(log(yPlus E))
=1/rho (mut+mu)Up/yP = ( y uTau kappa )/(log(yPlus E))
= (nut+nu)Up/yP = ( y uTau kappa )/(log(yPlus E))
and then from there by using the definition of yPlus, one obtains equation (2)
__________________
~roman

Last edited by romant; November 28, 2011 at 05:50. Reason: found solution to eliminate rho
romant is offline   Reply With Quote

Old   November 28, 2011, 06:40
Default
  #19
Senior Member
 
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 152
Rep Power: 6
Anne Lincke is on a distinguished road
Thank you so much!
Now I understand the derivation.
There is one typo on the site, the missing density at the definition of tau_w,
But there is no problem as the whole equation

rho (nut + nu) * up/yp = (rho *u_tau *kappa *u) / ln (E y+)

can simply be divided by rho.
Anne Lincke is offline   Reply With Quote

Old   November 28, 2011, 06:50
Default
  #20
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
Quote:
Originally Posted by Anne Lincke View Post
Thank you so much!
Now I understand the derivation.
There is one typo on the site, the missing density at the definition of tau_w,
But there is no problem as the whole equation

rho (nut + nu) * up/yp = (rho *u_tau *kappa *u) / ln (E y+)

can simply be divided by rho.
I changed it on the wiki page, it should be OK now
Anne Lincke and aerothermal like this.
__________________
~roman
romant 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
second order schemes marine OpenFOAM 66 October 26, 2014 12:10
nutWallFunction questions madad2005 OpenFOAM 9 October 28, 2011 06:36
nutWallFunction vs. nutkWallFunction romant OpenFOAM Running, Solving & CFD 1 October 21, 2011 10:12
A few questions about Spalart Allmaras model and y seb62 OpenFOAM Running, Solving & CFD 8 January 15, 2011 14:09
nutWallFunction boundary condition AirS OpenFOAM Pre-Processing 6 March 3, 2010 08:31


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