Równanie ciepła Fouriera: Dyfuzja, Warunki Brzegowe i Zanik Trybów

Rozpływ ciepła zależy od gradientu temperatury.

Równanie przewodzenia ciepła formalizuje prostą zasadę fizyczną: ciepło przepływa z gorącego do zimnego w tempie proporcjonalnym do nachylenia się gradientu temperatury. W jednym wymiarze, wzdłuż cienkiej prętowej osi, temperatura u(x,t) obeys równanie przewodzenia ciepła z stałą proporcyjna alfa, czyli współczynnik dyfuzyjności termicznej materiału – jak szybko wyrównuje się różnice temperatur, łącząc przewodność cieplną, gęstość i pojemność cieplną w jednym liczniku.

∂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)

Rozwiązanie za pomocą różnic skończonych

Metoda różnic skończonych zastępuje ciągłą pręt z łańcuchem dyskretnych punktów i zastępuje drugą pochodną przestrzenną różnicą między sąsiadami – drugiego rzędu centralną różnicę przedstawioną powyżej. Postęp tego algorytmu w czasie nazywany jest schematem jawnym (FTCS, forward-time central-space): tani, prosty w implementacji, ale warunkowo stabilny. Jeśli krok czasowy jest zbyt duży w stosunku do kroku przestrzennego, małe błędy numeryczne wzmacniają się na każdym kroku zamiast zanikać, a symulowany pręt generuje oscylujące, nieprawdopodobne temperatury zamiast gładkiego dyfuzyjnego rozpraszania. Granica stabilności jest bezpośrednim wynikiem analizy stabilności von Neumanna, która sprawdza, czy tryby błędów fourierowskich rosną lub maleją pod wpływem algorytmu aktualizacji.

Warunki brzegowe decydują o długoterminowym kształcie

To samo równanie wewnętrzne produkuje bardzo różne zachowanie w dłuższej perspektywie, w zależności od tego, co dzieje się na obu końcach pręta. Warunek Dirichlet ustala temperaturę na jednym końcu – jakby zamknąć koniec pręta w kąpieli lodowej – a ciepło ciągle odpływa lub napływa przez ten koniec. Warunek Neumann ustala przepływ ciepła (zerowy dla idealnie izolowanego końca) zamiast samej temperatury, więc nie ma przepływu ciepła przez tę granicę i temperatura końcowa jest wolna, by się poruszać, aż cały pręt osiągnie stan równowagi. Pręt ocieplony z obu stron zachowuje całkowitą energię cieplną i relaksuje się do pojedynczej, jednolitej temperatury – średniej wartości rozkładu początkowego, podczas gdy pręt zamknięty z obu stron przy różnych stałych temperaturach relaksuje się do profilu równowagi w linii prostej między nimi.

Dlaczego rozwiązanie analityczne jest sumą rozpadających się fal sinusoidalnych

Dla pręta o długości L, którego oba końce są utrzymywane w zerze (warunek Diraca), ogólne rozwiązanie równania ciepła ma postać szeregu Fouriera sinusów:

u(x,t) = ∑(n=1 do ∞) Bₙ sin(nπx/L) exp(-α(nπ/L)²t) Bₙ = współczynniki szeregu Fouriera sinusa dla profilu temperatury początkowego u(x,0) Każdy wyraz jest trybem sinusoidalnym przestrzennym, którego amplituda rozpadają się wykładniczo, i co najważniejsze, tryby o wyższej częstotliwości (większe n) rozpadają się szybciej – współczynnik zaniku skaluje się z n do kwadratu. To właśnie dlatego dyfuzja wygładza rzeczy: ostre, wysokiej częstotliwości szczyty temperatury (duże n) prawie natychmiast znikają, podczas gdy powolna, szeroka, o niskim n wariacja trwa, więc każdy początkowy profil temperatury – bez względu na to, jak nieregularny – zaokrągla się do gładkiego kłębusza i następnie niemal płaskiej linii w miarę wzrostu t. Jest to ta sama struktura matematyczna (podział zmiennych plus rozwinięcie Fouriera), która jest używana do rozwiązywania równania falowego i równania Schrödingera, dlatego równanie ciepła jest często pierwszym PME nauczanym: jest to najprostszy przypadek tej techniki, która powtarza się w całej matematycznej fizyce.

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)

Przestrzeń wielowymiarowa i przyczyny użycia schematów implicitnych

W dwóch lub trzech wymiarach ten sam równanie opisuje rozprzestrzenianie się ciepła przez płytę lub ciało, a ograniczenie stabilności dla metod różnic skończonych jest jeszcze bardziej restrykcyjne (zmniejszone o około czynnik równy liczbie wymiarów przestrzennych), co powoduje, że kody produkcyjne dla czegokolwiek poza prostymi demonstracjami zazwyczaj przechodzą na schematy implicitne (takie jak Crank-Nicolson). Te schematy rozwiązują liniowy układ równań w każdym kroku, ale zachowują stabilność dla dowolnej wielkości kroku czasowego – kosztem większego kosztu obliczeniowego na krok i możliwości użycia znacznie większych kroków całkowitych, a tym samym mniejszej liczby kroków.

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Fourier Heat Equation i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.