Skip to content
Snippets Groups Projects
Commit abfcafc8 authored by andy's avatar andy
Browse files

ENH: Use alpha in film standard radiation model instead of re-calculating wet/dry

parent ba717a55
No related branches found
No related tags found
No related merge requests found
......@@ -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;
......
......@@ -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_;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment