[gu1]

Hi,

I would like help answering a question.

The NicenoKEqn.C file calculates kSource() as follows:

template<class BasicMomentumTransportModel>
tmp<fvScalarMatrix> NicenoKEqn<BasicMomentumTransportModel>::kSource() const
{
    const alphaField& alpha = this->alpha_;
    const rhoField& rho = this->rho_;

    const phaseCompressibleMomentumTransportModel& gasTurbulence =
        this->gasTurbulence();

    const volScalarField phaseTransferCoeff(this->phaseTransferCoeff());

    return
        alpha*rho*bubbleG()
      + phaseTransferCoeff*gasTurbulence.k()
      - fvm::Sp(phaseTransferCoeff, this->k_);
}

Can I deduce that the last three lines of this code are replacing their respective lines in the KEqn model? Code below:

KEqn.C:

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(this->Ce_*alpha*rho*sqrt(k_)/this->delta(), k_)
      + kSource()
      + fvModels.source(alpha, rho, k_)
    );

Or am I inserting more tree terms into the k equation?

I'm trying to believe that the result of this combination will generate the code below, but I'm not sure... which would even make the code written by Niceno and referenced in OpenFOAM contrary to the article.

Expectation:

tmp<fvScalarMatrix> kEqn
    (
        fvm::ddt(alpha, rho, k_)
      + fvm::div(alphaRhoPhi, k_)
      - fvm::laplacian(alpha*rho*DkEff(), k_)
     ==
        alpha*rho*bubbleG()
      - fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_)
      + phaseTransferCoeff*gasTurbulence.k()
      - fvm::Sp(phaseTransferCoeff, this->k_);
      + fvModels.source(alpha, rho, k_)
    );

Thanks.