Recent interest in human-scale Mars missions has motivated research into the effects of dust particles that are suspended in the Martian atmosphere on reentry vehicles. For instance, these dust particles can enhance erosion of thermal protection systems and amplify surface heat fluxes. As such, in this work, we perform numerical simulations of hypersonic dusty flows over a blunt body using a two-way-coupled Lagrangian particle method in a discontinuous Galerkin framework. Due to the lack of a well-established physical model of the disperse phase appropriate for this multiphase flow regime, we aim to evaluate the sensitivities of numerically predicted dust-induced heat flux augmentation to various components and parameters of the model. Specifically, we investigate the drag coefficient correlation, the Nusselt number correlation, different terms in the particle momentum and energy equations, and particle size.