С использованием теории термодинамически согласованных гиперболических систем законов сохранения сформулирована система определяющих уравнений двухфазного течения сжимаемой газожидкостной смеси в длинных трубах. Для полученных определяющих уравнений разработаны методы Рунге-Кутты-TVD и Рунге-Кутты-WENO высокого порядка точности. Применение сформулированной модели и численных методов для решения ряда тестовых задач продемонстрировало их эффективность и возможность использования для моделирования пробковых режимов течений.