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

Контроль и испытания

УДК 621.22: 621.517: 532.528

ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ РАБОЧЕГО ПРОЦЕССА КАВИТАЦИОННОГО СМЕСИТЕЛЯ

Е,К. Спиридонов, А.А. Пантюхин

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

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

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

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

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

Рис. 1. Принципиальная схема струйного насоса

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

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

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

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

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

Для исследования рабочего процесса струйного кавитационного смесителя, конструкция и основные размеры проточной части которого представлены на рис. 2, в лаборатории кафедры «Гидравлика и гидропневмосистемы» Южно-Уральского государственного университета была проведена серия экспериментов.

Многоструйное сопло

Дроссельные шайбы

с/отв^О ММ Д?ол1<г=0,5 ММ

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

В состав испытательного гидравлического стенда входят: напорный бак, электронасосный агрегат ЦНСА 38-110-01, кавитационный смеситель, дроссельный расходомер и контрольноизмерительные приборы (рис. 3).

Рис. 3. Принципиальная схема испытательного стенда

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

При проведении экспериментальных исследований изменение подачи насоса осуществляется с помощью регулирующей задвижки 31; изменение противодавления на выходе исследуемого гидравлического устройства - с помощью регулирующей задвижки 33. Перепад давления в устройстве контролируется по показаниям манометров МНЗ и МН4 . Расход жидкости определяется величиной перепада давления на дроссельном расходомере, контролируемого по показаниям манометра МН1 и вакуумметра ВМ2 .

При экспериментальном исследовании характера рабочего процесса струйного кавитационного смесителя в качестве активной среды выступала вода, а эжектируемой среды - воздух. Конструкция кавитационного смесителя позволяет исключать полностью либо изменять расход эжектируемого потока с помощью набора сменных дроссельных шайб (см. рис. 2). Проводились две серии экспериментов: первая серия с нулевым коэффициентом эжекции - исключался подсос эжектируемой среды (входной патрубок конфузора герметично закрывался заглушкой); вторая серия с малыми коэффициентами эжекции - при наличии подсоса эжектируемой среды (во входной патрубок конфузора устанавливалась дроссельная шайба с диаметром центрального отверстия с1отв = 0,5 мм). Массовый расход воздуха в среднем составлял т = 4,84-10-5 кг/с, а объемный коэффициент эжекции изменялся в диапазоне а = 0,084...0,094.

В каждой серии испытаний снималась дроссельная характеристика аппарата, выражающая зависимость давления всасывания Р2 от противодавления Р5 при том или ином постоянном значении давления питания Р1. Методика проведения испытаний заключалась в следующем: вначале устанавливалось с помощью регулирующей задвижки 31 заданное значение давления Рх перед кавитационным смесителем, контролируемое по показаниям манометра МНЗ; перекрытием регулирующей задвижки 33 осуществлялось постепенное увеличение противодавления Р5, контролируемого по манометру МН4 , установленному на выходе кавитационного смесителя; одновременно с помощью вакуумметра ВМЗ измерялось давление Р2 в приемной камере. При увеличении противодавления Р5 до некоторого критического значения происходил срыв потока и затопление приемной камеры с резким увеличением давления Р,. Некоторые результаты экспериментальных исследований представлены в табл. 1, 2 и на графике (рис. 4).

Таблица 1

