Search

We present novel fully-symmetric quadrature rules with positive weights and strictly interior nodes of degrees up to 84 on triangles and 40 on tetrahedra. Initial guesses for solving the nonlinear systems of equations needed to derive quadrature rules are generated by forming tensor-product structures on quadrilateral/hexahedral subdomains of the simplices using the Legendre-Gauss nodes on the first half of the line reference element. In combination with a methodology for node elimination, these initial guesses lead to the development of highly efficient quadrature rules, even for very high polynomial degrees. Using existing estimates of the minimum number of quadrature points for a given degree, we show that the derived quadrature rules on triangles and tetrahedra are more than 95% and 80% efficient, respectively, for almost all degrees. The accuracy of the quadrature rules is demonstrated through numerical examples., Comment: 17 pages, 6 figures

We present an approach to construct efficient sparse summation-by-parts (SBP) operators on triangles and tetrahedra with a tensor-product structure. The operators are constructed by splitting the simplices into quadrilateral or hexahedral subdomains, mapping tensor-product SBP operators onto the subdomains, and assembling back using a continuous-Galerkin-type procedure. These tensor-product split-simplex operators do not have repeated degrees of freedom at the interior interfaces between the split subdomains. Furthermore, they satisfy the SBP property by construction, leading to stable discretizations. The accuracy and sparsity of the operators substantially enhance the efficiency of SBP discretizations on simplicial meshes. The sparsity is particularly important for entropy-stable discretizations based on two-point flux functions, as it reduces the number of two-point flux computations. We demonstrate through numerical experiments that the operators exhibit efficiency surpassing that of the existing dense multidimensional SBP operators by more than an order of magnitude in many cases. This superiority is evident in both accuracy per degree of freedom and computational time required to achieve a specified error threshold., Comment: 28 pages, 15 figures

High-order Hadamard-form entropy stable multidimensional summation-by-parts discretizations of the Euler and compressible Navier-Stokes equations are considerably more expensive than the standard divergence-form discretization. In search of a more efficient entropy stable scheme, we extend the entropy-split method for implementation on unstructured grids and investigate its properties. The main ingredients of the scheme are Harten's entropy functions, diagonal-$ \mathsf{E} $ summation-by-parts operators with diagonal norm matrix, and entropy conservative simultaneous approximation terms (SATs). We show that the scheme is high-order accurate and entropy conservative on periodic curvilinear unstructured grids for the Euler equations. An entropy stable matrix-type interface dissipation operator is constructed, which can be added to the SATs to obtain an entropy stable semi-discretization. Fully-discrete entropy conservation is achieved using a relaxation Runge-Kutta method. Entropy stable viscous SATs, applicable to both the Hadamard-form and entropy-split schemes, are developed for the compressible Navier-Stokes equations. In the absence of heat fluxes, the entropy-split scheme is entropy stable for the compressible Navier-Stokes equations. Local conservation in the vicinity of discontinuities is enforced using an entropy stable hybrid scheme. Several numerical problems involving both smooth and discontinuous solutions are investigated to support the theoretical results. Computational cost comparison studies suggest that the entropy-split scheme offers substantial efficiency benefits relative to Hadamard-form multidimensional SBP-SAT discretizations., Comment: 34 pages, 8 figures

We analyze the stability and functional superconvergence of discretizations of diffusion problems with the narrow-stencil second-derivative generalized summation-by-parts (SBP) operators coupled with simultaneous approximation terms (SATs). Provided that the primal and adjoint solutions are sufficiently smooth and the SBP-SAT discretization is primal and adjoint consistent, we show that linear functionals associated with the steady diffusion problem superconverge at a rate of $ 2p $ when a degree $ p+1 $ narrow-stencil or a degree $ p $ wide-stencil generalized SBP operator is used for the spatial discretization. Sufficient conditions for stability of adjoint consistent discretizations with the narrow-stencil generalized SBP operators are presented. The stability analysis assumes nullspace consistency of the second-derivative operator and the invertibility of the matrix approximating the first derivative at the element boundaries. The theoretical results are verified by numerical experiments with the one-dimensional Poisson problem., Comment: 21 pages, 8 figures

Several types of simultaneous approximation term (SAT) for diffusion problems discretized with diagonal-norm multidimensional summation-by-parts (SBP) operators are analyzed based on a common framework. Conditions under which the SBP-SAT discretizations are consistent, conservative, adjoint consistent, and energy stable are presented. For SATs leading to primal and adjoint consistent discretizations, the error in output functionals is shown to be of order $h^{2p}$ when a degree $p$ multidimensional SBP operator is used to discretize the spatial derivatives. SAT penalty coefficients corresponding to various discontinuous Galerkin fluxes developed for elliptic partial differential equations are identified. We demonstrate that the original method of Bassi and Rebay, the modified method of Bassi and Rebay, and the symmetric interior penalty method are equivalent when implemented with SBP diagonal-E operators that have diagonal norm matrix, e.g., the Legendre-Gauss-Lobatto SBP operator in one space dimension. Similarly, the local discontinuous Galerkin and the compact discontinuous Galerkin schemes are equivalent for this family of operators. The analysis remains valid on curvilinear grids if a degree $\le p+1$ bijective polynomial mapping from the reference to physical elements is used. Numerical experiments with the two-dimensional Poisson problem support the theoretical results., Comment: 44 pages, 11 figures