Connectivity provided from space is the ultimate method for a true and uninterrupted global coverage, where there are several massive satellite constellations that are currently undergoing rapid deployment. This paper lays a tractable analytic approach for modeling the downlink coverage probability in such dense satellite networks. The presented framework utilizes tools from stochastic geometry leveraging recent developments in Satellite-to-Ground path-loss modeling, and in Line-of-Sight probability formulation. Simulation results show a very good match to the presented analytic forms even for lower number of satellites, suggesting that the obtained insights can be used by system designers to analytically predict and optimize the expected downlink coverage performance.