Skip to content
Snippets Groups Projects
Commit 86af21c8 authored by mattijs's avatar mattijs
Browse files

radiation fixes

parent e5b647ed
Branches
Tags
No related merge requests found
...@@ -27,6 +27,9 @@ License ...@@ -27,6 +27,9 @@ License
#include "radiativeIntensityRay.H" #include "radiativeIntensityRay.H"
#include "fvm.H" #include "fvm.H"
#include "fvDOM.H" #include "fvDOM.H"
#include "constants.H"
using namespace Foam::constant;
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
...@@ -196,7 +199,7 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct() ...@@ -196,7 +199,7 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct()
fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)") fvm::div(Ji, ILambda_[lambdaI], "div(Ji,Ii_h)")
+ fvm::Sp(k*omega_, ILambda_[lambdaI]) + fvm::Sp(k*omega_, ILambda_[lambdaI])
== ==
1.0/Foam::mathematicalConstant::pi*omega_ 1.0/constant::math::pi*omega_
*( *(
k*blackBody_.bLambda(lambdaI) k*blackBody_.bLambda(lambdaI)
+ absorptionEmission_.ECont(lambdaI) + absorptionEmission_.ECont(lambdaI)
......
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