Under consideration is the finnite-size scaling of effective thermoelastic properties of random microstructures from a Statistical Volume Element(SVE) to a Representative Volume Element (RVE), without invoking any periodic structure assumptions, but only assuming the microstructure's statisticsto be spatially homogeneous and ergodic. The SVE is set up on a mesoscale,i.e. any scale finite relative to the microstructural length scale. The Hill condition generalized to thermoelasticity dictates uniform Neumann and Dirichletboundary conditions, which, with the help of two variational principles, lead toscale dependent hierarchies of mesoscale bounds on effective (RVE level) properties: thermal expansion and stress coefficients, effective stiffness, and specificheats. Due to the presence of a non-quadratic term in the energy formulas,the mesoscale bounds for the thermal expansion are more complicated thanthose for the stiffness tensor and the heat capacity. To quantitatively assessthe scaling trend towards the RVE, the hierarchies are computed for a planarmatrix-inclusion composite, with inclusions (of circular disk shape) located atpoints of a planar, hard-core Poisson point field. Overall, while the RVE isattained exactly on scales infinitely large relative to the microscale, depending on the microstructural parameters, the random fluctuations in the SVEresponse may become very weak on scales an order of magnitude larger thanthe microscale, thus already approximating the RVE.