A numerical investigation of conjugate natural convection in a 2D square enclosure with multiple solid blocks, considering the impact of volumetric thermal radiation, was conducted using the lattice Boltzmann method (LBM). The effect of Planck number wall surface emissivity extinction coefficient solid-to-fluid volume fraction and the number of solid blocks on the fluid-structure and heat transfer within the enclosure were analyzed. The numerical model was validated using benchmark problems from the literature, showing good agreement. The results revealed that at and increasing wall surface emissivity significantly increased the Nusselt number, ranging from approximately 51-116% for different wall emissivity values, as compared to scenarios without surface radiation (). On the other hand, reducing χ, Pl, and β improved the heat transfer rate. The results also demonstrated that subdividing the initial block reduced the heat transfer rate depending on the number of blocks and the value of the Planck number. For instance, when 64 solid blocks were present in the cavity and the Nusselt number decreased by 59.36% compared to the case with a lower solid-to-fluid volume fraction (χ = 15%). The parametric analysis provided useful insights for engineering design in the modern construction industry, with possible applications in building block design and energy efficiency enhancements.