CFD中文网

太有帮助了! 感谢

感觉就是不同的研究者对这个的规定不一样，论文里面推导过程正确就行。
当然如果按照张量的知识，向量并积计算的话，我觉得应该是
$$
\nabla {\bf{U}} = \frac{{\partial {U_j}}}{{\partial {x_i}}}
$$

😅我自己写的程序必须是simple的分解，scotch都不能用，更不用说metis了

不好讲故事是比较困境了，哎..

学流体的小明 在 interFoam计算气泡槽道流时的压力问题 中说：

应该是 并行 + 自适应网格 的问题？

更正一下，是 并行+自适应网格+压力参考点 三个要素的共同作用，去掉任何一个都可以算出来好的结果。

一个新的发现是壁面上会有奇异点，如下图，着色是p_rgh，这些奇异点都处于并行分区的各个界面上。更奇怪的是，上壁面并没有这样的问题。Note:参考点放置在流场中间。
01292e97-19c8-45be-889b-189c0ff6fe54-image.png

现在暂时放弃自适应加密了，直接全场加密算，压力就是正确的。
😂

@学流体的小明 不是，我的意思是说，对于处理器1的边界网格单元，他会有一个边界面，这个边界面，会跟处理器2的边界面对应，那么处理器2的边界面必然对应一个边界网格单元，去这么找。

这个问题我去年处理过，我就是这么搞的。不过我那个不是多processor，我那个是多计算域，也涉及到互相传递信息。下面这个图：

微信图片_20230721095200.jpg

感谢分享！！:146: :146: :146:

OpenFOAM-v2206/tutorials/incompressible/pimpleFoam/LES/planeChannel/setups.orig/common/system/controlDict，有一个计算Cf场的函数，可以参考下：

Cf { type coded; libs (utilityFunctionObjects); name Cf; timeStart $/timeStart; writeControl writeTime; codeExecute #{ auto* CfPtr = mesh().getObjectPtr<volScalarField>("Cf"); if (!CfPtr) { Info<< "Create skin-friction coefficient field" << nl; CfPtr = new volScalarField ( IOobject ( "Cf", mesh().time().timeName(), mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh(), dimless ); regIOobject::store(CfPtr); } auto& Cf = *CfPtr; Info<< "Computing skin-friction coefficient field\n" << endl; const auto& tau = mesh().lookupObject<volVectorField>("wallShearStress"); const dimensionedScalar Ubulk(dimVelocity, 17.55); Cf = mag(tau.component(0))/(0.5*sqr(Ubulk)); #}; }

更新一下进度
暂时不用可压缩求解器计算了，因为发现当前的网格算不到声波那种程度。不过我倒是实现了向温度方程内添加源项这件事，就是不太确定公式以及编程对不对。

系数B确实会随着雷诺数不同，有略微波动。《Kim J, Moin P, Moser R. Turbulence statistics in fully developed channel flow at low Reynolds number[J]. Journal of fluid mechanics, 1987,177:133-166.》文章中也提到了雷诺数比较低，B=5.5

908490eb-b40f-4968-ad79-3bc0ae94a791-image.png
a3903e77-0e01-4ba3-ab3b-dee1013905db-image.png

《Moin P, Kim J. Numerical investigation of turbulent channel flow[J]. Journal of fluid mechanics, 1982,118:341-377.》中的结果 B=5.0
c596e78d-8102-43aa-902f-d38d7fecd1cd-image.png

@李东岳 并行文件夹里面的2e-08文件夹是有pointProcAddressing文件的，但是后面算出来的2e-02，4e-02就没有这个文件了。所以reconstructPar说找不到这个文件，报错了。

更新一下
发现了问题的源头，是网格的问题。
我使用ICEM画网格，或者OpenFOAM的blockMesh工具画网格，都是先画背景网格，然后对近壁面进行多次加密。加密是八叉树类型的加密，所以在网格尺寸变化的地方，大网格的尺寸直接是小网格的二倍。这样的网格确实不好，可能一些拓扑结构上的不好导致了singleGraph程序进行插值时找不到插值的基点。
解决办法：
举例说距离壁面0~1的网格尺寸是1,1~2的网格尺寸是2，那么就使用singleGraph在0~1这个范围内提取一条线，在1~2再提取一条线，这时的type可以是uniform，也可以是midPoint或者midPointAndFace。一句话，就是分段提取信息。
我是重新画了网格， 老师说之前画的这网格根本不能用:mihu: