Optimal profiles in a phase-transition model with a saturating flux