Fully-Discrete Analysis of High-Order Spatial Discretizations with Optimal Explicit Runge–Kutta Methods
High-order unstructured methods have become a popular choice for the simulation of complex unsteady flows. Flux reconstruction (FR) is a high-order spatial discretization method, which has been found to be particularly accurate for scale-resolving simulations of complex phenomena. In addition, it ha...
Gespeichert in:
Veröffentlicht in: | Journal of scientific computing 2020-06, Vol.83 (3), p.63, Article 63 |
---|---|
Hauptverfasser: | , |
Format: | Artikel |
Sprache: | eng |
Schlagworte: | |
Online-Zugang: | Volltext |
Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
Zusammenfassung: | High-order unstructured methods have become a popular choice for the simulation of complex unsteady flows. Flux reconstruction (FR) is a high-order spatial discretization method, which has been found to be particularly accurate for scale-resolving simulations of complex phenomena. In addition, it has been shown to provide sufficient dissipation for implicit large-eddy simulation (ILES). In conjunction with an FR discretization, an appropriate temporal scheme must be chosen. A common choice is explicit schemes due to their efficiency and ease of implementation. However, these methods usually require a small time-step size to remain stable. Recently, the development of optimal explicit Runge–Kutta (OERK) schemes has enabled stable simulations with larger time-step sizes. Hence, we analyze the fully-discrete properties of the FR method with OERK temporal schemes. We show results for first, second, third, fourth and eighth-order OERK schemes. We observe that OERK schemes modify the spectral behaviour of the semidiscretization. In particular, dissipation decreases in the region of high wavenumbers. We observe that higher-order OERK schemes require a smaller time step than the low-order schemes. However, they follow the dispersion relations of the FR scheme for a larger range of wavenumbers. We validate our analysis with simple advection test cases. It was observed that first and second-degree temporal schemes introduce a relatively large amount of error in the solutions. A one-dimensional ILES test case showed that, as long as the time-step size is not in the vicinity of the stability limit, results are generally similar to classical RK schemes. |
---|---|
ISSN: | 0885-7474 1573-7691 |
DOI: | 10.1007/s10915-020-01243-8 |