Моделирование оптимальных режимов течения жидкой пленки в контакте с газовым потоком

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

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

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

С поверхностью стекания связана система координат OXY, направление оси ОХ совпадает с направлением стекания пленки, ось OY направлена по нормали в пленки жидкости.

Пленка обтекается газовым потоком, который движется параллельно оси ОХ и создает на свободной поверхности жидкости постоянное касательное напряжение т. Если направление движения газового потока совпадает с направлением стекания жидкой пленки, реализуется режим прямотока (г> 0), если направления движения пленки и газа противоположны - режим противотока (г < 0), если газовый поток отсутствует (г = 0) - режим свободного стекания.

Уравнение эволюции свободной поверхности жидкой пленки [1, 2] в линейном приближении имеет вид

5             Э*хр Э²хр _кахр

(О

Ь_хх --1 Н^ + ^д— v +   v + —*-=0

,nax ) at 4 эх* эх2 1 ax где хр(х, t) - отклонение свободной поверхности жидкой пленки от стационарного состояния, bxx=^Ke2F, b4 = -^|- , Ь^ = -^Re3F^F-т^, bx=-Re(F-TY Re - число Рейнольдса,

F = — + 1,5т - число Фруда, о- поверхностное натяжение.

Re

Из уравнения (1) получим дисперсионное соотношение оДЬ_ххк + \^л-Ь₄к* -b^k² + b_xik = 0,                              (2)

где со = co_r + ico,, co_r - частота, o_t - инкремент.

Из соотношения (2) в явном виде найдем:

(-b4k4 +b2k2Ybxxk-bxk

®r = "-----:—j ------->                        (3)

= b4k^(b2+bxbxx^

Результаты вычислительного эксперимента [2] для жидкой пленки воды для чисел Рейнольдса 7?ее[1; 10] приведены на рис. 1-4. Режим прямотока (т< 0) способствует уменьшению скоростей роста возмущений ®_! (рис. 1, кривая 1), режим противотока (т> 0) - увеличению скоростей роста возмущений <у₍ (рис. 1, кривая 3).

В режиме противотока фазовые скорости C-o>_rjk (рис. 2, кривая 3) возмущений меньше фазовых скоростей в режиме свободного стекания (рис. 2, кривая 1).

Рис. 1. Инкремент жидкой пленки (Re = 10): 1-т = -0,15; 2- т = 0; 3-т= 0,15

Рис. 2. Фазовые скорости (Re=10): Т — т = —0,15; 2-т=0; 3-т=0,15

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

Из условия нейтральной устойчивости

®,=0                                    (5)

определим волновое число нейтрально устойчивых возмущений к2 = .№ = Л Re2 f(F -т).                  (6)

DO"

В плоскости (7?е, к) выражение (6) задает кривую (рис. 5), которая отделяет область неустойчивости (щ > 0) от области устойчивости (со, < 0).

Течение жидкой пленки неустойчиво в рассматриваемых диапазонах изменения числа Рейнольдса 7?е. Режим противотока (т > 0) оказывает дестабилизирующее влияние на течение жидкой пленки, область неустойчивости расширяется по диапазону волновых чисел к. Режим прямотока (т<0) способствует стабилизации течения, наблюдается сокращение области неустойчивости по диапазону волновых чисел к. Воздействие газового потока позволяет расширить диапазон неустойчивых возмущений на 48 % в режиме противотока, в режиме прямотока - уменьшить на 57 %.

Режимы течения жидких пленок, для которых инкремент ®, достигает максимального значения (волны максимального роста), называют оптимальными. Параметры потока в целом близки к параметрам волн максимального роста [4].

Рис. 3. Инкремент жидкой пленки (Re = 10): 1-Дет = -10 %; 2-До=0; 3-Да=10%

Рис. 4. Фазовые скорости (Re = 10): 1 - Д<т = -10 %; 2-Дсг = 0; 3 - Д<т = 10 %

Прокудина Л. А., Саламатов Е.А.

Моделирование оптимальных режимов течения жидкой пленки в контакте с газовым потоком

Параметры волн максимального роста найдем из уравнения dk решив которое, получим

,2_ V         64          1 + V1 + (^И^Л?)

к2             ”         к2         ‘

"11                            "11

В плоскости (Re, к) выражение (8) задает кривую максимального роста (рис. 5).

Рис. 5. Области неустойчивости жидкой пленки:

1 - режим противотока (т= 0,15); 2-т= 0; 3- режим прямотока (т= -0,15); сплошные линии - кривые нейтральной устойчивости;

пунктирные линии - кривые максимального роста

Из зависимости(4) найдем

^ = -lbuk—. dk 11 dk

Для волн максимального инкремента (7) получим

^ = 0. dk

(Ю)

т.е. волны максимального роста обладают минимальной фазовой скоростью (рис. 2). Для групповой скорости

К =       (£С) = £— dk dkv ’ dk

Учитывая соотношение (9), получим равенство фазовой скорости С и групповой скорости V волн максимального роста (рис. 6)

7<к_т^С(к_тУ                       (И)

Увеличение |г| в режиме прямотока (т< 0) способствует уменьшению фазовых С и групповых V скоростей. При увеличении |г| в режиме противотока (г> 0) фазовые С и групповые V скорости увеличиваются, если Re < 8 (рис. 6, кривые 1-4). Для Re > 8 увеличение |т| в режиме противотока приводит к уменьшению фазовых и групповых скоростей возмущений, обладающих максимальной скоростью роста (рис. 6, кривые 5, 6).

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

Рис. 6. Фазовые С и групповые V скорости волн максимального роста: 1-Re = 5; 2-Re = 6; 3-Re = 7;

4 - Re = 8; 5 - Re = 9; 6 - Re = 10

Рис. 7. Инкременты ед волн максимального роста: 1 - Re = 5; 2-Re = 6; 3-Re= 7; 4-Re = 8; 5-Re = 9; 6-Re =10

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