research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

class template hePsiThermo modified method correct to update qc_

Browse source
This commit is contained in:
ignis 2015-11-30 13:32:07 +09:00
parent 5c9dbef47b
commit a8d3a17984
1 changed files with 9 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
src/thermophysicalModels/basic/psiThermo/hePsiThermo.C
View file

@ -37,6 +37,7 @@ void Foam::hePsiThermo<BasicPsiThermo, MixtureType>::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<BasicPsiThermo, MixtureType>::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<BasicPsiThermo, MixtureType>::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<BasicPsiThermo, MixtureType>::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<BasicPsiThermo, MixtureType>::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();
}
}
}