CFD中文网

你这个方法是对的，是不是拼写写错了，inelt.stl？

可以的

看看这个示例https://github.com/ancolli/kOmegaSSTCC/blob/master/kOmegaSSTCC/makeTurModel.C

速度矢量方向对吗 看起来没进去气体

可能是内存不够，用top看一下

TRPjet是啥，估计应该是fluid:592
solid2
solid1 之一

这就得找论文看看了

把出口映射到入口

“拷贝”指的是映射过去，openfoam有个mapFields用法

如题，两种粒子相互独立，不考虑两种粒子间的相互作用，是否要在createFields.H中创建两个kinematicCloud，如kinematicCloud1，kinematicCloud2，分别对应case/constant下的kinematicCloud1，kinematicCloud2，再分别指定物性

Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
basicKinematicTypeCloud kinematicCloud
(
    kinematicCloudName,
    rhoc,
    Uc,
    muc,
    g
);

@SHUKK
有点久远了额，湍流模型参数betaStar要取0.028，测试了下高度方向网格50个，simpleGrading高度方向取5，有如下结果：

感觉网格还能调试一下，低处会保持得很好

0/U boundary type movingWallVelocity

//k
rho_*epsilon_/k_
//epsilon
rho_*Cmu_*sqr(k_)/epsilon_

为k_和epsilon_接近0，导致除法崩溃，所以使用bound（）

可以试试增大epsilon 减小k；连续性残差也偏大，可以试试fvSchemes散度项二阶换一阶。

sample.lines 去掉符号“.”呢

@疏影横斜水清浅 在 大气边界层迎风坡台阶地形的风场模拟 中说：

@cccrrryyy 不是发散，是结果不对，非常明显。我已经找到错误原因了，我的压力边界设置错误了，设置成零梯度了，改成零压力计算结果就可以了。再次谢谢您，也谢谢各位大佬吧友的解答。

你好，请问你说的是上边界的压力零梯度改为固定值0就正常了吗

手动直接改那个时间步的边界条件就可以，只不过需要下拉到最下面一点

decomposePar -latestTime 可以不

@李东岳 在 请教一下标准k-epsilon湍流模型中的epsilon_()与k_()与epsilon_和k_有什么区别 中说：

目的：减少并行通信，从而降低并行计算的计算成本。（只有boundaryField才需要通信）

我不觉得是这个原因。我觉得应为这个是源项。源项是源项。源项不应该对boundary产生印象。

感谢李老师回复， 我再查查看看

https://www.cfd-online.com/Forums/openfoam-programming-development/224302-difference-between-rho-rho_-rho.html
@李东岳 李老师，这里。修正下：大部分通信发生在boundaryField
另外，两个湍流模型里的divU和G的定义也不一样：

//standard k-epsilon
    volScalarField::Internal divU
    (
        fvc::div(fvc::absolute(this->phi(), U))()
    );

    tmp<volTensorField> tgradU = fvc::grad(U);
    volScalarField::Internal G
    (
        this->GName(),
        nut()*(dev(twoSymm(tgradU().v())) && tgradU().v())
    );

//LaunderSharmaKE
    volScalarField divU(fvc::div(fvc::absolute(this->phi(), U)));
    tmp<volTensorField> tgradU = fvc::grad(U);
    
    volScalarField G(this->GName(), nut*(tgradU() && dev(twoSymm(tgradU()))));

可以看到，标准模型里都是InternalField

epsilon，k在声明时：

volScalarField epsilon_；
volScalarField k_;

这俩变量类型都是volScalarField，volScalarField就包含internalField 和boundaryField ，epsilon_()和k_()以及rho()和alpha()中的 operator()返回internalField ，类型相当于volScalarField::Internal

目的：减少并行通信，从而降低并行计算的计算成本。（只有boundaryField才需要通信）

注意：epsilon_()和k_()以及rho()和alpha()使用时仅在方程右边构造源项