Рівняння Фурієвого Тепла: Дифузія, Граничні Умови та Спадання Режиму

Тепло поширюється через градієнти

Рівняння теплопровідності формалізує просте фізичне правило: тепло тече від гарячого до холодного з швидкістю, пропорційною нахилу температурного градієнту. У одні dimensioni, вздовж тонкої жердини, температура u(x,t) підпорядковується рівнянню теплопровідності з константою пропорційності alpha, теплопередаваність матеріалу - як швидко воно вирівнює температурні відмінності, об'єднуючи теплопровідність, щільність і піжвідійсний капітальний прилад у одному числі.

∂u/∂t = α · ∂²u/∂x²          heat equation
finite-difference discretisation (explicit scheme):
u[i,t+1] = u[i,t] + α·Δt/Δx² · (u[i+1,t] - 2·u[i,t] + u[i-1,t])
stability requires: α·Δt/Δx² ≤ 0.5   (1D explicit / FTCS scheme)

Розв'язання за допомогою кінцевих різниць

Підхід з кінцевими різницями замінює безперервний стержень ланцюгом дискретних точок і замінює друге просторове похідне різницею між сусідами - другим порядком централізованої різниці, показаної вище. Поширення цього правила оновлення вперед у часі називається явним (FTCS, forward-time central-space) схемою: дешева, проста в реалізації, але лише обмежує стабільність. Якщо розмір кроку часу занадто великий відносно розміру кроку простору, невеликі чисельні помилки збільшуються на кожному кроці замість того, щоб зменшуватися, і зігноване охолодження змодельованого стержня виробляє нерівномірні, нефізичні температури замість плавної дифузії. Обмеження стабільності є прямим наслідком аналізу стабільності фон Неймана, який перевіряє, чи зростають чи зменшуються ф’юрієрові режими помилок під правилом оновлення.

Межальні умови визначають довгострокову форму

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

Почему аналитическое решение представляет собой сумму убывающих синусоидальных волн

Для стержня длиной L с обеими концами, удерживаемыми на нуле (Дирихле), общее решение уравнения теплопроводности выглядит следующим образом: u(x,t) = ∑(n=1 до ∞) Bₙ sin(nπx/L) * exp(-α(nπ/L)²*t)

Bₙ = коэффициенты Фурье синусной функции начального температурного профиля u(x,0)

Каждый член представляет собой пространственную синусоидальную моду с убывающей экспоненциально амплитудой и, что критически важно, более высокие частотные модули (большие n) затухают быстрее - скорость затухания масштабируется как n в квадрате. Это именно то, почему диффузия выравнивает вещи: острые, высокочастотные температурные всплески (большие n) почти мгновенно исчезают, а медленное, широкое, низко-n колебание сохраняется, поэтому любой начальный температурный профиль - каким бы неровным он ни был - сглаживается в плавную выпуклость и затем почти плоскую линию по мере увеличения t. Эта же математическая структура (разделение переменных плюс Фурье разложение) используется для решения уравнения волны и уравнения Шредингера, что объясняет, почему уравнение теплопроводности часто является первым PDE, которое преподают: это самый простой случай техники, которая повторяется во всей математической физике.

u(x,t) = Σ (n=1 to ∞)  Bₙ · sin(nπx/L) · exp(-α·(nπ/L)²·t)
Bₙ = Fourier sine coefficients of the initial temperature profile u(x,0)

Поза однією розмірністю та чому імпліцитні схеми є актуальними

У двох або трьох вимірах одна й та ж сама формула регулює поширення тепла через пластину чи твердий матеріал, і обмеження стабільності для явних кінцевих різниць стає ще більш суворим (на пару факторів більше, залежно від кількості просторових вимірів), тому код виробництва для будь-яких демонстрацій крім найпростіших зазвичай переходить на імпліцитні схеми (як Crank-Nicolson), які розв’язують лінійну систему на кожному кроці, але залишаються стабільними для будь-якого розміру часу кроку - обмінюючи більшу вартість обчислень за крок на свободу робити значно більші та, отже, менше кроків загалом.

Frequently asked questions

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

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