The paper is focused on numerical identification of 2D temperature fields in flow boiling of the liquid through a horizontal minichannel with a rectangular cross-section. The heat transfer process in the minichannel is described by a two-dimensional energy equation with the corresponding boundary conditions. Liquid temperature is determined using the homotopy perturbation method (HPM) with Trefftz functions for Laplace'a equation. The numerical solution to the energy equation found with the HPM is compared with the solution obtained for the simplified form of the energy equation. Considering that only the thermal sublayer is taken into account, both solutions give similar results.