View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0001108 | OpenFOAM | Bug | public | 2013-12-12 10:46 | 2013-12-12 12:31 |
Reporter | Assigned To | ||||
Priority | normal | Severity | minor | Reproducibility | N/A |
Status | resolved | Resolution | fixed | ||
Platform | Linux | OS | Debian | OS Version | (please specify) |
Summary | 0001108: Bug in externalWallHeatFluxTemperatureFvPatchScalarField | ||||
Description | I Believe something may have gone wrong in the updated implementation of externalWallHeatFluxTemperatureFvPatchScalarField (line 256-268 of the .C file). The heat flux is now calculated as: scalar totalSolidRes = 0.0; if (thicknessLayers_.size() > 0) { forAll (thicknessLayers_, iLayer) { const scalar l = thicknessLayers_[iLayer]; if (l > 0.0) { totalSolidRes += kappaLayers_[iLayer]/l; } } } q = (Ta_ - Tc)/(1.0/h_ + totalSolidRes); I believe the resistance should be calculated as kappa/l, rather than its inverse? In the previous version of the code (only single layer), this effectively has been the case, i.e. q = (Ta_ - Tc)/(1.0/h_ + thicknessLayer_/(kappaLayer_ + VSMALL)); | ||||
Tags | No tags attached. | ||||
|
Solved in commint 7d4b14579af9fb5dd8dd04dc9839432ab044e81a http://openfoam.org/mantisbt/view.php?id=1108 Thanks |