Laser-tissue interaction is a multi-physics phenomenon not yet mathematically describable and computationally predictable. It is a challenge to model the laser-tissue interaction for real time laser Benign Prostatic Hyperplasia (BPH) simulation which requires the laser-tissue interaction model to be computationally efficient and accurate. Under the consideration and enforcement of the thermodynamic first law and treating the laser-tissue interaction as a graybox, utilizing the sensitivity analysis of some key parameters that will affect the laser intensity on the tissue surface with respect to the tissue vaporization rate, a phenomenological model of laser-tissue interaction is developed. The developed laser-tissue interaction model has been implemented for a laser BPH simulator and achieves real time performance (more than 30 frames per second). The model agrees well with the available experimental data.