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

nut and nuTilda

Register Blogs Community New Posts Updated Threads Search

Like Tree41Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 27, 2010, 15:51
Angry nut and nuTilda
  #1
New Member
 
Join Date: Jan 2010
Posts: 1
Rep Power: 0
phile is on a distinguished road
Dear Foamers,

I am new to OpenFOAM, so my question may be silly. What's the difference between nut and nuTilda? Isn't nut for kinematic turbulence viscosity? But what for nuTilda?

Thanks,
Phile
phile is offline   Reply With Quote

Old   January 29, 2010, 00:43
Default
  #2
Member
 
Jinbiao Xiong
Join Date: Oct 2009
Location: China/Japan
Posts: 50
Rep Power: 16
JinBiao is on a distinguished road
Send a message via MSN to JinBiao
Quote:
Originally Posted by phile View Post
Dear Foamers,

I am new to OpenFOAM, so my question may be silly. What's the difference between nut and nuTilda? Isn't nut for kinematic turbulence viscosity? But what for nuTilda?

Thanks,
Phile
Hi,phile

I am also confused about this. And I have checked the OpenFOAM Programmer's C++ Reference, which you can find in the website of opencfd ltd. It seems nuTilda is a parameter for LES simulation. However, I donot know much about LES. If you can please check it and tell me about it if possible.

Jinbiao
JinBiao is offline   Reply With Quote

Old   January 29, 2010, 02:25
Default
  #3
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
perhaps it might help if you check in src/turbulenceModels/incompressible/RAS/SpalartAllmaras

and scroll to the bottom of the .C file
nut_.internalField() = fv1*nuTilda_.internalField();
niklas is offline   Reply With Quote

Old   January 29, 2010, 03:07
Default
  #4
Member
 
Jinbiao Xiong
Join Date: Oct 2009
Location: China/Japan
Posts: 50
Rep Power: 16
JinBiao is on a distinguished road
Send a message via MSN to JinBiao
Thanks. Niklas. It does help a lot.

It is great to find out the answer. nuTilda should b
e Spalart-Allmaras variable.
So it belongs to the
Spalart-Allmaras model which is one of the RAS model.
That is to say, when I am using the k-e model, I donot need it.

Am I right?

Jinbiao
nanavati, pela145, granzer and 3 others like this.

Last edited by JinBiao; January 29, 2010 at 03:24.
JinBiao is offline   Reply With Quote

Old   January 29, 2010, 03:44
Default
  #5
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 29
niklas will become famous soon enoughniklas will become famous soon enough
correct
niklas is offline   Reply With Quote

Old   February 1, 2010, 11:32
Default
  #6
New Member
 
Anonymous
Join Date: Jan 2010
Posts: 9
Rep Power: 16
Philer is on a distinguished road
I was really confused since the nuTilda file exists even with kEpsilon models in the tutorial. Now I know it is superfluous.

Thank you, guys!

Phile
Philer is offline   Reply With Quote

Old   September 30, 2010, 10:48
Unhappy nut is needed after V1.6
  #7
New Member
 
Jie (Jay) Zhang
Join Date: Sep 2010
Location: Phoenix, AZ, U.S.
Posts: 28
Rep Power: 15
jiez is on a distinguished road
nut (turb viscosity) is needed after V1.6
jiez is offline   Reply With Quote

Old   January 30, 2011, 19:25
Default
  #8
DLC
Member
 
Davide Lupo Conti
Join Date: Nov 2009
Posts: 34
Rep Power: 16
DLC is on a distinguished road
hello everyone!
can someone tell me how to initialize nut and nuTilda? I always used kOmegaSST and for those I found the formulas to initialize k and omega, but I couldn't find out how to initialize variables for spalartAlllmaras..
Thank you!
Endrew Barreto and watermelon like this.
DLC is offline   Reply With Quote

Old   January 31, 2011, 09:29
Default
  #9
