zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

multiphaseEulerFoam: sizeDistribution: Corrected spatial averaging of concentration fields

Browse Source 
Patch contributed by Institute of Fluid Dynamics,
Helmholtz-Zentrum Dresden - Rossendorf (HZDR)
This commit is contained in:
Will Bainbridge
2022-02-16 10:22:04 +00:00
parent fb4e237292
commit d0c63582a2
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
applications/solvers/multiphase/multiphaseEulerFoam/functionObjects/sizeDistribution/sizeDistribution.C
View File

@ -452,7 +452,7 @@ bool Foam::functionObjects::sizeDistribution::write()
const diameterModels::sizeGroup& fi = sizeGroups[i]; const diameterModels::sizeGroup& fi = sizeGroups[i];
resultValues[i] = resultValues[i] =
gSum(filterField(mesh_.V()*fi*fi.phase()/fi.x()))*this->V(); gSum(filterField(mesh_.V()*fi*fi.phase()/fi.x()))/this->V();
} }
if (normalise_ && sum(resultValues) != 0) if (normalise_ && sum(resultValues) != 0)
@ -469,7 +469,7 @@ bool Foam::functionObjects::sizeDistribution::write()
const diameterModels::sizeGroup& fi = sizeGroups[i]; const diameterModels::sizeGroup& fi = sizeGroups[i];
resultValues[i] = resultValues[i] =
gSum(filterField(mesh_.V()*fi*fi.phase()/fi.x()))*this->V(); gSum(filterField(mesh_.V()*fi*fi.phase()/fi.x()))/this->V();
} }
if (normalise_ && sum(resultValues) != 0) if (normalise_ && sum(resultValues) != 0)
@ -491,7 +491,7 @@ bool Foam::functionObjects::sizeDistribution::write()
const diameterModels::sizeGroup& fi = sizeGroups[i]; const diameterModels::sizeGroup& fi = sizeGroups[i];
resultValues[i] = resultValues[i] =
gSum(filterField(mesh_.V()*fi*fi.phase()))*this->V(); gSum(filterField(mesh_.V()*fi*fi.phase()))/this->V();
} }
if (normalise_ && sum(resultValues) != 0) if (normalise_ && sum(resultValues) != 0)
@ -508,7 +508,7 @@ bool Foam::functionObjects::sizeDistribution::write()
const diameterModels::sizeGroup& fi = sizeGroups[i]; const diameterModels::sizeGroup& fi = sizeGroups[i];
resultValues[i] = resultValues[i] =
gSum(filterField(mesh_.V()*fi*fi.phase()))*this->V(); gSum(filterField(mesh_.V()*fi*fi.phase()))/this->V();
} }
if (normalise_ && sum(resultValues) != 0) if (normalise_ && sum(resultValues) != 0)
@ -534,7 +534,7 @@ bool Foam::functionObjects::sizeDistribution::write()
( (
filterField(mesh_.V()*fi.a().ref()*fi*fi.phase()/fi.x()) filterField(mesh_.V()*fi.a().ref()*fi*fi.phase()/fi.x())
) )
*this->V(); /this->V();
} }
if (normalise_ && sum(resultValues) != 0) if (normalise_ && sum(resultValues) != 0)
@ -555,7 +555,7 @@ bool Foam::functionObjects::sizeDistribution::write()
( (
filterField(mesh_.V()*fi.a().ref()*fi*fi.phase()/fi.x()) filterField(mesh_.V()*fi.a().ref()*fi*fi.phase()/fi.x())
) )
*this->V(); /this->V();
} }
if (normalise_ && sum(resultValues) != 0) if (normalise_ && sum(resultValues) != 0)