/
Text
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ, МЕХАНИКИ И ОПТИКИ Г. Н. Дульнев, С. В. Тихонов ОСНОВЫ ТЕОРИИ ТЕПЛОМАССООБМЕНА Учебное пособие Санкт-Петербург 2010
Дульнев Г. Н., Тихонов С. В. Основы теории тепломассообмена, - СПб: СПбГУИТМО, 2010. - 93с. В пособии излагаются основы тепло- и массообмена, при этом основное внимание уделено кондуктивному процессу переноса тепловой энергии. Рассматриваются тепловые модели, позволяющие решать прикладные инженерные задачи. Приводятся примеры по применению теории к решению характерных задач. Рекомендовано Учебно — методическим объединением по университетскому политехническому образованию в качестве учебного пособия для студентов высших учебных заведений, обучающихся по направлению «Техническая физика». Протокол № УМК-30-10 от 10.06.2010г. САНКТ-ПЕТЕРБУРГ II В 2009 году Университет стал победителем многоэтапного конкурса, в результате которого определены 12 ведущих университетов России, которым присвоена категория «Национальный исследовательский университет». Министерством образования и науки Российской Федерации была утверждена Программа развития государственного образовательного учреждения высшего профессионального образования «Санкт- Петербургский государственный университет информационных технологий, механики и оптики» на 2009-2018 годы. ©Санкт-Петербургский государственный университет информационных технологий, механики и оптики, 2010 ©Г.Н.Дульнев, С.В.Тихонов, 2010 2
Введение Тепломассообмен — учение о процессе распространения тепла и массы в природе. Основы были заложены еще акад. М.В.Ломоносовым, он создал механическую теорию теплоты и установил законы сохранения массы и энергии. В дальнейшем учение о тепломассообмене развивалось как один из разделов технической физики. В XIX веке в связи с изобретением паровой машины, паровой турбины и двигателя внутреннего сгорания основное внимание уделялось превращению энергии в работу. В дальнейшем в связи с ростом мощности отдельных агрегатов тепловых машин стала возрастать роль процессов тепломассообмена. Этим процессам большое внимание стали уделять и в других отраслях техники - строительной, металлургической, холодильной, машиностроительной, электротехнической, в атомной энергетике, радиоэнергетике и др. Например, тепловые ограничения играют решающую роль при эксплуатации электрических машин и в радиоэлектронных устройствах. Их нагрев лежит в основе многообразных теплофизических процессов, которые могут угрожать жизнеспособности конструкции. При этом опасность повреждения связана не только с уровнем температуры, но и его распределением в пространстве и во времени. Отсюда вытекает необходимость в достоверной и подробной информации о распределении температур в данном агрегате. Охлаждение высокоэффективных реактивных и газотурбинных двигателей, «тепловой барьер» при больших скоростях движения летательных аппаратов, отвод малых тепловых потоков в атомных реакторах, генерация пара сверхвысоких параметров и другие проблемы новейшей техники расширили область практических приложений теории тепло- и массообмена. Во многих случаях знания процессов тепло- и массообмена являются определяющим фактором при выборе облика и параметров объекта. Современный уровень развития науки тепло- и массообмена предполагает использование системного подхода к их созданию в том числе и моделированию процессов переноса тепла и массообмена в различных средах и конструкциях. Правильное применение методов моделирования подразумевает увязку математического и физического (экспериментального) процессов областях на различных этапах проектирования и испытаний. Знакомство с учебной и монографической литературой показало, что усилиями зарубежных и отечественных ученых и педагогов созданы в течении XX столетия прекрасные пособия, в которых в строгой и доступной форме изложены основы теплообмена (теплопередача) и массообмена. Поэтому авторы считали естественным использовать имеющийся богатый научный и методический опыт. При этом возникает 3
обычный вопрос, зачем публиковать еще одно учебное пособие, когда можно воспользоваться имеющимися книгами. Прежде всего, для освоения предлагаемого материала студенту пришлось бы пользоваться примерно десятью книгами, которые к тому же малодоступны. Во-вторых, в данном учебном пособии использован и современный материал. В-третьих, отдельные разделы могут быть рекомендованы студентам других специальностей, которые сталкиваются с процессами теплообмена. В настоящем учебном пособии из-за ограниченного объема отсутствуют задачи и лабораторные работы, которые студенты решают на практических и лабораторных занятиях. По данным вопросам можно рекомендовать специальные учебные пособия. Надеемся, что наличие доступного для студентов учебного пособия поможет увеличить эффективность учебного процесса. 4
Обозначения t - температура (обычно по шкале Цельсия), °C Т — температура по шкале Кельвина, К Д/- перепад температур, К & - перегрев относительно некоторой температуры, К т - время, с п - нормаль к поверхности, Q - количество теплоты, Дж Ф — тепловой поток, Вт q - плотность теплового потока, Вт/м2 W - объемная плотность внутренних источников теплоты, Вт/м3 А - площадь (поверхности, сечения), м2 Я - теплопроводность материала, Вт/м.К а - коэффициент теплоотдачи, Вт/м2.К р - коэффициент объемного расширения, К'1; также температурный коэффициент изменения теплопроводности, К'1 Ср - удельная теплоемкость при постоянном давлении, Дж/кг.К 2/ v - кинематическая вязкость, м /с а - температуропроводность материала, м2/с g - ускорение силы тяжести, м/с2 и - скорость, м/с (р - угловой коэффициент F - тепловой коэффициент, К/Вт R - тепловое сопротивление, К/Вт 1 - расстояние, геометрический размер, м 6 - толщина, м V - объем, м3 U - периметр, м. Подстрочные индексы с - среда i,j - номер тела, поверхности пр - приведенное значение п — пластина ц - цилиндр ш - шар, сфера к - конвективная составляющая л - лучистая составляющая x,y,z - указывает соответствующее направление эф - эффективное значение 5
Оглавление Введение...................................................................4 Обозначения................................................................6 L Основные закономерности процессов тепло- и массообмена....................8 1.1. Процессы тепло- и массообмена в природе.............................8 1.2. Перенос тепловой энергии кондукцией. Основные понятия и определения..9 13. Перенос энергии конвекцией. Закон Ньютона-Рихмана...................14 1.4 Перенос энергии излучением..........................................17 1.5 Тепловое сопротивление и тепловой коэффициент.......................19 1.6. Тепловое сопротивление плоской, цилиндрической и сферической стенок..21 1.7. Составные стенки. Применение законов Кирхгофа......................22 1.8. Сложный теплообмен.................................................26 П. Коцдукция..............................................................28 2.1. Уравнение теплопроводности. Краевые условия........................28 2.1.1. Уравнение теплопроводности для анизотропного тела с источником энергии и переменными теплофизическими параметрами..............................28 2.1.2. Частные случаи уравнения Фурье.................................31 2.1.3. Краевые условия................................................33 22. Стационарное поле температур оболочек простейшей формы.............35 2.2.1. Плоская стенка.................................................35 2.2.2. Цилиндрическая стенка..........................................38 2.2.3. Шаровая стенка.................................................39 2.2.4. Тепловое сопротивление от сферы к неограниченному пространству.40 2.3. Стационарное поле температур тел с источниками тепла.................42 2.4. Нестационарный тепловой режим тела с равномерным полем температур....48 2.4.1. Дифференциальное уравнение процесса............................48 2.4.2. Интегрирование системы уравнений (2.61), (2.64)................50 2.4.3. Нагревание или охлаждение тела в среде с постоянной температурой.52 2.4.4. Нагревание или охлаждение в среде, температура которой изменяется во времени с постоянной скоростью........................................55 2.4.5. Температурный режим тела, помещённого в среду с гармонически меняющейся температурой...............................................57 2.4.6. Термическая инерция тела.......................................59 2.4.7. Внутренние источники энергии в теле............................61 6
2.5. Стационарное температурное попе стержней и пластин................61 2.5.1. Особенности теплообмена стержней и пластин. Дифференциальные уравнения...........................................................61 2.5.2. Дифференциальное уравнение теплопроводности для стержня.......62 2.5.3. Дифференциальное уравнение теплопроводности для пластины.....66 2.5.4. Д ифференциальное уравнение для диска........................68 2.6. Критерии неравномерности поля температур в теле...................69 2.7. Температурное попе стержня с источником тепла.....................73 2.8. Температурное поле пластины и диска...............................80 2.8.1. Обобщенное решение уравнения Бесселя..........................80 2.8.2. Круглое ребро постоянной толщины..............................82 2.8.3. Эффективность круглого ребра постоянной толщины..............84 2.9. Анализ ошибки измерения температуры..............................85 2.10. Принцип суперпозиции температурных полей........................88 Библиографический список................................................90 7
Основные закономерности процессов тепло- и массообмена 1.1. Процессы тепло- и массообмена в природе Тепломассообмен - раздел физики, в котором рассматриваются процессы переноса теплоты (энергии) и массы (вещества). Явления теплообмена связаны с необратимым переносом энергии из одной части пространства в другую и вызваны разностью температур. Явления массообмена связаны с необратимым перемещением вещества из одной части пространства в другую и вызваны разностью концентраций. Если оба явления сопутствуют друг другу и их приходится рассматривать во взаимосвязи, то соответствующие процессы называют тепломассообменом. Когда явления теплообмена и массообмена мало влияют друг на друга, их можно рассматривать порознь, иногда имеет место только какое-нибудь одно явление. Соответствующие процессы в этих случаях называют процессами тепло- и массообмена, теплообмена, массообмена. Перенос теплоты и вещества происходит благодаря молекулярному или конвективному процессу. Молекулярный перенос осуществляется микрочастицами (атомами, молекулами) в среде с неоднородным распределением температуры или концентрации. Конвективный перенос осуществляется макрообъёмами среды при их перемещении. Перенос энергии происходит с помощью механизма теплопроводности, конвекции и излучения. Теплопроводность - молекулярный перенос энергии в сплошной среде, вызванной разностью температур. Конвекция - перенос энергии совместно молекулярными и конвективными механизмами, вызванный также разностью температур. Излучение - перенос энергии электромагнитными волнами. Этот процесс обусловлен превращением внутренней энергии вещества в энергию излучения, переносом излучения и его поглощением веществом. Перенос массы происходит с помощью двух механизмов: диффузии и конвективного массообмена. Диффузия - молекулярный перенос вещества в среде, вызванный разностью концентраций (концентрационная диффузия), температур (термодиффузия) или давлений (бародиффузия). Конвективный массообмен - перенос массы, вызванный совместным действием конвективного переноса вещества и диффузии. Приведём примеры, связанные с переносом тепла и массы в природе. Теплообмен человека со средой (тепло- и массообмен); перенос тепла из жилища в окружающую среду и, наоборот, из среды в жилище (теплообмен); перенос энергии от Солнца к Земле (теплообмен); различные способы переработки вещества и продуктов — от приготовления 8
пищи до сложных производственных технологических процессов - также связаны с переносом энергии и вещества в пространстве (тепломассообмен). Такие процессы, как испарение, сушка, образование облаков, представляют собою целый комплекс явлений тепломассообмена, сопровождающихся фазовыми превращениями. В некоторых случаях имеют место какие-нибудь одни процессы. Например, процессы диффузии при производстве транзисторов, интегральных схем. В других случаях эти процессы настолько переплетаются, что трудно выделить один из них в качестве основного. Например, при взаимодействии мощных потоков лазерного излучения с веществом происходит нагревание последнего до температуры плавления, затем испарение; испарившееся вещество выбрасывается в окружающее пространство; дальнейшее поступление энергии приводит к ионизации паров, образованию плазмы и т.д. Благодаря сложности процессов тепломассообмена, целесообразно начать их изучение отдельно, т.е. рассмотреть порознь явления кондукции, конвекции, излучения, диффузии, конвективного массообмена, а затем обратить внимание на более глубокое изучение методов исследования таких явлений, в которых процессы могут встречаться в любой комбинации. 1.2. Перенос тепловой энергии кондукцией. Основные понятия и определения Две системы могут обмениваться теплом лишь в том случае, если они находятся при различных температурах, причём перенос тепла происходит в направлении от системы с более высокой к системе с более низкой температурой. Здесь не рассматривается вопрос о механизме переноса тепловой энергии кондукцией, конвекцией и излучением; основное внимание обращено на знакомство с феноменологическими законами переноса тепла. Количественные зависимости, связывающие разность температур систем с величиной теплового потока в случае кондукции, конвекции и излучения, имеют различный характер. Рассмотрим перенос тепла кондукцией (теплопроводностью) в различных телах и системах тел. Совокупность тел с различными теплофизическими параметрами (коэффициент теплопроводности, объёмная теплоёмкость) и явно выраженными границами раздела будем называть системой тел или неоднородным телом (рис. 1.1). 9
Рис. 1.1. Системы тел. Неоднородное тело (рис а) Каждая часть такой системы будет однородным телом. Однородные тела могут быть изотропными и анизотропными. В изотропном теле теплофизические параметры одинаковы во всех направлениях, в анизотропном - они различны в разных направлениях, но могут быть постоянными в выбранном направлении. Тепловое состояние тела или системы количественно характеризуется его температурным полем, т.е. совокупностью численных значений температуры в различных точках системы в данный момент времени. В том случае, когда температура во всех точках системы не изменяется с течением времени, поле температур называется стационарным, если же температуры в теле с течением времени изменяются, то поле температур называется нестационарным (рис. 1.2). € Рис. 1.2. Изменение температуры во времени (иллюстративный пример) Если температуры всех точек некоторого объёма равны между собою в любой момент времени, то это поле температур называют равномерным. Температурное поле в частном случае может зависеть только от одной координаты, тогда его называют одномерным. 10
Аналогичный смысл имеют термины двумерное и трёхмерное поле температур. Если тела находятся при различных температурах, то, как указывалось выше, возникает поток тепла, направленный от тела с более высокой температурой к телу с более низкой температурой. Для количественного описания этого процесса вводят два основных понятия: изотермическая поверхность и градиент температуры. Изотермической поверхностью называют геометрическое место точек, имеющих одинаковую температуру. По определению, через каждую точку внутри тела можно провести в данный момент времени только одну изотермическую поверхность. На рис. 1.3 линиями S] , $2 изображены следы на плоскости чертежа различных изотермических поверхностей в фиксированный момент времени. Рис. 1.3. Изотермические поверхности и линии теплового потока Вдоль изотермы температура не изменяется, в любом другом направлении — изменяется, причём в направлении нормали к изотермической поверхности наблюдается наибольшее изменение температуры на единицу длины. Возрастание температуры в направлении нормали к изотермической поверхности характеризуется отношением изменения температуры А/ между выбранными изотермами к расстоянию между ними по нормали А«. Предел этого отношения при устремлении Ан к нулю называется градиентом температуры (grad t). И
Градиент температуры есть вектор, направленный по нормали к изотермической поверхности в сторону возрастания температуры, т.е. lim __ = ц — = gradt К/м fl ° дп v 7 где п0 - единичный вектор, направленный по нормали в сторону возрастания температуры (рис. 1.3); — -производная температуры по дп направлению нормали. За положительное направление градиента температуры принимают направление, в котором температура возрастает. Линии, перпендикулярные к изотермическим поверхностям, называют линиями теплового потока, или, короче, линиями тока. На рис. 1.3 эти линии обозначены стрелками. В дальнейшем будем использовать следующие понятия: - количество теплоты, О. Дж; - мощность источников тепла или тепловой поток (если рассматривается перенос тепла через поверхность) Р - —, Вт; т - время, с; т - поток на единицу площади Р Вт называется удельным тепловым потоком или плотностью теплового потока. Основной закон, устанавливающий количественную связь между тепловым потоком и перепадом температур при кондуктивном теплообмене, называется законом Фурье по имени французского математика Ж.Фурье (1768 - 1830). Плотность теплового потока прямо пропорциональна градиенту температуры, т.е. — Qf q = -Xgradt = -AVt = -Лп0 —, (1.2) дп где Л - коэффициент пропорциональности, называемый коэффициентом теплопроводности; grad, V -математические символы, обозначающие градиент. 12
В математике положительное изменение функций направлено в сторону возрастания функции. В рассматриваемом случае тепловой поток направлен в сторону убывания температуры. Чтобы устранить противоречие между математическим и физическим определением положительного направления функции в зависимость (1.2) вводится знак «-». Соотношение (1.2) является обобщением опытных данных и лежит в основе всей современной теории теплопроводности. Выражение (1.2) можно представить в виде суммы составляющих градиента по осям декартовых координат , -dt - dt - dt gradt = i---F / — + к —, дх dy dz (1-3) одномерное поле где i, j,k -ортогональные между собою единичные векторы, направленные по координатным осям. n (dt Л Рассмотрим неизменное во времени — = 0 ^йт ) температур, изменяющееся лишь в одном направлении формулировка этих условий имеет вид x. Математическая dt dt п dt л dt л — =, — == 0; — ф 0; — - 0 йу & йх дт будет постоянной гт dt Пусть, кроме этого, градиент температур — дх величиной, т.е. температура линейно изменяется с координатой х, тогда dt t2 - — = ——- = const, где /у и t2- значения температур на поверхностях ху и х2 дх х2 -х, .причём tj> t2 , х2 > Xi. Величину плотности теплового потока q можно определить по формуле, которая следует из (1.2), Ат х2 -х. (1.4) Из (1.4) запишем выражение для коэффициента теплопроводности С *2 Ат tx-t2 т.е. коэффициент теплопроводности равен плотности теплового потока при перепаде температур на единице длины нормали, равном одному градусу. Размерность коэффициента теплопроводности в системе Си - Вт/м-К. 13
1.3. Перенос энергии конвекцией. Закон Ньютона-Рикмана Прежде всего, отметим, что капельные жидкости и газ здесь и в дальнейшем называем жидкостью. При конвективном теплообмене между твердым телом и жидкостью связь теплового потока Ф с температурами поверхности тела t и жидкости tc устанавливается формулой Ньютона-Рихмана P = a(t-tc)S, (1,5) где а - коэффициент теплообмена (теплоотдачи) конвекцией, А - площадь поверхности теплообмена. Если температура жидкой среды выше чем температура поверхности твердого тела, то Р = a(tc - t)S. Параметры, входящие в формулу (1.5), имеют следующие единицы измерения в системе СИ: [Ф] = Вт; [t] = [tj= К; [t - tj= К; [А] = м2; [а,] = Вт/м2К . Коэффициент теплообмена конвекцией численно характеризует тепловой поток, который рассеивает или воспринимает единица поверхности твердого тела при разности температур между твердым телом и средою в один Кельвин. Вся сложность теплообмена в формуле (1.5) концентрируется в одной величине - коэффициенте теплообмена а, который представляет собой сложную функцию большого числа параметров, существенно влияющих на процесс теплообмена. Так, для естественной конвекции a=/(t, tc ,Д Я,ср, v, a, g,0), (1.6) где /У - коэффициент объемного расширения, [1/К]', в частности, для газа п 11 - Z- ft = — =-----, для жидкостей берется из справочников; Т с Гс+273 А — коэффициент теплопроводности жидкости, Вт/(м-К); ср - удельная теплоемкость жидкости при постоянном давлении, Дж/(кг-К); v — коэффициент кинематической вязкости жидкости, м2/с; а - коэффициент температуропроводности жидкости, м2/с; g - ускорение силы тяжести, м/с2; Ф — символическое обозначение совокупности параметров, характеризующих форму, строение поверхности и ее размеры. К выводу о том, что в условиях естественной конвекции коэффициент теплообмена должен зависеть от перечисленных выше параметров, можно придти на основании анализа физической картины переноса. Предположим, что вертикально ориентированная стенка имеет 14
более высокую температуру, чем омывающий ее газ. В этих условиях объём газа, соприкасающегося со стенкой, получает от нее энергию благодаря столкновению молекул газа со стенкой; молекулы газа начинают двигаться с большей скоростью, первоначальный объем рассматриваемого количества газа увеличивается. Отсюда следует, что процесс переноса энергии должен быть связан с параметрами t, tc, Д ср . Увеличенный объем газа имеет меньшую плотность по сравнению с окружающими объемами, что приводит к движению этого объема вверх. Движение является результатом действия на объем газа нескольких сил: силы тяжести (появляется параметр g), силы Архимеда и сил внутреннего трения (последние связаны с коэффициентом кинематической вязкости у). Приведенные рассуждения показывают, что процессы движения объема газа связаны с его плотностью р, которая, однако, не фигурирует в перечне параметров. Плотность все-таки косвенно учитывается через соотношение между параметрами а, ср, р: aCpp = 'k (1-7) Так как в (1.6) есть параметры а, ср и 2, то тем самым косвенно учтено и влияние плотности р. Действительно, вместо параметра а можно было бы ввести р на основании соотношения (1.7). Такой сложный учет величины р имеет свои объяснения, которые будут приведены при более детальном знакомстве с процессом конвективного теплообмена. Перенос тепла от твердого тела к жидкости осуществляется не только благодаря движению макрообъемов жидкости, но также и за счет кондуктивного переноса тепла. В рассматриваемом примере кондуктивный перенос связан с движением и столкновением между собой молекул газа. Обычно конвективный и кондуктивный перенос тепла связаны с такой физической характеристикой, как теплопроводность газа, которая по указанной причине и фигурирует в функциональной зависимости (1.6). Конвективный перенос тепла связан с движением жидкости, а характер этого движения во многом зависит от геометрических характеристик тела и его ориентации в пространстве, чем объясняется присутствие в (1.6) параметра Ф. Часто движение жидкости или газа возле нагретой стенки является результатом не только температурного перепада t - tc, но и внешних воздействий (перепад давлений, создаваемый насосом, вентилятором или каким-либо другим путём), приводящих к перемещению жидкости или газа. Такой процесс называется вынужденной конвекцией. При малых скоростях движения наблюдается струйный характер течения (ламинарный режим), а при больших - неупорядоченно-вихревой (турбулентный режим). Этот факт был установлен О. Рейнольдсом с помощью следующих опытов. В протекающую по трубе воду вводилась тоненькая струйка 15
окрашенной жидкости. При скоростях течений, не превышающих некоторой критической UKp, отдельные части окрашенной струйки движутся только по направлению всего потока (рис. 1.4 а), но при U>UKp окрашенная струйка на небольшом расстоянии от входа в трубу размывалась и окрашивала всю воду. Из этого следует, что продольное движение частиц сопровождается движением, в котором частицы приобретают также значительные радиальные составляющие скорости (рис. 1.4 в). Напомним, что при переходе ламинарного режима в турбулентный сопротивление трения в трубе резко возрастает. Существует также весьма неустойчивый переходный режим движения жидкости (рис. 1.4 б). Рис. 1.4. Опытные результаты О. Рейнольдса а) ламинарное течение б) переходный режим течения в)турбулентное течение При вынужденной конвекции коэффициент теплоотдачи, в основном, зависит от следующих параметров а = а(С, v, 2, а, Ф) (1-8) Здесь через U обозначена скорость потока, омывающего поверхность тела. Молекулярный перенос тепла так же, как и в случае свободной конвекции, зависит от коэффициента теплопроводности жидкости и количества переносимой жидкости, которое, в свою очередь, связано с ее скоростью U. Силы противодействия движению потока жидкости связаны с ее вязкостью v и совокупностью геометрических параметров Ф. Меньшее влияние при вынужденной конвекции оказывают температуры поверхности и жидкости, поэтому в (1.8) температуры в явном виде не фигурируют. 16
Рис. 1.5. Характер теплообмена при наличии вынужденной и естественной конвекции Возможен конвективный перенос, в котором и вынужденная, и естественная конвекции вносят соизмеримые доли в общий процесс (рис. 1.5). Тогда коэффициент теплоотдачи зависит от всех параметров, фигурирующих в зависимостях (1.6) и (1.8). Приведём порядок величин коэффициентов конвективной теплоотдачи: - естественная конвекция: - для газов а = 2-10, Вт/м2 К; - для масла а = 200 - 300, Вт/м2 К; - для воды а = 400 - 600, Вт/м2 -К; - вынужденная конвекция: - для газов а = 10 - 100, Вт/м2 -К; - для масла а = 300 - 1000, Вт/м2 -К; - для воды а = 1000 - 10.000, Вт/м2 К, в зависимости от режима течения; - кипение воды - а = 20.000 - 40.000, Вт/м2 -К, однако для фреонов эта величина на порядок меньше; - конденсация паров воды - а = 30.000 - 50.000, Вт/м2 -К,. Данные значения являются сугубо ориентировочными и не могут служить основанием для расчетов. 1.4 Перенос энергии излучением В разделе 4 будет обоснована следующая зависимость между потоком энергии Фу, излучаемым через прозрачную среду с поверхности тела i с площадью А{, к поверхности тела j с площадью Aj (рис. 1.6): (1.9) 17
где Т., Tj - абсолютные температуры поверхностей тел, К; С=5,67 Вт/м2К4 - коэффициент излучения абсолютно черного тела. Коэффициент enpij носит название приведенной степени черноты и учитывает физические свойства поверхностей, их взаимное расположение, форму и размеры. Параметр ф - коэффициент облученности j -го тела i -ым - учитывает долю потока от тела i, попадающего на тело]. Представим (1.9) в форме, аналогичной закону Ньютона-Рихмана (1.5) для конвективного теплообмена: = <Мб -Q4 (1.Ю) где a3ij - коэффициент теплообмена излучением между поверхностями i и j. Если тело i находится в неограниченной среде, то tj равно температуре среды tc. Рис. 1.6. К расчету теплообмена излучением между двумя телами В формуле (1.10) вся сложность процесса теплообмена излучением сконцентрирована в одной величине - aMj, структуру которой нетрудно определить, приравнивая правые части (1.9) и (1.10) (1-Й) = 5.67-10-8 У tj (1-12) Если температуры и близки, то /О,,^) можно рассчитывать по приближенной формуле 18
3 ж,?.) = 0.227 — ; Т = 0.5(7] + Г) J v1 з' |i qq ’ v 1 j ' (1-13) 1.5 Тепловое сопротивление и тепловой коэффициент Обозначим через Ф = qA = — тепловой поток, проходящий через т поверхность Л, и преобразуем уравнение (1.4): ti —12 = ФР; (1.14) (1.15) ЛА где параметр F - тепловой коэффициент. Структура параметра F в выражении (1.15) справедлива для плоского тела с одномерным стационарным полем температур. В общем случае, когда градиент температуры зависит от координат х, у, z, выражение для F можно получить из (1.2) и (1.14). Представим (1.2) в виде 5г =—— дп и найдем разность температур АЛ между изотермами ti и t2 , расположенными на расстоянии 7/ и 12 от начала отсчета (рис. 1.3): t1-t2=-(t2-t1) = (1.16) Сравнивая последнее выражение с (1.14), находим общее выражение для F: 1 } Ф(Г)дп где дп - элемент длины пути теплового потока; А (I) - аналитическое выражение площади изотермической поверхности на расстоянии I от начала отсчета; Ф(1) - величина теплового потока через изотермическую поверхность A(l)’t Ф1 - величина теплового потока через изотермическую поверхность А(Ц); Ц и 12 - расстояния от начала отсчета изотермических поверхностей А] и А2 (рис. 1.3). Если на пути теплового потока между изотермами Aj и А2 отсутствуют источники или стоки энергии как в теле, так и на его 19
границах, то поток Ф в этой области не меняет своей величины, т.е. Ф = Ф1 и (1.17) приобретает более простой вид: WA(l) (1.18) При неизменном значении коэффициента теплопроводности: 17 1 дп F=— ------. Я J Л(/) (1.19) Значительное число задач стационарной теплопроводности по существу состоит в определении теплового коэффициента F. Рассмотрим аналогию между переносом тепла через твердое тело и протеканием электрического тока через проводник. Сила электрического тока I и разность потенциалов (Ui - U2) связаны между собой законом Ома (1.20) где R3 - электрическое сопротивление между эквипотенциальными поверхностями А} и Аг. Сопоставляя выражения (1.14) и (1.20), видим, что аналогом температуры является электрический потенциал, теплового потока - электрический ток, а теплового коэффициента - электрическое сопротивление. Если между эквипотенциальными поверхностями отсутствуют источники тока, границы области имеют абсолютную изоляцию, а удельная электрическая проводимость и является постоянной величиной, то электрическое сопротивление также можно представить в форме, аналогичной (1.19): (1.21) Приведенная аналогия объясняет, почему параметр F часто называют тепловым сопротивлением. Укажем границы применимости рассмотренной здесь аналогии. Зависимость (1.20) между током и разностью потенциалов в проводнике справедлива при указанных выше условиях, а именно: в проводнике отсутствуют стоки электрического тока (например, утечка с боковой поверхности) и источники, создающие дополнительный ток внутри 20
проводника. Эти условия можно обобщить следующим образом: если между двумя изо потенциальными поверхностями (изотермическими или эквипотенциальными поверхностями) отсутствуют стоки и источники энергии, то поток энергии (тепловой поток, электрический ток) остаётся неизменным; в этом случае справедлива аналогия между процессами переноса тепла и электричества в цепях с распределёнными параметрами. Аналогом электрического сопротивления Яэ является тепловой коэффициент F который будем обозначать в этом случае через R и называть тепловым сопротивлением. Величины R3 и R определяются из одинаковых по структуре выражений (1.19) и (1.21). Эту аналогию можно продолжить и применить разработанные в электротехнике приёмы, основанные на законах Кирхгофа для электрических цепей, для расчета теплового сопротивления сложной цепи. Если сформулированное выше условие не выполняется, то электрическое сопротивление R3 не является аналогом теплового коэффициента F. 1.6. Тепловое сопротивление плоской, цилиндрической и сферической стенок На рис. 1.7 изображены однородные стенки различной конфигурации, поверхности которых х=1] и х=12 являются изотермическими с температурами tj и t2 , а торцы плоской и цилиндрической стенок являются адиабатическими; внутренние источники тепла в стенке отсутствуют; коэффициент теплопроводности материала - Л. Найдем выражение для стационарного теплового потока Ф через эти стенки. Воспользуемся зависимостями (1.14) и (1.15), связывающими разность температур (ti - ti) с тепловым потоком Ф, а значение теплового коэффициента F найдем для каждого конкретного случая с помощью выражения (1.19). Элемент длины dn пути теплового потока для плоской (п), цилиндрической (ц) и шаровой (ш) стенок равен dn = dx, а аналитические выражения А(х) изотермических поверхностей имеют вид Ап =L] L2 , Ац =2яхЬц, Аш=4т1х2, где Li и Ь2-длина и ширина плоской стенки; Ьц— длина цилиндрической стенки. Поскольку по условиям задачи между изотермическими поверхностями отсутствуют источники и стоки энергии, коэффициент F имеет смысл теплового сопротивления, которое для плоской, цилиндрической и шаровой стенок обозначим через Rn, R4, Rw. Учитывая неизменность потока тепла 0(x)=const и подставляя значения dn и А(х) в выражение (1.19), получим 21
г дх _ $ J ALL, ~ ALL, ~ ЛА ’ г дх 1 ’г clx 1 , /2 • Л2лхЬа 2лАЬц х 2яйХц 12 г дх _ 1 ‘ г с/х 1 ' Л4лх2 4л2’ х2 а) б) (1.22) в) 11 4л2 I. L } где 5 = I2 -1]- толщина плоской стенки. Рис. 1.7. К расчёту теплового сопротивления стенок различной конфигурации 1.7. Составные стенки. Применение законов Кирхгофа Рассмотрим теперь последовательно составленную плоскую стенку, состоящую из п разнородных, ориентированных перпендикулярно потоку тепла слоев, толщины и коэффициенты теплопроводности которых - с)г и 2г .Температура наружных поверхностей стенок равны t} и tn+1 (см. рис. 1.8 а). Изотермическими поверхностями в этом случае являются плоскости, параллельные поверхности стенок. Между изотермическими поверхностями отсутствуют стоки и источники энергии и тепловой поток, не изменяясь, проходит через все стенки. Следовательно, каждой стенке 22
Рис. 1.8. Последовательное соединение стенок можно приписать тепловое сопротивление Яг ,все тепловые сопротивления, как это видно согласно рис. 1.8 б, соединены последовательно, т.е. суммарное тепловое сопротивление: 8, 1 8> д. л” (1-23) На основании (1.14) между температурами и t„+1 и тепловым потоком Ф справедлива зависимость =™ = у2д-- (1.24) Проведём аналогичные рассуждения для тепловых сопротивлений последовательно составленных цилиндрической и сферической неоднородных стенок, состоящих из п неоднородных слоёв, расположенных перпендикулярно потоку тепла: in А. 1 1 1 1 г 1 _ 1 ' (1-25) (1-26) где Ц - радиус z-ro цилиндрического или сферического слоя; Ьц - длина цилиндра. Найдём выражение для теплового сопротивления параллельно- составной стенки, образованной системой плоских стенок, расположенных параллельно потоку тепла (рис. 1.9 а). 23
Рис. 1.9. Параллельное соединение стенок Ограничивающие поверхности с температурами ti и Ь являются равноотстоящими плоскостями, ориентированными перпендикулярно тепловому потоку. Для определения теплового сопротивления параллельно-составной стенки введем допущение, что разнородные стенки отделены одна от другой бесконечно тонкими адиабатическими (не проводящими тепла) прослойками. Тогда температурное поле в каждой стенке становится одномерным, и ее тепловое сопротивление может быть рассчитано по формуле (1.22 а). Тепловые сопротивления отдельных стенок соединены параллельно (см. рис. 1.9 в), поэтому их общее тепловое сопротивление потоку тепла равно (1-27) Принятое допущение о наличии адиабатических прослоек позволило существенно упростить вывод выражения для теплового сопротивления параллельно-составной плоской стенки. Нетрудно получить выражение и для теплового сопротивления параллельно-составной цилиндрической стенки (рис. 1.9 б). Полагая, что разнородные стенки отделены одна от другой адиабатическими прослойками, из формулы (1.22 б) для Кц и схемы соединений сопротивлений (рис. 1.9 в) получим (1-28) 24
где Lt- длина z-ro участка цилиндрической стенки. Сложное соединение стенок. При расчете тепловых сопротивлений различных технических устройств чаще встречаются случаи сложного соединения, в которых тела в системе соединены как последовательно, так и параллельно. Примером такого соединения может служить крепление некоторых деталей к шасси с помощью болта (рис. 1.10 а). В этом случае корпус детали электрически изолирован от шасси шайбами 2 и 4 из электроизоляционных материалов, которые одновременно являются и теплоизоляцией. Рис. 1.10. Пример сложного соединения стенок Болт 2 отделен от шасси воздушной прослойкой 8, поэтому теплообмен через прослойку между болтом и шасси практически отсутствует. Тепло от детали к шасси поступает двумя потоками Ф/ и Ф2 Поток Ф/ идет непосредственно через изоляцию 2, поток Ф2 - более сложный путем: от детали 7 через болт 1, гайку 6, шайбу 5 и слой изоляции 4. На рис. 1.10,а пути тепловых потоков обозначены стрелками. В данном случае тепловые потоки Ф/ и Ф2 движутся параллельно, преодолевая тепловое сопротивление R2 участка 2 и тепловое сопротивление нескольких последовательно соединенных элементов. На рис. 1.10,6 представлена общая схема соединения тепловых сопротивлений для рассматриваемого случая. Если предположить, что тепловой поток, проходящий через отдельные участки системы, не рассеивается, то результирующее тепловое сопротивление R крепёжного соединения найдем на основании закона Кирхгофа, а именно: 25
„ /?'/?" nl п 82 R =--------, R = R2 =----------, R'+R" 2 Я2Л2 fl" = + + ^5 + ^1^1 ^4'^4 ^5^5 ^“6'^6 Параметры 8^ и А6 в данном примере довольно условны, приближенные методы их определения устанавливаются в каждом конкретном случае. 1.8. Сложный теплообмен Зависимости между тепловым потоком и разностью температур (ti~ ti) можно представить для кондуктивного, конвективного и лучистого механизмов переноса в единой форме: 6 t2 „ (1-29) где т - индекс, характеризующий механизм переноса (кондуктивный (т=Т), конвективный (т=К) и лучистый (т=Л))’, Фпт - тепловой поток между изотермическими поверхностями 1 и 2 для механизма m; tjia. t2- температуры изотермических поверхностей 1 и 2; В-12т - тепловое сопротивление потоку между изотермическими поверхностями 1 и 2 для механизма переноса т. Если сопоставить формулу (1.29) с формулами для различных механизмов переноса, например, с формулами (1.5), (1.10), (1.14), то получаем следующие выражения для теплового сопротивления: кондуктивный перенос R12T =R, где R- тепловое сопротивление, определяемое по формуле (1.16), а для частных случаев плоской, цилиндрической и шаровой стенок - по формулам (1.22); конвективный перенос (130) где ак - конвективный коэффициент теплообмена; лучистый перенос Я,2Л =(«лЛ)-', (1-3 где ал - коэффициент теплообмена излучением, структура которого дана формулой (1.11). Если все три механизма переноса тепла существуют одновременно, 26
то тепловой поток Ф/2 между изотермическими поверхностями 1 и 2 равен Ф12 = Ф12К + Ф12Т + Ф12Л', (1-32) Заметим, что все приведенные здесь зависимости получены в предположении отсутствия стоков или источников энергии между изотермическими поверхностями 1 и 2. Расширим понятие теплопередачи от одной жидкой среды к другой через твердую стенку. С этой целью рассмотрим перенос тепла в трубчатом теплообменнике, где жидкость с температурой tci, протекающая по трубе длиной L, нагревает жидкость с температурой tc2, омывающую наружную поверхность трубы (см. рис. 1.11). Рис.1.11. К расчету теплового сопротивления трубчатого теплообменника На основании закона Кирхгофа тепловой поток от среды с температурой tcl проходит через стенку к жидкости с температурой tc2, преодолевая следующие последовательно соединенные тепловые сопротивления: среда-стенка (Rci), стенка (Rn) и стенка-среда (Rc2); общее сопротивление потоку R равно: R Rc} + Rj2 + Rc2, (1-33) где R]2 - тепловое сопротивление простой или многосоставной стенки, определяется по формулам (1.22), (1.25), a Rc/, Rc2 равны: RC1 =—, Re2 а{Ах а2А2 (1.34) 27
Здесь ос/ и а.2 - полные коэффициенты теплообмена (конвективно- кондуктивный и лучистый) между средой и поверхностями стенок; А; =2яЫ1 и А2 =2tiL12 - площади поверхностей 1 и 2. Из формул (1.33), (1.34) следует, что 11 i,z2 1 ---- ----4- Jn 2 н------ Л / а212 (1-35) Анализ формулы (1.35) приводит к выводу, что тепловой поток через отдельные сопротивления для рассматриваемого цилиндрического тела уменьшается непропорционально увеличению толщины изоляции. При росте /2и неизменном 1} тепловое сопротивление —-—In— увеличивается, а 2я£2 /, тепловое сопротивление —?—уменьшается. Такого рода двойной эффект '2.тА220и2 означает, что для цилиндрической стенки существует определенный критический радиус 1кр, при котором потеря тепла является максимальной. Дифференцируя значение R из (1.35) по 12 и приравнивая производную нулю, найдем выражение для 1кр - критического радиуса изоляции: 4Р » (1.3< при котором будет наименьшее сопротивление R потоку. II. Кондукция 2.1. Уравнение теплопроводности. Краевые условия 2.1.1. Уравнение теплопроводности для анизотропного тела с источником энергии и переменными теплофизическими параметрами Из обширного класса анизотропных тел рассмотрим такие, в которых тело имеет различные коэффициенты теплопроводности Лх, /у, 2z в трех взаимно перпендикулярных направлениях, принятых за оси координат х, у, z. Если принять это ограничение, то останется справедливой запись закона Фурье в форме (1.2). Предположим, что в объеме AV= AxAyAz могут находиться источники тепла, удельная мощность которых в момент времени т равна W(x, у, z, т), Вт/м3 (рис.2.1). Рассмотрим изменение теплового состояния объема AV за промежуток времени dr. При этом будем базироваться на двух законах: законе сохранение энергии и законе Фурье. 28
Количество тепла, поступающего в объем слева через поверхность ЛуЛг за время dr, обозначим dQx, а количество тепла, выходящее из объема AV через правую поверхность dydz, обозначим dQ'x. Аналогично введем обозначения dQy , dQ’y , dQz , dQ’z. Пусть за время dz температура t( х, у, z, т ) объема А V повысилась на dt градусов, тогда количество тепла, поглощенное объемом, равно cpAVdt; действие источников тепла вызовет дополнительное поступление энергии, равное WA Vdz. На основании закона сохранения энергии: dQx + dQy + dQz + = dQ'x+dQ'y +dQ'z+cpA Vdt, (2.1) Рис. 2.1. К выводу дифференциального уравнения теплопроводности На основании закона Фурье имеем [ ] dQx = qjXy&zdz = —2хДуДг — I dz; (бх А dQ'x =q'xb.yhzdz = 2.v/\v/Vl — dz, (2.2) где q» q'x - удельные тепловые потоки в направлении оси х. Аналогично записываем выражения для направлений у и z: dt dO'v =q'vAxAzdz = AxAz — dz; /у + &у (2.3) 29
dQ = q AxAydr = —Л AxAy — dr, ( dQ'z =q'-/\xA.ydT = -ДАхАу — dr. Подставим (2.2) и (2.3) в равенство (2.1): 1 Л Л ( I J 1 Л Л I I J 1 * A I 1 J ТТЛЛ Л * J ДДуАг — dr —Л AxAz — dr - Д Ах Ду I — d‘F + Wax&y&zdT = {dx)r у \ду ) ':5zj . I qI а ЛхAyAz — dr — Л ArAz (6xLAv у dt А <2Д) [ ot । dr ~ ДАхАу — dv + cp AxAyAz dt. \ rfe / / у+Ду 7z+Az । ot Заметим, что t есть функция х, у, z, т, а следовательно, произведение Д— I di есть функция х, у, z, т; здесь i = х, у, z. Обозначим последнее произведение через Фг, т.е. , dt] ~ . dt) , =Ф-> А- = di л \ di /1 X S 1У Д1 (2-5) Используя обозначение (2.5), перепишем равенство (2.4): ф„д.-А . ф^-фу, I I I Г Г t-. f-j • Дх Av Az дт (2.6) Переходя в (2.6) к пределу при Ai—^O и используя определение производной Получим на основании определения (2.5) д дх (2.7) Д'->о Qi Qi Это и есть дифференциальное уравнение теплопроводности для рассматриваемого класса анизотропных тел с источниками тепла и теплофизическими параметрами, зависящими от температуры. Рассмотрим класс анизотропных тел в цилиндрической системе координат. Пусть коэффициенты теплопроводности в направлениях х и z равны Лх иД. Тело ограничено плоскостями, перпендикулярными оси z. 30
Для таких тел уравнение теплопроводности с учётом осевой симметрии имеет вид I х L* vV I I I XL* I I У У О JC-X • х дхI дх) dz\ dz) дт (2-8) Замечание. Как уже отмечалось, в общем случае для анизотропной среды направление вектора теплового потока в какой-либо точке, вообще говоря, не совпадает с направлением нормали к изотерме, проходящей через эту точку, а закон Фурье записывается в ином виде: где Лп -коэффициенты теплопроводности (компоненты тензора второго ранга). Общая теория переноса тепла через анизотропные среды изложена в [10]. 2.1.2. Частные случаи уравнения Фурье а) Пусть коэффициенты теплопроводности 2Л и л- не зависят от координат, но зависят от направления; тогда уравнение (2.7) примет вид б) Если Лх = ку = =2, то тело изотропное и уравнение (2.9) становится более простым: о л ’ 'Л-1’ л а дт ср Здесь через а обозначен коэффициент температуропроводности материала. Разъясним физический смысл этого коэффициента. Температуропроводность характеризует способность материала , „ .. dt повышать свою температуру с большей или меньшей скоростью — под действием притекающего тепла. Если в каком-либо слое материала происходит повышение температуры, то скорость этого повышения будет, 31
прежде всего, определяться тем, какое количество тепла передается этому слою в единицу времени от соседних с ним слоев. Это последнее количество тепла пропорционально теплопроводности, а, следовательно, скорость повышения температуры пропорциональна 2, С другой стороны, чем больше теплоемкость единицы объема материала cp=Cvoi (Cvoi - объемная теплоемкость), тем меньше будет повышаться температура рассматриваемого слоя; следовательно, скорость повышения температуры слоя обратно пропорциональна Cvoi. Отсюда приходим к заключению, что свойство, характерное для поведения материала при изменяющейся температуре, должно определяться параметром а = .Единица С vol СР измерения температуропроводности в системе Си — м2/с. Покажем, что путем тождественных преобразований уравнение (2.9) можно привести к форме, аналогичной (2.10). Введем новые координаты х'=хЛ-у,=у& z'=z& (2-и) где 2 - так называемая базовая теплопроводность, выбор которой произволен; обычно за Л принимают одно из трех значений коэффициентов теплопроводностей Лх, Луи 2г. Подставим новые значения координат (2.11) в уравнение (2.9), для этого предварительно проделаем следующие выкладки: dt _ dt dx' _ Гл~ dt d2t 2 d2t дх дх' дх у Лх дх’’ дх2 Лх д(х')2 ’ ( • ) Аналогично d2t _ Л d2t _ d2t _ Л d2t ду1 4 3(У)2 ’ dz2 Лх d(z')2 ’ Уравнение (2.9) принимает после преобразования вид d2t d2t d2t W 1 dt .. (2-12) аналогичный по форме уравнению (2.10). Преобразованию типа (2.11) подвергаются и граничные условия, так как в них всегда входят координаты и геометрические параметры. Итак, преобразования (2.11) позволяют свести решение задачи для анизотропных твёрдых тел к решению соответствующих задач для изотропных тел. Напомним, что этот вывод справедлив только для рассматриваемого здесь 32
класса анизотропных тел, определение которого указано в начале раздела. Для цилиндрической системы координат уравнение (2.8) с помощью преобразований (2.13) приводится к виду, характерному для изотропных тел d2t 1 dt d2t W _ 1 dt c>(x')2 + / дх' + d(z')2 + Л ~ адт' (2.14) в) Если в теле отсутствуют источники энергии W=0, то уравнение (2.10) переходит в следующее „2 «V7 = — дт (2-15) Большинство исследований по аналитической теории теплопроводности связано с решением этого уравнения. г) В стационарных условиях —=0, и уравнения (2.7), (2.9), (2.10) дт записываются в более простой форме: , d2t . d2t 3 d2t а? 'а/ -а? ; ™ = 0. Л (2.16) (2.17) (2.18) 2.1.3. Краевые условия Задача теории теплопроводности состоит в определении поля температур в теле в данный момент времени. Для решения этой задачи, кроме дифференциального уравнения, необходимо знать поле температур для какого-нибудь предшествующего момента времени (начальное условие), а также форму тела и закон взаимодействия между окружающей средой и поверхностью тела (граничные условия). Начальное и граничные условия в совокупности называются краевыми условиями. Начальное условие определяется заданием закона распределения температуры в теле в начальный момент времени, т.е. 33
t( x, у, z, 0) = ^(х, у, z). (2.19) Граничные условия можно представить в различной форме в зависимости от характера явлений, протекающих на границе тела. а) Например, может быть задано распределение температуры на поверхности тела в любой момент времени (задача Дирихле, или условие первого рода) Ь(т) =f(r), (2.20) где ts(z) — температура на поверхности тела в момент времени т. б) Задано распределение плотности теплового потока в любой момент времени (задача Неймана, или условие второго рода) qs(r) =f(r), (2.21) где qs(r) - плотность теплового потока на поверхности тела. в) Задан закон теплообмена между поверхностью тела и окружающей жидкой или газообразной средой (условие третьего рода). Математически выражение этого закона может быть получено на основании следующих рассуждений. По закону Ньютона удельный тепловой поток на границе тело-среда равен: qs(r) = a(ts -tc), (2.22) где а - коэффициент теплообмена. По закону Фурье к поверхности изнутри тела подходит удельный поток СQ ["] q's=-A\ — , (2.23) где п - внешняя нормаль к поверхности тела. Если на границе тело-среда отсутствуют стоки или источники энергии, то qs(r) = q's(r) и граничное условие принимает вид: Л (д Л ---- 4" СС \dn)s -Zt) = 0 (2.24) Если на границе тела имеется источник энергии, производительность Вт которого равна q, —- , то граничное условие нетрудно получить, м основываясь на законе сохранения энергии 34
q's+q= qs (2.25) В последнее равенство следует вместо q's и qs подставить их выражение (2.22) и (2.23). г) Условия теплообмена на границе двух твёрдых тел (условия четвертого рода): - отсутствие температурного скачка на границе (при отсутствии контактного теплового сопротивления) <,|s = <2|s . (2.26) где 5 - символическое обозначение границы. - если на границе двух тел имеется источник (или сток) энергии, поверхностная мощность которого q, то на основании закона сохранения энергии (2.27) где п - нормаль к поверхности раздела тел. 2.2. Стационарное поле температур оболочек простейшей формы 2.2.1. Плоская стенка Полагая в уравнении (2.16) W=0 и принимая во внимание, что б? _ dt _ бу oz получим d ( , dt} _ = 0 и зададим граничные условия ?(0) = , /(/) = t2 Будем считать, что 2 =f(t), тогда, интегрируя (2.28), получим f(t)dt = C\dx, 35
Постоянная Ci имеет простой физический смысл. Сравним выражение , dt „ . „ dt 2— = С, с законом Фурье Л— = -q, dx dx тогда Ci = - q, т.е. Ci численно равна удельному тепловому потоку через стенку, взятому с обратным знаком. fE Следовательно, f(t)dt = qx. t Рассмотрим частные случаи. Пусть l=const. t —t fAdt = A(tl-t) = qx; q = AJ—. (2.30) При x = /, t=t2 и q = /J} Пусть 2=2o(7 +/?(). t2 _ /2 t-t ?x = A>f(i + ^=A)ai-0 + A£LyL = A)0i-0(i + ^I1y-). (2.31) Если p>0, то распределение температур в плоской стенке изобразится параболой с выпуклостью, обращённой в сторону (+/) (рис. 2.2). Если fl<0, то распределение температур в плоской стенке будет иметь вогнутый вид, представленной на рис. 2.2 Рис. 2.2. Характер изменения температуры плоской стенки 36
При X=l, t=t2 И q —уЦ1 + /7г)01 -t2), где t = 0,5(t! + t2) Обозначим A = Ao (1 + fit}, тогда _ A(tx -t2) (2.32) Тепловое сопротивление плоской стенки равно j? = ~ = " Sq AS ’ где 5 - площадь плоской изотермической поверхности. Отклонение распределения температур от прямой линии можно найти из формулы (2.31). Полагая в этой формуле х = I, получим ql ~ Afitx -tfi ! + /?( Zi +t2 2 Отношение равно 1 + Д 2 Ошибки в оценке расстояния х, на котором наблюдается заданная температура t, дана вторым множителем (дробью) в правой части последней формулы. Пример 1. Для большинства технических материалов-диэлектриков можно принять: А=А0(1 +fit), fi>0. Например, теплопроводность кирпича следующим образом зависит от температуры Ао = 1,0 + 10 м-К Пусть I = 0,4 м, /;=900°С, /2=100°С , тогда 37
Л =1,0 + 500 1000 = 1,5-?^- м К 900-100 Вт д = 1,5 = 3000-^ 0,4 м Если бы пренебрегли зависимостью X=k(t), то в определение теплового потока внесли бы ошибку в 1,5 раза. Пример 2. Пробковая стенка: коэффициент теплопроводности изменяется с температурой по закону 2 = 0,04(1 + 4,5-10’4), 7=0,1 м, ti=30 °C, t2= -20 °C, +tl =5°C 2 = 0,04(1 + 4,5 -10 1 • 5) = 0,0409 2 м • К g=0,0409 30"("2Q)=20,45g| 0,1 M2 Если не учитывать зависимость k=k(t), то ошибка составила бы 2,25%. Отсюда видно, что, смотря по обстоятельствам, зависимостью Л от t можно иногда пренебрегать, особенно если рассматривается неширокий диапазон температур. 2.2.2. Цилиндрическая стенка R=o Из уравнения (2.8) для стационарного режима dz неограниченного цилиндра получаем уравнение 1 d dt. _ ---(Лх—) = 0. х dx dx Проинтегрируем это уравнение: Лх— = С, dx без источников тепла {№=()), (2.33) (2.34) Найдём значение постоянной интегрирования С; из следующих рассуждений. Выражение для теплового потока через цилиндрическую поверхность имеет вид: Р = -Л—2яхЬ. dx (2-35) 38
Сравнивая (2.34) и (2.35), находим с =___ 1 Проинтегрируем уравнение (2.34) от изотермической поверхности радиуса X] до поверхности радиуса х, температуры этих поверхностей соответственно равны ti и t: f Adt = f — = С, in—- -С, In—— ' ' Л' X X] Если X=Xq(1 +flt), то согласно (2.31) г Л 4“ и температура t связана с тепловым потоком Р выражением: -Ain—= 2о(?1-О(1 + ^^). (2.36) 2л£ X] 2 ' В частности, при (3 = 0, Pin — t^t =__(Ш) 1 2л£Л0 Найдём выражение для теплового сопротивления цилиндрической стенки 6 G Р Для этого положим х=х2 и t=t2 и введём средний коэффициент теплопроводности (2.32), тогда 1пх/ Al 2я£20 Итак, выражение для получено иным методом, чем в разделе 1.5, с использованием решения уравнения теплопроводности. 2.2.3. Шаровая стенка Тепловой поток через сферическую изотермическую поверхность радиуса х равен 39
Р = -Л—4ях2, dx откуда Принимая во внимание выражение (2.31) для правого интеграла и вычисляя интеграл в левой части последнего равенства, получим р ( j i^i t +t --------+ (2.38) 4л-|<х1 xj 2 В частности, при Р=0 получаем Р 1 1 t =----(----- 4я20 %] х (2.39) Положим теперь в (2.38) х=Х2 vit=t2^ найдём выражение для потока через шаровую стенку _ 4л‘2(/1 -Т2) т"2"т.......... (2.40) где Л = Л)(1 + pt}, f = 0,5^+t,) 2.2.4. Тепловое сопротивление от сферы к неограниченному пространству Пусть в неограниченном пространстве расположена сфера радиуса X], поверхность которой поддерживается при температуре tf, на большом удалении от сферы температура пространства равна t2 (рис. 2.3 а). Определить сопротивление R тепловому потоку Р от сферы к неограниченному пространству. По определению Р следовательно, задача сводится к нахождению величины теплового потока Р, который на основании закона Фурье равен 40
а) б) Рис. 2.3. К определению теплового сопротивления от сферы к неограниченному пространству dt dx 4лх2. (2.41) Итак, необходимо найти температурное поле в неограниченном пространстве, окружающем сферу. Дифференциальное уравнение теплопроводности без источников тепла в установившемся режиме имеет в сферических координатах вид d2t 2 dt Л 1 d1 . . п — +-------= о или---------(xt) = Q, dx х dx х dx откуда ^ = C., xt = C{x + C2 (2.42) dx Запишем условия на границах: t(x}) = tj, t(x>) = t2 (2.43) Из (2.42) и (2.43) следует, что Cj = t2, C2=(ti~t2)xi (2.44) Итак, t = t2 +(/j -r2)— x Из (2.41) и (2.42) находим величину Р: P = 4W,-',)x„ (2.45) 4лЛх{ 41
Найдём выражение для теплового сопротивления полупространства с адиабатической границей, в поверхность которого вдавлена полусфера с постоянной температурой поверхности (рис. 2.3 б). Этот случай отличается от предыдущего только тем, что весь поток Р направлен вглубь полупространства. В предыдущей задаче этот поток поровну распределился в верхнюю и нижнюю части полупространства вследствие симметричности тела. Следовательно, для полупространства тепловое сопротивление должно быть вдвое больше, чем для пространства, т.е. R = _Л_ 2/гЯх (2.46) 2.3. Стационарное поле температур тел с источниками тепла Введём следующие ограничения: тела рассматриваются с постоянным тепловыделением по объему; на границах тел теплообмен подчиняется закону Ньютона (условия третьего рода); теплопроводность и теплоотдача не зависят от температуры. Неограниченная пластина. (^t Из уравнения (2.9) при —F г = 0 следует, что ду dz Я^-± + ^ = 0, 0<х<8. (2.47) dx Запишем условия на границах (рис. 2.4 а) . (dt । . ч —Л] —— — OI=o 5 \dxjx=o .(dt\ , ч — Л] -у- “ ^2 0 ~~ С ' \dxjx=s (2.48) (2.49) Решение d ( dt dx\dx dt dx w t = ~—x +C]% + C2. (2.50) W Л ’ Постоянные интегрирования Cj и C2найдём из уравнений (2.48) - (2.49): 1 ВЛ 1 2Л ЛВ 1 £| " ' , W8 аг8. a{t а2<Р a2t ВЛ 2Л В Л В 42
В = а{ + (Х}(Х2 S а2. Рис. 2.4. К определению стационарных полей температур тел с источниками энергии Рассмотрим частные случаи ах=а2 = а. Тогда W8. „ W8 ---, с, — t, ч-- 22 2 L 2а Подставим значения Ci и С2 в (2.50): .2 W8 ,1 t — t ч-(— . 2 а Л Л8 Последнее выражение преобразуем, введя новую систему координат: WL WL t — tc —-1- а Л WL WL( , т х а Л 2L 8 2 7 2L {L-x') WL W 2 2 -------—-----1---1 L- X )• а 22 8 г г S х = х---= x-L; L = —. 2 WL WL . , T.r2L-x'-L = + — (x'+Z)[- а Л 2L (2-51) Найдём максимальную температуру tmax , которая соответствует центру пластины x=L, а х'=0. 43
_ WL WL2 ^max , • а 2л Неограниченный сплошной цилиндр. Дифференциальное уравнение и граничные условия имеют вид d2t ! 1 dt W. dx2 х dx А ’ - =0, dt a "I 1 — £ ) — (J, t/x т=0 dx Л L Jx=i (2-52) Представим (2.52) в виде d z dr W — (x—) =---- dx dx 2 Решение этого уравнения dt _ W х — — “—л* ч- С|, dx 22 W 2 f X ч- О, in X Ч" О-,. 42 1 2 dt _ Wx C] dx 2Л x Найдём постоянные интегрирования С] и С2. При х = О, абсурдно, поэтому полагаем С/ = 0 и запишем решение в виде /=-ОО, что 42 Из условий на границе x=L определим С2. dt Wx WL 2Л ' dx 2Л x=L 2Л Л 4Л O = o. Окончательно температурное поле цилиндра имеет вид W t~tc 42 (L2 -х2) + WL 2а (2-53) Полный неограниченный цилиндр (рис. 2.4 б). Дифференциальное уравнение для неограниченного цилиндра имеет вид (2.52), его решение было получено выше. 44
w t =--x1 + C. In A* + C2 42 1 2 На границах x = Rjhx = R2 теплообмен происходит по закону Ньютона со средами, температуры которых равны г и г , а коэффициенты теплообмена - ах и а2 . Граничные условия имеют вид (2.54) Используя эти условия, найдём постоянные Cj и С> 1 W а а 1'Г С. + +ffa2) + ^^(/?2-/?2)]; 1 -q >-izx£?2 t?i' 2 х I z 4/L 4 *" J' "* Q = I[ Y (уг"++ (7г‘1п *2 - *21п я.) - у («Л1п *1 + М.1п *2) + л 4 л2 4Л 2. WX R„ R, a.tr a2tf. +----(—-----) + a}a2(tc 1пЛ2-? 1п7?.) + Я(-L +--)]; X jFJ JZXC| Z £?2 ' X JTJ jrj ’ уГ^> В = a,a, In—-+A а, а7 Пример. Полый электрический проводник охлаждается водой, текущей в его внутренней полости и естественной конвекцией на наружной поверхности, находящейся в воздухе. Размер проводника R} = 2 мм, /Г =3 мм, удельное сопротивление материала р=0,1 Ом мм2 / м, 2= 15 Вт/м-К, сила тока /=1000 A, «/ =1000 Вт/м2К, а2 =10 Вт/м2 К, -30°С, =20°С. Требуется определить максимальную температуру проводника. Решение: Найдём координату х0, соответствующую максимальной температуре dt dx г-х0 -п х0 22 " ’ J2 2С -^1 = 0,00298м W Поскольку Ri<xo<R2 , то максимальная температура достигается внутри проводника и можно воспользоваться полученными выше формулами (в противном случае максимальная температура наблюдалась бы на одной из поверхностей). С2 =1304,9; С. =120,5; Гтах =543,9°С, 45
при такой высокой температуре следует рассмотреть вопрос кипения воды. Шар с источником энергии. Дифференциальное уравнение имеет вид: d2t ^2 dt _ W dx2 x dx A’ 0<x<L (2-55) На границе x=L теплообмен происходит по закону Ньютона, поэтому граничное условие имеет вид dt Л---- x-L v с lx=L Максимальная температура должна быть в центре шара, что позволяет записать dt dx = 0. х=0 Покажем решение приведённой выше системы уравнений. Представим (2.55) в форме и2и_ W U(xt)\ W г— — — а -- — — Хол ♦ х dx-------------------2 \ dx ) 2 Интегрируя дважды последнее уравнение, получим d(xt) W х2 ч W х2 , J -----_ _-----1- С], d(xt)-~----------dx + С\ dx, dx Л 2 Л 2 t W х2 ~ Wx2 С2 tx —-----f С, х + С-,, t = —-1- (7i ч-, Л 2-3 1 2 62 1 х где С/и С2- постоянные интегрирования. Устремим х—>0, тогда t —>00, что физически абсурдно, поэтому Wx2 полагаем С2 = 0 и t =---+ С,. 62 1 На основании условия (2.56) определим Сц dt WL WL г WL2 „ — ------------л = а\------1-С-Ц1; dx ЗЛ ЗЛ 6Л 1 с запишем выражение для искомой температуры: 46
w w t-t = !L(l2-x2)+---L 6Л 3a (2-57) Приведём сводку формул для трех тел простейшей формы: 2л а W /Г2 2Ч W т ~ + — / / V 1 I I • I------\ Т_- A, I 1 Т_-, L 4Л 2а W 1Т1 2ч W г L “Л =— (ъ ~х )н----------L. ш с 6Л За t цил Во всех случаях температурное поле тела представляется параболическим законом. При ос=оо и х=0: _ WL2 Гпл ^~^2Л’ WL2 4Л WL2 6Л ' Найдём среднюю объёмную температуру пластины, цилиндра и шара; по определению (2-58) где V и dV- объём и элемент объёма тела. dVn -dx-1'l, КП=2Л.7Ч Кц=^£2-7, L<х<L• <7Рц = lnxdx'1, Q<x<L- <7ГШ = 4лх2<7х, о < х < L Подставляя в (2.58) значения разности температур из (2.51), (2.53) и (2.57) и производя интегрирование, получим после простых преобразований выражение для среднеобъемного перегрева п Р aV. Г = &v =--(7 + п--) v aS ЛЗ (2.59) где P=W-V - полная мощность источников в теле; V и S - объём и теплоотдающая поверхность в теле; п - коэффициент, равный для пластины, цилиндра и шара 47
1 1 3 пп ~ з ’ Пч ~ 2 ’ И,я ~ 5 (2.60) Формулу (2.59) можно также применять для оценки средней объёмной температуры тел иных конфигураций с распределёнными источниками энергии. Для этого следует все тела разбить на три группы: тела группы шара имеют все три измерения одного порядка; тела группы цилиндра имеют два конечных измерения одного порядка и третье измерение - неограниченно большое ( z «z «z3); тела группы пластины обладают одним измерением конечной величины с двумя другими неограниченно большими измерениями (z, <<Z2 « Z3 )• При оценке для тел указанных групп следует выбирать соответствующие значения п из (2.60) или промежуточные значения этих величин. 2.4. Нестационарный тепловой режим тела с равномерным полем температур 2.4.1. Дифференциальное уравнение процесса Рассмотрим тело произвольной конфигурации, которое имеет мощность источника W(r). В начальный момент времени т-0 тело имеет температуру t0 ~ ^0 5 (2-61) и вносится в среду, температура tc(r) которой изменяется во времени; теплообмен тела со средой подчиняется закону Ньютона (2.24). Задача состоит в определении температуры тела в любой момент времени. Температурное поле такого тела описывается дифференциальным уравнением (2.10) Я Л дт граничным условием (2.24) и начальным условием (2.61). Приведём эту систему уравнений к иному виду, применив ко всем членам системы уравнений следующую операцию: Рассмотрим сначала (2.10): 48
V2tdV 1 гW 1 rep dt + — —dV = — ——dV. V{ A V{ A dr (2.62) Используем теорему Грина J V~tdV = j gradtdS V s и преобразуем первый член (2.62), при этом будет использовано граничное условие (2.24): j gradtdS = -J ^dS = ~— (tt -tc), с с А А t = i to Тогда — [v2tdV = -~—(t -t). vl AV ’ ‘ По определению j WdV = P, [tdV — tv, ' V где P - полная мощность источников тепла в теле; t - средняя объёмная температура тела. В этих преобразованиях предположено, что коэффициенты а, А, с не зависят от температуры. Используем полученные результаты и перепишем (2.62) в виде dtv z ч -Р(г) — + moft = mQtc (г) + ——, dr С С = cpV , . (2.63) Если предположить, что поле температур в теле равномерное, то tv = ts = t, и уравнение (2.63) примет вид dt z ч Р(т) 4" THqI — dr----------------------С (2.64) 49
Система уравнений (2.61) и (2.64) полностью описывает температуру тела с равномерным полем температур. 2.4.2. Интегрирование системы уравнений (2.61), (2.64) Используем следующий метод определения t(r): представим эту величину в виде произведения двух неизвестных пока величин V и W: t(r) = V(r)-W(T) (2.65) и подставим (2.65) в (2.64): w(dV ! y^ydW. т{ t Р<г) dr 0 dr 0 c C Выберем функцию v = v(t) так, чтобы сумма в скобках обратилась в нуль: dV Г/ А ---h mt,V = 0; dr ° решение этого уравнения даёт выражение для V в виде: V = С^т°т (2.66) где Cj -постоянная интегрирования. Из (2.64) и (2.65) получаем дифференциальное уравнение для W: dW 1 r z . Р(г)1 ^=?[^(г)+—]. решение, которого имеет вид W = J е'""г IX (г) + ^-]dr + С2 т0С где С2 - новая постоянная интегрирования. Принимая во внимание это выражение, а также (2.64) и (2.66), получаем t(T) = e-^{D + m0\e’notFdi}, (2.67) 50
р где D = C{C2, F = rf.+ mQC Неопределённый интеграл в (2.67) заменим определённым, взяв его в пределах от г0 до т. При г = г0 последнее выражение примет вид to=e-m°r°[D + morfem°rFJr]. г0 Интеграл в последнем выражении равен нулю, отсюда Подставив значение D в (2.67), получим г(т) = + тое-щт Jеи"г Fdz. (2.68) T<i Далее применим к интегралу в правой части (2.68) формулу интегрирования по частям m^e^F^dr - jF(r)t/(ew’r) = [F(r)em“rXo - ^evF'(r)dr, го го после чего формула (2.68) примет вид t = + е {Fev - Faem°ra - je^F'dv}, (2.69) где FQ=tc0+-^-, mnC F = tc + P m0C Последнее выражение перепишем в несколько ином виде: Т t-F = (/0 -Fo)e‘m°(r’ro) -e~v ^F'dr. (2.70) Дальнейший анализ можно проводить, если задан вид функциональных зависимостей tc =tc(r) и р = р(т) . Здесь возможны различные сочетания, поэтому рассмотрим случаи, наиболее часто встречающиеся на практике. 51
2.4.3. Нагревание или охлаждение тела в среде с постоянной температурой Пусть Р=0 , и тело помещено в среду с постоянной температурой tc = const. Тогда Fo = F = tc = const, и зависимость (2.70) примет вид 0 = = *о ~ tc (2.71) Из последнего выражения следует, что разность температур тела и среды изменяется по закону экспоненты (рис. 2.5 а). Прологарифмируем формулу (2.71): In® - ~т0(т ~т0) = ~тйт + тото = ~тот + const. (2.72) Рис. 2.5. Охлаждение или нагревание тела в среде с постоянной температурой 52
На рис. 2.5, б дано графическое представление этой зависимости в полулогарифмических координатах. Из рис. 2.5, б и формулы (2.72) следует, что in 0. - In 02 1п(т - Г) - 1п(/ - /) Wo =----1----1 = ----L2—lZ (2.73) тг-тх тг-т, Последнее выражение позволяет определить параметр тй , который в дальнейшем будем называть темпом охлаждения (нагревания) тела, опытным путём. Пусть из опыта получена зависимость (t~tc) = /(г); построив её в полулогарифмических координатах Int/---/) — /Jr) и выбрав два каких-либо момента времени г, и г2, находим по формуле (2.73) величину т0. Полезно аналитическое представление темпа охлаждения тела т0 в иной форме: дифференцируя (2.72), получим d In 0 ----= -m„ или dr 0 d® d(t-te') (2-74) ---=------- Qdr (t-tc)dr Рассмотрим теперь простое нагревание тела в среде с температурой / > /. Вычтем из правой и левой частей уравнения (2.71) по единице. После преобразований получим = щи (2.75) г-/=(ге-^0)[7-^ММ£1>] Графическая зависимость представлена на рис. 2.5. Сделаем следующие тождественные преобразования с выражением (2.75): t-tQ=t-t0 +te -tc= (te -t); (/ — tQ) — (/ — t) = (/ — t0 )[1 — e ]; Представим, наконец, это выражение в виде , / -t In----------------------------= —тот + const tc ~ t0 53
и изобразим последнюю зависимость графически (рис. 2.5). Если имеется экспериментальная зависимость t = /(г) по нагреванию тела в среде с температурой tc >t0, то можно определить темп нагревания тела по аналогии с темпом охлаждения тела. В этом случае результаты опытов следует обрабатывать в форме ln(^ -fj-lnfc -t2) Г -г. (2-76) При простом нагревании или охлаждении тела с равномерным полем температур темпы нагревания и охлаждения, как это видно из всех предыдущих выводов, равны между собой. 2.4.4. Нагревание или охлаждение тела в среде, температура которой изменяется во времени с постоянной скоростью Рассмотрим следующий случай (рис. 2.6 е): te=b(r-T0) + tc0, р ~ Ро = 0. (2.77) Согласно (2.69) Fo = tc0, F =tc, и из формулы (2.70) следует, что =(f0 -be~m°r После преобразований получим ? — — ——г ’0 т —)^ 0 , (278) m0 c »l0 V • 7 Графическое представление разновидностей рассматриваемых режимов дано на рис. 2.6, из которого видно, что возможны случаи пересечения кривых и tc (г), однако, такое пересечение может быть только в одной точке, соответствующей максимуму (рис.2.6 г) или минимуму (рис. 2.6 б). Запишем условие пересечения t(T^ = tc(T,) (2.79) и найдём время т*, когда происходит пересечение кривых. Из (2.79) и (2.78) следует, что 54
Рис. 2.6. Охлаждение или нагревание тела в среде, температура которой меняется по линейному закону ь . ь ~~ Vo ~ ^сО X м0 т0 Откуда г. = г() +—1п(1 + ——— т0 Ытй (2.80) 55
Если t0 = te0, то пересечение происходит в точке т, =т0 (рис. 2.6 а). При t0 > tc0 и b>0 точка пересечения соответствует времени т, > г0(рис. 2.6 б). Рассмотрим отдельно случай С) •> 7) 1 Ь>0, Как следует из (2.80), аргумент In будет меньше единицы, и второй член в правой части этого уравнения всегда будет отрицательным, т.е. уравнение имеет решение лишь для г, < г0 . Последнее означает, что пересечение никогда не произойдёт, так как рассматриваются процессы, протекающие при г>т0, а г,>г0 (рис. 2.6 в). Для Ь<0 и tc0 > t0 выражение ——— >0 и, Ы тй следовательно, г«>г0, т.е. возможно пересечение кривых ?(г) и tc(r) (рис. 2.6 г). Наконец, для Ь<0 и tc0 <t0 выражение —— <0, т.е. г. < г0, что Ы т0 означает невозможность пересечения кривых ^(г)и гс (г) (рис. 2.6 д). Рассмотрим теперь, как изменяется ход кривой t(t) с течением времени. Со временем второе слагаемое в (2.78), содержащее экспоненту е~Мо(г~Го) в качестве множителя, становится пренебрежимо малым по сравнению с первым, т.е. разность температур тела и среды стремится стать постоянной b т0 (2.81) при больших значениях (г-т0). Из рис. 2.6 видно, что с течением времени кривые г(г) и гс(г) становятся практически параллельными. Обозначим время, начиная с которого с заданной степенью точности можно не учитывать второе слагаемое в (2.78), через rui тепловой режим тела при г > тР назовём регулярным режимом второго рода, а при г < тР - иррегулярным (дорегулярным) режимом. Слово регулярный {regular) дословно означает «правильный», «закономерный». Здесь оно употребляется в смысле «упорядоченный» и в рассматриваемом конкретном случае означает, что спустя некоторое время начальное температурное состояние перестаёт влиять на температурное поле тела. Иногда такой режим называют «квазистационарным». Из сопоставления общего выражения для температурного поля тела (2.78) с выражением (2.81), а также анализа рисунков 2.6 следует, что для всех случаев соотношений t0 и tc0 при т>тр разность температур (tp-tc) изменяется по одинаковому закону (2.81). В развёрнутом виде уравнение (2.81) примет вид 56
tP(r^tc^b(,T-T^-—, b<Q (2.82) Выражение (2.78) является точным решением задачи, а (2.82) — приближенным. Для некоторых технических задач часто бывает удобно использовать такого рода приближенные выражения. В заключение заметим, что процесс нагревания или охлаждения тела в среде, температура которой изменяется во времени с постоянной скоростью, в чистом виде встречается редко, а для своего осуществления требует применения специальных средств регулирования. Однако некоторые явления, протекающие в естественных или лабораторных условиях, можно в первом приближении свести к рассмотренному здесь процессу или к их совокупности. Например, охлаждение метеозонда при его подъёме можно с довольно грубым приближением описать с помощью приведённых выше уравнений. 2.4.5. Температурный режим тела, помещённого в среду с гармонически меняющейся температурой Простейший закон периодического (гармонического) изменения температуры среды: — — tc =tc + Лсоз(——) =tc + ЛсоэбУг, (2.83) где tc -среднее значение температуры среды, около которой происходят её колебания, Т- период колебаний, <у = —. Применим для отыскания разности температур (t-tc) общую формулу (2.70), в которой F = tc, F0=tc0, F'= —Act>sin а>т . Тогда t -tc = (f0 - tL.o)+ AaxT^ jev sin axdv. Интеграл в правой части является табличным: г е'и°г sin<yrc/r = —--------------™(zn0 sintwr-cycoscyr); J m0 +су 57
применив его, получим t-tc=(t0- tc0 )е " Го) + А* ое m°r [e";,:‘r sin от cos р - e"*1 cos от sin Р -(етйТй sin суг0 cos Р - е44 cos от^ sin /?)] или t-tc =(z0-r£0)e-M“<r-r°) + Z(ye-Mor[sin((yr-^)-sin((yr0-^)e’Mo<r’r“)] . Здесь использованы следующие обозначения и формулы: tgp = —, p-arctg — , А* = г = Л sin р, mo mo М + О2 , V ° (2.84) 1 _ 1 1 _ cos2 р ml + о2 ml 1 + tg2 Р ml Разность температур представим теперь в следующем виде: i-rc = [/0 -tc0 - A" sin(6yr0 -Р)]е т<,{г ‘°-’ + A* sin(oT-p) (2.85) По прошествии некоторого времени от начала процесса экспоненциальный сомножитель в (2.85) становится столь малым, что первым членом в (2.85) можно будет пренебречь по сравнению со вторым, т.е. = Л* sin(«?r—/?) * (2.86) где индекс «р» означает, что рассматривается только та часть процесса, в которой начальное температурное состояние (?0,/с) , а также момент фиксации начала процесса тоуже не играют роли, т.е. изучаемый процесс вступил в наиболее простую, упорядоченную (регулярную) стадию. Этот температурный режим называют в литературе регулярным режимом третьего рода, или квазистационарным режимом. Более наглядное представление о колебаниях t получим, если сравним колебания обоих температур tc и t около среднего значения te (рис. 2.7). Для этого в последней формуле заменим tc с использованием выражения (2.83), а Л* - через 4sin/7; тогда получим / - - Л cos от = A sin р sin(tyr - Р) или (2.87) 58
t-tc= В cos(cyr - 0), В - A cos 0 Рис. 2.7. Температурный режим тела, помещённого в среду с гармонически изменяющейся температурой Из (2.87) следует, что амплитуда В колебаний температуры тела в cos 0 раз меньше амплитуды А колебаний температуры среды; отставание по фазе дано величиной 0, которая согласно (2.84) заносит от периода колебаний Т и параметра т0. 2.4.6. Термическая инерция тела Рассмотренные выше три случая изменения температуры тела в среде с переменной во времени температурой нашли широкое применение в задаче о тепловой инерции различных термоприемников. Из анализа полученных ранее формул видно, какую важную роль во всех типичных случаях играет параметр т0. Величину, обратную т0, ) с т0 aS (2.88) называют показателем тепловой инерции. В том случае, когда tc = const, параметр е всецело определяет 59
быстроту приближения системы к тепловому равновесию со средой. В двух других случаях тепловая инерционность выражается более сложным комплексом, который, однако, по истечении некоторого времени становится прямо пропорциональным s. В табл. 2.1. приведены значения комплекса U, определяющего тепловую инерцию тела для трёх вышеприведённых случаев. Таблица 2.1 Закон изменения Характеристика тепловой инерционности tc - const £ tc =4о+Хг-го) be tc =tc+A cos cot eco Изложенная выше теория справедлива не только для однородного тела, но и для системы тел, если выполняется основная предпосылка: температурное поле системы тел равномерно. Пусть полная теплоёмкость С системы: т С = £ с> i = 1 где С, - полная теплоёмкость г-ой части системы, п - число тел в системе. где S - наружная площадь системы, участвующая в теплообмене со средой. Укажем ещё один важный параметр, служащий для характеристики инерционных свойств тела - это время установления Z температуры тела или системы тел. Временем установления Z системы называют время, по истечении которого разность ( t - tc ) температур системы и среды станет меньше заданной величины А (например, 0,05 град). Часто под А понимают ту небольшую разность температур, которая находится на пределе точности измерений, осуществляемых данной аппаратурой. Если tc неизменна во времени, то Z определяется из формулы (2.71), в которой положим (г-т0) = Z, |/-^| = А, а Оо “О будем считать заданными, тогда (2.89) 60
2.4.7. Внутренние источники энергии в теле Рассмотрим неизменные значения мощности источников и температуры среды: tc = const, Р = const Из (2.68) найдём температуру t тела в любой момент времени: р t = +---) [ e^d? m«c ИЛИ p p t = t 4"-+ —t—----° ( °) C -r v и C ' m0C m0C При т - co наступит стационарный режим, и температура тела станет равной ?™=^со) = ^+^7- тпС Учитывая это, перепишем формулу для t: t -t = (t -t lcm 1 Vcfn l0/^ (2.90) Эта формула аналогична формуле (2.71) для простого нагревания тела в среде с более высокой температурой. Таким же способом, как это было выполнено выше, возможно рассмотреть другие случаи определения температуры тела при изменении температуры среды tc или внутренних источников Р по иным законам. 2.5. Стационарное температурное поле стержней и пластин 2.5.1. Особенности теплообмена стержней и пластин. Дифференциальные уравнения Характерной чертой стержней и пластин является малый градиент температуры в поперечном сечении этих тел, обычно его считают равным нулю. Примером стержней могут служить проводники радиотехнических деталей, термоэлектроды термопар, проводники термометров сопротивлений и т.д. Стержень вытянут в одном измерении, в котором движется поток тепла, остальные два размера тела по сравнению с первым 61
малы. В пластинах, толщина которых по сравнению с другими размерами тела мала, тепловой поток движется в плоскости. Если температурное поле пластин обладает круговой симметрией, то их называют дисками. К ним можно отнести рабочий элемент полупроводникового выпрямителя радиаторного типа, отдельное ребро радиатора; шасси, на котором производится монтаж детали и т.д. Однако следует помнить, что «стержень» и «пластина» - это физические модели и связь с геометрическими размерами не всегда однозначная. Процесс распространения тепла в стержнях и пластинах существенно отличается от процесса распространения тепла в стенках: тепловой поток, протекающий через любую изотермическую поверхность стенки без внутренних источников тепла, будет в установившемся режиме неизменным, а через различные изотермические поверхности стержней и пластин проходит разный по величине тепловой поток. Это объясняется тем, что при передаче тепловой энергии кондукцией в стержне или пластине происходит непрерывное рассеяние или приток тепловой энергии с поверхности этих тел в окружающую газообразную или жидкую среду благодаря конвекции и излучению. 2.5.2. Дифференциальное уравнение теплопроводности для стержня Температурное поле стержня или пластины может быть получено на основании решения соответствующего дифференциального уравнения для стержня или пластины вместе с граничными условиями. Такое уравнение обычно выводят специально, анализируя процесс переноса тепла в стержне или пластине. Однако оно может быть получено из общего уравнения Фурье, так как стержень или пластина представляют собою частные случаи тел с двух- и трёхмерным полем температур. Рассмотрим наиболее общий случай - анизотропное тело, все три размера которого имеют один порядок; пусть в теле равномерно распределены источники тепла удельной мощностью W Вт/м3 . Стационарное поле температур в таком теле описывается дифференциальным уравнением (2.16) для параллелепипеда Для ограниченного цилиндра в стационарном режиме можно использовать уравнение (2.8), положив в нём — = 0: дт 15,, dt . д dt ---у—) -I —у щ = о х дх дх dz-dz (2.91) 62
Рассмотрим сначала температурное поле параллелепипеда, на границах у = ±1 и z = ±/ которого справедливы условия (2.24) третьего рода (рис. 2.8): Обозначим <9 и, в силу симметрии температурного поля, перепишем граничные условия в виде а) (2.92) в) Условия на границах х = ±1х пока не оговариваем. Предположим, что градиент температур в данном теле в направлениях у и z мал и покажем, как можно учесть эту особенность температурного поля в дифференциальном уравнении (2.16), которое представим в форме + + + = О (2.93) дх ду dz Проведем над уравнением (2.93) следующую операцию: 1 h = (2.94) tyt? О О т.е. умножим почленно уравнение на dydz и проинтегрируем его в пределах от 0 до 1У и lz, а затем разделим результат на произведение lylz. Физически по отношению к ,9(x,y,z) эта операция означает осреднение температурного поля в теле по направлениям у и z. Обозначим осредненное в направлениях у и z значение температуры (рис. 2.8 а) I l>L- &(х) = = — J J&(x,y,z)dydz. (2.95) у z 0 О Итак, 63
д2& 1 У z О О ^-^-dydz = Л [— f f 3dydz] ~Л ^-у-; &2 ^х2 Z/JJ0 Л * дх2 д2& гду lylz OOV/ У z О О Рис. 2.8. К выводу дифференциальных уравнений параллелепипеда, пластины и стержня Поскольку в последнем выражении интегрирование и дифференцирование происходит по одной переменной у, то эти операции нельзя менять 64
местами. Рассмотрим интеграл О 53 дЗ_ & 1у дУ о Последнее преобразование сделано на основании условия (2.92 а). Следовательно, —f &(x,l ,z)dz = —-<9г(х), UVo 1у (2.96) где 1 — I &(x,ly,z)dz - значение осреднённой по z температуры на 4 О грани у=1у . Аналогично проведём операцию / над третьим членом уравнения (2.93). Используем условия (2.92): С/^ I (2.97) £/x) = -~p(x,y,/JJy, 1у О &у(х) - значение осреднённой по у температуры на грани z = /г. 1 1у1~ Iyz[W(x,y,z)\ - —jJW(x,у,z)dydz - JF*(x). lyh о 0 Если W(x,y,z) = const, TO W = W. Перепишем уравнение (2.93) с учётом преобразований . о2>.9(х) а а п _ 5х ly lz Введём обозначения «ы=т-; №Г=Ч'' (298) 65
Поясним физический смысл безразмерных параметров Ч' и т (см. рис. 2.8 а). Ч\ -отношение осреднённой по z температуры на грани у = I к средней по сечению yz температуре. В общем случае и Т должны изменяться с х. Используя (2.98), перепишем последнее уравнение в следующем виде: дх2 + Лх ’ b = —-------— ил а) (2.99) б) Если градиент температур по сечению yz стремится к нулю, то это означает, что средняя по сечению температура и температура на гранях iv и lz стремятся друг к другу, т.е. Ч* ->1 и ->1. Тогда выражение для Ь2 примет вид = (2.100) Уравнение (2.99), в котором Ь определяется по формуле (2.100), является дифференциальным уравнением для стержня с источником тепла. Если стержень имеет периметр U и площадь поперечного сечения 5, а теплообмен на границах происходит по закону Ньютона с коэффициентом теплообмена а , то можно обосновать следующее выражение для b: . 2 aU о --- (2.101) 2.5.3. Дифференциальное уравнение теплопроводности для пластины Рассмотрим температурное поле параллелепипеда, в котором градиент температур в направлении z мал, и составим для такого тела дифференциальное уравнение теплопроводности. Проведём над всеми членами уравнения (2.93) операцию /г: z о (2.102) 66
где f(x, у, z)- функция координат x,y,z. Обозначим (2.103) здесь <9](х,у) - осреднённая по оси z температура в точке х,у параллелепипеда (рис. 2.8 б). Итак, дх дх оу Оу ^дэ I [1 -- ] = 3 -- (/2 = ---) =^L-- Л г&2 Z?o г&2 zj ^dzJ lz dz (2.104) = -у-«9(х,у,/2). о h Последнее преобразование сделано на основании условия (2.92). Iz[W{x,y,z)] = y\w{x,y,z)dz = Wx{x,y). (2.105) 4 0 В частности, при W = const имеем W, = W. Запишем теперь уравнение (2.93) после проведения операции Д над всеми его членами: д2&{ д2Ц а Л ----7- + Л.--7- —- &(х, у, 1„) + W, = 0. х дх2 у ду2 i2 1 Обозначим (2.106) где Т, — отношение температуры в точке х,у на грани z = lz параллелепипеда к средней по направлению z температуре на отрезке 0< z < Zz с теми же координатами (х,у) (рис. 2.8 б). Учитывая (2.106), перепишем последнее уравнение , д2% . д2&{ ,2п nr п г? + л’~ву" ^=^11 к (2.107 ) 67
Если в направлении оси z градиент температур стремится к нулю, то -> 1 и -> . Этот случай соответствует пластине с двумерным полем 4 температур. 2.5.4. Дифференциальное уравнение для диска Приведём теперь операцию Iz над всеми членами дифференциального уравнения (2.91) для ограниченного цилиндра с источниками энергии. Обозначим при этом (рис. 2.8 в) I92(x) = p|I9(x,z)Jz, (2.108) о Тогда Г1 д 6<9Л1 1 d , д&2 ЛЕ-—(2^ —)] = я )• х дх дх х dx дх Операции lz над вторым и третьим членами уравнения (2.91) приведут к выражениям (2.104), а 1 /Z[f7(x,z)] = — J W(x, z)dz = W2(x). h о После проведения указанной операции уравнение (2.91) примет вид х dx dx л. Ы _ ^2 ф _ ^Л) 2 2Л ’ 2 .92(х) а) (2.109) б) Здесь Т2 — отношение температуры в точке х на грани z = lz ограниченного цилиндра к средней по длине цилиндра температуре в той же точке х (рис. 2.8 в). При стремлении градиента температур по длине цилиндра к 0 &(x,lz) -> <92(х) и . Для Т2=1 приходим к дифференциальному уравнению (2.109) для температурного поля пластины с круговой симметрией; такие пластины называются дисками. Температурное поле пластины с двумерным полем температур описывается уравнением (2.107), а для диска - уравнением (2.109). Прежде 68
чем переходить к решению этих уравнений, проведём некоторое их обобщение, позволяющее рассмотреть следующий более широкий класс задач: температура среды с разных граней параллелепипеда или со стороны торцов цилиндра разная, тогда dt ос' dt ос" F-+=0, =°- oz л - oz Л Если при проведении преобразований, приведших к зависимости (2.104), учесть условия на границах в приведённой здесь форме, то уравнение (2.107) сохранило бы свой вид 21 Z но значение & было бы другим, а именно: a' t + ОС" t" 9(x,y) = t(x,y)-t (2.110) a z+a \ v 7 где гсэф- эффективная температура среды. Если ос \ = а\, то ^=0’5(С+С)> (2.111) т.е. эффективная температура в этом случае будет равна среднеарифметической температуре двух сред. Тот же вывод остается справедливым и для уравнения (2.109). 2.6. Критерии неравномерности поля температур в теле Критерии Ч'2,Ч'у,'Г1,Ч'2, определенные выражениям (2.98), (2.106), (2.109), характеризуют неравномерность температурного поля в поперечном сечении тела. Можно попытаться найти точное значение для этих критериев. Для этого необходимо иметь точные аналитические выражения температурного поля и произвести с ними операции, указанные формулами (2.98), (2.106), (2.109). В общем случае эти критерии будут изменяться с координатами, что затруднит интегрирование дифференциальных уравнений (2.99), (2.107), (2.109). Для случая Ч'2=Ч' =Ч'1=Ч'2= 1 (стержни и пластины) интегрирование этих уравнений, как будет показано ниже, не представляет труда и приводит к сравнительно простым выражениям. Однако для того, чтобы знать 69
величину критерия Ч* в каждом конкретном случае и быть уверенным в возможности свести задачу к пластине или к стержню, необходимо получить хотя бы приближенные зависимости для оценки 4х . Предположим, что критерий ч* не зависит от координат; численные значения этого критерия для многомерного и одномерного температурного поля незначительно отличаются друг от друга. Эти два предположения существенно упрощают задачу отыскания аналитического выражения для 4х. Изложим один из путей решения этой задачи. Рассмотрим процессы переноса тепла в твёрдом теле и покажем физический смысл приведённых выше операции осреднения. На рис. 2.9, а показан элементарный объём dxdy, в левую грань которого входит удельный поток 0[(х,у),а из правой - выходит 02(х,у), с верхней и нижней граней рассеиваются удельные потоки 0з(х, у) и 04(х,у) . Обычно при составлении уравнения баланса энергии рассматриваем малый объём dxdy и рассчитываем ту часть потока QT, которая за счёт теплопроводности пришла (или ушла) в данный объём Qt = (<h -<k)dy + (03 - 04)rfx, а затем, учитывая другие источники или стоки энергии Q, на основе закона сохранения энергии QT = Q. Если теперь произвести операцию осреднения по одной из осей, например у, 1 /[у] = — j fdy, 1у О то от локальных потоков 9](х,у)И 02(х,у) перейдём к осреднённым потокам 1 2i(*) =—J Ч[(х,у)ф^ 02 (*) = у- f <72 (х, y)dy, > О ‘у а а вместо локальных потоков 03(х,у) и 04(х,у) будем рассматривать их частные значения при у = 0 и у = у, т.е. 04(аО). Тогда пришедший (или ушедший) в элемент I dx поток выразится в виде Qt = [21 (*) - 2 2 + [2з (А I у) - 2 4 (x,0)]Jx. 70
Другие источники или стоки для объёма lydx дадут поток Q = |Qdy, и закон О сохранения примет вид QT = Q. Рис. 2.9. Анализ физического смысла операций осреднения б) Из этих рассуждений следует, что при осреднении по одной из осей (например, у), разность между вошедшим и вышедшим из объёма потоком, во-первых, рассматривается для осредненных значений потоков (рис. 2.9 б), а во-вторых, эта разность потоков может быть представлена в виде внутреннего (в данном элементе 11х ) источника энергии постоянной мощности. Следовательно, для оценки критерия неравномерности температурного поля в каком-либо направлении возможно рассмотреть простейшую задачу об одномерном температурном поле тела с постоянным источником энергии. Заметим, что здесь идёт речь не о точном вычислении критериев неравномерности поля, а только об их оценке. Например, для неограниченной пластины толщиной 21х с постоянным источником энергии и условиями третьего рода на границах х = ~1Х и х = 1х на основании (2.15) имеем . w п1 2. Wlx t — t — iz-(I ~ X ) 4-. c 2Л a Найдём отношение средних на поверхности ( л = 1х) и по объёму перегревов t9s и &v. Из приведённого выражения следует, что 71
Wlx a a &v =— &dx=—-(1ч-----). IA a 3 Я Отношение температуры <7 и &r по определению является критерием неравномерности поля температур. Для пластины равно 1Ч--В/ 3 (2.112) Температурное поле неограниченного цилиндра с источником тепла описывается выражением (2.53). Определяя температуры и для цилиндра: = 6'lx=R ’ &V = f &(X)xdx Я О и беря их отношение, найдём критерий ч'ц: =___1= ЕЗ. " 1Да-’ ^ИЗ) 4 Предположим теперь, что критерии ¥ не зависят от координат, и численные значения этих критериев для одномерного и многомерного полей температур незначительно отличаются друг от друга. Специально проведённый для многочисленных конкретных случаев анализ показал, что для стационарного режима сделанные допущения, а также способ оценки критерия Т по формулам (2.112) и (2.113) для простейших тел с равномерно-распределёнными источниками тепла приводят к удовлетворительным результатам. Зависимость (2.112) позволяет оценить критерий Ч' и Ч'_ для приведённых выше случаев. у 1а/. 1 + — 3 2 П/ _ П/ _ \Т/ __1 z 1 2 , 1 «Л ’ 3 Я (2.114) 72
Эти формулы позволяют решить вопрос о возможности использовать уравнения (2.99), (2.107), (2.109) для нахождения температурного поля стержня и пластины. Пример. Параллелепипед, размеры граней которого равны 21х = 20 см, 21 = 2/г = 8 см , обменивается теплом с окружающей средой по закону Ньютона, коэффициенты теплообмена на границах параллелепипеда = 15 Вт/(м2-К), ay=az=5 Вт/(м2-К). Коэффициенты теплопроводности параллелепипеда в направлениях осей х, у, z: Л, = 0,2 Вт1(м-К), Лх =Л2 =0,4 Вт!(м-К). Оценить возможность рассматривать температурное поле такого тела по формулам для стержня. Решение. Найдём значение критериев Bi и вычислим по формулам (2.114) критерии Тх, ¥у, Tz: «А=1ЛЛ02 = 75; ^5^ 1 Л 2-10 у 4-Ю’1 ^=^ = 0,286; 1 + — 3 =-4у = 0,86. 1 + — 3 Следовательно, по осям у и z отношение среднего поверхностного перегрева к среднеобъёмному близко к единице, т.е. градиент температур в этих направлениях мал. В направлении оси х существенна неравномерность поля температур. В данном случае параллелепипед можно рассматривать как стержень. Если в дифференциальных уравнениях (2.99), (2.107), (2.109) полагать Д' не равным единице, а подставлять их значения, вычисленные по формулам (2.113), (2.114), получим более точные выражения для температурного поля стержней и пластин. Естественно, такие поправки следует вводить для значений Д' , близких к единице, например, Т > 0,8 . Для Т < 0,8 , по-видимому, следует искать точное решение задачи. 2.7. Температурное поле стержня с источником тепла Пусть стержень имеет постоянное поперечное сечение S и периметр U. Теплообмен стержня со средой описывается законом Ньютона, средний коэффициент теплообмена равен а. Материал стержня имеет постоянный коэффициент теплопроводности Л. На одном торце стержня х = 0 задан тепловой поток Р , а на другом торце х = 1х теплообмен происходит по закону Ньютона с коэффициентом теплообмена аТ . Температура окружающей среды всюду постоянна и равна te. 73
Полагая в (2.99) W* = const = w и принимая Ь2 по формуле (2.103), получим дифференциальное уравнение, описывающее температурное поле стержня (рис. 2.10): W d23 dx2 -b2& + = 0, 2 aU b =7s’ & = t(x) - tc. A Условия на границах имеют вид Рис. 2.10. Стержень произвольного сечения Общий интеграл равен &{x) = C{chbx + C2shbx + -—^. (2.115) Из граничных условий находим aTW — ASb ASb A^ibshblt + ^chbl^ Л 74
chblx+—shblx \ + — thblx л_ x Ab Ab shbL + — chblx thbL + — Ab Ab Окончательное решение примет вид P W & =----(Achbx - shbx) + —- (1 - /LSb f—chhx Ab shbl + — chbl. Ab (2.116) Рассмотрим частные случаи. а) Теплообмен с торца х = 1х отсутствует означает что ат = 0 и (2.116) примет вид (ад иабатич еский торец). Это или Р W & =----(Achbx - shbx) ч--, ASb Ab А = cthbl P chb(lx - л) И' “ ASb shbl + ~AtP’ (2.117) б) То же, что и (а), но отсутствуют внутренние источники энергии W - 0. $ _ P chb(lx - x) ~ AbS shbl (2.118) Обозначим температуру торца х = 0 в этом случае через <90 и найдём связь между <90и Р ASb P = 9flASbihblx. (2.119) 3 = ? cthbl Условие (2.119) позволяет записать формулу (2.118) для случая, когда задана на границе х = 0 не мощность Р , а перегрев ; действительно, chb(lx - х) (2.120) Из сопоставления (2.116) и (2.117) следует, что расчётная формула существенно упрощается, когда соблюдено условие аг = 0. Покажем, как 75
для ат^0 можно решение (2.116) привести к приближенной формуле, имеющей структуру (2.117), из условия сохранения общего теплового баланса. Для того, чтобы в (2.117) учесть теплоотдачу с торца стержня (снять условие = 0 ), следует условно увеличить площадь боковой поверхности стержня путем увеличения длины стержня. Тогда условная длина стержня станет равной . Эти параметры связаны очевидными соотношениями: (Jx -lx)-U а = SaT Откуда (2.121) Подставив в уравнение (2.117) вместо 1Х условную длину /’ , учтём теплообмен торца с окружающей средой. в) Полуограниченный стержень 1х=<х> . Выражение для А примет вид А\ = 1, а 1^=00 ----в Ч KSb Xb2 (2.122) Здесь использованы известные зависимости chx - — (е* +е х); shx = — (е* ~ е х). Пример 1. Для увеличения рассеяния тепловой энергии с поверхности 50 торца функционального микроминиатюрного элемента выведено п медных (Z = 400 Вт/м-К) жил диаметром 27? = 10' м, длиной / = 10 мм (рис. 2.11). Пусть 50=1 см2, коэффициент теплообмена между поверхностью и средой при отсутствии жил ао=1О Вт/м2К, коэффициент теплообмена одной жилы с окружающей средой а = 20 Вт/м2К. Найти, во сколько раз увеличится отток тепла с поверхности при наличии жил. При решении задач считать, что наличие жил не изменяет значений коэффициентов теплообмена а и а,,. Решение. Обозначим мощность, рассеиваемую с поверхности жил, -Р; мощность, рассеиваемую оставшейся поверхностью (при наличии жил), -Д. Р + Р Требуется вычислить отношение---L. 76
Рис. 2.11. Увеличение интенсивности теплообмена с помощью медных жил Найдём Здесь через О0 обозначен перегрев над средой поверхности So . Если пренебречь количеством тепла, рассеиваемым с открытого торца жилы в среду, то на основании (2.119) P~n&aZSbihbl Значение Pt найдем из очевидного равенства = а,Р0(50-ял/12) Итак, tiR2 nkbSthbl - an.S’(1 (I -/z-) So nXbS , nR2 —=----------thbl + l-n-----; ao*^o = 200-b 6 = 14,1-; M M Po Ct о 2 _ ct(7 _ 2a _ 2 20 " IS ~ kR ~ 400 0,5 10’3 bl = 0,14; thbl « 0,14; P + P — 1 sl^z-0,62. Po Если n = 10, то рассеяние увеличивается в 7 раз. Замечание. Вывод (2.119) сделан в предположении, что основание жилы 77
имеет температуру 90, а коэффициент теплоотдачи на поверхности жилы а остается неизменным. Однако с увеличением количества жил, если не принимать специальных мер, интенсивность теплообмена на их боковых поверхностях будет снижаться и последняя формула дает завышенный результат. Дальнейшее увеличение количества жил приведет к противоположному эффекту. Пример 2. Полупроводниковое термосопротивление типа бусинки имеет следующие параметры (рис. 2.12): диаметр бусинки d{=\ мм, диаметр медных проводников d2 = 0,1 мм, их коэффициент теплопроводности Х2 = 400 Вт/м-К, длина проводников во многом превышает их диаметр, т.е. можно считать 12 = со . Бусинка находится в воздухе, коэффициент теплообмена проводов и бусинки со средой а = 40 Вт/м2К. Рис. 2.12. Полупроводниковое сопротивление бусинкового типа Необходимо найти величину измерительного тока, при котором погрешность измерения температуры tc воздуха из-за перегрева термосопротивления не превысит 0,1% ; tc = 40 °C. Электрическое сопротивление R =2500 Ом при t = 40 °C. Решение. Обозначим допустимый измерительный ток через I# , он связан с рассеиваемой допустимой мощностью Р$ зависимостью = откуда где R - электрическое сопротивление термистора. Погрешность измерения 5 за счёт перегрева термосопротивления равна 8 = LiLl = 2. = 0,001 = Ю”3,^ -И0 “Згс (а) Это значение 9 будем считать допустимым и обозначать &д. Установим 78
связь между i9( и Pd. Для этого представим термистор в виде бусинки с равномерным полем температур и двумя стержнями (рис. 2.12 б). В стержни отводится поток Р2 , а со свободной поверхности ПТС рассеивается поток Pt. Полная мощность Р() равна р = р +Р 1 д 1 1 1 2 (б) Найдём : Р{ = a&.S, = а&д (nrf — 2) = &д 1,25 Ю’4 Вт Полагая, что Р2/2-поток, входящий в один торец бесконечного стержня с температурой основания по формуле (2.119) найдём ^2 п 2 L С .7 LI Z. I i 4tz Ttd2 — &dA2b2S2thbl1, b2 - I - J , S2 = — 2 у A2R2 у Л2а2 4 Поток, входящий в оба проводника, равен J4-40 77-10 s ----= 4 • 10’4 Вт. 400 10 4 4 Итак, Рд = (1,25 + 4) 10 4 = 5,25 10 4 &д Вт. Из (а) и (б) находим Р, =5,25-10 4-10 ''-40 = 2,1-10 s Вт. Из (а) также следует, что перегрев термосопротивления при этой мощности будет равен ,9 = Ю-3^ = 0,04К, т.е. электрическое сопротивление следует брать при температуре t = 40,04 °C; будем считать его мало отличающимся от заданного R =2500 Ом при t = 40 °C. Тогда 79
f2,110 5 2,5 IO3 -92 мкА 2.8. Температурное поле пластины и диска 2.8.1. Обобщенное решение уравнения Бесселя Постановка задач теплопроводности в цилиндрической системе координат часто приводит к уравнению Бесселя. В частности, уравнение теплопроводности (2.109) для диска является уравнением Бесселя. Уравнение этого типа путём преобразований обычно приводится к канонической форме, для которой известны решения. Иногда бывает довольно сложно найти такую замену независимой переменной, которая позволила бы преобразовать (если это вообще возможно) заданное дифференциальное уравнение в обычное уравнение Бесселя или же такое уравнение, общее решение которого содержало бы линейные комбинации функции Бесселя. Однако простое сравнение заданного дифференциального уравнения с приводимым ниже обобщённым уравнением Бесселя позволяет просто ответить на вопросы: сводится ли данное дифференциальное уравнение к уравнению типа Бесселя и каково решение этого уравнения. Одна из форм обобщенного уравнения Бесселя имеет вид [7]: х2 + [(1 - 2А) х- 2Bx2]^~- + [C2D2x2c +В2х2 - В(1 - 2А)х + А2 - С2п2]у - 0, (2.123) dx dx где и определяет порядок уравнения Бесселя. Обобщённое решение имеет вид [7]: у = хАев*\С{ JK(Dxc) + С2У„ (Dx с )], (2.124) где Ci и Сг -постоянные интегрирования; Jn и У - функции Бесселя первого и второго рода порядка и. Рассмотрим простейшую форму уравнения Бесселя вида dx2 dx Его решение имеет вид у = СЛ(х) + С2У0(х), (2.125) 80
где J0(x) и У0(х) - функции Бесселя вещественного аргумента нулевого порядка, вторую из них иногда называют функцией Неймана. Напомним некоторые свойства функций Бесселя. Первые производные функций нулевого порядка JQ(x) и У0(х) являются функциями Бесселя первого порядка (2.126) У'0(х) = ^(х), У’0(х) = ^(х) Рис. 2.13. Функции Бесселя и Неймана вещественного аргумента нулевого и первого порядков На рис. 2.14 представлен вид функций Бесселя /0(х),^0(х) нулевого порядка и первого порядка Х’Дх) и Д(х) . Функция A?v(z) называется функцией Макдональда. Функции А?0(х) и К(х) при малых значениях аргумента могут быть представлены в виде ^0(х)»1п(-)™ 0,577 X К} (х) » х 1 х < 0,5 (2.127) 81
5 4 3 2 1 О 2 3 X Рис. 2.14. Функция Бесселя мнимого аргумента К0(х), Ki(x), Iq(x), Ii(x) 2.8.2. Круглое ребро постоянной толщины На рис. 2.15 изображено круглое ребро, радиусы которого Т?7и R2, а толщина 5 ; поверхности ребра находятся в теплообмене со средой по закону Ньютона, причём со стороны верхней поверхности коэффициент теплообмена а’г, температура среды t'cz , а со стороны нижней поверхности эти параметры равны a"., t”cz . На границе х = R{ в ребро входит постоянный поток Р, на границе х = R2 теплообмен со средой настолько мал, что им можно пренебречь. Теплопроводность материала диска 1. Требуется найти температурное поле диска. Рис. 2.15. Круглое ребро постоянной толщины 82
Задача сводится к решению уравнения (2.109), которое в нашем случае примет вид гс12& d& ,2 2q х —г + х--------Ь х+ dx dx &(x) = t(x)-tc3ip, a'+a\ ^-0. л _ a’z+a" лЗ (2.128) Граничные условия Р = -Л— -2nRx-8, х=1> d& dx x-R = 0 Сопоставив соответствующее (2.128) однородное уравнение с обобщённым уравнением Бесселя (2.123), находим А=В=0, п=0) С=1 d2 =_ Ь2 Используя обобщённое решение (2.124), получаем для соответствующего (2.125) однородного уравнения решение &0 = С, JQ(ibx) + C2Y0(ibx) = C3I0(bx) + C4K0(bx). Частное решение (2.128) имеет вид ,9; =—- тогда общее решение уравнения (2.128) & = С3Ц (Ьх) + С4К0 (Ьх) + —у; Ли Сз и Сд найдём из граничных условий. Окончательно Р Кх (bR2 )/0 (bx) + Ко (fec)7, (bR2) W 2xR^Ab' Il(bR2)Kl(bR1)-Il(bRl)Kl(bR2) + Ab2' (2.129) Приведём последнее выражение к форме, в которой фигурируют только безразмерные параметры: 83
I. (у)^0 (/ —) + K{ (y)/0 (/ —) 1 1 ° R2 ‘ ° /?, w^r2 + “psT 2.8.3. Эффективность круглого ребра постоянной толщины Рассмотрим дальше частный случай w = 0, т.е. круглый диск без источников тепла, и найдём эффективность ребра е: отношение полного количества тепла, рассеянного ребром, к тому количеству, которое рассеивалось бы в случае, если бы вся поверхность ребра находилась при температуре to основания ребра (или X = ooj, т.е. р0 = («X, >(й22 - )(г0 - гсэф). Эффективность ребра является важной технической характеристикой при выборе параметров ребра. Обозначим в = 11(/)^1(/1)-^1(х)11(Г1) тогда эффективность ребра /?,2) 5 (l + /?2//?,)£• la' + d" £ = (R2-Ri)lEb_^ V do (2.131) Зависимость e = e(s,—) приведена на рис. 2.16. Если /?2->сю , то Л1 (7?7?2) -> оо и решение (2.129) принимает вид Р К0(Ьх) 2nR{A6b K^bR') (2.132) Покажем в заключение, что решение (2.129) можно обобщить, учтя приближенно отвод тепла с торца х = R2. Пусть толщина ребра значительно меньше (R2 - /?,), тогда некоторым увеличением радиуса R2 до /?20 можно учесть дополнительный отвод тепла с торца в среду с температурой гсэф. 84
Величину R16 найдём из равенства 2ж/?20 = 2жК2 + 2xR.,d’, Откуда R2o=R,(1 + 8 2.9. Анализ ошибки измерения температуры Ошибки за счёт оттока тепла по проводам термопары могут быть полностью устранены только в случае равенства температур пластины и окружающей среды. Такой предельный случай встречается редко, поэтому рассмотрим реальный процесс измерения температуры термопарой (рис. 2.17 а). Пусть поверхности 1 и 2 тонкой пластины омываются средами с температурами гс1 и гс2 с коэффициентами теплообмена ах и а2 ; термопара прикреплена к поверхности 1 пластины и находится в среде 1. Рассмотрим случай, когда /С1 > ?с2 , тогда термопара становится источником, а разность между тепловым потоком , подводимым к спаю, 85
и потоком Р2, рассеянным на лицевой стороне круглого источника, равна р = р{ - р2 мощности источника. Выразим Р{ и Р2 через термические проводимости (рис. 2.17 б) Рис. 2.17. К оценке погрешностей измерения температуры с помощью термопары — (Тц. (/с1 Z[ ) , Р,— (У2 (^1 , где ст1Г- тепловая проводимость термопары между спаем и средой 1, а2- тепловая проводимость между лицевой поверхностью 2 источника и средою 2. Тогда Р — (o']+ СГ2^е2 ) — ((Гц- + (Т2X (2.133) Для вычисления температуры можно воспользоваться соотношением (2.132) t =t I Р _ а^с1 + a2tc2 , 2 _ ’ ~~хГ + а-. Подставим в последнее выражение значение р: ^2^:2 ) (^ir *-^2 )^1 ^0 (^-^1) 2kR,дЛЬ K^bRJ Заметим, что - температура, зарегистрированная термопарой; 1езф - характеризует или температуру пластины в точке х = <®, или температуру пластины при отсутствии термопары. Разность -1сзф) представляет собой абсолютную ошибку измерения температуры контактным способом 86
с помощью термопары. Обычно а1Г »<г2, поэтому можно положить а2 =0, и последняя формула примет вид П 1сэф _ 1 (cl ~ $сэф | + 2П |}R ^1(^1) СТ1Г 1 K0(bR{) (2.134) На рис. 2.18 приведена зависимость hT> ч ~ J I >УЛ1 /’ М "Чэф ^17 5 = °\т где по оси абсцисс отложено отношение проводимости диска к проводимости термопары. Для определения <т1Г можно воспользоваться результатами решения задачи о стержне (2.119). Обозначим теплопроводность термопары , площадь её поперечного сечения - ST, параметр ь для термопары - через Ьт, её длину - 1Т; учтём также, что термопара состоит из двух электродов 1 и 2, поэтому (2.119) запишем в виде 87
Ру — — (7'1 1 )('^'7'1^7'1 ту^ЬЬ-ру!yj T 2 у 2/l у 2 *7у 2 2 (г 2 )" Далее примем во внимание, что длина термоэлектродов значительна, так что thbTlT~\. Тогда О ।у b'j'IЭ” 2 aTl^T1 Ьт'=ТТ~ Лтгп Г2’ 2 _ aT1UT2 Окончательно получим СГ|У Если провода термопары изолированы, то изоляцию можно представить по формуле для плоской стенки. В последнюю формулу в этом случае следует вместо aTi подставить a'Ti, найденное из формулы 1 1 Ai .......™.......•+ -- Ои -р- 0=1,2), где Az и Л- толщина и теплопроводность изоляции электрода. Вычисления по формуле (2.134) можно упростить, если вместо K{(bRy) и K0(bRy) подставить асимптотические выражения, справедливые для малых х = bR{: f 2'i /С,(х) = х”1, (х < 0,05); К0(х) = In - -0,577, (х < 0,05). I х) ^2 ~д дд [ j ду дд д д л д дд дд УД дд ^у И д И И ПГ1 J^y ДУ Ду д^|у,'уг ду д д ДЛ jjy Для линейных задач теплопроводности справедлив принцип суперпозиции (сложения) температурных полей: температурное поле тела, которое формируется в результате нескольких тепловых воздействий, может быть представлено в виде алгебраической суммы температурных полей, вызванных каждым из воздействий в отдельности. Под «тепловыми воздействиями» понимается любая неоднородность в уравнении теплопроводности и краевых условиях. Так, для уравнений (2.8), (2.14), это величина W, для граничных условий (2.20), (2.21) - t(r), для (2.24 - 2.25) - величины tc и q, а для начального условия (2.19) - T(x,y,z). Покажем на примере задачи (2.47) - (2.49) применение принципа суперпозиции. Поле температур в этой задаче формируется в результате 88
следующих тепловых воздействий: w,t ,t . Два последних воздействия одинаковы по своей природе, и их логично объединить (хотя можно рассмотреть порознь). Таким образом, искомая температура будет складываться из двух составляющих: где /,(х)- поле температур под действием только внутреннего источника =/ =0), a t2(x) - поле температур только под действием температур окружающих сред tc. и (JF=O). Соответствующие математические постановки задач для /,(х) и t2(x) имеют вид: т2. dx ( dt. А ( dt, —'Ч ~Г~ +^1^1(л=0) = 0’ + ^2^2(л=(?) =0, k dx )х=о к dx )x=s = "^-о’ =a2(^2“^)J=lS, к dx Л=о к dx )x=s Принцип суперпозиции полезен при качественном анализе температурных полей и тестировании программ. 89
Библиографический список 1. Михеев М.А., Михеева И.М. Основы теплопередачи. М., «Энергия», 1973. 2. Исаченко В.П., Осипова В.А., Сукомел А.С., Теплопередача. М., «Энергия», 1975. 3. Болгарский А.В., Мухачева Г.А., Щукин В.К. Термодинамика и теплопередача. М., «Высшая школа», 1964. 4. Жуковский В.С. Основы теории теплопередачи. М., «Энергия», 1969. 5. Кондратьев Г.М. Тепловые измерения. М.-Л., Машгиз, 1957. 6. Дульнев Т.Н., Тарновский Н.Н. Тепловые режимы электронной аппаратуры. Л., «Энергия», 1971. 7. Шнейдер П. Инженерные проблемы теплопроводности. М., ИИЛ, 1960 8. Спэрроу Э.М., Сесс Р.Д. Теплообмен излучением. Л., «Энергия», 1971. 9. Эккерт Э.Р., Дрейк Р.М. Теория тепло- и массообмена. М.-Л., Госэнергоиздат, 1961. Ю.Карслоу Г., Егер Д. Теплопроводность твердых тел. М., «Наука», 1964. П.Варгафтик Н.Б. Справочные по теплофизическим свойствам газов и жидкостей. М., «Наука», 1972. 12. Кутателадзе С.С. Основы теории теплообмена. Новосибирск, «Наука», 1970 90
Дульнев Геннадий Николаевич Тихонов Сергей Васильевич Основы теории тепломассообмена Учебное пособие В авторской редакции Редакционно-издательский отдел Санкт-Петербургского государственного университета информационных технологий, механики и оптики Зав. РИО Н.Ф. Гусарова Лицензия ИД № 00408 от 05.11.99 Подписано к печати Заказ № Тираж Отпечатано на ризографе 93