# Single-material DG hydrodynamics

### Contents

Inciter supports multiple hydrodynamics schemes. This page describes the DG method for single-material flows.

The discontinuous Galerkin (DG) method implemented in the code is a high order finite element method developed for solving the conservation laws. Similar to classical finite element method, the DG method can achieve high order accuracy by approximating the numerical solution within the element as high order polynomials and admitting discontinuties at the cell interfaces.

## Governing equations of the compressible flows

The governing equations used in the DG algorithm can be represented as

For the scalar advection equation, the scalar variable and are defined as

Here represents the wave speed.

While for the compressible Euler equations, the conservative vector and the flux vector are defined as

Here , , and denote the density, pressure, and specific total energy of the fluid respectively, and is the velocity of the flow in the coordinate direction .The pressure for ideal gas can be computed from the equation of state

where is the ratio of specific heats.

## Discontinuous Galerkin discretization

In order to apply the discontinuous Galerkin discretization, a weak formulation of the governing equations over the computational domain is obtained as

where denotes the boundary of , is the unit outward vector to the boundary and is the test function.

By subdividing the domain into a collection of non-overlapping elements and considering functions U and W defined within each element, we obtain the following semi-discrete formulation,

where denotes the boundary of , and represent the piecewise polynomial approximations to the analytical solution and test function. The Galerkin method assumes the test function to be equal to the basis function. Then the above equation becomes the following system with equations:

The unknows of the conservative variables on each element are defined as

where is the basis functions.

## Basis function

The Dubiner basis function designed for tetrahedron grid is implemented in the DG algortithm. This set of basis function is defined in the reference tetrahdral domain, which helps to achieve the orthogonality features of the basis function. The projection procedure between the physical domain and reference domain is linear. Therefore, each point in the either domain can be found its corresponding point using linear transport relations. Let us consider a third-order approximation in three-dimensional space. In such case, the basis functions are shown as,

## Numerical flux

Due to the discontinuous function approximations, flux terms are not uniquely defined at element interfaces. In this respect, the flux function is replaced by a numerical flux function . In the current code, the approximate Lax-Friedrich, HLL and HLLC type of Riemann solvers are implemented.

## Temporal integration

The semi-discrete form of the governing equations is able to construct a system of ODE equations.

The above system is solved using the TVD-RK3 method with the form of

The TVD-RK3 method is linearly stable for a CFL number less than or equal to ( is the order of the polynomial solution) and the global time step is obtained by finding the minimum value for all the elements in the computational domain.