КУРСОВАЯ РАБОТА

"Исследование теплоотдачи цилиндра в закрученном потоке"

Введение

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

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

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

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

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

1. Теоретическое обоснование

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

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

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

, (1)

где - число Нуссельта ( – коэффициент теплоотдачи конвекцией от цилиндра к закрученному потоку; – внешний диаметр цилиндра; – коэффициент теплопроводности среды); – критерий, определяющий движение жидкости и теплоотдачу в закрученном потоке (j_m – максимальное значение центростремительного ускорения j в системе; l_m – характерный линейный размер, определяющий положение j_m относительно поверхности теплоотдачи; ∆T– разность температур в точках потока, где инерционное ускорение достигает максимального и минимального значений; β, ν – коэффициенты объемного расширения и кинематической вязкости среды; – критерий Прандтля, – коэффициент температуропроводности).

В рассматриваемой задаче минимальное ускорение имеет место на поверхности теплоотдачи (цилиндра) и равно нулю. Физические константы в уравнении (1) определяются по средней температуре потока, практически равной температуре Е точке потока, где j=j_m .

При сравнительно небольших температурных напорах ∆Т (в условиях, близких к изотермическим) сомножитель (1+Т) незначительно отличается от единицы, а при небольших изменениях ∆Т, даже если и отличается от единицы, то изменяется сравнительно мало. Поэтому изменение критерия Р полностью определяется безразмерным комплексом

. (2)

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

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

.

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

. (3)

Знак плюс (+) соответствует нагреву среды; знак минус (–) – охлаждению.

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

. (4)

Анализ критерия S с использованием результатов исследований аэродинамики циклонных камер и созданной на их основе методики аэродинамического расчета позволил установить, что он может быть заменен на число Рейнольдса Re_φm , рассчитанное по максимальной вращательной скорости в рабочем объеме циклонной камеры w_φm . Эта связь определяется уравнением

(5)

Здесь ,

где ν_m – коэффициент кинематической вязкости среды в точке, где вращательная скорость равна максимальной; D=D(η_я , В) – безразмерный комплекс, зависящий от аэродинамических характеристик потока; – безразмерный радиус осесимметричного ядра потока (r_Я – радиус ядра потока; r_Ц – радиус цилиндра; r_{φ m} – радиус, определяющий положение w_φm ); – безразмерный параметр.

Величины η_Я , В зависят от основных геометрических параметров циклонной камеры и могут быть рассчитаны, либо определены по экспериментально найденному профилю вращательной скорости в объеме камеры. Значения D в зависимости от η_Я , В приведены в табл. 1 [1].

Замена S на Re_{φ m} дает определенные преимущества в обработке опытных данных и использовании полученных уравнений подобия, так как в число Рейнольдса входит величина w_{φ m} , определяемая по методике аэродинамического расчета циклонных камер. Максимальная вращательная скорость циклонного потока для рассматриваемой задачи является фактически скоростью набегающего потока, которая обычно принимается в качестве характерной во внешних задачах конвективного теплообмена.

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

, (6)

где A, m, n – постоянные коэффициенты, определяемые из опыта.

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

, (7)

В условиях опытов величина k_β сохраняется примерно постоянной и равной 0,885 (влияние коэффициента на расчетные уравнения невелико и уменьшается с понижением показателя n).

В ряде случаев представляет интерес использовать в уравнении (6) не сам комплекс D, а заменяющие его более простые характеристики, определяемые методикой аэродинамического расчета циклонных камер. В частности хорошие результаты дает использование коэффициента крутки в ядре потока ε_Я :

,

где W_{φ Я} – линейная вращательная скорость на радиусе r_Я определяющем границу ядра потока.

2. Описание экспериментальной установки и методики измерений

Исследовательская часть работы выполняется на специальном экспериментальном стенде. Циклонно-вихревая камера (рис. I) представляет собой гладкостенную вертикальную металлическую модель внутренним диаметром D_К =2R_К = 310 мм.

Ввод воздуха в камеру осуществляется тангенциально расположенными к внутренней поверхности ее рабочего объема входными каналами (шлицами) 6 с двух диаметрально противоположных сторон. Ширина прямоугольных шлицев (сопел) l_ВХ и их высота h_ВХ могут варьироваться специальными вкладышами. При этом соответственно меняется и суммарная площадь входа потока , где – число входных каналов. Отвод газа из модели осуществляется через плоский торец с круглым осесимметричным выходным отверстием, безразмерный диаметр которого может меняться в диапазоне от 0,2 до 0,6.

Схематический чертеж цилиндра – калориметра приведен на рис. 3. Длина калориметра – 400 мм, толщина стенки рабочего участка – 2 мм. Наружный диаметр цилиндра-калориметра при проведении опытов может изменяться от 45 до 140 мм. Верхний торец рабочего участка калориметра теплоизолирован текстолитовым диском, нижний – охранным участком.

Принципиальная схема тепловых измерений приведена на рис. 4.

Как видно из рис. 4, греющий пар из электрокотла через электрические основной и выносной перегреватели по подводящей трубке поступает в рабочий участок калориметра. Для исключения возможности попадания в калориметр жидкой фазы на входе пара в рабочий участок поддерживается и непрерывно контролируется протарированной медь – константановой термопарой небольшой перегрев ()°С. Отсчет ЭДС термопары производится переносным потенциометром типа ПП-63. Для обеспечения стока конденсата нижняя часть парового пространства калориметра имеет небольшую конусность (рис. 3). Отвод излишнего пара и паровоздушной смеси осуществляется в охранный участок калориметра, а оттуда через отводящий штуцер в дренажную систему.

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

3. Порядок проведения опытов

Пуск экспериментальной установки производится путем включения в сеть электрокотла, предварительно напитанного водой. Уровень воды в водомерном стекле должен быть не ниже средней отметки. Продувочные краники дренажных магистралей пароперегревателя, паропровода и калориметра устанавливаются в открытое положение. Последовательным включением нагревательных элементов котла устанавливается электрическая мощность, близкая к максимальной (сила тока по показаниям амперметра должна составлять около 30А). Одновременно с котлом включается и основной пароперегреватель на 50% его мощности. Процесс разогрева ведется в данном режиме до достижения избыточного давления в электрокотле порядка (0,2…0,3) кгс/см ((19,6…29,4) кПа). После этого производится включение стенда по воздушной стороне. Для этого необходимо закрыть заслонку на. воздухопроводе установки и осуществить запуск воздуходувки с электрощита управления. Изменение частоты вращения, а следовательно, и производительности осуществляется вручную реостатами ступенчатой и плавной регулировки. Максимальная нагрузка вентилятора устанавливается таким образом, чтобы показания амперметра на щите управления не превышали 45–50 А. Полный напор, развиваемый воздуходувкой, при этом составит (420…480) мм вод. ст. ((4,12…4,71) кПа). Далее по указанию преподавателя устанавливается соответствующий режим работы, определяемый расходом воздуха через установку. Регулировка расхода может быть осуществлена как заслонкой, так и изменением частоты вращения привода вентилятора. Прежде чем приступить к производству замеров, необходимо вывести установку на стационарный гидродинамический и тепловой режимы. Для этого обычно требуется не менее 40–60 мин. В это время необходимо тщательно продуть паровое пространство установки, после чего закрыть краники на пароперегревателе и паропроводе. С помощью реостатов плавной регулировки мощности и продувочного краника калориметра установить необходимое избыточное давление в нем путем совмещения нижнего уровня жидкости гидрозатвора с отметкой на смотровом стекле, включенном в магистраль сбора конденсата. После этого приступить к регулированию перегрева пара реостатами выносного, а если необходимо, и основного пароперегревателей, ведя контроль по показаниям потенциометра и сравнивая их с градуировочной характеристикой термопары. Убедившись в достижении стационарного режима по паровой и воздушной сторонам, приступают к проведению эксперимента.

В начале и конце каждого опыта производится запись показаний следующих приборов: термометров в воздушных коллекторах перед измерительным участком трубопровода t и циклонной камерой t_вх , термометра, измеряющего температуру окружающей среды t_ос , микроманометра, измеряющего перепад давления в сужающем устройстве (нормальной диафрагме) ∆Р, дифференциальных водяных манометров, измеряющих избыточные статические давления в контрольных сечениях перед сужающим устройством Р_сд , во входных каналах Р_с.вх и на боковой поверхности циклонной камеры Р_с.ст , барометра-анероида В.

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

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

Для изучения гидродинамической обстановки в рабочем объеме камеры производится также измерение скоростей и давлений в одном или нескольких радиальных сечениях. Для тепловой задачи основной интерес представляет величина максимальной вращательной скорости w_{φ m} , достигаемой на радиусе r_{φ m} . Поэтому снятие полей скоростей и давлений производят на одном режиме (обычно максимальном), устанавливая величину r_{φ m} , а на остальных режимах замеры осуществляют лишь на данном радиусе.

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

1) ввернуть ножку координатника зонда в бобышку измерительного сечения;

2) установить зонд в горизонтальное положение по уровню;

3) поворотом лимба координатника относительно указательной стрелки в соответствующую сторону установить величину угловой поправки к показаниям зонда по данным его аэродинамической тарировки;

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

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

При снятии полей скоростей и давлений производится отсчет и запись следующих величин: радиального расстояния y = R_K – r (r – радиус точки замера) от стенки рабочего объема камеры до точки измерения (по рейке координатника); угла скоса потока φ (угла между направлением вектора полной скорости потока и направлением горизонтали), отсчитываемого по лимбу координатника; полного избыточного давления потока A₂ перепада давления между центральным и боковым отверстиями A₁ пропорционального скоростному напору в точке замера.

4. Обработка результатов наблюдений

Определение расхода воздуха через установку

Модуль сужающего отверстия

.

Плотность воздуха перед сужающим устройством

кг/м³ ,

где g – ускорение силы тяжести, g = 9,80665 м/с² .

Поправочный множитель на расширение измеряемой среды ε для нормальных диафрагм:

,

где k – показатель адиабаты измеряемой среды (для воздуха k = 1,4).

Величина перепада давления:

мм вод. ст.

где ρ_ж – плотность рабочей жидкости микроманометра (спирта) при 20 °С, г/см³ ; 0,8095 – плотность эталонного спирта при 20 °С, г/см³ ; 0,0011 – коэффициент объемного расширения спирта, 1/°С; t_ж – температура рабочей жидкости, принимаемая равной t_ОС ; k_изм – множитель, равный синусу угла наклона измерительной трубки микроманометра; Н_пок – показание прибора, отсчитываемое по измерительной шкале, мм.

Теоретический расход воздуха через экспериментальную установку:

м³ /с,

где α_и – исходный коэффициент расхода; k₂ – поправочный множитель на шероховатость измерительного трубопровода; k₃ – поправочный множитель на неостроту входной кромки диафрагмы.

Значение α_и зависит от типа сужающего устройства и модуля m:

при

.

Произведение поправочных множителей k_2,3 = k₂ k₃ находят из таблицы 2 [1] по значениям m, D_ТР .

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

м² /с.

Число Рейнольдса, отнесенное к диаметру измерительного трубопровода:

Действительный расход воздуха:

м³ /с,

где k₁ – поправочный множитель на число Рейнольдса, определяемый по номограмме, представленной на рис. 5 [1].

Для основной диафрагмы, используемой на стенде (m= 0,259), k₁ можно рассчитать также по приближенному аппроксимационному уравнению

Плотность воздуха во входных каналах циклонной камеры

кг/м³ .

Среднерасходная скорость воздуха в шлицах

м/с.

Безразмерное избыточное статическое давление газа на боковой поверхности циклонной камеры

Безразмерное (отнесенное к динамическому давлению на входе) статическое давление во входных каналах

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

где Р_П – перепад полного давления в циклоне, мм вод. ст.;

Р_ДИН – динамическое давление во входных каналах, мм вод. ст.

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

Входное число Рейнольдса вычисляется по формуле

Расчет распределений скоростей и давлений в объеме циклонной камеры (по данным замеров пневмометрическим цилиндрическим зондом)

Избыточное статическое давление потока в произвольной точке

мм вод. ст.,

где А₁ , А₂ – показания пневмометрического насадка, мм вод. ст.;

k_Ц , k_Ц-Б – поправочные коэффициенты, определяемые тарировкой насадка.

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

кг/м³ .

Полная скорость потока в точке замера

м/с.

Безразмерная осевая составляющая полной скорости

Безразмерная вращательная составляющая полной скорости потока

.

Избыточное безразмерное статическое давление в точке замера

.

Безразмерное избыточное полное давление в точке замера

Размерная величина максимальной вращательной скорости

м/с

Обсчет опытных данных по конвективному теплообмену

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

Вт,

где G_К – масса конденсата, кг;

τ – время опыта, с;

r_n – теплота парообразования, Дж/кг.

Температуры насыщения

⁰ С,

где Р_ИЗБ – избыточное давление в калориметре, создаваемое столбом воды в гидрозатворе, Р_ИЗБ =Н_ГИДР , мм вод. ст.

Дж/кг.

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

где – приведенная степень черноты системы; ε_Ц , ε_СТ – соответственно степени черноты калориметра и поверхности циклонной камеры(, ); C₀ =5,77 – постоянная Стефана-Больцмана, ; Т_Ц – абсолютная температура внешней поверхности калориметра, К; T_СТ – абсолютная температура внутренней поверхности камеры (); F_Ц – площадь поверхности теплообмена калориметра, F_Ц =πd_Ц l_Ц , м² ; F_СТ - площадь боковой поверхности циклонной камеры, , м² .

м² ,

м² .

Конвективный тепловой поток:

Вт.

Приращение температуры воздуха, охлаждающего калориметр

⁰ С,

где С_Р – средняя массовая теплоемкость воздуха при постоянном давлении, равная 1005 Дж/(кг °С) в диапазоне температур от 0 до 60 ⁰ С.

Средняя температура циклонного потока

⁰ С.

Коэффициент теплоотдачи конвекцией от калориметра к закрученному воздуху

Вт/(м^{2 0} С).

Число Нуссельта

где – коэффициент теплопроводности воздуха при средней температуре потока.

Вт/(м^{2 0} С).

Плотность воздуха на радиусе r_{φ m}

кг/м³ ,

где Р_СТ – избыточное статическое давление воздуха на радиусе r_{φ m} , мм вод. ст.

Коэффициент кинематической вязкости воздуха на радиусе r_{φ m}

Число Рейнольдса

Статистическая обработка данных

После обработки опытных данных по теплоотдаче для всех исследованных в работе режимов по числу Рейнольдса устанавливается функциональная связь между числами Nu и Re в виде зависимости (6).

Наиболее обоснованным и широко распространенным в практике научных исследований видом аппроксимации опытных данных является метод «наименьших квадратов», связанный со статистическим законом распределения случайных ошибок эксперимента.

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

Применим данный метод для окончательной обработки результатов экспериментов, используя значение показателя степени m при комплекса. D, полученное ранее. (Следует заметить, что m, как правило, значительно меньше n и погрешность в определении m не сказывается существенно на результатах расчетов чисел Nu.).

Запишем зависимость в виде:

где

,

,

Все опытные данные обобщим в координатах , (рисунок 1 приложения).

где z – число опытов.

5. Расчет статистических показателей

Общая дисперсия (рассеивание) точек относительно среднеарифметического значения lgKo:

Дисперсия точек относительно линии регрессии:

где – дисперсия относительно среднего значения lgKo

Среднеквадратичное отклонение точек от прямой линии регрессии:

Граница доверительного интервала по Ко:

где t_α – коэффициент Стьюдента, определяемый по табл. 4 [1] в зависимости от величины выборки (числа опытов);

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

Коэффициент корреляции:

Среднеквадратичная погрешность коэффициента корреляции.

Оценка достоверности прямолинейной связи:

, t_α =2,78; – условие выполняется

6. Расчет погрешностей измерений

При можно принять, что

.

Учитывая, что

.

Если исследуемая величина является функцией нескольких независимых переменных а, b, …, c, измеряемых в процессе эксперимента, относительная ошибка измерения определяется следующим образом:

.

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

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

,

где – средняя квадратичная относительная погрешность измерения расхода. Последняя определяется по формуле:

,

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

Средняя квадратичная относительная погрешность коэффициента расхода для диафрагм:

,

где – среднеквадратичные относительные погрешности соответственно исходного коэффициента расхода, а также поправочных коэффициентов К₁ , К₂ , К₃ ; Определяются по номограммам на рис. 8–11 [1].

Среднеквадратичная относительная погрешность поправочного множителя на расширение измеряемой среды ε:

,

где n – множитель, принимаемый равным 0,02 для диафрагм с .

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

Средняя квадратичная погрешность измерения перепада давления на диафрагме:

где ∆h – основная допускаемая абсолютная погрешность дифференциального микроманометра, мм вод. ст.

Средняя квадратичная относительная погрешность значения плотности ρ:

Предельная относительная погрешность определения погрешность числа Нуссельта

Предельная относительная ошибка определения коэффициента теплоотдачи конвекцией:

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

где

,

,

где

где

где , φ – рад.