A new multiscale method of primary atomization is proposed. The liquid blobs larger than the grid volumes were captured by the volume of fluid algorithm, the droplets comparable with the grid volumes or smaller were simplified to particles and tracked by a Lagrangian particle tracking model on a virtual mesh. Large eddy simulation with one equation sub-grid scale turbulent energy transport model was used to describe the turbulent flow. A new developed code was validated by several cases, and some key parameters were investigated to improve the precision. The primary atomization of a single jet was calculated, and the results of instantaneous and average characteristics are in good agreement with the experimental results.