Near Earth Objects (NEOs) enter the Earths atmosphere on a regular basis. Depending on the size, object and entry parameters; these objects can burn-up through ablation (complete evaporation), undergo fragmentation of varying nature, or impact the ground unperturbed. Parameters that influence the physics during entry are either unknown or highly uncertain. In this work, we propose a probabilistic approach for simulating entry. Probabilistic modeling typically requires an expensive Monte Carlo approach. In this work, we develop and present a novel engineering approach of developing surrogate models for simulation of the atmospheric entry accounting for drag, ablation, evaporation, fragmentation, and ground impact.