@李东岳 :wocao:
大哥，你这样我压力好大= =

@veen 不清楚，cfl只是说保证计算结果收敛，精度的问题可能重点不在这上面了，稳态计算的话cfl要求没那么高

最近也基于pisoFoam植入了一下RKprojection foam，单个计算步不限制残差，projection method的确快很多，但是限制残差之后速度就变得特别慢，并且该算法内存占用量是pisoFoam的20倍。
不知道有没有大佬看到过把这个算法植入到多相流求解器的呢？想尝试一下把这个方法放到多相流计算里边去。

感谢分享！！:146:

十分感谢

@星星星星晴 明白了.我这个方法确实直接套用不上.并且碰撞模型从同一个cell扩展到周围cells之间的碰撞,还会有其他的问题.我再想想,谢谢了.以后常联络!

@CarmeloSun 谢谢前辈

改写了coalChemistryFoam求解器，添加了volume fraction。现在想把pairCollision模型也添加到coalChemistryFoam中，有些一筹莫展，请求大神指点！！！
目前想到2个思路，一是把CollidingParcel加入到coalParcel中，已经试过了，有问题，应该不可行吧；二是改写pairCollision模型，现在卡壳了，不知道该怎么改。

#!/bin/sh # Source tutorial run functions . WM_PROJECT_DIR/bin/tools/RunFunctions # Get application name application="potentialFoam" runApplication blockMesh \cp save/boundary constant/polyMesh/boundary \rm -rf 0 \cp -r 0_org 0 runApplication $application -writep runApplication writeIbMasks

这个是cylinderInChannelPotential中Allrun文件代码，就是比较想知道save文件夹里边有block Mesh，constant文件夹下的ployMesh里也有blockMesh，所以Allrun中的blockMesh命令运行的是哪个文件。然后将save/boundary 拷贝到constant/polyMesh/boundary的意义是什么呢

也是实测。OpenFOAM自带算例。256万网格。虚拟机16核跑了1772s，真系统跑了548s。

虚拟机大约要慢3倍

@队长别开枪 专业 学习了学习了:140: 我看朋友用blender做过一些美术动画，看起来还是挺舒适的:papa: 目测不太适合工程:xinlei:

@同学博 是的，基于 freeCAD，专门给 OF 开发的前处理工具

@李东岳 谢谢李老师！我发现问题了，这个报错其实就是最后的correctADDStress的问题，所以看这个报错提示就是看最后一个信息，错误的原因是我用wclean去移除编译的库的时候其实没有移除，所以后面不管我怎么改，都是调用的原来的库，解决办法就是手动删除openfoam-7的platform中生成的自己编译的库，重新编译就好了:136:

可能是超算上有不只一个节点，所以节点间通过IP地址交流信息。
而linux系统里有一套“权限”和“环境变量”的继承规则。你加载 openfoam 环境变量的的那句应该是放在家目录的.bashrc了（就是source ~/OpenFOAM/etc/bashrc，大概是这句）

而 miprun 再通过IP调用其他节点的时候，这个网络通信的用户权限和环境变量，应该没有继承你的用户名的环境变量；而是系统的环境变量。所以不包含你的用户目录下.bashrc里的东西。

把你家目录那句 source 放到 /etc/profile 试试 。这是系统范围的环境变量。得是绝对路径：/home/<用户名>/OpenFOAM/……

当然也可能得在各个节点上都配上环境变量、装上 openfoam，这细致活得慢慢调……

部分求解器出现上述状况（编译成功，但是不能算），还有一部分求解器完全正常。

没用过多节点的，但问题是出在环境变量上没跑了。
到这个路径下
~/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin
这都是可运行的二进制文件。

直接./<文件名>，比如
./icoFoam
就知道是没编译好还是环境变量没配置好了。

http://c.biancheng.net/view/2226.html

this就是当前对象的别称，有哪些成员，就看它所在的函数是哪个对象的成员函数了。
去定义这个对象的头文件查，与之同名的.c文件有各成员函数的定义。

感谢回复，我研究研究。我把体积分转化成了面积分，然后用python calculator进行了求和。感觉结果不太理想

@dzw05 你遇到的极值局部点距离给定压力条件的边界接近吗？我今天下午算的时候发现时间步过大时，距离fixedValue p_rgh边界最近的cell中心也会出现温度极值....调小了时间步后就没有了。我算的这段时间一直感觉温度和压强的联系太密切了...

不好判断 但感觉就是设置的问题 如果可以上传算例并且网格数量小的话 我可以给你看看

电势是相对的，上下机板电势差不变，就像电池正负两极保持电势差一样。不影响是否绝缘。
电荷密度零梯度，仅在极板与空气交界线上。看这种情况：紧挨着交接线两边各有10个电荷，所以零梯度。但是一边诗歌电荷后面站着100个电荷，另一边十个电荷后面站着20个电荷。

不同于流场的变化，仅分析相邻区域即可。电荷之间的相互作用是“超距”的，不相邻也可以用场去作用。
总的来看，电极板之间已有电势，只要有载流子就有电流。就像在电容间加水，就存不住电势差了。
交界处电荷梯度为零可以两边同时增加或者减少，或者增加减少动态平衡。

如果要绝缘，按照物理意义的绝缘最靠谱：没有载流子。
如果强行绝缘，应该是电势梯度为0，没有电势差，就没有电荷受力。但是电势场是场内所有电荷分布造成的结果，应该不能直接设定。

具体怎么实现得看源代码是怎么算电场的了

@veen 多谢大佬！^O^/