Abstract: A reliable multiphase flow simulator is an important tool to improve wellbore integrity and production decision-making. To develop a multiphase flow model with high adaptability and high accuracy, we first build a multiphase flow database with 3561 groups of data and developed a drift closure relationship with stable continuity and high adaptability. Second, a high-order numerical scheme with strong fault capture ability is constructed by effectively combining MUSCL technology, van Albada slope limiter and AUSMV numerical scheme. Finally, the energy equation is coupled into the AUSMV numerical scheme of the drift flow model in the form of finite difference. A transient non-isothermal wellbore multiphase flow model with wide applicability is formed by integrating the three technologies, and the effects of various factors on the calculation accuracy are studied. The accuracy of the simulator is verified by comparing the measurement results with the blowout experiment of a full-scale experimental well.