FEMM 5

Граничные условия теплообмена и тепловые схемы

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

Пластина с конвекцией на правой границе. Аналитический эталон. b и h — размеры выделенной геометрии между крайними поверхностями по координатам; пропорции сечения сохранены.
Пластина с конвекцией на правой границе1 — Расчётная областьАналитический эталон. b и h — размеры выделенной геометрии между крайними поверхностями по координатам; пропорции сечения сохранены.Открыть SVG целикомРезультат отдельно
Предварительные понятия

Зачем нужны граничные условия

Уравнение теплопроводности описывает, что происходит внутри тела. Но нагрев обмотки или остывание стенки определяются прежде всего тем, как тело обменивается теплом с окружением: омывается ли оно воздухом или водой, излучает ли в холодное помещение, прижато ли к радиатору. Все эти процессы — конвекция, излучение, контакт с другим телом — в расчётной модели теплопроводности сводятся к граничным условиям на поверхности тела. Решение краевой задачи единственно, только если на каждой части границы задано ровно одно условие (теорема единственности); выбор правильного типа условия — главное инженерное решение при постановке тепловой задачи.

Обозначения: TsT_s — температура поверхности тела, n\un — внешняя нормаль к ней, qn=q⋅n=−λ ∂T/∂nq_n=\uq\cdot\un=-\lambda\,\partial T/\partial n — плотность теплового потока, выходящего из тела через поверхность.

Три рода граничных условий

Граничное условие первого рода (Дирихле)

На поверхности задана температура: T∣S=Ts(r,t)T\big|_S=T_s(\vec r,t). Физически это означает, что внешняя среда настолько интенсивно обменивается теплом с поверхностью, что навязывает ей свою температуру: кипящая жидкость, конденсирующийся пар, массивный термостатированный радиатор. Плотность потока через поверхность при этом заранее неизвестна и находится из решения.

Граничное условие второго рода (Нейман)

На поверхности задана нормальная плотность потока: −λ ∂T/∂n∣S=qn(r,t)-\lambda\,\partial T/\partial n\big|_S=q_n(\vec r,t). Важнейший частный случай qn=0q_n=0 — теплоизолированная (адиабатная) поверхность. Такое же условие ставят на плоскости симметрии задачи: если распределение температуры симметрично относительно плоскости, то на ней ∂T/∂n=0\partial T/\partial n=0 и тепло через неё не переходит; это позволяет рассчитывать половину или четверть тела. Условие второго рода задаёт также известный подвод тепла — например, от плёночного нагревателя на поверхности.

Граничное условие третьего рода (конвекция, закон Ньютона–Рихмана)

Поверхность омывается средой с температурой T∞T_\infty вдали от неё; плотность потока, отдаваемого поверхностью среде, пропорциональна разности температур: −λ∂T∂n∣S=α (Ts−T∞).-\lambda\frac{\partial T}{\partial n}\Big|_S=\alpha\,(T_s-T_\infty). Коэффициент α\alpha — коэффициент теплоотдачи, Вт/(м²·К). Условие связывает температуру и её производную на границе и содержит два параметра среды: α\alpha и T∞T_\infty.

Коэффициент теплоотдачи — не свойство материала, а характеристика процесса обтекания: он зависит от скорости среды, её вязкости и теплопроводности, от ориентации и размера поверхности, от самой разности Ts−T∞T_s-T_\infty при свободной конвекции. Типичные порядки величины [Исаченко]:

Условия теплоотдачи α\alpha, Вт/(м²·К)
Воздух, свободная конвекция 5–25
Воздух, вынужденная конвекция (обдув) 25–250
Вода, вынужденная конвекция 500–10 000

Разброс в три порядка объясняет, почему водяное охлаждение позволяет снимать с той же поверхности в сотни раз большую мощность, чем естественное охлаждение воздухом.

Конвекция и излучение здесь — модели границы, а не расчёт течения

Условие третьего рода не решает задачу о движении жидкости и переносе тепла в ней: вся сложная гидродинамика пограничного слоя спрятана в одном числе α\alpha, взятом из опытных зависимостей (критериальных уравнений). Точно так же условие излучения, о котором речь ниже, не решает задачу переноса излучения между многими поверхностями с переотражениями — оно описывает обмен малой поверхности с большим окружением. Точность теплового расчёта определяется, как правило, не точностью решения уравнения теплопроводности, а точностью задания α\alpha: ошибка в α\alpha на 3030 % — обычное дело, и к ней нужно относиться как к неустранимой неопределённости модели.

Излучение

Любая поверхность с абсолютной температурой TsT_s излучает в окружение, имеющее температуру TокрT_{\text{окр}}, и поглощает излучение окружения. Для малого тела в большом окружении (стенки помещения, небо) результирующая плотность потока qизл=εи σСБ(Ts4−Tокр4),q_{\text{изл}}=\varepsilon_{\mathrm и}\,\sigma_{\text{СБ}}\left(T_s^4-T_{\text{окр}}^4\right), где εи\varepsilon_{\mathrm и} — степень черноты поверхности (0<εи≤10<\varepsilon_{\mathrm и}\le1; окисленный металл 0,60{,}6–0,90{,}9, краска 0,90{,}9, полированный металл 0,050{,}05–0,10{,}1), σСБ=5,670⋅10−8\sigma_{\text{СБ}}=5{,}670\cdot10^{-8} Вт/(м²·К⁴) — постоянная Стефана–Больцмана [CODATA]. Температуры входят в четвёртой степени и только в кельвинах: подстановка градусов Цельсия даёт бессмыслицу, вплоть до отрицательного потока от более горячего тела к более холодному.

Нелинейное условие неудобно для схем замещения, поэтому его линеаризуют, представляя разность четвёртых степеней как произведение: Ts4−Tокр4=(Ts2+Tокр2)(Ts+Tокр)(Ts−Tокр),T_s^4-T_{\text{окр}}^4=\left(T_s^2+T_{\text{окр}}^2\right)\left(T_s+T_{\text{окр}}\right)\left(T_s-T_{\text{окр}}\right), откуда qизл=αизл (Ts−Tокр),αизл=εи σСБ(Ts2+Tокр2)(Ts+Tокр).q_{\text{изл}}=\alpha_{\text{изл}}\,(T_s-T_{\text{окр}}),\qquad \alpha_{\text{изл}}=\varepsilon_{\mathrm и}\,\sigma_{\text{СБ}}\left(T_s^2+T_{\text{окр}}^2\right)\left(T_s+T_{\text{окр}}\right). Величина αизл\alpha_{\text{изл}} — «коэффициент теплоотдачи излучением» — зависит от обеих температур, поэтому её приходится уточнять итерациями: задаться TsT_s, вычислить αизл\alpha_{\text{изл}}, решить линейную задачу, получить новое TsT_s и повторить. Так же поступают и в численных расчётах; лаборатория ниже решает нелинейное условие итерационно и приводит Rизл=1/(αизлS)R_{\text{изл}}=1/(\alpha_{\text{изл}}S) при найденной TsT_s.

Коэффициент теплоотдачи излучением

Окрашенный корпус, εи=0,9\varepsilon_{\mathrm и}=0{,}9, при Ts=350T_s=350 К в помещении с Tокр=300T_{\text{окр}}=300 К: αизл=0,9⋅5,670⋅10−8⋅(3502+3002)(350+300)=0,9⋅5,670⋅10−8⋅2,125⋅105⋅650≈7,05 Вт/(м2⋅К).\alpha_{\text{изл}}=0{,}9\cdot5{,}670\cdot10^{-8}\cdot(350^2+300^2)(350+300)=0{,}9\cdot5{,}670\cdot10^{-8}\cdot2{,}125\cdot10^5\cdot650\approx7{,}05\ \text{Вт/(м}^2\cdot\text{К)}. Проверка по исходной формуле: qизл=0,9⋅5,670⋅10−8⋅(1,5006⋅1010−8,1⋅109)=352q_{\text{изл}}=0{,}9\cdot5{,}670\cdot10^{-8}\cdot(1{,}5006\cdot10^{10}-8{,}1\cdot10^{9})=352 Вт/м², и 352/50=7,05352/50=7{,}05 Вт/(м²·К). При свободной конвекции в воздухе (α≈5\alpha\approx5–1010) излучение переносит столько же или больше тепла, чем конвекция, — пренебрегать им при естественном охлаждении нельзя.

Если тело одновременно омывается воздухом с T∞T_\infty и излучает в окружение с TокрT_{\text{окр}}, потоки складываются: qn=α (Ts−T∞)+εиσСБ(Ts4−Tокр4).q_n=\alpha\,(T_s-T_\infty)+\varepsilon_{\mathrm и}\sigma_{\text{СБ}}\left(T_s^4-T_{\text{окр}}^4\right). При T∞=TокрT_\infty=T_{\text{окр}} это условие третьего рода с суммарным коэффициентом αΣ=α+αизл\alpha_\Sigma=\alpha+\alpha_{\text{изл}}.

Условия на границе раздела и контактное сопротивление

На границе двух тел с теплопроводностями λ1\lambda_1, λ2\lambda_2 выполняются два условия (см. границы раздела сред): непрерывность нормальной плотности потока — энергия на границе не накапливается — и непрерывность температуры при идеальном контакте: λ1∂T1∂n=λ2∂T2∂n,T1=T2.\lambda_1\frac{\partial T_1}{\partial n}=\lambda_2\frac{\partial T_2}{\partial n},\qquad T_1=T_2 . Реальные поверхности шероховаты, соприкасаются лишь отдельными выступами, а в промежутках заключён воздух с его ничтожной теплопроводностью. Поэтому при одинаковом потоке температура на границе испытывает скачок: T1−T2=Rк qn,T_1-T_2=R_{\mathrm к}\,q_n , где RкR_{\mathrm к} — удельное контактное тепловое сопротивление, К·м²/Вт. Его типичные значения 10−510^{-5}–10−310^{-3} К·м²/Вт; оно уменьшается при росте давления прижима, шлифовке поверхностей, применении теплопроводящих паст и прокладок. При qn=104q_n=10^4 Вт/м² и Rк=2⋅10−4R_{\mathrm к}=2\cdot10^{-4} К·м²/Вт скачок составляет 22 К — для силового полупроводникового прибора на радиаторе это заметная величина.

Тепловые сопротивления и схемы замещения

В стационарном одномерном случае через каждый слой и каждую граничную поверхность идёт один и тот же тепловой поток PP, а температура падает на величину, пропорциональную PP. Это полная аналогия с электрической цепью: PP играет роль тока, перепад температур — роль напряжения.

Тепловое сопротивление

Rт=ΔT/PR_{\mathrm т}=\Delta T/P, К/Вт. Из одномерных решений предыдущего урока и граничных условий: Rпл=dλS,Rцил=ln⁡(r2/r1)2πλℓ,Rконв=1αS,Rизл=1αизлS,Rконт=RкS.R_{\text{пл}}=\frac{d}{\lambda S},\qquad R_{\text{цил}}=\frac{\ln(r_2/r_1)}{2\pi\lambda\ell},\qquad R_{\text{конв}}=\frac{1}{\alpha S},\qquad R_{\text{изл}}=\frac{1}{\alpha_{\text{изл}}S},\qquad R_{\text{конт}}=\frac{R_{\mathrm к}}{S}. Здесь SS — площадь, через которую идёт поток (для цилиндра — площадь той поверхности, на которой задано условие), ℓ\ell — длина трубы. Сопротивления слоёв, стоящих на пути одного потока, складываются (последовательное соединение); сопротивления параллельных путей — например, конвекции и излучения с одной поверхности — складываются по обратным величинам (параллельное соединение).

Для многослойной плоской стенки между двумя средами полное сопротивление на единицу площади 1k=1α1+∑idiλi+∑jRкj+1α2,q=k (T∞1−T∞2),\frac1k=\frac1{\alpha_1}+\sum_i\frac{d_i}{\lambda_i}+\sum_j R_{\mathrm к j}+\frac1{\alpha_2},\qquad q=k\,(T_{\infty1}-T_{\infty2}), где kk — коэффициент теплопередачи, Вт/(м²·К). Для цилиндрической стенки удобнее работать с сопротивлениями на единицу длины R′=R ℓR'=R\,\ell, К·м/Вт, и потоком P′P', Вт/м. Температура на каждой границе находится последовательным вычитанием P RiP\,R_i: это и есть «расчёт цепи».

Проверка теплового баланса

В стационарном режиме для любого замкнутого объёма — слоя, всей стенки, всего устройства — выполняется Pвх−Pвых+Pвыд=0,P_{\text{вх}}-P_{\text{вых}}+P_{\text{выд}}=0 , где Pвыд=∫VqV dVP_{\text{выд}}=\int_Vq_V\,\dV — мощность внутренних источников. Без источников тепловой поток через все сечения одинаков. Этот баланс — первая проверка любого расчёта, аналитического или численного: если он нарушен, ошибка либо в граничных условиях, либо в постановке задачи. В численных решениях невязка баланса, отнесённая к PвхP_{\text{вх}}, служит мерой погрешности.

Трёхслойная стена с конвекцией с обеих сторон

Стена: бетон d1=0,15d_1=0{,}15 м (λ1=1,5\lambda_1=1{,}5), минеральная вата d2=0,10d_2=0{,}10 м (λ2=0,04\lambda_2=0{,}04), бетон d3=0,05d_3=0{,}05 м (λ3=1,5\lambda_3=1{,}5). Внутри воздух T∞1=20T_{\infty1}=20 °C, свободная конвекция α1=8\alpha_1=8; снаружи T∞2=−20T_{\infty2}=-20 °C, ветер, α2=25\alpha_2=25 Вт/(м²·К). Сопротивления на 1 м²: 1α1=0,125,d1λ1=0,100,d2λ2=2,500,d3λ3=0,033,1α2=0,040  К⋅м2/Вт;\frac1{\alpha_1}=0{,}125,\quad \frac{d_1}{\lambda_1}=0{,}100,\quad \frac{d_2}{\lambda_2}=2{,}500,\quad \frac{d_3}{\lambda_3}=0{,}033,\quad \frac1{\alpha_2}=0{,}040\ \ \text{К·м}^2\text{/Вт}; сумма 1/k=2,7981/k=2{,}798, k=0,357k=0{,}357 Вт/(м²·К). Плотность потока q=40/2,798=14,3q=40/2{,}798=14{,}3 Вт/м². Температуры границ, начиная от внутреннего воздуха: 20−14,3⋅0,125=18,220-14{,}3\cdot0{,}125=18{,}2 °C (внутренняя поверхность), далее 16,816{,}8 °C, −19,0-19{,}0 °C, −19,4-19{,}4 °C и, наконец, −19,4−14,3⋅0,040=−20,0-19{,}4-14{,}3\cdot0{,}040=-20{,}0 °C — вернулись к температуре наружного воздуха: баланс сошёлся. Почти весь перепад (35,735{,}7 К из 4040) приходится на утеплитель, который даёт 8989 % сопротивления; бетонные слои и конвекция вместе — лишь 1111 %.

Изолированная труба

Стальная труба r1=25r_1=25 мм, r2=30r_2=30 мм (λ=50\lambda=50), покрыта минеральной ватой до r3=60r_3=60 мм (λи=0,04\lambda_{\mathrm и}=0{,}04). Внутри вода 150150 °C, α1=2000\alpha_1=2000; снаружи воздух 2020 °C, α2=10\alpha_2=10 Вт/(м²·К). Сопротивления на 1 м длины: R1′=12π⋅0,025⋅2000=0,0032,Rст′=ln⁡(30/25)2π⋅50=0,0006,Rи′=ln⁡(60/30)2π⋅0,04=2,758,R2′=12π⋅0,06⋅10=0,265  К⋅м/Вт.R'_1=\frac1{2\pi\cdot0{,}025\cdot2000}=0{,}0032,\quad R'_{\text{ст}}=\frac{\ln(30/25)}{2\pi\cdot50}=0{,}0006,\quad R'_{\mathrm и}=\frac{\ln(60/30)}{2\pi\cdot0{,}04}=2{,}758,\quad R'_2=\frac1{2\pi\cdot0{,}06\cdot10}=0{,}265\ \ \text{К·м/Вт}. Сумма 3,0273{,}027 К·м/Вт, P′=130/3,027=42,9P'=130/3{,}027=42{,}9 Вт/м. Температуры: вода 150150, внутренняя стенка 149,9149{,}9, наружная стенка трубы 149,8149{,}8, поверхность изоляции 31,431{,}4 °C. Без изоляции R2′=1/(2π⋅0,03⋅10)=0,531R'_2=1/(2\pi\cdot0{,}03\cdot10)=0{,}531, сумма 0,5340{,}534, P′=243P'=243 Вт/м — изоляция снизила потери в 5,75{,}7 раза, а стальная стенка и внутренняя теплоотдача практически не сопротивляются потоку.

Критический радиус изоляции

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

Критический радиус

Пусть труба с наружным радиусом r1r_1 покрыта изоляцией (λ\lambda) до радиуса rr, снаружи — теплоотдача α\alpha. Сопротивление изоляции и теплоотдачи на единицу длины: R′(r)=ln⁡(r/r1)2πλ+12πrα.R'(r)=\frac{\ln(r/r_1)}{2\pi\lambda}+\frac{1}{2\pi r\alpha}. Дифференцируем: dR′dr=12πλr−12παr2=12πr2(rλ−1α)=0  ⇒  rкр=λα.\frac{\dd R'}{\dd r}=\frac{1}{2\pi\lambda r}-\frac{1}{2\pi\alpha r^2}=\frac{1}{2\pi r^2}\left(\frac r\lambda-\frac1\alpha\right)=0 \;\Rightarrow\; r_{\text{кр}}=\frac{\lambda}{\alpha}. При r<rкрr<r_{\text{кр}} производная отрицательна: сопротивление убывает с ростом толщины изоляции, и потери растут. При r>rкрr>r_{\text{кр}} изоляция начинает работать как изоляция. Таким образом, rкрr_{\text{кр}} — радиус наибольших потерь. Критерий простой: если наружный радиус трубы r1r_1 уже больше λ/α\lambda/\alpha, любая изоляция уменьшает потери; если меньше — сначала потери возрастут.

Для минеральной ваты в воздухе при α=10\alpha=10: rкр=0,04/10=4r_{\text{кр}}=0{,}04/10=4 мм. Для труб это неактуально — их радиусы заведомо больше. Но для электрических проводов в изоляции из полимера (λ≈0,2\lambda\approx0{,}2, rкр=20r_{\text{кр}}=20 мм при α=10\alpha=10) эффект работает в полезную сторону: слой изоляции на тонком проводе улучшает его охлаждение, увеличивая поверхность теплоотдачи. Безразмерная форма условия, αr/λ<1\alpha r/\lambda<1, — это число Био, меньшее единицы; с ним мы встретимся в нестационарных задачах.

Лаборатория

Лаборатория рассчитывает многослойную стенку — плоскую или цилиндрическую (труба с изоляцией) — до трёх слоёв. Для каждого слоя задаются толщина, λ\lambda и объёмное тепловыделение qVq_V; между слоями можно ввести контактное сопротивление RкR_{\mathrm к}. На левой и правой границах выбирается тип условия: заданная температура, конвекция (α\alpha, T∞T_\infty), излучение (εи\varepsilon_{\mathrm и}, TокрT_{\text{окр}}), конвекция вместе с излучением (при T∞=TокрT_\infty=T_{\text{окр}}), заданная плотность потока или теплоизоляция. Все температуры вводятся в °C; в формулу излучения подставляются кельвины. Выводятся температуры на всех границах слоёв, qq на наружной границе, тепловые потоки PP через обе границы (положительные — наружу) при заданной площади SS (для трубы — длине ℓ\ell), выделяемая мощность, баланс «сумма потоков наружу − выделенное = 0», тепловые сопротивления слоёв, контактов и граничных условий, схема замещения и графики T(x)T(x), q(x)q(x). Решение аналитическое по слоям; нелинейное условие излучения решается итерационно, а RизлR_{\text{изл}} вычисляется по линеаризованному αизл\alpha_{\text{изл}} при найденной TsT_s.

Что стоит проделать:

Распространённые заблуждения

  • Градусы Цельсия в формуле излучения. T4T^4 требует абсолютной температуры. (50 °C)4(50\ °\text{C})^4 не имеет физического смысла, а разность четвёртых степеней в градусах Цельсия может дать даже неверный знак. В разности Ts−T∞T_s-T_\infty закона Ньютона–Рихмана — любые единицы.
  • «Коэффициент теплоотдачи — свойство материала стенки». α\alpha характеризует среду и режим её движения; материал стенки на него почти не влияет. Два одинаковых корпуса в спокойном воздухе и под вентилятором имеют α\alpha, различающиеся в десять раз.
  • Условие первого рода «вместо» конвекции. Задать на поверхности Ts=T∞T_s=T_\infty — значит предположить α→∞\alpha\to\infty. Для водяного охлаждения это иногда приемлемо, для воздушного — грубая ошибка, занижающая температуру тела на десятки кельвинов.
  • Два условия на одной границе. Нельзя одновременно задать и температуру, и поток: задача станет переопределённой. Поток на границе первого рода — результат решения, а не входная величина.
  • Сложение сопротивлений параллельных путей. Конвекция и излучение с одной поверхности действуют параллельно: складываются коэффициенты α+αизл\alpha+\alpha_{\text{изл}} (проводимости), а не сопротивления.
  • Пренебрежение излучением «потому что температуры низкие». При естественном охлаждении в воздухе αизл≈5\alpha_{\text{изл}}\approx5–77 Вт/(м²·К) сопоставим с конвективным α\alpha уже при 5050 К перегрева.

Вопросы и упражнения

  1. Пластина толщиной 2020 мм из стали (λ=50\lambda=50) с одной стороны обогревается потоком qn=2000q_n=2000 Вт/м², с другой омывается воздухом T∞=20T_\infty=20 °C, α=20\alpha=20 Вт/(м²·К). Найдите температуры обеих поверхностей. Какая доля перепада приходится на саму пластину?
  2. Для стены из примера найдите, во сколько раз уменьшится qq, если толщину минеральной ваты удвоить. Почему не вдвое?
  3. Поверхность Ts=400T_s=400 К излучает в окружение Tокр=300T_{\text{окр}}=300 К, εи=0,8\varepsilon_{\mathrm и}=0{,}8. Найдите qизлq_{\text{изл}} точно и через αизл\alpha_{\text{изл}}; убедитесь, что результаты совпадают. Какую ошибку даст подстановка градусов Цельсия (127127 и 2727)?
  4. Силовой транзистор выделяет P=50P=50 Вт и прижат к радиатору через площадь S=4S=4 см² с Rк=3⋅10−4R_{\mathrm к}=3\cdot10^{-4} К·м²/Вт. Найдите скачок температуры на контакте и тепловое сопротивление контакта в К/Вт. Как его уменьшить вдвое?
  5. Выведите критический радиус для сферического слоя изоляции и покажите, что он равен 2λ/α2\lambda/\alpha.
  6. Трёхслойная стенка без источников рассчитана численно, и программа выдала потоки на левой и правой границах 103,2103{,}2 и 98,798{,}7 Вт. Что это означает и какой параметр расчёта следует проверить в первую очередь?

Ответы к пунктам 1, 3 и 4

1. В стационарном режиме весь подводимый поток уходит в воздух: Ts2=T∞+qn/α=20+2000/20=120T_{s2}=T_\infty+q_n/\alpha=20+2000/20=120 °C. Перепад в пластине qnd/λ=2000⋅0,02/50=0,8q_nd/\lambda=2000\cdot0{,}02/50=0{,}8 К, Ts1=120,8T_{s1}=120{,}8 °C. На пластину приходится 0,8/100,8≈0,80{,}8/100{,}8\approx0{,}8 % перепада — всё определяется теплоотдачей.

3. Точно: q=0,8⋅5,670⋅10−8⋅(2,56⋅1010−8,1⋅109)=794q=0{,}8\cdot5{,}670\cdot10^{-8}\cdot(2{,}56\cdot10^{10}-8{,}1\cdot10^{9})=794 Вт/м². Через αизл=0,8⋅5,670⋅10−8⋅(1,6⋅105+9⋅104)⋅700=7,94\alpha_{\text{изл}}=0{,}8\cdot5{,}670\cdot10^{-8}\cdot(1{,}6\cdot10^5+9\cdot10^4)\cdot700=7{,}94 Вт/(м²·К), q=7,94⋅100=794q=7{,}94\cdot100=794 Вт/м². В градусах Цельсия: 1274−274=2,60⋅108−5,3⋅105≈2,6⋅108127^4-27^4=2{,}60\cdot10^8-5{,}3\cdot10^5\approx2{,}6\cdot10^8, что даёт q≈11,8q\approx11{,}8 Вт/м² — в 6767 раз меньше истинного.

4. qn=50/(4⋅10−4)=1,25⋅105q_n=50/(4\cdot10^{-4})=1{,}25\cdot10^5 Вт/м²; скачок ΔT=Rкqn=3⋅10−4⋅1,25⋅105=37,5\Delta T=R_{\mathrm к}q_n=3\cdot10^{-4}\cdot1{,}25\cdot10^5=37{,}5 К; Rконт=Rк/S=0,75R_{\text{конт}}=R_{\mathrm к}/S=0{,}75 К/Вт. Уменьшить вдвое — удвоить площадь контакта или вдвое снизить RкR_{\mathrm к} (теплопроводящая паста, большее усилие прижима).

Итоги

Далее — нестационарная теплопроводность: как быстро тело прогревается и когда его можно считать телом с одной температурой.

Применение в FEMM 5