zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Merge branch 'master' of ssh://dm/home/dm4/OpenFOAM/OpenFOAM-dev

Browse Source
This commit is contained in:
Henry
2012-05-04 16:20:14 +01:00
parent e5cd8fff64 060a2b0f42
commit 8b5826bdc7
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/thermophysicalModels/thermalPorousZone/thermalModel/fixedTemperature/fixedTemperature.C
View File

@ -82,20 +82,20 @@ void Foam::porousMedia::fixedTemperature::addEnthalpySource
scalarField& hDiag = hEqn.diag(); scalarField& hDiag = hEqn.diag();
scalarField& hSource = hEqn.source(); scalarField& hSource = hEqn.source();
tmp<volScalarField> Cp = thermo.Cp(); const scalarField T(hDiag.size(), T_);
// TODO: generalize for non-fixedTemperature methods
const scalar rate = 1e6; const scalar rate = 1e6;
forAll(zones, zoneI) forAll(zones, zoneI)
{ {
const labelList& cells = mesh.cellZones()[zones[zoneI]]; const labelList& cells = mesh.cellZones()[zones[zoneI]];
tmp<scalarField> h = thermo.h(T, cells);
forAll(cells, i) forAll(cells, i)
{ {
hDiag[cells[i]] += rate*V[cells[i]]*rho[cells[i]]; hDiag[cells[i]] += rate*V[cells[i]]*rho[cells[i]];
hSource[cells[i]] += hSource[cells[i]] +=
rate*V[cells[i]]*rho[cells[i]]*Cp()[cells[i]]*T_; rate*V[cells[i]]*rho[cells[i]]*h()[cells[i]];
} }
} }
} }