|
[Sponsors] |
![]() |
![]() |
#1 |
Senior Member
Farzad Faraji
Join Date: Nov 2019
Posts: 206
Rep Power: 7 ![]() |
Dear Foamers
I have a question about MPPICFoam. Why before pimple loop we put cloud.source = Zero?? cccccccccccccccccccccccccccccccccccccccccccccccccc cccccccccccccccccccc // Update continuous phase volume fraction field alphac = max(1.0 - kinematicCloud.theta(), alphacMin); alphac.correctBoundaryConditions(); alphacf = fvc::interpolate(alphac); alphaPhic = alphacf*phic; fvVectorMatrix cloudSU(kinematicCloud.SU(Uc)); volVectorField cloudVolSUSu ( IOobject ( "cloudVolSUSu", runTime.timeName(), mesh ), mesh, dimensionedVector ( "0", cloudSU.dimensions()/dimVolume, Zero ), zeroGradientFvPatchVectorField::typeName ); cloudVolSUSu.primitiveFieldRef() = -cloudSU.source()/mesh.V(); cloudVolSUSu.correctBoundaryConditions(); cloudSU.source() = Zero; // --- Pressure-velocity PIMPLE corrector loop while (pimple.loop()) { #include "UcEqn.H" // --- PISO loop while (pimple.correct()) { #include "pEqn.H" } if (pimple.turbCorr()) { continuousPhaseTurbulence->correct(); } } cccccccccccccccccccccccccccccccccccccccccccccccccc cccccccccccccccccccc |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 844
Rep Power: 18 ![]() |
This is the explicit term, it is included in phiForce instead and goes to the pressure equation.
__________________
My OpenFOAM algorithm website: http://dyfluid.com By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html |
|
![]() |
![]() |
![]() |
Tags |
cloudsu.source, mppicfoam, pimpleloop |
|
|