Результаты первой серии опытов (с нулевым

Таблица 2

Результаты второй серии опытов (при малых

/^■ІО5, Па /-Гс-105,Па ^■10?, Па е. м3/с

11 0,09 1,75 0,007

0,09 2,48 0,007

0,09 3,27 0,007

0,08 4,74 0,007

0,08 5,19 0,007

Срыв потока - 5,23

9 0,09 2,48 0,006

0,09 3,46 0,006

0,09 4,25 0,006

0,09 4,35 0,006

0,10 4,44 0,006

Срыв потока - 4,49

7 0,09 2,24 0,005

0,09 2,83 0,005

0,10 3,56 0,005

0,10 3,41 0,005

0,10 3,54 0,005

Срыв потока - 3,56

^а&-105,Па Р2вбе-105,Па Р"г,с • 105, Па <2, м3/с

11 0,07 3,98 0,007

0,06 4,57 0,007

0,07 5,13 0,007

0,07 5,18 0,007

0,08 5,23 0,007

Срыв потока - 5,23

9 0,07 2,97 0,006

0,07 3,66 0,006

0,07 3,95 0,006

0,07 4,16 0,006

0,08 4,35 0,006

Срыв потока - 4,39

7 0,08 1,99 0,005

0,08 2,92 0,005

0,08 3,37 0,005

0,09 3,41 0,005

0,10 3,50 0,005

Срыв потока - 3,50

Р2 • 10, Па 6,0 т

5.0

4.0 +

3.0

2.0 +

1,0-

0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 5,5

Р5обс-10, Па

/>а&=7-10Па ?“&=9-105Па ^Р^МнОПа

Рис. 4. Зависимость давления всасывания от противодавления при нескольких значениях давления питания

Данные, полученные в серии экспериментов с нулевым коэффициентом эжекции, представлены на графике светлыми точками, а при малых коэффициентах эжекции - темными. Из представленных на графике рис. 4 и в табл. 1 и 2 результатов опытов следует, что при постоянном давлении Рх на входе устройства и изменении противодавления Р5 параметры течения парогазожидкостной смеси: объемный расход <2 и давление Р2 за срезом сопла не изменяются. Это является подтверждением наличия режима «запирания» в рабочей камере кавитационного смесителя, который формируется при сверхзвуковом парожидкостном течении до прыжка перемешивания. Изменение противодавления до момента срыва потока и затопления приемной камеры определяет только положение прыжка перемешивания в рабочей камере, что подтверждает предположение о том, что переход сверхзвукового течения в дозвуковое с одновременной конденсацией паров происходит в прыжке перемешивания, за которым образуется течение капельной жидкости. При повышении давления прыжок перемешивания смещается против течения к началу рабочей камеры; в момент срыва потока происходит резкое повышение давления в приемной камере и течение в струйном смесителе становится однофазным течением капельной жидкости.

Базовое уравнение расчета кавитационного смесителя, полученное в процессе математического анализа рабочего процесса [2] имеет вид:

р\ ~Р5

Р\~Рнп

0)

1 + 0 + СООП

Здесь Рнп- давление насыщенных паров жидкости; £соп, £диф - коэффициенты сопротивления многоструйного сопла и диффузора. В работе [2] приведены предполагаемые значения коэффициентов сопротивления £со„ = 0,1; ^диф = 0,25; СрК ~ коэффициент потерь на трение по длине рабочей камеры; €1 - относительная площадь сопла; о - число кавитации.

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

Выразим Р5 из уравнения (1) для того, чтобы сопоставить результаты расчета оптимального противодавления с опытными значениями противодавления срыва потока (табл. 1,2):

(** \ Д; + 0 ■- п); - Т" ] (2)

1 + а + Ъсоп

Для расчета оптимальных значений Р5 необходимо определить фактические значения относительной площади сопла, коэффициентов сопротивления, а также числа кавитации.

Коэффициент потерь на трение по длине рабочей камеры [3]

С=Л-г- = °.57>

>рк

рк

где Л - коэффициент гидравлического трения; 1рк - длина рабочей камеры; <Зрк - диаметр проходного сечения рабочей камеры (см. рис. 2)

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

/ \

Л : 0,0255,

2 = 0,11-

где А - эквивалентная абсолютная шероховатость, А = 0,75 мм .

Коэффициент сопротивления диффузора, согласно Идельчику [4], равен

ідиф

8-8Іп(/?/2)

Г £„„1

1 Рк + віп р -

1 *5 ] [_ 1 Ъ ]

= 0,34.

Здесь £5 - площадь проходного сечения проточной части в сечении 5-5, £5 =5^; Р - угол конусности диффузора, р = Т.

Значения оставшихся членов уравнения (2) определяются не только геометрией проточной части устройства, но и параметрами потока, так как в связи со сжатием струи при истечении из круглого отверстия площадь проходного сечения сопел и действительная площадь струй различны. Степень сжатия струи оценивается коэффициентом сжатия £:

р — с /с ° и струи °соп 9

где ^ - действительная площадь поперечного сечения струи; 5сои - площадь проходного се-

чения многоструйного сопла.

Основываясь на опытных данных, определим значения коэффициента сопротивления соплового устройства, относительной площади сопла и числа кавитации исследуемого кавитационного смесителя, расчетная схема которого представлена на рис. 5.

1. 2,

Рис. 5. Расчетная схема кавитационного смесителя

В соответствии с уравнением Бернулли потери полного давления, вызванные внезапным сужением проточной части, при движении жидкости между сечениями 1 и 2:

