From a8d3a17984876c4ba154cb219d67fbbe3f97dabc Mon Sep 17 00:00:00 2001 From: ignis Date: Mon, 30 Nov 2015 13:32:07 +0900 Subject: [PATCH] class template hePsiThermo modified method correct to update qc_ --- src/thermophysicalModels/basic/psiThermo/hePsiThermo.C | 9 +++++++++ 1 file changed, 9 insertions(+) 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(); } } }