def METHOD_NAME(run_in_tmpdir, void_coincident_geom_model):
    # update source to a number of points on the outside of the sphere
    # with directions pointing toward the origin
    n_sources = 20
    phi_vals = np.linspace(0, np.pi, n_sources)
    theta_vals = np.linspace(0, 2.0*np.pi, n_sources)
    bbox = void_coincident_geom_model.geometry.bounding_box
    # can't source particles directly on the geometry boundary
    outer_r = bbox[1][0] - 1e-08
    sources = []
    energy = openmc.stats.Discrete([14.06e6], [1])
    for phi, theta in zip(phi_vals, theta_vals):
        src = openmc.IndependentSource()
        src.energy = energy
        pnt = np.array([np.sin(phi)*np.cos(theta), np.sin(phi)*np.sin(theta), np.cos(phi)])
        u = -pnt
        src.space = openmc.stats.Point(outer_r*pnt)
        src.angle = openmc.stats.Monodirectional(u)
        # set source strengths so that we can still expect
        # a tally value of 0.5
        src.strength = 0.5/n_sources
        sources.append(src)
    void_coincident_geom_model.settings.source = sources
    sp_filename = void_coincident_geom_model.run()
    _check_void_spherical_tally(sp_filename)