Nut和nuTilda【翻译】

翻译自：CFD-online

帖子地址：http://www.cfd-online.com/Forums/openfoam-solving/72176-nut-nutilda.html

phile：

亲爱的Foamer们，

我是一个OpenFOAM的新手，因此我的问题可能很愚蠢。nut和nuTilda有什么区别？nut是否就是湍流运动粘度？那么nuTilda又是什么呢？

谢谢

Phile

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

你好，phile

我对此也很困惑。我查看了OpenFOAM C++编程指南，你能在opencfd找到它。似乎nuTilda是大涡模拟的一个参数。然而，我对大涡模拟了解不多。如果可能，请你检验一下然后告诉我。

Jinbiao

niklas：

如果你查看src/turbulenceModels/incompressible/RAS/SpalartAllmaras目录可能会对你有帮助。滚动.C文件的底部，nut_.internalField() = fv1*nuTilda_.internalField();

JinBiao：

感谢Niklas。这对我非常有帮助。

很高兴找到答案。nuTilda应该是Spalart-Allmaras模型的变量。因此它属于平均雷诺模型中的Spalart-Allmaras模型。

也就是说，当我使用k-e模型时，我不需要它。

我正确吗？

Jinbiao

niklas：

正确

Philer：

nuTilda文件存在于tutorial的kEpsilon模型下让我很困惑。

现在我知道它是多余的。

谢谢你，朋友！

Phile

jiez：

nut（湍流粘性）在版本1.6以后是需要的。

DLC：

大家好！

谁能告诉我怎么初始化nut和nuTilda?我经常使用kOmegaSST模型并且我发现了初始化k和omega的公式，但是我没有找到如何初始化spalartAlllmaras模型的变量

tcarrigan：

如果我理解正确的话，如果在使用kOmegaSST模型时，你想要初始化和记录nut和nuTilda，简单的将它们放在0目录中并且设置他们为0...或者你基于k和omega计算它们。

当使用KOmegaSST模型时（只需要k和omega），nut和nuTilda是不需要的。通过将nut和nuTilda目录文件放置在0目录下，你就强制OpenFOAM每次迭代都记录这些值。当你完成模拟之后你不必自己去计算这些值，这样有利于你的后处理。

注：我认为nut和nuTilda仅在使用SpalartAllmaras湍流模型时需要。

DLC：

你好，Trevis

感谢你的快速回复...我意识到这个发布的帖子英语很糟糕（深夜...）...

我很清楚的知道当使用KOmegaSST模型时无需nut和nuTilda...

现在我尝试测试spalartAllmaras模型而且我还是不知道如何设置0目录下的nut和nuTilda...

（除了薄片周围我都是用的freestream边界，薄片的nut我使用的壁面函数而nutTilda使用的是fixedValue）。

Akuji：

大家好

我有一个问题...我使用spalartAllmaras模型，因此我需要nut和nuTilda。根据我的理解，nut是湍流粘度和nuTilda是湍流运动粘度。我是否正确？

我应该如何计算它们？

rob3rt 0ng：

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();

大家好

对于fv1取什么值好？有些人说0.1而其他人说0.5？

感谢并且祝好

Robert

rob3rt 0ng：

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

我已经得到答案了。对不起，打扰大家了。

感谢并且祝好

Robert

rafamusura：

rob3rt 0ng

能请你回答一下你原来问的问题吗？我还不知道怎么使用该值！

谢谢

rob3rt 0ng：

你好，Rafael

D. Wilcox的CFD的书籍当中找到湍流模型，fv1被定义为x^3/(x^3 + cv1^3)，其中x是nut/nu，cv1是封闭常数7.1，在我的算例当中，fv1增加到1.00

祝好

Robert

rafamusura：

好的，但是我很困惑x被定义为nutilda/nu。Nut仅被定义为nutilda 的"任意"值。从其他方面，我读到的nutilda=3*nut=nu*0.1是否正确？

谢谢

rob3rt 0ng：

x被nut/nu定义。Nut是可以通过你设置的free-stream湍流参数通过计算得到的湍流运动粘度。

祝好

Robert

Tobi：

大家好

