Pytania otagowane jako advection-diffusion

1
Jak debugować kod numeryczny, co może być źródłem tego błędu oscylacyjnego?
Cichą wiedzę można uzyskać z doświadczenia, zastanawiałem się tylko, czy ktoś wcześniej widział coś podobnego do tego. Wykres pokazuje początkowy warunek (zielony) dla równania rada-dyfuzja, następnie roztwór przy iteracji 200 (niebieski), a następnie ponownie przy iteracji 400 (czerwony). Rozwiązanie równania dyfuzyjno-doradczego pojawia się po kilku iteracjach. Liczba Pécleta i warunek …
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 / …
3
Jakie są podstawowe zasady generowania ruchomej siatki?
Jestem zainteresowany wdrożeniem ruchomej siatki w celu rozwiązania problemu dyfuzyjnego. Metody adaptacyjnej siatki ruchomej stanowią dobry przykład tego, jak to zrobić dla równania Burgera w 1D przy użyciu skończonej różnicy. Czy ktoś byłby w stanie podać praktyczny przykład rozwiązania równania doradczego z dyfuzją 1D przy użyciu różnicy skończonej z ruchomą …
3
Metody rozwiązywania nieliniowych układów doradczo-dyfuzyjnych poza Newton-Raphson?
Pracuję nad projektem, w którym mam dwie domeny sprzężone za pomocą adv-diff za pomocą odpowiednich terminów źródłowych (jedna domena dodaje masę, druga odejmuje masę). Dla zwięzłości modeluję je w stanie ustalonym. Równania to standardowe równanie transportowo-dyfuzyjne z terminem źródłowym wyglądającym tak: ∂c1∂t=0=F1+Q1(c1,c2)∂c2∂t=0=F2+Q2(c1,c2)∂c1∂t=0=F1+Q1(c1,c2)∂c2∂t=0=F2+Q2(c1,c2) \frac{\partial c_1}{\partial t} = 0 = \mathcal{F}_1 + …
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.