Mean flow modeling in high-order nonlinear Schrödinger equations