openfoam添加湍流入口方法?

[dyb]

@星星星星晴 您好，您这边有没有关于OpenFoam 求解器详细介绍的资料呢？谢谢您

[星星星星晴]

@dyb 没有，你可以google一下

[cccrrryyy]

@dyb 这个提示你的Foam::Random类里面没有叫做vector的函数或者参数，你需要搜一下源代码了。我记得我编译的时候没这个问题，不过我用的是4.x。可能不同版本之间它改了名字，一般这种随机数产生器里面肯定会有叫做什么vector之类的函数给你调用，比如这个出错误提示里面的ranGen.vector01()，我猜大概率是产生一个0到1之间的随机数？

[dyb]

问题找到了,低版本的Random和高版本的有区别，但是，我不知道怎么 改

[cccrrryyy]

你看看它哪里用到了这个vector函数，搞清楚用这个是想实现什么，然后在你现有OpenFOAM里面的随机数产生器找一下实现同样功能的函数。当然了，前提是你现有的OpenFOAM里面确实包含了他要用到的功能，这个我没法百分百确定，但这种涡方法里面产生随机数大概率就是调用一个函数然后产生一个0到1之间的随机数。

或者更直接点的，这个湍流入口好像是针对2.3.x（不一定准确）开发的吧，你直接找到2.3.x的源码，看一下这个湍流入口调用的是2.3.x版本里面的哪个函数。这样你就知道在你的版本下应该调用什么函数了。

[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里面提到的那个公式即可。