openfoam添加湍流入口方法?

[sjlouie91]

@cccrrryyy 请问一下，你在安装到OF4时是直接应用的他们的代码还是修改后才可以在新版本上编译的？

[cccrrryyy]

@sjlouie91 可能改过吧，不记得了。即使是改过应该也就是很小的地方改了一下，按照编译时候出现的错误提示改的，没有额外的东西。

[dyb]

@cccrrryyy 老师，您好，您这边有InflowGenerator的算例吗，这个值我不知道怎么设置

[cccrrryyy]

@dyb 好久不登录了- -它需要的就是速度、长度尺度和雷诺应力。简单的话速度你给个均一值，长度尺度取水力直径的10%，雷诺应力（RField）我一般都按照 k epsilon 里面 k 的算法（3/2*(refU*turbulentIntensity)^2，turbulentIntensity取个百分之五啥的），RField格式是 uniform (Rxx Rxy Rxz Ryy Ryz Rzz)，我一般就把算出来的这个值给Rxx Ryy Rzz。这个雷诺应力的顺序你自己再查一下，应该是这样的。

边界条件的设定其实没有很重要，大致差不多就行，比如长度尺度我看有人取水力直径的7%的，turbulentIntensity取到百分之二十的。雷诺应力，严格来说其他分量也应该有值，但我一直就给Rxx Ryy Rzz好像也没啥问题。

这个边界条件只是给一个扰动，真正的湍流应该是靠你场里面的流动产生，所以在边界条件设置上不要过于纠结具体的数值。

[dyb]

@cccrrryyy 谢谢老师

[dyb]

@cccrrryyy 老师，这个速度可不可以给定一的速度梯度

[cccrrryyy]

给速度梯度是什么意思？

[dyb]

@cccrrryyy 我的入口是大气边界，有一个风廓线，入口速度沿着高度方向有个速度梯度，按指数变化

[cccrrryyy]

@dyb 那就是速度有一个profile？我一般的做法是直接把这个profile用fixedValue赋值到每个cell上。这个需要你先用writeCellCentres命令导出每个边界的坐标，然后找到入口对应的坐标。我一般拿到这些值后用MATLAB根据位置坐标把速度值赋过去。

OpenFOAM是有codedFixedValue这种东西的，就可以直接给一个profile的，这个应该操作比较容易。

[dyb]

@cccrrryyy 在open foam中atmBoundaryLayerInletVelocity这个边界是带有风廓线，而InflowGenerator只有扰流没有风廓线，如果这两个边界可以编译到一起用作入口是非常不错。
atmBoundaryLayerInletVelocity这个边界的源码在OpenFOAM-2.4.0\src\turbulenceModels\incompressible\RAS\derivedFvPatchFields\atmBoundaryLayerInletVelocity这个位置，我虽然在用openfoam但是我不懂编程，我解决不了这个问题
非常感谢老师的解答，

[星星星星晴]

@dyb 和你说一个更硬核的方法，对进口的每个cell进行赋值 你可以参照这个，这个方法最大的问题就是你改变一次网格就要重新赋值。但是当你在为边界条件搞得头疼不已，就用这个吧

http://www.cfd-china.com/topic/2089/分享-搬运-自定义非均匀inlet-u

[dyb]

@星星星星晴 谢谢老师

[cccrrryyy]

@星星星星晴 我上面回复的做法跟你这个一样，你这个很详细~我觉得这样自由度更高一些，实际上写好一个赋值速度的程序之后可以一直用，很快很方便。

[dyb]

@cccrrryyy atmBoundaryLayerInletVelocity这个就是一个可赋值的，这是他生成的场很方便，就是不懂怎把这个入口碧边界和InflowGenerator入口边界结合起来，

[cccrrryyy]

@dyb 把两个结合起来肯定是有办法的，就相当于自己写一个边界条件了。但是我感觉这样需要的代码工作量可能很大，因为两部分代码里面都包含一些实际上你不想要管的东西，比如什么构造函数啊之类的，还要去好好的想每个类里面哪些数据和函数应该做成protected，哪些做成public之类的，种种问题。除非你是一个编程高手很熟悉这些东西，不然真的可能工作量太大了。实际上对inflowGenerator来说你就是需要给一个速度的profile，我觉得这个profile你按照之前那个post里面的方法去给就好了，根据位置去确定每个cell的速度是多少的时候你用atmBoundaryLayerInletVelocity里面提到的那个公式即可。

[dyb]

@cccrrryyy 谢谢老师

[dyb]

@星星星星晴 谢谢老师，我按照老师的方法已经弄出来了

[zxx]

@cccrrryyy https:Eugene de Villiers的扰动在这里//github.com/wyldckat/perturbU

[cccrrryyy]

@zxx 感谢！

[浪迹天大]

我把代码移植到OF-上了，并附上了我的算例以供参考：
https://github.com/ZhangYanTJU/boundaryConditions