输出壁面摩擦速度

Samuel-Tu

OpenFOAM的postProcess有计算壁面剪力或者壁面摩擦速度的命令么？因为我的计算已经完成了，想看看壁面摩擦速度，所以需要后处理一下。在paraFoam中利用Cmu^0.25*k^0.5这个公式计算了摩擦速度，但是分布明显不对，我觉得是这个公式过于经验了。所以我想postProcess输出下壁面剪力或者摩擦速度看看。
望大家不吝赐教。

benqing

有一个叫wallshearstress，2.3.0里面可以直接输入这个命令，不知道高版本会不会不一样

benqing

准确的应该是wallShearStress，计算的是壁面剪切应力tau_wall

Samuel-Tu

@benqing 感谢

Samuel-Tu

成功了，在system里加一个文本文件，包含如下代码即可：

functions
{
    wallShearStress1
    {
        type  wallShearStress;
        libs        ("libfieldFunctionObjects.so");
    }
}

Samuel-Tu

wallShearStress倒是输出出来了，但是有两个问题

1. wallShearStress是一个vector，按理说就应该放心应该是切向的（如X向），但是我发现X和Y向都有数值，说明壁面剪力不是切向的，让我很奇怪
1. 输出后能否可视化？本来想导入paraview然后作图，但是发现没有选项，如何可视化呢?

Samuel-Tu

@Samuel-Tu 或者不加这个文件，直接用命令

pisoFoam -postProcess -func wallShearStress -noZero -noFunctionObjects

Samuel-Tu

@Samuel-Tu 可以可视化，可以导入paraview

vbcwl

@benqing 请问下这个是计算utau=sqrt(tauw/密度)里的tauw么？我求出来怎么全是负数