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:
)/(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!
Have you looked up in alphaEqn.H ?? As far as I know, there are the equation are stored.
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.
|All times are GMT -4. The time now is 10:34.|