These html pages are based on the PhD thesis "Cluster-Based Parallelization of Simulations on Dynamically Adaptive Grids and Dynamic Resource Management" by Martin Schreiber.

There is also more information and a PDF version available.

There is also more information and a PDF version available.

So far, we discretized our problem in space by decomposing the domain into non-overlapping cells _{i}.
The simulation itself can be executed with cell-local computations and computations based on
adjacent cells. However, for the domain boundary ∂Ω, particular boundary conditions have to be
constructed and applied. Considering our flux computations, we need DoF at the edges shared by both
adjacent cells. Our approach is based on the construction of a ghost edge layer based on the data of
the cell which is sharing the ghost edge.

We can apply Dirichlet boundary conditions directly by setting the DoF on the edge to the Dirichlet parameters instead of recomputing them based on cell information adjacent to the boundary edge. Those boundary conditions are of particular interest for simulations of analytic benchmarks with changing inflow conditions in each time step (see Section 6.2).

One possible way of handling outflow boundary conditions is to set the conserved quantities at the ghost edges to be equal to those from the adjacent cell. However, one severe drawback exists by directly applying this outflow boundary condition. In case of a flow into the domain (e.g., due to a negative momentum), the non-directional dependent conserved quantity such as mass or density would be reconstructed to be of equal quantity. This might lead to possible undesired mass inflow. We circumvent this issue by a combination of Dirichlet and outflow boundary conditions: the non-direction related components are set to Dirichlet values, whereas the direction related components orthogonal to the edge are set to zero.

Combining in- and outflow for conserved quantities in edge space q^{L} of the cell at the boundary,
the conserved quantities q^{R} on the right side of the Riemann problem are respectively set to

With bounce back boundaries, we can simulate walls streaming back the flow with components
orthogonal to the domain boundary. Let the DoF on an edge in edge space be given by u_{i}^{L}. We can
then reconstruct the DoF on the ghost layer boundary with U^{R} := (u_{1}^{L},…,-u_{k}^{L},…,u_{n}^{L}) in edge space
and u_{k} the direction-dependent component perpendicular to the edge. Among others, this boundary
condition is used in our analytic benchmark (see Sec. 6.2).