Научная статья на тему 'Теоретический анализ обтекания решетки профилей малорасходного центробежного насоса потоком вязкой жидкости'

Теоретический анализ обтекания решетки профилей малорасходного центробежного насоса потоком вязкой жидкости Текст научной статьи по специальности «Механика и машиностроение»

CC BY
201
78
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГИДРОДИНАМИКА / МАЛОРАСХОДНЫЙ ЦЕНТРОБЕЖНЫЙ НАСОС / РАБОЧЕЕ КОЛЕСО / ПОГРАНИЧНЫЙ СЛОЙ / СТЕПЕНЬ ДИФФУЗОРНОСТИ / HYDRODYNAMICS / LOW-FLOW CENTRIFUGAL PUMP / IMPELLER / BOUNDARY LAYER / DIFFUSER DEGREE

Аннотация научной статьи по механике и машиностроению, автор научной работы — Краева Екатерина Михайловна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Theoretical analysis of flow of cascade of low-flow centrifugal pump by viscid fluid stream

In view of design features of air-gas channel of impeller of low-flow centrifugal pump the analysis off low of cascade by viscid fluid flow is executed. Calculation of parameters of a boundary layer is executed.

Текст научной работы на тему «Теоретический анализ обтекания решетки профилей малорасходного центробежного насоса потоком вязкой жидкости»

E. N. Garin

THE INFLUENCE OF CONDITIONS FOR PASSING GLONASS AND GPS SATELLITE RADIO NAVIGATION SIGNALS ON THE ERRORS OF DEFINING RELATIVE COORDINATES

The influence of the conditions for passing GLONASS and GPS satellite radio navigation signals on the errors of defining relative coordinates is overlooked.

Keywords: signal delay, multipath propagation, error, frequency, navigation spacecraft.

© Гарин Е. Н., 2010

О

WKО* = J (W - Wx )dn,

(1)

где Wx - текущее значение скорости потока в пограничном слое.

Так как при п ; Wx ® Wк подынтегральное выражение стремится к нулю, тогда получим

УДК621.45 - 181.4:629.78

Е. М. Краева

ТЕОРЕТИЧЕСКИЙ АНАЛИЗ ОБТЕКАНИЯ РЕШЕТКИ ПРОФИЛЕЙ МАЛОРАСХОДНОГО ЦЕНТРОБЕЖНОГО НАСОСА ПОТОКОМ ВЯЗКОЙ ЖИДКОСТИ*

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

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

В анализе принято значение степени диффузорности канала рабочего колеса (РК) насоса, при котором на рабочей (напорной) стороне лопатки реализуется безотрывное обтекание профиля, а отрыв с развитием спутного потока реализуется на тыльной стороне профиля лопатки. Эта особенность гидродинамики РК МН была экспериментально изучена при визуализации и испытаниях [1].

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

К параметрам, характеризующим пограничный слой, как правило, относят [2] не саму толщину пограничного слоя 5 , в пределах которой величина скорости потока Wx отличается на 1 % от скорости ядра Wк, а толщину вытеснения динамического пограничного слоя 5 , соответствующую расстоянию, на которое сдвигаются от профиля лопатки линии тока внешнего течения вследствие образования пограничного слоя. Это позволяет при определении расхода вязкой жидкости принимать поток невязким, но считать, что он оттеснен от поверхности лопатки на величину 5*. Значение толщины вытеснения находим из выражения [2]

dn

или

WKО* = J W - Wx )dn = J (WK - Wx )dn. (2)

Аналогично для толщины потери импульса

Wt2 o'* = J Wx (Wx - Wx )dn.

(3)

Для анализа, в силу осесимметричности степени сдвига потока в радиальной плоскости по шагу решетки, действительный трехмерный поток в РК МН условно разделяем на два: внешнее ядро потока (вне пограничного слоя) и пограничный слой, в котором проявляются силы вязкости.

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

' Работа выполнена при поддержке грантов: АВЦП РПН ВШ № 2.1.2/802.

Применительно к РК лопаток одинарной кривизны, образующая которых параллельна вектору угловой скорости ю . Профили круговых решеток в полярных координатах Я, 9 и их конформное отображение на плоскостях х,у приведены на рис. 1. Используя граничные условия х = 0, Я = Я1 и х = £, Я = Я2, получим связь между текущими координатами прямой решетки х, у и полярными координатами Я, 9 профиля круговой решетки:

Я = Я1 ехр Х 1п—;

І Яі

9 = • _!2!_ =

Ґ ’ Ґ

пр пр

І„

і Я2 г 1п 2

Яі

(4)

2—

Рис. 1. Схемы решеток профилей при их конформном перестроении: а - круговые; 6 - прямые

Сформулируем основные положения принятой модели течения при расчете закрытого РК насоса.

