Direct Monte Carlo and multifluid modeling of the circumnuclear dust coma