zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

removed unnecessary construction of vol fields

Browse Source
This commit is contained in:
andy
2009-03-20 18:57:37 +00:00
parent 763817a3db
commit 0428874b2d
1 changed files with 6 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/thermophysicalModels/radiation/radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
View File

@ -172,26 +172,24 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct()
// reset boundary heat flux to zero // reset boundary heat flux to zero
Qr_ = dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0); Qr_ = dimensionedScalar("zero", dimMass/pow3(dimTime), 0.0);
scalar maxResidual = 0.0; scalar maxResidual = -GREAT;
forAll(IWave_, lambdaI) forAll(IWave_, lambdaI)
{ {
const volScalarField& k = dom_.aj(lambdaI); const volScalarField& k = dom_.aj(lambdaI);
volScalarField E =
absEmmModel_.ECont(lambdaI)/Foam::mathematicalConstant::pi;
surfaceScalarField Ji = dAve_ & mesh_.Sf(); surfaceScalarField Ji = dAve_ & mesh_.Sf();
volScalarField Ib =
blackBody_.bj(lambdaI)/Foam::mathematicalConstant::pi;
fvScalarMatrix IiEq fvScalarMatrix IiEq
( (
fvm::div(Ji, IWave_[lambdaI], " div(Ji,Ii_h)") fvm::div(Ji, IWave_[lambdaI], " div(Ji,Ii_h)")
+ fvm::Sp(k*omega_, IWave_[lambdaI]) + 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(); IiEq.relax();