skip to main content

Title: A high-order shock capturing discontinuous Galerkin–finite difference hybrid method for GRMHD
Abstract We present a discontinuous Galerkin (DG)–finite difference (FD) hybrid scheme that allows high-order shock capturing with the DG method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited DG candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust FD methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the DG and FD solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1D, 2D, and 3D relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.  more » « less
Award ID(s):
1931280 2207342 2209655 1931266 2209656 2011961
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Classical and Quantum Gravity
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Summary

    A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.

    more » « less
  2. Abstract

    This article presents high order accurate discontinuous Galerkin (DG) methods for wave problems on moving curved meshes with general choices of basis and quadrature. The proposed method adopts an arbitrary Lagrangian–Eulerian formulation to map the wave equation from a time‐dependent moving physical domain onto a fixed reference domain. For moving curved meshes, weighted mass matrices must be assembled and inverted at each time step when using explicit time‐stepping methods. We avoid this step by utilizing an easily invertible weight‐adjusted approximation. The resulting semi‐discrete weight‐adjusted DG scheme is provably energy stable up to a term that (for a fixed time interval) converges to zero with the same rate as the optimal error estimate. Numerical experiments using both polynomial and B‐spline bases verify the high order accuracy and energy stability of proposed methods.

    more » « less
  3. High power sources of electromagnetic energy often require complicated structures to support electromagnetic modes and shape electromagnetic fields to maximize the coupling of the field energy to intense relativistic electron beams. Geometric fidelity is critical to the accurate simulation of these High Power Electromagnetic (HPEM) sources. Here, we present a fast and geometrically flexible approach to calculate the solution to Maxwell’s equations in vector potential form under the Lorenz gauge. The scheme is an implicit, linear-time, high-order, A-stable method that is based on the method of lines transpose (MOLT). As presented, the method is fourth order in time and second order in space, but the A-stable formulation could be extended to both high order in time and space. An O(n) fast convolution is employed for space-integration. The main focus of this work is to develop an approach to impose perfectly electrically conducting (PEC) boundary conditions in MOLT by extending our past work on embedded boundary methods. As the method is A-stable, it does not suffer from small time step limitations that are found in explicit finite difference time domain methods when using either embedded boundary or cut-cell methods to capture geometry. This is a major advance for the simulation of HPEM devices. While there is no conceptual limitation to develop this in 3D, our initial work has centered on 2D. The extension to 3D requires validation that the proposed fixed point iteration will converge and is the subject of our follow-up work. The eventual goal is to combine this method with particle methods for the simulations of plasma. In the current work, the scheme is evaluated for EM wave propagation within an object that is bounded by PEC. The consistency and performance of the scheme are confirmed using the ping test and frequency mode analysis for rotated square cavities—a standard test in the HPEM community. We then demonstrate the diffraction Q value test and the use of this method for simulating an A6 magnetron. The ability to handle both PEC and open boundaries in a standard device test problem, such as the A6, gives confidence on the robustness of this new method.

    more » « less
  4. Aims . In this work, we study the structure of neutron stars under the effect of a poloidal magnetic field and determine the limiting largest magnetic field strength that induces a deformation such that the ratio between the polar and equatorial radii does not exceed 2%. We consider that, under these conditions, the description of magnetic neutron stars in the spherical symmetry regime is still satisfactory. Methods . We described different compositions of stars (nucleonic, hyperonic, and hybrid) using three state-of-the-art relativistic mean field models (NL3 ω ρ , MBF, and CMF, respectively) for the microscopic description of matter, all in agreement with standard experimental and observational data. The structure of stars was described by the general relativistic solution of both Einstein’s field equations assuming spherical symmetry and Einstein-Maxwell’s field equations assuming an axi-symmetric deformation. Results . We find a limiting magnetic moment on the order of 2 × 10 31 Am 2 , which corresponds to magnetic fields on the order of 10 16 G at the surface and 10 17 G at the center of the star, above which the deformation due to the magnetic field is above 2%, and therefore not negligible. We show that the intensity of the magnetic field developed in the star depends on the equation of state (EoS), and, for a given baryonic mass and fixed magnetic moment, larger fields are attained with softer EoS. We also show that the appearance of exotic degrees of freedom, such as hyperons or a quark core, is disfavored in the presence of a very strong magnetic field. As a consequence, a highly magnetized nucleonic star may suffer an internal conversion due to the decay of the magnetic field, which could be accompanied by a sudden cooling of the star or a gamma ray burst. 
    more » « less
  5. Trefftz methods are high-order Galerkin schemes in which all discrete functions are elementwise solution of the PDE to be approximated. They are viable only when the PDE is linear and its coefficients are piecewise-constant. We introduce a “quasi-Trefftz” discontinuous Galerkin (DG) method for the discretisation of the acoustic wave equation with piecewise-smooth material parameters: the discrete functions are elementwise approximate PDE solutions. We show that the new discretisation enjoys the same excellent approximation properties as the classical Trefftz one, and prove stability and high-order convergence of the DG scheme. We introduce polynomial basis functions for the new discrete spaces and describe a simple algorithm to compute them. The technique we propose is inspired by the generalised plane waves previously developed for time-harmonic problems with variable coefficients; it turns out that in the case of the time-domain wave equation under consideration the quasi-Trefftz approach allows for polynomial basis functions. 
    more » « less