请问你运行成功了吗？我按照源代码给的示例这样设置，程序读取的时候报错

@小菜鸟 在 OpenFOAM python PINN tensorflow gpu (小白，莫笑） 中说：

@李东岳 李老师，deepxde库属于基于pytorch tensorflow做了封包，形成的高级科学计算库，把计算物理方程的部分模块化掉了，在计算性能方面可能有做过优化，但是降低了自定义的灵活性。还有jax，这个会更加完善，我习惯用pytorch，用jax不太习惯，jax的风格有类似于tensorflow吧(具体不太了解)。国内的有百度系的paddle science和华为系的MindSpore 的MindSpore Flow。

老兄能否留个联系方式交流一下，我也在整这个，正在入门。

@xpqiu

十分感谢老师的指点！

我检查了所有异常算例，发现问题的确都与0文件中的初始设置有关。

在这个算例里，命令subsetMesh -overwrite frameCellSet -patch frame会将frameCellSet网格区域与frame链接，但是由于没有显式划分frame，因此，subsetMesh会默认生成一个empty类型的边界。

六楼的分析存在大问题，没有参考的必要了。

再次感谢两位老师的指点:146: ~

@KrisIrving 留个档吧，这样操作是可以的。只不过模型尺寸在微米，我用magics查看会出错，用Microsoft商店的3D viewer可以正常查看，也可以利用Blender进行查看和修改stl文件

后来这个怎么样了

在一维CFD密度基可压缩流动中，使用亚音速（总温总压）进口边界条件，进口静压使用线性外推，其余进口物理量靠总温总压和总温静压的关系给出，为什么算出来的进口处附近有震荡，有什么好的解决方法吗？

@李东岳 确实是我驱动段给的压力和温度不够，我错把激波后的压力和温度当成驱动段的压力和温度了，现在重新设置驱动段的压力和温度后，确实能够达到高马赫数了

大佬们，求问有没有什么方法可以简单的修改边界条件的值呢。感觉计算不同条件的case，一个个设置P、T、U边界条件的值很麻烦。能不能在边界条件类型已知情况下，直接输入一个txt文档的边界条件值，就能完成所有边界条件值得修改呢。用openFOAM7

@wangfei9088 非常感谢

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

@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就是冲突的东西。所以如果做这方面研究，自圆其说是最重要的。

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

@先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了，你设置的自适应加密等级也基本到位了。这个时候再运行一次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 ，顶大佬一波！！

感谢大佬回复@coolhhh

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

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

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

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

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

@五好青年 大佬，后来者大胆问下这个有解决的好办法或者好思路了么