openfoam中添加湍流激励

[Rachel0096]

请教各位老师，参考了相关文献，我想在动量方程中添加如下的湍流激励项（方程最右边）：

其中A为湍流forcing参数，k为瞬时计算的湍流动能，k0为期望的稳态湍流动能。请问这里需要怎么添加代码？

[李东岳]

@Rachel0096 最简单的就是直接去动量方程里面加个源项。

fvVectorMatrix uEqn
(
  ...
   ==
  A*(k0/k)*U
);

[Rachel0096]

@李东岳 请问李老师，这里添加的A，k0，k需要去哪里声名或者定义呢，是在求解器的C文件当中吗？

[李东岳]

@Rachel0096 是的。在求解器层面去定义。

[杨英狄]

@李东岳 请教一下李老师，怎么在creatFields或者求解器中声明湍流的相关信息呢，如湍动能k，谢谢李老师

[李东岳]

@杨英狄 你要在求解器里面调用湍流变量么？ openfoam这面一般把湍流相关的剥离出来了。一般我们不在createFields里面去操作。你要做什么？

[杨英狄]

@李东岳 李老师您好，我在做空化，需要考虑湍动能对空化的影响，其实是在饱和蒸气压上加一个修正项：pturb=0.39*（密度*湍动能）。所以就需要在求解器里面声明一下湍动能。
感谢李老师。

[李东岳]

@杨英狄

pturb=0.39*（密度*湍动能）

类似这种，可以通过这个实现，不需要复杂的操作

const volScalarField& k = mesh_.lookupObject<volScalarField>("k");