diff --git a/src/thermophysicalModels/basic/psiThermo/hePsiThermo.C b/src/thermophysicalModels/basic/psiThermo/hePsiThermo.C index 970fc82f..0c60ac25 100644 --- a/src/thermophysicalModels/basic/psiThermo/hePsiThermo.C +++ b/src/thermophysicalModels/basic/psiThermo/hePsiThermo.C @@ -37,6 +37,7 @@ void Foam::hePsiThermo::calculate() scalarField& psiCells = this->psi_.internalField(); scalarField& muCells = this->mu_.internalField(); scalarField& alphaCells = this->alpha_.internalField(); + scalarField& qcCells = this->qc_.internalField(); forAll(TCells, celli) { @@ -54,6 +55,8 @@ void Foam::hePsiThermo::calculate() muCells[celli] = mixture_.mu(pCells[celli], TCells[celli]); alphaCells[celli] = mixture_.alphah(pCells[celli], TCells[celli]); + + qcCells[celli] = mixture_.QC(); } forAll(this->T_.boundaryField(), patchi) @@ -67,6 +70,8 @@ void Foam::hePsiThermo::calculate() fvPatchScalarField& pmu = this->mu_.boundaryField()[patchi]; fvPatchScalarField& palpha = this->alpha_.boundaryField()[patchi]; + fvPatchScalarField& pqc = this->qc_.boundaryField()[patchi]; + if (pT.fixesValue()) { forAll(pT, facei) @@ -79,6 +84,8 @@ void Foam::hePsiThermo::calculate() ppsi[facei] = mixture_.psi(pp[facei], pT[facei]); pmu[facei] = mixture_.mu(pp[facei], pT[facei]); palpha[facei] = mixture_.alphah(pp[facei], pT[facei]); + + pqc[facei] = mixture_.QC(); } } else @@ -93,6 +100,8 @@ void Foam::hePsiThermo::calculate() ppsi[facei] = mixture_.psi(pp[facei], pT[facei]); pmu[facei] = mixture_.mu(pp[facei], pT[facei]); palpha[facei] = mixture_.alphah(pp[facei], pT[facei]); + + pqc[facei] = mixture_.QC(); } } }