CC BY
394
УДК 533.517.4
МАТЕМАТИЧЕСКИЕ МОДЕЛИ И РАСЧЕТ ГИДРОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ПОГРАНИЧНОГО СЛОЯ
Лаптев Анатолий Григорьевич д.т.н., профессор
ФГБОУ ВПО «Казанский государственный энергетический университет», Казань, Россия
Фарахов Тимур Мансурович к.т.н.
инженер, ООО Инженерно Внедренческий Центр «Инжехим», Казань, Россия
Рассмотрены процессы переноса импульса в ламинарном и турбулентном пограничном слое.
Получены оригинальные модели для расчета касательного напряжения, турбулентной вязкости, коэффициентов трения, толщины пограничного слоя и коэффициентов переноса импульса (импульсоотдачи) в пограничных слоях при различных условиях движения сред. Показаны примеры расчета данных характеристик и согласование с известными экспериментальными данными. Полученные уравнения и методы определения характеристик пограничного слоя могут использоваться при расчетах промышленных тепло- и массообменных аппаратов различных конструкций и другого оборудования
Ключевые слова: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, МАССОПЕРЕДАЧА,
ТЕПЛОПЕРЕДАЧА
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ
К явлениям переноса в пограничном слое не ослабевает интерес многих исследователей, как с точки зрения развития различных теорий, так и для прикладных расчетов гидравлических и тепломассообменных характеристик [1-3].
Движение вязких сред практически всегда связано с явлениями переноса в пограничном слое, где локализуются сопротивления трения, тепло- и массоотдачи. Одним из первых концепцию пограничного слоя в 1904 году предложил Л. Прандтль, исходя из гипотезы о том, что сопротивление, вызванное трением и вязкостью при движении среды существенно зависит от пристеночной зоны (слоя).
UDC 533.517.4
MATHEMATICAL MODELS AND CALCULATION OF THE HYDRODYNAMIC CHARACTERISTICS OF A BOUNDARY LAYER
Laptev Anatoliy Grigorievich Dr.Sci.Tech., professor
Kazan State Power Engineering University (KSPEU), Kazan, Russia
F arakhov T imur Mansurovich Cand.Tech.Sci., engineer
LLC Engineering-Promotional Center “Inzhekhim ”, Kazan, Russia
The processes of momentum transfer in laminar and turbulent boundary layers are considered. We have obtained the proprietary models for calculation of shear stresses, turbulent viscosity, friction factors, boundary layer thickness and momentum transfer coefficients in boundary layers under different conditions of motion of media. The examples of the calculation of the above characteristics and their agreement with known experimental data are shown. The obtained equations and the methods of calculation of the boundary layer characteristics can be used for calculations of industrial heat and mass transfer apparatuses of various designs as well as of some other types of equipment
Keywords: MATHEMATICAL MODELING, MASS TRANSFER, HEAT TRANSFER
Теоретические методы моделирования и исследования явлений переноса в пограничном слое условно можно классифицировать на точные, асимптотические, численные и приближенные. В связи с большим разнообразием конструкций аппаратов химической технологии и одновременно происходящих процессов переноса импульса, массы и энергии получить точные аналитические решения практически
невозможно, особенно для аппаратов большого масштаба. Поэтому для математического моделирования явлений переноса и определения эффективности проводимых процессов в промышленных аппаратах чаще используются приближенные и численные методы. К приближенным методам относятся, например, однопараметрические интегральные методы в теории пограничного слоя, пленочная и пенетрационная модели, модель диффузионного пограничного слоя, методы линерализации уравнений и др. Преимущество приближенных методов заключается в возможности получать необходимые формулы для выполнения конкретных инженерных расчетов [3-7].
Приближенное математическое описание процессов переноса в пограничном слое в первую очередь связано с моделями Прандтля и Кармана, гидродинамической аналогией Рейнольдса и Чилтона-Кальборна и далее моделью диффузионного пограничного слоя Ландау-Левича. Кроме этого, важное значение имеют консервативные свойства пограничного слоя - т.е. весьма слабые зависимости некоторых характеристик осредненного течения по отношению к внешним возмущениям, на которые одними из первых обратили внимание С.С. Кутателадзе и А.И. Леонтьев.
Теоретическая основа рассматриваемых далее подходов заключается в использовании известных свойств консервативности законов трения к возмущающим воздействиям (например, продольному градиенту давления), т.е. структура математического описания элементарных актов
переноса инвариантна к различным возмущениям и масштабу аппарата. Влияние возмущений не изменяет форму математического описания пограничного слоя, а учитывается параметрически. Определение корректирующих параметров выполняется путем удовлетворения баланса импульса [3, 5-7].
В данной статье рассматриваются различные модели пограничного слоя для вычисления напряжения трения и других характеристик на поверхности различных обтекаемых тел при ламинарном и турбулентном режиме. Касательное напряжение трения очень важно при определении коэффициентов гидравлического сопротивления обтекаемых тел, а также при решении тепломассообменных задач с использованием различных вариантов гидродинамической аналогии.
Рассматриваются пограничные слои при движении несжимаемой вязкой жидкости с постоянными физическими свойствами применительно к задачам химической технологии. В статье широко используется понятие коэффициента переноса импульса (импульсоотдачи), которое позволяет получить в аналитическом виде расчетные формулы для коэффициента трения и толщины пограничного слоя. Даны результаты расчета этих характеристик и сравнение с экспериментальными данными. Особый теоретический и практический интерес представляет расчет аппаратов с насадочными (зернистыми) слоями, гидроциклонов, пучков труб в теплообменниках и аппаратов с подводом внешней энергии в контактирующие фазы.
ОПРЕДЕЛЕНИЕ КАСАТЕЛЬНОГО НАПРЯЖЕНИЯ ТРЕНИЯ
Для расчета коэффициентов переноса в тепло- и массообменных аппаратах различных конструкций необходимо определить среднее значение касательного напряжения на стенке тст (или динамической скорости ш =у]тст/р ), где р - плотность среды, кг/м3.
Рассмотрим общий метод определения касательного напряжения трения (локального потока импульса) на стенке - поверхности обтекаемого тела. Пусть ось Оу направлена по нормали к стенке. В этом случае касательное напряжение на ось Оу в любой точке пограничного слоя будет определяться по известному выражению с учетом молекулярного и турбулентного переноса
ёи ! \ёы ! \ёы (л \
т = -ц-— Цт ш—=-р^+^ Н~, (1)
ау ау ау
где
ц , цт - динамическая молекулярная и турбулентная вязкость, Па- с;
V, V т - коэффициенты кинематической молекулярной и турбулентной вязкости, м2/с;
— - градиент скорости.
ау
Касательное напряжение на стенке, то есть при у = 0, может быть выражено, используя уравнение импульсоотдачи (потока импульса) в виде:
Тст = РУ(и„- ист )= РУи^ , (2)
где
у - коэффициент переноса импульса (импульсоотдачи), м/с;
и¥, ист - скорость в ядре потока и на стенке, м/с.
В уравнении (2) учитывая условие прилипания, имеем ист = 0 .
В литературе [8] у имеет размерность кг/(м • с). Однако, как показано ниже удобнее использовать размерность у - м/с.
Еще в середине XIX столетия изучением процесса трения в трубах занимались сначала Дарси, а потом Фаннинг, которые и ввели понятие о коэффициенте трения:
с, = %=, (3)
ри0
при и¥ = иср, где иср - средняя скорость потока, м/с.
Комплекс (3) соответствует критерию Стантона Бт для потока теплоты и критерию Стантона Бд для потока компонента при массообмене.
Полный поток импульса (сила р!,) в единицу времени, который переходит от среды к стенке, записывается в виде:
р = '1р = РУ^~, (4)
где
F - площадь поверхности канала (тела), м .
Аналогично записываются потоки массы компонента и теплоты:
М = pFАC, (5)
б = aFAT, (6)
где
а, р - коэффициенты массоотдачи и теплоотдачи;
АС, АТ - движущие силы соответствующих процессов.
Размерности а - Вт(м • К), а р зависит от способа выражения концентрации С компонента. Если С - кг/м , то р - м/с. Чтобы получить размерность м/с, то а надо разделить на рср, где ср - удельная
теплоемкость, Дж/(кг-К).
При полной гидродинамической аналогии (аналогия Рейнольдса) коэффициенты переноса, приведенные к одинаковой размерности, равны
а = Р = а/рср. (7)
При Бс ф Рг ф 1 используются различные модификации гидродинамической аналогии [3, 7-9], где Бс, Рг - числа Шмидта и Прандтля.
При использовании любых модификаций гидродинамической аналогии необходимо вычисление или экспериментальное определение коэффициента трения с, (3), который связан с касательным напряжением
тст (или динамической скоростью трения и*).
Например, для трубы динамическую скорость можно определить из уравнения баланса сил
АрБ = тст F, (8)
и*
=^=иСРт (9)
где
Ар - перепад давления канала, Па;
5 - площадь поперечного сечения, м2;
2
¥ - поверхность канала, м ;
X - коэффициент гидравлического сопротивления;
Ар записано в форме выражения Дарси-Вейсбаха.
Для пластины из выражения (3) имеем:
и* = и»^ С^2. (10)
Из (9) и (10) следует связь Сг = X/ 4 .
Запишем значение динамической скорости, используя выражение (2)
и* =у]уи». (11)
Рассмотрим подход определения коэффициента переноса импульса
[3,7].
Введем относительный поток импульса т * (у) = т(у)/тст . Тогда из выражения (1) запишем
ист
О <$и
: 5 т *
и» г т
тСт 0 V + (у)
где
5 - толщина гидродинамического пограничного слоя, м.
Сравнивая (12) и (2) приходим к следующему выражению для определения коэффициента импульсоотдачи
т * (у }1у
5 т * о V + (у)
У
Выражение (13) является достаточно общим и позволяет рассчитать коэффициент переноса импульса на основе известного коэффициента молекулярной вязкости, а также известного характера изменения коэффициента турбулентной вязкости V,,(у) и относительного потока импульса т * (у) в пограничном слое. Применение выражения (13) рассмотрено ниже.
В тех случаях, когда значения Сг или X неизвестны (например, в
аппаратах с вводом внешней энергии в контактирующие фазы), находит применение подход вычисления динамической скорости, с использованием средней диссипируемой энергии в единице объема среды [3, 5-7,10-12].
В общем случае для несжимаемой жидкости количество диссипирующейся энергии, т.е. количество механической энергии, превращающейся в тепловую, отнесенное к единице времени и единице объема, определяется по выражению:
8 = Ц
ЭуX 2 + 2 ( Эу у ' 2+2 Эу 2 I2 + ( Эу у + ЭУ 2 "
1 Эх 1 Эу V 1 Эг V 1 Эг Эу V
Эу 2 Эух ^2 ( Эу х ду у ^
+
У X +__________У
Эу Эх
у
(14)
Эх Эг
\ / \ у
где
ц - динамическая вязкость, Па- с; х , у, г - пространственные координаты, м;
у х , V у, у г - составляющие вектора скорости среды по
соответствующим координатам, м/с.
Применение данного выражения для прикладных расчетов сопряжено со значительными трудностями, поэтому скорость диссипации энергии 8 в жидкости, обусловленной вязкостью, при градиенте скорости ёп/ёу записывают в виде [13]
8=де = тёп, (15)
ёу ёу ’
где
е - плотность потока энергии, е = тп; http://ej.kubagro.ru/2012/08/pdf/52.pdf
2
т - касательное напряжение, Па; и - скорость среды, м/с; у - поперечная координата, м.
Ниже рассмотрены различные примеры применения выражения (15) для расчета касательного напряжения (и динамической скорости) при ламинарном и турбулентном режимах.
ЛАМИНАРНЫЙ РЕЖИМ
Касательное напряжение (поток импульса) на стенке по уравнению Ньютона имеет следующий вид:
. (16)
у=0
Из (15), (16) следует уравнение, характеризующее диссипацию
ёи
тст _ М ,
ёу
энергии для є 0 для ламинарного подслоя [11]:
є 0 = М
(17)
Исключив из уравнений (16) и (17) величину производной — и
ёу
подставив найденное при этом значение
?ст = (^80 )1/2 (18)
записывают [11]:
/ \1/4
и. ^ . (19)
I р ^
Первоначально рассмотрим ламинарное течение на пластине и в
трубе.
Диссипация энергии на пластине:
є =
тст^и~ тсти~ (20)
^сл 8 ’
и в трубе:
є _ тстРи¥ 4тсти¥ (21)
^сл ё ’
где
гсл - объем пристенного слоя, м ; ё - диаметр трубы, м.
Отсюда на участке стабилизированного течения в трубе 5 = 0,25ё.
Из (19) - (20) для пластины получим
/ \0,25 , ч0,5
и* —
(22)
Р8
Локальная толщина, пограничного слоя на пластине [14]:
8 » 5 ™ , (23)
\и„
где
х - продольная координата пограничного слоя, м.
Из (22) - (23) значение динамической скорости по координате х:
и*(х)_ 0,447иУ5(VX)0,25. (24)
Или используя выражение (10) с локальным коэффициентом трения
[14]:
и*(х)_ и¥^Ср./2 _ 0,407и¥,75(V/х)0,25, (25)
где
_ 0,332
С,х ‘Ж ’
Яе х _ и¥х/V - число Рейнольдса.
Из выражений (23) - (25) следует, что полное совпадение уравнений (24) и (25) обеспечивается при толщине пограничного слоя (23) с
численным коэффициентом 6,04.
Если в выражении (20) использовать среднее значение толщины пограничного слоя по длине пластины Ь:
8 »1Ь 1^ёх _ 3,33 (26)
Ь0\ите \и„
то из формулы (20) получим
и* _ 0,548иУ5^/Ь)°,25, (27)
а из выражения (10) со средним коэффициентом трения пластины [14] Сг — 1,328Яе-0’5 следует
и* — 0,576м¥,?5 (у/X)0,25. (28)
Расхождение численных коэффициентов в (27) и (28) около 5%.
Для трубы из (21) запишем:
I ¥
и* —
4и¥V ^5 (29)
(30)
й
Или по выражению (9):
г— (8и V ^|0,5
и* — исрЛ/У8 — —р- .
V й
В ламинарном потоке х — 64/Яе; Яе — исрй/V.
Учитывая, что для трубы и¥ — 2иср, получаем полную идентичность
выражений (29) и (30).
Приведенные выражения и сопоставления подтверждают
справедливость уравнения (19) для ламинарного течения с погрешностью не более ±5%.
ТУРБУЛЕНТНЫЙ РЕЖИМ
Получим значение динамической скорости на основе уравнения (15) и известного характера переноса импульса в турбулентном пограничном слое.
В уравнении (15) значение производной найдем из выражения:
йи тст (31)
/ \ йи
тст — - Р^ + ^ )-Т
йу
у—0
где
2
V, Vт коэффициенты кинематической и турбулентной вязкости, м /с. Из (15) и (31) получим:
е — Тст V, — V (у). (32)
p(v+ Ут Г т т Л
Среднее значение скорости диссипируемой энергии в пограничном слое толщиной 5 запишем в виде:
(33)
где подынтегральная функция определяет сопротивление переносу импульса в турбулентном пограничном слое (13), и тогда уравнение (33) получит форму:
Применение выражения (35) рассмотрено в следующих разделах.
ТУРБУЛЕНТНАЯ ВЯЗКОСТЬ
В промышленных аппаратах химической технологии и энергетических установках движение газов (паров) и жидкостей практически всегда происходит при турбулентном режиме. Поэтому для поиска путей совершенствования технологического и энергетического оборудования важное значение имеет достоверный расчет турбулентных течений.
Первые фундаментальные исследования турбулентных течений были выполнены еще в 19 веке Г.Гагеном, О.Рейнольдсом, Ж.Буссинеском, Г. Лоренцем и в начале 20 века развиты многими учеными (Л.Прандтль, Т.Кармен, В.Экман, И.Бюргерс, А.Колмогоров, Х.Драйден, Г.Клаузер, Дж.Тейлор, Г.Шлихтинг, Л.Лойцянский и др). Разработаны различные полуэмпирические теории пристенной турбулентности и отрывных течений. Осреднение по Рейнольдсу уравнений динамики вязкой жидкости приводит к системе уравнений Рейнольдса, называемые уравнениями для первых моментов. В общем случае проблема замыкания этих уравнений сводится к установлению связи тензора турбулентных напряжений с тензором осредненных скоростей деформаций.
(34)
5ру 5у
Отсюда найдем среднее значение динамической скорости:
_ " Y’25
(35)
Одна из первых попыток упростить расчет турбулентных течений принадлежит Буссинеску, который ввел понятие турбулентной вязкости. Далее эта гипотеза была развита Прандтлем и Карманом. При этом делались различные допущения о характере изменения коэффициента v,r (у).
Опытами было установлено, что ут изменяется по поперечному сечению потока, возрастая по мере удаления от стенки до некоторого максимума и затем убывая на оси потока. Часто с достаточной для практических расчетов принимают среднее значение ут »const в ядре потока, т.е. за пределами пристенного слоя.
Рассмотрим теоретический подход определения коэффициента турбулентной вязкости в ядре потока, связанный с диссипацией энергии.
Для энергии турбулентности е изотропного движения за пределами пограничного слоя следует соотношение [15]:
15
е _— v 2
( ЭиЛ
2
. Эу .
\ J j
(36)
Используя приближенное соотношение
Э и
Л 2
Э у
\ У
получим
Э и Э и (37)
Эу Эу
15 Эи Эи (38)
2 ду ду
Принимая турбулентность в канале, близкой к изотропной, Таунсенд
2 3
записал скорость диссипации энергии (м /с ) в виде [16]
е~ И! (39)
I ’
где
V' - пульсационная скорость, м/с;
I - масштаб турбулентности, м.
По определению турбулентная вязкость записывается в виде vт = VI, а касательное напряжение т = рv'м', где V' и и - поперечная и продольная составляющие пульсационной скорости. При изотропной турбулентности имеем V' = и = и*. Тогда из выражений (36) - (39) получим среднее значение коэффициента турбулентной вязкости в ядре потока (т.е. при у > 5):
V, = 7,5 £, (4°)
8
где
_ 2 3
8 - средняя объемная диссипация энергии, м /с .
По модели Прандтля средняя турбулентная вязкость при у > 5 :
vт = хи*5. (41)
Рассмотрим применение выражений (40) и (41) для турбулентного потока в трубе.
Динамическую скорость в трубе вычислим по выражению (9). Средняя скорость диссипации энергии в пристенном слое равна
- = тисрЕ = исри* = иср ^ (42)
рЕ5 5 85
где
2
Е - площадь поверхности трубы, м .
Среднюю толщину турбулентного пограничного слоя найдем, используя логарифмический профиль скорости в трубе [17]:
и = 5,751п -Уи* + 5,5 = 2,51п^* + 5,5. (43)
и*
Отсюда при у = 5 и и = иср получим
X V
5 =—ехр
и*
1 (и
ср - 5,5
V и
2,5
Из (40) с (42) получаем средние значения vт при у > 5
(44)
vт = 0,93иср5Х. (45)
Из (41) с и* (9) запишется:
V
V
= 0,14м ср 5л/х.
'-У
(46)
На рисунке 1 представлены зависимости vт от числа Re по формулам (45) и (46).
Рисунок 1. Зависимость средней турбулентной вязкости в трубе от числа Рейнольдса. 1 - расчет по формуле (45); 2 - расчет по формуле (46).
V
т
Из представленных расчетов видно удовлетворительное согласование результатов расчета средних V,, в ядре потока с использованием теории изотропной турбулентности и модели пограничного слоя Прандтля.
КОЭФФИЦИЕНТЫ ПЕРЕНОСА ИМПУЛЬСА
ЛАМИНАРНЫЙ РЕЖИМ
Для плоского ламинарного пограничного слоя значение у можно определить из выражения (3) с локальным или полным сопротивлением пластины
с,= ^ СГЬ = А&. (47)
Отсюда из (3) и (47) локальный и средний коэффициент переноса импульса
у; = 0,332^5 у = 0,664^. (48)
Данные выражения справедливы для плоского ламинарного пограничного слоя без возмущений. Рассмотрим подход [3], определения у для градиентных течений.
В прикладной аэрогазовой динамике [18, 19] получил применение метод «эффективной длины». В этом случае влияние градиента давления учитывается соответствующим подбором эффективной длины при условии равенства толщины потери энергии пограничного слоя. При известной эффективной длине характеристики пограничного слоя рассчитываются с помощью соответствующих формул для пластины. Следовательно, чтобы получить выражение для у, воспользуемся известными свойствами консервативности законов трения к различным возмущениям. Для этого осредним параметры градиентного потока по длине обтекаемого тела и приведем их к плоскому пограничному слою без гидродинамических возмущений [3, 20].
Предположим, что средний поток импульса на поверхности тела известен. Тогда в рамках плоского ламинарного пограничного слоя запишем [14]:
Из данных выражений найдем осредненную скорость и¥, приведенную к плоскому безградиентному пограничному слою:
В данном случае эквивалентными параметрами градиентного и безградиентного потоков являются касательное напряжение тст и характерный размер телаь, а влияние градиента давления учитывается путем расчета эффективной скорости и¥, на основе известного
(49)
где
и¥= 1,62(тсТ/р)2/3 (ь/-V )1/3.
(50)
1
( тс^ ^3
у = 0,62
коэффициента трения плоской поверхности, т.е. удовлетворяя балансу импульса.
Это дает возможность применить выражение (3)
у = Дс^. (51)
РМ¥
Отсюда с учетом (50) запишем
. (52)
Рь )
Выражение (52) является достаточно общим для ламинарных течений и в частном случае из него следует уравнение (48) с тст = и*.р (10) и С/ (47).
ТУРБУЛЕНТНЫЙ РЕЖИМ
Для турбулентного пограничного слоя на плоской пластине
коэффициенты импульсоотдачи найдем аналогично, как для ламинарного течения, только с локальным и средним коэффициентом трения при турбулентном режиме [21]:
с, = 0058. с = 0074 (53)
Сг. „0,2 ’ С / " • 47
ь
Из (3) с Сг и Сг (53), получим
у х = 0,029м¥ 0,8 (V; )0,2, (54)
у = 0,037М¥0,8 (V ь )0,2. (55)
В общем случае в выражении (35) коэффициент переноса импульса
определим путем интегрирования зависимости (13):
1 = |т*(^¥у . (56)
у 0 V + ^ (у У
Для наглядности получим формулу для у на основе применения двухслойной модели пограничного слоя Прандтля. Запишем интеграл (56) в виде сопротивлений переносу импульса в вязком подслое (V т = 0) и турбулентной области (vт = хи*у):
I = , (57)
У 0 V 51 %и*у’ получим
и*
у =-----? (58)
Я1 + -1п(Л5/Л1) ' 7
%
где
я5 = и* 5/у - безразмерная толщина турбулентного пограничного слоя; я1 = и* 5^ V = 11,6 - безразмерная толщина вязкого подслоя на пластине;
% = 0,4 - константа Прандтля.
Результаты расчета у для пластины по (55) и (58) согласуются с расхождением ± 7%.
Далее с применением трехслойной модели Кармана выражение (56) получит вид:
1 = + 5 / + 2,5 5 А (59)
У 0 у 51 и* у 52 и* У
где
51 = 5у/и*;
52 = 30у/и* .
Первое слагаемое в выражении (59) определяет сопротивление
переносу импульса в вязком подслое толщиной 51 , второе - в переходной
(буферной) области толщиной (52 - 5!), а третье - в турбулентной области
толщиной (5 - 5 2).
Для удобства далее используем безразмерные числа
~ и*51 ~ и*5 2 Л1 = ; ^-2 = ,
V V
где согласно модели Кармана ^ = 5, ~ = 30.
После интегрирования (59) получим
1 = А. + ^Цпб + —1п(д5/Я2) (60)
у и* и* и*
или
Y 2,792R1 + 2,5ln(Rs/ R2)'
(б1)
Как известно, в моделях Прандтля и Кармана в вязком подслое vT = о.
Согласно теории Ландау и Левича, подтвержденной Дайслером, турбулентность в вязком подслое описывается функцией
Vl = (о,124у+)4, при у+< 6, (62)
тогда после интегрирования аналогичного выражения (59) с vт(у) (62) получим
Y =----~ Ы* ( ), (63)
2,288R1 + 2,5ln(R5/ R2)
где
R1 = 6 для плоского пограничного слоя.
Как следует из расчета по выражениям (61) и (63) учет затухания турбулентных пульсаций в вязком подслое в модели Кармана практически слабо влияет на значение коэффициента переноса импульса. Наиболее значительное влияние наблюдается для двухфазных сред [3, 5, 6].
В работах [3,7] выполнено интегрирование выражения (13) с различными функциями турбулентной вязкости в пограничном слое. Например, с функциями vт (у) Owen P. получено:
Y =
U*
11,7З + 2,5 ln[(0,4u*5 - Зv)/(5v)]'
(б4)
Или с функциями Дайсслера и Ханратти
Y
_____и*____ (65)
5,309 + 2,51п Д5
Ниже показано, что использование полученных выражений для коэффициента импульсоотдачи позволяют получать в явном виде ряд важных характеристик турбулентного пограничного слоя.
ТОЛЩИНА ПОГРАНИЧНОГО СЛОЯ
Локальная и средняя толщина, турбулентного пограничного слоя на пластине имеет вид [21]:
0,37х _ 1 Ч 037 хйх 0,205!
5 х=•8=Ч ^
Х
Ке/0,2
(66)
где
Ч - длина пластины, м;
х - продольная координата;
числа РейнольдсаЯе х = и¥х/у; ЯеI = и¥Ч/у.
С применением выражения для у, толщину пограничного слоя найдем из потокового соотношения (2):
тст = и* р = руАи ¥,
(67)
где
Аи ¥ - средняя движущая сила переноса импульса; Аи ¥ = и ¥ ± и гр, где
и¥ - скорость на оси потока, м/с;
гр
- скорость на границе раздела, м/с. Для твердых поверхностей
и гр = 0.
На основе (67) с учетом коэффициента переноса импульса (65) получено [7, 22]:
- 5,309 . (68)
) _
Аналогичные выражения для 5 или к5 можно получить из (67) с коэффициентами у (58) и (61):
К5 = ехр
0,4
К5 = К1 ехР
*5 = *2 ехр
где
*1 = 5;
К2 = 30;
0,4
(
0,4
и ~
— - 2,792*
и*
(69)
(70)
К = 11,6.
При использовании вышеприведенных выражений для каналов необходимо учитывать, что скорость на оси потока [21]:
и¥ » 1,15и или и¥ = иср + 3,75и*.
В итоге уравнение (35) решается итерационным методом с одним из выражений для 5. Результатом решения являются средние значения динамической скорости и толщины турбулентного пограничного слоя.
Для проверки полученных уравнений выполнены сравнения результатов расчета толщины пограничного слоя для пластины (рис. 2, 3).
На рисунке 2 представлена зависимость толщины пограничного слоя на пластине от продольной координаты х .
Толщина пограничного слоя вычислялась по уравнению (66) и выражению (69).
В уравнении (69) динамическая скорость определена по уравнению (10), где локальный коэффициент трения на пластине С^ = 0,058/яех2 .
Й. м
0,03
ода
0.02
0,0 [ 5
0.01
0,005
------1-----1-----4______I_____I-__
0.2 0,4 0.& КО
Рисунок 2. Зависимость толщины пограничного слоя на пластине от продольной координаты. Рабочая среда - воздух; иср = 7,5 м/с; 1 - расчет по уравнению (69); 2 - расчет по уравнению (66).
На рисунке 3 представлены зависимости толщин пограничного слоя на пластине от числа Рейнольдса.
0.005
й, м
0,0035
0,0045
0.0025
1
Не
I _
2104 6-104 105 4-105
Рисунок 3. Зависимость толщины пограничного слоя на пластине от числа Рейнольдса. Рабочая среда - воздух; х = 0,1 м; 1 - расчет по уравнению (66); 2 - расчет по уравнению (69).
В результате проверки уравнения (69) получены удовлетворительные согласования с известными выражениями. Аналогичные результаты дают и уравнения (68) и (70).
Т.к. полученные выражения (68)-(70) дают близкие результаты (±5), то для вычислений удобнее пользоваться более компактной формулой (69), полученной на основе применения модели Прандтля.
РАСЧЕТ ДИНАМИЧЕСКОЙ СКОРОСТИ
В качестве примера найдем по выражению (35) значение и* для турбулентного потока на пластине.
Средняя диссипация энергии на пластине:
¥5 5 ’
где
2
Е - площадь поверхности пластины, м ;
- = Тст ^¥ = и*?ри¥
(71)
3
У5 = 5Е - объем пограничного слоя, м .
Из (35) и (71) следует выражение (11).
Сравнение результатов расчета значения и* для пластины по выражениям (10) и (11) с у (58), (61) или (65) показывает их согласование с расхождением не более 5-7%.
Аналогичную проверку сделаем для потока в трубе, где средняя диссипация энергии имеет вид (42). Тогда получим:
| 3 Л 0,25
и, = |и^1 (72)
8
V у
На рисунке 4 представлена графическая зависимость динамической скорости от числа Рейнольдса в трубе.
Динамическая скорость определяется решением системы уравнений (72) с у (69) - кривая 2 и известного уравнения (10) (кривая 1), где коэффициент сопротивления определяется по формуле Блазиуса.
Рисунок 4. Зависимость динамической скорости от числа Рейнольдса в трубе. Рабочая среда - воздух; аэ = 0,05 м; 1 - расчет по (9), 2 - расчет по (72); 3 - экспериментальные данные [21].
Далее рассмотрим выражение для коэффициентов импульсоотдачи (58), (61), (63) и (65). Например, из выражения (58) следует, что
сопротивление переносу импульса сосредотачивается в вязком подслое
Я
и*
и в турбулентной области пограничного слоя
2,51п(Я5/ Я)
и *
Многочисленные расчеты у показывают, что сумма слагаемых в знаменателе выражения (58) примерно составляет я1 + 2,51п(я5/я1 ) = 15 +18. Если приближенно взять среднее значение, то получим у»и*/16,5. Аналогичная оценка следует и из выражений (61), (63) и (65). Тогда в формуле (72) можно записать:
г 3 Л 0,25
иср1и*] = 0,196иср X0,333. (73)
и*
*ср
132
Как показывают расчеты значений и* по (73) данное выражение дает удовлетворительные результаты (± 20) для различных контактных устройств при турбулентном режиме движения потоков.
ПОЛУЭМПИРИЧЕСКИЕ МЕТОДЫ. ТУРБУЛЕНТНЫЙ РЕЖИМ
Значение динамической скорости можно получить из решения системы уравнений (35) с у по (58), (61) или (65) и 5 по (68) или (70). Для сокращения времени расчета рассмотрим полуэмпирический подход определения динамической скорости на основе использования средней диссипируемой энергии.
Для турбулентного движения потоков используют соотношение:
8 0 = к 4 е.
(74)
где
к — эмпирический коэффициент пропорциональности.
Тогда выражение для динамической скорости представляют в виде:
г Л0,25
и* = к
^ (75)
.р J ’
где коэффициент к находят экспериментальным путем для каждой
конструкции аппарата [11,23].
Далее рассмотрим турбулентный режим движения среды на пластине и в трубе.
Выражение (75) для динамической скорости справедливо с коэффициентом пропорциональности к.
Для пластины из (71) и (75) получим:
и* — к
( 2 1 и* vu¥
0,25
5
к —-
Отсюда запишем коэффициент пропорциональности
/ п г\0,25
[ и¥Сг5 1
^и„/ 5) , где
0,25
2v
— ^О/2)0,25,
(76)
(77)
д5 — и¥5^ - безразмерная толщина пограничного слоя.
При средних значениях для пластины (66) 5 — 0,205Ь/Яе0,2 и
Сг — 0,074/Яе0,2 из (77) получим к — 0,34Яе0,2.
В трубе:
к—
и¥\ /X,
0,25
— 0,42(Яе X)
0,25
(78)
0,187
(Хи¥ V/ 2ё )0'
При X — 0,31б/яе0,25 из (78) следует к — 0,31ЯеС
Полученные выражения для коэффициента пропорциональности к — /(яеп), (п — 0,187 и п — 0,2) в зависимости (75) объясняют полученные экспериментальным путем значения к —1,6+2,7 для аппаратов различных конструкций [11, 23].
Для определения коэффициента к в (75) теоретическим путем запишем поток импульса (67) с динамической скоростью (75). Получим
(~ ^ 8У
0,5
уЛи¥ — к
Отсюда найдем коэффициент
к — (уЛи ¥),
(- V0,25
(79)
(80)
Из (75) с коэффициентом (80) получаем динамическую скорость, аналогичную (11). Таким образом, используя различные подходы, приходим к единой формуле для расчета и* в виде (11).
Ьйр://е].киЬа§го.гц/2012/08/рёГ/52.рёГ
КОЭФФИЦИЕНТЫ ТРЕНИЯ
Рассмотрим применение выражений для коэффициентов импульсоотдачи в форме коэффициентов трения.
Из (11) с у (58) получим динамическую скорость в виде:
и* —---^---------? (81)
*1 + -1п(Л5/Л1) ' '
%
Выражение (81) связывает динамическую скорость трения и* с
характеристиками турбулентного пограничного слоя в двухслойной модели Прандтля - *5, я1, х. Расчет и* по (81) удовлетворительно согласуется с и*, полученному по выражению (10) и* — /(С/).
С использованием выражения (81) можно записать коэффициент трения
С/ — 2*1 + 11п(Л5/ *1)
1_ х
Аналогично с выражением (61)
Сг — 2[2,792*1 + 2,51п(*5/ *2 )]"2. (83)
и с формулой (65)
Сг — 2 [5,3 09 + 2,51п *5 ]-2. (84)
Полученные выражения для коэффициента трения (82)-(84) позволяют исследовать влияние толщины вязкого подслоя (Я1 или *1), буферной области *2 и всего пограничного слоя *5 на перенос импульса с учетом различных возмущений.
Так, например, на шероховатой поверхности с выступами, соизмеримыми с толщиной вязкого подслоя, выражение (82) получит вид
С/ — 2
Д. + ^А. I-2, (85)
X *1
1ш у
где
*1ш - толщина, вязкого подслоя с учетом шероховатостей,
*1ш — (51 - кV*/V ;
Ьйр://е].киЬа§го.ги/2012/08/рёГ/52.рёГ
к - высота выступов шероховатостей, м (прин < 5!).
Из выражения (85) следует, что при увеличении к коэффициент трения возрастает, что соответствует известным экспериментальным данным.
ЗЕРНИСТЫЕ НАСАДОЧНЫЕ СЛОИ
Рассмотрим применение полученных выражений для и*, у и V,, при движении газов или жидкостей в насадочных или зернистых слоях.
Насадочные слои тепломассообменных аппаратов подразделяются на регулярные и нерегулярные (засыпанные внавал). Нерегулярные насадки состоят из элементов с размерами 50мм и менее и имеют различную форму. Турбулентность в таких слоях возникает значительно раньше, чем в гладких трубах и границе ламинарного режима соответствует значение Яе —15-40. Полностью развитый турбулентный режим наступает при значениях Яе от 2 103 до 6 103, где Яе — исрёэ/V . При
обычно встречающихся на практике значениях Яе от 40 до 2 103 соответствует переходному режиму [24].
Для ламинарного режима движения газа или жидкости (Яе < 40) динамическую скорость найдем, используя выражение (19). Тогда из (19),
(86) и (87) получим
и* —
{к V/4 (\ з Л1/4
Ари V уа
pH
ср_
У
^исруау
8е3
(86)
где
н - высота насадочного слоя, м;
2 3
ау - удельная поверхность насадки, м /м ; єсв - удельный свободный объем;
йэ — 4Єсв/ау .
Коэффициент гидравлического сопротивления насадочного или зернистого слоя находится экспериментально для каждого типа элемента
(засыпки) и обычно дается в виде функции х — / (Яе) [12,24,25,26]. В работе [26] связь между коэффициентами имеет вид х — 4 /э.
Коэффициент переноса импульса для ламинарного режима в насадке запишем из выражения (52) и (88)
Современные насадки, которые по форме значительно отличаются от колец и седел, способствуют развитию интенсивной турбулентности в ядре потока. Однако, элементы насадок, образованными изогнутыми полосами металлической ленты шириной 3-5 мм, имеют особенности по формированию пограничного слоя. Известно, что при контакте турбулентного потока с обтекаемой поверхностью на передней кромке образуется ламинарный пограничный слой. При достижении некоторого критического размера ламинарное движение в пограничном слое становится неустойчивым и развивается турбулентность. Для элементов насадок, расположенных далее первого ряда, начальный участок ламинарного слоя подвержен турбулентным воздействиям из ядра потока. Такой пограничный слой классифицируется как псевдоламинарный. В первом приближении среднюю толщину псевдоламинарного пограничного слоя можно вычислить по выражению (26), а затем уточнить с учетом воздействия турбулентных пульсаций. Тогда турбулентную вязкость при у > 5 найдем с использованием модели локальной изотропной турбулентности по уравнениям (40), (42), получим
где динамическую скорость для насадочного слоя можно записать, используя приближенное выражение (73), тогда получим
(87)
(88)
, Л 0,333
2 V1 1
и* — 0,196исрХ°,333 — 0,196и,
ср Н 2й
V
У
У
где коэффициент сопротивления записан в известной форме [24] для насадочного слоя, используя перепад давления.
Выражение (88) с динамической скоростью (89) запишется в виде:
V,, = 0,2885исрХ °’666. (90)
В неупорядоченных насадочных слоях, образованных из металлических насадок, развивается высокая интенсивность турбулентности (Яе > 2 • 103), что дает возможность использовать модель локальной изотропной турбулентности со средней объемной скоростью диссипации энергии выражение (40), где е запишем в виде:
- Лр5иср Лриср = Чр_ (91)
р^сл-св рНесв 2^э-св
где
3
гсл - объем слоя, м ;
5 - площадь поперечного сечения, м2; н - высота слоя, м.
Перепад давления слоя насадки в выражении (91) записан, используя аналог выражения Дарси-Вейсбаха
н ри2
Лр = X
н риср Ч
(92)
где
dэ = 4-св/ау - эквивалентный диаметр насадки, м. Из выражения (40) с и* (89) и - (91) получим:
7,5(0,196исрХ0,333 )4
= 7,5<12-^0,1964м X0,333 = 0,113и*^е, (93)
, , / “ср'^ ^^u*dэ£cв
ЛРиср/ (^-свР)
или
V
= 0,022Яе -свХ0,333. (94)
св
V
В итоге из (93) следует аналог выражения Прандтля (41), где
5 » 0,Ш^-св/х .
Динамическая скорость при турбулентном режиме в насадке находится из решения системы уравнений (35) и (58). Многочисленные http://ej.kubagro.ru/2012/08/pdf/52.pdf
расчеты для разных типов насадок показывают, что значение и* можно вычислить по простой зависимости (75), где для насадочного слоя к »1,8 - 2,0. Удовлетворительные результаты так же дает выражение (73) с коэффициентом сопротивления насадочного слоя.
АППАРАТЫ С МЕХАНИЧЕСКИМ ПЕРЕМЕШИВАНИЕМ И ГИДРОЦИКЛОНЫ
Величину средней диссипируемой энергии в аппарате с мешалкой выразим в виде отношения мощности N, затрачиваемой на перемешивание, к объему жидкости в аппарате
где
в - диаметр аппарата, м; н ж - высота уровня жидкости, м.
При известном значении - динамическую скорость на стенке можно вычислить из решения системы уравнений (35), (58) и (69). Расчеты показывают с достаточной для прикладных расчетов точностью, что можно использовать выражение (75) с к » 2,0. Тогда получим
(96)
Аналогичное выражение получено для гидроциклона
(97)
где
Лр - перепад давления в гидроциклоне, Па;
5 - площадь входного патрубка, м ; иж - скорость жидкости в патрубке, м/с; уж - объем жидкости в гидроциклоне, м .
Таблица - ОСНОВНЫЕ ПОЛУЧЕННЫЕ УРАВНЕНИЯ
№ формулы Показатель
Динамическая скорость трения, м/с
(11) в общем виде: и* =у1 чи¥.
(35) и* — Г е8у V Р ) 0,25
(22) на пластине (ламинарный режим) / \ 0,5 -1 и¥У 1 и* — 1 1 . 1 8 )
(24) и* (х) — 0,447и¥,75 (V х)0,25.
(27) и* — 0,548и^,75 (V ^)0,25.
(29) в трубе (ламинарный режим) — Г 4и¥ V10,5 и* — 1 1 . 1 d )
(72) в трубе (турб и* — V улентн 3л 1 иср ^У 8 ) ый режим) 0,25
(73) обобщенная приближенная формула и* — 0,196иср 10’333. 7 ир
(75) обобщенная форму: и* — Для трубы к — 0,31Яе0,187 Для пластины к — 0,34Яе0,2 Для насадочного слоя к —1,8 - 2,0 В аппарате с мешалкой к — 2,0 В гидроциклоне к —1,55 Шахматные и коридорные пучки а (ту к ( ^' V р) труб рбулентный режим) 0,25 к — 2,0
(86) для насадочного сл и* — V для насадочного сл и* —1,8 оя (лам DPuсрV ^ рн ) )я (тур б Г DPuсрV V Рн инарный режим) 0,25 улентный режим) 10,25 )
Продолжение таблицы 1
№ формулы Показатель
(89) и* = 0,196иср ’ ср ( Ар Н 2' риср 2ёэ V у Г1 л / ) 0,333
(96) аппара ( и* = 2 V іт с мешалкой \ 0,25 4 | ржЛ 2 Нж )
(97) гид и* = 1,55 роциклон ( Ар5ижуж'' V рж^ж ) 0,25
Коэффициенты переноса импульса
(48) Ламинарный режим на пластине ух = 0,332^^5 у = 0,664^^.
(52) обобщенная формула у = 0,62 для ламинарного режима Г л 0,333 Тст^ V Р1 ) '
(58) турбулентный режим и*
У і . Я + -1п(Я5/ Я1) %
(61) и*
1 2,792Я + 2,51п(Я5/Я2 )
Толщина пограничного слоя - турбулентный режим
(69) Я5 = ^1 ехр 0, (и Л \ — -Я1 Vи* ) -
(70) Я5 = Я2 ехр 0, ( и ~ Л1 4 — - 2,792Я1 . V и*
(44) Я5 = ехР " ( 0,4 V и¥-5,5 Л и* )
Турбулентная вязкость
(40) в общем виде
4
_ , vт = 7,5—.
8
Продолжение таблицы 1
№ формулы Показатель
(45) для трубы ут = 0,93нср5Х.
(88) в насадочном слое 2 ут = 7,55—, 40 < Яе < 2 103 иср
(94) ^ = 0,022 Яе8свх0,333, Яе > 2 103 V
Коэффициенты трения (турбулентный режим)
(82) ^=2 *1 + *1) _ х т-2
(83) = 2[2,792.~ + 2,51п(* 5 / *2 )]"2.
(85) на шерох СТ = 2 ^оватой поверх] ( 1 * ^ (*1ш + -1П --5V х *1ш у ости -2
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ
Выполненные в данной работе аналитические исследования моделей пограничного слоя и сделанные сопоставления с экспериментальными данными и известными расчетными зависимостями позволяют сделать следующие обобщения и выводы:
1. Теоретически обоснована форма записи динамической скорости (19) для ламинарных течений с погрешностью ±5%.
2. Для турбулентного режима получено аналитическое выражение
(35), где основными параметрами являются средняя скорость диссипации энергии, толщина, пограничного слоя и коэффициент импульсоотдачи.
3. На основе применения теории локальной изотропной турбулентности получено выражение (40) для расчета средней турбулентной вязкости. Показано согласование с известными экспериментальными данными.
4. Получено выражение для коэффициента импульсоотдачи (52) для ламинарных течений с различными возмущениями, где основным параметром является касательное напряжение на стенке.
5. Путем использования моделей Прандтля, Кармана, Ландау, Левича и Дайсслера и интегрирования уравнения (13), получены выражения для коэффициентов импульсоотдачи (58), (61) и (63) для турбулентного режима.
6. На основе баланса импульса в пограничном слое с использованием выражений для коэффициентов импульсоотдачи получены уравнения для расчета средней толщины гидродинамического пограничного слоя. Показано согласование с известными результатами.
7. С использованием выражения (35) показаны примеры расчета динамической скорости на пластине и в трубе и сравнение с известными выражениями.
8. Получены выражения для коэффициента трения, которые содержат характеристики пограничного слоя и позволяют учесть шероховатость поверхности при турбулентном режиме.
9. Рассмотрены ламинарный и турбулентный режимы движения потоков в неупорядоченных насадочных и зернистых слоях и получены выражения для расчета динамической скорости, коэффициента импульсоотдачи и турбулентной вязкости. Аналогичный подход рекомендован и для аппаратов с механическим перемешиванием и гидроциклонов.
Приведенные в данной работе уравнения и методы расчета гидромеханических характеристик пограничного слоя могут использоваться в различных вариантах расчета промышленных аппаратов и другого оборудования в зависимости от решаемых задач.
Литература
1. Ломов С. А., Ломов И.С. Основы математической теории пограничного слоя. -М.: Изд-во Моск. ун-та, 2011. - 456с.
2. Репик Е.У., Соседко Ю.П. Турбулентный пограничный слой. Методика и результаты экспериментальных исследований. - М.: Физматлит, 2007. - 312с.
3. Лаптев А.Г. Модели пограничного слоя и расчет тепломассобенных процессов. - Казань: Изд-во Казанск. ун-та, 2007. - 500с.
4. Дильман В.В. Методы модельных уравнений и аналогий в химической технологии / В.В. Дильман, А. Д. Полянин. - М.: Химия, 1988. - 304 с.
5. Дьяконов С.Г. Теоретические основы и моделирование процессов разделения веществ / С.Г. Дьяконов, В.И. Елизаров, А.Г. Лаптев. - Казань.: Изд-во Казанск. ун-та, 1993. - 437с.
6. Лаптев А.Г. Модели переноса и эффективность жидкостной экстракции / А.Г. Лаптев. - Казань: Казан. гос. энерг. ун-т, 2005. - 299 с.
7. Лаптев А.Г., Николаев Н.А., Башаров М.М. Методы интенсификации и моделирования тепломассообменных процессов. - М.: Теплотехник, 2011. - 287с.
8. Бенедек П., Ласло А. Научные основы химической технологии. Под ред. П. Г. Романкова и М.И Курочкиной. - Изд. «Химия», 1970. - 376с.
9. Деменюк С.Л. Теплообмен и гидравлическое сопротивление при турбулентном течении в каналах / С. Л. Деменюк . - Санкт - Петербург: Бюро hi Пром, 2010. - 133с.
10. Брагинский Л.В. Перемешивание в жидких средах: физические основы и инженерные методы расчета / Л.Н.Брагинский, В.И.Бегачев, В.М. Барабаш. - Л.: Химия, 1984.
11. Соколов В.Н. Г азожидкостные реакторы / В.Н.Соколов, И.В. Доманский. -Л.: Машиностроение, 1976. - 216с.
12. Лаптев А.Г. Разделение гетерогенных систем в насадочных аппаратах / А.Г Лаптев, М.И.Фарахов. - Казань: Казан. гос. энерг. ун-т, 2006. - 342с.
13. Ландау Л.Д. Теоретическая физика / Л.Д.Ландау, В.М. Лившиц: учеб. пособие: в 10 т. - Т.6. Гидродинамика. - 5-е изд. - М.: Физматлит, 2003. - 560с.
14. Лойцянский Л.Г. Ламинарный пограничный слой / Л.Г. Лойцянский - М.: Физматгиз, 1962. - 410с.
15. Хинце И.О. Турбулентность, ее механизм и теория. - М.: Наука, 1963 -
680 с.
16. Таунсенд А.А. Структура турбулентного потока с поперечным сдвигом /
A. А. Таунсенд: пер. с англ. / под ред. А.Н. Колмогорова. - М., 1959. - 360с.
17. Лойцянский Л.Г. Механика жидкости и газа / Л.Г. Лойцянский. - М.: Наука, 1987. - 840с.
18. Авдуевский В.С. Основы теплопередачи в авиационной и ракетно -космической технике / В.С. Авдуевский, Б.М. Галицейский, Г.А.Глебов и др. под ред.
B.К. Кошкина. - М.: Машиностроение, 1975. - 330с.
19. Краснов Н. Ф. Основы прикладной аэрогазодинамики. Кн. 2. Обтекание тел вязкой жидкостью / Н.Ф Краснов, В.Н Кошевой, В.Ф Захарченко и др..; под ред. Н.Ф Краснова. - М.: Высшая школа, 1991. - 350с.
20. Лаптев А.Г. Модель массоотдачи в жидкой фазе насадочных колонн / А.Г. Лаптев, Р.Р. Насыров, Т.М. Фарахов // Тепло - массообменные процессы и аппараты химической технологии: межвуз. темат. сб. науч. тр. - Казань, 2005. - С.4-8.
21. Шлихтинг Г. Теория пограничного слоя. / Г. Шлихтинг. - М.: Наука, 1974. - 714с.
22. Лаптев А.Г., Башаров М.М., Фарахова А.И. Эффективность
турбулентной сепарации мелкодисперсной фазы в тонкоотслойных отстойниках // Энергосбережение и водоподготовка, 2011, - №5(73). - С. 43-46.
23. Доманский И.В. Обобщение различных случаев конвективного
теплообмена с помощью полуэмпирической теории турбулентного переноса / И.В.
Доманский, В.Н. Соколов. // Теор. Основы хим. технол. - 1968. - Т.2, № 5. С.761-767.
24. Рамм В.М. Абсорбция газов. 2-е изд. / В.М. Рамм. - М.: Химия, 1976. -
655с.
25. Аэров М.Э. Аппараты со стационарным зернистым слоем:
Гидравлические и тепловые основы расчета / М. Э. Аэров, О. М. Тодес, Д. А. Наринаский. - Л.: Химия, 1979.
26. Фарахов Т.М. Многофункциональные контактные устройства смешения котельного топлива с присадками и очистки газовых выбросов ТЭС: дис. канд. тех. наук / Т.М. Фарахов - Казань: КГЭУ, 2011.
