Text this: Phase-field modeling of coupled bulk photovoltaic effect and ferroelectric domain manipulation at ultrafast timescales