From a4d76986c7ae01ff10b7fea5fe349a61b7e2a594 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Tue, 6 Oct 2009 14:44:05 +0100
Subject: [PATCH] radiation fixes
---
.../fvDOM/absorptionCoeffs/absorptionCoeffs.C | 2 +-
.../radiation/radiationModel/fvDOM/fvDOM/fvDOM.C | 5 ++---
.../radiativeIntensityRay/radiativeIntensityRay.C | 10 ++--------
.../greyMeanAbsorptionEmission.C | 7 ++++---
4 files changed, 9 insertions(+), 15 deletions(-)
diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/absorptionCoeffs/absorptionCoeffs.C b/src/thermophysicalModels/radiation/radiationModel/fvDOM/absorptionCoeffs/absorptionCoeffs.C
index eeec89cef8f..041d0f7b007 100644
--- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/absorptionCoeffs/absorptionCoeffs.C
+++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/absorptionCoeffs/absorptionCoeffs.C
@@ -99,7 +99,7 @@ void Foam::radiation::absorptionCoeffs::initialise(const dictionary& dict)
{
dict.lookup("Tcommon") >> Tcommon_;
dict.lookup("Tlow") >> Tlow_;
- dict.lookup("Tlow") >> Thigh_;
+ dict.lookup("Thigh") >> Thigh_;
dict.lookup("invTemp") >> invTemp_;
dict.lookup("loTcoeffs") >> lowACoeffs_;
diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C b/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
index 81f8457ba46..aaf5b544966 100644
--- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
+++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
@@ -26,7 +26,6 @@ License
#include "fvDOM.H"
#include "addToRunTimeSelectionTable.H"
-#include "fvm.H"
#include "absorptionEmissionModel.H"
#include "scatterModel.H"
@@ -65,7 +64,7 @@ Foam::radiation::fvDOM::fvDOM(const volScalarField& T)
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
- IOobject::NO_WRITE
+ IOobject::AUTO_WRITE
),
mesh_,
dimensionedScalar("G", dimMass/pow3(dimTime), 0.0)
@@ -347,7 +346,7 @@ Foam::radiation::fvDOM::Ru() const
const DimensionedField<scalar, volMesh> a =
a_.dimensionedInternalField(); //absorptionEmission_->aCont()()
- return a*G - 4.0*E;
+ return a*G - E;
}
diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C b/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
index cb94ca5f9da..79029be3b80 100644
--- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
+++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
@@ -28,12 +28,6 @@ License
#include "fvm.H"
#include "fvDOM.H"
-#include "absorptionEmissionModel.H"
-#include "scatterModel.H"
-#include "mathConstants.H"
-#include "radiationModel.H"
-#include "Vector2D.H"
-
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
Foam::label Foam::radiation::radiativeIntensityRay::rayId(0);
@@ -202,9 +196,9 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct()
fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)")
+ fvm::Sp(k*omega_, ILambda_[lambdaI])
==
- 1.0/constant::math::pi
+ 1.0/Foam::mathematicalConstant::pi*omega_
*(
- k*omega_*blackBody_.bLambda(lambdaI)
+ k*blackBody_.bLambda(lambdaI)
+ absorptionEmission_.ECont(lambdaI)
)
);
diff --git a/src/thermophysicalModels/radiation/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C b/src/thermophysicalModels/radiation/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C
index cdb57c0f539..6204188f1aa 100644
--- a/src/thermophysicalModels/radiation/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C
+++ b/src/thermophysicalModels/radiation/submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C
@@ -259,10 +259,11 @@ Foam::radiation::greyMeanAbsorptionEmission::ECont(const label bandI) const
)
);
- if (mesh_.foundObject<volScalarField>("hrr"))
+ if (mesh_.foundObject<volScalarField>("dQ"))
{
- const volScalarField& hrr = mesh_.lookupObject<volScalarField>("hrr");
- E().internalField() = EhrrCoeff_*hrr.internalField();
+ const volScalarField& dQ =
+ mesh_.lookupObject<volScalarField>("dQ");
+ E().internalField() = EhrrCoeff_*dQ;
}
return E;
--
GitLab