Was ist die theoretische Konvergenzrate für einen FFT-Giftlöser? Ich löse eine Poisson-Gleichung: ∇2VH(x,y,z)=−4πn(x,y,z)∇2VH(x,y,z)=-4πn(x,y,z)\nabla^2 V_H(x, y, z) = -4\pi n(x, y, z) mit in der Domäne[0,2]×[0,2]×[0,2]mit periodischer Randbedingung. Diese Ladungsdichte ist netto neutral. Die...
Kürzlich habe ich verschiedene nichtlineare Löser von scipy verglichen und war besonders beeindruckt vom Newton-Krylov-Beispiel im Scipy-Kochbuch, in dem sie eine Differentialgleichungsgleichung zweiter Ordnung mit nichtlinearem Reaktionsterm in etwa 20 Codezeilen lösen. Ich habe den Beispielcode...
Ich bin daran interessiert, die Poisson-Gleichung mit dem Finite-Differenzen-Ansatz zu lösen. Ich möchte besser verstehen, wie man die Matrixgleichung mit Neumann-Randbedingungen schreibt. Würde jemand das Folgende überprüfen, ist es richtig? Die Finite-Differenz-Matrix Die...
Betrachten Sie ein symmetrisches positives definitives tridiagonales lineares System wobei A ∈ R n × n und b ∈ R n sind . Wenn wir bei drei Indizes 0 ≤ i < j < k < n nur Gleichungszeilen annehmen, die streng zwischen i und k gelten, können wir Zwischenvariablen eliminieren, um eine...
Ich würde gerne wissen, wie Dirichlet-Bedingungen normalerweise angewendet werden, wenn die Methode des endlichen Volumens auf einem zellzentrierten ungleichmäßigen Gitter angewendet wird. Meine aktuelle Implementierung legt einfach die Randbedingung fest, dass ich den Wert der ersten Zelle...
Ich habe gehört, dass eine schnelle Fourier-Transformation verwendet werden kann, um das Poisson-Problem zu lösen, wenn die Randbedingungen alle ein Typ sind ... Sinusreihen für Dirichlet, Cosinus für Neumann und beide für periodische. Angenommen, zwei gegenüberliegende Seiten haben periodische...
Ich habe in den letzten Tagen versucht, diesen Fehler zu beheben. Ich habe mich gefragt, ob jemand Ratschläge zum weiteren Vorgehen hat. Ich löse die Poisson-Gleichung für eine Stufenladungsverteilung (ein häufiges Problem in der Elektrostatik / Halbleiterphysik) auf einem ungleichmäßigen...
Problemstellung Ich habe ein geometrisches Multigrid für implementiert, wobei f = 3 π 2 ist- ∇2= f- -∇2=f-\nabla^{2}=f aufΩ∈[0,1]auf einemEinheitswürfel. Die Dirichlet-Grenzen auf der linken Seite, der Unterseite und der Vorderseite sind0. Neumann-Grenzen oben, rechts und hinten sind∂uf= 3 π24s i n...