[sharonyue]

Hi guys,

In twoLiquidMixingFoam, at first we solve alpha by MULES, then rhoPhi is calculated by the following:

rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;

This makes sense. Then alpha is solved in alphaDissEqn.H, rhoPhi is updated as following:

rhoPhi += alpha1Eqn.flux()*(rho1 - rho2);

If we put them together, finally, rhoPhi will be:

rhoPhi = (alphaPhi + alpha1Eqn.flux())*(rho1 - rho2) + phi*rho2;

I think we need to delete the old alphaPhi in the code above, so in the alphaDissEqn.H, it should be:

rhoPhi = alpha1Eqn.flux()*(rho1 - rho2) + phi*rho2;

Am I right?