[pbachant]

I'm trying to create a turbulence source to accompany a momentum source (think a bluff body or stalled airfoil). I have some reference values for k and epsilon from a simulation with a body-fitted grid that I am trying to match, but I am having trouble figuring out how to translate those values into source values (or constraints or corrections). Here is the k equation syntax:

    // 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()
      + fvOptions(alpha, rho, k_)
    );

    kEqn.ref().relax();
    fvOptions.constrain(kEqn.ref());
    solve(kEqn);
    fvOptions.correct(k_);
    bound(k_, this->kMin_);

It seems intuitively that `correct` is the method to use to set the values of the field directly, but since the equation has already been solved, does this mean the field won't necessarily be a solution of the equation? Does this matter? When I try to set values using the `correct` method, they seem to end up 100 times larger than they should be.

Is the`addSup` method more correct to use if this is a turbulence source, and if so, how would you go about estimating a rate from a target value?

[Tibo99]

Did you find a solution about your problem?

I got similar issue.

Regards