We present a detailed analysis and simulation of solenoidal, magnetically confined electron bombardment ion sources, aimed at molecular beam detection. The aim is to achieve high efficiency for singly ionized species while minimizing multiple ionization. Electron space charge plays a major role and we apply combined ray tracing and finite element simulations to determine the properties of a realistic geometry. The factors controlling electron injection and ion extraction are discussed. The results from simulations are benchmarked against experimental measurements on a prototype source.