We propose high-resolution numerical schemes for two-phase compressible flow simulations to reproduce dynamically created gas/vapor-liquid interfaces with phase change. In the Godunov-type finite volume framework, suppressing numerical dissipation errors in numerical schemes is crucial for capturing the discontinuous solutions. The MUSCL scheme has second-order accuracy for smooth solutions and non-oscillatory behavior near discontinuous solutions. However, the MUSCL scheme introduces excessive numerical dissipation and diffuses the discontinuities nonphysically, leading to difficulties in distinguishing gas/vapor and liquid phases, and thus causing a blur in the interfaces during the multi-phase flow simulations. The hybrid-type boundary variation diminishing (BVD) scheme in this paper combines the MUSCL scheme and the THINC scheme to reduce the numerical dissipation errors near the discontinuities. The MUSCL-THINC-BVD scheme applies the MUSCL scheme for smooth solutions and the THINC scheme for discontinuous solutions, resulting in the successful capture of the discontinuities including the dynamically created gas/vapor-liquid interfaces. The Adaptive THINC-BVD scheme, which switches two types of THINC schemes with different values of gradient parameter, also captures the discontinuities clearly. The numerical results of the benchmark tests show that the proposed BVD schemes can lucidly reproduce the vapor-liquid interfaces newly created during the dynamical process of phase change. •Hybrid BVD schemes are presented to compute gas/vapor interfaces in 6-equation model for multi-phase flows.•Interfaces newly generated between liquid and vapor during the evaporation process are well-reproduced for the first time.•Basic benchmark validations are conducted for two-phase flows associated with phase changes.•The numerical methods presented are well-suited for simulating multi-phase flows involving phase changes.