Розподіл температури: Як тепло поширюється через тверде тіло

Рівняння, що визначає поширення тепла

Тепло передається через провідність у твердих тілах: енергія перескакує від більш гарячих, високоенергетично вібруючих атомів до холодніших сусідніх без будь-якого руху матеріалу. Закон Ньютона-Фор’є (1822) встановлює, що місцевий тепловий потік пропорційний місцевому температурному градієнту — тепло тече від гарячого до холодного, і швидше там, де температура змінюється крутіше:

q = -k · ∇T (Закон Ньютона-Фор’є: тепловий потік проти температурного градієнта) ∂T/∂t = α · ∇²T (Рівняння тепла, що об'єднує закон Ньютона-Фор’є з законами збереження енергії) k = теплопровідність (Вт/м·К) α = k / (ρcp) = тепловідмінність (м2/с) демонстрація в реальному часі · поле температури, що розслаблюється до стаціонарного стану вздовж стержня● LIVE Об'єднуючи закон Ньютона-Фор’є з законами збереження енергії, отримуємо рівняння тепла, частковий диференціальний рівняння, яке описує еволюцію поля температури T(x,t) як у просторі, так і в часі. Воно математично ідентичне формулюванняму дифузійного рівняння, що керує розповсюдженням барвника через нерухому воду або розширенням газу в кімнаті — тепло, в цьому сенсі, буквально дифундує.

q = -k · ∇T          (Fourier's law: heat flux vs. temperature gradient)

∂T/∂t = α · ∇²T      (the heat equation, combining Fourier's law with energy conservation)

k = thermal conductivity (W/m·K)
α = k / (ρc_p) = thermal diffusivity (m²/s)

Стабільний стан проти тимчасового: дві дуже різні задачі

Коли кінці стержня утримуються на фіксованих температурах достатньо довго, температурний розподіл перестане змінюватися з часом — термін ∂T/∂t зникає — залишаючи значно простішу рівняння Лапласа, ∇²T = 0. Для простого 1D стержня без внутрішніх джерел тепла це змушує стабільний стан температури змінюватися лінійно між двома кінцевими температурами, незалежно від провідності матеріалу — провідність впливає на те, як швидко потік тепла та скільки часу потрібно для досягнення рівноваги, але не формує кінцевої стабільної форми. Проблема тимчасового характеру — як змінюється температурний розподіл у дорозі до цієї рівноваги, починаючи з деякого початкового розподілу — це складніша, часова залежна версія, і саме тут важливе повне рівняння теплопровідності з її терміном ∂T/∂t.

Аналітичні рішення та їх межі

Для простих геометрій – нескінченної стерні, плити з фіксованими температурами на межах, сфери, що охолоджується в однорідному середовищі – рівняння теплопровідності має точні аналітичні розв’язки, часто будуються за допомогою ряду Фур’є (іронічно повний цикл, оскільки Джозеф Фуріє розробив свій знаменитий ряд саме для вирішення цього рівняння) або з використанням функції помилки для напівнескінченних областей. Ці рішення є безцінними для валідації, але працюють лише для ідеалізованих форм із простими, незмінними граничними умовами. Реальні компоненти – теплообмінник із змінною перерізною площиною, стерня, що також передає тепло до навколишнього повітря, об’єкт із змінним матеріалом – потребують числових методів.

Метод кінцевих різниць: перетворення числення в арифметику

Метод кінцевих різниць (МЕР) замінює безперервний стержень на дискретну ланцюжок точок, розділених кроком Δx, і замінює другий похідний член у рівнянні теплопровідності простим середнім значенням із трьох точок:

∂²T/∂x² ≈ ( T[i+1] - 2*T[i] + T[i-1] ) / Δx² явне оновлення: T_new[i] = T[i] + α * Δt / Δx² * ( T[i+1] - 2*T[i] + T[i-1]) стабільність вимагає: α * Δt / Δx² ≤ 0.5 (1D явне МЕР) Це межеве значення стабільності не є рекомендацією — цей явний метод справді перетворюється на осцилюющі безглуздість, якщо розмір кроку часу занадто великий відносно просторового розрізу, як той самий CFL-умова, що обмежує кроки часу в симуляціях рідини та хвиль. Удвоєння просторового розрізу зменшує максимальний стабільний крок часу вчетверо, оскільки Δx з'являється в квадраті в межі, тому точні транзиентні теплові симуляції можуть швидко стати дорогими. Імплицитні схеми обмінюються більш дорогою лінійною розв’язкою на кожному кроці за безперечних умов стабільності, дозволяючи інженерам робити значно більші кроки часу лише тоді, коли важлива остання стабільна поведінка та точний часовий розріз транзиєнту не є проблемою.

∂²T/∂x²  ≈  ( T[i+1] - 2·T[i] + T[i-1] ) / Δx²

explicit update:
T_new[i] = T[i] + α · Δt/Δx² · ( T[i+1] - 2·T[i] + T[i-1] )

stability requires:  α · Δt / Δx²  ≤  0.5   (1D explicit FDM)

Рашотки: чо теплоносії мають таку форму

Охолодна рашотка втрачає тепло як шляхом проведення його вздовж власної довжини, так і конвекцією від поверхні в навколишній повітряний простір по всій довжині, що робить рівняння для рашки додатковий термін, який є більш складним, ніж звичайний термічний розрахунок. Розв’язування показує, що температура падає приблизно експоненційно від основи рашки, що і пояснює, чому теплоносії роблять тонкими та численними, а не кількома та товстими: більша площа поверхні на одиницю об’єму переносить більше тепла в повітря, хоча тонша рашотка менше проводить тепло вздовж власної довжини.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Temperature Distribution і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.