Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
157
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite volume (FV) approach. We demonstrate that the DG technique offers distinct advantages over FV formulations on both static and moving meshes. The DG method is also easily generalized to higher than second-order accuracy without requiring the use of extended stencils to estimate derivatives (thereby making the scheme highly parallelizable). We implement the technique in the AREPO code for solving the fluid and the magnetohydrodynamic (MHD) equations. By examining various test problems, we show that our new formulation provides improved accuracy over FV approaches of the same order, and reduces post-shock oscillations and artificial diffusion of angular momentum. In addition, the DG method makes it possible to represent magnetic fields in a locally divergence-free way, improving the stability of MHD simulations and moderating global divergence errors, and is a viable alternative for solving the MHD equations on meshes where Constrained-Transport (CT) cannot be applied. We find that the DG procedure on a moving mesh is more sensitive to the choice of slope limiter than is its FV method counterpart. Therefore, future work to improve the performance of the DG scheme even further will likely involve the design of optimal slope limiters. As presently constructed, our technique offers the potential of improved accuracy in astrophysical simulations using the moving mesh AREPO code as well as those employing adaptive mesh refinement (AMR).
rate research
Read More
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.
In this paper, a high order quasi-conservative discontinuous Galerkin (DG) method using the non-oscillatory kinetic flux is proposed for the 5-equation model of compressible multi-component flows with Mie-Gruneisen equation of state. The method mainly consists of three steps: firstly, the DG method with the non-oscillatory kinetic flux is used to solve the conservative equations of the model; secondly, inspired by Abgralls idea, we derive a DG scheme for the volume fraction equation which can avoid the unphysical oscillations near the material interfaces; finally, a multi-resolution WENO limiter and a maximum-principle-satisfying limiter are employed to ensure oscillation-free near the discontinuities, and preserve the physical bounds for the volume fraction, respectively. Numerical tests show that the method can achieve high order for smooth solutions and keep non-oscillatory at discontinuities. Moreover, the velocity and pressure are oscillation-free at the interface and the volume fraction can stay in the interval [0,1].
In this paper we present energy-conserving, mixed discontinuous Galerkin (DG) and continuous Galerkin (CG) schemes for the solution of a broad class of physical systems described by Hamiltonian evolution equations. These systems often arise in fluid mechanics (incompressible Euler equations) and plasma physics (Vlasov--Poisson equations and gyrokinetic equations), for example. The dynamics is described by a distribution function that evolves given a Hamiltonian and a corresponding Poisson bracket operator, with the Hamiltonian itself computed from field equations. Hamiltonian systems have several conserved quantities, including the quadratic invariants of total energy and the $L_2$ norm of the distribution function. For accurate simulations one must ensure that these quadratic invariants are conserved by the discrete scheme. We show that using a discontinuous Galerkin scheme to evolve the distribution function and ensuring that the Hamiltonian lies in its continuous subspace leads to an energy-conserving scheme in the continuous-time limit. Further, the $L_2$ norm is conserved if central fluxes are used to update the distribution function, but decays monotonically when using upwind fluxes. The conservation of density and $L_2$ norm is then used to show that the entropy is a non-decreasing function of time. The proofs shown here apply to any Hamiltonian system, including ones in which the Poisson bracket operator is non-canonical (for example, the gyrokinetic equations). We demonstrate the ability of the scheme to solve the Vlasov--Poisson and incompressible Euler equations in 2D and provide references where we have applied these schemes to solve the much more complex 5D electrostatic and electromagnetic gyrokinetic equations.
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.
We investigate the applicability of curvilinear grids in the context of astrophysical simulations and WENO schemes. With the non-smooth mapping functions from Calhoun et al. (2008), we can tackle many astrophysical problems which were out of scope with the standard grids in numerical astrophysics. We describe the difficulties occurring when implementing curvilinear coordinates into our WENO code, and how we overcome them. We illustrate the theoretical results with numerical data. The WENO finite difference scheme works only for high Mach number flows and smooth mapping functions whereas the finite volume scheme gives accurate results even for low Mach number flows and on non-smooth grids.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
