Entropy Stable Summation-by-Parts Methods for Hyperbolic Conservation Laws on h/p Non-Conforming Meshes

In this work we present high-order primary conservative and entropy stable schemes for hyperbolic systems of conservation laws with geometric (h) and algebraic (p) non-conforming rectangular meshes. Throughout we rely on summation-by-parts (SBP) operators which discretely mimic the integration-by-parts rule to construct stable approximations. Thus, the discrete proofs of primary conservation and entropy stability can be done in a one-to-one fashion to the continuous analysis. Here, we consider different SBP operators based on finite difference as well as discontinuous Galerkin approaches. We derive non-conforming schemes by extending ideas of high-order primary conservative and entropy stable SBP methods on conforming meshes. Here, special attention is given to the coupling between non-conforming elements. The coupling is instructed to entropy stable projection operators. However, these projection operators suffer from a suboptimal degree. Therefore, we develop degree preserving SBP operators where the norm matrix has a higher degree compared to classical SBP operators. With these operators it is possible to construct entropy stable projection operators which have the same degree as the SBP differentiationmatrix. Typically, high-order primary conservative and entropy stable schemes are semi-discrete methods with a discretized spatial domain and assuming continuity in time. Therefore, temporal errors are introduced when integrating the conservation laws in time with standard methods, e.g. Runge-Kutta schemes, for which the entropy can have an unpredictable temporal behaviour. Thus, we extend high-order primary conservative and entropy stable semi-discrete methods to fully-discrete schemes on conforming and non-conforming meshes. This results in an implicit space-time method. We introduce a simple mesh generation strategy to obtain quadrilateral meshes surrounding two dimensional complex geometries. Finally, with the generated meshes we simulate a flow around a NACA0012 airfoil to demonstrate the benefits of considering non-conforming elements for a practical simulation

Encapsulated generalized summation-by-parts formulations for curvilinear and non-conforming meshes

We extend the construction of so-called encapsulated global summation-by-parts operators to the general case of a mesh which is not boundary conforming. Owing to this development, energy stable discretizations of nonlinear and variable coefficient initial boundary value problems can be formulated in simple and straightforward ways using high-order accurate operators of generalized summation-by-parts type. Encapsulated features on a single computational block or element may include polynomial bases, tensor products as well as curvilinear coordinate transformations. Moreover, through the use of inner product preserving interpolation or projection, the global summation-by-parts property in extended to arbitrary multi-block or multi-element meshes with non-conforming nodal interfaces