摩擦速度比论文偏小



  • 各位CFDer,我最近做了一个高低屋面的流场计算,流线如下图:
    jietu.png
    计算稳定后(低屋面上的分离点位置不再变化),提取低屋面上的摩擦速度,发现摩擦速度比论文里说的偏小,但分离点的位置(摩擦速度为0的点)对应上了。
    论文摩擦速度如下(看u*=0.8的线):
    Screenshot from 2020-01-15 13-53-43.png

    我算出来的摩擦速度如下:
    Screenshot from 2020-01-15 14-00-09.png
    想不通是什么原因,我的yPlus如下:
    Screenshot from 2020-01-15 14-02-49.png
    使用的realizedKE(和论文一致,但论文使用的是FLUENT,未说明壁面模型),我的k壁面使用的是kqRWallFunction(高低雷诺数都适用),我的epsilon壁面使用的是epsilonWallFunction(高低雷诺数都适用),我的nut壁面使用的是nutkWallFunction(高低雷诺数都适用)
    请大家帮忙分析下原因



  • 另外,我的摩擦速度是用壁面剪切应力wallShearStress的X分量开平方计算出来的,我的模型是不可压流体,因此没有除以密度,这个计算应该没问题吧。



  • @Samuel-Tu摩擦速度比论文偏小 中说:

    kqRWallFunction是零梯度,应该只适合高雷诺数,kLowReWallFunction才是高低雷诺数都适合。epsilonWallFunction应该也只适合高雷诺数,epsilonLowReWallFunction才是高低雷诺数都适用。



  • @我是河滩 kqr确实是我记错了,不过epsilon我印象中已经合并了,只有一个epsilonwallfunction了


Log in to reply