From d94744e9f767d684c809566784eb28c71260c489 Mon Sep 17 00:00:00 2001
From: Andrew Heather <>
Date: Tue, 30 May 2023 12:07:33 +0100
Subject: [PATCH] BUG: parcel injection models - corrected volume calculations.
Fixes #2708
---
.../InjectedParticleDistributionInjection.C | 2 +-
.../InjectedParticleInjection/InjectedParticleInjection.C | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleDistributionInjection/InjectedParticleDistributionInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleDistributionInjection/InjectedParticleDistributionInjection.C
index 3cdad615aba..8bc762c114a 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleDistributionInjection/InjectedParticleDistributionInjection.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleDistributionInjection/InjectedParticleDistributionInjection.C
@@ -179,7 +179,7 @@ void Foam::InjectedParticleDistributionInjection<CloudType>::initialise()
sumPow3 += pow3(diameters[particlei]);
}
- const scalar volume = sumPow3*mathematical::pi/16.0;
+ const scalar volume = sumPow3*mathematical::pi/6.0;
sumVolume += volume;
volumeFlowRate_[injectori] = volume/dTime;
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleInjection/InjectedParticleInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleInjection/InjectedParticleInjection.C
index 76e05f08fa2..cd3f5ae3447 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleInjection/InjectedParticleInjection.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectedParticleInjection/InjectedParticleInjection.C
@@ -119,7 +119,7 @@ void Foam::InjectedParticleInjection<CloudType>::initialise()
scalar sumVolume = 0;
forAll(volume, i)
{
- scalar vol = pow3(diameter_[i])*mathematical::pi/16.0;
+ scalar vol = pow3(diameter_[i])*mathematical::pi/6.0;
volume[i] = vol;
sumVolume += vol;
}
--
GitLab