参考：https://zhuanlan.zhihu.com/p/427686653
另外注意网格的尺寸，of里面是米的单位，一般需要使用transformPoints从mm转到m

@李东岳 静态网格也是有这个问题的。按理说不应该啊...我再跑跑其他类似的原生算例

p：

#include "../varsSettings" //Include global variables dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value $internalField; } frontAndBack { type empty;//zeroGradient; } topAndBottom { type zeroGradient; } frame { type zeroGradient;//fixedFluxPressure; } }

U：

#include "../varsSettings" //Include global variables dimensions [0 1 -1 0 0 0 0]; internalField uniform $flowVelocity; boundaryField { inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; inletValue uniform (0 0 0); value $internalField; } frontAndBack { type empty;//slip; //value uniform (0 0 0); } topAndBottom { type slip; value uniform (0 0 0); } frame { type noSlip;//movingWallVelocity; //value uniform (0 0 0); } }

算例文件 OpenfOAM2012 网格数4w
functionalTest_NoPZRefine2_2D_staticMesh.zip

现在收费了么？我还不知道。好久没关注了。我手头也有个免费版。好像起码得3 5年之前的了。我当时还好奇，为啥免费下的大家还留邮箱要。我一直以为跟我那个书一样，我那个书网上免费下，还有人付费买。原来是收费了。那要是收费了确实可以联系原作者购买，支持一下。

不知道卖多少钱。我手头有一本I do Like CFD。也是类似个人写的。我买的原版，220元人民币。电子版也有。

@李东岳 李老师好，这里我请教一个问题->链接

为了简化算例，我把网格三维转化为二维，简化了结构形式，以便排错。

【算例云图】
6101495a-adbe-455d-ab41-45e2b1b7b415-image.png
bb59f1c1-9dde-4169-83d9-cf00c0daf433-image.png

【简化版算例下载】
functionalTest_NoPZRefine2_2D.zip
生成网格createMesh，
跑算例runCaseOnly，
全流程runCase

感谢~

@李东岳 嗯嗯，谢谢老师！

@hy1112006 您好，如果方便的话可以有您的邮箱或者联系方式交流一下使用pisoFOAMopenfast吗？

想尝试与一篇文献中实验结果作对比，但根据文中所提供的条件算出的入口速度很小（0.18m/s），以此作为入口的边界条件；同时入口边界的alpha.water也固定为1，计算一段时间后如图所示，发现入口处几乎没有水相的存在了。这是为什么？

ae0a7d49-651d-4636-a27b-c8a0ff4b9334-image.png

入口处的边界条件如下：

U { type fixedValue; value uniform (0.18 0 0); } p_rgh { type zeroGradient; } alpha.water { type fixedValue; value uniform 1; } alpha.air { type zeroGradient; }

我觉得这个是一个自圆其说的东西。如果看sci，肯定有人对比LES和DNS的结果，证明他们的LES或者DNS的结果更好。

但我认为，本身欧拉拉格朗日里面存在了很多模化，在粒子受力的角度，存在更多的模化，尤其是湍流分散力那面，基本都是参数。本身DPMFoam和DNS就是冲突的东西。所以如果做这方面研究，自圆其说是最重要的。

https://cfd-china.com/topic/8018

能做的话可以上传到这里我给你看看

您好，我最近研究这方面的问题，请问您参考的这篇文章标题可以发一下吗，我想学习一下

@余正东 这个我关注了，未来可期，目前只是初步移植了一些

@先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了，你设置的自适应加密等级也基本到位了。这个时候再运行一次setFields，给加密后的较密的网格赋值，再去分块。
分块的话也得想想加密的部分会不会在你流场里面移动，你的不同核之间的负载要基本平衡，可能还需要你手动去划分区域。

OpenFOAM中使用multiDimAMR编译和运行正常，使用reconstructParMesh也正常，但是之后再使用reconstructParMesh会报错。
第一个时刻还是正常的，第二个时刻发生报错，还有一个算例是第5个时刻开始报错，报错内容一样。使用paraview中的decomposePar也能查看结果，为什么会出现这个问题呢~
Reconstructing fields for mesh region0
Time = 0.01
Reconstructing FV fields
Reconstructing volScalarFields
alpha.water
alpha.water_0
p
p_rgh
rAU
Reconstructing volVectorFields
U
Reconstructing surfaceScalarFields
alphaPhi0.water
phi
Reconstructing surfaceVectorFields
Uf
Reconstructing point fields
No point fields
No lagrangian fields
No FA fields

Time = 0.02

--> FOAM Warning :
From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 210
Number of patches has changed. This may have unexpected consequences. Proceed with care.
--> FOAM Warning :
From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 210
Number of patches has changed. This may have unexpected consequences. Proceed with care.

**--> FOAM FATAL ERROR: (openfoam-2012)
Processor 2 has a different polyMesh at time 0.02 compared to any previous processors.
Please check time 0.02 directories on all processors for consistent mesh files.

From Foam::polyMesh::readUpdateState Foam::processorMeshes::readUpdate() in file processorMeshes.C at line 192.**

FOAM exiting

刚好在用kinematicParcel ，顶大佬一波！！

@学流体的小明 是的 想通过设置入口湍流速度边界来发展出槽道内湍流。刚接触这一块，觉得是不是可以 用贴子这个初始扰动的方法，设置成入口边界，结果试了一下单纯加上去不太行。就想问问有没有类似的设置入口边界的方法。

大佬们，这个算例用rhoCentralFoam能跑通。因为后面想加入颗粒计算，现在在同样的边界条件下我用rhoPimpleFoam计算时，却一直卡在这里，入口的压力也不对，总压给定的是3e6。这是什么问题呢？
image.png

image.png

image.png

感谢大佬回复@coolhhh

这两个图分别对应 方法1 和 方法2 的计算结果吗，还是不同方向的结果？

不同方向的结果。$\xi$是$x$方向的间距，$\zeta$是$z$方向的间距。上面我画的都是相关系数，也就是相关函数进行了归一化，所以在间距$\xi$或$\zeta$为0时，相关系数都为1。

看起来对于槽道流这种展向和顺流向是均匀+时间平稳的流动是可以用时间平均之后的。问AI，AI也是这么说的。

我的数据在时间上和空间上也比较多，$\phi(x,z,y=y_0,t)$每个计算网格的值都输出了，时间上也以很小的时间间隔输出了一段时间。无论是时间上的相关系数还是空间上的，计算起来都没什么难度。

东岳流体虚拟机4月更新, 配套ISO镜像同样上传至百度云. 安装方法不再赘述.