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

DPMFoam: Changed the cloud source splitting to handle symmetric semi-implicit sources consistently

Browse Source 
Resolves bug-report https://bugs.openfoam.org/view.php?id=3385
This commit is contained in:
Henry Weller
2019-11-13 17:15:12 +00:00
parent 80bb6e2b0f
commit b22545b931
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
applications/solvers/lagrangian/DPMFoam/DPMFoam.C
View File

@ -133,9 +133,10 @@ int main(int argc, char *argv[])
zeroGradientFvPatchVectorField::typeName
);
cloudVolSUSu.primitiveFieldRef() = -cloudSU.source()/mesh.V();
cloudVolSUSu.primitiveFieldRef() =
(cloudSU.diag()*Uc() - cloudSU.source())/mesh.V();
cloudVolSUSu.correctBoundaryConditions();
cloudSU.source() = Zero;
cloudSU.source() = cloudSU.diag()*Uc();
// --- Pressure-velocity PIMPLE corrector loop
while (pimple.loop())