اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدMulti-Frequency General Relativistic Radiation-Hydrodynamics with $bf{M}_1$ Closure
104
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++, which expands the two-moment, $bf{M}_1$, radiation treatment from grey to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose-Einstein) or neutrino (Fermi-Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation-matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interactions terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.
قيم البحث
اقرأ أيضاً
We present a new radiative transfer method (SPH-M1RT) that is coupled dynamically with smoothed particle hydrodynamics (SPH). We implement it in the (task-based parallel) SWIFT galaxy simulation code but it can be straightforwardly implemented in oth
er SPH codes. Our moment-based method simultaneously solves the radiation energy and flux equations in SPH, making it adaptive in space and time. We modify the M1 closure relation to stabilize radiation fronts in the optically thin limit. We also introduce anisotropic artificial viscosity and high-order artificial diffusion schemes, which allow the code to handle radiation transport accurately in both the optically thin and optically thick regimes. Non-equilibrium thermo-chemistry is solved using a semi-implicit sub-cycling technique. The computational cost of our method is independent of the number of sources and can be lowered further by using the reduced speed of light approximation. We demonstrate the robustness of our method by applying it to a set of standard tests from the cosmological radiative transfer comparison project of Iliev et al. The SPH-M1RT scheme is well-suited for modelling situations in which numerous sources emit ionising radiation, such as cosmological simulations of galaxy formation or simulations of the interstellar medium.
We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, thre
e neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_odot$ progenitor in full unconstrained 3D and in octant symmetry for $gtrsim$$ 380,mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion.
We present an analysis of the general relativistic Boltzmann equation for radiation, appropriate to the case where particles and photons interact through Thomson scattering, and derive the radiation energy-momentum tensor in the diffusion limit, with
viscous terms included. Contrary to relativistic generalizations of the viscous stress tensor that appear in the literature, we find that the stress tensor should contain a correction to the comoving energy density proportional to the divergence of the four-velocity, as well as a finite bulk viscosity. These modifications are consistent with the framework of radiation hydrodynamics in the limit of large optical depth, and do not depend on thermodynamic arguments such as the assignment of a temperature to the zeroth-order photon distribution. We perform a perturbation analysis on our equations and demonstrate that, as long as the wave numbers do not probe scales smaller than the mean free path of the radiation, the viscosity contributes only decaying, i.e., stable, corrections to the dispersion relations. The astrophysical applications of our equations, including jets launched from super-Eddington tidal disruption events and those from collapsars, are discussed and will be considered further in future papers.
We present a method for analyzing the interaction between radiation and matter in regions of intense, relativistic shear that can arise in many astrophysical situations. We show that there is a simple velocity profile that should be manifested in reg
ions of large shear that have lost memory of their boundary conditions, and we use this self-similar velocity profile to construct the surface of last scattering, or $tau simeq 1$ surface, as viewed from any comoving point within the flow. We demonstrate that a simple treatment of scattering from this $tau simeq 1$ surface exactly conserves photon number, and derive the rate at which the radiation field is heated due to the shear present in the flow. The components of the comoving radiation energy-momentum tensor are calculated, and we show that they have relatively simple, approximate forms that interpolate between the viscous (small shear) and streaming (large shear) limits. We put our expression for the energy-momentum tensor in a covariant form that does not depend on the explicit velocity profile within the fluid and, therefore, represents a natural means for analyzing general, radiation-dominated, relativistic shear flows.
Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant applica
tion using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20%$, the large-scale electromagnetic jet efficiency is of order $10%$, and the total radiative efficiency that reaches large distances remains low at only order $1%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs. disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $zsim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
