A set of codes for solving 1D/2D/3D time-dependent hyperbolic PDEs using high order DG methods on conforming unstructured meshes of triangular, quadrilateral, or hexahedral elements.
These codes are inspired by the Matlab codes for the book Nodal Discontinuous Galerkin methods by Hesthaven and Warburton (2007). While unstructured meshes are supported, all demos use uniform meshes, and the codebase is intended mainly for experimentation and method development.
- The simplest demo is "dg1D_advec.jl", which solves the 1D advection equation on a uniform 1D mesh.
- The demo "dg2D_advec_tri.jl" simulates the advection equation on a triangular mesh.
- The demo "dg2D_wave_quad.jl" computes solutions to the acoustic wave equation on quadrilateral meshes.
- The demo "dg3D_advec_hex.jl" solves the advection equation on a hexahedral mesh.
- The file "dg2D_euler_quad.jl" and "dg3D_euler_hex.jl" provide entropy stable DG methods on quadrilateral and hexahedral meshes for the compressible Euler equations, with a sparsity-optimized implementation of the Hadamard sum step in flux differencing.
The discretizations used are based on the following references: