我想要一些超级简单的计算代码来解决磁流体动力学问题。高精度和性能不是我关心的问题。我不希望它只是直观地探索等离子体的定性行为。对我来说重要的是：

• 通用性 - 它不应该专门针对 MHD 的任何特定子问题
• 易于使用- 不需要任何难以安装的库。输入是直截了当的，并有一些例子有据可查
• 易于阅读和修改- 您可以阅读代码并从中了解 MHD 解决方案。您可以将此代码用作您自己的 MHD 代码开发的起点

像这个Navier stokes 求解器这样简单的东西是最好的

最好是 python-numpy 或 java 中的东西，但 C++ 或 Fortran 对我来说不是一个大问题。

到目前为止，我发现：

一些库（MHD-Hermes和python-mhd）在 python 中，但似乎没有任何文档或示例如何使用它。

这 是一个免费、快速、简单和高效的 TVD MHD 代码，它在 frotran 中，但只有 400 行

我正在考虑自己写一些东西，但我想确定没有类似的东西已经可用。我认为这很容易，只需从上面的链接中获取 Navier-Stokes 求解器，并通过 FFT 泊松求解器（或 Particle Mesh Ewald ）引入电势和磁势。