Режим аккреционной колонки возникает при аккреции на замагниченную нейтронную звезду (НЗ) при больших темпах аккреции. Давление излучения внутри колонки приводит к формированию стоячей ударной волны над поверхностью НЗ. Стационарное аналитическое решение для течения ниже фронта ударной волны известно давно (Basko & Sunyaev 1976). Однако пределы применимости аналитической модели, процесс установления динамического и теплового равновесия в колонке, устойчивость ударной волны, моды переменности -- всё это изучено мало. \r\nПоэтому мы решили рассмотреть аккреционную колонку численно как одномерное гидродинамическое течение, зависящее от времени. Для этого написан оригинальный код, использующий упрощенную разностную схему годуновского типа. Для широкого диапазона параметров получается хорошее согласие с предсказаниями аналитической модели Баско-Сюняева. Наиболее важные различия возникают при больших темпах аккреции. При этом давление излучения начинает превышать магнитные напряжения на конечной высоте над поверхностью НЗ, что, по-видимому, должно приводить к потере массы через боковые стенки колонки. При сверхкритических темпах аккреции это, вероятно, будет приводить к формированию оттоков со скоростями ~0.1-0.5c. \r\nУстановление стационарного течения в наших моделях неизменно сопровождается затухающими колебаниями с периодами, соответствующими времени распространения звука от ударной волны до поверхности. В спектре мощности реального рентгеновского пульсара такая переменность должна приходиться на гектогерцовый диапазон.