Optimizing Computational Process of High-Order Taylor Discontinuous Galerkin Method for Solving the Euler Equations

Solving the Euler equations often requires expensive computations of complex, high-order time derivatives. Although Taylor Discontinuous Galerkin (TDG) schemes are renowned for their accuracy and stability, directly evaluating third-order tensor derivatives can significantly reduce computational efficiency, particularly for large-scale, intricate flow problems. To overcome this difficulty, this paper presents an optimized numerical procedure that combines Taylor series time integration with the Discontinuous Galerkin (DG) approach. By replacing cumbersome tensor derivatives with simpler time derivatives of the Jacobian matrix and finite difference method inside the element to calculate the high-order time derivative terms, the proposed method substantially decreases the computational cost while maintaining accuracy and stability. After verifying its fundamental feasibility in one-dimensional tests, the optimized TDG method is applied to a two-dimensional forward-facing step problem. In all numerical tests, the optimized TDG method clearly exhibits a computational efficiency advantage over the conventional TDG method, therefore saving a great amount of time, nearly 70%. This concept can be naturally extended to higher-dimensional scenarios, offering a promising and efficient tool for large-scale computational fluid dynamics simulations.