仅有一个问题

SpalartAllmaras模型仅是一个分析描述为了得到湍流的长度尺度和为了计算mu_t的湍流的动能。我参考的Ferziger and Peric。

因此，这个模型在替代雷诺平均当中的双方程模型的优势在什么地方？在我的印象当中，应该是双方称能得到更好的结果比起使用单方程或者分析描述，不对吗？可能在过去，计算机资源非常昂贵，一方程模型适用那个时候而不是现在？

Tobi

rafamusura：

Robert Ong

我不确定你的那个方法是否能够运行，因为你定义了nut=fv1*nutilda, fv1=x^3/(x^3 + cv1^3)，正如你所述的，x是nut的一个函数，如果你定义nut是一个nut的函数，nut=nutida*f(nut)。

在wiki网页上http://www.cfd-online.com/Wiki/Spalart-Allmaras_model，x被定义为nutlinda/nu，因为哪一个才是正确的解释？

Tobi，我和你有相同的问题，我尝试使用k-e模拟，但是它总是崩溃，因此退后一步首先使用一个简单的模型。如果你知道问题的答案请让我知道。

Tobi：

Rafamusura你好

这个时候，我正在撰写我的硕士论文并且在理论上更深入。所有的纳维斯托克斯方程以微分方程呈现，因此它的逻辑性和意义在于使用微分方程来计算湍流粘度。因此在简单的理论中-无，一方程或者两方程-，你应该参考两方程，特别是在你使用复杂几何或更好的雷诺模型的时候，参考Ferziger and Peric.的应力张量模型。但是你必须确定你的模拟需要达到什么精度，并不是每个模型都能很好的针对所有的几何和问题。很可能一方程的解与两方程对某一个几何的解是一样的而对于其他几何则不是。对于雷诺平均湍流模型我不是专家然而对于粘性模型，我认为你应该使用双方程模型封闭纳维斯托克斯方程。

你的k-e问题是边界问题或者你初始化的问题。但是那是另外一个主题了。

如果我错了，请纠正我！

vishal_s：

大家好，你们是否得到了x的正确解释？

aerogt3：

Quote:

Originally Posted by vishal_s 

Hello Everyone, did u got the correct interpretation for X??

一样的问题，如何从freestream湍流特性计算它

vonboett：

X被nutilda/nu指定，nutilda在壁面和freestream当中都为0，但是为了稳定原因，nutilda < 0.5，nu可以被运用。在算例中，你仅需要将它作为一个起始值（在Hybrid URANS_LES DES, DDES使用Spalart-Allmaras）。Wikipedia建议使用nutilda = 5nu在freestream，用来造成湍流发展为剪切流。（参阅http://www.cfd-online.com/Wiki/Turbu...ary_conditions）

SCG：

在提到的出口的nut文件应该设置什么值？

tareqkh：

如果你知道了，请分享

祝好

alyucel：

Quote:

Originally Posted by Philer 

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

这个问题还没有答案

Nut>动态湍流粘度

nutilda>湍流运动粘度

对吗？

他们之间有什么不同？

msuaeronautics：

Quote:

Originally Posted by alyucel 

This question was not answered yet 
nut> dynamic turbulence viscosity
nutilda> kinamatic turbulence viscosity

aren't they? 

and what is the different ?

"nuTilda"场是一种用来求解所有的基于Spalart-Allmaras湍流模型一个输运方程的伪湍流粘度（SA-RANS或者SA-DES模型）。

"nut"场是真实的"湍流粘度"并且通过fv1函数与"nutTilda"成比例，其中fv1是一种nuTilda/nu的局部函数。

人们总是应该合理初始化"nuTilda"场。

aerospain：

亲爱的Albrecht

这里有一个关于OpenFOAM3.1从nuTilda初始化nut的简单问题（SA模型）。我是否应该手动的从nuTilda计算0/目录下包含的nut文件内的值？（使用fv1和Cv1），或者OF在模拟开始的时候自动生成nut文件？

常常，为了保险起见，我会手动初始化0/目录下的nut文件，但是我怀疑是否会与OF当中的某些算法相冲突当不同版本改变的时候。

谢谢！