We propose a novel breakup model to simulate the catastrophic breakup regime in a supersonic cross flow. A developed model has been extended from an existing Kelvin-Helmholtz/Rayleigh-Taylor (K-H/R-T) hybrid model. A new mass reduction rate equation, which has critical effects on overall spray structure, is successfully adopted, and the breakup length, which is an important parameter in existing model, is replaced by the breakup initiation time. Measured data from the supersonic wind tunnel with a dimension of 762×152×127 mm was employed to validate the newly developed breakup model. A nonaerated injector with an orifice diameter of 0.5 mm is used to inject water into a supersonic flow prescribed by the momentum flux ratio of the liquid jet to free stream air, q 0 . The conservation-element and solution-element (CE/SE) method, a novel numerical framework for the general conservation law, is applied to simulate the supersonic compressible flow. The spray penetration height and average droplet size along with a spray penetration axis are quantitatively compared with data. The shock train flow structures induced by the presence of a liquid jet are further discussed.