Gdzie mogę znaleźć dobre odniesienie do właściwości stabilności kilku metod rozwiązywania parabolicznych PDE?

10

W tej chwili mam kod, który używa algorytmu Crank-Nicholson, ale myślę, że chciałbym przejść do algorytmu wyższego rzędu w celu pomiaru czasu. Wiem, że algorytm Crank-Nicholson jest stabilny w dziedzinie, w której chcę pracować, ale martwię się, że niektóre inne algorytmy mogą nie być.

Wiem, jak obliczyć region stabilności algorytmu, ale może to być rodzaj bólu. Czy ktoś zna jakieś dobre odniesienie do właściwości stabilności wielu algorytmów pomiaru czasu dla parabolicznych PDE?

Dan
źródło

Odpowiedzi:

10

Bardzo krótka odpowiedź: dla wyczerpującego opisu nie można pokonać Hairer i Wanner tom II .

Krótka odpowiedź: Oto kilka skryptów MATLAB do wykreślenia regionu stabilności liniowej metody wieloetapowej lub metody Runge-Kutta , biorąc pod uwagę współczynniki. Możesz także użyć pakietu nodepy pakietu Python (zastrzeżenie: to mój pakiet i nie jest to najbardziej dopracowane oprogramowanie, ale wykreślanie regionów stabilności to jedna rzecz, którą robi bardzo dobrze). Instrukcje dotyczące kreślenia regionów stabilności znajdują się tutaj .

Dłuższa odpowiedź: są trzy klasy metod, którymi możesz się tutaj zainteresować.

  • AA-stabilność. Niektóre przykłady takich metod to metody Gaussa-Legendre'a, Radau i Lobatto. Wszystkie są w pełni dorozumiane, a zatem dość drogie.

  • A(α)ode15s()α

  • Jawne metody , które będą konieczne, obejmują tylko skończony przedział na ujemnej osi rzeczywistej. Istnieją specjalne „ustabilizowane” wyraźne metody (w szczególności metody Runge-Kutta-Chebyshev ), które mają duże ujemne regiony stabilności osi rzeczywistej i są odpowiednie dla lekko sztywnych problemów, ale zwykle nie dla problemów parabolicznych. Dobrym wejściem do tej literatury jest ten artykuł , który zawiera wiele informacji o regionach stabilności.

LL

Aktualizacja : Jeśli naprawdę potrzebujesz wiedzieć wszystko na ten temat, zdobądź kopię monografii Dekkera i Verwera . Ma jedno z najlepszych istniejących pojęć takich jak jednostronne stałe Lipschitza, norma logarytmiczna i kilka głębszych koncepcji stabilności. Jest wyczerpany, ale zwykle możesz znaleźć używane kopie na Amazon (za cenę!)

David Ketcheson
źródło
Hairer II jest zdecydowanie najlepszy. Jest to prawdopodobnie jedyne miejsce, w którym można znaleźć wzmiankę na temat adaptacji PI w krokach. Brakuje jednak ważnych szczegółów, takich jak dodatkowe warunki zamówienia dla metod Rosenbrocka na parabolicznych PDE. Oczywiście żadna książka nie może mieć wszystkiego, ale powinno być coś lepszego konkretnie na temat parabolicznych PDE.
Chris Rackauckas,