wallShearStress使用的一些问题

[刘旭阳]

各位老师好，想请教一个关于wallShearStress的问题。我想要计算一个表面的剪应力，用wallShearStress的代码算出来的只是一个表面的最小值和最大值，但是我想得到此表面上的每个点的剪应力分布。
如：
类似于这样的图，简单说来就是在圆柱绕流的计算上得到表面剪应力。那么如何用openFoam来实现呢？
我用simpleFoam求解器。作为一个地理人，对偏力学的知识确实有些欠缺，请各位老师不吝赐教，谢谢。

[bestucan]

https://www.openfoam.com/documentation/guides/latest/doc/guide-fos-field-wallShearStress.html

wallShearStress 输出就是一个场呀，你用的是别人修改过的 wallShearStress 么

[刘旭阳]

@bestucan 不好意思，老师，刚看到。我用的是您发的那个链接上的原代码，没有进行修改。那如何才能计算得到每个坐标点对应的剪应力呢？我想计算一个面（如长方体计算域底面）上的所有坐标点的剪应力数据，那么这个用openfoam如何实现呢？我用我能想到的办法计算的结果都不对。请老师赐教，谢谢！

[bestucan]

@刘旭阳 我没搞过，不过这个帖子好像和你相干的事一样：帖子

[刘旭阳]

@bestucan 收到，谢谢老师~

[李东岳]

@刘旭阳 wallShearStress是计算不出来你贴的那种云图的。这是“壁面剪切力”，只有壁面附近有。你可以用下面的命令来看wallShearStress

simpleFoam -postProcess -func wallShearStress

[刘旭阳]

@李东岳 嗯嗯，东岳老师，我用wallShearStress计算出来的结果是一个壁面上的剪应力大小范围，并不能计算出我想要的结果。bestucan老师附的帖子我正在研究。请问东岳老师，您有没有过类似的经验呢？我使用了各种方法，都得不到正确的结果，无法与风洞实验进行对比。

[李东岳]

@刘旭阳

wallShearStress计算出来的结果是一个壁面上的剪应力大小范围

为什么是个范围，不是一个值？不理解。

[李东岳]

我跑了一下，理解你的意思了，你是看log的这个输出吧：

No MRF models present

No fvModels present
No fvConstraints present
wallShearStress wallShearStress:
    processing all wall patches

wallShearStress wallShearStress write:
    writing object wallShearStress
    min/max(upperWall) = (-0.470645 -0.0433369 -2.56057e-17), (-0.155955 0.0424452 1.57794e-17)
    min/max(lowerWall) = (-0.613171 -0.129826 -5.78407e-18), (0.0915054 0.0731431 4.48414e-18)


End

不要看这个，看你时间步文件下面的文件，云图

[刘旭阳]

@李东岳 我试了一下，可以的，是我想要的结果。谢谢东岳老师！