This study focuses on the analysis of an approach to the simulation of the phase transition in porous media when hot steam is injected into the oil reservoir. The reservoir is assumed to consist of a porous medium with homogeneous thermal properties. Its porous space is filled with a three-phase mixture of steam, water, and oil. The problem is considered in a non-stationary and non-isothermal formulation. Each phase is considered to be incompressible, with constant thermal properties, except for the dynamic viscosity of oil, which depends on the temperature. The 1D mathematical model of filtration, taking into account the phase transition, consists of continuity, Darcy, and energy equations. Steam injection and oil production in the model are conducted via vertical or horizontal wells. In the case of horizontal wells, the influence of gravity is also taken into account. The Lee model is used to simulate the phase transition between steam and water. The convective terms in the balance equations are calculated without accounting for artificial diffusion. Spatial discretization of the 1D domain is carried out using the finite volume method, and time discretization is implemented using the inverse (implicit) Euler scheme. The proposed model is analyzed in terms of the accuracy of the phase transition simulation for various sets of independent phases and combinations of continuity equations. In addition, we study the sensitivity of the model to the selected independent phases, to the time step and spatial mesh parameters, and to the intensity of the phase transition. The obtained results allow us to formulate recommendations for simulations of the phase transition using the Lee model.