收藏一下，感觉有用

@李东岳 好的李老师，感谢感谢

@五好青年 老哥，我也会有这个问题，单核正常，并行发散，这是不是与分配核数有关呀？

@李东岳 是可以不选择的，但怕设置错了，这玩意用的人确实不多啊，老师还有啥推荐的网格软件吗

VIV.zip 这是我配置的openFoam动网格算VIV的算例，从目前的结果来看，程序能正常运行，但程序运行结束以后，圆柱不动，想请各位大佬帮我看看，哪儿出问题了，谢谢大家。

@疏影横斜水清浅 重力以及体积力

@sgggggg 硕士论文我也想拜读一下，可以加个联系方式吗2132306535

@李东岳 请问这可能是是网格太小导致吗？

各位大佬好，我是CFDEM新手，最近尝试参照twoSpheresGlowinskiMPI构建自己的算例，采用dynamicRefine，在计算过程中断构建的算例基本与参照算例一致，但试跑算例出现中断，不知道是否是哪里没有注意到，求大佬指导
报错描述
当不使用自适应网格时，算例正常计算；当采用自适应网格功能interDyMFoam，网格未加密前（未达到设置加密条件），正常计算;网格开始加密计算中断。
报错代码如下：

Time = 0.0001 Selected 0 cells for refinement out of 1613716. Selected 0 split points out of a possible 0. Courant Number mean: 0.000199582 max: 0.0804916 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' Setting up run at Fri Jul 26 12:59:28 2024 Memory usage per processor = 6.75634 Mbytes Step Atoms KinEng rke Volume dragtota dragtota dragtota 1 1 2.5261073e-08 0 0.0260508 0 0 0 CFD Coupling established at step 10 11 1 3.5918126e-08 0 0.0260508 0 0 0 Loop time of 0.000962037 on 16 procs for 10 steps with 1 atoms, finish time Fri Jul 26 12:59:28 2024 Pair time (%) = 4.70437e-06 (0.489001) Neigh time (%) = 0 (0) Comm time (%) = 4.65044e-06 (0.483395) Outpt time (%) = 6.52643e-05 (6.78397) Other time (%) = 0.000887418 (92.2436) Nlocal: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Total # of neighbors = 0 Ave neighs/atom = 0 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished Foam::cfdemCloudIB::reAllocArrays() nr particles = 1 evolve done. DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.06984e-07, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.99282e-07, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 6.95219e-10, No Iterations 2 DICPCG: Solving for p, Initial residual = 1, Final residual = 9.26695e-07, No Iterations 573 time step continuity errors : sum local = 3.69904e-10, global = -1.28591e-13, cumulative = -1.28591e-13 DICPCG: Solving for p, Initial residual = 0.0476412, Final residual = 9.26405e-07, No Iterations 513 time step continuity errors : sum local = 3.16922e-08, global = -3.49546e-12, cumulative = -3.62406e-12 DICPCG: Solving for p, Initial residual = 0.0055667, Final residual = 9.51869e-07, No Iterations 500 time step continuity errors : sum local = 3.3706e-08, global = -6.10801e-12, cumulative = -9.73206e-12 DICPCG: Solving for p, Initial residual = 0.00178549, Final residual = 9.94251e-07, No Iterations 473 time step continuity errors : sum local = 3.53268e-08, global = -5.71629e-12, cumulative = -1.54484e-11 No finite volume options present particleCloud.calcVelocityCorrection() DICPCG: Solving for phiIB, Initial residual = 1, Final residual = 9.63924e-07, No Iterations 501 ExecutionTime = 30.7 s ClockTime = 32 s Time = 0.0002 Selected 784 cells for refinement out of 1613716. Refined from 1613716 to 1619204 cells. Selected 0 split points out of a possible 784. Courant Number mean: 0.126921 max: 9.46381 - evolve() timeStepFraction() = 2 Starting up LIGGGHTS Executing command: 'run 10 ' Setting up run at Fri Jul 26 12:59:56 2024 Memory usage per processor = 6.75634 Mbytes Step Atoms KinEng rke Volume dragtota dragtota dragtota 11 1 3.5918126e-08 0 0.0260508 0 0 0.017347056 CFD Coupling established at step 20 21 1 3.6972707e-08 0 0.0260508 0 0 0.017347056 Loop time of 0.00149142 on 16 procs for 10 steps with 1 atoms, finish time Fri Jul 26 12:59:56 2024 Pair time (%) = 3.75399e-06 (0.251706) Neigh time (%) = 0 (0) Comm time (%) = 3.65471e-06 (0.24505) Outpt time (%) = 8.7527e-05 (5.86871) Other time (%) = 0.00139648 (93.6345) Nlocal: 0.0625 ave 1 max 0 min Histogram: 15 0 0 0 0 0 0 0 0 1 Nghost: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 16 0 0 0 0 0 0 0 0 0 Total # of neighbors = 0 Ave neighs/atom = 0 Neighbor list builds = 0 Dangerous builds = 0 LIGGGHTS finished nr particles = 1 -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 15 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. --------------------------------------------------------------------------

@李东岳 谢谢您

postProcess -func yPlus

--> FOAM FATAL ERROR:
Unable to find turbulence model in the database

From function virtual bool Foam::functionObjects::yPlus::execute() in file yPlus/yPlus.C at line 180.

非常感谢李老师！

请教各位大佬：
openfoam有没有自带的插值可以将volField转化成pointField，使用tecplot进行后处理时插值操作会将cell data转化成point data，操作比较麻烦于是自行编写openfoam程序，很多操作是基于volField，然后结果和tecplot有比较大的区别，不知道是不是插值的影响

我给你看了一下，openfoam2012这个没有植入颗粒动理学理论，算颗粒流不太准了。你可以使用openfoam8。

@李东岳 好的我简要描述下问题，就是我按照完整的形式推导了不可压缩流动的能量守恒公式。
4422b726-0b92-4704-82cd-c145d75dcc1c-image.png 但是案例结果并不满足。
我对这一项有疑问4326e809-66fe-472b-a04c-686c03c3cb2b-image.png
我在程序里输出是有值的，但是之前的回复中提到不可压缩流动中应当为零。
这一项是按照下面这个式子，考虑了速度梯度为零，去掉右边第一项积分来的
e470d379-bbaf-47b1-ba05-c1afb5d0c8d9-image.png
我发现这个式子中的三项，积分后的形式都是一致的，感到很困惑。

在openfoam 5.x下编译你试试 我那个虚拟机里面有
如果编译你都不知道怎么编译的话 你需要先学一下openfoam 然后再弄ibScourFoam

谢谢李老师，问题解决了，他们的程序调数据都要走一下虚拟硬盘，不是直连的，中间文件乱码了。

@李东岳 好的，谢谢东岳老师！

第一次在这个网站上求助各位大佬，报错图片不知道怎么看不了了，不好意思![IWE)}@6(NH98E9LRH@PDEK.png

@wangfei9088 不好意思哈，我这边没有描述清楚。是这样的，我用of2012当中reactingTwoPhaseEulerFoam里面的 fluidisedbed，用reactingTwoPhaseEulerFoam求解器计算双流体，这个文件里面第22行到后面都没有报错，以及我用这个算例，增加一个颗粒相的文件，用multiPhaseEulerFoam求解器，计算多流体，也没有报错，但是我用reactmultiPhaseEulerfoam求解器计算就报错了，说是不识别RASModel kineticTheory.不知道什么原因，如何修改。大神有遇到这种情况嘛？屏幕截图 2024-07-23 113729.jpg