| Authors : | Kushagra Sharma1, Abhishek Kumar Srivastava1, José Juan González-Aviles2
1 Department of Physics, Indian Institute of Technology (BHU), Varanasi 221005, India.
2 Cátedras CONACYT, SCIESMEX-LANCE, Instituto de Geofísica, Unidad Michoacán, Universidad Nacional Autónoma de México. Morelia, Michoacán, México
|
| Abstract : | We have developed a python-based finite volume numerical code to solve ideal Magnetohydrodynamic(MHD) equations. The code uses dimensionally unsplit CTU method of Collela (1990) for integration of MHD equations and maintaining the divergence free condition of the magnetic field. Also it consists of a multi state HLL approximate Riemann solver (HLL, HLLC, HLLD) as adopted from Miyoshi & Kusano (2005) for the flux calculations, with piecewise constant reconstruction. The code has successfully passed the Orszag-Tang Vortex test as demonstrated firstly by Orszag & Tang (1998), as well as shock-tube test. The aim of this code is to model and study the dynamics of various solar and astrophysical plasma systems in the frame-work of MHD. In the present paper, we will demonstrate the basic structure of the MHD code, it’s test runs, and applications to the solar atmosphere.
References -
Colella P., 1990, Journal of Computational Physics, 87, 171
Miyoshi T., Kusano K., 2005, Journal of Computational Physics, 208, 315
Orszag S. A., Tang C.-M., 1979, Journal of Fluid Mechanics, 90, 129 |
