songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Revert "BUG: Fixing ray dAve and omega for 1D and 2D cases"

Browse Source 
This reverts commit 5848b0afd5.
This commit is contained in:
Kutalmis Bercin 2023-05-05 13:23:41 +01:00
parent 324ce9adec
commit fe97c8fee9
1 changed files with 0 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
View File

@ -140,17 +140,6 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
if (mesh_.nSolutionD() == 2) if (mesh_.nSolutionD() == 2)
{ {
// Omega for 2D
omega_ = deltaPhi;
// dAve for 2D
dAve_ = vector
(
2*sinPhi*Foam::sin(0.5*deltaPhi),
2*cosPhi*Foam::sin(0.5*deltaPhi),
0
);
vector meshDir(Zero); vector meshDir(Zero);
if (dom_.meshOrientation() != vector::zero) if (dom_.meshOrientation() != vector::zero)
{ {
@ -172,7 +161,6 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
dAve_ = coordRot & dAve_; dAve_ = coordRot & dAve_;
d_ = coordRot & d_; d_ = coordRot & d_;
} }
else if (mesh_.nSolutionD() == 1) else if (mesh_.nSolutionD() == 1)
{ {
@ -195,9 +183,6 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay
dAve_ = (dAve_ & normal)*meshDir; dAve_ = (dAve_ & normal)*meshDir;
d_ = (d_ & normal)*meshDir; d_ = (d_ & normal)*meshDir;
// Omega normalization for 1D
omega_ /= 2;
} }
autoPtr<volScalarField> IDefaultPtr; autoPtr<volScalarField> IDefaultPtr;