OpenFOAM版本：OpenFOAM v9
求解器 ：multiphaseEulerFoam
模拟对象：圆柱状鼓泡床反应器，多相流，液相反应，无液相进出，气相不反应，底部进气，顶部出气。
算例来源：在tutorials的bubbleColumnEvaporatingReacting算例基础上修改而来，各相初始分布由setFields控制。

问题：即使关掉反应，从数值迭代上看，多相流计算第一个时间步的PIMPLE循环的相分数计算发散，见下日志。存下第一个PIMPLE迭代结果看，气相和液相的基础相莫名其妙地在原本该是0的地方，变成了1，见下图。
确认过setFields正确地初始化了相分布，上半部气相，下半部液相。由于第一次做多相流的算例，无法定位发散的原因是边界设置？constant的多相参数？相分数求解？还是别的问题？求各位大神指点一二，不胜感激。

尝试过的操作：

加松弛 改矩阵迭代求解器，GAMG换PCG，smoothSolver换PBiCGStab 调整气相液相和压力的出口边界类型为zeroGradient 关掉反应 加压力限制limitPressure 将气相从multiComponentMixture改成pureMixture

/*------------------------------------------------------------------------------
算例：
运行meshGen脚本生成网格，allrun脚本运行。
f21a20e8-830b-41ca-8175-7667fab9ab9b-p2.tar.gz

/*------------------------------------------------------------------------------
计算日志，第一个时间步

PIMPLE: Iteration 1 MULES: Solving for alpha.gas MULES: Solving for alpha.liquid gas fraction, min, max = 0.4999984327 0.009998944 0.990001056 liquid fraction, min, max = 0.5000016007 0.01 0.990001056 Phase-sum volume fraction, min, max = 1.000000033 1 1.000001067 MULES: Solving for alpha.gas MULES: Solving for alpha.liquid gas fraction, min, max = 0.4999968323 0.009997888113 0.9900010559 liquid fraction, min, max = 0.500003201 0.01 0.9900021119 Phase-sum volume fraction, min, max = 1.000000033 1 1.000001067 MULES: Solving for alpha.gas MULES: Solving for alpha.liquid gas fraction, min, max = 0.499995232 0.009996832338 0.9900010559 liquid fraction, min, max = 0.5000048013 0.01 0.9900031677 Phase-sum volume fraction, min, max = 1.000000033 1 1.000001067 smoothSolver: Solving for NH4+.liquid, Initial residual = 1.484075337e-05, Final residual = 1.159761323e-16, No Iterations 2 smoothSolver: Solving for NO3-.liquid, Initial residual = 0, Final residual = 0, No Iterations 1 smoothSolver: Solving for O2.liquid, Initial residual = 0.2795617027, Final residual = 6.228131728e-16, No Iterations 2 smoothSolver: Solving for H+.liquid, Initial residual = 0, Final residual = 0, No Iterations 1 Constructing momentum equations Tf.gasAndLiquid: min = 299.9947355, mean = 300.0339633, max = 300.0732112 DILUPBiCGStab: Solving for e.gas, Initial residual = 0.9999999986, Final residual = 7.01286441e-16, No Iterations 1 DILUPBiCGStab: Solving for e.liquid, Initial residual = 0.1121489206, Final residual = 1.679510257e-16, No Iterations 1 gas min/max T 299.9857036 - 300.0249711 liquid min/max T 299.9714439 - 300.032703 DICPCG: Solving for p_rgh, Initial residual = 0.9950749712, Final residual = 0.0007476549794, No Iterations 58 DICPCG: Solving for p_rgh, Initial residual = 0.01716581647, Final residual = 1.54813425e-05, No Iterations 40 DICPCG: Solving for p_rgh, Initial residual = 0.977211124, Final residual = 0.0007350869486, No Iterations 58 DICPCG: Solving for p_rgh, Initial residual = 0.01014511789, Final residual = 6.818727967e-12, No Iterations 107 PIMPLE: Iteration 2 MULES: Solving for alpha.gas MULES: Solving for alpha.liquid DICPBiCGStab: Solving for alpha.gas, Initial residual = 5.821421844e-17, Final residual = 8.934357652e-34, No Iterations 1 DICPBiCGStab: Solving for alpha.liquid, Initial residual = 4.308481516e-17, Final residual = 7.155231175e-34, No Iterations 1 gas fraction, min, max = 0.4310004415 -8.01486001 2.514951181 liquid fraction, min, max = 0.571470277 -0.01414169716 9.01486001 Phase-sum volume fraction, min, max = 1.002470718 0.9165471615 2.502164998 MULES: Solving for alpha.gas MULES: Solving for alpha.liquid DICPBiCGStab: Solving for alpha.gas, Initial residual = 4.38709213e-17, Final residual = 6.868609435e-34, No Iterations 1 DICPBiCGStab: Solving for alpha.liquid, Initial residual = 4.105585202e-17, Final residual = 6.57320396e-34, No Iterations 1 gas fraction, min, max = 0.3970086899 -84.20635147 74.29030785 liquid fraction, min, max = 0.632777971 -73.29030785 85.20635147 Phase-sum volume fraction, min, max = 1.029786661 0.9005093894 99.99686523 MULES: Solving for alpha.gas MULES: Solving for alpha.liquid DICPBiCGStab: Solving for alpha.gas, Initial residual = 2.799140256e-17, Final residual = 4.373723282e-34, No Iterations 1 DICPBiCGStab: Solving for alpha.liquid, Initial residual = 3.729944455e-17, Final residual = 5.877805245e-34, No Iterations 1 gas fraction, min, max = -7.192486694 -84273.15366 35913.14461 liquid fraction, min, max = 0.7041286022 -35912.14461 84274.15366 Phase-sum volume fraction, min, max = -6.488358091 -37604.99677 580.9315927 smoothSolver: Solving for NH4+.liquid, Initial residual = 0.0112574308, Final residual = 4.168831596e-13, No Iterations 13 smoothSolver: Solving for NO3-.liquid, Initial residual = 0, Final residual = 0, No Iterations 1 smoothSolver: Solving for O2.liquid, Initial residual = 0.2270999707, Final residual = 5.714438306e-13, No Iterations 12 smoothSolver: Solving for H+.liquid, Initial residual = 0, Final residual = 0, No Iterations 1 Constructing momentum equations Tf.gasAndLiquid: min = 299.9711527, mean = 299.9985045, max = 300.0301728

强行存下的第一次PIMPLE循环的结果，下半部应该全都0。
5e042630-92c4-4ba5-ba00-ba9f28ba6b20-1653619214(1).png

请教大家，看了一些论文，为什么大家在研究圆柱绕流涡旋脱落问题时候，是研究不同雷诺数下的涡旋脱落形态？因为我通过数值模拟发现，在圆柱直径一定的时候，只有特定范围内速度的时候才会有涡旋脱落现象发生，在特定流速以外圆柱尾流区域的涡流是几乎不变的没有交替脱落现象。所以一些论文里所说的不同雷诺数下的涡旋脱落形态是既改变了速度也改变了圆柱的直径吗？51eb44d1-55c6-4bfb-ac4a-87d58f2140be-image.png7f0490c7-2424-4dfc-95d0-b73b4cd46b41-image.png 比如图片上这个论文所研究的

@百善孝为先0626 还有一些没有用到的变量也加了注释，然后就成功了，可能删除的那次没有删除干净吧

启动的时候在启动命令后面加上后缀试试

-driver x11

patchInternalField表示你说的这些量

不好意思时间太久远了，我记得我那个算例的后处理的u+和y+的图，做的还挺费事的，是通过paraview导出来的，但是具体是180还是395还是8.89，我已经记不清了

没收敛吧...

:chitang:

好奇老师最后问题是怎么解决的？今天我也遇到了类似问题？

@韬智tz 请问您找到ccm向OF9网格转换的方法了嘛

求推荐OpenFOAM编程方面的资料或课程。
最好是整个流程讲解得比较全面的。

目前想实现自定义格式的流场数据输出。后续打算开发自定义dsmcFoam求解器。

已解决，个人发现在controldict中的 “ writeformat ” 为binary，更改为ascii即可

@shrine 你多试几个，我也是新手，之前遇到过类似的，换几个总有能用的

李老师，我个人觉得公式3.7，3.8的表述有点问题，如果在某一个位置（固定坐标）求导数，那么就是偏导数，而不是全导数，全导数下面是不应该加坐标限制的，也许是理解不太一样

@number44 好的，已经按照这个思路在做了

@bestucan老师我这个算发散了，后续算一个正确的再上传！

@李东岳 东岳老师，我了解了以下结构化网格的排序，然后刚才通过postProcess -func writeCellCentres 写出了channel395的网格，这里网格是按照先x在y后z的顺序排列的，时间步中的结果应该是和这个对应的请问下东岳老师，有办法把结果输出成三维列表来么？

@bestucan 谢谢前辈，版本是一样的，后来和学校超算人员联系发现是少了一个环境配置加载：module load gcc/7.4.4，加上以后就可以了

@冠竹 按道理来说不应该啊 你这个都属于是神bug级别了 如果有时间你就多debug一下 确定是bug可以去官方提bug 没时间就算了这样吧

三台电脑已经通过网线建立了小集群

最开始看到这个话给我整蒙了，我说3个电脑连起来不通过交换机咋搞呢？

就是普通的网线，水晶头的那种，都连到普通交换机上，网卡什么级别我不太清楚，他们的型号都一样

我仔细看了一下，原来楼主通过交换机连起来的。

各位大神请问一下，我最近在模拟气泡与自由液面之间的影响，用的Depthcharge2d改编的，我想问我设置了初始压力和初始半径，还定了气泡位置，怎么得到最大气泡半径呀，还有有没有什么什么办法知道计算结果中的每个时间步的射流速度和气泡半径，像下面图所示。
微信图片_20220523110946.jpg 微信图片_20220523111112.jpg

@石油cfder 理论上越小越好，但也不一定，可以参考文献的设置或者文献中的结果来选择

@fubianhanshu
在particleIO.C下添加一行代码，注释原先的输出可以输出颗粒全局坐标。

os << position() << token::SPACE << celli_;

先采用最笨的方法，把保存时间步设置小，遍历读取坐标，其他的再慢慢实现，看着lagrangian粒子程序太复杂了:136:

好的，非常感谢

@浪迹天大 请问下这个怎么使用啊

管长375mm，管径25mm，叶片厚度2.5mm,结构如下：
static mixer-six element.png

局部尺寸设置
tubewall和yepian的尺寸一样
局部.JPG

表面网格设置
表面.JPG

边界层设置
指定在管道壁面上画边界层
边界.JPG

体网格尺寸设置
体网格.JPG

运行后出现错误
错误.JPG
提示：
Error at Node 1: SIGSEGV(segmentation violation)

Error at Node 0: SIGSEGV(segmentation violation)

Error: Poly meshing not completed.
Error Object: #f
done.
Warning: Deleting empty objects (geom-1)

Error at Node 1: SIGSEGV(segmentation violation)

Error at Node 0: SIGSEGV(segmentation violation)

请问有大佬遇到过这种情况吗？可以给点建议或啥不？提前感谢~~
不在管道壁面添加边界层，网格生成顺利。

@队长别开枪 您好，我在OF8中编译了ccm26ToFoam，想模拟波物相互作用，转化后的网格checkMesh也是ok的，转化之后的网格我放到OF v1912里面进行计算的(还没编译ccmToFoam)，但是相同的设置用OpenFOAM自带的sHM划分的网格可以计算，ccm26ToFoam转化的网格却在0.02s就突然发散了。想问一下有什么解决办法吗？:135:

sed -e 's|^mirrorlist=|#mirrorlist=|g' \ -e 's|^#baseurl=http://mirror.centos.org|baseurl=http://mirrors.tuna.tsinghua.edu.cn|g' \ -i.bak \ /etc/yum.repos.d/CentOS-*.repo yum clean all && yum makecache

@李东岳 @bestucan
感谢两位老师，我大概明白了，马上就去详细看一看:xiexie:

求助各位，目前of的动网格技术可不可以实现下图功能：在计算域顶端和计算域内的物块之间连接弹簧，同时给整个计算域施加简谐运动？
9dbaecc8-35ff-49cb-a874-3db65bdce331-1653113557(1).png
目前学习of五个月有余，从我的了解来看，of自带动网格库中的弹簧只支持一端连接固定点，一端连接运动的物体。想问问大家有没有做过类似的东西，能提供些思路，不动底层代码能够解决这个问题，不胜感激。

paraview 的操作可以写 python 脚本。直接写比较麻烦，可以宏录制之后复制粘贴改参数

tools -> start trace

@coolhhh 感谢回复！

@李东岳 李老师，我用的是extend中的fsifoam求解器，在open foam中画出网格用extend计算，会这样报错。
433cf721-c2ed-4e1f-8490-c3cbeba84ec8-ac9b3169084b49719ca8976178b2dc4.png
但是我的几何边界中没有命名global的边界。请问这是什么问题呢。

@coolhhh 厉害厉害，我就觉着是变量定义格式问题，但是自己不太懂。感谢～

现在高版本只有Schnerr-Sauer模型和ZGB模型

@yfleung 请问下用paraview该怎么做流向和展向的平均呢？

请教一下各位，我做的圆环柱的传热，流向为周期性边界条件，内壁面恒定热流，外壁面绝热，内外壁面无滑移，fvOptions加了一个meanVelocityForce，（0.7 0 0），0.7就是流向，可是最终的速度云图，产生了一块速度较大的区域，并且该区域还会沿环向流动

如图，可以明显看到红色高速区在缓慢的逆时针移动，请问下为什么会产生这个高速区，以及为什么高速区会移动？

压力与速度是耦合的 重力是个源项啊 没法取梯度吧

最近老师让做一个垂直轴风力机的模拟，但是没做过类似的，想请问一下各位前辈，怎么在运动中设置旋转之后，增加一个旋转时的阻尼用来模拟有负载下的风机转动,麻烦前辈们了

@zhe :146: :146:

@hhh 解决了，需要现在ICEM中设置周期性边界，之后就可以了
我参照的是https://blog.csdn.net/muyangzixue/article/details/108631943

@tens 感谢您的提醒，我再去看看

@星星星星晴 感觉还是不太对，网格需要再细一些。

把 return 挨个屏蔽，看看屏蔽了哪个报错。。。:yes:
就知道调用哪个了

我想知道它是什么意思呢？

好像是用 u 和 dim 初始化了一个新的 fvVectorMatrix 类型的变量

@李东岳 东岳老师您好，感谢您的回复。我重新安装gcc 9.4.0和gsl 2.7之后，使用root账户尝试编译发现错误，如下：

make[1]: *** [/home/user1/OpenFOAM/user1-v2012/platforms/linux64GccDPInt32Opt/bin/faceSetToSTL] Error 1 make[1]: Leaving directory `/home/user1/OpenFOAM/user1-v2012/applications/utilities/waves2Foam/applications/utilities/preProcessing/faceSetToSTL' make: *** [faceSetToSTL] Error 2

想询问一下这种错误是来自 gsl 版本问题吗？或者是否需要重新编译OpenFOAM-v2012（之前使用的gcc 4.8.5 编译的v2012）？

在高版本Paraview中我通常是采用contour读取界面，然后选中该表面save data，进行输出界面上的各种数据，但是由于实验室的服务器过于老旧，只能使用paraview4.4版本，同样的操作在输出数据时，得到的数据文件及其繁琐乱糟糟的，一个时刻下竟然生出了5个文件，请教一下大家有没有更好的方法简化输出数据的方法

你有热通量输入到计算域里，但是没有边界让输入的热量离开，系统的温度无法达到统计稳态那就自然随着时间不断提高了

@vien锁 您好，dsmcFoam输出哪些流场文件是由什么控制的。注意到像icoFoam这种求解器是通过fvSolution来设置的，但是dsmcFoam好像并不是。

而且我在使用dsmcFoamPlus时发现，提供的几个案例输出的流场参数也是不一样的。按理说相同求解器输出的哪些流场文件应该是相同的，除非设置了字典文件，然而我在字典文件中并未发现这种设置。

@hellofoam 请问一下,自动画图的温度范围区间不合适，是300-1200，我想把他调整到800-1200该怎么整？

@zhanghan 您好，我最近在用buoyantPimpleFoam做管道的传热，可是温度分布云图却几乎和层流的一样，请问您方便把您之前的算例文件借我借鉴借鉴么？或者说一下您的温度是否出现过我这种问题以及如何解决的？