A reduced model of neutral-plasma interactions in the edge and scrape-off-layer: Verification comparisons with kinetic Monte Carlo simulations