Fluid flow with phase change heat transfer in a three-dimensional porous channel with asymmetrically heating from one side is numerically studied in this paper. The “modified” Kirchhoff method is used to deal with the spatial discontinuity in the thermal diffusion coefficient in the energy equation. The velocity and temperature fields, as well as the liquid saturation field on the heated section of the wall with different Peclet and Rayleigh numbers are investigated. The results show that the liquid flow bypasses the two-phase zone, while the vapor flows primarily to the interface between the sub-cooled liquid zone and the two-phase zone. An increase in the Peclet number decreases the two-phase region while an increase in the Rayleigh number helps to spread the heat to a larger region of the domain. The distribution of the liquid saturation on the heated section of the wall indicates that the minimum liquid saturation increases with the increase of both the Peclet and Rayleigh numbers