Discontinuous Galerkin discretisation with embedded boundary conditions
The purpose of this paper is to introduce discretisation methods of discontinuous Galerkin type for solving second order elliptic PDEs on a structured, regular rectangular grid, while the problem is defined on a curved boundary. The methods aim at high-order accuracy and the difficulty arises since the regular grid cannot follow the curved boundary. Starting with the Lagrange multiplier formulation for the boundary conditions, we derive variational forms for the discretisation of 2-D elliptic problems with embedded Dirichlet boundary conditions. Within the framework of structured, regular rectangular grids, we treat curved boundaries according to the principles that underlie the discontinuous Galerkin method. Thus, the high-order DG-discretisation is adapted in the cells with embedded boundaries. We give examples of approximation with tensor products of cubic polynomials. As an illustration, we solve a convection dominated boundary value problem on a complex domain. Although, of course, it is impossible to accurately represent a boundary layer with a complex structure by means of a cubic polynomial, the boundary condition treatment appears quite effective in handling such complex situations.