Method and non-transitory computer readable medium thereof for thermal analysis modeling

A non-transitory computer readable medium storing a computer program having code to execute a method for thermal analysis modeling, wherein the method comprises:

• establishing an electrothermal network π model on the basis of a plurality of electronic modules of an electronic system to define a heat source, at least one propagation path and a common base of the electronic system;
• defining a plurality of observation points in the electronic system, wherein each of the observation points is located at an isothermal surface enclosing a volume surrounding a reference point, where the reference point is the heat source or one of the observation points;
• calculating a heat conduction temperature difference and a heat convection temperature difference according to a power density function of the reference point, a thermal conductivity coefficient and a distance vector between the reference point and each of the observation points;
• and
• establishing a temperature distribution of the electronic system according to the heat conduction and the heat convection temperature differences and a defined temperature of the common base,
• wherein when the reference point is the heat source, the heat conduction temperature difference is calculated according to the power density function, the thermal conductivity coefficient and the distance vector,
• wherein the step of calculating the heat conduction temperature difference further comprises performing volume integration on the volume surrounding the reference point and surface integration on the isothermal surface enclosing the volume according to the power density function, the thermal conductivity coefficient and the distance vector,
• wherein the power density function comprises a volume power density function qv and a surface power density function qs, the thermal conductivity coefficient comprises a volume thermal conductivity coefficient kv and a surface thermal conductivity coefficient ks, the location of each of the observation points is r, the location of the heat source is r′, a unit normal vector of the isothermal surface is {circumflex over (n)}, the surface is a, and the heat conduction temperature difference is expressed as:
• ΔT=(¼π)*∫v(qv/kv| r− r′|)dv′−(¼π)*s(qs|{circumflex over (r)}−{circumflex over (r)}′|·{circumflex over (n)}/ks| r− r′|)da′.