**Abstract**

Surface heat flux is an important parameter in various industrial applications, which is often estimated based on measured temperature by solving inverse heat conduction problem (IHCP). In this chapter, the available IHCP methods including sequential function specification (SFS), transfer function (TF) and Duhamel's theorem were compared, taking the example of surface heat flux estimation during spray cooling. The Duhamel's theorem was improved to solve 1D multi-layer ICHP. Considering the significant nonuniformity of heat transfer, the 2D filter solution method was proposed to estimate surface heat flux for 2D multi-layer mediums. The maximum heat flux calculated by the 1D method was underestimated by 60% than that calculated by 2D filter solution, indicating that the lateral heat transfer cannot be ignored. The cooling performances based on 2D filter solution demonstrated that substituting the environment friendly R1234yf for R134a can remarkably reduce global warming potential to <1, but its cooling capacity is insufficient. The effective heat flux of R1234yf can be enhanced by 18.8% by reducing the nozzle diameter and decreasing the back pressure, providing the theoretical basis for the clinical potential substitution of R1234yf with low global warming potential (GWP) for commercial R134a with high GWP in laser dermatology.

**Keywords:** surface heat flux, inverse heat conduction problem, improved Duhamel's theorem, 2D filter solution, spray cooling
