The outwardly-opening piezoelectric injector is gaining popularity as a high efficient spray injector due to its precise control of the spray. However, few modeling studies have been reported on these promising injectors. Furthermore, traditional linear instability sheet atomization (LISA) model was originally developed for pressure swirl hollow-cone injectors with moderate spray angle and toroidal ligament breakups. Therefore, it is not appropriate for the outwardly-opening injectors having wide spray angles and string-like film structures. In this study, a new spray injection modeling was proposed for outwardly-opening hollow-cone injector. The injection velocities are computed from the given mass flow rate and injection pressure instead of ambiguous annular nozzle geometry. The modified Kelvin-Helmholtz and Rayleigh-Taylor (KH-RT) breakup model is used with adjusted initial Sauter mean diameter (SMD) for modeling breakup of string-like structure. Spray injection was modeled using a Lagrangian discrete parcel method within the framework of commercial CFD software CONVERGE, and the new model was implemented through the user-defined functions. A Siemens outwardly-opening hollow-cone spray injector was characterized and validated with existing experimental data at the injection pressure of 100 bar. It was found that the collision modeling becomes important in the current injector because of dense spray near nozzle. The injection distribution model showed insignificant effects on spray due to small initial droplets. It was demonstrated that the new model can predict the liquid penetration length and local SMD with improved accuracy for the injector under study.