Pytania otagowane jako time-integration

2
Co oznacza „symplektyczny” w odniesieniu do integratorów numerycznych i czy używa ich odeint SciPy?
W tym komentarzu napisałem: ... domyślny integrator SciPy, który, jak zakładam, używa tylko metod symplektycznych. w którym mam na myśli SciPy odeint, który używa albo „metody niesztywnej (Adamsa)”, albo „metody sztywnej (BDF)”. Według źródła : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, …
2
Co to jest pseudo stepping?
Czytając literaturę o rozwiązaniach PDE, spotkałem się dziś z terminem pseudo-stepping . Wydaje się, że jest to powszechny termin, jednak nie znalazłem dla niego dobrej definicji ani artykułu wprowadzającego. Stąd: Co to jest pseudo stepping i jak się go zwykle stosuje?
1
BDF vs niejawne przyspieszenie czasu Runge Kutta
Czy są jakieś powody, dla których należy wybierać domyślny Runge Kutta (IMRK) wyższego rzędu zamiast skokowego czasu BDF? BDF wydaje mi się znacznie łatwiejszy, ponieważ etap IMRK potrzebuje qqqqqqq rozwiązań liniowych na krok czasowy. Stabilność BDF i IMRK wydaje się być kwestią sporną. Nie mogę znaleźć żadnych zasobów porównujących / …
1
Jaki jest właściwy sposób włączenia do symulacji astronomicznych?
Tworzę prosty symulator astronomiczny, który powinien wykorzystywać fizykę newtonowską do symulacji ruchu planet w układzie (lub dowolnych obiektach, jeśli o to chodzi). Wszystkie ciała są okręgami na płaszczyźnie euklidesowej, które mają takie właściwości, jak położenie, prędkość, masa, promień i wynikowa siła. Chcę aktualizować wszechświat w krótkich odstępach czasu, zwykle kilka …
5
Przykłady obliczeń PDE z wykorzystaniem równoległości w przestrzeni i czasie
W numerycznym rozwiązaniu początkowych wartości granicznych PDE bardzo często stosuje się równoległość w przestrzeni . O wiele rzadziej stosuje się jakąś formę paralelizmu w dyskretyzacji czasu , a paralelizm ten jest zwykle znacznie bardziej ograniczony. Jestem świadomy rosnącej liczby kodów i opublikowanych prac wykazujących równoległość czasową, ale żaden z nich …
1
Jak sformułować macierz bryły skupionej w MES
Przy rozwiązywaniu PDE zależnych od czasu za pomocą metody elementów skończonych, na przykład powiedzmy równanie cieplne, jeśli zastosujemy wyraźne stopniowanie czasu, wówczas musimy rozwiązać układ liniowy z powodu macierzy masy. Na przykład jeśli trzymamy się przykładu równania cieplnego, ∂u∂t= c ∇2)u∂u∂t=do∇2)u\frac{\partial{u}}{\partial{t}} = c\nabla{}^{2}u następnie za pomocą forward Euler otrzymujemy M.( …
2
Test integratora symetrycznego trzeciego rzędu vs czwartego rzędu z dziwnym wynikiem
W mojej odpowiedzi na pytanie dotyczące MSE dotyczące symulacji fizyki Hamiltoniana 2D zasugerowałem użycie integratora symplektycznego wyższego rzędu . Potem pomyślałem, że dobrym pomysłem może być wykazanie wpływu różnych kroków czasowych na globalną dokładność metod z różnymi zamówieniami. Napisałem i uruchomiłem skrypt Python / Pylab w tym celu. Do porównania …
1
Jakie są różnice między Parareal, PITA i PFASST?
Algorytmy Parareal, PITA i PFASST są technikami obejmującymi całą domenę , służącymi do równoległego rozwiązywania problemów zależnych od czasu w czasie. Jakie są główne zasady tych metod? Jakie są główne różnice między nimi? Czy mogę powiedzieć, że jedna opiera się na innej? W jaki sposób? Co z ich aplikacjami? Wiem, …
3
Czy mogę użyć jawnego schematu krokowego w celu ustalenia liczbowego, czy ODE jest sztywny?
Mam ODE: u′=−1000u+sin(t)u′=−1000u+sin(t)u'=-1000u+sin(t) u(0)=−11000001u(0)=−11000001u(0)=-\frac{1}{1000001} Wiem, że ten konkretny ODE jest sztywny, analitycznie. Wiem również, że jeśli użyjemy jawnej (krokowej) metody krokowej (Euler, Runge-Kutta, Adams itp.), Metoda powinna zwrócić bardzo duże błędy, jeśli krok czasowy jest zbyt duży. Mam więc dwa pytania: Czy w ten sposób określa się sztywne wartości ODE …
2
Dyskretyzacja czasoprzestrzennych elementów skończonych dla PDE zależnych od czasu
W literaturze MES metody półwariacyjne są zwykle stosowane w rozwiązaniu zależnych od czasu PDE. Nie widziałem podejścia w pełni wariacyjnego, tj. Gdzie MES dyskretuje przestrzeń i czas, być może pozwalając na użycie nieustrukturyzowanych siatek czasoprzestrzennych. Chociaż metody pomiaru czasu mogą być łatwiejsze do wdrożenia, czy istnieje szczególny powód, dla którego …
1
Algorytm obliczania wykładniczej macierzy Hessenberga
Interesuje mnie obliczenie rozwiązania systemu lage ODE przy użyciu metody krylova jak w [1]. Taka metoda obejmuje funkcje związane z wykładniczym (tzwφφ\varphi-Funkcje). Zasadniczo polega na obliczeniu działania funkcji macierzowej przez zbudowanie podprzestrzeni Kryłowa za pomocą iteracji Arnoldiego i rzutowanie funkcji na tę podprzestrzeń. Zmniejsza to problem obliczania wykładniczej znacznie mniejszej …
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.