From an analytical standpoint, a majority of calculations use known boundary conditions (temperature or flux) and the so-called direct route to determine internal temperatures, strains, and/or stresses. For such problems where the thermal boundary condition is known a priori, the analytical procedure and solutions are tractable for the linear case where the thermophysical properties are independent of temperature. On the other hand, the inverse route where the boundary conditions must be determined from remotely determined temperature and/or flux data is much more difficult mathematically, as well as inherently sensitive to data errors (i.e., ill-posed). When solutions are available, they are often restricted to a harsh, albeit unrealistic step change in temperature or flux and/or are only valid for relatively short time frames before temperature changes occur at the far boundary. While the two approaches may seem to be at odds with each other, a generalized direct solution based on polynomial temperature or strain-histories can also be used to determine unknown boundary conditions via least-squares determination of coefficients. Once the inverse problem (and unknown boundary condition) is solved via these coefficients, the resulting polynomial can then be used with the generalized direct solution to determine the thermal- and stress-states as a function of time and position. When used for both thick slabs and tubes, excellent agreement was seen for various test cases. In fact, the derived solutions appear to be well suited for many thermal scenarios, provided the analysis is restricted to the time interval used to determine the polynomial and the thermophysical properties that do not vary with temperature. Since temperature dependent properties can certainly be an issue that affects accuracy in these types of calculations, some recent analytical procedures for both direct and inverse solutions are also discussed.