From abfcafc8470a1d9c5d5e06b72490b8f15a4b0736 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 12 Nov 2012 12:03:03 +0000
Subject: [PATCH] ENH: Use alpha in film standard radiation model instead of
re-calculating wet/dry
---
.../standardRadiation/standardRadiation.C | 7 +++----
.../standardRadiation/standardRadiation.H | 6 ------
2 files changed, 3 insertions(+), 10 deletions(-)
diff --git a/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.C b/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.C
index be5460757fe..1ff32747161 100644
--- a/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.C
+++ b/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.C
@@ -85,8 +85,6 @@ standardRadiation::standardRadiation
dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0),
zeroGradientFvPatchScalarField::typeName
),
- delta_(owner.delta()),
- deltaMin_(readScalar(coeffs_.lookup("deltaMin"))),
beta_(readScalar(coeffs_.lookup("beta"))),
kappaBar_(readScalar(coeffs_.lookup("kappaBar")))
{}
@@ -129,9 +127,10 @@ tmp<volScalarField> standardRadiation::Shs()
scalarField& Shs = tShs();
const scalarField& QinP = QinPrimary_.internalField();
- const scalarField& delta = delta_.internalField();
+ const scalarField& delta = owner_.delta().internalField();
+ const scalarField& alpha = owner_.alpha().internalField();
- Shs = beta_*(QinP*pos(delta - deltaMin_))*(1.0 - exp(-kappaBar_*delta));
+ Shs = beta_*QinP*alpha*(1.0 - exp(-kappaBar_*delta));
// Update net Qr on local region
QrNet_.internalField() = QinP - Shs;
diff --git a/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.H b/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.H
index 9b96fb60fe1..799bb6c6493 100644
--- a/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.H
+++ b/src/regionModels/surfaceFilmModels/submodels/thermo/filmRadiationModel/standardRadiation/standardRadiation.H
@@ -65,15 +65,9 @@ private:
//- Remaining radiative flux after removing local contribution
volScalarField QrNet_;
- //- Reference to the film thickness field / [m]
- const volScalarField& delta_;
-
// Model coefficients
- //- Minimum thickness to apply radiation model
- scalar deltaMin_;
-
//- Beta coefficient
scalar beta_;
--
GitLab