Determination of failure envelope for faceted snow through numerical simulations