Термодинаміка та теплопередача — Карно, випромінювання абсолютно чорного тіла та фазові діаграми

Чому термодинаміка неминуча

Кожен процес у Всесвіті — хімічний, біологічний, механічний, електромагнітний — зрештою керується термодинамікою. Термодинаміку не можна вивести із законів Ньютона; вона стоїть на власному логічному фундаменті, побудованому на спостереженнях про тепло, роботу та стрілу часу. Її закони прості, але їхні наслідки глибокі: машини не можуть бути ідеально ефективними, безлад зростає, а абсолютний нуль назавжди залишається недосяжним.

Термодинаміку унікально складно викладати саме тому, що її ключові поняття — ентропія, температура, теплоємність — це статистичні середні за величезною кількістю частинок. Жодна окрема частинка не має температури. Інтерактивна симуляція, що дозволяє спостерігати, як розподіли Максвелла-Больцмана перебудовуються, коли ви пересуваєте повзунок температури, або спостерігати, як спонтанно утворюються валики конвекції в нагрітій рідині, робить ці середні значення відчутними так, як цього не можуть самі лише рівняння.

Цей огляд проведе вас крізь шість термодинамічних симуляцій платформи, показуючи, чого навчає кожна з них і як вони пов'язані з базовою фізикою.

Шар 1: закон охолодження Ньютона

Симуляція охолодження кави

Найпростіший термодинамічний експеримент є водночас одним із найпоказовіших: гаряче тіло, вміщене в холодніше середовище. Закон охолодження Ньютона стверджує, що швидкість втрати тепла пропорційна різниці температур між об'єктом і навколишнім середовищем. Результат — експоненційний спад — та сама математична форма, що керує радіоактивним розпадом, зарядкою RC- кола та зростанням популяції до межі ємності середовища.

Охолодження Ньютона:
  dT/dt = -k(T - T_середовища)
  Розв'язок: T(t) = T_середовища + (T₀ - T_середовища)·e^{-kt}
  k = h·A / (m·c_p)       [константа охолодження, с⁻¹]

Теплопровідність Фур'є (1D):
  q = -λ · dT/dx           [Вт/м², тепловий потік]
  ∂T/∂t = α · ∂²T/∂x²     [рівняння теплопровідності, α = λ/(ρ·cₚ)]

Температуропровідність:
  α = λ / (ρ·cₚ)           [м²/с]
  Мідь: α ≈ 1.17×10⁻⁴ м²/с
  Вода:  α ≈ 1.43×10⁻⁷ м²/с (у 820 разів повільніше)

Закон випромінювання Ньютона (Стефан-Больцман):
  P = ε·σ·A·(T⁴ - T_сер⁴)   [домінує при високій T]
  σ = 5.67×10⁻⁸ Вт/(м²·К⁴)

Симуляція охолодження кави дозволяє налаштувати початкову температуру, температуру середовища, матеріал чашки (теплопровідність) і площу поверхні. Спостерігайте, як крива охолодження зсувається від майже лінійної при малому ΔT до крутого експоненційного спаду при великому ΔT. Перемикайтеся між чистою теплопровідністю та охолодженням, де домінує випромінювання, щоб побачити, чому член T⁴ Стефана-Больцмана має значення лише при дуже високих температурах.

Шар 2: межа теплового двигуна

Симуляція циклу Карно

У 1824 році Саді Карно довів дещо неймовірне без жодних знань про атоми: жоден тепловий двигун, що працює між двома резервуарами з температурами T_Г та T_Х, ніколи не може перевищити ефективність Карно η = 1 − T_Х/T_Г. Ця межа абсолютна. Вона не залежить від робочої рідини, конструкції двигуна чи винахідливості інженера — вона випливає лише з другого закону.

Процес           Тип                Тепло      Робота
─────────────────────────────────────────────────────
1→2  Ізотермічне розширення     Q_Г > 0    W₁₂ = Q_Г
2→3  Адіабатичне розширення     Q = 0      W₂₃ = -ΔU
3→4  Ізотермічне стиснення      Q_Х < 0    W₃₄ = Q_Х
4→1  Адіабатичне стиснення      Q = 0      W₄₁ = -ΔU

Чистий вихід роботи:
  W_нетто = Q_Г - |Q_Х|

Ефективність Карно:
  η_К = W_нетто/Q_Г = 1 - T_Х/T_Г   (температури в Кельвінах)

Обмеження другого закону:
  ΔS_всесвіту ≥ 0
  |Q_Х|/T_Х ≥ Q_Г/T_Г  (нерівність Клаузіуса)

