Relatively accurate component (and board) thermal simulations?

After playing around with SkiDL+NgSpice on SkiDL+NgSpice (via PySpice): simulating a simple parallel resistors circuit's power outputs, I wonder what’s the recommended/best/suitable way to go about simulating thermodynamics on a circuit. Talking with an expert colleague he recommended the following theoretical approach:

Essentially I would like to have a way to “heat map” a circuit pretty much like a FLIR camera does but via simulation/software… is it possible with NgSpice alone? Or would this require more ammo with FEA modelling or the like?

I would like to know if somebody has attempted this successfully and has good experience/code to share? :slight_smile:

Just found this, sweet:

Which goes to: