The present work focuses on accurately estimating heat flux variations on the hot faces of the billet and slab molds during continuous casting using the Metropolis Hastings–Markov Chain Monte Carlo (MH-MCMC) Bayesian inverse method and experimental data. Accurate estimation of heat flux variation is crucial for operational efficiency and ensuring high-quality steel products. The variation in the hot face heat flux (qn(y)) of the billet mold is estimated using the proposed MH-MCMC Bayesian inverse method and experimental data for various casting speeds. The study extends this methodology to the slab mold for estimating the variations in the heat fluxes (qnr(y), qw(x,y)) on the hot narrow and wide faces of the slab mold. The estimated heat fluxes for both molds are used in respective forward models to obtain temperatures. The temperatures obtained excellently match experimental temperatures and show superior accuracy in estimating the heat fluxes. From the present methodology, the errors between these simulated and experimental temperatures are much less compared with the literature. The predicted heat flux variations on the hot faces of the billet and slab molds are also compared with the heat flux obtained from the literature.This study provides valuable insights into accurately estimating heat fluxes in continuous casting molds (billet and slab) using the MH-MCMC Bayesian inverse method. The accurately estimated heat fluxes of the molds are essential for ensuring the quality, efficiency, and innovation of continuous casting processes.