[tehache]

Hello there,

I am having problems trying to understand mixedEnergyFvPatchScalarField::updateCoeffs(), code below.
My problem is this part:

refValue() = thermo.he(pw, Tw.refValue(), patchi);

My feeling is that the refValue for the enthaly should be computed using the thermo from the neighbouring region, at least in the case were Tw is e.g. a turbulentTemperatureCoupledBaffleMixedFvPatchScala rField. Could someone enlighten me, perhaps?


================================================== ========

void Foam::mixedEnergyFvPatchScalarField::updateCoeffs( )
{
if (updated())
{
return;
}

const basicThermo& thermo = basicThermo::lookupThermo(*this);
const label patchi = patch().index();

const scalarField& pw = thermo.p().boundaryField()[patchi];
mixedFvPatchScalarField& Tw = refCast<mixedFvPatchScalarField>
(
const_cast<fvPatchScalarField&>(thermo.T().boundar yField()[patchi])
);

Tw.evaluate();

valueFraction() = Tw.valueFraction();
refValue() = thermo.he(pw, Tw.refValue(), patchi);
refGrad() =
thermo.Cpv(pw, Tw, patchi)*Tw.refGrad()
+ patch().deltaCoeffs()*
(
thermo.he(pw, Tw, patchi)
- thermo.he(pw, Tw, patch().faceCells())
);

mixedFvPatchScalarField::updateCoeffs();
}