现在明白了，谢谢东岳老师讲解

@东岳 李老师，关掉考虑自转之后出现倒流现象更严重了

@星星星星晴 好的谢谢，我目前网格有7500W，质量0.72以上，结构是简单的圆管，在算DNS，ncorrector设置的2，您觉得可以吗？

大家好，我参考openfoam自带的求解器，在求解器中添加自适应网格的功能。我的算例网格总数是50000，运行算例的时候，遇到下面的错误
1.png
以下是reactMinMax.H 17-30行

17 const label tautableDim = 3; 18 19 forAll(p,i) 20 { 21 List<scalar> lookupCoordinates(tautableDim,0.0); 22 lookupCoordinates[0]=fH[i]; 23 lookupCoordinates[1]=pHI[i]; 24 lookupCoordinates[2]=THI[i]; 25 26 const List<scalar> &LOG_t_ign_temp = tignTable.lookUpBounded(lookupCoordinates,false); 27 tIgnHI[i] = min(max(1e-14,pow(10,LOG_t_ign_temp[0])),1e9); 28 29 } 30 tIgnHI.correctBoundaryConditions();

请教大家问题在哪呢

请问PBiCGStab在更稳定，收敛速度更快的情况下，会不会牺牲精度呢？还是两者都会给出正确结果？我目前正在算DNS，对于计算准确性比较关心。谢谢。

大家好，问一个很弱智的问题。。。用了openfoam的sampling surface，输出格式为raw。请问用什么软件处理呢？

此贴基本完结，按着源代码琢磨了一下午，基本弄清楚了是如何计算的。nearestPoint()主要使用了实时碰撞检测算法，可以计算出网格中心点到一个面的最短距离。这个最短距离是通过寻找面上离网格中心点的最近点来计算的。等我过段时间总结好了，到时候再补充到这个贴子下面。

@东岳 谢谢李老师！我去试试

@anubis 哦那我记错了，talphaPhi1Un是高阶的，alphaPhi10是低阶的。

在头文件找到了描述：The low-Re correction is activated by setting the entry lowReCorrection to 'on'; in this mode the model switches between laminar and turbulent functions based on the laminar-to-turbulent y+ value derived from kappa and E. When the lowReCorrection is inactive, the wall function operates in high-Re mode.可能是在constant文件里使用这epsilonWallFunction的时候能够选择打不打开LowRe

之前用interFoam，VOF做液滴在气体中的碰撞（湍流LES模型），现在想做水滴在油里的碰撞，除了要改油的物性参数（transportproPerties）以外还需要改哪些文件？
另外，我想对比层流中液滴碰撞和湍流中的液滴碰撞，但是没有看出区别是什么原因

https://github.com/OpenFOAM/OpenFOAM-2.3.x

有区别，最好下载这个官方版

if (Pstream::parRun()) { if (Pstream::master()) { Info << "balabala..." << endl; } } else { Info << "balabala..." << endl; }

这样串行并行就都OK了，可能需要#include "Pstream.H"。希望能帮到你。

@麦迪文 那您觉得应该从哪方面入手，小白我还没有改过边界

equation.png
大气边界层风力机的控制方程为上图（比较前沿的文章基本都采用上式）。我觉得OpenFOAM中的boussinesq浮力项与文章的并不一致，因此我并不敢直接用，有知道其中原理的朋友可以解释一下吗？假如要自己写一个浮力项，<potential temperature>这个带了“<>”的势温按照文献中的意思为"horizontal average"，在求解器应如何实现？是同一水平面的势温取平均还是邻近网格取平均？

请问大家 我用如下的并行命令
mpirun -np 20 mapFields sourceDirectory -parallelTarget
但是在执行该命令时出现下面的情况
Screenshot from 2019-10-22 12:49:34.png
但是map过程没有结束, 请问这种情况如何解决?还是说并行命令有问题?

@Hungryandfool 在 有没有哥们做破碎波浪的，可以交流下么？我用REEF3D做，可以和OpenFOAM的哥们讨论 中说：

2106936

OK,

@东岳 多谢李老师了，是时间步长的原因，我开始按库伦数取得 听了您的建议往小了试，没有出现发散

@东岳 感谢感谢

那位OpenFOAM经验丰富的前辈能讲讲 OpenFOAM的算例之间如何耦合？