PL 
EnglishDomyślny
Українська🌡️ Równanie przewodnictwa cieplnego Fouriera — dyfuzja i warunki brzegowe
Rozwiąż jednowymiarowe równanie przewodnictwa cieplnego Fouriera metodą różnic skończonych. Ustaw warunki brzegowe Dirichleta lub Neumanna i obserwuj dyfuzję ciepła w pręcie oraz wykładnicze zanikanie modów.
O tej symulacji
Ta symulacja numerycznie rozwiązuje jednowymiarowe równanie przewodnictwa cieplnego ∂T/∂t = α·∂²T/∂x² na pręcie o 200 punktach, wykorzystując jawny schemat różnic skończonych. Wybierasz początkowy profil temperatury i warunek brzegowy, a następnie obserwujesz, jak ciepło rozprasza się i wygładza w czasie rzeczywistym. Aktualizowany na żywo wykres słupkowy rozkłada krzywą temperatury na pierwsze pięć sinusoidalnych modów Fouriera, pokazując wprost, dlaczego cechy o wyższej częstotliwości zanikają najszybciej.
🔬 Co pokazuje
Pręt złożony ze 200 punktów siatki ewoluuje według wzoru T_i^(n+1) = T_i^n + r·(T_{i+1} − 2T_i + T_{i-1}), gdzie r = α·dt/dx², a dt jest dobierane automatycznie na poziomie 80% granicy stabilności (r ≤ 0,5). Pasek kolorów i wykres profilu śledzą tę samą tablicę temperatur; osobny wykres słupkowy pokazuje amplitudę sinusoidalnych modów Fouriera n = 1 do 5, wyodrębnionych z bieżącego profilu.
🎮 Jak korzystać
Ustaw suwak Thermal Diffusivity α (0,01–0,5), wybierz warunek początkowy (Gaussian Pulse, Step Function, Sine Wave, Hot End) oraz typ warunku brzegowego (Dirichlet — stała temperatura, Neumann — izolacja, lub Periodic), a następnie naciśnij Run. Cztery przyciski presetów wczytują gotowe kombinacje — spróbuj Hot End z warunkiem Dirichleta, aby zobaczyć front ciepła rozchodzący się od jednego końca, lub Periodic z falą sinusoidalną, aby zaobserwować czysty zanik pojedynczego modu.
💡 Czy wiesz, że
Wykres modów Fouriera to nie ozdoba — to ten sam rozkład sinusoidalny, który Joseph Fourier wprowadził w 1822 roku, aby rozwiązać dokładnie to równanie. Ponieważ każdy mod n zanika jak e^(−(nπ/L)²·α·t), wyższe mody znikają w przybliżeniu n² razy szybciej niż mod podstawowy, dlatego każdy nieregularny profil początkowy zawsze wygładza się ostatecznie do jednego gładkiego garbu.
Najczęściej zadawane pytania
Jakie równanie faktycznie rozwiązuje ta symulacja?
Rozwiązuje jednowymiarowe równanie przewodnictwa cieplnego (dyfuzji) ∂T/∂t = α·∂²T/∂x² na pręcie o 200 punktach siatki, wykorzystując jawny schemat różnic skończonych: każda nowa wartość temperatury to stara wartość powiększona o poprawkę proporcjonalną do krzywizny sąsiadów, przeskalowaną przez r = α·dt/dx². Krok czasowy dt jest przeliczany na podstawie dyfuzyjności α i ustalany automatycznie na poziomie 80% granicy stabilności (r ≤ 0,5), dzięki czemu symulacja nigdy się nie „wysypuje”, niezależnie od ustawienia suwaka dyfuzyjności.
Co fizycznie oznaczają trzy warunki brzegowe?
Dirichlet utrzymuje oba końce pręta w stałej temperaturze — 1,0 na gorącym końcu w presecie „Hot End”, 0,0 w pozostałych przypadkach — jak pręt z końcami przyciśniętymi do grzejnika i kąpieli lodowej. Neumann ustala zerowy gradient temperatury na każdym końcu (T_new[0]=T_new[1], T_new[N-1]=T_new[N-2]), modelując pręt doskonale izolowany, który nie może wymieniać ciepła z otoczeniem. Periodic zwija pręt w pierścień, tak że ciepło, które opuściłoby prawy koniec, ponownie wchodzi z lewej strony.
Jak suwak Thermal Diffusivity zmienia zachowanie symulacji?
Suwak α (zakres od 0,01 do 0,5) ustala, jak szybko rozprzestrzenia się ciepło: pojawia się bezpośrednio w aktualizacji metodą różnic skończonych oraz w automatycznie obliczanym kroku czasowym dt = 0,8 × 0,4 × dx² / α. Większe α wygładza profil temperatury szybciej w czasie symulowanym, a ponieważ dt skaluje się odwrotnie proporcjonalnie do α, przy wysokiej dyfuzyjności symulacja przechodzi też przez więcej czasu fizycznego na każdą renderowaną klatkę.
Co reprezentują cztery warunki początkowe?
Gaussian Pulse zaczyna od dzwonowego gorącego punktu wyśrodkowanego na środku pręta; Step Function zaczyna jako płaski gorący obszar między x=0,3 a x=0,7; Sine Wave zaczyna jako pojedyncza półokresowa krzywa sinusoidalna, będąca czystym modem Fouriera, dzięki czemu jest najprostszym przypadkiem do obserwowania zaniku; Hot End ustala maksymalną temperaturę tylko dla lewych 5% pręta, symulując jeden koniec przyciśnięty do źródła ciepła.
Dlaczego wykres modów Fouriera ma znaczenie i jak jest obliczany?
Dla każdego modu n od 1 do 5 symulacja numerycznie całkuje bieżącą tablicę temperatur względem sin(nπx/L), aby wyodrębnić amplitudę tego modu — dokładnie dyskretny odpowiednik współczynnika sinusoidalnego szeregu Fouriera. Ma to znaczenie, ponieważ dokładne analityczne rozwiązanie równania przewodnictwa cieplnego to T(x,t) = Σ aₙ·sin(nπx/L)·e^(−(nπ/L)²·α·t): każdy mod zanika wykładniczo, przy czym wyższe n zanika znacznie szybciej, co jest właśnie tym zachowaniem wygładzającym widocznym na wykresach pręta i profilu.