CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   MULES problem (http://www.cfd-online.com/Forums/openfoam/94163-mules-problem.html)

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 =
(
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 November 25, 2011 06:10

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

Cheers

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 10:16.