Abstract

Abstract

We present an analysis of the performance of some standard and optimized explicitly Runge­ Kutta schemes that are equipped with CFL-based and error-based time step adaptivity when they are coupled with the relaxation procedure to achieve fully-discrete entropy stability for complex compressible flow simulations. We investigate the performance of the temporal integration algorithms by simulating the flow past the NASA juncture flow model using the in-house KAUST SSDC hp-adaptive collocated entropy stable discontinuous Galerkin solver. In addition, we present a preliminary analysis of the performance of the SSDC framework on the Amazon web service cloud computing. The results indicate that SSDC scales well on the most recent and exotic computing architectures available on the Amazon cloud platform. Our findings might help select a more robust and efficient temporal integration algorithm and guide the choice of the EC2 AWS instances that give the best price and wall-clock-time performance to simulate industrially relevant turbulent flow problems.

Abstract

In the context of the numerical treatment of convective terms in compressible transport equations, general criteria for linear and quadratic invariants preservation, valid on uniform and non-uniform (Cartesian) meshes, have been recently derived by using a matrix-vector approach, for both finite-difference and finite-volume methods ([1, 2]). In this work, which constitutes a follow-up investigation of the analysis presented in [1, 2], this theory is applied to the spatial discretization of convective terms for the system of Euler equations. A classical formulation already presented in the literature is investigated and reformulated within the matrix-vector approach. The relations among the discrete versions of the various terms in the Euler equations are analyzed and the additional degrees of freedom identified by the proposed theory are investigated. Numerical simulations on a classical test case are used to validate the theory and to assess the effectiveness of the various formulations.

Abstract

The aim of this work is to model compressible flows involving shock waves past a solid obstacle using a non-conformal mesh. An Immersed Boundary Method (IBM) with feedback forcing and a volume penalization method are considered and compared. Both methods are validated on various test-cases. Accuracy and computational cost are discussed.

Abstract

We review and compare two techniques to get entropy stability for nodal Discontinuous Galerkin Spectral Element Methods (DG) in compressible flows. One technique is based on entropy split-forms, e.g., [8, 15, 5] and one is based on a direct algebraic correction [1]. We have implemented the Flux Differencing methodology for both, Legendre-Gauss-Lobatto (Lobatto) and Legendre-Gauss (Gauss) based spectral element basis functions. While the Lobatto operators belong to the class of diagonal norm summation-by-parts (SBP) operators, the Gauss operators belong to the generalized class of SBP operators, where it is not necessary that the boundary nodes are included. To reach entropy stability, respectively guaranteed entropy dissipation, a key ingredient is an entropy conserving numerical flux function. With this ingredient, only the volume integral term of the DG method has to be modified accordingly. We have also implemented an alternative technique, which is in general applicable for a wide range of discretizations. Abgrall [R. Abgrall, A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes. Journal of Computational Physics, vol 372, 2018] introduced an algebraic correction term that retains conservation of the primary quantities and is furthermore constructed such, that an entropy (in-) equality can be shown. The second technique is at first sight a simpler alternative to the split-form based approach. Hence, questions regarding its advantages and disadvantages naturally come up.