Проведено исследование развития магниторотационной неустойчивости (МРН) в
аккрецирующей оболочке протозвезды в неизотермическом случае, т.е. с учетом
уравнения энергии, RKDG методом второго порядка с использованием распараллеливания по технологии MPI. Показано, что учет уравнения энергии в общей системе МГД-уравнений не оказывает существенного влияния на моделирование процесса образования и эволюции крупномасштабной структуры турбулентности,
приводящей к отводу углового момента на периферию оболочки и, как следствие, к
росту темпа аккреции на протозвезду. Результаты моделирования качественно
сравниваются с наблюдениями области звездообразования G31.41+0.31 и другими
наблюдательными данными.