PL 
EnglishDomyślny
Українська🌊 Fale chaotyczne — łańcuchy oscylatorów nieliniowych
Symuluj fale chaotyczne w sprzężonych oscylatorach nieliniowych. Poznaj problem Fermiego-Pasty-Ulama-Tsingou, podział energii między mody i chaos falowy w modelach Duffinga i sieci Tody.
O łańcuchu Fermiego–Pasty–Ulama–Tsingou
Ten symulator modeluje łańcuch N punktowych mas połączonych identycznymi nieliniowymi sprężynami, z obydwoma końcami zamocowanymi do nieruchomych ścian. Każda masa i spełnia mẍⁱ = V′(xᵢ₊₁−xᵢ) − V′(xᵢ−xᵢ₋₁), całkowane symplektycznym schematem velocity Verlet, dzięki czemu całkowita energia jest dobrze zachowana przez długie przebiegi. Do wyboru są trzy potencjały sprężyste: α-FPUT (V = ½x² + ⅓αx³, mała poprawka sześcienna do prawa Hooke'a), β-FPUT (V = ½x² + ¼βx⁴, poprawka czwartego rzędu) oraz sieć Tody (V = (eᵃˣ − αx − 1)/α², dokładnie całkowalna sprężyna wykładnicza). Ustawienie suwaka nieliniowości na 0 przywraca czysto harmoniczny łańcuch, w którym energia wprowadzona do jednego modu normalnego pozostaje w nim na zawsze — mody normalne łańcucha liniowego nigdy nie wymieniają energii.
Przy starcie (i po Reset) łańcuch jest wzbudzany w pojedynczym modzie normalnym k, xᵢ(0) = A·sin(kπi/(N+1)), za pomocą suwaków modu wzbudzenia, amplitudy początkowej i długości łańcucha N. Górny canvas rysuje chwilowe wychylenie każdej masy jako łańcuch sprężyn i kulek, z kolorem kulki kodującym lokalną energię kinetyczną (niebieski = zimny, czerwony = gorący). Dolny canvas to widmo energii modów: w każdej klatce bieżący stan jest rzutowany na N dyskretnych sinusowych modów normalnych łańcucha liniowego (dyskretna transformata sinusowa), a energia harmoniczna każdego modu Eₖ = ½pₖ² + ½ωₖ²qₖ² jest gromadzona w czasie, pokazując, jak — i czy w ogóle — energia wycieka z początkowo wzbudzonego modu k do jego sąsiadów. Panel Statystyk śledzi upływający czas symulacji, całkowitą (zachowaną) energię, energię pozostającą w modzie k oraz największe wychylenie w łańcuchu.
Najczęściej zadawane pytania
Czym jest problem FPUT?
W 1953 roku Fermi, Pasta, Ulam i (nieuznana aż do znacznie później) Tsingou przeprowadzili jeden z pierwszych eksperymentów komputerowych, symulując łańcuch 64 mas z małą nieliniowością sześcienną (α) dodaną do w innym wypadku liniowych sprężyn. Spodziewali się, że początkowe wzbudzenie pojedynczego modu szybko rozłoży swoją energię równomiernie na wszystkie mody normalne (termalizacja/ekwipartycja). Zamiast tego zaobserwowali niemal-rekurencję: energia dzielona krótko z kilkoma sąsiednimi modami, po czym wracająca niemal w całości do pierwotnego modu, wielokrotnie, bez oznak termalizacji. Ten symulator odtwarza dokładnie ten sam układ — wzbudź jeden mod, obserwuj widmo energii modów i zobacz, czy energia wraca.
Jaka jest różnica między modelami α-FPUT, β-FPUT i Toda?
Wszystkie trzy zastępują idealną sprężynę liniową V=½x² lekko nieliniową. α-FPUT dodaje wyraz sześcienny (asymetryczny, miększy w jednym kierunku i sztywniejszy w drugim); β-FPUT dodaje wyraz czwartego rzędu (symetryczny, usztywniający w obu kierunkach); sieć Tody wykorzystuje potencjał wykładniczy, który jest dokładnie rozwiązywalny (całkowalny) metodą odwrotnego rozpraszania i wspiera dokładne rozwiązania solitonowe. Zwiększanie suwaka nieliniowości (α lub β) wzmacnia sprzężenie między modami normalnymi; model Tody rekurencjonuje niemal idealnie przy dowolnej amplitudzie, podczas gdy α- i β-FPUT ostatecznie przechodzą w prawdziwy chaos i termalizację, gdy nieliniowość jest wystarczająco silna.
Dlaczego energia wraca do pierwotnego modu zamiast się rozprzestrzeniać?
Rozwiązanie, znalezione przez Zabusky'ego i Kruskala w 1965 roku, polega na tym, że łańcuch FPUT w granicy continuum zachowuje się jak równanie Kortewega–de Vriesa (KdV), którego rozwiązania są zdominowane przez solitony — samotne fale przechodzące przez siebie bez utraty tożsamości. Niemal-rekurencja jest dyskretnym echem dynamiki solitonowej: układ jest wystarczająco bliski układowi całkowalnemu (Toda, KdV), że przeszkoda twierdzenia KAM dla termalizacji nadal w dużej mierze obowiązuje. Dopiero przy silnej nieliniowości, gdy wystarczająco dużo sąsiednich modów nakłada się i rezonuje, łańcuch traci tę niemal-całkowalną strukturę i naprawdę się termalizuje w chaotyczny, ekwipartycjowany ruch.
Co dokładnie pokazuje panel widma energii modów?
W każdej klatce animacji symulator rzutuje bieżące położenia i prędkości na N modów normalnych łańcucha liniowego za pomocą dyskretnej transformaty sinusowej (qₖ = Σᵢ xᵢ·sin(kπi/(N+1)), i podobnie dla pędów), a następnie oblicza energię harmoniczną każdego modu Eₖ = ½pₖ² + ½ωₖ²qₖ² z relacją dyspersji ωₖ = 2sin(kπ/2(N+1)). Te energie poszczególnych modów są gromadzone w bieżącą sumę dla każdego z pierwszych 24 modów i pokazywane jako słupki, z podświetlonym początkowo wzbudzonym modem k — jest to dokładnie diagnostyka, której Fermi, Pasta, Ulam i Tsingou użyli w swoim oryginalnym raporcie z 1955 roku.
Dlaczego stosuje się całkowanie velocity Verlet zamiast prostszej metody?
Velocity Verlet to integrator symplektyczny: w przybliżeniu zachowuje strukturę hamiltonowską równań ruchu, więc całkowita energia oscyluje wokół stałej wartości zamiast dryfować przez długi czas symulacji. Dla układu chaotycznego lub niemal-całkowalnego, jak łańcuch FPUT, metoda niesymplektyczna (np. zwykły Euler) sztucznie traciłaby lub zyskiwała energię i mogłaby wytworzyć fałszywą termalizację lub fałszywą rekurencję. Statystyka Energii całkowitej w tym symulatorze powinna pozostawać zasadniczo płaska w czasie — to bezpośredni sprawdzian jakości integratora.
Co się stanie, jeśli ustawię suwak nieliniowości na zero?
Przy α=β=0 każdy model redukuje się do czysto harmonicznego łańcucha, którego dokładnym rozwiązaniem jest suma niezależnych modów normalnych, które nigdy nie wymieniają energii. Wzbudzenie modu k pozostawia widmo energii modów pokazujące całą energię trwale zaparkowaną w słupku k, przy każdym innym słupku równym zero — nie ma dudnień, nie ma rekurencji i nie ma chaosu, ponieważ liniowe mody normalne po prostu nie oddziałują. To jest punkt odniesienia, wobec którego zaprojektowano eksperyment FPUT.
Jaką rolę odgrywa tu twierdzenie KAM?
Twierdzenie Kołmogorowa–Arnolda–Mosera (KAM) opisuje, jak przy dodaniu małego nieliniowego zaburzenia do całkowalnego układu hamiltonowskiego (jak łańcuch liniowy czy całkowalna sieć Tody) większość quasi-okresowych niezmienniczych tori przetrwa, zamiast natychmiast rozpuścić się w chaos. Dlatego łańcuchy α- i β-FPUT wykazują długotrwałą niemal-rekurencję przy słabej do umiarkowanej nieliniowości zamiast natychmiastowej termalizacji: układ pozostaje bliski układowi całkowalnemu. Dopiero gdy suwak nieliniowości zostanie podniesiony wystarczająco wysoko, wystarczająco dużo rezonansów nakłada się, aby zniszczyć ocalałe tori i wytworzyć prawdziwe chaotyczne dzielenie energii między modami.
Jak długość łańcucha N i mod wzbudzenia k wpływają na dynamikę?
Większe N daje drobniejszą rozdzielczość przestrzenną i widmową oraz pozwala istnieć wyższym numerom modów, ale też zwiększa liczbę niemal-rezonansowych trójek/czwórek modów, które wyrazy nieliniowe mogą sprzęgać, na ogół ułatwiając termalizację przy danej nieliniowości. Wzbudzenie niskiego modu (małe k, długa fala) zwykle pokazuje najczystszą rekurencję FPUT, ponieważ jego najbliżsi sąsiedzi w widmie modów są dobrze rozdzieleni częstotliwościowo; wzbudzenie wysokiego modu (k blisko N) stawia układ bliżej krawędzi strefy Brillouina, gdzie relacja dyspersji spłaszcza się i sąsiednie mody łatwiej wchodzą w rezonans.
Dlaczego kulki na górnym canvasie zmieniają kolor?
Kolor każdej kulki koduje jej chwilową energię kinetyczną, ½mvᵢ², skalowaną względem najgorętszej kulki w łańcuchu w danym momencie (niebieski → niska energia kinetyczna, czerwony → wysoka energia kinetyczna). Obserwowanie wzorca kolorów jest szybkim wizualnym odpowiednikiem tego samego pytania o dzielenie energii, na które liczbowo odpowiada widmo modów: w reżimie niemal-rekurencyjnym „gorący punkt” przelewa się tam i z powrotem w dość zorganizowanym wzorcu związanym z długością fali wzbudzonego modu, podczas gdy w reżimie chaotycznym/stermalizowanym gorące i zimne kulki stają się przestrzennie nieuporządkowane i w przybliżeniu jednorodne w średniej czasowej.
Czy to jest ten sam układ, który pierwotnie symulowali Fermi, Pasta, Ulam i Tsingou?
Tak, strukturalnie: jednowymiarowy łańcuch mas na identycznych nieliniowych sprężynach z nieruchomymi (zamocowanymi) końcami, wystartowany z pojedynczego sinusoidalnego wzbudzenia modu normalnego i ewoluujący w czasie na komputerze — dokładnie taki jak układ z raportu Los Alamos z 1955 roku (pierwotnie z N=32 lub 64 masami i małym wyrazem α lub β). Ten symulator pozwala dodatkowo wybrać dokładnie całkowalny potencjał Todego, który wprowadzono później (1967), specjalnie po to, aby wyjaśnić — poprzez teorię solitonów — dlaczego oryginalny łańcuch FPUT odmawiał termalizacji.
Symuluj chaotyczne fale w sprzężonych nieliniowych oscylatorach. Odkrywaj problem Fermiego-Pasty-Ulama-Tsingou, dzielenie energii modów i chaos falowy poprzez modele Duffinga i sieci Todego.
3D · Renderer Three.js / WebGL · cel 60 FPS · działa w całości po stronie klienta, bez instalacji