A meshless particle code based on the Lagrangian Smoothed Particle Hydrodynamics method has been developed with the objective of simulating primary atomization. It is capable of predicting compressible and incompressible fluids. Single phase flows as well as multiphase flows have been simulated successfully. Furthermore, surface tension effects can be considered. Preliminary results will be presented, which underline the promising performance of this approach and which are the pre-stage towards simulating primary atomization.