Жидкости в трубопроводах

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

Умножим второе уравнение системы уравнений (12.14) на и сложим результат с первым уравнением. Получим:

.

Аналогичным образом после вычитания второго уравнения системы (12.14), умноженного на , из первого, получим

.

Если на плоскости переменных рассмотреть прямые линии, которые определяются уравнениями

и которые называются характеристиками системы дифференциальных уравнений (12.13), можно заметить, что для любого параметра :

Это означает, что выражение, стоящее слева, есть производная от функции по направлению первой характеристики (или, как говорят, производная вдоль характеристики ). Аналогично

,

т.е. выражение, стоящее слева, есть производная от функции по направлению второй характеристики (или вдоль направления ).

Используя понятие производной по направлению, полученные уравнения можно записать в следующем виде:

или

. (12.16)

Система (12.16) называется характеристической формой системы уравнений (12.14).

Если , то правые части уравнений в (12.15) равны нулю. Это означает, что вдоль характеристики положительного наклона ( ) сохраняется величина , а вдоль характеристики отрицательного наклона ( ) сохраняется величина . Величины и называются инвариантами Римана.

Итак, при , т.е. при отсутствии силы трения и скатывающей составляющей силы тяжести, вдоль характеристик положительного наклона сохраняется первый инвариант Римана, а вдоль характеристики отрицательного наклона - второй инвариант Римана.

При величины и не сохраняются на соответствующих характеристиках. Однако формулы (12.16) могут служить для расчета различных неустановившихся течений в трубопроводе, особенно если использовать численные методы.

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

Рис. 12.6. Расчетная схема метода характеристик

Через узлы получившейся сетки проведем характеристики и положительного и отрицательного наклонов, соответственно. Непрерывное распределение искомых функций и заменим дискретными значениями и сеточных функций в узлах построенной сетки. Предположим, что все значения и известны в каком-нибудь слое и требуется найти значения сеточной функции при , т.е. и . Покажем, как это сделать на примере произвольной точки .

Заменим производные по направлению в уравнениях (16) конечными разностями вдоль характеристик и . Получим:

где

.

Отсюда получаем систему уравнений для определения давления и скорости жидкости в точке через известные величины этих параметров в точках и :

или

где и значения функции , вычисленные по параметрам точек и , соответственно. Из последней системы вычисляем значения давлений и скоростей потока нефти в трубопроводе в момент времени через значения этих же параметров в момент времени :

(12.17)

Таким образом, рекуррентные формулы (12.17) в принципе решают поставленную задачу о расчете неустановившихся течений слабо сжимаемой жидкости в трубопроводе, поскольку позволяют рассчитать значения давлений и скоростей течения в последующий момент времени по известным значениям этих параметров в предыдущий момент времени. Так как за первый “предыдущий” момент времени можно взять начальное состояние потока (т.е. значения давлений и скоростей течения в момент времени, принимаемый за начальный ), то вычисляя по формулам (12.17) шаг за шагом значения этих величин в последующие моменты времени, можно рассчитать параметры потока в произвольный момент времени и затем найти все интересующие нас технологические параметры нестационарного режима.