CFD中文网

@先先贝 课题组的求解器，rhoCentral耦合化学反应的，，但是我这个算例没开化学反应，相当于原生的rhoCentral

从去年开始导师开始对AI For Fluid产生了超级浓厚的兴趣，于是本人也开始做了一些工作，主要是纯数据驱动这一块（就是把CFD计算出来的物理场当成数据集，按照图像处理那一套东西，一顿卷积，然后预测出一个结果）。经过了半年多的研究，我发现这个东西非常不靠谱，存在特别多难以调和的问题，刷到国内外的一些帖子，发现大伙基本都在给纯数据驱动for fluid判死刑。 于是最近我看了一些文章，发现一些其他思路，比如PINN，这部分岳哥很慷慨的在论坛里开源许多代码。除此之外，好像还有DeepONet（深度学习算子）和FNO（傅立叶神经算子）这两个关注度比较高的思路，于是想来论坛问问诸位大哥，对这一块有什么看法。 另外，本人研究方向主要是气固两相流，不知道在Gas-Solid flow for AI这方面，国内外有没有一些比较不错的工作，希望诸位大佬不吝赐教。

这里面除了无滑移还有剪切应力，但是我那选择moving wall的时候就剩一个无滑移,

@heyi999 偶尔还在用之前的代码跑算例，有些细节可能有所以往，不过可以交流下

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

@旺财R @HIT_lly @youhaoyu @1064168551 已发

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

@先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了，你设置的自适应加密等级也基本到位了。这个时候再运行一次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镜像同样上传至百度云. 安装方法不再赘述.

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

@李东岳 是的。我又从头到尾排了一边，好像是求解马赫数的方式有点问题，现在还在验证，希望是这个原因

@李东岳 感谢李老师，我准备转向org版本进行研究，采用全新求解器重新实操一下

求助，目前在使用hy2Foam求解器计算类似乘波体外形，已经计算过钝头算例可以得到收敛解。但目前计算外形相对复杂，二维计算，计算域为1*0.55m，网格数量42000。计算过程中经常报错，求问有没有什么调试建议。下面是报错信息。

Mean and max Courant Numbers = 5.29417e-05 0.0491796 deltaT = 3.99074e-11 Time = 2.017990284e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000161695, Final residual = 3.0494e-13, No Iterations 1 DILUPBiCG: Solving for O2, Initial residual = 0.000148597, Final residual = 1.8029e-13, No Iterations 1 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 1.07385e-06, Final residual = 1.49309e-06, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 8.83461e-07, Final residual = 1.19425e-06, No Iterations 1000 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 2.94288e-05, Final residual = 3.15326e+222, No Iterations 400 Attempt to use rho2ReactionThermo out of temperature range 11152 times during this iteration. -1.64159e+236 < Tlow: 200 Phase no 1.0 ExecutionTime = 3564.14 s ClockTime = 3564 s Iteration no 4641 (0.79 s) Mean and max Courant Numbers = 5.41311e-05 0.0489379 deltaT = 4.07806e-11 Time = 2.018031064e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000162296, Final residual = 1.67236e-16, No Iterations 2 DILUPBiCG: Solving for O2, Initial residual = 0.000148383, Final residual = 6.70388e-18, No Iterations 2 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:? #9 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #10 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #11 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"

嗯嗯，还是谢谢老师了，我再比较比较文献和我程序实现里面的不同看看

这两个不能叠加使用。