Приклади реальних двигунів (η_К проти реальної):
  Парова турбіна: T_Г=600К, T_Х=300К → η_К=50%, реально ~42%
  Двигун автомобіля: T_Г=900К, T_Х=300К → η_К=67%, реально ~25-35%

Симуляція циклу Карно малює повну P-V діаграму та T-S діаграму одночасно. Перетягуйте повзунки T_Г і T_Х та спостерігайте, як обмежена площа — чиста робота — розширюється чи стискається. Порівняйте з циклом Отто (ідеалізований бензиновий двигун), щоб побачити, чому дизельні двигуни з вищим ступенем стиснення досягають кращої реальної ефективності.

Шар 3: швидкість молекул

Розподіл Максвелла-Больцмана

Температура — це середня кінетична енергія молекул, але не всі молекули рухаються з однаковою швидкістю. Джеймс Клерк Максвелл і Людвіг Больцман вивели розподіл ймовірностей швидкостей молекул в ідеальному газі. Цей розподіл не гаусів: він має довгий хвіст високих швидкостей, що пояснює, чому хімічні реакції можливі навіть тоді, коли середня молекула не має достатньої енергії.

Густина ймовірності:
  f(v) = 4π·n·(m/2πk_BT)^{3/2}·v²·exp(-mv²/2k_BT)

Характерні швидкості:
  v_p  = √(2k_BT/m)          найімовірніша
  v̄   = √(8k_BT/πm)         середня
  v_rms = √(3k_BT/m)         середньоквадратична

Співвідношення: v_p < v̄ < v_rms  (пропорція ≈ 1 : 1.128 : 1.225)

Теорема про рівнорозподіл:
  ½m⟨v_x²⟩ = ½m⟨v_y²⟩ = ½m⟨v_z²⟩ = ½k_BT
  Одноатомний ідеальний газ: U = (3/2)Nk_BT
  Двоатомний ідеальний газ:  U = (5/2)Nk_BT  (+ 2 обертальні ступені вільності)

Фактор Больцмана:
  P(E) ∝ exp(-E/k_BT)         Арреніус: k = A·exp(-E_a/k_BT)

Симуляція Максвелла-Больцмана дозволяє встановити температуру газу й молекулярну масу, спостерігаючи, як зсувається крива розподілу, а жива гістограма змодельованих швидкостей молекул узгоджується з теоретичною кривою. Хвіст високих швидкостей стає відчутним: навіть при помірних температурах невелика частка молекул має достатньо енергії, щоб покинути рідку фазу (випаровування) або подолати бар'єри реакції (хімія).

Шар 4: світло від гарячих об'єктів

Симуляція випромінювання абсолютно чорного тіла

Кожен об'єкт вище абсолютного нуля випромінює електромагнітне випромінювання. Спектр цього випромінювання залежить лише від температури об'єкта — факт, який класична фізика не могла пояснити, породивши "ультрафіолетову катастрофу". Макс Планк вирішив цю кризу 1900 року, припустивши, що енергія квантована в дискретні пакети — революційна гіпотеза, що започаткувала квантову механіку.

Спектральна яскравість Планка:
  B(λ,T) = (2hc²/λ⁵) · 1/(exp(hc/λk_BT) - 1)

Закон зсуву Віна:
  λ_max · T = 2.898×10⁻³ м·К
  Сонце (T≈5778К):  λ_max ≈ 502 нм (зелений)
  Тіло людини (310К): λ_max ≈ 9.35 мкм (інфрачервоний)

Закон Стефана-Больцмана:
  P = ε·σ·A·T⁴
  σ = 5.670×10⁻⁸ Вт/(м²·К⁴)
  Подвоєння T → у 16 разів більше потужності

Класичний Релей-Джинс (не працює на коротких λ):
  B_РД(λ,T) = 2ck_BT/λ⁴    → ∞ при λ → 0 (ультрафіолетова катастрофа)

Виправлення Планка: квантовані осцилятори
  ⟨E⟩ = hν/(exp(hν/k_BT) - 1)  замість k_BT

Симуляція випромінювання абсолютно чорного тіла показує спектр Планка поряд із класичним передбаченням Релея-Джинса. Пересуньте температуру від 300 К до 30 000 К і дивіться, як пікова довжина хвилі проходить від інфрачервоної через видимий червоний, жовтий і білий — точно ті самі кольори, що й у розжареного металу. Класична крива катастрофічно розходиться на коротких довжинах хвиль, тоді як формула Планка залишається фізичною, роблячи квантову революцію відчутно неминучою.

