From fe97c8fee9792c261f75cae21e82e4f1170acd50 Mon Sep 17 00:00:00 2001 From: Kutalmis Bercin Date: Fri, 5 May 2023 13:23:41 +0100 Subject: [PATCH] Revert "BUG: Fixing ray dAve and omega for 1D and 2D cases" This reverts commit 5848b0afd5c931a3e85d673638bae6dab29ce62d. --- .../radiativeIntensityRay/radiativeIntensityRay.C | 15 --------------- 1 file changed, 15 deletions(-) diff --git a/src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C b/src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C index cc1a15a1e8..4fea2866a9 100644 --- a/src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C +++ b/src/thermophysicalModels/radiation/radiationModels/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C @@ -140,17 +140,6 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay 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); if (dom_.meshOrientation() != vector::zero) { @@ -172,7 +161,6 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay dAve_ = coordRot & dAve_; d_ = coordRot & d_; - } else if (mesh_.nSolutionD() == 1) { @@ -195,9 +183,6 @@ Foam::radiation::radiativeIntensityRay::radiativeIntensityRay dAve_ = (dAve_ & normal)*meshDir; d_ = (d_ & normal)*meshDir; - - // Omega normalization for 1D - omega_ /= 2; } autoPtr IDefaultPtr;