[shanyeyun]

Hi, Foamers

In the class chtRcTemperatureFvPatchScalarField'' there is a member function:


Source(), defined as follows.


{ const fvPatch& p = patch(); const magLongDelta& mld = magLongDelta::New(p.boundaryMesh().mesh());
const scalarField TcOwn = patchInternalField(); const scalarField TcNei = patchNeighbourField(); const scalarField& Tw = *this;
const chtRegionCoupleBase& K = dynamic_cast<const chtRegionCoupleBase&> ( p.lookupPatchField<volScalarField, scalar>(kName_) );
const scalarField k = K*p.deltaCoeffs();
const scalarField kOwn = K.originalPatchField()/(1 - p.weights())/mld.magDelta(p.index());
return kOwn*(Tw - TcOwn) - k*(TcNei - TcOwn); }


I am confused about the last line of code. As far as I understand, kOwn = thermalConductivity at cellCenter divided by the distance from the inrerface to faceCellCenter.
k = thermalConductivity at the interface divided by the distance from Owner to Neighbour cell center. Then I should return only - k*(TcNei - TcOwn);


Can anybody help me? Thanks a lot.