We present a framework to simulate nuclear magnetic resonance (NMR) signal enhancement experiments using parahydrogen at ultralow magnetic fields. Our approach captures the full quantum evolution of the system while reducing computational complexity for multi-spin systems through symmetry-based dimensionality reduction. This enables efficient prediction of optimal polarization transfer fields and simulation of the resulting ultralow-field NMR spectra.