改变每个网格的粘度



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


  • 管理员

    听起来不太符合物理:chigua2:



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


  • 副教授

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



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


  • 副教授

    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里面做过,稳态情况下(流体流过一个圆柱)收敛效果还是不错的,守恒性也挺好。瞬态的时候有点小问题,算是没问题的,也没去细究。我没做很深,所以没有什么成体系的方法论或者结论。
    我个人其实还蛮喜欢这种能解决工程问题的简单处理方式,更复杂的沉浸边界法也有很多,有的也做得很精细。



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


  • 副教授

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



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


  • 副教授

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


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1