twoPhaseEulerFoam & Source term
I'm willing to include a mass source term in twoPhaseEulerFoam solver :
S_Mass = - rho1 * alpha1 / dT, with dT the time step.
It should extract the mass of the phase1 (dispersed phase) in a zone I define with a function :
F = 1 if I activate the source term,
That's for the theory. I'm stuck with the OpenFOAM code now.
Since the equation is in volumic form, I wrote :
fvm::laplacian(Dp,p) == fvc::div(phiHByA) + F*alpha1/(runTime.deltaT().value)
Does it seems correct to you ?
|All times are GMT -4. The time now is 19:04.|