在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。
