怎么在k-e湍流模型中加源项

[史浩]

计算中，需要在k-e模型两个方程中添加源项，形式为-Kk和 -Ke，但是目前来看k-e模型是封装好的，好像改不动。这个有修改的方法么？

[东岳]

https://github.com/OpenFOAM/OpenFOAM-3.0.x/blob/master/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C

去这里加个源项不就好了？

// Turbulent kinetic energy equation
    tmp<fvScalarMatrix> kEqn
    (
        fvm::ddt(alpha, rho, k_)
      + fvm::div(alphaRhoPhi, k_)
      - fvm::laplacian(alpha*rho*DkEff(), k_)
     ==
        alpha*rho*G
      - fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_)
      - fvm::Sp(alpha*rho*epsilon_/k_, k_)
      + kSource()
      // ?
      + K*k_
    );