Eu gostaria de um código computacional super simples para resolver problemas de magneto-hidrodinâmica . Alta precisão nem desempenho não é minha preocupação. Não quero apenas explorar visualmente o comportamento qualitativo do plasma. Importante para mim é:
- Generalidade - não deve ser especializada em nenhum subproblema específico de MHD
- Fácil de fazer funcionar - não requer bibliotecas difíceis de instalar. A entrada é direta e bem documentada por alguns exemplos
- Fácil de ler e modificar - você pode ler o código e aprender sobre a solução MHD. Você pode usar esse código como ponto de partida para o seu próprio desenvolvimento de código MHD
Algo tão simples como este Navier Stokes Solver seria o melhor
melhor seria algo em python-numpy ou java, mas C ++ ou Fortran não é um grande problema para mim.
Até agora eu encontrei:
algumas bibliotecas ( MHD-Hermes e python-mhd ) que estão em python, mas parece não ter nenhuma documentação ou exemplo de como usá-lo.
e este Um código MTV TVD MHD gratuito, rápido, simples e eficiente que está em frotran, mas tem apenas 400 linhas
Eu estava pensando em escrever algo sozinho, mas quero ter certeza de que não há nada semelhante já disponível. Eu acho que poderia ser bem fácil, basta pegar o solucionador Navier-Stokes no link acima e introduzir o potencial elétrico e magnético pelo solucionador de poisson FFT (ou Particle Mesh Ewald).