To study the temperature field distribution of engine cylinder head more precisely, based on the Euler multi-phase model, a new two-phase flow model combining the Rohsenow core boiling heat transfer model was established, which can be used for the numerical simulation of subcooled boiling. Applying BDL and Rohsenow transfer models separately, the cylinder head fluid-solid coupling heat transfer system was simulated and analyzed with the CFD technology. The results have shown that the Rohsenow core boiling heat transfer model is more precise than the BDL model in calculating the temperature field of engine cylinder head. The highest temperature of the cylinder head is 192.22 ℃, and the maximum of stress is 235 MPa, both of which can meet the strength requirements of cylinder head materials.