A Computational Approach for Approximating Center Manifolds of Black-Box Simulators