[dschmidt]

Dear Foamers,

I'm trying to implement a passive scalar transport into multiphaseEulerFoam (2.2.x).
The scalar is only transported by one of the phases, so I need to access the corresponding phi.

I got a little confused by phases.phiAlpha(), which is used in the UEqns.H and is calculated within the multiphaseSystem
and phases.phi(), which is used for calculating the overall phi (Peqn.H) and is calculated as:

PEqn.H
...
 phase.phi() =
                    phiHbyAs[phasei]
                  + rAlphaAUfs[phasei]*mSfGradp/phase.rho();
...

Both seem to be described in the phaseModel.H, but that didn't help me either

        //- Volumetric flux of the phase
        surfaceScalarField phiAlpha_;

        //- Volumetric flux for the phase
        autoPtr<surfaceScalarField> phiPtr_;

My scalar transport equation with the two variations looks like this:

 fvm::ddt(alpha, Tracer)     
+ fvm::div(fvc::interpolate(alpha)*phase.phi(),Tracer)
 //+ fvm::div(phase.phiAlpha(), Tracer)  
 - fvm::laplacian(fvc::interpolate(alpha)*DTotal, Tracer)

I would be glad if anyone could clarify this for me.


Thanks,
Dominik