MULES problem
 

Hello every foamer:
I have a problem about MULES, and I want to add source to my VOF equation, just like below:
MULES::explicitSolve(geometricOneField(), alpha1, phi, phiAlpha, source , geometricZeroField() , 1, 0);
but the result of the surface is inconsistent, so I look deep into the MULES, and find the quation below:
psiIf =
(
rho.oldTime().field()*psi0/deltaT
+ Su.field()
- psiIf
)/(rho.field()/deltaT - Sp.field());

I really can not find any relationship with the VOF equation in Rusche's and Jasak's theises, can anyone explain it? Thank you in advance!

bojiezhang

Have you looked up in alphaEqn.H ?? As far as I know, there are the equation are stored.

Cheers

Hi bojiezhang,
just take a look in the file MULESTemplates.C starting with the line:
00218 fvScalarMatrix psiConvectionDiffusion

Good luck ;)

The interPhaseChangeFoam solver could give you a good example of how a source can be implemented.