1. Абсолютная скорость потока на входе в РК постоянна и направлена параллельно оси вращения РК.

2. Сдвиговые течения в плоскости вращения по дискам РК существенны и их необходимо учитывать.

3. Средняя скорость потока в решетке постоянна и с учетом сдвига равна скорости на выходе РК ^к = W2).

4. Форма входной кромки лопатки представлена на рис. 2, а обтекание ее соответствует натеканию потока вязкой жидкости окрестности клинообразного тела с критической точкой Е на вершине.

5. Влиянием пограничных слоев вдоль каналов и начальным сдвигом потока в меридиональной плоскости на основной поток пренебрегаем.

6. Характеры течений на покрывном и ведущем дисках РК при IV <5 отличаются незначительно, что допускает принять плоскость симметрии потока по средней линии (по ширине) межлопаточного канала (рис. 2). Поэтому при расчете пограничных слоев на лопатках (с учетом сдвиговых течений по дискам) учитывается симметричность параметров вторичных течений.

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

Наличие градиента давления по шагу межлопаточного канала приводит к поперечному сдвиговому течению, что вызывает изменение характера потока в пограничном слое вдоль дисков РК. Большая относительная толщина пограничного слоя на дисках РК с учетом сдвигового течения ядра потока в радиальной плоскости, предопределяет условное разделение потока в канале РК на ядро со скоростью (I ’ и вязкий пограничный слой, рассчитываемый с использованием интегрального уравнения импульсов [2]. Этот метод достаточно точен для инженерной практики и позволяет сравнить результаты расчета с опытом, учитывая более надежное определение в РК МН опытных интегральных, а не локальных параметров потока.

Течение на входном участке в РК рассмотрим с учетом особенностей поворота потока в окружности крити-

Рис. 2. Расчетная схема течения в РК центробежного насоса 97

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

Точное решение такой задачи для осесимметричного течения ввиду его аффинно-подобия вблизи критической точки выполнено Н. Фреслингом [2]. Расчет параметров осесимметричного течения на входе в РК сводится к расчету пограничного слоя плоского течения около кли-

л

на с углом раствора у = лу = [2]. В этом случае ско -

рость потенциального течения вдоль ограничивающего диска имеет вид [2]

Wx = схт,

где с и m - постоянные величины.

Тогда m рассчитывается по формуле

У

У

2 -у

2я-у

(5)

(6)

По оси г соответствующее значение толщины пограничного слоя определяется по выражению [2]

= 4,4 -.

(7)

скоростей

wx

w2

= 1 при x = 1пр имеем Wx

c = -

l;

Уравнение (8) запишем следующим образом:

W„ cosє„

w = -

і;

(9)

x=

l( x )(l„R - lR) l R^

(1О)

wx=wK п;

lR - lR

cos є

R

(11)

7 = l (x) где 'пр = /

l ; l(x) - текущее значение длины канала.

/ пр

При g = 0 имеем случай продольного обтекания пластины со сдвигом. Тогда из уравнения (5) m = 0 и величина относительно скорости на выходе из РК МН имеет вид

W(x)2 = W C0S eR = W2. (12)

Будем считать, что по всей поверхности дисков РК пограничный слой ламинарный (в последствии это докажем). Решение уравнения пограничного слоя проводим с учетом градиентного течения, используя для оценки

0'*

закона изменения толщины потери импульса 0zl вдоль ограничивающих дисков межлопаточного канала выражение из работы [2], сведенное к простой квадратуре:

\2 О,47п

(dz) w 6

n( x ) О

(13)

С учетом уравнения (11), решая относительно 5 имеем

5” = О,686 х

піп;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

R,

(5; +1) [lR - lR ];WK cos єr2 l (x);

(14)

Из уравнения (7) следует, что вблизи критической точки величина, на которую распространяется влияние трения, пропорциональна корню квадратному из кинематической вязкости и с учетом ю составляет малое значение. Следует заметить, что безразмерное распределение

При m = О; єr = О и l(x) = l уравнение (14) принимает вид

5z1 = О,686(

vk

w, ''

(15)

и значение 5z1 не зависят от координа-

ты х, т. е. не изменяются вдоль радиуса диска вплоть до начала лопаток на радиусе R1. Это позволяет проводить расчет течения вдоль ограничивающих дисков межлопа-точных каналов РК МН от радиуса R1 и далее по каналу, принимая Oz0 = Ог1 = const.

Для каналов с гидродинамической степенью диффу-

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

Для оценки толщины вытеснения пограничного слоя воспользуемся вспомогательной функцией [2]

8*

f1(C) =

(16)

в которой значение формпараметра (с), согласно Г Хольштейну и Т. Болену [2], выразим следующим образом:

** 2

С учетом сдвига по дискам относительно средней линии каналов будет W(x) cos e = Wк cos eR . Тогда из уравнения (5) получим

Wx W, cose R

. (8)

c =

dW (x)

(17)

ёх (10,64т + 2,128)

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

Л(Х) = -%ХЬ. (18)

Откуда

При переходе от круговой к прямой решетке, с учетом выражения (5), получим уравнение средней линии канала РК в виде

где

vpW (x)

5*1 ’

Л1 f 37 Л Л2

(19)

f2(c) = l 2+ -jl—- —-

тогда выражение (9) будет выглядеть следующим образом:

(2О)

6 Д315 945 9О72;

С учетом постоянного значения гидродинамической степени диффузорности РК ( WK = W2), согласно

** jw і—

К. Польгаузену в работе [2], 5z1J= yCo = О,278

уравне-

ние (19) запишем

; =

t> z = f2(l)pvW^\~ =

vc

Wк f2(c)pvWK

5z1

или, с учетом значения f2 (c) = 0,235, имеем 0,235pWк v

(21)

(22)

где значение 5*1 определяется по уравнению (14).

Решение уравнений пограничного слоя К. Польгаузе-ном для случая продольного обтекания плоской пластины (т = 0) дает выражение

** V/

5*1 = 0,686 I —^.

УК

Тогда уравнение (21) преобразуем к виду

х21 = 0,343рп0,5^к1,5 или с учетом сдвигового течения

т-1

(23)

(24)

tz1 = 0,343pl 2 yjWl3v(5; +1) cos єR

f 12 ln R2 - ln R1

, Л R2 /

V lnpln 7R 0

V /1 0

(25)

С учетом уравнений (14) и (16) для приведенных выше параметров РК МН при угловой скорости ю = 800 рад/с для ряда реальных рабочих тел число Рейнольдса ^екр5;; изменялось от 34,2 (п = 10-4 м2/с) до 342 (V = 0,7-10-6 м2/с), что существенно ниже его критического значения (рис. 3).

If'

Ілр- tof

Л*

Результаты расчета параметров пограничного слоя 5*, 5** вдоль ограничивающих дисков каналов РК в зависимости от скорости набегающего потока W(х), с учетом сдвига в радиальной плоскости на угол е со следующими параметрами: Яг = 8,25 мм; Я2 = 20,25 мм; шаг решетки ?пр = 14,95 мм (г = 6); /пр = 18,6 мм; ширина канала й1 = 4,3 мм и угловой скорости ю = 800 рад/с (рис. 3); для ряда рабочих жидкостей: вода (п = 1.10-6 м2/с; р = 103 кг/м3); водный раствор глицерина (п = 1. 10-4 м2/с; р = 1,23. 103 кг/м3); ЛЗ-ТК-2 при t = +20 °С (V = 0,7. 10-6 м2/с; р = 691 кг/м3) и при t = -100 °С (V = 1,55.10-5 м2/с; р = 784 кг/м3).

Из графиков видна существенная зависимость роста 5** и 5* с увеличением вязкости рабочей среды. Это приводит к оттеснению области течения, которое можно рассматривать как потенциальное, от стенок канала и сужению проходного сечения межлопаточных каналов на выходе РК до 50 %, что необходимо учитывать при расчете МН. Следует отметить, что значения толщины вытеснения и потери импульса возрастают, достигая на половине расчетного участка 80 % своего значения, соответствующего параметрам на выходе из РК.

Рис. 3. Результаты расчета параметров пограничного слоя для ряда рабочих жидкостей вдоль ограничивающих дисков канала РК (Я1 = 8, 25 мм; Я2 = 20,25 мм; Ь1 = 4,3 мм; ґпр = 14,95 мм (г = 6); Іпр = 18,6 мм; ю = 800 рад/с)

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

Библиографические ссылки

1. Отчет по НИР Б 1.6.08 «Математическое моделирование и гидродинамический анализ трехмерного пространственного течения (I этап). N° 0120.0804367 ; Сиб. гос. аэрокосмич. ун-т. Красноярск, 2008.

2. Шлихтинг Г. Теория пограничного слоя. М. : Наука, 1969.

3. Степанов Г. Ю. Гидродинамика решеток гидромашин. М. : Физматгиз, 1962.

E. M. Kraeva

THEORETICAL ANALYSIS OF FLOW OF CASCADE OF LOW-FLOW CENTRIFUGAL PUMP BY VISCID FLUID STREAM

In view of design features of air-gas channel of impeller oflow-flow centrifugal pump the analysis offlow of cascade by viscid fluid flow is executed. Calculation ofparameters of a boundary layer is executed.

Keywords: hydrodynamics, low-flow centrifugal pump, impeller, boundary layer, diffuser degree.

©Краева Е. М., 2010

i Надоели баннеры? Вы всегда можете отключить рекламу.