Numerical solver for supercontinuum generation in multimode optical fibers

We present an approach for numerically solving the multimode generalized nonlinear Schrödinger equation (MM-GNLSE). We propose to transform the MM-GNLSE to a system of first-order ordinary differential equations (ODEs) that can then be solved using readily available ODE solvers, thus making modeling of pulse propagation in multimode fibers easier. The solver is verified for the simplest multimode case in which only the two orthogonal polarization states in a non-birefringent microstructured optical fiber are involved. Also, the nonlinear dynamics of the degree and state of spectral polarization are presented for this case.