\(\left(\rho c_{p}\right)_{\text{hnf}}=\left(1-\phi_{2}\right)\left[\left(1-\phi_{1}\right){(\rho c_{p})}_{\text{bf}}+\phi_{1}{(\rho c_{p})}_{np,1}\right]+\phi_{2}{(\rho c_{p})}_{np,2}\) (2)