اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدConservative Moment Equations for Neutrino Radiation Transport with Limited Relativity
128
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We derive conservative, multidimensional, energy-dependent moment equations for neutrino transport in core-collapse supernovae and related astrophysical systems, with particular attention to the consistency of conservative four-momentum and lepton number transport equations. After taking angular moments of conservative formulations of the general relativistic Boltzmann equation, we specialize to a conformally flat spacetime, which also serves as the basis for four further limits. Two of these---the multidimensional special relativistic case, and a conformally flat formulation of the spherically symmetric general relativistic case---are given in appendices for the sake of comparison with extant literature. The third limit is a weak-field, `pseudo-Newtonian approach citep{kim_etal_2009,kim_etal_2012} in which the source of the gravitational potential includes the trace of the stress-energy tensor (rather than just the mass density), and all orders in fluid velocity $v$ are retained. Our primary interest here is in the fourth limit: `$mathcal{O}(v)$ moment equations for use in conjunction with Newtonian self-gravitating hydrodynamics. We show that the concept of `$mathcal{O}(v)$ transport requires care when dealing with both conservative four-momentum and conservative lepton number transport, and present two self-consistent options: `$mathcal{O}(v)$-plus transport, in which an $mathcal{O}(v^2)$ energy equation combines with an $mathcal{O}(v)$ momentum equation to give an $mathcal{O}(v^2)$ number equation; and `$mathcal{O}(v)$-minus transport, in which an $mathcal{O}(v)$ energy equation combines with an $mathcal{O}(1)$ momentum equation to give an $mathcal{O}(v)$ number equation.
قيم البحث
اقرأ أيضاً
Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy t
o extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.
Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for neutrino-driven explos
ions. Due to the multiscale nature of neutrino transport in CCSN simulations, an implicit treatment of neutrino-matter interactions is desired, which requires solutions of coupled nonlinear systems in each step of the time integration scheme. In this paper we design and compare nonlinear iterative solvers for implicit systems with energy coupling neutrino-matter interactions commonly used in CCSN simulations. Specifically, we consider electron neutrinos and antineutrinos, which interact with static matter configurations through the Bruenn~85 opacity set. The implicit systems arise from the discretization of a non-relativistic two-moment model for neutrino transport, which employs the discontinuous Galerkin (DG) method for phase-space discretization and an implicit-explicit (IMEX) time integration scheme. In the context of this DG-IMEX scheme, we propose two approaches to formulate the nonlinear systems -- a coupled approach and a nested approach. For each approach, the resulting systems are solved with Anderson-accelerated fixed-point iteration and Newtons method. The performance of these four iterative solvers has been compared on relaxation problems with various degree of collisionality, as well as proto-neutron star deleptonization problems with several matter profiles adopted from spherically symmetric CCSN simulations. Numerical results suggest that the nested Anderson-accelerated fixed-point solver is more efficient than other tested solvers for solving implicit nonlinear systems with energy coupling neutrino-matter interactions.
We analyze the mixed frame equations of radiation hydrodynamics under the approximations of flux-limited diffusion and a thermal radiation field, and derive the minimal set of evolution equations that includes all terms that are of leading order in a
ny regime of non-relativistic radiation hydrodynamics. Our equations are accurate to first order in v/c in the static diffusion regime. In contrast, we show that previous lower order derivations of these equations omit leading terms in at least some regimes. In comparison to comoving frame formulations of radiation hydrodynamics, our equations have the advantage that they manifestly conserve total energy, making them very well-suited to numerical simulations, particularly with adaptive meshes. For systems in the static diffusion regime, our analysis also suggests an algorithm that is both simpler and faster than earlier comoving frame methods. We implement this algorithm in the Orion adaptive mesh refinement code, and show that it performs well in a range of test problems.
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. Th
e code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our codes capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
