Co analiza statyczności von Neumanna mówi nam o nieliniowych równaniach różnic skończonych?


9

Czytam artykuł [1], w którym rozwiązują następujące równanie nieliniowe przy użyciu metod różnic skończonych. Analizują także stabilność schematów za pomocą analizy stabilności von Neumanna. Jednak, jak zdają sobie sprawę autorzy, dotyczy to tylko liniowych PDE. Autorzy więc to poprzez „zamrożenie” nieliniowego terminu, tzn. termin , gdzie „uważa się za reprezentujące lokalnie stałe wartości ”.

ut+ux+uuxuxxt=0
uuxUuxUu

Więc moje pytanie jest dwojakie:

1: jak interpretować tę metodę i dlaczego (nie) działa?

2: czy moglibyśmy również zastąpić termin terminem , gdzie jest „uważany za reprezentujący lokalnie stałe wartości ”?uuxuUxUxux

Bibliografia

  1. Eilbeck, JC i GR McGuire. „Studium numeryczne znormalizowanego równania długofalowego I: metody numeryczne”. Journal of Computational Physics 19.1 (1975): 43-57.

1
Źle wpisałeś równanie. Równanie w artykule jest równaniem RLW.
Ömer

3
Powiązane pytania, bez pełnych odpowiedzi: scicomp.stackexchange.com/q/8717/713 , mathoverflow.net/q/186760 , scicomp.stackexchange.com/q/16142 , scicomp.stackexchange.com/q/6863 . Myślę, że heurystycznie rzecz biorąc, powinno to działać, ponieważ interesuje Cię stabilność trybów o bardzo wysokiej częstotliwości (przy których występują błędy, długość fali w kolejności odstępów między siatkami), podczas gdy samo rozwiązanie zmieniałoby się natomiast ze znacznie niższą częstotliwością, więc można zamrażać współczynniki i badać stabilność zamrożonych współczynników PDE.
Kirill,

2
Udzieliłem odpowiedzi na niektóre pytania powiązane z Kirillem. Niestety nie znam żadnych wyników dla równania RLW, ale prawdopodobnie stabilność można udowodnić, o ile rozwiązanie jest wystarczająco gładkie.
David Ketcheson

Odpowiedzi:


1

To, co mówisz, nazywa się linearyzacją. Jest to powszechna technika stosowana w analizie nieliniowych PDE. Wykonane jest rzutowanie równań w formacie,

ut+Au=0

Tutaj A jest macierzą wynikającą z linearyzacji równania.

Teraz do twoich pytań,

  1. Jak myślisz, działa do pewnego stopnia, ale do pewnego stopnia nie. Użyteczność polega na tym, że stabilność można udowodnić w przypadku systemów liniowych, ale niełatwo w przypadku systemów nieliniowych. Zatem wyniki liniowe są rozszerzone na układy nieliniowe. Często w poszczególnych przypadkach przyjmuje się różne metody. Na przykład,

uux=12(u2)x

która jest formą ochrony. Więc,

ut+12(u2)x=0

gdy jest reprezentowany w sensie skończonej objętości, ogranicza ograniczenia ewolucji u.

  1. Jaka jest użyteczność wymiany. Usuniesz równanie z postaci równania falowego. Co oznaczałoby, że rozwiązania nie zachowywałyby się jak równanie falowe. Zatem w analizie stabilności rozwiązania testowe musiałyby być zupełnie inne i niefizyczne.

2

Aby rozwinąć argument dotyczący linearyzacji, w uu_x chcesz założyć, że u jest lokalnie stały, a nie u_x, z dwóch powodów: a) u zmienia się wolniej niż jego pochodna, i b) w tym konkretnym przypadku, jeśli założymy, że u_x jest lokalnie stały , z definicji zakładasz także, że u Ciebie jest lokalnie liniowy, co oznacza, że ​​pochodne wyższych przestrzeni są zerowe, co nie tylko wprowadza dodatkowy błąd aproksymacji, ale może oznaczać, że możesz wylewać dziecko z kąpielą, w zależności od twojego równania.

Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.