CFD中文网

@xpqiu
1161114073@qq.com
感谢老师！！

求助：如果我需要将热膨胀系数编译进类rhoThermo中，我需要怎么操作。仅仅在rhoThermo文件中添加该项会出错，可能是这是一个基类，在运行时会影响其他子类，请问这个该怎么处理呢？

@hery 图片未上传上，如下：

用buoyantPimpleFoam求解器在kEqn亚格子LES计算进口速度0.21m/s的圆管流动，管径在3mm。div（phi，u）用的是Gauss LUST grad（U）格式，有壁面热流。入口处用的是mappedPatch，计算的泰勒尺度在0.0004左右，时间步设置为4e-5，开始计算收敛较好，计算一阵子后总会发散，而且查看云图会发现整根管的流体除入口处全部回流，请问这是什么原因啊，是离散格式的原因，还是偏离实际情况吗，我用稳态buoyantsimpleFoam计算可以很好收敛？
细节图片如下：

好的，谢谢老师

各位老师好，请问做大涡模拟时进口条件像设置成循环赋值的形式，比如将流场0.3m处的值从新copy到入口处这样的该怎么设置啊，用cyclic边界条件吗？
形式类似于下图这样：

请问在fluent2020R2中怎么设置v2f湍流模型，我在viscous中没有找到？

@东岳 用的komega湍流模型，有时在outlet处设置了intOut边界，value会设置成0，也尝试过zerogradient，没有value
，0文件夹的omega设置如下：

internalField   uniform 53.6;

boundaryField
{
    INLET
    {
        type            turbulentMixingLengthFrequencyInlet;
	mixingLength	8.52e-4;
        value           uniform 53.6;
    }
    OUTLET
    {
        type            zeroGradient;
    }
    BOTWALL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    UPWALL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    HEATWALL
    {
        type            omegaWallFunction;
        value           $internalField;
    }
}

我对直径为2mm的圆管中液态二氧化碳换热进行计算，进口温度298.15K，粘度为0.00007左右，用的是fluent里面获得的komega模型的初始值，k=0.000628，omega=713，进口的质量流量是0.001022，速度大概是0.4098，其他的稳态条件设置的与fluent默认值差不多，但是一计算第一步就发散，发散情况如下：

而如果用下面公式：
计算获得omega初始值大概为50左右，这样的初始化又可以计算，没有在第一步就发散。
我想请教下大家，这是什么原因呢，是因为openfoam的omega定义的和其他cfd软件不同吗？而且我用微信中的cfd工具箱计算的比耗散率甚至达到了265626，这个数值又与前面两个数值有什么差别呢？