[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

[MMNCH]

dear ziyang

I suggest you to read this thesis:
An Extended Mixture Model for the Simultaneous Treatment of Short and Long Scale Interfaces.

I hope you find it useful

[wangziyang]

Originally Posted by MMNCH

dear ziyang

I suggest you to read this thesis:
An Extended Mixture Model for the Simultaneous Treatment of Short and Long Scale Interfaces.

I hope you find it useful

thank you very much!

i will read it carefully!

ziyang

[Elol]

Hello Foamers,

So, I have some problems regarding the the functions in Mules as I understand the theory from Santiago thesis but I cant understand the difference betweent the MULES functions. as follow

MULES::limit()
MULES::limitSmm()
MULES::correct()
MULES::explicitSolve()

I went back to source code but I couldn't find out the differnces

what each of above functions do ? and when we can use them?

Best Regards
Mohamed Elerian