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