ENH: Updated film radiation model

This commit is contained in:
andy
2012-10-12 15:47:19 +01:00
parent 643f4cb906
commit eab309a838
2 changed files with 9 additions and 9 deletions

View File

@ -57,11 +57,11 @@ standardRadiation::standardRadiation
) )
: :
filmRadiationModel(typeName, owner, dict), filmRadiationModel(typeName, owner, dict),
QrPrimary_ QinPrimary_
( (
IOobject IOobject
( (
"Qr", // same name as Qr on primary region to enable mapping "Qin", // same name as Qin on primary region to enable mapping
owner.time().timeName(), owner.time().timeName(),
owner.regionMesh(), owner.regionMesh(),
IOobject::NO_READ, IOobject::NO_READ,
@ -103,7 +103,7 @@ standardRadiation::~standardRadiation()
void standardRadiation::correct() void standardRadiation::correct()
{ {
// Transfer Qr from primary region // Transfer Qr from primary region
QrPrimary_.correctBoundaryConditions(); QinPrimary_.correctBoundaryConditions();
} }
@ -128,13 +128,13 @@ tmp<volScalarField> standardRadiation::Shs()
); );
scalarField& Shs = tShs(); scalarField& Shs = tShs();
const scalarField& QrP = QrPrimary_.internalField(); const scalarField& QinP = QinPrimary_.internalField();
const scalarField& delta = delta_.internalField(); const scalarField& delta = delta_.internalField();
Shs = beta_*(QrP*pos(delta - deltaMin_))*(1.0 - exp(-kappaBar_*delta)); Shs = beta_*(QinP*pos(delta - deltaMin_))*(1.0 - exp(-kappaBar_*delta));
// Update net Qr on local region // Update net Qr on local region
QrNet_.internalField() = QrP - Shs; QrNet_.internalField() = QinP - Shs;
QrNet_.correctBoundaryConditions(); QrNet_.correctBoundaryConditions();
return tShs; return tShs;

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -59,8 +59,8 @@ private:
// Private data // Private data
//- Radiative flux mapped from the primary region / [kg/s3] //- Radiative incident flux mapped from the primary region / [kg/s3]
volScalarField QrPrimary_; volScalarField QinPrimary_;
//- Remaining radiative flux after removing local contribution //- Remaining radiative flux after removing local contribution
volScalarField QrNet_; volScalarField QrNet_;