Ich hätte gerne einen supereinfachen Rechencode zur Lösung von Problemen der Magnetohydrodynamik . Hohe Genauigkeit und Leistung sind nicht mein Anliegen. Ich möchte nicht nur das qualitative Verhalten von Plasma visuell untersuchen. Wichtig für mich ist:
- Allgemeines - Es sollte nicht auf ein bestimmtes Teilproblem von MHD spezialisiert sein
- Einfach zum Laufen zu bringen - erfordert keine schwer zu installierenden Bibliotheken. Die Eingabe ist unkompliziert und anhand einiger Beispiele gut dokumentiert
- Einfach zu lesen und zu ändern - Sie können den Code lesen und sich über die MHD-Lösung informieren. Sie können diesen Code als Ausgangspunkt für Ihre eigene MHD-Codeentwicklung verwenden
Etwas so Einfaches wie dieser Navier Stokes Solver wäre am besten
Am besten wäre etwas in Python-Numpy oder Java, aber C ++ oder Fortran ist für mich kein großes Problem.
Bisher habe ich gefunden:
Einige Bibliotheken ( MHD-Hermes und Python-mhd ) befinden sich in Python, scheinen jedoch keine Dokumentation oder Beispiele für deren Verwendung zu haben.
und dies ist ein kostenloser, schneller, einfacher und effizienter TVD-MHD-Code, der in Frotran vorliegt, aber nur 400 Zeilen enthält
Ich habe darüber nachgedacht, etwas selbst zu schreiben, aber ich möchte sichergehen, dass nichts Ähnliches bereits verfügbar ist. Ich denke, es könnte ganz einfach sein, nehmen Sie einfach den Navier-Stokes-Löser aus dem obigen Link und führen Sie das elektrische und magnetische Potential durch den FFT-Poisson-Löser (oder Particle Mesh Ewald) ein.