In this work, we carry out three-dimensional mesoscale simulations of heterogeneous solid propellant combustion. We solve the reactive low-Mach-number equations in the gas phase with complete coupling to the solid phase. The model takes into account thermal expansion and deformation in the solid phase by using a hypoelastic law in the quasi-static limit. To account for morphology, we select two different propellant formulations with different particle size distributions and present the results as a function of pressure. The thermomechanical behavior is accessed by examining quantities such as strain and stress in the propellant as a function of pressure and propellant morphology. We also show that temperature and velocity fluctuations exist in the far field above the propellant surface and that these fluctuations can be significant. To better understand the nature of these fluctuations, we vary the pressure and make relevant plots of normal velocity and temperature probability density functions, as well as time autocorrelations. Such descriptions are necessary to account for the coupling between the mesoscale and the macroscale, where the fluctuations at the mesoscale can affect quantities at the macroscale, such as head-end pressure, trigger parietal vortex shedding, and aeroacoustics.