PrincetonUniversity / athena-public-version

(MOVED) Athena++ GRMHD code and adaptive mesh refinement (AMR) framework. Active repository --->
https://github.com/PrincetonUniversity/athena
BSD 3-Clause "New" or "Revised" License
160 stars 118 forks source link

A request for a relativistic turbulence driver #47

Closed jianfuzhangxtu closed 3 years ago

jianfuzhangxtu commented 4 years ago

Dear Developers and Users, I am trying to employ Athena++ to simulate a driven relativistic MHD turbulence.

On the basis of the problem generator file turb.cpp (Athena++/src/pgen/), I added an initialization of magnetic fields, which can correctly evolve in the case of the non-relativistic MHD turbulence. When configuring with “configure.py -fft -b -mpi --prob turb -b –s -hdf5 –mpi”, my attempting run with changing input parameters can only evolve a very short time ~0.1, printing “[perturb]: s is NaN!”. Although I am not sure whether an initialization of a gas pressure is necessary, this was done as well in order for describing an adiabatic state.

A bigger issue is that the turbulence driver module (Athena++/src/fft/turbulenc.cpp) was only written for Newtonian physics, which needs to be changed to work in the case of special relativity.

This post issue is to seek a possible help. Whether does the developer or user have the relativistic MHD turbulence driver? Are you willing to provide assistance to change the driver module, or share it.

With my thanks and best regards,

Jianfu