Commit 2ec0a1df authored by Henry's avatar Henry
Browse files

Corrected calculation of multi-layer resistance

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1447
parent c347c80f
......@@ -117,14 +117,12 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
if (thicknessLayers_.size() > 0)
{
// Calculate effective thermal resistance by harmonic averaging
forAll (thicknessLayers_, iLayer)
{
const scalar l = thicknessLayers_[iLayer];
if (l > 0.0)
{
contactRes_ += kappaLayers_[iLayer]/l;
}
contactRes_ += thicknessLayers_[iLayer]/kappaLayers_[iLayer];
}
contactRes_ = 1.0/contactRes_;
}
}
......
......@@ -123,16 +123,12 @@ turbulentTemperatureRadCoupledMixedFvPatchScalarField
if (thicknessLayers_.size() > 0)
{
// total thermal transmittance by harmonic averaging
// Calculate effective thermal resistance by harmonic averaging
forAll (thicknessLayers_, iLayer)
{
const scalar l = thicknessLayers_[iLayer];
if (l > 0.0)
{
contactRes_ += l/kappaLayers_[iLayer]; // inverse sum
}
contactRes_ += thicknessLayers_[iLayer]/kappaLayers_[iLayer];
}
contactRes_ = 1.0/contactRes_; // new total inverse
contactRes_ = 1.0/contactRes_;
}
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment