книги / Тепловые и гидродинамические процессы в колеблющихся потоках
..pdfQ/2^;fiA
Рис. |
105. |
Коэффициент |
ослабления |
Рис. |
106. |
|
Зависимость |
коэффициента |
||||||||
при малоамплитудных колебаниях дав |
ослабления от числа Рейнольдса и |
|||||||||||||||
ления |
жидкости в |
канале: |
|
частоты |
колебания: |
|
|
|
|
|||||||
О — Re = 1,Ь 103; |
д — Re = 104; |
П — |
• |
— |
f = |
200 Гц; |
X |
— |
f = |
300 Гц; |
V |
— |
||||
Re = |
2- 104; |
# - |
Re = |
5- 104; |
X - |
f |
= |
400 Гц; |
о — |
f |
= |
900 |
Гц; □ — |
f |
= |
|
Re = |
10* |
|
|
|
|
= |
600. Гц |
|
|
|
|
|
|
|
|
Рис. 107. |
Зависимость |
безразмерного |
Рис. 108. |
Зависимость |
безразмерного |
||||||
коэффициента ослабления |
рd0 от ча |
коэффициента ослабления рd0 от ча |
|||||||||
стоты колебания f и числа М0 для |
стоты f |
и числа М для |
шероховатой |
||||||||
гладкой трубы d0 = 75 мм: |
|
|
трубы (h = |
242 мкм) d0 = |
55 мм: |
||||||
Д — М 0 = |
0; X — М0 = |
0,1; |
□ |
— м 0= |
• |
_ м 0 = |
— |
0 , 1 ; Д —М0 = |
0,15; О — М0= |
||
= 0,2; А |
— м 0 = 0,3; О |
— М0 = |
0,4 |
= |
0,2; |
□ |
М0 = 0,3; А |
— |
М 0 = 0,4 |
ционально и число Рейнольдса. Таким образом, в данных опытах увеличение коэффициента (5 от числа М0 обусловлено увеличе нием числа Рейнольдса, которое изменялось в диапазоне 105— 7,5 -105.
Шероховатость на поверхности канала существенно увеличи вает коэффициент ослабления. На рис. 108 приведены значения безразмерного коэффициента ослабления $d0 от числа М0 и ча стоты колебания / для трубы диаметром 55 мм с шероховатостью
221
Рис. 109. Зависимость безразмер ного коэффициента ослабления pdQ от динамической скорости v* в шероховатой трубе d0 = 75 мм для различной степени шероховатости:
О |
— h = |
242 мкм; |
□ |
— h = 96 мкм; |
д |
— h =вв |
43 мкм; • |
— |
h =*в 4 мкм |
поверхности 242 мкм. В этом случае коэффициент ослаб ления существенно увеличи вается. На рис. 109 приве дена зависимость коэффици ента ослабления $d0 от динамической скорости v* =
= К т^/р для трубы диамет ром 75 мм для различной степени шероховатости при частоте колебания 1000 Гц. Как и в случае стационар
ного режима течения, влияние шероховатости выявляется только в том случае, если высота бугорков шероховатости превышает вязкий слой. Это следует из приведенного рис. 109: чем больше степень шероховатости, тем при меньшем значении v* наблюдается влияние динамической скорости на коэффициент ослабления.
В приведенных выше экспериментальных данных влияние амплитуды колебания на коэффициент ослабления не учитыва лось, и предполагалось, что амплитуда колебания мала (Аи0/и0
«1).
Коэффициент поглощения высокоинтенсивных колебаний боль
шой амплитуды при неизотермическом течении газа в канале экспериментально исследован авторами данной монографии. Ис следования проводились в области наибольшего влияния колеба ний на теплообмен на экспериментальной установке (рис. 110). Экспериментальная установка состоит из экспериментального участка, системы подачи рабочего тела, источника колебаний давления и системы измерения. В качестве рабочего тела исполь зовался воздух, который поступал к экспериментальному участку из системы баллонов высокого давления через предварительный подогреватель и понижающий редуктор. Экспериментальный участок представлял собой специально калиброванную трубку из нержавеющей стали, с толщиной стенки 0,15 мм и внутренним диаметром 12 мм. Общая длина экспериментального участка составляла 2045 м. Обогревалась экспериментальная трубка переменным током низкого напряжения, непосредственно про пускаемым по экспериментальной трубке.
На экспериментальной установке измерялись следующие пара метры: температура воздуха на входе и выходе из эксперимен-
222
тального участка, температура наружной поверхности трубки, электрическая мощность, подводимая к экспериментальной трубке, частота и амплитуда колебания давления на входе и выходе из экспериментального участка, осредненный по времени расход воздуха и перепад давления на экспериментальном участке.
Средний по времени расход воздуха измерялся нормальной диафрагмой, которая предварительно тарировалась по мерному баку. Чтобы исключить влияние колебания давления на измере ние расхода, диафрагму устанавливали на значительное расстоя ние от пульсатора. Э. д. с. термопар измерялась полуавтомати ческим потенциометром Р2/1 и цифровым милливольтметром по стоянного тока Ш-15-13. Для измерения температур использова лись хромель-алюмелевые термопары с диаметром провода 0,1 мм.
Амплитуда колебания давления в экспериментальном участке регулировалась изменением давления подачи и сменных сопел на вы>оде из экспериментального участка.
Опыты проводились в диапазоне изменения основных пара* метров; давление составляло (1—20) 104 кгс/м2; температура поверхности экспериментальной трубки 0—600° С; число Рей нольдса 104—5 -105; относительная амплитуда колебания давления 0—0,5; частота колебаний 50—1000 Гц. Основное исследование
Рис. 110. Схема экспериментальной установки для исследования высокочастот ных колебаний давления и их влияние на теплообмен в канале:
1 |
— |
подогреватель; |
2 |
— |
усилитель; |
3 — |
частотомер; |
4 — усилитель-преобразователь; |
|||
5 |
— |
осциллограф; |
6 |
— |
спектрометр; 7 — |
милливольтметр; 8 — потенциометр; 9 — пере |
|||||
ключатель; 10 — |
экспериментальный |
участок; |
11 — |
индуктивные датчики давления; |
|||||||
12 |
— |
пульсатор; |
13 |
— |
электродвигатель; |
1 4 |
— |
измерительный трансформатор; 15 — |
|||
силовой трансформатор; |
17 — вольтметр; |
1 6 |
— |
вакуумметр |
223
|
‘______ 1 — |
_____ |
Рис. 111. |
Зависимость коэффициента |
ос |
||||||
/,2 |
лабления |
от амплитуды |
колебания |
|
мас |
||||||
2Ч0Гц |
х X |
совой скорости в канале |
d — 12 мм |
для |
|||||||
|
|||||||||||
|
х Х |
о А з 0 |
различных |
резонансных |
гармоник |
|
(80, |
||||
0,8 |
|
O r |
|
160 и 240 |
Гц) |
|
|
|
|
||
0 ° |
|
|
|
|
|
|
|
|
|
||
|
|
|
было проведено в области резо |
||||||||
1- |
}£ь80Г >/н -------- |
|
нансных |
частот |
колебания, |
|
где |
||||
ОЛ |
|
|
|
|
|
|
|
|
|
|
|
_____ |
|
|
|
наблюдается максимальная ампли |
|||||||
|
|
|
|
туда |
колебания |
давления, |
т. е. |
||||
|
100 . 200 |
300 |
400 |
Д (ры) |
|
(ры)0. |
Коэффициент |
|
ос |
||
|
лабления определялся по методу |
||||||||||
|
|
а - м |
|||||||||
|
|
^9uWc |
двух датчиков. На |
рис. 111 |
при |
||||||
|
|
|
|
ведена |
типичная |
характеристика |
|||||
изменения коэффициента ослабления |
р в канале (диаметр 12 мм) |
от амплитуды колебания массовой скорости для первых трех резонансных гармоник (fs =*■80; 160 и 240 Гц). С увеличением амплитуды и частоты колебания коэффициент ослабления увели чивается. Причем при сравнительно больших значениях Д (ры) коэффициент ослабления достигает 0,8—1,0, что на порядок пре вышает соответствующие квазистационарные значения. Такие большие коэффициенты ослабления обусловлены тем, что при ин тенсивных колебаниях измеряемый коэффициент ослабления вклю чает в себя потери в волне, на стенке канала, а также и потери, обусловленные нелинейными и газодинамическими эффектами. С увеличением числа Рейнольдса коэффициент ослабления также увеличивается.
Согласно экспериментальным данным, полученным для глад ких цилиндрических каналов различных размеров (d0 = 9,8; 12; 20; 30 мм), коэффициент ослабления р в области резонансных частот колебания и при Д (ры)0/(ры)0 ^ 1 в диапазоне исследо ванных параметров (приведенных выше) обобщается критериаль
ной |
зависимг. |
ь |
|
|
|
|
к |
- с«■.-[> + XAS'-•>]"■ ( x f |
{ЦТ Re«‘’ |
<507> |
|
где |
С = |
8,2- 10-8; п1 = 0,517; п2 = 0,517; |
п3 = 0,4; п4 = |
0,1; |
|
«5 = |
0,167; /\ |
(P*V0 / — средняя j по длине |
канала относитель- |
||
ная |
амплитуда массовой скорости; Rea = a0d/v. |
|
Измерение коэфициента гидравлического сопротивления
Для экспериментального измерения осредненного по времени коэффициента гидравлического сопротивления в первом приближении можно воспользоваться стационарной методикой, измеряя в эксперименте перепад давления на определенной длине экспериментального участка. В этом случае коэффициент гидра-
224
влического сопротивления целесообразно определять как отно шение касательного напряжения на стенке канала к полной кине тической энергии потока. Если предположить, что статические отборы давления измеряют среднее квадратическое по времени давление, то, интегрируя уравнение (39) по длине измеряемого участка, получим
(Рг) — (Pi) + Ро«о (Мог — м01) + \ Р о |
Ди\ |
|
- \P o ~ 2 |
||
П |
л , |
(508) |
= — y |
xw Aft, |
откуда согласно уравнению (36) на длине канала ДL среднее зна чение коэффициента гидравлического сопротивления
|
/ |
Aw? 4 |
Д«1' |
xw |
l<Pi> — <р2>1 — Ро«о (“ ог — “ 01) + Ро < |
—к - |
— < Ро |
|
|
|
|
Ро»о |
|
|
|
8 |
|
|
(509) |
|
|
|
где (pi)—(р2) — средний по времени перепад статического давле ния на длине канала AL; черта сверху означает среднее по длине значение кинетической энергии потока; индексы 1 и 2 относятся к параметрам потока соответственно на входе и выходе из канала. При такой методике естественно возникает вопрос, на сколько точно измеряется осредненное по времени давление посредством стационарных методов и в какой степени одномерная модель тече ния описывает движение реальной жидкости. Изменение осредненного по времени статического перепада давления в цилиндри ческом канале (диаметр 12 мм) на экспериментальной установке (см. рис. 110) представлено на рис. 112. Из рисунка видно, что
Рис.. 112. Изменение осредненного по времени перепада давления в канале d0 — 12 мм в зависимости от частоты колебания:
О —•р 0 = 15* 104 кгс/м2; Д — р 0 = 8* 104 кгс/м*
15 Б. М. Галилейский
Рис. И З . Зависимость относительного коэффициента гидравлического сопро
тивления Къ = £ / |о от параметра Г =
= 0 ,1 /|0 (Ди0/« 0)2 При бо/бк = 0,75:
О — экспериментальные точки
225
с изменением частоты перепад давления изменяется довольно сложно. Вблизи резонансных частот, соответствующих первой и второй резонансным гармоникам (80, 160 Гц), наблюдается отри цательный перепад давления. По-видимому, это объясняется тем, что при резонансных колебаниях отборы давления расположены не в одинаковых условиях, поскольку амплитуда колебания ско рости переменна по длине канала. Обработка результатов опытов показала, что относительный коэффициент гидравлического со противления К%= |/£ 0 зависит от относительной амплитуды колебания массовой скорости А (ры)0/(ри)0, числа Рейнольдса Re
осредненного движения |
и соотношения |
толщин |
колебательного |
|
и стационарного вязкого слоя 8к/80, т. е. |
|
|||
M |
t ) |
! = ' ( w ' |
Re' - | ) - |
<510> |
На рис. 113 приведены результаты обработки эксперимен тальных данных по относительному коэффициенту сопротивления трения.
С увеличением относительной амплитуды колебания скорости А (ры)0/(ры)0 величина увеличивается при прочих равных условиях.
5. ТЕПЛООБМЕН ПРИ ВЫСОКОЧАСТОТНЫХ КОЛЕБАНИЯХ ПРИ ТУРБУЛЕНТНОМ РЕЖИМЕ ТЕЧЕНИЯ
Теоретический анализ теплообмена в колеблющихся потоках при турбулентном режиме течения значительно услож няется, поэтому существующие полуэмпирические теории, кото рые достаточно удовлетворительно с практической точки зрения описывают стационарные осредненные по времени турбулентные потоки, вряд ли могут быть использованы для исследования колеблющихся турбулентных потоков. Как известно, турбулент ные стационарные потоки включают в себя большой диапазон частот колебаний параметров потока. Возможно, что колеблю щийся поток определенной частоты может избирательно вклю чаться в энергетический спектр турбулентных пульсаций. Тур булентные пульсации, частоты которых близки к частоте выну жденных колебаний, могут возбуждаться, тогда как турбулент ные пульсации с другой частотой, наоборот, могут подавляться под действием вынужденных колебаний.
В случае высокочастотных колебаний (при больших со), когда толщина колеблющегося пограничного слоя много меньше, чем толщина стационарного вязкого подслоя (6к/80 <С 1), в первом приближении можно считать, что влияние колеблющегося потока сказывается только в вязком подслое, тогда как в турбулентном ядре профиль скорости квазистационарный. В этой области воз можен приближенный теоретический анализ, основанный на ме тоде расчета для ламинарного колеблющегося пограничного слоя.
226
Вслучае средних частот колебаний (при сравнительно малых значениях <о) толщина колеблющегося пограничного слоя больше, чем толщина вязкого подслоя (6к/60 Г> 1), и вынужденным ко лебаниям будет подвергаться не только вязкий подслой, но и про филь скорости в турбулентном ядре потока. Поскольку турбулент ная вязкость в ядре потока значительно больше, чем в пристеноч ной области, то деформация профиля скорости в основном будет наблюдаться вблизи поверхности тела. При ускорении потока профиль скорости будет более заполненным, а при замедлении, наоборот, менее заполненным по сравнению со стационарным случаем. Колебания профиля продольной скорости будут сопро вождаться колебаниями поперечной скорости потока.
Таким образом, влияние колеблющегося потока на теплообмен при турбулентном режиме течения обусловлено, во-первых, изме нением коэффициента турбулентной теплопроводности посредством воздействия вынужденных колебаний на энергетический спектр турбулентных пульсаций, -во-вторых, вынужденные колебания продольной скорости генерируют колебания поперечной скорости потока, что приводит к дополнительному поперечному переносу количества тепла.
Вслучае, когда колеблющийся поток возбуждается стоячей волной, влияние этих колебаний на теплоотдачу, как и в случае ламинарного режима течения, переменно по длине стоячей волны.
Это объясняется особенностью стоячих волн. Поскольку ампли
туда колебания скорости изменяется по длине стоячей волны от нуля в узлах до максимума в пучностях, то теплоотдача в пуч ностях скорости стоячей волны будет максимальной, а в узлах — минимальной. Если не вдаваться детально в механизм передачи тепла в высокочастотных колеблющихся потоках, то критериаль ное уравнение для относительного коэффициента теплоотдачи
можно представить |
в виде |
|
|
* = |
= |
= / ( - “ -. Re, |
(511) |
Используя соотношение (469) для течения в трубе, получим
Следовательно, критериальное уравнение (511) можно записать в виде
К = Nu*Nu f ( ~ ° l , Re, R e»,^ , м). |
(512) |
Таким образом, критериальное уравнение (512) является наиболее общей функциональной зависимостью, определяющей основные параметры, по которым следует моделировать процесс теплообмена в области высокочастотных колебаний.
15* |
227 |
Теплообмен в условиях низкочастотных колебаний
Результаты исследования влияния низкочастотных ко лебаний на средний по времени и по длине коэффициент тепло отдачи в теплообменнике приведены в работе [60]. Исследование проводилось в трубчатом теплообменнике, состоящем из 36 или 12 медных труб с внутренним диаметром 13,5 мм и длиной 940 мм. В качестве теплоносителя использовалась вода (t = 17° С), обогрев которой осуществлялся паром, подаваемым в межтрубное про-, странство. Число Рейнольдса изменялось в пределах 2,9 • 103— 17,5-Ю3. Колебания расхода воды создавались емкостью, вклю ченной в систему подачи воды, в результате изменения давления в ней. Частота колебаний составляла от 0,4— 1,1 Гц. Амплитуда колебания скорости воды определялась по уровню жидкости
в мерном баке |
= 2—9^. Результаты опытов для различных |
чисел Рейнольдса приведены на рис. 114. Из графиков ри сунка видно, что с увеличением относительной амплитуды коле бания скорости теплоотдача увеличивается и при Аи/и0 ^ 8 составляет для развитого турбулентного потока К = 1,4. На этом же рисунке приведена кривая, соответствующая квазистационарной зависимости, определенной по осредненной по вре мени скорости колебания в степени 0,8, соответствующей квазистационарному мгновенному изменению коэффициента теплоот
дачи:
т
(Аы) = J (1 -f- Л«о/ sin (at)0’9 dt. 0 о
Рис. 114. Зависимость относительного среднего .по длине теплообменника .коэф фициента теплоотдачи (К) от относительной амплитуды колебания скорости
(Ди/цр):
a) |
Re = 2900, |
f = |
0,5-* 1 Гц; б) |
Re = |
6,6-10»; f = |
0,47 |
+ 0,75 Гц; в) Re = 1,75104; |
/ = |
0,48^0,93 |
Гц; |
ф — давление |
пара |
1,62 KFC/CM8; |
0- —. |
р =в 1,1-25 кгс/см2 |
228
к
Рис. 115. Зависимость относительного коэффициента теплоотдачи К (сред него по длине теплообменника) от относительной амплитуды колебаний
скорости |
Аи0/ и 0: |
|
|
|
X — |
охлаждение А Т = |
75° С; |
• — наг |
|
ревание |
Д Т = 75° С, |
О — |
нагревание |
|
А Т = |
15° С |
|
|
Рис. 116. Зависимость относительного коэффициента теплоотдачи К (сред него по длине теплообменника) от числа Рейнольдса при частоте коле бания / = 2,68 Гц:
О — экспериментальные точки
Модель квазистационарной теплоотдачи удовлетворительно описывает эксперименты только в области сравнительно малых значений относительных амплитуд колебания скорости. Аналогич ные исследования были проведены в трубчатом теплообменнике «труба в трубе» [73]. Диаметр внутренней трубы составлял 51 мм, длина 600 мм. Обогрев воды (до t = 75° С), поступающей в теплообменник, осуществлялся паром, подаваемым в рубашку теплообменника. Колебания воды создавались посредством порш невого насоса, включенного в систему подачи воды. Амплитуда колебания измерялась, как и в предыдущей работе, по колебанию
уровня |
воды в мерном |
баке. Частота |
колебаний |
составляла |
|||
1,75 |
Гц, |
относительная |
амплитуда |
колебания скорости |
1— 1,6 |
||
при |
числах Рейнольдса |
3 *104—8,5 • |
104. |
Результаты |
опытов при |
||
ведены на рис. 115, из которого видно, что при Ди/и0 #=» |
1,4 на |
блюдается максимальное увеличение теплоотдачи (К *=» 2,0). В ра боте [76 ] приведены результаты исследований на трубчатых тепло обменниках с диаметром внутренней трубы 12,7 мм и длиной тепло обменника в пределах 1,2—4,8 м. В качестве рабочего тела исполь зовалась вода, нагреваемая паром до температуры 95° С. Число Рейнольдса изменялось в пределах 8 • 103—8 • 104. Колебания рас хода воды генерировались поршневым насосом. Частота колеба
ний |
составляла 0,67; 1,35; |
2,68 Гц, амплитуда колебания — |
1-4-2. |
На рис. 116 приведены |
результаты опытов по теплоотдаче |
в зависимости от числа Рейнольдса. Максимальное увеличение теплоотдачи соответствует малым числам Рейнольдса и наиболь шей частоте (в данным опытах f — 2,68 Гц). В работе [34] экспе риментально исследовался теплообмен в теплообменнике с диа метром трубы 27 мм и длиной 1 м. В качестве теплоносителя
229
использовались вода и водные растворы сахара с концентрацией 30 и 60%. Обогрев теплообменника осуществлялся паром. Колеба ния создавались поршневым насосом. Частота колебаний состав ляла 0,17— 17 Гц, относительная амплитуда колебания скорости изменялась в пределах 0,027—82,7. Число Рейнольдса состав ляло 8*102—3 • 104; число Прандтля 2,5— 112. Результаты опытов по теплоотдаче в зависимости от амплитудного числа Рейнольдса приведены на рис. 117. Их можно обобщить следующей зави симостью:
Nu = |
0,57 • 10-6ReJ'6Re°Af-f,574.io-‘Re„ Pr°.« |
, (513) |
где ReA„ = |
Au0d/v — амплитудное число Рейнольдса. |
|
Аналогичные эффекты увеличения интенсивности теплообмена при турбулентном режиме течения наблюдаются и в случае вибра ции поверхности теплообменного аппарата. В этих случаях, как и при ламинарном режиме течения, в качестве характерной ско рости выбирается амплитуда скорости вибрации, определяемая как Avof = Af (А — амплитуда вибрации поверхности; f — ча стота вибраций). Результаты исследования влияния поперечных вибраций на теплообмен в кольцевом зазоре и в канале при пере ходном и турбулентном режимах течения приведены в работе [71 ]. В работе [70] исследовалась теплоотдача как при ламинар ном, так и при турбулентном режиме (рис. 118) течения воды в кольцевом зазоре с внешним диаметром 36 мм и в канале труб
Ни |
чатого |
теплообменника |
диаметром |
|
12,7 мм при вибрациях внутренней |
||||
Рг°'Чеп0 |
||||
трубы с частотой 5—45 Гц и амплитудой |
||||
|
||||
|
вибраций |
1—5,27 мм. На рис. 118 при |
||
|
ведены результаты опытов по тепло- |
|||
|
отдаче’для «чисел^Рейнольдса |
в канале |
||
|
Re = 7970; 9050 и в кольцевом зазоре |
к
Рис. 117. Зависимость среднего относительного коэффициента теплоотдачи в теплообменнике от амплитудного числа Рей нольдса ReAu:
• — экспериментальные точки
1,3
V
V
ол
о го W 6 0 8 0 10 0 12 0 М 16 0
Рис. 118. Зависимость относительного коэф фициента теплоотдачи в цилиндрическом и кольцевом канале от амплитуды поперечной скорости вибрации стенок канала Ап = A0f:
X |
— труба Re = |
7,87* 108; зазор Re = 1,54» 104; |
О |
— труба Re = |
9,63* 108, зазор Re ==. 1,9* I04 |
230