Nice, it turned out to scale almost linearly with the hardcoded values I made up.
Had to tweak the linear relation to make powerful engines stand out more.
powerFactor = max_energy_usage / 10K -- 10K means 600kW
if isShip(....) then powerFactor /= 3.33
powerFactor = math.pow(powerFactor, 1.25)