CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   MULES problem (

bojiezhang November 7, 2011 23:32

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 =
+ 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!


idrama November 25, 2011 06:10

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


sabin.ceuca November 28, 2011 08:45

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

Good luck ;)

romant November 28, 2011 09:47

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

All times are GMT -4. The time now is 04:44.