Rozkład Temperatury: Jak Ciepło Rozprzestrzenia Się W Ciele Stałym

Równanie opisujące rozprzestrzenianie się ciepła

Ciepło przemieszcza się w materiale stałym przez przewodzenie cieplne: energia przeskakuje od gorących, bardziej energetycznie drgujących atomów do chłodniejszych sąsiadujących bez jakiejkolwiek ruchomości materiału. Prawo Fouriera dotyczące przewodzenia ciepła (1822) stanowi prawo rządzące tym procesem, mówiące, że lokalny przepływ ciepła jest proporcjonalny do lokalnego gradientu temperatury – ciepło płynie z gorącego do zimnego i płynie szybciej tam, gdzie zmiana temperatury jest bardziej stroma:

q = -k ċ ∇T (Prawo Fouriera: strumień ciepła w funkcji gradientu temperatury) ∂T/∂t = α ċ ∇²T (Równanie Fouriera, łączące prawo Fouriera z zasadą zachowania energii) k = przewodność cieplna (W/m·K) α = k / (ρc_p) = złożona przenikalność cieplna (m²/s) demo na żywo · pole temperatury dążące do stanu równowagi w długotrwałym pręcie● NA ŻYWO Połączenie prawa Fouriera z zasadą zachowania energii daje równanie Fouriera, równanie różniczkowe cząstkowe opisujące ewolucję pola temperatury T(x,t) w przestrzeni i czasie. Jest ono matematycznie identyczne w formie do równania dyfuzji, które rządzi rozprzestrzenianiem się barwnika przez nieruchomą wodę lub rozpraszaniem gazu w pokoju – ciepło, w tym sensie, dosłownie się dyfuzyjnie rozprasza.

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)

Stan staczy się w porównaniu z stanem przejściowym: dwa zupełnie różne problemy

Gdy końce prutu są utrzymywane w stałych temperaturach przez wystarczająco długi czas, pole temperatury przestaje zmieniać się w czasie – zanika wyraz ∂T/∂t – pozostawiając znacznie prostsze równanie Laplace’a: ∇²T = 0. Dla prostego pruta 1D bez wewnętrznych źródeł ciepła to zmusza temperaturę stanu staczy do zmiany liniowej między dwoma końcowymi temperaturami, niezależnie od przewodności materiału – przewodność wpływa na to, ile ciepła przepływa i jak długo trwa osiągnięcie równowagi, ale nie wpływa na ostateczny kształt stanu staczy. Problem przejściowy – czyli jak pole temperatury ewoluuje w drodze do tej równowagi, zaczynając od pewnej początkowej rozkładu – jest trudniejszą, zależną od czasu wersją i to tam faktycznie ma znaczenie pełne równanie ciepła z wyrazem ∂T/∂t.

Analiza rozwiązań i ich ograniczenia

Dla prostych geometrii – nieskończony pręt, płyta o ustalonych temperaturach brzegowych, kula chłodząca w jednorodnym środowisku – równanie ciepła ma dokładne, zamknięte rozwiązania analityczne, często opierające się na szeregach Fouriera (ironicznym pełnym cyklu, ponieważ Józef Fourier opracował swój nazwany jego imieniem szereg specjalnie do rozwiązania tego równania) lub z funkcji błędu dla półnieskończonych domen. Te rozwiązania są bezcennym narzędziem weryfikacji, ale działają tylko dla idealnych kształtów o prostych, niezmiennych warunkach brzegowych. Elementy rzeczywiste – rzadzeń o zmiennej przekrojach poprzecznych, pręty oddające ciepło konwekcyjnie do otaczającego powietrza, obiekty z rozkładem temperatury w przestrzeni – wymagają metod numerycznych.

Metoda różnic skończonych: przekształcanie rachunku różniczkowego w arytmetykę

Metoda różnic skończonych (MRŚ) zastępuje ciągłą pręt z dyskretową łańcuchą punktów oddzielonych odstępem Δx, a terminu drugiego pochodnego we wzorze na przewodzenie ciepłem zamienia się na prosty trzyprzedni uśredniony:

∂²T/∂x² ≈ ( T[i+1] - 2·T[i] + T[i-1] ) / Δx² aktualizacja jawna: T_nowe[i] = T[i] + α · Δt / Δx² · ( T[i+1] - 2·T[i] + T[i-1]) stabilność wymaga: α · Δt / Δx² ≤ 0.5 (MRŚ w 1D) To ograniczenie stabilności nie jest sugestią – powyższy schemat jawny rzeczywiście ulega oscylacjom i chaotycznym rozwiązaniom, jeśli krok czasowy jest zbyt duży w stosunku do odstępu przestrzennego, stanowiąc formę tego samego warunku CFL ograniczającego kroki czasowe w symulacjach płynów i fal. Podwojenie odstępu przestrzennego podwaja dokładność, a zatem zmniejsza maksymalny stabilny krok czasowy o czwarty, ponieważ Δx pojawia się do kwadratu w tym ograniczeniu, dlatego też precyzyjne symulacje transjentne przewodzenia ciepła mogą szybko generować wysokie koszty. Schematy ukryte (implicitne) oferują bardziej kosztowne rozwiązanie liniowe na każdym kroku w zamian za bezwarunkową stabilność, pozwalając inżynierom przyjmować znacznie większe kroki czasowe, gdy tylko interesuje nas ostateczne zachowanie stałe, a drobna rozdzielczość czasowa transjentnego procesu nie jest istotna.

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

Płytki chłodne: dlaczego radiatory mają taki kształt

Płytka chłodząca traci ciepło zarówno przez przewodzenie wzdłuż własnej długości, jak i przez konwekcję, czyli przenoszenie tego ciepła na zewnątrz na powierzchnię do otaczającego powietrza. To sprawia, że równanie opisujące działanie płytki chłodnej jest bardziej złożone niż równanie dla prostej belki – uwzględnia ono balans między przewodzeniem wzdłuż płytki a stratą ciepła przez konwekcję ze ścianek.

Rozwiązanie tego równania pokazuje, że temperatura spada w przybliżeniu wykładniczo od podstawy płytki. Dlatego też radiatory mają cienkie i liczne płytki zamiast kilku grubych – większa powierzchnia na jednostkę objętości przenosi więcej ciepła do powietrza, nawet jeśli cieńsza płytka przewodzi mniej ciepła wzdłuż swojej długości.

Frequently asked questions

Wypróbuj na żywo

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