Skip to content
Snippets Groups Projects
Commit 8bee77c7 authored by Henry Weller's avatar Henry Weller Committed by Andrew Heather
Browse files

BUG: blackBodyEmission: correct non-coupled boundary values

Resolves bug-report https://bugs.openfoam.org/view.php?id=2709
parent 141e219f
Branches
Tags
No related merge requests found
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
......@@ -168,9 +168,9 @@ Foam::radiation::blackBodyEmission::blackBodyEmission
IOobject::NO_WRITE
),
physicoChemical::sigma*pow4(T)
)
);
}
}
......@@ -188,7 +188,7 @@ Foam::scalar Foam::radiation::blackBodyEmission::fLambdaT
const scalar lambdaT
) const
{
return table_(lambdaT*1.0e6);
return table_(1e6*lambdaT);
}
......@@ -226,11 +226,32 @@ Foam::radiation::blackBodyEmission::EbDeltaLambdaT
forAll(T, i)
{
scalar T1 = fLambdaT(band[1]*T[i]);
scalar T2 = fLambdaT(band[0]*T[i]);
const scalar T1 = fLambdaT(band[1]*T[i]);
const scalar T2 = fLambdaT(band[0]*T[i]);
Ebif[i] *= T1 - T2;
}
volScalarField::Boundary& EbBf = Eb.ref().boundaryFieldRef();
forAll(EbBf, patchi)
{
fvPatchScalarField& EbPf = EbBf[patchi];
if (!EbPf.coupled())
{
const scalarField& Tpf = T.boundaryField()[patchi];
forAll(EbPf, facei)
{
const scalar T1 = fLambdaT(band[1]*Tpf[facei]);
const scalar T2 = fLambdaT(band[0]*Tpf[facei]);
EbPf[facei] *= T1 - T2;
}
}
}
return Eb;
}
}
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment