Инженерия / Теплопередача
Температура контакта двух тел через тепловую эффузивность
При идеальном быстром контакте двух массивных тел температуры их соприкасающихся поверхностей переходят к общей контактной температуре. Она является средним начальных температур, взвешенным тепловыми эффузивностями тел, а не просто их средним арифметическим.
Формула
При равных эффузивностях граница имеет температуру 60 °C. С ростом отношения эффузивностей результат приближается к 100 °C.
Обозначения
- $T_c$
- общая температура поверхности на идеальной границе сразу после установления теплового контакта, K или °C при единой шкале температур
- $T_1,T_2$
- начальные равномерные температуры тел до соприкосновения, K или °C, одна шкала
- $e_1,e_2$
- тепловые эффузивности материалов, e=√(kρc_p), Вт·с^(1/2)/(м²·К)
- $k$
- теплопроводность материала, Вт/(м·К)
- $ρ$
- плотность материала, кг/м³
- $c_p$
- удельная теплоёмкость при постоянном давлении; для твёрдых тел во многих оценках используют табличную теплоёмкость, Дж/(кг·К)
Условия применения
Модель предполагает два однородных изотропных полупространства, каждое из которых до контакта имеет равномерную температуру. Контакт возникает практически мгновенно, передача идёт одномерной теплопроводностью, а контактное сопротивление, поверхностные покрытия, тепловыделение и теплообмен с окружающей средой отсутствуют. Теплопроводность, плотность и теплоёмкость считаются постоянными в диапазоне температур. Оценка описывает начальный интерфейс в полубесконечной модели: характерная глубина прогрева должна быть мала относительно толщин тел.
Ограничения
Это не формула конечного равновесия всей пары тел: она описывает границу двух полубесконечных сред в модели идеального контакта. При конечных размерах спустя время температурное поле зависит от толщин, теплоёмкостей и внешних условий. Воздушная прослойка, шероховатость и контактное сопротивление могут заметно изменить мгновенный скачок температуры; при фазовом переходе, температурно-зависимых свойствах или сильном конвективном/радиационном теплообмене требуется другая модель. Не следует путать эффузивность e=√(kρc_p) с температуропроводностью α=k/(ρc_p): первая задаёт вес поверхности в этой контактной задаче, вторая — скорость распространения температурного возмущения.
Подробное объяснение
При соприкосновении двух тел тепло начинает проникать в каждый материал от границы. Температура границы может заметно отличаться от среднего двух начальных температур: высокоэффузивный материал способен принять или отдать больше тепла в раннем переходном процессе и сильнее удерживает контактную температуру возле своего исходного значения. Поэтому формула использует взвешенное, а не обычное среднее.
Для полубесконечного тела решение уравнения теплопроводности после ступенчатого изменения температуры поверхности содержит профиль функции ошибок. Поток через поверхность пропорционален произведению теплопроводности на перепад температуры, делённый на корень из температуропроводности и времени. Приравнивание потоков, уходящих в два тела, даёт равенство e₁(T₁−T_c)=e₂(T_c−T₂), где e=√(kρc_p); решение этого баланса и есть контактная температура.
Эффузивность имеет единицы Вт·с^(1/2)/(м²·К). Она возрастает как с теплопроводностью, так и с объёмной теплоёмкостью ρc_p. Если e₁=e₂, контактная температура лежит ровно посередине между T₁ и T₂. Если e₁≫e₂, T_c приближается к T₁; при этом холодное или горячее тело с меньшей эффузивностью меняет температуру сильнее. Это объясняет ощущение прохлады у металла при комнатной температуре по сравнению с древесиной, но ощущение само по себе не является измерением температуры.
Формулу применяют как оценку начальной граничной температуры в контактном охлаждении, при оценке взаимодействия капли с твёрдой поверхностью и при интерпретации теплового отклика материалов. Она помогает проверить направление эффекта до численного моделирования: результат обязан лежать между исходными температурами и быть ближе к температуре материала с большей эффузивностью.
Важное отличие от калориметрического среднего: для конечных тел после полного теплообмена используется сохранение общей энергии с массами и теплоёмкостями, если система изолирована. Контактная оценка для полупространств вместо этого уравнивает мгновенные потоки через границу. Обе записи могут давать разные числа, потому что отвечают на разные вопросы и относятся к разным временным и геометрическим режимам.
В расчёте используйте согласованные свойства при температуре, близкой к контактной. Цельсий допустим для взвешивания температурных значений, если обе температуры и результат выражены в °C: добавление одной и той же константы к T₁ и T₂ не меняет среднего. Если требуется термодинамическая запись абсолютных температур, ведите расчёт в кельвинах. Округляйте только итог, иначе отношение эффузивностей может быть искажено.
Как пользоваться формулой
- Определите две стороны границы и их равномерные начальные температуры до контакта; запишите обе температуры в одной шкале, например в градусах Цельсия.
- Для каждого материала найдите k, плотность ρ и теплоёмкость c_p в подходящем диапазоне температур и вычислите e=√(kρc_p).
- Проверьте режим: контакт должен быть близок к идеальному, а глубина прогрева за интересующее время — мала по сравнению с толщиной каждого тела.
- Подставьте начальные температуры и две эффузивности в взвешенное среднее T_c=(e₁T₁+e₂T₂)/(e₁+e₂).
- Проверьте, что знаменатель положителен, ответ лежит между T₁ и T₂ и расположен ближе к температуре тела с большей эффузивностью.
- Для конечных тел или заметного контактного сопротивления замените оценку моделью с геометрией, граничными условиями и соответствующими контактными свойствами.
Историческая справка
Математическая основа задачи появилась вместе с теорией распространения тепла. В 1822 году Жозеф Фурье опубликовал «Théorie analytique de la chaleur», где систематизировал уравнение теплопроводности и методы анализа температурных полей. Для двух тел в контакте взвешенная контактная температура следует из решений нестационарной задачи для полубесконечных сред: профиль температуры после резкого изменения граничного условия выражается через функцию ошибок, а равенство тепловых потоков на общей поверхности задаёт её значение. Это не отдельный эмпирический закон, приписываемый одному автору, а следствие модели Фурье при конкретных предпосылках. Учебники по теплопроводности XX века развили этот подход для контактов, полубесконечных тел и тепловых датчиков. Термин «тепловая эффузивность» стал использоваться для комбинации √(kρc), связывающей теплопроводность и объёмную теплоёмкость; в современных работах эта величина служит коэффициентом в выражении для потока на поверхности. Позднее идея применялась при исследовании контакта капель с нагретыми стенками, сенсорных методов и термического ощущения поверхности. Исторически формула отражает переход от стационарного потока к нестационарному процессу: в самый ранний момент важны свойства материала у поверхности и рост прогретого слоя во времени.
Историческая линия формулы
Модель основана на уравнении теплопроводности Фурье. Контактная температура получается при равенстве тепловых потоков через идеальную границу двух полупространств. В современной литературе коэффициент, который определяет вклад каждого материала в это среднее, называют тепловой эффузивностью.
Пример
Две большие пластины быстро соприкасаются; требуется оценить начальную температуру границы. До контакта первая пластина имеет T₁=100 °C, k₁=40 Вт/(м·К), ρ₁=7800 кг/м³ и cₚ₁=500 Дж/(кг·К). Для второй заданы T₂=20 °C, k₂=0,15 Вт/(м·К), ρ₂=600 кг/м³ и cₚ₂=1400 Дж/(кг·К). Предположим идеальный контакт и полубесконечные тела. Дано: T₁=100 °C, T₂=20 °C. Сначала считаем эффузивности: e₁=√(40·7800·500)=√156 000 000≈12 490 Вт·с^(1/2)/(м²·К). e₂=√(0,15·600·1400)=√126 000≈355 Вт·с^(1/2)/(м²·К). Подстановка: T_c=(12 490·100+355·20)/(12 490+355)≈1 256 100/12 845≈97,8 °C. Ответ: начальная температура границы около 97,8 °C. Проверка: она лежит между 20 и 100 °C и близка к 100 °C, потому что эффузивность первой пластины примерно в 35 раз выше. Эффузивность имеет размерность Вт·с^(1/2)/(м²·К), поэтому в отношении взвешенного среднего она сокращается и температура сохраняет единицу °C.
Частая ошибка
1) Подменяют контактное среднее обычным средним температур. Так получилось бы только при равных эффузивностях. 2) Подставляют температуропроводность вместо эффузивности: α=k/(ρc_p) отвечает за диффузию температуры, но не является весом в балансе начальных потоков. 3) Используют формулу для конечных тел и времени полного выравнивания; здесь требуется другая энергетическая модель. 4) Игнорируют контактное сопротивление, покрытия и шероховатость, хотя они могут управлять реальным интерфейсом. 5) Берут свойства, сильно меняющиеся с температурой, как постоянные без проверки диапазона. 6) Считают, что найденная температура сохраняется бесконечно: в полубесконечной идеализации она постоянна на границе, но температуры внутри распространяются вглубь, а конечные размеры нарушают эту картину.
Практика
Задачи с решением
Отношение эффузивностей по контактной температуре
Условие. Первое тело имеет температуру 100 °C, второе 20 °C. При идеальном контакте граница сразу нагрелась до 80 °C. Найдите e₁/e₂.
Решение. Из равенства потоков получаем e₁/e₂=(T_c−T₂)/(T₁−T_c). Подставляем температуры: (80−20)/(100−80)=60/20=3. Температурные разности заданы в одинаковых единицах, поэтому отношение безразмерно.
Ответ. Эффузивность первого тела в 3 раза больше.
Одинаковые эффузивности
Условие. Два полубесконечных тела имеют T₁=90 °C, T₂=10 °C и равные эффузивности e₁=e₂. Найдите контактную температуру.
Решение. При e₁=e₂=e формула сокращается: T_c=(e·90+e·10)/(2e)=(90+10)/2=50 °C. Ответ лежит между исходными температурами и равноудалён от каждой.
Ответ. T_c=50 °C.
Дополнительные источники
- https://openaccess.city.ac.uk/id/eprint/13580/3/2010_IJTS_Strotos_Non-dimensionalisation.pdf
- https://www.e-rara.ch/doi/10.3931/e-rara-19706
- https://www.enpress-publisher.com/files/journals/34/articles/2487/public/2487-9841-2-PB.pdf
Связанные формулы
Инженерия
Закон Фурье для теплопроводности
Закон Фурье для теплопроводности показывает, как получить расчетную величину из проверяемых исходных данных. Формула полезна для предварительного инженерного расчета, потому что сразу связывает результат с единицами измерения и областью применимости.
Инженерия
Конвективный тепловой поток
Конвективный тепловой поток показывает, как получить расчетную величину из проверяемых исходных данных. Формула полезна для предварительного инженерного расчета, потому что сразу связывает результат с единицами измерения и областью применимости.
Инженерия
Число Нуссельта для конвективного теплообмена
Число Нуссельта показывает, во сколько раз конвективный теплообмен у поверхности эффективнее простой теплопроводности через характерный слой среды.
Инженерия
Теплопроводность через плоскую стенку
Теплопроводность через плоскую стенку: формула Q=\lambda A\frac{\Delta T}{\delta} помогает величины Q, lambda, A, DeltaT заданы для одной и той же ситуации, периода или объекта. В тексте есть условия, пример, ошибки и проверка результата.
Физика
Температуропроводность вещества
Температуропроводность a характеризует скорость выравнивания температуры теплопроводностью. Она равна отношению теплопроводности λ к объёмной теплоёмкости ρc и измеряется в м²/с.