Щ-2 = Ссоп ^уи = РХ ~ Р2 +^1 ~ К1руи ).

где V) - скорость течения жидкости в сечении 1; V,

струи

стии.

Выразим коэффициент сопротивления сопла в виде

Ґ А

Ссоп

(р1-рг)

+

V2 у струи

-1

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

(3)

Р ^струи

Согласно уравнению (3) для определения коэффициента сопротивления сопла необходимо определить действительную скорость течения струи Уструи . Для этого запишем уравнение коли-

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

V = Ул + -

струи 1

\Рх-Рг)

Р-6

(4)

Зная значение

■ струи» можно рассчитать действительную площадь сечения струи и определить, в первом приближении, ее коэффициент сжатия

е = 0/1рсоп -уструи)- (5)

Подставляя в уравнения (3) и (5) значения действительной скорости течения струи, рассчитаем значения коэффициентов сопротивления сопла и сжатия струи и представим результат на графике (рис. 6).

Р"&=7*Ю^1а

♦ Р1й&’=9*105Па

*Р,абМ1-1071а

Рис. 6. Зависимость коэффициента сопротивления сопла от коэффициента сжатия струй при нескольких значениях давления питания

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

Рх =1 МО5 Па до Р] = 7• 105 Па величина коэффициента сопротивления снижается в диапазоне от Ссоп = 0;50 до С,соп - 0,41, и это в значительной степени отличается от значения, приведенного в работе [2], где предполагается, что £соп * 0,1 носит универсальный характер. В действительности коэффициент соплового устройства определяется как геометрическими параметрами проточной части, так и динамическими характеристиками потока.

Относительная площадь сопла

П =

I*.

струи

^рк ^ рк

где - площадь проходного сечения рабочей камеры (см. рис. 2).

Число кавитации <т в зависимости от изменения параметров течения можно определить по следующему уравнению:

Л-Р„

СГ =

струи

Зная действительную скорость течения струй и величину давления насыщенных паров воды в диапазоне изменения температур Т = 25...30 °С, РИП «2500 Па, рассчитаем число кавитации,

подставляя в уравнение осредненные значения давления Р2 из табл. 1. Результаты расчетов представим на графике (рис. 7).

/ГЧО' Па

Рис. 7. Зависимость числа кавитации от давления питания

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

В работе [2] приводится следующая формула расчета числа кавитации: а = 0,07 +1,36 • О ■ (1 - □), (6)

согласно которой в исследованном диапазоне давлений изменение расчетных значений составляет 0,367 < сг < 0,378, что значительно больше фактических величин числа кавитации (см. рис. 7). Из чего следует, что при расчете по формуле (6) недооценивается степень развития кавитационных процессов в исследуемом устройстве.

Рассчитаем оптимальные значения противодавления Р5, подставляя в уравнение (3) фактические значения (табл. 3) относительной площади, коэффициента сопротивления сопла, а также числа кавитации. Результаты расчетов и экспериментальные данные представим на графике.

Таблица 3 Сводная таблица результатов расчета

105,Па О Ссоп а Р?бс-10?, Па

11 0,35 0,50 0,008 3,83

9 0,33 0,44 0,011 3,18

7 0,32 0,41 0,015 2,50

РГЧ05, Па

Рис. 8. Зависимость величины оптимального противодавления и противодавления срыва потока от давления питания

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

Спиридонов Е.К., Пантюхин А.А. Экспериментальные исследования рабочего процесса

________________________________________________________________кавитационного смесителя

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

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

Данная работа выполнена при финансовой поддержке Правительства Челябинской области.

Литература

1. Спиридонов, Е.К. Энергетический анализ жидкостногазовых течений / Е.К. Спиридонов //Вестник ЮУрГУ. Сер. «Машиностроение». - 2003. - Вып. 3.№ 1 (17). - С. 141-151.

2. Спиридонов, Е.К Рабочий процесс и характеристики гидродинамических кавитационных эмульгаторов /Е.К. Спиридонов, Л.С. Прохасъко //Вестник ПГТУ. Серия «Аэрокосмическая техника». - Пермь, 2000. —№ 5. - С. 93-98.

3. Прандтль, Л. Гидроаэромеханика / Л. Прандтлъ. - М. : Изд-во иностранной литературы, 1951.-575с.

4. Идельчик, И.Е. Справочник по гидравлическим сопротивлениям / И.Е. Иделъчж. - М.: Машиностроение, 1992.