12. Течение в ступени центробежного насоса. Уравнения и методы решения. Сравнение эксперимента с численными результатами по интегральным характеристикам.

Цель работы : построение математической модели, описывающей течение флюида в ступени центробежного насоса.

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

• учесть вязкость и сжимаемость перекачиваемого флюида, турбулентность потока;

• учесть потери на трение и вихреобразование в элементах ротора и статора ступени насоса;

• получить в общем виде систему дифференциальных уравнений для описания течения сжимаемой вязкой жидкости в ступени ЦН.

На основе численных реализаций построенной математической модели :

• используя полученные распределения, сравнить теоретические и экспериментальные данные;

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

П остановка задачи:

• Ротор вращается со скоростью Ω

• Среда – однофазная, сжимаемая, вязкая

• Течение – развитое турбулентное

Численное моделирование турбулентности:

Турбулентность - это "…завихреное течение с очень большим числом возбужденных степеней свободы и с хаотическим распределением дисперсионных соотношений и фазовых сдвигов”*.

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

2. Моделирование методом крупных вихрей.

Гипотеза: статистические характеристики крупномасштабных турбулентных движений не зависят от молекулярной вязкости.

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

- Фильтрация мелкомасштабных пульсаций.

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

3. Решение систем уравнений, усредненных по Рейнольдсу.

Статистический подход к исследованию турбулентности - осреднение уравнений движения по ансамблю состояний. Все параметры движения разлагаются на среднюю и турбулентную составляющие.

Рассмотрим уравнения Навье-Стокса для несжимаемой вязкой среды .

Входящие в них величины представляются в виде сумм средних полей и пульсаций .

Подставив разложения и проведя осреднение, получим уравнения Рейнольдса

где - тензор напряжений Рейнольдса.

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

- зависимость Буссинеска

где - коэффициент турбулентной вязкости

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

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

Методы исследования течения среды в ступени насоса:

• Исходные уравнения, представляющие собой усреднённые по Рейнольдсу уравнения Навье-Стокса, записываются в дивергентной форме для цилиндрической системы координат (х, r, φ), вращающейся с угловой скоростью ротора Ω.

• Среда предполагается однофазной сжимаемой вязкой, подчиняющейся двухпараметрическому уравнению состояния.

• Течение – турбулентное.

• При расчетах течения в не вращающихся областях проточной части насоса (элементов статоров, направляющих устройств и др.) удобно использовать абсолютную цилиндрическую систему координат, а для вращающихся областей проточной части (элементов ротора)–относительную систему координат, вращающуюся с соответствующей угловой скоростью ротора .

Уравнения Навье-Стокса в трехмерной нестационарной постановке:

В безразмерном виде система уравнений во вращающейся цилиндрической системе координат имеет вид:

t – время, ρ – плотность, p – давление, T – температура, S – энтропия, е* – внутренняя энергия,

- –компоненты относительной скорости вдоль осей х, r и φ

– коэффициенты молекулярной и турбулентной вязкости,

и - числа Прандтля; Re - число Рейнольдса

Модель турбулентности

В безразмерном виде система уравнений во вращающейся цилиндрической системе координат имеет вид:

где - кинетическая энергия пульсаций

Здесь d_n расстояние вдоль нормали до стенки. Значения констант, входящих в последние формулы, следующие:

Начальные условия:

Начальными условиями задачи являются заданные в момент времени t=0 приближенные распределения гидродинамических параметров по всему рабочему пространству ступени насоса.

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

Необходимым условием решения задачи является задание полной 3D геометрии всей проточной части.

Граничные условия:

Граничные условия задачи:

• на неподвижных стенках все компоненты скорости задаются равными нулю – условие прилипания,

• на подвижных стенках скорость = Ω , остальные компоненты скорости равны нулю;

• на входе в насос задаются полные параметры набегающего потока, а также направление скорости;

• на выходе - статическое давление окружающей среды.

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

Исходные данные:

Численный эксперимент проводился для ступени электроцентробежной погружной насосной установки для нефтедобычи - ЭЦН АКМ.

Параметры расчета:

• Рабочее тело – вода : плотность ρ = 997 кг/м3, вязкость m = 0.001003 кг/м×с, температура Т = 300 К.

• Рассматривается изотермический процесс.

• Частота вращения ротора насоса: ω = 10000 об/мин.

В результате расчетов обнаружены существенные различия в характере течения при малых и больших расходах через ступень.

В связи с этим анализ результатов будет приведен для двух предельных значений расходов, при которых проводился численный эксперимент: 30 м3/сут и 160 м3/сут.

Интегральные характеристики:

Расходно-напорная характеристика ступени ЦН

1 - эксперимент, 2 – расчет на алгебраической одномерной модели, 3 - расчет на математической модели.

Относительная погрешность эксперимента составляет 15%.

Анализ результатов моделирования:

Рабочий агент: дегазированная нефть,

• плотность: 921 кг/м3,

• сжимаемость: 1,3*10^-3 1/МПа,

• вязкость:

• 0.005 Па*с,

• 0.015 Па*с,

• 0.025 Па*с.

Расходно-напорные характеристики при различных значениях вязкости

перекачиваемого флюида

Выводы:

1. Предложено математическое описание течения флюида, основанное на уравнениях Навье-Стокса с учетом вязкости флюида и развитой турбулентности потока.

2. Разработан алгоритм решения полученной системы дифференциальных уравнений, который позволил улучшить сходимость численного процесса.

3. Для проведения численных экспериментов создана твердотельная модель ступени насоса, который используется в нефтяной промышленности,.

4. Сравнение интегральных характеристик ступени, полученных численными методами, с экспериментальными результатами продемонстрировало их хорошее согласование.

5. Проведенный анализ многовариантных расчетов течения флюида в диапазоне значений расходов 30 – 160 м3/сут показал, что

• в случае больших величин объемных расходов основные потери идут на вихреобразование, вызванное инжекцией высокоскоростной струи в полость статора, а в случае малых объемных расходов потери связаны с возникновением большого количества малых вихрей во входном канале полости диффузора;

• особенности геометрии входного канала полости статора неминуемо приводят к его запиранию при малых значениях объемного расхода;

• особенности геометрии полости диффузора приводят к возникновению системы присоединенных вихрей, занимающих большую часть объема полости и существенно уменьшающих площадь проходного сечения;

• применив комплекс мер по перепрофилированию статорной части ступени центробежного насоса, можно существенно улучшить его гидродинамические характеристики и увеличить эффективность работы насоса в целом. 6.

6. Проведен анализ зависимости течения флюида в ступени от вязкости в диапазоне значений 0.001 – 0.025 мПа*с, который показал, что с увеличением вязкости:

• высота напора уменьшается;

• коэффициент полезного действия снижается;

• область оптимальных подач смещается в область меньших расходов.