CFD中文网

超算中的参数有节点数、CPU数、任务数，OpenFOAM自身有个划分网格的数量。怎么设置可以最大化计算效率呢？

@学流体的小明
1、谢谢大佬！就是这个原因，mesh.Cf()改成mesh.C()就没问题了！
2、当时注释这个是因为我发现有问题后用Info输出坐标和mesh的size想看看问题出在哪，但是只输出了processer0的网格信息，所以有了这个疑问。后来我查到了是Info只会输出master processor的信息，我改成Pout就可以输出所有processor的信息。

@学流体的小明
红线包围的是我限制的区域，黑线表示的是网格，这样看应该有满足限制条件的网格，但计算得到的结果是bForce全是0，没有任何赋值的情况。

@学流体的小明

volVectorField bForce
(
    IOobject
    (
        "bForce",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedVector("bForce", dimensionSet(1, -2, -2, 0, 0, 0, 0), Zero) // dimless
);

for (label cellI = 0; cellI < mesh.C().size(); cellI++) // only process 0?
{
    const scalar x0 = 0;
    const scalar y0 = 0.0245;
    const scalar x = mesh.Cf()[cellI][0] - x0;
    const scalar y = mesh.Cf()[cellI][1] - y0;

    if (x > 0 && x < 0.012 && y > 0 && y < 0.0005) {
        bForce[cellI][0] = 50;
    }

@学流体的小明 谢谢大佬，我昨天写好了。但是目前的问题是体积力场添加的位置不对，我现在改成了最简单的模型来检查错误，希望在一个长方形区域中添加固定数值的体积力场，还是不能实现。

体积力场是随时间变化的量，是(x, y)的函数

@李东岳 不对的，实验结果约为1.2。

@李东岳

@李东岳 圆柱尾涡不是周期性的。maxco从1.91——0.19。

@李东岳 老师我遇到个问题，使用y+=0.2的网格和laminar模型进行计算，目前算的1e-4、5e-5、2e-5、1e-5的流场都没有收敛。这是为什么呢？是网格不合适吗？

@李东岳 好的老师，我算完发上来。

@李东岳 老师，层流的话雷诺数定在多少呢？网格还需要换吗？

@李东岳

    front     
    {
	type            empty; 
    }
   
    back 
    {
	type            empty;  	
    }

    wall
    {
	type            nutkWallFunction;
        value           uniform 0;
    }
	
    inlet
    {
	type            calculated;
        value           uniform 0;
    }
	
    outlet
    {
	type            calculated;
        value           uniform 0; 
    }

@李东岳 李老师，这是我换v2306的结果。

背景：
二维圆柱绕流

网格拓扑结构：

时间无关性的研究：

对于同一套网格，阻力系数随时间步长的变化如下图所示。结果表明，当最大库朗数的值不过分小时，时均阻力系数随着时间步长的减小有收敛的趋势。当最大库朗数的值的量级到1e-2时，时均阻力系数发散。

• 疑问：是否可认为在openfoam的计算中，并不是时间步长越小，计算结果越精确。
  

网格无关性的研究：

• 主要针对y+
  对于相同时间步长，阻力系数随y+的变化如下图所示。结果表明，当y+的值不过分小时，时均阻力系数几乎收敛到同一个值。当y+远小于1时，时均阻力系数发散。
  

结论：

根据目前的结果来看，y+过小或者dt过小都使得阻力系数发散。

• 疑问：是openfoam针对小的y+和dt有什么额外的设置吗？应该如何理解目前的结果呢？谢谢各位老师同学。

好的老师，我算完给您回复。

@李东岳 老师抱歉，最近一直没登录论坛，我现在参照一个参考文献更改了网格，还是存在和之前一样的问题，即减小时间步长后的Cd没有收敛到一个固定值。

@李东岳 老师，请问现在关于时间步长无关性的研究有结论了吗？

@李东岳 cylinder_1e-4_last.7z 老师这是1e-4包含第5s结果的文件夹。

@李东岳 好的谢谢老师！