[bojiezhang]

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

[idrama]

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

Cheers

[sabin.ceuca]

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

Good luck

[romant]

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