In this paper we demonstrate how a computational physics approach based on a discrete particle model can be employed for the visualization of magnetic field lines. This can be regarded as an illustration technique for the underlying vector field. To accomplish this, a model for ellipsoidal shaped particles consisting of ferromagnetic materials is described in combination with a parallel simulation algorithm based on an electrostatically enhanced Molecular Dynamics approach. For interactive viewing of the time-dependent simulation data, a hardware-based projection algorithm for the perspectively correct rendering of large numbers of illuminated ellipsoids is presented.