An intrusive hybrid method for discontinuous two-phase flow under uncertainty

The discretization of the non-smooth region is based on a generalization of the HLL flux, and have many properties in common with its deterministic counterpart. It is simple and robust, and captures the statistics of the shock. The discretization of the smooth region is carried out with high-order finite-difference operators satisfying a summation-by-parts property.