How to include the multiple heat sources in constraints

Dear Developer,

First of all, thank you very much for providing Modulus-sym! I have encountered a problem that I am unable to solve on my own, and I would like to seek your advice as well as insights from friends on the forum:

I have studied several examples in the tutorials, and I noticed that they all involve only a single heat source. In such cases, the heat flux density or heat flow can be converted into a temperature gradient on the heat source surface using Fourier’s law of heat conduction. However, my model involves multiple surface heat sources. Since each heat source is influenced by the others, the heat flux passing through each surface heat source is unknown. As a result, I am unable to calculate the temperature gradient for each heat source surface. I only know the original heat flux density for each individual heat source. In this situation, where multiple heat sources are coupled, how should I handle the heat sources? In what form should they be included in the constraints?

I would greatly appreciate your and the community’s feedback. Thank you in advance!