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

reactingEulerFoam: phaseForces: Corrected calculation of accumulated force

Browse Source 
Patch contributed by Institute of Fluid Dynamics,
Helmholtz-Zentrum Dresden - Rossendorf (HZDR)
This commit is contained in:
Will Bainbridge
2019-07-25 08:39:39 +01:00
parent 824b243ee8
commit d0579252a0
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
applications/solvers/multiphase/reactingEulerFoam/functionObjects/phaseForces/phaseForces.C
View File

@ -225,7 +225,7 @@ bool Foam::functionObjects::phaseForces::read(const dictionary& dict)
bool Foam::functionObjects::phaseForces::execute()
{
forAllIter
forAllConstIter
(
HashPtrTable<volVectorField>,
forceFields_,
@ -267,17 +267,17 @@ bool Foam::functionObjects::phaseForces::execute()
if (type == "liftModel")
{
force = nonDragForce<liftModel>(pair);
force += nonDragForce<liftModel>(pair);
}
if (type == "wallLubricationModel")
{
force = nonDragForce<wallLubricationModel>(pair);
force += nonDragForce<wallLubricationModel>(pair);
}
if (type == "turbulentDispersionModel")
{
force = nonDragForce<turbulentDispersionModel>(pair);
force += nonDragForce<turbulentDispersionModel>(pair);
}
}
}
@ -289,7 +289,7 @@ bool Foam::functionObjects::phaseForces::execute()
bool Foam::functionObjects::phaseForces::write()
{
forAllIter
forAllConstIter
(
HashPtrTable<volVectorField>,
forceFields_,