CFD中文网

0_1460680809776_捕获.JPG

@李东岳 在 有关侵入边界法。 中说：

Hi，

OpenFOAM-extend里面有IBM方法，Extend和官方的OpenFOAM侧重不同，几个我了解到的不同主要是IBM，粘弹流体，动网格表面跟踪。

还有可压缩流动高分辨率格式。https://github.com/Unofficial-Extend-Project-Mirror/Unofficial-Extend-Project-Mirror.github.io/wiki

@搬运工不好当

这种情况很常见，如果使用第三方网格生成程序比如ICEM做的网格需要使用periodic严格定义才可以。如果你确定你的面试匹配的话，可以加大tolerance来忽略这个提醒，也不会有错误。这就是网格的原因。

这里的errors是误差的意思，而不是错误。可参考cfd-online上的这个帖子。

omegaWallFunction 这个边界条件继承自 fixedValueFvPatchField类，其构造函数如下
0_1460115657665_upload-47e7e05b-ec74-4419-8ff9-d55730af0dfa

这个构造函数里，需要注意两点，第一是传给基类fixedValueFvPatchField的参数，也就是 fixedValueFvPatchField<scalar>(p, iF, dict)，这一句里，基类会通过参数iF 去读取 类似 value uniform 2; 这样的语句，来讲当前边界的值初始化为 2。不过，看下面最后一句，this->operator==(patchInterField());，这里将当前边界的值重新赋值为临近边界的第一层网格中心的值。

所以，综合起来，value uniform 2; 这一句不起实际作用。

@cfd-china 我觉得不是，并不是大于$\Delta$的窦唯大涡。

ｏｐｅｎｆｏａｍ　边界条件引入到离散控制方程求解中，实际是把边界节点的值引入到了源项来看待的吧？最想不明白的就是再靠近壁面的体积单元Ｐ，再其通过惯性力/压力差/剪切力形成动量方程的时候，是不是还要把壁面剪切应力考虑进来。

@搬运工不好当

CFD界：相对于商业代码，OpenFOAM的特性是什么？

Henry：OpenFOAM是完全免费并且开源的（CFD界注：开源的意思即为使用者可以随意更改软件的内核）！这提供了HUGE的方便，不仅仅在费用上，也表现在灵活性、以及从软件本身来学习CFD上。OpenFOAM不仅仅让CFD工程师可以模拟算例，它也允许CFD工程师来学习CFD方法、开发CFD算法并完善，甚至某一天创造一个可以替代OpenFOAM的软件。
并且，许多CFD问题需要一些特别的功能、算法、方程、边界条件等。在商业软件中植入这些功能可能令人沮丧并且很可能失败。有些情况下可能有必要完全的重新写一个求解器、边界条件等。OpenFOAM提供了最大的灵活以及扩展性。这些对使用CFD做研究的人非常重要，因此我觉得OpenFOAM可能在某一天会替代闭源软件。

@搬运工不好当 说：

求速度方程之后为什么没有求压力方程，而

Hi，

我觉得能量方程放在压力方程之后和之前结果是相同的。就像湍流求解，标量传输放在压力方程之前和压力方程之后也是相同的。

你可以自己测试一下并反馈？

@wwzhao 在OpenFOAM-2.1的时候，程序里写的就是“CrankNicholson” ，所以controlDict里也应该写“CrankNicholson”，虽然这个拼写确实是不正确的。但是从OpenFOAM-2.3开始，这个名字改成了正确的拼写：“CrankNicolson”，相差一个字母。所以这就导致了2.1和2.3的算例会在这里不兼容。

@搬运工不好当 这种问题不是特别想回答了，监控你关注的量即可。 设置动画和自动保存。

我也是相同的问题啊！！！

0_1462072825790_14.jpg
调大此值，如果还出现警告，建议改为标准初始化

你得把错误内容贴出来，不然一句话看不出来。出错的原因太多了，而且不知道你修改了求解器，甚至case的哪部分...............

@搬运工不好当 说：

OpenFOAM2.3.0中有buoyantPressure这个边界条件吗?没有的话要用哪个边界条件来代替 buoyantPressure呢？

OpenFOAM 2.3.0 中没有buoyantPressure，替代的边界条件为fixedFluxPressure，详细可参考cfd-online的这个帖子。

Hi，

这和你的求解方程有关，例如：

fvm::ddt(rho) + fvc::div(U)
和
fvm::laplacian(nu, p) = fvc::div(phi)

是对称矩阵。更进一步的，第一个是对角阵。

fvm::ddt(U) + fvm::div(U)

则是非对称矩阵。

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

Hi，

我遇到的使用三菱柱的情况就是四面体的边界层网格。当然需要使用三棱柱来做边界层，有利于使用三棱柱的上下面把误差抵消掉。

计算效率是指？
目前出了纯结构网格，其他的网格在做索引的时候编程都要比结构网格复杂。
如果值得是精度，三棱柱在边界层肯定要比4面体好，这个在教材中都有提及。

“\”好像是Linux的语法。 我的理解，是换行的意思，表示后面还有代码需要执行。你看最后一行就没有""了

你这个是三维的LES单相流模拟，应该不难分析，常见的不收敛的原因有很多，上面大神们说的时间步长、网格质量、松弛因子都是很不错也很常见的调试首选项，还有就是边界类型也挺关键，再就是几何域合不合适，上下游延伸段够不够，比如存在很狭小的流动域之类的、动网格更新策略，都可以作为改进收敛性的切入点。