Проведено численное моделирование нелинейных спиновых волн в графеновых структурах псевдоспектральным методом Чебышёва и методом Ритца-Галеркина. Кроме того, реализовано численное моделирование нелинейных спиновых конфигураций плотности валентных электронов на поверхности графеновой пленки, таких как кинки, антикинки, бризеры и их взаимодействия. Разработаны численные методы моделирования потенциальной энергии взаимодействующих метастабильных состояний (бризеров) псевдоспиновых волн на поверхности графеновой пленки. Численно исследован дискретный спектр взаимодействия бризеров с целью установления инверсной заселенности уровней метастабильных состояний псевдоспиновых волн на поверхности графеновых пленок различной топологии.