In the framework of further development of our FEATFLOW based in-house simulation tools we have achieved improvements on the solver level such that reasonably large jumps of physical quantities can be tackled without causing serious convergence problems. The underlying cell-averaging technique of physical properties plays a key role in the design of the embedded linear multigrid driver in order to suppress numerical oscillations in the solution process. The Level Set method - serving as our interface capturing method - has been successfully enriched with a robust hyperbolic PDE based reinitialization technique in order to keep the Level Set indicator function smooth and close to a distance function. Furthermore, taking advantage of a Discontinuous Galerkin FEM approach for the Level Set equation offers at one hand the possibility of reduction of the computational domain only to the subset of elements located in the vicinity of the interface thus reducing the computational costs, and on the other hand no smoothing of the discontinuous sign function is required. The mentioned treatments have been realized according to the so called fictitious boundary method. Initially, the newly designed simulation tool has been subjected to validation against experimental results in case of moderate operating conditions. Subsequently, more realistic nozzle configurations and more dramatic operating conditions have been analyzed. Finally, the influence of geometry and inflow modulation on the achieved droplet sizes and droplet generation frequencies has been investigated.