[farzadmech]

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

[sharonyue]

This is the explicit term, it is included in phiForce instead and goes to the pressure equation.