interCondensatingEvaporatingFoam中的rhoCpPhi

[BlookCFD]

在v2012版本中，TEqn.H中有rhoCpPhi的定义为：

     const surfaceScalarField rhoCpPhi
    (
         "rhoCpPhi",
         rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2
     );

这个定义是如何得到的? 如果Cp1 = Cp2，那rhoCpPhi = phi*rho2*Cp2，这个和 phi*rho*Cp有差距。

如果直接把中心点的rhoCp插值到面上和phi相乘应该表示为 fvc::interpolate(rho1*alpha1*Cp1 + rho2*alpha2*Cp2)*phi。