Stable numerical methods for conservation laws with discontinuous flux function