Senior Member
 
Travis Carrigan
Join Date: Jul 2010
Location: Arlington, TX
Posts: 161
Rep Power: 15
tcarrigan is on a distinguished road
If I understand correctly, If you'd like to initialize and record nut and nuTilda when using kOmegaSST, simply place them in the 0 directory and set them to 0...or you could calculate them based on k and omega.

Because nut and nuTilda aren't required when using kOmegaSST (only k and omega are required), by placing the nut and nuTilda dictionary files in the 0 directory, you're forcing OpenFOAM to record those values at each iteration. This is good for postprocessing as you don't have to calculate them yourself after the simulation is over.

Note: I believe that nut and nuTilda are only required for the SpalartAllmaras turbulence model.

Hope this helps.
tcarrigan is offline   Reply With Quote

Old   January 31, 2011, 18:55
Default
  #10
DLC
Member
 
Davide Lupo Conti
Join Date: Nov 2009
Posts: 34
Rep Power: 16
DLC is on a distinguished road
Hi Trevis,
thanks for the quick reply.. I realized now that the post was written in an awful English...(late night...)

I well know that nut and nuTilda are not requested when using kOmegaSTT...
Now I am testing spalartAllmaras model and I don't know which values assign in the 0/nut and 0/nuTilda...
(i'm using freestream bc everywhere except for on the foil, where I'm using wallfunctions for nut and fixedvalue for nuTilde )...

Thaks!
DLC is offline   Reply With Quote

Old   February 9, 2011, 09:55
Default
  #11
Member
 
Arina
Join Date: Oct 2009
Location: Belarus
Posts: 71
Rep Power: 16
Akuji is on a distinguished road
Send a message via ICQ to Akuji
Hello!
I have a question... I use Spalart-Allmaras model, so I need nut and nuTilda. As I understand, nut is turbulent viscosity and nuTilda is kinematic turbulenr viscosity. Am I right?
And how can I calculate them?
Endrew Barreto and Gang Wang like this.
Akuji is offline   Reply With Quote

Old   February 16, 2011, 03:11
Default
  #12
Member
 
Robert Ong
Join Date: Aug 2010
Posts: 86
Rep Power: 15
rob3rt 0ng is on a distinguished road
Quote:
Originally Posted by niklas View Post
perhaps it might help if you check in src/turbulenceModels/incompressible/RAS/SpalartAllmaras

and scroll to the bottom of the .C file
nut_.internalField() = fv1*nuTilda_.internalField();
Hi Nicklas or others,

What is the appropriate value for fv1? some say 0.1 the others say 0.5??

Thanks and regards,
Robert
rob3rt 0ng is offline   Reply With Quote

Old   February 16, 2011, 03:42
Default
  #13
Member
 
Robert Ong
Join Date: Aug 2010
Posts: 86
Rep Power: 15
rob3rt 0ng is on a distinguished road
Quote:
Originally Posted by rob3rt 0ng View Post
Hi Nicklas or others,

What is the appropriate value for fv1? some say 0.1 the others say 0.5??

Thanks and regards,
Robert
I've already got the answer. Sorry for bothering you.

Thanks and regards,
Robert
rob3rt 0ng is offline   Reply With Quote

Old   August 6, 2012, 16:08
Default
  #14
New Member
 
Rafael Valenzuela Musura
Join Date: Feb 2012
Posts: 27
Rep Power: 14
rafamusura is on a distinguished road
rob3rt 0ng

Can you answer the question you've asked please I don't know what value should I use!

thanks!
rafamusura is offline   Reply With Quote

Old   August 6, 2012, 23:49
Default
  #15
Member
 
Robert Ong
Join Date: Aug 2010
Posts: 86
Rep Power: 15
rob3rt 0ng is on a distinguished road
Hi Rafael,

From Turbulence modelling for CFD book by D. Wilcox, fv1 is defined as x^3/(x^3 + cv1^3), where x is nut/nu and cv1 is a closure coefficient of 7.1. In my case, fv1 adds up to 1.00.

