[wangziyang]

hello, everyone!

i am read the InterPhaseChangeFaom code recently,(the version of of is 4.1)

in alphaEqn.H, the MULES algorithm is used to solve the equation, but i don't konw what the meaning of the follow code:

# if (MULESCorr)
{
talphaPhiCorr.ref() -= talphaPhi();

volScalarField alpha100("alpha100", alpha10);
alpha10 = alpha1;

MULES::correct
(
geometricOneField(),
alpha1,
talphaPhi(),
talphaPhiCorr.ref(),
vDotvmcAlphal,
(
divU*(alpha10 - alpha100)
- vDotvmcAlphal*alpha10
)(),
oneField(),
zeroField()
);

// Under-relax the correction for all but the 1st corrector
if (aCorr == 0)
{
talphaPhi.ref() += talphaPhiCorr();
}
else
{
alpha1 = 0.5*alpha1 + 0.5*alpha10;
talphaPhi.ref() += 0.5*talphaPhiCorr();
}
}
else // explict MULES, i have know it
{
MULES::explicitSolve
(
geometricOneField(),
alpha1,
phi,
talphaPhiCorr.ref(),
vDotvmcAlphal,
(divU*alpha1 + vDotcAlphal)(),
oneField(),
zeroField()
);

talphaPhi = talphaPhiCorr;
}

alpha2 = 1.0 - alpha1;
}#

the code seems the semi-MULES algorithm in alphaEqn.H, but i don't know it how works

can some one help me? and can some one recommend some referene?

best wishes!

ziyang