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.
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
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
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,
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.
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.
Tab / T to search, Esc to close
Search for symbols, directories, files, pages or modules. You can omit any
prefix from the symbol or file path; adding a : or /
suffix lists all members of given symbol or directory. Navigate through the
list using ↓ and
Enter to go.