From 008d015689b3fb90ee9277e19947da4105a90991 Mon Sep 17 00:00:00 2001
From: Kutalmis Bercin <kutalmis.bercin@esi-group.com>
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 12bd0c2f69d..55aed1f96a4 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<volScalarField> IDefaultPtr;
--
GitLab