Qc and Qc2 field calculation dimension correctly set

This commit is contained in:
ignis 2018-04-30 05:09:10 +09:00
parent 8670f39b9f
commit 6b3e8c972b

View file

@ -96,15 +96,19 @@ Foam::tmp<Foam::volScalarField> Foam::basicSpecieMixture::Qc() const
Y[0].mesh()
),
Y[0].mesh(),
dimensionedScalar("zero", dimless, 0)
dimensionedScalar("zero", dimCurrent*dimTime/dimMass, 0)
)
);
volScalarField& Qc = tQc.ref();
dimensionedScalar Qci("Qci", dimCurrent*dimTime/dimMass, 0);
forAll(Y, i)
{
Qc += Y[i]*this->Qc(i);
Qci.value() = this->Qc(i);
Qc += Y[i]*Qci;
// Qc += Y[i]*this->Qc(i);
}
return tQc;
@ -126,15 +130,19 @@ Foam::tmp<Foam::volScalarField> Foam::basicSpecieMixture::Qc2() const
Y[0].mesh()
),
Y[0].mesh(),
dimensionedScalar("zero", dimless, 0)
dimensionedScalar("zero", sqr(dimCurrent*dimTime)/dimMass, 0)
)
);
volScalarField& Qc = tQc.ref();
dimensionedScalar Qc2i("Qc2i", sqr(dimCurrent*dimTime)/dimMass, 0);
forAll(Y, i)
{
Qc += Y[i]*this->Qc2(i);
Qc2i.value() = this->Qc2(i);
Qc += Y[i]*Qc2i;
//Qc += Y[i]*this->Qc2(i);
}
return tQc;