externalWallHeatFluxTemperature in turbulence models
I wish to apply a nonuniform heat-transfer-coefficient boundary coefficient to a wall. I have used the following boundary condition for temperature:
h nonuniform List<scalar>
If I specify 'h' with units [W/m2/K] I get far too much heat leaving the boundary.
I'm having trouble following the units through externalWallHeatFluxTemperatureFvPatchScalarField. C
00210 q = (Ta_ - *this)*h_;
I think this implies [q] = [K - K]*[W/m2/K] = [W/m2].
00224 if (q[i] > 0) //in
00226 this->refGrad()[i] = q[i]/K(*this)()[i];
00227 this->refValue()[i] = 0.0;
00228 this->valueFraction()[i] = 0.0;
00230 else //out
00232 this->refGrad()[i] = 0.0;
00233 this->refValue()[i] = KDelta[i]*q[i] + patchInternalField()()[i];00234 this->valueFraction()[i] = 1.0;
For the heat input case, I see [refGrad]=[W/m2] / [W/m/K] = [K/m], which is good.
For the heat output case, I think the units of kDelta are [W/m/K]*[1/m]=[W/m2/K] based on mesh.surfaceInterpolation.c and
00202 scalarField KDelta(K(*this)*patch().deltaCoeffs());
Doesn't this imply the units of refValue are [refValue] = [W/m2/K] * [W/m/K] + [K] ?
Is this a bug or have I mis-understood something?
first i want to say that i am interested in getting the problem solved from the first post "If I specify 'h' with units [W/m2/K] I get far too much heat leaving the boundary" too. Can anyone explain the question from the first post?
My question deals with the problem how to read out data for Ta from a .dat file; is it in general possible?
K solidThermo; // solidThermo or lookup
q uniform 1000; // Heat flux / [W/m2]
Ta uniform 300.0; // Tambient temperature /[K]
h uniform 10.0; // Heat transfer coeff /[W/Km2]
value uniform 300.0; // Initial temperature / [K]
Can i combine this type with:
Thank you very much!
|All times are GMT -4. The time now is 22:51.|