Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRealizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
62
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Building on the framework of Zhang & Shu cite{zhangShu_2010a,zhangShu_2010b}, we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function $f$. The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky & Bludman cite{cernohorskyBludman_1994} and Banach & Larecki cite{banachLarecki_2017a}. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Paulis exclusion principle, which demands a bounded distribution function (i.e., $fin[0,1]$). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. Numerical results demonstrate the realizability-preserving properties of the scheme. We also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.
rate research
Read More
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 explosions. 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.
Recently, a 4th-order asymptotic preserving multiderivative implicit-explicit (IMEX) scheme was developed (Schutz and Seal 2020, arXiv:2001.08268). This scheme is based on a 4th-order Hermite interpolation in time, and uses an approach based on operator splitting that converges to the underlying quadrature if iterated sufficiently. Hermite schemes have been used in astrophysics for decades, particularly for N-body calculations, but not in a form suitable for solving stiff equations. In this work, we extend the scheme presented in Schutz and Seal 2020 to higher orders. Such high-order schemes offer advantages when one aims to find high-precision solutions to systems of differential equations containing stiff terms, which occur throughout the physical sciences. We begin by deriving Hermite schemes of arbitrary order and discussing the stability of these formulas. Afterwards, we demonstrate how the method of Schutz and Seal 2020 generalises in a straightforward manner to any of these schemes, and prove convergence properties of the resulting IMEX schemes. We then present results for methods ranging from 6th to 12th order and explore a selection of test problems, including both linear and nonlinear ordinary differential equations and Burgers equation. To our knowledge this is also the first time that Hermite time-stepping methods have been applied to partial differential equations. We then discuss some benefits of these schemes, such as their potential for parallelism and low memory usage, as well as limitations and potential drawbacks.
We introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.
We design, analyze and numerically validate a novel discontinuous Galerkin method for solving the coagulation-fragmentation equations. The DG discretization is applied to the conservative form of the model, with flux terms evaluated by Gaussian quadrature with $Q=k+1$ quadrature points for polynomials of degree $k$. The positivity of the numerical solution is enforced through a simple scaling limiter based on positive cell averages. The positivity of cell averages is propagated by the time discretization provided a proper time step restriction is imposed.
We present a finite element based variational interface-preserving and conservative phase-field formulation for the modeling of incompressible two-phase flows with surface tension dynamics. The preservation of the hyperbolic tangent interface profile of the convective Allen-Cahn phase-field formulation relies on a novel time-dependent mobility model. The mobility coefficient is adjusted adaptively as a function of gradients of the velocity and the order parameter in the diffuse interface region in such a way that the free energy minimization properly opposes the convective distortion. The ratio of the convective distortion to the free energy minimization is termed as the convective distortion parameter, which characterizes the deviation from the hyperbolic tangent shape due to the convection effect. The mass conservation is achieved by enforcing a Lagrange multiplier with both temporal and spatial dependence on the phase-field function. We integrate the interface-preserving and conservative phase-field formulation with the incompressible Navier-Stokes equations and the continuum surface tension force model for the simulation of incompressible two-phase flows. A positivity preserving scheme designed for the boundedness and stability of the solution is employed for the variational discretization using unstructured finite elements. We examine the convergence and accuracy of the Allen-Cahn phase-field solver through a generic one-dimensional bistable diffusion-reaction system in a stretching flow. We quantify and systematically assess the relative interface thickness error and the relative surface tension force error with respect to the convective distortion parameter. Two- and three-dimensional rising bubble cases are further simulated to examine the effectiveness of the proposed model on the volume-preserving mean curvature flow and the interface-preserving capability.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
