Sparkling identifies spherical cosmic voids in numerical N-body simulations and galaxy surveys. The algorithm uses a Voronoi tessellation to estimate the density field and selects underdense candidates by iteratively computing the integrated density contrast inside spheres of increasing radius until a specific threshold is reached. Void centers are refined through repeated random displacements that are accepted only when they increase the candidate radius, yielding well-defined centers for subsequent analysis. Overlapping spheres are rejected to produce a non-redundant void catalog suitable for computing stacking statistics such as correlation functions and radial density profiles.