Шар 5: спонтанне утворення візерунків

Симуляція конвекції Бенара

Нагрійте плаский шар рідини знизу, і він залишиться провідним — тепло дифундує вгору без руху рідини — поки не буде перевищено критичний поріг, число Релея. Вище цього порога стаціонарний стан стає нестійким, і спонтанно виникають валики конвекції, утворюючи прекрасні шестикутні клітини Бенара, що спостерігаються в нагрітій олії або гранулах поверхні Сонця.

Число Релея:
  Ra = g·β·ΔT·d³ / (ν·α)
  g  = гравітація [м/с²]
  β  = коефіцієнт теплового розширення [1/К]
  ΔT = різниця температур [К]
  d  = товщина шару [м]
  ν  = кінематична в'язкість [м²/с]
  α  = температуропровідність [м²/с]

Число Прандтля:
  Pr = ν/α   (відношення дифузійності імпульсу до теплової)
  Вода: Pr ≈ 7    Повітря: Pr ≈ 0.71    Силіконова олія: Pr ≈ 1000

Початок конвекції:
  Ra_к ≈ 1708  (для жорстко-жорстких граничних умов)

Довжина хвилі валика на початку:
  λ ≈ 2d  (горизонтальні клітини приблизно тієї самої ширини, що й товщина шару)

Число Нуссельта (посилення конвекції):
  Nu = Q_конв/Q_пров = 1 + C·(Ra/Ra_к - 1)^n   для Ra > Ra_к

Симуляція конвекції Бенара використовує розв'язувач рівнянь Нав'є-Стокса методом скінченних різниць із наближенням Бусінеска. Збільшіть градієнт температури понад критичне число Релея й дивіться, як безлад раптово кристалізується в структури, що обертаються. Збільшіть далі, і валики стануть турбулентними — яскрава демонстрація того, як порядок виникає з граничних умов, а не з проєкту.

Шар 6: матеріали та фазові переходи

Симуляція фазової діаграми сплаву

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

Області фазової діаграми:
  Рідина          — однорідний розплав
  Твердий розчин L (α)   — твердий розчин, багатий на A
  Твердий розчин R (β)   — твердий розчин, багатий на B
  Рідина + α      — співіснування двох фаз
  Рідина + β      — співіснування двох фаз
  α + β           — субсолідусна двофазна область

Правило важеля (частка фаз у двофазній області):
  f_α  = (C_β - C₀) / (C_β - C_α)
  f_β  = (C₀ - C_α) / (C_β - C_α)
  C₀ = загальний склад, C_α, C_β = межі фаз

Евтектична точка:
  Мінімальна температура плавлення на фазовій діаграмі
  В евтектиці: рідина → α + β одночасно
  Евтектика припою Sn-Pb: 63% Sn, 37% Pb, T_e = 183°C

Правило фаз Гіббса:
  F = C - P + 2   (C компонентів, P фаз, F ступенів вільності)
  В евтектиці:    F = 2 - 3 + 2 = 1  (зафіксувавши P, T визначено)

Симуляція фазової діаграми сплаву відображає повну бінарну евтектичну діаграму та відстежує ваш курсор через фазові області з живим показанням правила важеля. Встановіть початковий склад рідини, задайте швидкість охолодження й дивіться, як формується мікроструктура: первинні дендрити α кристалізуються першими, збагачуючи решту рідини компонентом B, поки вона не досягне евтектики й не замерзне одночасно.

Повна термодинамічна колекція

Шість симуляцій охоплюють основні розділи класичної термодинаміки:

Зв'язки з іншими колекціями

Термодинаміка проникає майже в кожну іншу категорію на платформі. Той самий розподіл Максвелла-Больцмана, що керує кінетикою газу, також визначає константи швидкості в симуляції ферментативної кінетики через рівняння Арреніуса. Рівняння теплопровідності Фур'є математично ідентичне рівнянню дифузії, що керує симуляцією дифузії в кристалах — замініть температуру на концентрацію, і фізика буде та сама. Конвекція Бенара — близька родичка симуляції атмосферної конвекції, де ефект Коріоліса порушує симетрію конвекційних клітин, перетворюючи їх на циклони. Навіть модель Ізінга — це термодинамічна система: її фазовий перехід при температурі Кюрі — підручниковий приклад фазового переходу другого роду, керованого фактором Больцмана.

Алгоритми та методи в цій колекції