Question

我正在尝试关注this OpenFOAM tutorial使用Gmsh生成轴对称网格。文件提供here。但是，当我尝试使用pimpleFoam解算器解决问题时，我收到错误：

- ＆GT;泡沫致命错误：   未知的TurbulenceModel类型RASModel

有效的TurbulenceModel类型：

3   （   LES   RAS   层   ）

来自函数static Foam :: autoPtr＆gt; Foam :: TurbulenceModel :: New（const alphaField＆amp;，const rhoField＆amp;，const volVecto rField＆amp;，const surfaceScalarField＆amp;，const surfaceScalarField＆amp;，const transportMo del＆amp;，const Foam :: word＆amp;）[with Alpha = Foam :: geometricOneField; Rho = Foam :: geom etricOneField; BasicTurbulenceModel = Foam :: incompressibleTurbulenceModel; Trans portModel = Foam :: transportModel; Foam :: TurbulenceModel :: alphaField = Foam :: geometricOneField; Foam :: Turbulenc eModel :: rhoField = Foam :: geome tricOneField; Foam :: volVectorField = Foam :: GeometricField; Foam :: surfaceScalarField = Foam :: GeometricFi eld; Foam :: TurbulenceModel :: transportModel = Foam :: transportMo del]   在文件/opt/CFDSupportFOAM4.0/beta/OpenFOAM-dev/src/TurbulenceModels/turbul enceModels / lnInclude / TurbulenceModel.C第113行。

正如this page中所解释的那样turbulenceProperties中case/constant的语法已经发生了变化。所以我编辑了turbulenceProperties文件：

simulationType RASModel;

到

simulationType RAS;

RAS
{
RASModel kEpsilon;

turbulence      on;

printCoeffs     on;
}

然后我得到了一个不同的错误：

泡沫致命IO错误：   试图超越EOF

file：... / Axisymmetric2D / case / system / fvSchemes.divSchemes.default at 29。

来自功能虚拟泡沫:: Istream＆amp;泡沫:: ITstream ::阅读（泡沫::令牌放大器;）   在第82行的文件db / IOstreams / Tstreams / ITstream.C中。   FOAM退出

本教程似乎适用于较旧版本的OpenFOAM。如果你能帮助我知道问题是什么以及如何解决它，我将不胜感激。

我的目标是学习如何使用Gmsh制作轴对称网格。所以开箱即用的解决方案或新版OpenFOAM als的教程都可以。

P.S。我已经报告了问题here in the Github repo

Answer 1

通过查看official OpenFOAM github repo（version 2.3.x和version 5.x）中提供的不同版本的axisymmetricJet模板，我能够解决问题。要做出的改变：

在case/constant/RASProperties中添加以下内容：

kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3; }

case/constant/turbulenceProperties

将行simulationType RASModel;替换为：

simulationType RAS;

RAS {

RASModel            kEpsilon;

turbulence          on;
printCoeffs         on;

}

case/system/fvSchemes

将行div((nuEff*dev(T(grad(U))))) Gauss linear;更改为div((nuEff*dev2(T(grad(U))))) Gauss linear;

解算器按预期收敛。我仍然不知道这些变化意味着什么以及它们如何运作。我一想好就会加上它们。我已经分叉the github repo here，包括所需的修改。

P.S。亲爱的主持人，出于某种原因我的代码格式化不起作用！

OpenFOAM：FOAM致命错误：未知的TurbulenceModel类型RASModel

1 个答案: