Yttrium Aluminum Garnet (YAG, Y3Al5O12) and its varieties have applications in thermographic phosphors, lasing mediums, and thermal barriers. In this work, thermal properties of crystalline YAG where aluminum atoms are substituted with gallium atoms (Y3(Al1−xGax)5O12) are explored with molecular dynamics simulations. For YAG at 300K, the simulations gave values close to experimental values for constant-pressure specific heat, thermal expansion, and bulk thermal conductivity. For various values of x, the simulations predicted no change in thermal expansion, an increase in specific heat, and a decrease in thermal conductivity for x = 50%. Furthermore, the simulations predicted a decrease in thermal conductivity with decreasing system size.