Me gustaría un código computacional súper simple para resolver problemas de magnetohidrodinámica . La alta precisión ni el rendimiento no son mi preocupación. Solo quiero explorar visualmente el comportamiento cualitativo del plasma. Importante para mí es:
- Generalidad : no debe estar especializado en ningún subproblema particular de MHD
- Fácil de hacer que funcione : no requiere ninguna biblioteca difícil de instalar. La entrada es sencilla y está bien documentada por algunos ejemplos.
- Fácil de leer y modificar : puede leer el código y aprender sobre la solución MHD a partir de él. Puede usar este código como punto de partida para su propio desarrollo de código MHD
Algo tan simple como este Navier aviva el solucionador sería lo mejor
lo mejor sería algo en python-numpy o java, pero C ++ o Fortran no es un gran problema para mí.
Hasta ahora encontré:
algunas bibliotecas ( MHD-Hermes y python-mhd ) que están en python, pero no parecen tener documentación ni ejemplos de cómo usarlo.
y este código TVD MHD gratuito, rápido, simple y eficiente que está en frotran pero tiene solo 400 líneas
Estaba pensando en escribir algo por mí mismo, pero quiero estar seguro de que ya no hay nada similar disponible. Creo que podría ser bastante fácil, solo tome el solucionador Navier-Stokes del enlace de arriba e introduzca el potencial eléctrico y magnético mediante el solucionador de poisson FFT (o Particle Mesh Ewald).