Lista zajętych projektów:
2 Wahadło z siłą wymuszającą
Wykonaj zadanie 12.8 i 12.10 z podanej
literatury, ale dla różnych amplitud siły wymuszającej
(parametru kontrolnego) tak, by zaobserwować zachowanie chaotyczne i
okresowe. Porównaj wyniki uzyskane dwoma (dowolnymi) metodami
całkowania różnych rzędów. Wyniki opisz w raporcie
końcowym.
Literatura:
- Mechanika Klasyczna 2: rozdział 12. Mechanika nieliniowa i chaos, John R. Taylor
3 Symulacja periodycznej trajektorii 3 ciał
Przeprowadzenie symulacji układu 3 ciał
oddziaływujących grawitacyjnie, z warunkami początkowymi podanymi z
literaturze do zadania.
Zbadanie energii i momentu pędu układu w funkcji czasu. Porównanie metod całkowania : Explicit
Euler, Runge-Kutta i Verlet.
Wyniki opisz w raporcie końcowym.
4 Symulacja periodycznej trajektorii 4 ciał
Przeprowadzenie
symulacji układu 4 ciał oddziaływujących grawitacyjnie,
z warunkami początkowymi podanymi z literaturze do zadania. Zbadanie
energii i momentu pędu układu w funkcji czasu. Porównanie
metod całkowania : Explicit Euler, Runge-Kutta i Verlet. Wykonanie
animacji.
Wyniki opisz w raporcie końcowym.
Literatura:
5 Symulacja periodycznej trajektorii 5 ciał
Przeprowadzenie
symulacji układu 5 ciał oddziaływujących grawitacyjnie,
z warunkami początkowymi podanymi z literaturze do zadania. Zbadanie
energii i momentu pędu układu w funkcji czasu. Porównanie
metod całkowania : Explicit Euler, Runge-Kutta i Verlet. Wykonanie
animacji.
Wyniki opisz w raporcie końcowym.
Literatura:
6 Symulacja układu Lorenza
Analiza układu Lorenza danego układem równań
Rozważ
kilka zestawów parametrów (w tym używany przez
Lorenza) i opisz zachowanie układu. Wykonaj animację
atraktora (motyla). Znajdz i opisz punkty krytyczne (obiegane przez
atraktor). Porównaj wyniki uzyskane dwoma (dowolnymi)
metodami całkowania różnych rzędów. Wyniki opisz w
raporcie końcowym.
Literatura:
11 Temperatura zwilżania w 2D modelu Isinga z polem zewnętrznym (pomysł: Piotr Nowakowski)
Układ dwuwymiarowy, sieć kwadratowa, rozmiar M na N. W kierunku M okresowe (periodyczne) warunki brzegowe, M>>N Stała sprzężenia
między spinami J. W kierunku N na brzegach spiny są swobodne i działa na nie zewnętrzne pole brzegowe. Na spiny z jednej strony pole +H1, na
spiny z drugiej strony -H1. W takim układzie są dwie charakterystyczne temperatury - przemiana objętościowa TC, dla J/kT=1/2 Log(1+Sqrt{2})
oraz temperatura zwilżania TW zależna od H1.
Znajdz temperaturę zwilżania oraz profil magnetyzacji dla T<TW algorytmem Metropolisa.
Literatura