Regards,
Robert
Gang Wang and mostafa kareem like this.
rob3rt 0ng is offline   Reply With Quote

Old   August 7, 2012, 21:09
Default
  #16
New Member
 
Rafael Valenzuela Musura
Join Date: Feb 2012
Posts: 27
Rep Power: 14
rafamusura is on a distinguished road
Yes but the X is defines as nutilda/nu so i'm confused. Nut is dependent only for the "arbitrary" value of nutilda. In other thread I've read that nutilda=3*nut=nu*0.1 is that correct?

Thanks
rafamusura is offline   Reply With Quote

Old   August 7, 2012, 23:25
Default
  #17
Member
 
Robert Ong
Join Date: Aug 2010
Posts: 86
Rep Power: 15
rob3rt 0ng is on a distinguished road
x is defined as nut/nu. Nut is turbulent kinematic viscosity which can be readily work out based on your free-stream turbulence parameters.

Regards,
Robert
Alisa_W likes this.
rob3rt 0ng is offline   Reply With Quote

Old   August 8, 2012, 02:56
Default
  #18
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

just one question.
The SpalartAllmaras model is just a analytical description to get the turbulent length scale and the kinetic energy of the turbulence to calculate mu_t. I refer to Ferziger and Peric.

So whats the advantage of that model instead of using a two differential equation model in RANS? In my human mind you should get the better results with a two equation model instead of using one equation or just analytical description for that, or not? Maybe in the past computational resources were very expensive and it was better to use
a none, or one equation model for that but today ?


Tobi
nimasam likes this.

Last edited by Tobi; August 8, 2012 at 03:31.
Tobi is offline   Reply With Quote

Old   August 8, 2012, 08:00
Default
  #19
New Member
 
Rafael Valenzuela Musura
Join Date: Feb 2012
Posts: 27
Rep Power: 14
rafamusura is on a distinguished road
Robert Ong
I'm not shure that's the way it work, because you are defining nut=fv1*nutilda, fv1=x^3/(x^3 + cv1^3), and like you say X is a function of nut, if so you're defining nut as a function of nut, nut=nutida*f(nut).

in the wiki page http://www.cfd-online.com/Wiki/Spalart-Allmaras_model, X is defined as nutlinda/nu. So what's the rigth interpretation?.

Tobi, I have the same question, I was trying to get a k-e simulation but it always crashed, so I took a step backward to a more simple model first. If you have the answer for that just let me know.
rafamusura is offline   Reply With Quote

Old   August 8, 2012, 08:14
Default
  #20
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by rafamusura View Post

Tobi, I have the same question, I was trying to get a k-e simulation but it always crashed, so I took a step backward to a more simple model first. If you have the answer for that just let me know.

Hey rafamusura,
well I am writing my masterthesis at the moment and getting deeper in the theory. All your navier stokes equation are in differential equations so its logical and makes sence to use differential equation for calculating the turbulent viscosity. So in the simple methods - none, one, or two equation -, think you should prefer two equation models, especially if you are using
complex geometry or better reynolds-models, couse of the modelling of the stress tensor refering to Ferziger and Peric. But you have do decide how accurate your simulation should be and not every model is working good for all geometry`s and problems. Its possible that a solution of reynolds-models compared to two equation models is the same for the first geometry and could be very different for another geometry. I am not a pro in turbulence modelling but for RANS and the viscosity models i think you should use two equation models for closure the navier stokes equations.

And your k-eps. problem is a BC problem or a problem with your initialization. But that `s a other topic

If I am wrong, please correct me!
Blanco and openfoam_aero like this.
Tobi is offline   Reply With Quote

Reply

Tags
nut nutilda


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
SpalartAllmaras question egp OpenFOAM Running, Solving & CFD 45 October 28, 2010 03:30


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