能问下关于湍流方程是否可以添加源项吗？

mohui

最近想添加一个湍动能的源项，不知道是否可以直接在k方程添加源项呢？

wwzhao

可以用 fvm::susp 添加。

mohui

@wwzhao 准备怼of湍流模型的调用机制，导致扫了一下K-equation 是可以添加源项的，但是具体的还不是很清楚，比如说如何引入源项（因为源项是通过某个数学模型计算得到的）。

wwzhao

@mohui 打开源码，找到kEqn，添加源项，完成。

mohui

@wwzhao 谢谢，我看到这个方程了，现在我在想如何引进源项。

wwzhao

@mohui 你可以把源项 IOobject 注册到内存中，再在湍流模型里通过 lookupObject 查找得到源项。

mohui

@wwzhao 好的，谢谢。现在稍微了解到3.0以后版本的湍流模型的编译较为复杂，需要仔细看看。

zym604

用pimpleFoam的话可以直接在fvoptions里添加momentum source项。

wwzhao

@zym604 动量方程和湍动能方程是两个不同的概念，fvOptions只能在动量方程中使用。因为只有动量方程的源项中加了fvOptions(U)，而湍动能方程没有。

zym604

@wwzhao 啊，sorry，没看清题目。

yuan_neu

@mohui 也不复杂吧，比如对于kepsilon，用incompressible turbulence模型，
1 在ras目录下直接在kepsilon做更改，
2 回到incompressible目录下，编译
3 回到求解器，在编译一下
ok 还算挺方便的
当然也可以添加新的文件 mykepsilon，需要在incompressible目录下修改某子文件夹下的 xxxmodels 添加两句代码，就可以了

yuan_neu

@wwzhao 其实很多模型里面已经有ksource项了，当然默认下这项为零，可以直接将这个source给激活就行了。虽然我没用过

Cloud1990

of5.0里面已经有fvOption了，

    tmp<fvScalarMatrix> epsEqn
    (
        fvm::ddt(alpha, rho, epsilon_)
      + fvm::div(alphaRhoPhi, epsilon_)
      - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_)
     ==
        C1_*alpha()*rho()*G*epsilon_()/k_()
      - fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_)
      - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_)
      + epsilonSource()
      + fvOptions(alpha, rho, epsilon_)
    );

或者复制模仿realizableKE的格式，赋值于epsilonSource() 或者 kSource()．