The main objective of this article is the exploitation of a stochastic hybrid mesh-free method based on stochastic generalized finite difference (SGFD), Newmark finite difference (NFD) methods and Monte Carlo simulation for thermoelastic wave propagation and coupled thermoelasticity analysis based on GN theory (without energy dissipation). A thick hollow cylinder with Gaussian uncertainty in mechanical properties is considered as an analyzed domain for the problem. The effects of uncertainty in mechanical properties with various coefficients of variations on thermo-elastic wave propagation are studied in details. Also, the time histories and distribution on thickness of cylinder of maximum, mean and variance values of temperature and radial displacement are studied for various coefficients of variations (COVs).