Discontinuous Galerkin Finite Element Method with Anisotropic Local Grid Refinement for Inviscid Compressible Flows

A new discretization method for the three-dimensional Euler equations of gas dynamics is presented, which is based on the discontinuous Galerkin finite element method. Special attention is paid to an efficient implementation of the discontinuous Galerkin method that minimizes the number of flux calculations, which is generally the most expensive part of the algorithm. In addition a detailed discussion of the truncation error of the presented algorithm is given. The discretization of the Euler equations is combined with anisotropic grid refinement of an unstructured, hexahedron-type grid to achieve optimal resolution in areas with shocks, vortices, and other localized flow phenomena. The data structure and searching algorithms necessary for efficient calculation on highly irregular grids obtained with local grid refinement are discussed in detail. The method is demonstrated with calculations of the supersonic flow over a 10� ramp and the ONERA M6 wing under transsonic flow conditions.

[1]  Stanley Osher,et al.  Upwind schemes and boundary conditions with applications to Euler equations in general geometries , 1983 .

[2]  Jaap J. W. van der Vegt,et al.  Experiences with Advanced CFD Algorithms on NEC SX-4 , 1996, VECPAR.

[3]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[4]  Van Der Vegt.  Higher-order accurate Osher schemes with application to compressible boundary layer stability , 1993 .

[5]  J. Bonet,et al.  An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems , 1991 .

[6]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[7]  Directionally adaptive finite element method for multidimensional Euler and Navier-Stokes equations , 1993 .

[8]  A. Stroud.  Approximate calculation of multiple integrals , 1973 .

[9]  J va,et al.  Anisotropic grid refinement using an unstructured discontinuous Galerkin method for the three-dimensional Euler equations of gas dynamics , 1995 .

[10]  Nigel P. Weatherill,et al.  Adaptivity techniques for compressible inviscid flows , 1993 .

[11]  Kim S. Bey,et al.  A Runge-Kutta discontinuous finite element method for high speed flows , 1991 .

[12]  V. Venkatakrishnan.  Convergence to steady state solutions of the Euler equations on unstructured grids with limiters , 1995 .

[13]  San-Yih Lin,et al.  Discontinuous Galerkin Finite Element Method for Euler and Navier-Stokes Equations , 1993 .

[14]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[15]  M. V. Wilkes,et al.  The Art of Computer Programming, Volume 3, Sorting and Searching , 1974 .

[16]  M. J. Aftosmis,et al.  Upwind method for simulation of viscous flow on adaptively refined meshes , 1994 .

[17]  Donald E. Knuth,et al.  Sorting and Searching , 1973 .

[18]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[19]  Harold L. Atkins,et al.  QUADRATURE-FREE IMPLEMENTATION OF DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS , 1996 .

[20]  J. W. Boerstoel,et al.  Test Cases for Inviscid Flow Field Methods. , 1985 .

[21]  P. Raviart,et al.  On a Finite Element Method for Solving the Neutron Transport Equation , 1974 .

[22]  P. Lesaint.  Finite element methods for the transport equation , 1974 .

[23]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[24]  Russ D. Rausch,et al.  Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations , 1993 .