• Nie Znaleziono Wyników

(Równanie dyfuzji) Zapisz jawny schemat Eulera dla równania dyfuzji, a następnie korzystając z zasady maksimum określ warunki dla których schemat będzie stabilny

N/A
N/A
Protected

Academic year: 2021

Share "(Równanie dyfuzji) Zapisz jawny schemat Eulera dla równania dyfuzji, a następnie korzystając z zasady maksimum określ warunki dla których schemat będzie stabilny"

Copied!
2
0
0

Pełen tekst

(1)

Przykładowe pytania egzaminacyjne.

1. (RRZ) Zapisz schemat różnicowy wykorzystując niejawny schemat trapezów dla poniższego rów- nania różniczkowego

d2u(t)

dt2 = s(t)− g(t) · u(t) (1)

gdzie: s(t) i g(t) są znanymi funkcjami.

2. (RRZ) Dla równania oscylatora harmonicznego z tlumieniem d2x

dt2 =−ω2x− βdx

dt (2)

zapisz schemat różnicowy korzystając z jawnej metody punktu pośredniego (metoda RK) okre- ślonej tablicą Butchera

0 0 0

1 2

1

2 0

0 1

3. (Równanie adwekcji) Dla jednowymiarowego równania adwekcji ∂u/∂t = −v∂u/∂x zaproponuj jawny i niejawny schemat Eulera pozwalający rozwiązać je numerycznie. Z kombinacji obu me- tod (jawny i niejawny Euler) utwórz schemat Crancka-Nicolson i określ, stosując analizę von Neumanna, obszar jego bezwzględnej stabilności.

4. (Równanie dyfuzji) Zapisz jawny schemat Eulera dla równania dyfuzji, a następnie korzystając z zasady maksimum określ warunki dla których schemat będzie stabilny.

5. (Równanie Poissona) Dane jest jednowymiarowe równanie Poissona z warunkami brzegowymi d2u

dx2 =−ρ(x), x ∈ [0, L], u(0) = uA, du dx

x=L

= C (3)

Dokonaj jego dyskretyzacji stosując odpowiedni iloraz różnicowy na siatce zbudowanej z 4 rów- noodległych węzłów, następnie zapisz równanie w postaci macierzowej A⃗u = ⃗b (z uwzględnieniem warunków brzegowych). Podaj postać A i ⃗b.

6. (RRZ) Dla poniższego równania różniczkowego dy(t)

dt + A· y2(t) + g(t) = 0 (4)

gdzie: A jest stałą, a g(t) jest pewną znaną funkcją czasu, proszę zapisać schemat różnicowy (szukamy y(t) dla chwili czasowej tn+1) wykorzystując metodę RK2 określoną przez tablicę Butchera

0 0 0 1 1 0

1 2

1 2

7. (RRZ) Dla problemu autonomicznego du/dt = λ· u proszę wyznaczyć współczynnik wzmocnie- nia i określić region bezwzględnej stabilności niejawnej metody RK punktu pośredniego, którą opisuje tablica Butchera

1 2

1 2

1

8. (Równanie adwekcji) Zapisz schemat różnicowy Laxa-Friedrichsa dla równania adwekcji a na- stępnie, stosując analizę von Neumanna, określ współczynnik wzmocnienia tego schematu. Jaki warunek musi być spełniony, aby metoda była bezwzględnie stabilna?

(2)

9. (Dyfuzja numeryczna) Wyznacz współczynniki dyfuzji numerycznej w schematach Laxa-Friedrichsa i Laxa-Wendroffa dla równania adwekcji.

10. (Równanie falowe) Stosując rozwinięcie w szereg Taylora skonstruuj schematy: położeniowy i prędkościowy Verleta dla równania falowego. Określ rząd dokładności obu schematów.

Cytaty

Powiązane dokumenty

Wyznacz wartości parametru m, dla których równanie ma dwa różne pierwiastki o różnych znakach.. Wyznacz wartości parametru m, dla których równanie ma dwa różne pierwiastki

1. Przeanalizuj graf obrazujący ciąg reakcji, jakim ulegają związki organiczne i wykonaj polecenia znajdujące się pod grafem. Używaj jedynie wzorów

Tego typu mapki nieba, można sobie wygenerować za pomocą internetowych programów astronomicznych (np. Heavens-Above) – wystarczy podać współrzędne geograficzne

→ jeśli rozwiązanie startowe jest „bliskie” dokładnemu to ilość iteracji może być mała (rel. Poissona nie trzeba jej nawet tworzyć (zysk w postaci ograniczenia

(Równanie falowe) Skonstruuj niejawny schemat różnicowy Eulera dla równania falowego i określ rząd dokładności poszczególnych rozwiązań względem ∆t i

(Równanie adwekcji) Dla równania adwekcji zastosuj schemat numeryczny z centralną pochodną przestrzenną, a następnie użyj analizy von Neumanna do określenia

→ jeśli M jest macierzą rzadką to koszt jednej iteracji jest rzędu O(n), dla pełnej macierzy O(n 2 ). → jeśli rozwiązanie startowe jest „bliskie” dokładnemu to

Rysunek 1: Geometria układu, w którym rozkład gęstości u(x, y, t) zmienia się ze względu na adwekcję w polu prędkości i pod wpływem dyfuzji... , 20 jest numerem