Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity
484
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A considerable amount of attention has been given to discontinuous Galerkin methods for hyperbolic problems in numerical relativity, showing potential advantages of the methods in dealing with hydrodynamical shocks and other discontinuities. This paper investigates discontinuous Galerkin methods for the solution of elliptic problems in numerical relativity. We present a novel hp-adaptive numerical scheme for curvilinear and non-conforming meshes. It uses a multigrid preconditioner with a Chebyshev or Schwarz smoother to create a very scalable discontinuous Galerkin code on generic domains. The code employs compactification to move the outer boundary near spatial infinity. We explore the properties of the code on some test problems, including one mimicking Neutron stars with phase transitions. We also apply it to construct initial data for two or three black holes.
rate research
Read More
We present a discontinuous Galerkin internal-penalty scheme that is applicable to a large class of linear and non-linear elliptic partial differential equations. The scheme constitutes the foundation of the elliptic solver for the SpECTRE numerical relativity code. As such it can accommodate (but is not limited to) elliptic problems in linear elasticity, general relativity and hydrodynamics, including problems formulated on a curved manifold. We provide practical instructions that make the scheme functional in a production code, such as instructions for imposing a range of boundary conditions, for implementing the scheme on curved and non-conforming meshes and for ensuring the scheme is compact and symmetric so it may be solved more efficiently. We report on the accuracy of the scheme for a suite of numerical test problems.
Understanding fundamental kinetic processes is important for many problems, from plasma physics to gas dynamics. A first-principles approach to these problems requires a statistical description via the Boltzmann equation, coupled to appropriate field equations. In this paper we present a novel version of the discontinuous Galerkin (DG) algorithm to solve such kinetic equations. Unlike Monte-Carlo methods we use a continuum scheme in which we directly discretize the 6D phase-space using discontinuous basis functions. Our DG scheme eliminates counting noise and aliasing errors that would otherwise contaminate the delicate field-particle interactions. We use modal basis functions with reduced degrees of freedom to improve efficiency while retaining a high formal order of convergence. Our implementation incorporates a number of software innovations: use of JIT compiled top-level language, automatically generated computational kernels and a sophisticated shared-memory MPI implementation to handle velocity space parallelization.
We study a class of nonlinear eigenvalue problems of Scrodinger type, where the potential is singular on a set of points. Such problems are widely present in physics and chemistry, and their analysis is of both theoretical and practical interest. In particular, we study the regularity of the eigenfunctions of the operators considered, and we propose and analyze the approximation of the solution via an isotropically refined hp discontinuous Galerkin (dG) method. We show that, for weighted analytic potentials and for up-to-quartic nonlinearities, the eigenfunctions belong to analytic-type non homogeneous weighted Sobolev spaces. We also prove quasi optimal a priori estimates on the error of the dG finite element method; when using an isotropically refined hp space the numerical solution is shown to converge with exponential rate towards the exact eigenfunction. In addition, we investigate the role of pointwise convergence in the doubling of the convergence rate for the eigenvalues with respect to the convergence rate of eigenfunctions. We conclude with a series of numerical tests to validate the theoretical results.
We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $fin[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.
We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for regularization of shocks. High order multi-stage forward Euler and strong stability preserving Runge-Kutta time integration options complement high order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, though we note an equivalent capability currently also exists in CosmosDG for Newtonian systems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
