Analysis of hyperbolic conservation laws with random discontinuous flux functions and their efficient simulation