Based on dissipative particle dynamics (DPD) methods and experimental data, we used an empirical relationship between the DPD temperature and the real temperature to build a model that describes the viscosity of molten TNT fluids. The errors in the predicted viscosity based on this model were no more than 2.3 %. We also studied the steady-state shear rheological behavior of molten TNT fluids containing nanoparticles ("nanofluids"). The dependence of the nanofluid viscosity on the temperature was found to satisfy an Arrhenius-type equation, η = Ae (B/T) , where B, the flow activation energy, depends on particle content, size, and shape. We modified the Einstein-type viscosity model to account for the effects of nanoparticle solvation in TNT nanofluids. The resulting model was able to correctly predict the viscosities of suspensions containing nano- to microsized particles, and did not require any changes to the physical background of Einstein's viscosity theory. Graphical Abstract The revised Einstein viscosity model that correctly predict the viscosity of TNT suspensions containing nanoparticles.