The paper presents an analytical two-dimensional model of two-phase turbulent jets with focus on fuel sprays in internal combustion engines. The developed model allows prediction of the fuel spray parameters including local fuel concentration and mixture velocity. The model proposed in this paper is based on the single-phase steady-state laminar axisymmetric jet flow field solution by Schlichting. This solution is amended to include transport of the discontinuous fuel phase in a stagnant air in the limit of a dilute fuel concentration. This two-phase jet flow model admits a closed form analytical solution for the fuel concentration distribution. This solution is then applied to turbulent jet flow as per the approach described by Schlichting and in other studies, and used to predict point-wise properties of fuel sprays in internal combustion engines. The results of model simulations are compared with the available experimental data. It was found that the analytical model predicts satisfactorily spray properties without additional assumptions or fitting coefficient. The advantage of the developed model is in its simplicity and rigor in comparison with empirical and numerical models.