In this work, we extend a separate magnitude and phase regularization framework for Phase Contrast MRI by incorporating the divergence-free condition.