CFD中文网

@izumi

//const tmp<volScalarField> tmuEff = turbModel.muEff(); // 错误，智能指针不能以引用方式返回 const tmp<volScalarField> tmuEff = turbModel.muEff(); // 正确 // 获得被封装对象的引用 const volScalarField& muEff = tmuEff(); // OpenFOAM 4 以前的用法 //const volScalarField& muEff = tmuEff.ref(); // OpenFOAM 4 以后的用法

@东岳
另一个alphatJayatillekeWallFunction和Fluent中的壁面函数一样，根据$y^+$使用不同的公式计算，打算试一下这个

应该研究下为什么transonic需要提升对角占优

@李东岳 回答一下自己的问题。transonic下的压力方程存在对流项，对流项不利于对角占优矩阵的产生，因此transonic的压力方程要增加方程松弛。

http://dyfluid.com/rhoSimpleFoam.html

@东岳
谢谢！用Surface Vectors filter就直接画出来了

抱歉，找到原因了，是我忘记wclean了

虚拟机肯定是没问题，跟Ubuntu真系统感觉没什么区别。
但是BSL下面的OpenFOAM能感觉出来有点慢，之前对比过，主要是来自于磁盘读写。如果你有兴趣的话，可以试试对比下，谢谢。

我看的一篇文献，在CFX中使用SSG模型，网格y+~1，但是我在OpenFOAM中网格y+较大（>11）时，能够得到收敛的结果，但是加密网格到y+~1时，就一直发散。
请问：SSG模型要求的y+是多大？此时nut、k、epsilon和R的边界条件应该怎么给？谢谢！

我大略看了一下：

#include "createFluidMeshes.H" #include "createSolidMeshes.H"

划分两个网格的时候就存在了边界场，在求解fluid之后，边界场信息会更新，你不需要显性的传递边界场。

fvMesh& mesh = solidRegions[i]; solidThermo& thermo = thermos[i]; tmp<volScalarField> trho = thermo.rho(); const volScalarField& rho = trho(); tmp<volScalarField> tcp = thermo.Cp(); const volScalarField& cp = tcp();

这几行代码自动包含边界场，其中给固体的是网格共享的边界场信息。

应该是这个思路， 你可以确认下。

看论文发现有作者说fluent稳态计算残差降到1e-7则认为收敛停止计算，但是我用of的simpleFoam的时候发现残差根本降不到那么低，连降到１e-６都已经是极限了。我都不敢在论文里写残差判断标准了，一般你们怎么说收敛标准呢

@exthan 您好，我现在还没找到正确的方法，还在尝试中。

选用 kappaMethod fluidThermo;时，kappa是根据alphaEff*thermo.Cp来计算的，在我的例子里面kappa不是恒定的，那么用fluidthermo是否会存在问题？我算出的热流密度跟原来设定的热流密度有点差别，所以有些疑惑。

@vbcwl constant/thermoPhy 文件里面去设置。

thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; }

@李东岳 哈哈哈，物理隔绝了:mihu: