Модели расчета турбулентного пограничного слоя у проницаемой поверхности

39

Размещено на http://www.allbest.ru/

УДК 517.925.8;517.928

Модели расчета турбулентного пограничного слоя у проницаемой поверхности

Лаптева Елена Анатольевна

к.т.н., старший преподаватель

Саитбаталов Марат Викторович

cоискатель

ФГБОУ ВПО «Казанский государственный энергетический

университет», Казань, Россия

Фарахов Тимур Мансурович

к.т.н.

инженер, ООО Инженерно Внедренческий Центр «Инжехим»,

Казань, Россия

Аннотации

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

Ключевые слова: КОНДЕНСАЦИЯ, ИСПАРЕНИЕ, ТЕПЛООТДАЧА, МАССООБМЕН, НАСАДКА, ИМПУЛЬСООТДАЧА

UDС 517.925.8;517.928

THE MODELS OF CALCULATION OF TURBULENT BOUNDARY LAYER AT PERMEABLE SURFACE

Lapteva Elena Anatolevna

Cand.Tech.Sci., senior lecturer

Saidbatalov Marat Viktorovich

competitor for degree

Kazan State Power Engineering University (KSPEU), Kazan, Russia

Farakhov Timur Mansurovich

Cand.Tech.Sci., engineer

LLC Engineering-Promotional Center “Inzhekhim”, Kazan, Russia

In the article we have considered the mathematical models of turbulent boundary layer corresponding with condensation and evaporation. The expressions for accounting that effect have been derived. We have also compared them with the known experimental data

Keywords: CONDENSATION, EVAPORATION, HEAT TRANSFER, MASS EXCHANGE, NOZZLES, IMPULSE TRANSFER

Введение

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

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

Перенос импульса

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

Уравнение для коэффициента переноса импульса выводится из потоковых соотношений и в общем виде [1]:

(1)

Уравнение (1) является интегральным выражением потерь импульса потока в пограничном слое посредством трения. С другой стороны величина касательного напряжения трения входит в известное уравнение Рейнольдса:

(2)

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

Таким образом, решение уравнения (2) при 0 асимптотически сводится к решению уравнения:

(3)

Интегрируя уравнение (3) по y от нуля до текущей координаты y в рамках рассматриваемой модели безградиентного (или квазибезградиентного) течения имеет:

(4)

где .

Решение уравнения (4), очевидно, является асимптотическим и действительным лишь в области . Несмотря на кажущуюся грубость такого приближения оно с успехом используется в алгебраических моделях турбулентности при распространении их на случаи обтекания проницаемых поверхностей: [4-10].

В работе [11] была предложена форма аппроксимации распределения касательного напряжения по высоте пограничного слоя у проницаемой поверхности на основе (4) в виде:

(5)

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

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

(6)

Коэффициент пропорциональности в данном случае полагается неизменным ввиду консервативности структуры турбулентного пограничного слоя. Впрочем, для градиентных турбулентных пограничных слоёв у проницаемой поверхности многими авторами [5-7] отмечается некоторое повышение коэффициента . В первых двух приведенных работах для случая обтекания гладкой поверхности с положительным перепадом давления рекомендуется . В работе [7] значение определялось исходя из соотношений двузонного пограничного слоя и интегральных характеристик гидродинамических условий промышленного аппарата.

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

(7)

Уравнение (5) с учетом (7) может быть переписано в виде:

(8)

Здесь индекс «0» будет употребляться как к внешнему потоку (вне пограничного слоя) и внешнему потоку при стандартных условиях, так как полагается, что влияние вдува ограничивается лишь областью пограничного слоя.

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

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

(9)

Учитывая успешный опыт авторов [5,6,8] можно полагать, что форма демпфирующего множителя в виде , показавшая наилучшие результаты для простых случаев турбулентного течения может быть применена и для описания обтекания проницаемых поверхностей. В рамках соотношений Таунсенда демпфирующий фактор как функция от запишется в следующем виде:

(10)

Для случая обтекания полупроницаемой пластины турбулентным потоком, исходя из вышеописанного, рекомендуются значения коэффициентов в формуле (10) близкие к турбулентному невозмущенному течению: , (последнее значение коэффициента A предлагается авторами [6] в функции турбулентной вязкости для течений с совокупным воздействием трех возмущающих факторов: поперечной кривизны обтекаемой поверхности, градиента давления, и вдува).

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

Для проверки этой модели были выбраны экспериментальные данные полученные в работах [10,11] как наиболее достоверные. Сравнивались экспериментальные и вычисленные согласно приведенным выше формулам профили скорости по высоте пограничного слоя. В указанных источниках результаты экспериментов были обработаны в безразмерных координатах

, .

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

Уравнение (1), переписанное в дифференциальном виде:

(11)

В безразмерных координатах и уравнение (11) выразится в следующем виде:

(12)

Если использовать в (12) функцию относительного трения , предложенную авторами [4], то (12) выразится как:

(13)

Функция относительной турбулентной кинематической вязкости (9) и демпфирующего множителя (10) в этих координатах предстанут соответственно в виде:

(14)

(15)

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

(16)

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

В качестве функции коэффициента трения при стандартных условиях была использована эмпирическая формула из [6]:

(17)

Единственным неопределенным параметром в системе уравнений (13)-(17) является функция относительного трения . В работе [9] была предложена аппроксимационная форма закона относительного трения в виде:

(18)

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

В работе [4] были предложены две формулы для нахождения критического значения параметра b при заданном числе Рейнольдса:

(19)

(20)

В формуле (20) - параметр неизотермичности, в рассматриваемой задаче .

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

Рисунок 1. Графики функций критических значений параметра вдува. от числа соответственно a-(19), b -(20).

Как видно на рисунке 1 наблюдается довольно большое расхождение между этими двумя законами связи и .

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

Очевидно, что для точки отрыва пограничного слоя должны соблюдаться условия:

(21)

Относительная толщина пограничного слоя находится из решения уравнения:

(22)

Таким образом, нахождение значения для каждого заданного сводится к краевой задаче для обыкновенного дифференциального уравнения (13) и входящих в его состав функций (14)-(17) при условиях (21) и (22).

Задача решалась численными методами в свободно распространяемом пакете Scilab. Поскольку дифференциальное уравнение не является жестким, для его решения был применен стандартный алгоритм Рунге-Кутты- Фельдберга 4-5-го порядка точности, для численного интегрирования использован алгоритм Гаусса-Кронрода. Общий алгоритм вычислений организован следующим образом: задавалось число Рейнольдса , первоначальные значения и ; интегрирование дифференциального уравнения проводилось от верхнего предела (начальные условия - ), в случае, если получаемый профиль скорости не удовлетворял условию (22) требуемое значение искалось путем градиентного приближения до достижения условия . Аналогично, если при соблюдении предыдущего условия при заданном не соблюдалось условие, подбиралось новое значение до достижения условия (пересчет производился для каждого нового ). Были выбраны следующие допустимые величины невязки , .

Рисунок 2. Расчетный предотрывный профиль скорости, полученный по уравнениям (13), (21)

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

Результаты расчетов аппроксимируются формулой:

(23)

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

Рисунок 3. График результатов численного исследования краевой задачи (13), (21) (точки) и аппроксимирующей их зависимости (23).

С помощью формулы (23) система уравнений (13)-(18) оказывается замкнутой. Для сравнения в расчет включены выражения для безразмерного критического параметра вдува (19) и (20). Расчет производился для каждой пары и , указанных работах [10,11]. Начальным условием для (13) было . Интегрирование производилось до достижения условия (22). Как и в решении, только что рассмотренной краевой задачи для интегрирования обыкновенного дифференциального уравнения использовался алгоритм Рунге-Кутты Фельдберга 4-5-го порядка точности. Численное интегрирование получавшегося профиля скорости (после достижения указанного в экспериментальных данных верхнего уровня) проводилось по двадцатиодноточечному алгоритму Гаусса-Кронрода.

На рисунке 4-8 показаны результаты расчетов профилей скорости по высоте пограничного слоя в безразмерном виде, вычисленные согласно уравнениям (13)-(20), (22), (23), а также экспериментальные данные источников [10] и [11]. Как, видно, приведенные модели уравнения с использованием функции (23) дает наилучшие результаты для всего рассматриваемого диапазона, даже для больших значений .

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

Рисунок 4. Графики распределения относительной скорости по высоте турбулентного пограничного слоя, отнесённой к толщине потери импульса при параметре вдува и числе Рейнольдса: a - экспериментальные данные [11]; b - расчет с использованием функции (19); c - расчет с использованием функции (20); d - расчет с использованием функции (23).

Рисунок 5. Графики распределения относительной скорости по высоте турбулентного пограничного слоя, отнесённой к толщине потери импульса при параметре вдува и числе Рейнольдса : a - экспе-риментальные данные [11]; b - расчет с использованием функции (19); c - расчет с использованием функции (20); d - расчет с использованием функции (23)

Рисунок 6. Графики распределения относительной скорости по высоте турбулентного пограничного слоя, отнесённой к толщине потери импульса при параметре вдува и числе Рейнольдса: a - экспериментальные данные [10]; b - расчет с использованием функции (19); c - расчет с использованием функции (20); d - расчет с использованием функции (23).

Рисунок 7. Графики распределения относительной скорости по высоте турбулентного пограничного слоя, отнесённой к толщине потери импульса при параметре вдува и числе Рейнольдса: a - экспериментальные данные [10]; b - расчет с использованием функции (19); c - расчет с использованием функции (20); d - расчет с использованием функции (23)

Рисунок 8. Графики распределения относительной скорости по высоте турбулентного пограничного слоя, отнесённой к толщине потери импульса при параметре вдува и числе Рейнольдса: a - экспериментальные данные [11]; b - расчет с использованием функции (19); c - расчет с использованием функции (20); d - расчет с использованием функции (23)

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

Перенос тепла и массы компонента

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

Так же как при выводе уравнения (4) (в окрестности ) запишем асимптотическое решение уравнения сохранения энтальпии:

(24)

где .

Функция распределения относительной плотности теплового потока по высоте пограничного слоя, от координаты :

(25)

Величина имеет ту же размерность что и i является аналогом динамической скорости для теплообмена. Очевидно также, что при тепломассообмене с постоянной плотностью член .

Коэффициент турбулентной температуропроводности в случае необходимости может быть выражен в форме предложенной Себеси [12] и Мейером и Роттой [13] с учетом закона стенки Таунсенда (6):

(26)

где демпфирующий коэффициент определяется из выражения:

(27)

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

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

(28)

где

функция вычисляется согласно (25).

Для, собственно, массопереноса в пограничном слое балансовое дифференциальное уравнение (при условии размерности концентрации в долях массы - (кг/кг) запишется аналогично (2):

(29)

Откуда массовый поток в окрестности аналогично (4) и (24) можно записать: турбулентный импульс перенос теплота

(30)

После преобразования (30) в безразмерный вид и с сохранением формы выражения относительной плотности потока аналогично (28) и (25) имеем:

(31)

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

Так, например, при использовании в качестве функции относительного потока импульса в невозмущенном турбулентном пограничном слое полинома третьей степени имеем

, ,

где равно толщине диффузионного пограничного слоя для рассматриваемой задачи (течение у полупроницаемой пластины).

И соответственно коэффициент турбулентной диффузии может быть выражен аналогично (9):

(32)

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

Наконец, интегральное уравнение для определения коэффициента массоотдачи в турбулентном пограничном слое:

(33)

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

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

Окончательно система уравнений для расчета характеристик тепломассообмена в пограничном слое запишется в следующем виде:

(34)

(35)

(36)

Порядок решения может быть построен следующим образом:

1) задается массив первоначального распределения безразмерных величин и по высоте пограничного слоя, вычисляются, в первом приближении, относительные плотности переноса тепла и массы: , ;

2) из уравнения (34) находится распределение безразмерной разности скоростей по высоте пограничного слоя - и толщина эквивалентного пограничного слоя - ;

3) итерационно решается уравнение (35) до достижения сходимости по параметру , в результате вычисляются массив распределения безразмерной концентрации - и относительный параметр вдува (отсоса) - ;

4) проводится итерационная процедура из пунктов 2 и 3 до достижения сходимости по величине , в результате уточняются профили и , а также безразмерный параметр вдува ;

5) решается уравнение (36), в результате чего определяется профиль безразмерной энтальпии по высоте пограничного слоя - и относительной плотности потока теплоты через поверхность раздела - .

6) В случае если учитывается неизотермичность по высоте пограничного слоя, расчет проводится повторно с пункта 2 до уменьшения невязки по параметру до приемлемого уровня.

Первоначальное приближение в пункте 1 не обязательно, и призвано лишь, ускорить сходимость итерационной процедуры. Рекомендуется в качестве первого приближения расчетных профилей использовать форму гидродинамической аналогии предложенной в [14]. Следует все же отметить, что полное сведение задачи расчета переноса в пограничном слое к интегрально-параметрическому виду нежелательно виду неизбежного ограничения области применения полученных таким образом формул. В качестве первоначального приближения коэффициента массоотдачи рекомендуется использовать классическую форму гидродинамической аналогии Чильторна и Кольборна [1,2].

Для решения системы уравнений (34)-(36) использовался алгоритм Рунге-Кутты Фельдберга 4-5-го порядка точности.

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

Учет возмущающих факторов в турбулентном пограничном слое

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

Используемая здесь функция турбулентной вязкости записана в следующем виде:

(37)

При использовании ее в рамках предлагаемого метода должно выполняться условие:

(38)

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

(39)

где - скорость турбулентных пульсаций.

Разложив функцию (37) в ряд Тейлора по переменной в окрестности нуля, получим:

Таким образом, с точностью до члена порядка степень n в уравнениях (38) и (37) может быть определена из приближенного равенства:

Как показывают расчеты, основное влияние на степень m в формулах (38) и (39) при использовании (37) оказывает степень n. Коэффициент A в большей степени влияет на значение суммарной толщины вязкого и переходного подслоя,

:

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

Нахождение второго коэффициента демпфирующего фактора в (37) возможно из условия консервативности длины пути смешения к возмущающим воздействиям показанному Кутателадзе С.С. [4]. Иными словами, в расчетах предлагается руководствоваться гипотезой, что все действие возмущающих факторов ограничивается изменением толщины вязкого подслоя. В алгебраических моделях турбулентности такие условия реализуются путем изменения коэффициента B [1,4] в профиле скорости в логарифмической области:

(40)

Очевидно, что

(41)

В приближении эквивалентного внутреннего (логарифмического) подслоя можно также записать:

(42)

и для течения осложненного возмущающими факторами:

(43)

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

(44)

Тогда для любой безразмерной координаты , отвечающей требованиям , справедливо равенство:

(45)

Уравнение (45) относительно является нелинейным и может быть решено только численно.

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

(46)

Откуда, задаваясь погрешностью O можно при известном n вычислить коэффициент A.

Проверка на модели невозмущенного течения показала, что влияние коэффициента O на значение A может практически не приниматься во внимание при .

Проверка по результатам экспериментов

В качестве проверки адекватности модели применительно к описанию массопереноса в газовой фазе при движении совместно со стекающей пленкой жидкости проведено сравнение с критериальными формулами полученных различными авторами в формах:

(47)

(48)

Здесь - число Рейнольдса газа, вычисляемое по эквивалентному диаметру и осредненной по поперечному сечению скорости.

Коэффициенты в формулах (47) и (48) указаны в таблице 1.

Также для сравнения использовалась формула, полученная авторами [1,7] на основании - модели диффузионного пограничного слоя:

(49)

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

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

Толщина пленки воды с учетом поперечной кривизны канала рассчитывалась по уравнению, полученному в [27]:

(50), где

; .

Коэффициент сопротивления для газовой фазы вычислялся по зависимости [26]:

(51)

Число Рейнольдса пленочного течения:

(52)

На рисунке 9 и 10 показаны графики зависимостей числа Шервуда в газовой фазе от числа Рейнольдса на основе средней скорости в свободном сечении канала - при числах Шмидта и соответственно. Расчеты осуществлялись при использовании в критериальных формулах числа Рейнольдса по абсолютной скорости газа - . Для уравнений, основанных на гидродинамической аналогии, расчет диффузионного числа Стантона также осуществлялся на основе . Такой метод анализа зависимостей по массообмену, показанный в [26] позволяет унифицировать данные, полученные в опытах при противоточном и прямоточном движении фаз.

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

Таблица 1. ДАННЫЕ ПО ФОРМУЛАМ ДЛЯ ФОРМУЛ (47) И (48).

№	Форма уравнения	A	m	n	Литература	Усл. обозн. на рис. 9 и 10
1	(47)	0,023	0,83	0,44	16	a
2	(47)	0,026	0,83	0,33	17	b
3	(47)	0,03	0,8	0,33	18	c
4	(47)		0,75	0,63	19	d
5	(47)	0,0208	0,8	0,44	20	e
6	(47)	0,025	0,8	0,33	12,1	f
7	(47)	0,033	0,77	0,33	22	g
8	(47)	0,023	0,8	0,4	23	h
9	(48)	-	-	0,33	24	i
10	(48)	-	-	0,33	25	j
11	(4)	-	-	0,33	14	k

`Рисунок 9. Сравнительный график зависимостей числа Шервуда от числа Рейнольдса газовой фазы, полученных по формулам различных авторов при ; условные обозначения: a-k - формулы, (48) и (49), см. табл. 1., , и расчет по модели при кг/кг, кг/кг и кг/кг соответственно.

`Рисунок 10.Сравнительный график зависимостей числа Шервуда от числа Рейнольдса газовой фазы, полученных по формулам различных авторов при ; условные обозначения: a-k - формулы (47)-(49), см. табл. 1, , и расчет по модели при кг/кг, кг/кг и кг/кг соответственно.

Как видно на рисунке 9 и рисунке 10, в целом можно говорить о удовлетворительном согласовании полученной модели с данными других авторов во всем интервале практически значимых чисел Шмидта. Наблюдаемый довольно большой разброс результатов используемых формул, по всей видимости, связан с различными условиями проведения экспериментов. Наилучшее совпадение результатов при малых значениях движущих сил переноса приходиться на различные формы гидродинамической аналогии между трением и массоотдачей (уравнения (48) и (49)). Практически важно, хорошее совпадение с результатами авторов [1, 23], чьи формулы часто рекомендуются для практических инженерных расчетов [16,26]. Влияние разности концентраций долгое время остается незначительной и лишь при сравнительно больших значениях сказывается на результатах.

В качестве поверки при моделировании тепломассообмена в насадочном слое проводилось сравнение коэффициентов массоотдачи, полученных экспериментально при испытаний макета вентиляторной градирни с сетчатыми насадками [25].

Экспериментальный макет представлял собой цилиндрическую трубу диаметром 0,211 м, заполненную 18-ю цилиндрическими блоками сетчатой насадки на высоту 0,4 м. Осуществлялось противоточное движение фаз. Вода с помощью форсунки подавалась сверху навстречу потоку охлаждающего воздуха. Скорость воздуха (на полное сечение макета) варьировалась в диапазоне 0,72 - 1,07 м/с; плотность орошения 7,61 м³/ м² •ч; температура воздуха изменялась в диапазоне 26-28 єС; температура воды - 28 - 40 єС.

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

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

Коэффициент сопротивления для потока воздуха вычислен с использованием зависимости (51).

На рисунке 11 показано сравнение значений коэффициента массоотдачи

представленной модели с экспериментальными данными.

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

W_Г , м/с.

Рисунок 11. Сравнительный график расчетов с экспериментом.

Как показали расчеты, сделанные допущения могут применяться лишь при относительно высоких плотностях орошения и нагрузках по газу ниже точки подвисания. Так при понижении нагрузки по воде до 4,93 м³/м²•ч существенно нарушается качественное согласование модели с экспериментом. Это может косвенно свидетельствовать о значительной неравномерности обтекания поверхности насадки при низких плотностях орошения и, как следствие, большом снижении коэффициента активности поверхности насадки.

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

Список использованной литературы

1. Лаптев А.Г. Модели пограничного слоя и расчет тепломассообменных процессов. - Казань: Изд-во Казанск. гос. ун-та, 2007.

2. Александров А.И. Массопередача при ректификации абсорбции многокомпонентной смесей. Л., «Химия»,1975

3. Долгов А. Н., А.Г. Лаптев Математическая модель дегазации в насадочных аппаратах. // Изв. ВУЗов. Проблемы энергетики. №5-6. 2012.

4. Кутателадзе С.С., Леонтьев А.И. Тепломассообмен и трение в турбулентном пограничном слое - 2-е изд., перераб. - М.: Энергоатомиздат, 1985.

5. Патанкар С.В., Сполдинг Д.Б. Тепло- и массообмен в пограничных слоях / пер. Лыков А.В. - М.: Энергия, 1971.

6. Гарбарук А.В. Л.Ю.В., Стрелец М.Х. Турбулентный пограничный слой при одновременном влиянии продольного градиента давления, вдува (отсоса) и поперечной кривизны поверхности // Теплофизика высоких температур. 2002. Том 40. N № 3. с. 436-442.

7. Дьяконов С.Г., Елизаров В.И., Лаптев А.Г. Теоретические основы и моделирование процессов разделения веществ. - Казань: Изд-во Казанского университета, 1993.

8. Гарбарук А.В., Лапин Ю.В., Стрелец М.Х. Простая алгебраическая модель турбулентности для расчета турбулентного пограничного слоя положительным перепадом давления // Теплофизика высоких температур. 1999. Т. 37. № 1. с. 87-91.

9. Squire H.B., Young A.D. The calculation of the profile drag of aerofoils. R. & M. No. 1838, British A.R.C., 1938.

10. Леоньтев А.И., Миронов Б.П., Луговский П.П. Экспериментальное определение критического параметра вдува на пористой пластине // Инженерно-физический журнал. 1966. Том 10. N № 4. с. 447-451.

11. Кутателадзе С.С., Леоньтев А.И. Тепломассообмен и трение в турбулентном пограничном слое. - М.: Энегрия, 1972.

12. Себеси Т. Метод расчета коэффициента турбулентной теплопроводности и турбулентного числа Прандтля // Теплопередача. 1973. Том 95. N № 2.

13. Мейером и Роттой .

14. Лаптев А.Г., Саитбаталов М.В. Аналогия переноса импульса, массы и теплоты в насадочных элементах градирен// Изв. Вузов. Проблемы энергетики,2009,-№1-2,-С.140-144.

15. Gilliland E.R., Sherwood Т.К. Diffusion of Vapors into Air Streams // Ind. Eng. Chem. 1934. v. 26. NN 5. 516-523.\

16. Войнов Н.А., Николаев Н.А. Пленочные трубчатые газожидкостные реакторы (гидродинамика тепло и массообмен). - Казань: "Отечество", 2008.

17. Van_Krevelen D.W., Hoftijzer P.J. Studies of gas absorption // Rec. Trav. Chim. 1949. v. 68. N N 4. p. 221-223.

18. Кузнецов М.Д. Определение коэффициентов скорости абсорбции по методу подобия // Журнал прикладной химии 1948. т. 26. N 1. с. 48-57.

19. Турхан Э.Я., Жиделева К.П. Скорость абсорбции-десорбции паров воды и теплопередача в системе газ - серная кислота // Хим. пром. 1946. 11. N с. 6-9.

20. Barnet W.I., Kobe K.A. Heat and Vapor Transfer in a Wetted-Wall Tower // Ind. Eng. Chem. 1941. v. 33. N 33. p. 436-442.

21. Johnstone H.F., Pigford R.L. Distillation in a Wetted Wall Column. // Trans. Am. Inst. Chem. 1942. v. 38. N 1. p. 25--51.

22. Жаворонков Н.М., Малюсов В.А. Исследование гидродинамики и массопередачи в процессах абсорбции и ректификации при высоких скоростях потоков // Теор. основы хим. технологии. 1967. т. 1. N 5. c. 562-577.

23. Дытнерский Ю.И., Борисов Г.С. Исследование массообмена в жидкой фазе // Процессы хим. технол. (Гидродинамика, тепло- и массопередача) / Под ред. М.Е. Позина. - М.: Наука, 1965. - С. 266-270.

24. Олевский В.М. Пленочная тепло- и массообменная аппаратура/ под. ред. Олевского Под . ред. В. М. - М.: Химия, 1988.

25. Лаптев А.Г., Ведьгаева И.А. Устройство и расчет промышленных градирен. - Казань: Изд-во КГЭУ, 2004.

26. Рамм В.М. Абсорбция газов. Изд. 2-е, переработ, и доп. - М.: Химия, 1976.

27. Лаптев А.Г., Бажиров Т.С., Саитбаталов М.В. Модель гидродинамических характеристик пленочного течения // Вестник Казанского государственного энергетического университета 2010. №3 - С. 18-23

References

1. Laptev A.G. Modeli pogranichnogo sloja i raschet teplomassoobmennyh processov. - Kazan': Izd-vo Kazansk. gos. un-ta, 2007.

2. Aleksandrov A.I. Massoperedacha pri rektifikacii absorbcii mnogokomponentnoj smesej. L., «Himija»,1975

3. Dolgov A. N., A.G. Laptev Matematicheskaja model' degazacii v nasadochnyh apparatah. // Izv. VUZov. Problemy jenergetiki. №5-6. 2012.

4. Kutateladze S.S., Leont'ev A.I. Teplomassoobmen i trenie v turbulentnom pogranichnom sloe - 2-e izd., pererab. - M.: Jenergoatomizdat, 1985.

5. Patankar S.V., Spolding D.B. Teplo- i massoobmen v pogranichnyh slojah / per. Lykov A.V. - M.: Jenergija, 1971.

6. Garbaruk A.V. L.Ju.V., Strelec M.H. Turbulentnyj pogranichnyj sloj pri odnovremennom vlijanii prodol'nogo gradienta davlenija, vduva (otsosa) i poperechnoj krivizny poverhnosti // Teplofizika vysokih temperatur. 2002. Tom 40. N № 3. s. 436-442.

7. D'jakonov S.G., Elizarov V.I., Laptev A.G. Teoreticheskie osnovy i modelirovanie processov razdelenija veshhestv. - Kazan': Izd-vo Kazanskogo universiteta, 1993.

8. Garbaruk A.V., Lapin Ju.V., Strelec M.H. Prostaja algebraicheskaja model' turbulentnosti dlja rascheta turbulentnogo pogranichnogo sloja polozhitel'nym perepadom davlenija // Teplofizika vysokih temperatur. 1999. T. 37. № 1. s. 87-91.

9. Squire H.B., Young A.D. The calculation of the profile drag of aerofoils. R. & M. No. 1838, British A.R.C., 1938.

10. Leon'tev A.I., Mironov B.P., Lugovskij P.P. Jeksperimental'noe opredelenie kriticheskogo parametra vduva na poristoj plastine // Inzhenerno-fizicheskij zhurnal. 1966. Tom 10. N № 4. s. 447-451.

11. Kutateladze S.S., Leon'tev A.I. Teplomassoobmen i trenie v turbulentnom pogranichnom sloe. - M.: Jenegrija, 1972.

12. Sebesi T. Metod rascheta kojefficienta turbulentnoj teploprovodnosti i turbulentnogo chisla Prandtlja // Teploperedacha. 1973. Tom 95. N № 2.

13. Mejerom i Rottoj .

14. Laptev A.G., Saitbatalov M.V. Analogija perenosa impul'sa, massy i teploty v nasadochnyh jelementah gradiren// Izv. Vuzov. Problemy jenergetiki,2009,-№1-2,-S.140-144.

15. Gilliland E.R., Sherwood T.K. Diffusion of Vapors into Air Streams // Ind. Eng. Chem. 1934. v. 26. NN 5. 516-523.\

16. Vojnov N.A., Nikolaev N.A. Plenochnye trubchatye gazozhidkostnye reaktory (gidrodinamika teplo i massoobmen). - Kazan': "Otechestvo", 2008.

17. Van_Krevelen D.W., Hoftijzer P.J. Studies of gas absorption // Rec. Trav. Chim. 1949. v. 68. N N 4. p. 221-223.

18. Kuznecov M.D. Opredelenie kojefficientov skorosti absorbcii po metodu podobija // Zhurnal prikladnoj himii 1948. t. 26. N 1. s. 48-57.

19. Turhan Je.Ja., Zhideleva K.P. Skorost' absorbcii-desorbcii parov vody i teploperedacha v sisteme gaz - sernaja kislota // Him. prom. 1946. 11. N s. 6-9.

