CFD中文网

各位好，我最近用OpenFOAM也在做这方面的模拟，现在使用SHM划分网格遇到了一些问题，主要是网格质量、非正交性、凹度和负体积网格方面的，能否交流一下

算不准大概率是网格问题，不仅仅是数量，还有质量，另外，系数值的获取是如何计算，自己手动计算还是直接通过软件监测？若通过软件监测，则需注意参考值的设定。

参考文献
Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.

@李东岳 模型缺陷吗.....那就有点麻烦了。我不单纯写文章，是老师的项目。我再考虑考虑

@李东岳 多谢李老师鼓励，还得感谢李老师的讨论答疑和提供的交流平台:ok: 。

@evensun 嗯，涡粘系数是二阶张量这个按我理解其实就是非线性湍流模型了

目前这个课题做了半年了，纯粹的是算法类研究，网格1000足矣，1维就能跑。今天下午和3个大老板碰头了一下，聊了一个小时，很遗憾三个大老板一个没懂。

这个算法在8月份就已经植入好了并且和解析解对比完美符合。不过即使是我自己都觉得这个算法太复杂。毕竟大老板们肯定不会去太关注算法的细节，但是这个算法有点挑战传统的观念，9月份给他们讲了一次，没懂，搁置了。。。今天又给他们讲了一次，还是没懂。

于是写了一下午的文档，打算邮件再给他们说一说..

TIRED

世界上最难搞的算法，就是连你老板都不知道你做的是啥

...

如果你认为动量方程增量很小，可以尝试把动量方程的松弛因子设为1，看看稳定性:duang: 好像问题不是这么简单

@李东岳 是的，现在bug还是比较多，估计迭代个一年效果会不错

@ycx1042341217 @huningdong 请问两位，壁面第一层网格高度可以做到什么程度？我发现最小只能设置到大约5e-5m:136:

按照文档里的说法，有时换一下meshShrinker也许会有效果，感觉文档里那个例子和你的也有点像。

addLayersControls { ... meshShrinker displacementMotionSolver; solver displacementLaplacian; displacementLaplacianCoeffs { diffusivity quadratic inverseDistance 1(wall); } ... }

相应的还需要调整下fvSolution和fvScheme，具体文档里都有。

话说新版文档虽然变好看了，但图居然是糊的，可能还需要参考下旧版文档。

另外这里也请教下各位老师，按我的理解snappyHexMesh做完snap之后，在layer这步应该是从几何表面出发，向外挤出来一块空间把边界层塞进去。怎么在楼主的这个例子中反而是向内把几何都挤变形了

@李东岳 老师我OF-4.1版本算例constant文件夹下得有turbulenceProperties字典，里面写的是laminar，不涉及湍流。好像OF-8以后都不需要turbulenceProperties字典了。

:136: 好吧找到原因了，是这个版本的link edge factor默认是0

@李东岳 好的博主 但是为什么计算后在paraview里的质量分数大于1呢

IMG_20231127_193223.jpg
在计算单相流作为初场，可以收敛，但是湍流粘度比超限，继续设置为多相流空化计算的时候，开始发散，不知道有哪些解决办法，已经重新画过网格。
边界条件是速度入口，压力出口，湍流sst k-w，稳态计算

@如果随风 pointswise可以画

@bestucan hhhh谢谢您，我发现我哪里错了，我画网格的时候是xoz面上的，正确的应该是在xoy画个2D的网格，然后openfoam会在z轴增加一个面。

0_1510344430862_草图.png

这本书我2015年买的，花了118欧，还让外导签了名。今天双11，琢磨亚马逊有没有什么书可以买，忽然发现这本书65欧了，便宜了一半啊！

这几年买国外书估计都花了1000欧了 :angry: 因为喜欢写写画画，打印的话A4纸太重太厚。

@ibelief 谢谢，我再看一看。