改变每个网格的粘度

[Samuel-Tu]

OF里面能够对每个网格赋予不同的粘度么？不知道该调用什么

[李东岳]

听起来不太符合物理

[Samuel-Tu]

@东岳 我发现用UEqn.setValues，把一个区域的速度设为0，但是求解后，那个地方的速度并不完全为0，只是比较小。我觉得可能simple算法中，U最后是用p重构的，UEqn只是提供了预测速度。
所以我想通过改变粘度的方法，设置一个很大的粘度，使得那个区域速度完全为0。这需要对每个网格的粘度都进行控制，所以很麻烦。。不知道调用什么函数。因为我看很多有关于nu的都是const，改不了。。

[cccrrryyy]

@Samuel-Tu 能不能直接在每一个迭代步最后强行把速度设成0呢？之前在FLUENT里面植入沉浸边界法的时候这么弄过。

[Samuel-Tu]

@cccrrryyy 应该可以，只是不知道收敛性如何。你用的direct forcing 的浸入边界法吗。为什么会直接把速度设为0呀。
有没有发表论文我参考下( ˘•ω•˘ )，我也在做浸入边界法。。

[cccrrryyy]

Interaction between particle clusters and particle-induced turbulence (International Journal of Heat and Fluid Flow, Volume 23, Issue 5, October 2002, Pages 639-646)，Unsteady three-dimensional simulation of interactions between flow and two particles(International Journal of Multiphase Flow, Volume 29, Issue 9, September 2003, Pages 1431-1450)。
记得没错我当时参考的是这个，不过它里面是用不同相的速度差来构建流场的源项，然后加入动量方程。我当时想，既然加入这个源项的作用是把计算域的一部分变成想要的速度（固体相体积分数的线性函数，全是固体相那就是0），那不如直接在每一步最后强制让速度等于想要的速度。在FLUENT里面做过，稳态情况下（流体流过一个圆柱）收敛效果还是不错的，守恒性也挺好。瞬态的时候有点小问题，算是没问题的，也没去细究。我没做很深，所以没有什么成体系的方法论或者结论。
我个人其实还蛮喜欢这种能解决工程问题的简单处理方式，更复杂的沉浸边界法也有很多，有的也做得很精细。

[Samuel-Tu]

@cccrrryyy 噢噢，这样啊啊。我应该试过。对于圆柱绕流就是把圆柱内部的速度全部设为0。但是我现在做的是高雷诺数下的屋面剪切力。我发现直接设为0，剪切力不太对，太小了。。

[cccrrryyy]

@Samuel-Tu 完全在圆柱内部的网格点设为0，但是圆柱边界跨过的那些网格处的速度是需要重构的吧？

[Samuel-Tu]

@cccrrryyy 我试过重构，发现壁面剪力曲线毛刺非常多。。。有可能是网格比较大。。

[cccrrryyy]

@Samuel-Tu 原来这样，简单的重构方法好像确实有一个代价就是网格要比较细，不然不准确。