期刊
ASTROPHYSICAL JOURNAL
卷 837, 期 1, 页码 -出版社
IOP PUBLISHING LTD
DOI: 10.3847/1538-4357/aa5f55
关键词
accretion, accretion disks; magnetohydrodynamics (MHD); plasmas; relativistic processes
资金
- NSF
- NASA through an Einstein Postdoctoral Fellowship [PF4-150122]
- Chandra X-ray Center [NAS8-03060]
- Illinois Distinguished Fellowship from the University of Illinois
- NSF [AST-333612]
- visiting fellowship at All Souls College, Oxford
- Direct For Mathematical & Physical Scien [1333682] Funding Source: National Science Foundation
- Division Of Astronomical Sciences [1333682] Funding Source: National Science Foundation
- Division Of Astronomical Sciences
- Direct For Mathematical & Physical Scien [1333091] Funding Source: National Science Foundation
- Division Of Astronomical Sciences
- Direct For Mathematical & Physical Scien [1333612] Funding Source: National Science Foundation
Hot, diffuse, relativistic plasmas such as sub-Eddington black-hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects such as heat conduction and viscosity can be important in a collisionless plasma and will potentially alter the dynamics and radiative properties of the flow from that in ideal fluid models; we refer to models that include these processes as Extended GRMHD. Here we describe a new conservative code, grim,. that enables all of the above and additional physics to be efficiently incorporated. grim. combines time evolution and primitive variable inversion needed for conservative schemes into a single step using an algorithm that only requires the residuals of the governing equations as inputs. This algorithm enables the code to be physics agnostic as well as flexibility. regarding time-stepping schemes. grim. runs on CPUs, as well as on GPUs, using the same code. We formulate a performance model. and use it to show that our implementation runs optimally on both architectures. grim. correctly captures classical GRMHD test problems as well as a new suite of linear and nonlinear test problems with anisotropic conduction and viscosity in special and general relativity. As tests and example applications, we resolve the shock substructure due to the presence of dissipation, and report on relativistic versions of the magneto-thermal instability and heat flux driven buoyancy instability, which arise due to anisotropic heat conduction, and of the firehose instability, which occurs due to anisotropic pressure (i.e., viscosity). Finally, we show an example integration of an accretion flow around a Kerr black hole, using Extended GRMHD.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据