Częściowe równania różniczkowe (PDE) to równania, które odnoszą się do pochodnych cząstkowych funkcji więcej niż jednej zmiennej. Ten znacznik jest przeznaczony do pytań o modelowanie zjawisk za pomocą PDE, rozwiązywanie PDE i innych powiązanych aspektów.
Silna postać PDE wymaga, aby nieznane rozwiązanie należało do . Ale słaba forma wymaga tylko, aby nieznane rozwiązanie należało do .H 1H.2)H.2)H^2H.1H.1H^1 Jak to pogodzisz?
Załóżmy, że miałem problem z wartością graniczną: dud2udx2+dvdx=f in Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega u=h in ∂Ωdudx+d2vdx2=g in Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u=h in ∂Ωu=h in ∂Ωu=h \text{ in } \partial\Omega Moim celem jest rozłożenie rozwiązania tego sprzężonego problemu na sekwencję …
Próbuję symulować podstawowe modele półprzewodnikowe do celów pedagogicznych - zaczynając od modelu dyfuzyjnego Drifta. Chociaż nie chcę używać gotowego symulatora półprzewodników - nauczę się innych (powszechnych, najnowszych lub niejasnych) modeli, ale chcę używać gotowego solvera PDE. Ale nawet w prostym przypadku 1D model dyfuzyjno-dyfuzyjny składa się z szeregu sprzężonych nieliniowych …
Aby użyć szybkiej transformacji Fouriera (FFT) na danych o jednakowym próbkowaniu, np. W połączeniu z rozwiązaniami PDE, dobrze wiadomo, że FFT jest algorytmem ). Jak dobrze skala FFT jest przetwarzana równolegle dla n → ∞ (tj. Bardzo duża)?O (nlog( n )O(nlog(n)\mathcal{O}(n\log(n)n → ∞n→∞n\to\infty
Przepływy o dużej liczbie Reynoldsa wytwarzają bardzo cienkie warstwy przyścienne. Jeśli rozdzielczość ściany jest używana w symulacji dużego wiru, współczynnik kształtu może być rzędu . Wiele metod staje się niestabilnych w tym systemie, ponieważ stała inf-sup degraduje się jako pierwiastek kwadratowy współczynnika kształtu lub gorzej. Stała inf-sup jest ważna, ponieważ …
Chciałbym rozwiązać niektóre PDE na rozmaitościach, powiedzmy na przykład równanie eliptyczne na kuli. Gdzie zaczynam? Chciałbym, aby znaleźć coś, że używanie wcześniej istniejących bibliotek / kod 2D, nic tak nadzwyczajnego (na razie) Dodano później: Artykuły i raporty są mile widziane.
Próbuję zrozumieć, jak działa metoda optymalizacji oparta na sprzężeniu dla optymalizacji ograniczonej przez PDE. W szczególności staram się zrozumieć, dlaczego metoda łączenia jest bardziej wydajna w przypadku problemów, w których liczba zmiennych projektowych jest duża, ale „liczba równań jest niewielka”. Co rozumiem: Rozważ następujący problem optymalizacji ograniczonej przez PDE: minβ …
W przypadku niektórych prostych domen wypukłych w 2D mamy spełniające następujące równanie: z pewnymi warunkami brzegowymi Dirichleta i / lub Neumanna. Według mojej wiedzy zastosowanie metody Newtona w przestrzeni elementów skończonych byłoby względnie prostym sposobem numerycznego rozwiązania tego równania.ΩΩ\Omegau ( x )u(x)u(x)- d i v ( A ∇ u ) …
Czy ktoś może mi pomóc znaleźć książki na temat rozwiązań numerycznych (różnic skończonych i metod Cranka – Nicolsona) równań Poissona i dyfuzji, w tym przykładów na nieregularnej geometrii, takich jak dziedzina składająca się z obszaru między prostokątem a okręgiem (zwłaszcza książek lub linków na przykładach kodu MATLAB w tym przypadku)?
Wdrażam artykuł „ Optymalny transport masowy do rejestracji i wypaczania ”, moim celem jest umieszczenie go w Internecie, ponieważ po prostu nie mogę znaleźć żadnego kodu eulerowskiego transportu masowego w Internecie, co byłoby interesujące, przynajmniej dla społeczności naukowej zajmującej się przetwarzaniem obrazu. Artykuł można podsumować w następujący sposób: - znajdź …
Chciałbym się dowiedzieć, jak działa element Raviart-Thomas (RT). W tym celu chciałbym analitycznie opisać, jak funkcje podstawowe wyglądają na kwadracie odniesienia. Celem nie jest wdrożenie go osobiście, ale po prostu intuicyjne zrozumienie elementu. W dużej mierze opieram tę pracę na omawianych tutaj trójkątnych elementach , być może rozszerzenie jej na …
Jestem początkującym z FE. Moje zastosowanie to wycena pochodnych instrumentów finansowych, w których przestrzeń jest pięciowymiarowa. Tak więc, dodając czas, problem ma sześć wymiarów. Próbowałem się rozejrzeć (Fenics, escript, deal.II, ...), ale rozumiem, że te programy są ograniczone do 3 + 1 (przestrzeń 3D + czas 1d). Czy to jest …
Używam schematu różnic skończonych Cranka-Nicolsona do rozwiązania równania cieplnego 1D. Zastanawiam się, czy zasada maksimum / minimum równania ciepła (tj. Że maksimum / minimum występuje w stanie początkowym lub na granicach) również obowiązuje dla rozwiązania dyskretnego. Prawdopodobnie wynika to z faktu, że Crank-Nicolson jest stabilnym i zbieżnym schematem. Wygląda jednak …
Rozwiązuję układ dyfuzyjny reakcji Turinga za pomocą następującego kodu C ++. Jest zbyt wolny: dla tekstury 128 x 128 pikseli dopuszczalna liczba iteracji wynosi 200 - co powoduje 2,5 sekundy opóźnienia. Potrzebuję 400 iteracji, aby uzyskać interesujący obraz - ale 5 sekund oczekiwania to za dużo. Ponadto rozmiar tekstury powinien …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.