Жидкости в трубопроводах
Расчет неустановившихся течений жидкости в трубопроводах более сложен, чем расчет установившихся течений. Не существует простых алгебраических формул для вычисления параметров течения в неустановившихся течениях хотя бы потому, что таких течений существует бесчисленное множество. Поэтому расчеты неустановившихся течений, например тех, которые были поименованы в начале этого параграфа, осуществляют численно с использованием компьютеров. Приведем один из эффективных алгоритмов расчета, основанный на математическом методе, получившем название метода характеристик.
Умножим второе уравнение системы уравнений (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) шаг за шагом значения этих величин в последующие моменты времени, можно рассчитать параметры потока в произвольный момент времени и затем найти все интересующие нас технологические параметры нестационарного режима.
Дата добавления: 2016-05-16; просмотров: 799;