High-order discrete ordinate transport in non-conforming 2D Cartesian meshes

We present in this paper a numerical scheme for solving the time-independent first-order form of the Boltzmann equation in non-conforming 2D Cartesian meshes. The flux solution technique used here is the discrete ordinate method and the spatial discretization is based on discontinuous finite elements. In order to have p-refinement capability, we have chosen a hierarchical polynomial basis based on Legendre polynomials. The h-refinement capability is also available and the element interface treatment has been simplified by the use of special functions decomposed over the mesh entities of an element. The comparison to a classical SN method using the Diamond Differencing scheme as spatial approximation confirms the good behaviour of the method. (authors)