CFD中文网

@黒洛克 提取完之后就是foam格式的，算是ascii编码的数据，文本文件就能读。
我是用matlab处理的。
用tecplot的话，你试试vtk格式的？再不成，就自己拿ascii的写一个tecplot能处理的文件。

我不太懂网格
只是感觉还是得从网格着手吧

我自己的经验：同一块代码，可以看看不同进程的执行耗时，看看是哪个进程比较慢，也许就拖累了进度。
我也就使用reduce、gather、scatter这些进程间通信的函数。执行到这些函数的时候，快的进程会停下来等慢的进程就绪吧。

在controlDict中使用#includeFunc wallShearStress就可以得到wallShearStress的场数据，应该不是东岳老师说的那个原因。
wallShearStress只在边界的patch上面有数据，提取一条线的数据容易选到没有数据的内部网格上？所以没有文件输出吧。
用surfaces提取面上的数据好了

#includeEtc "caseDicts/postProcessing/visualization/surfaces.cfg"

fields       (wallShearStress);

surfaces
(
    bottomWall
    {
        type patch;
        patches ("bottomWall");
        interpolate     false;
        invariant       true;  // Unaffected by mesh motion
        surfaceFormat   foam; //ensight vtk foam
    }
)

迭代到1000次是没收敛，和网格、边界条件、时间步长都有关系。
10mm量级的网格，要画出10um的狭缝，网格感觉就很难画了。得保证狭缝那里有足够的分辨率吧，看看网格？

@风云5091 我们这些都是槽道流里面的实际坐标，除以特征长度，得到的无量纲化坐标y+。
你如果只做验证的话，OF输出的yPlus足够了，应该可以在paraview显示出来。
OF输出的yPlus确实只存在于那个边界patch上。

是的，praview这里会有问题，我也遇到过。应该是paraview没有正确识别noSlip边界条件的原因。

h大小会影响到槽道的无量纲参数，比如摩擦雷诺数$Re_\tau=u_\tau h / \nu$，你要算湍流的话，网格要跟着$Re_\tau$去画。当然不能简单地缩放。

我不理解这句话

你说的这个是绝对速度，看相对速度的话，和转轴没多大关系

旋转中心设置为-2 m，计算域还是0到0.04m，计算域距离旋转中心就很远了啊，旋转产生的速度$v=\omega \times r$也应该很大的吧。

@Afr1yne 我还读在读博😂
我是拿pimpleFoam算的，网格都是精心设计过，和前人文献建议的一致，即使这样，OpenFOAM算出来的结果都不是很好。还是使用了更加精细的网格，得到的才和文献的DNS结果一致。

DNS计算要考虑到网格尺寸的，你改变计算域大小，网格尺寸也改变了，$ \Delta x^+$这些都改变了。
可能是因为大计算域时候网格尺寸不满足DNS的条件？就是计算域尺度方面的一些问题。

@bcp OKK，可以把怎样解决的发出来吗？大家有相同问题的话也可以参考你的方案。

取一个或多个平行于壁面的面（Kim的图也是不同高度的，y^+=5.39和149.23），使用surfaces工具输出这个面上的速度，然后每个面计算速度的展向、流向空间相关函数就行了。
可以直接计算二维的相关。
也可以固定z坐标，一条线一条线地计算流向的相关函数，再做个平均。

@bcp 我搞错了，之前没看对你计算rField的公式。我以为是你rField的计算公式不对。

你把报错代码发一下。
debug，分别输出rho、U、rField的量纲看一看

UEqn当中各项的量纲是 $[\mu][U]/L^2=M/(L^2T^2)$
rField的量纲应该是$T^{-1}$才对，但看起来你写出来的并不是这样
你可以在定义场的时候显式地指定量纲

volScalarField B                                            
(
    IOobject
    (
        "B",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar(dimless, 0.0)
);

平均流速就是用瞬时的流场速度U进行三维空间上的积分。可以看看代码
https://develop.openfoam.com/Development/openfoam/blob/OpenFOAM-v2012/src/fvOptions/sources/derived/meanVelocityForce/meanVelocityForce.C

和重力不挂勾的，重力产生的加速度是它自己的效果，记为A；这个体积源项本身会产生一个加速度，记为B。
A不变，B会一直变化以使得平均速度保持为你设定的那个值。

@东方白杨

东方白杨 在 展向平均 中说：
这一段意思就是对流场同时X、Z方向平均，得到沿Y方向的分布吗？

是的

东方白杨 在 展向平均 中说：
如果我只想得到流场沿X方向平均，得到一个YZ面，就是流向平均的一个截面，这该怎么操作呀？

这两个工具都是同时对两个方向进行平均的，如果想得到一个截面，我觉得你可以研究一下代码然后自己改改程序，或者再找找，可能别人写过相关的工具。

我自己一般在本地服务器做计算，超算上传下载比较麻烦。
内存1TB算1000万网格有点浪费吧，李老师卖的服务器，内存256G，我也能算2700万的网格，8000万网格的也不会爆。大网格情况下，CPU限制了计算速度，内存是确保程序能运行下去的。

还有1000万网格时需要开到64核以上计算吗？

可以往上开，计算速度应该会更快。可以先研究一下加速比再做决定。

感觉是topoSet选择网格，提取各种面的时候有些问题了，我没用过这方面的功能，不太了解。
你如果只是想把边界面提取出来的话，在blockMesh当中就可以设置吧。