research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

turbulentHeatFluxTemperature: replaced Cp0 with rhoCp0 and added comment concerning units to clarify usage

Browse source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1433
This commit is contained in:
Henry 2015-01-18 22:01:55 +00:00
parent f7a485069c
commit 7e448b3295
2 changed files with 8 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/turbulenceModels/incompressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.C
View file

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -194,19 +194,19 @@ void turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
// retrieve (constant) specific heat capacity from transport dictionary
const IOdictionary& transportProperties =
db().lookupObject<IOdictionary>("transportProperties");
const scalar Cp0(readScalar(transportProperties.lookup("Cp0")));
const scalar rhoCp0(readScalar(transportProperties.lookup("rhoCp0")));
switch (heatSource_)
{
case hsPower:
{
const scalar Ap = gSum(patch().magSf());
gradient() = q_/(Ap*Cp0*alphaEffp);
gradient() = q_/(Ap*rhoCp0*alphaEffp);
break;
}
case hsFlux:
{
gradient() = q_/(Cp0*alphaEffp);
gradient() = q_/(rhoCp0*alphaEffp);
break;
}
default:

8
src/turbulenceModels/incompressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H
View file

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -54,9 +54,9 @@ Description
\endverbatim
Note
- it is assumed that the units of \c alphaEff are [kg/m/s]
- the specific heat capcaity is read from the transport dictionary entry
\c Cp0
- The units of \c alphaEff are [kg/m/s]
- The density*specific heat capacity \c rhoCp0 (units [kg/m3]*[J/kg/K])
is read from the transport dictionary entry
SourceFiles
turbulentHeatFluxTemperatureFvPatchScalarField.C