A discontinuity-enriched finite element method for the computational design of phononic crystals