From 0428874b2dc015046bfb389925803465f9b4f6b7 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Fri, 20 Mar 2009 18:57:37 +0000
Subject: [PATCH] removed unnecessary construction of vol fields
---
.../radiativeIntensityRay/radiativeIntensityRay.C | 14 ++++++--------
1 file changed, 6 insertions(+), 8 deletions(-)
diff --git a/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C b/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
index b90835db567..23db640feaf 100644
--- a/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
+++ b/src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
@@ -172,26 +172,24 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct()
// reset boundary heat flux to zero
Qr_ = dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0);
- scalar maxResidual = 0.0;
+ scalar maxResidual = -GREAT;
forAll(IWave_, lambdaI)
{
const volScalarField& k = dom_.aj(lambdaI);
- volScalarField E =
- absEmmModel_.ECont(lambdaI)/Foam::mathematicalConstant::pi;
-
surfaceScalarField Ji = dAve_ & mesh_.Sf();
- volScalarField Ib =
- blackBody_.bj(lambdaI)/Foam::mathematicalConstant::pi;
-
fvScalarMatrix IiEq
(
fvm::div(Ji, IWave_[lambdaI], " div(Ji,Ii_h)")
+ fvm::Sp(k*omega_, IWave_[lambdaI])
==
- k*omega_*Ib + E
+ 1.0/Foam::mathematicalConstant::pi
+ *(
+ k*omega_*blackBody_.bj(lambdaI)
+ + absEmmModel_.ECont(lambdaI)
+ )
);
IiEq.relax();
--
GitLab