1. A computer-implemented method for enhancing thermal design of a system having a plurality of boundary values, said method comprising the steps of: representing thermal response of said system to said plurality of boundary values as a superposition of temperature fields associated with given ones of said boundary values; obtaining at least one constraining parameter; and determining at least one of spatial and temporal distribution of said plurality of boundary values, based on said thermal response represented in said representing step, to satisfy said at least one constraining parameter.