20. Barnet W.I., Kobe K.A. Heat and Vapor Transfer in a Wetted-Wall Tower // Ind. Eng. Chem. 1941. v. 33. N 33. p. 436-442.

21. Johnstone H.F., Pigford R.L. Distillation in a Wetted Wall Column. // Trans. Am. Inst. Chem. 1942. v. 38. N 1. p. 25--51.

22. Zhavoronkov N.M., Maljusov V.A. Issledovanie gidrodinamiki i massoperedachi v processah absorbcii i rektifikacii pri vysokih skorostjah potokov // Teor. osnovy him. tehnologii. 1967. t. 1. N 5. c. 562-577.

23. Dytnerskij Ju.I., Borisov G.S. Issledovanie massoobmena v zhidkoj faze // Processy him. tehnol. (Gidrodinamika, teplo- i massoperedacha) / Pod red. M.E. Pozina. - M.: Nauka, 1965. - S. 266-270.

24. Olevskij V.M. Plenochnaja teplo- i massoobmennaja apparatura/ pod. red. Olevskogo Pod . red. V. M. - M.: Himija, 1988.

25. Laptev A.G., Ved'gaeva I.A. Ustrojstvo i raschet promyshlennyh gradiren. - Kazan': Izd-vo KGJeU, 2004.

26. Ramm V.M. Absorbcija gazov. Izd. 2-e, pererabot, i dop. - M.: Himija, 1976.

27. Laptev A.G., Bazhirov T.S., Saitbatalov M.V. Model' gidrodinamicheskih harakteristik plenochnogo techenija // Vestnik Kazanskogo gosudarstvennogo jenergeticheskogo universiteta 2010. №3 - S. 18-23

Список условных обозначений

b - безразмерный параметр вдува;

- безразмерный параметр вдува в критическом, предотрывном, состоянии турбулентного пограничного слоя;

c - концентрация, кг/кг;

d - внутренний диаметр трубы, м;

D - коэффициент диффузии, ;

g - коэффициент ускорения свободного падения, м/с²;

i - удельная энтальпия, Дж/кг;

- плотность потока массы поперек пограничного слоя, ;

L - длина трубы, м;

q - плотность потока тепла поперек пограничного слоя, ;

- эффективная высота вязкого диффузионного подслоя турбулентного пограничного слоя, в координатах «закона стенки»;

- эффективная высота турбулентного пограничного слоя, в координатах «закона стенки»;

u - продольная составляющая скорости потока (вдоль координаты x) газовой фазы, м/с;

- осредненная по поперечному потока продольная составляющая скорости, м/с;

- продольная составляющая скорости на внешней границе пограничного слоя газовой фазы, м/с;

- динамическая скорость, м/с;

- осредненная скорость движения пленки жидкости, м/с;

- координата по высоте пограничного слоя, м;

- безразмерная координата по высоте пограничного слоя (в масштабах «закона стенки»);

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

- коэффициент теплоотдачи (переноса теплоты), отнесенный к максимальной по высоте пограничного слоя разности энтальпий,

- коэффициент массоотдачи (переноса массы), м/с;

- коэффициент переноса импульса, м/с;

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

- коэффициент турбулентности Прандтля-Кармана;

- функции движущих сил переноса (разности концентраций, энтальпий и скоростей соответственно) между границей раздела фаз и текущей координатой по высоте пограничного слоя, [кг/кг, Дж/кг, м/с];

- максимальная по высоте пограничного слоя движущая сила массопереноса (разность концентраций), кг/кг;

- толщина динамического пограничного слоя, м;

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

- эффективная толщина турбулентного пограничного слоя, м;

- осредненная толщина стекающей пленки жидкости, м;

- приведенная толщина пленки жидкости, м;

- коэффициент трения для закона перепада давления Вейсбаха-Дарси;

- коэффициент динамической вязкости, Па•с;

- коэффициент кинематической вязкости, ;

- плотность, ;

- коэффициент поверхностного натяжения, Н/м;

- напряжение касательного трения (плотность потока импульса поперек пограничного слоя), ;

- относительная скорость газовой фазы в пограничном слое;

- коэффициент относительного трения;

;

- число Прандтля.

Размещено на Allbest.ru