Skip to content
Snippets Groups Projects
Commit e0c1ac3b authored by Andrew Heather's avatar Andrew Heather
Browse files

cosmetics

parent 98363e58
Branches
Tags
No related merge requests found
......@@ -37,19 +37,19 @@ Description
i.e. dictionary
fvDOMCoeffs
{
Nphi 1; // azimuthal angles in PI/2 on X-Y.(from Y to X)
Ntheta 2; // polar angles in P1 (from Z to X-Y plane)
nPhi 1; // azimuthal angles in PI/2 on X-Y.(from Y to X)
nTheta 2; // polar angles in PI (from Z to X-Y plane)
convergence 1e-4; // convergence criteria for radiation iteration
}
nFlowIterPerRadIter 1; // Number of flow iterations per radiation
solverFreq 1; // Number of flow iterations per radiation
iteration
The total number of solid angles is 4*Nphi*Ntheta.
The total number of solid angles is 4*nPhi*nTheta.
In 1D the direction of the rays is X (Nphi and Ntheta are ignored)
In 2D the direction of the rays is on X-Y plane (only Nphi is considered)
In 3D (Nphi and Ntheta are considered)
In 1D the direction of the rays is X (nPhi and nTheta are ignored)
In 2D the direction of the rays is on X-Y plane (only nPhi is considered)
In 3D (nPhi and nTheta are considered)
SourceFiles
fvDOM.C
......@@ -127,9 +127,6 @@ class fvDOM
//- Disallow default bitwise assignment
void operator=(const fvDOM&);
//- Update Absorption Coefficients
// void updateAbsorptionCoeffs(void);
//- Update nlack body emission
void updateBlackBodyEmission();
......
......@@ -113,11 +113,11 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
0.5*deltaPhi*Foam::sin(2.0*theta)*Foam::sin(deltaTheta)
);
forAll(ILambda_, i)
forAll(ILambda_, lambdaI)
{
IOobject IHeader
(
"ILambda_" + name(rayId) + "_" + name(i),
"ILambda_" + name(rayId) + "_" + name(lambdaI),
mesh_.time().timeName(),
mesh_,
IOobject::MUST_READ,
......@@ -129,7 +129,7 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
{
ILambda_.set
(
i,
lambdaI,
new volScalarField(IHeader, mesh_)
);
}
......@@ -150,7 +150,7 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
ILambda_.set
(
i,
lambdaI,
new volScalarField(IHeader, IDefault)
);
}
......@@ -201,7 +201,6 @@ Foam::scalar Foam::radiation::radiativeIntensityRay::correct()
).initialResidual();
maxResidual = max(eqnResidual, maxResidual);
}
return maxResidual;
......
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