Hey there,
You applied those boundary conditions to the same face, so that’s not going to work out very well. I copied your project and applied the surface heat flux to the other face. Note that heatflux has to be negative if heat is leaving the domain.
The result is exactly -4ºC as you predicted:
