There is the ``externalWallHeatFluxTemperature`` which allows to directly define a power or heat flux on your patch: https://doc.openfoam.com/2306/tools/processing/boundary-conditions/rtm/derived/thermal/externalWallHeatFluxTemperature/ Also, do you have under-relaxation on energy (e or h) in ``fvSolution``?