We present a novel simulation approach combined with pulsed laser experiments, spectroscopic ellipsometry, and Raman spectroscopy to comprehensively analyze phase transformation dynamics in thin films. The simulations apply to any thin film stack and incorporate critical factors, such as thin film interference, heat transfer, and temperature-dependent optical properties during heating and melting. As a case study, we investigate the picosecond laser-induced amorphization of antimony sulfide (Sb2S3) thin films, a promising alternative to traditional phase-change materials in photonic applications to validate the simulation model. The computational efficiency of our simulations enables not only the investigation of the laser-induced phase transformation but also the optimization of key process parameters and parameter fitting. The simulations identified optimal film thickness and laser fluence parameters that maximize energy efficiency, melting effectiveness, and quenching rate while ensuring high reflectivity contrast between the amorphous and crystalline states. By constructing a wide-ranging, high-resolution parameter map of the laser fluence and film thickness dependence of the melting process, we demonstrate how this model guides the understanding of phase transformation dynamics. Raman spectroscopy confirms the polycrystalline to amorphous transition of Sb2S3 and provides a semiquantitative estimate of the amorphous fraction as a function of laser fluence, which is qualitatively consistent with the simulation predictions of the model. The open-source simulation framework, experimentally validated, provides valuable insights into laser-induced amorphization dynamics in Sb2S3 and related phase-change material thin films, enabling rapid optimization of photonic devices.