Discontinuous Galerkin and Petrov Galerkin methods for compressible viscous flows

Abstract In this paper, high-order finite-element discretizations consisting of discontinuous Galerkin (DG) and streamline upwind/Petrov Galerkin (SUPG) methods are investigated and developed for solutions of two- and three-dimensional compressible viscous flows. Both approaches treat the discretized system fully implicitly to obtain steady state solutions or to drive unsteady problems at each time step. A modified Spalart and Allmaras (SA) turbulence model is implemented and is discretized to an order of accuracy consistent with the main Reynolds Averaged Navier–Stokes (RANS) equations using the present high-order finite-element methods. To accurately represent the real geometry configurations for viscous flows, high-order curved boundary meshes are generated via a Computational Analysis PRogramming Interface (CAPRI), while the interior meshes are deformed subsequently through a linear elasticity solver. The mesh movement procedure effectively avoids the generation of invalid elements that can occur due to the projection of curved physical boundaries and thus allows high-aspect-ratio curved elements in viscous boundary layers. Several numerical examples, including large-eddy simulations of viscous flow over a three-dimensional circular cylinder and turbulent flows over a NACA 4412 airfoil and a high-lift multi-element airfoil at high angles of attack, are considered to show the capability of the present high-order finite-element schemes in capturing typical viscous effects such as flow separation and to compare the accuracy between the high-order DG and SUPG discretization methods.

[1]  William Gropp,et al.  Skjellum using mpi: portable parallel programming with the message-passing interface , 1994 .

[2]  Abrar H. Mohammad.  Large Eddy Simulation of Flow over a Cylinder Using High-Order Spectral Difference Method , 2010 .

[3]  T. Hughes,et al.  Stabilized finite element methods. I: Application to the advective-diffusive model , 1992 .

[4]  Ralf Hartmann,et al.  Error estimation and anisotropic mesh refinement for 3d laminar aerodynamic flow simulations , 2010, J. Comput. Phys..

[5]  Todd A. Oliver.  A High-Order, Adaptive, Discontinuous Galerkin Finite Element Method for the Reynolds-Averaged Navier-Stokes Equations , 2008 .

[6]  W. K. Anderson,et al.  Three-Dimensional Stabilized Finite Elements for Compressible Navier–Stokes , 2013 .

[7]  Li Wang,et al.  Adjoint-based h-p adaptive discontinuous Galerkin methods for the 2D compressible Euler equations , 2009, J. Comput. Phys..

[8]  W. K. Anderson,et al.  Petrov-Galerkin and discontinuous-Galerkin methods for time-domain and frequency-domain electromagnetic simulations , 2011, J. Comput. Phys..

[9]  David L. Darmofal,et al.  Shock Capturing with Higher-Order, PDE-Based Artificial Viscosity , 2007 .

[10]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[11]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[12]  Robert Haimes,et al.  Adjoint-Based Design Optimization Using CAD Parameterization Through CAPRI , 2012 .

[13]  Derek M. Causon,et al.  On the Choice of Wavespeeds for the HLLC Riemann Solver , 1997, SIAM J. Sci. Comput..

[14]  Dimitri J. Mavriplis,et al.  Adjoint Sensitivity Formulation for Discontinuous Galerkin Discretizations in Unsteady Inviscid Flow Problems , 2010 .

[15]  W. K. Anderson,et al.  High-Order Finite-Element Method for Three-Dimensional Turbulent Navier-Stokes , 2013 .

[16]  Dimitri J. Mavriplis,et al.  High-order discontinuous Galerkin methods using an hp-multigrid approach , 2006, J. Comput. Phys..

[17]  J. Peraire,et al.  Sub-Cell Shock Capturing for Discontinuous Galerkin Methods , 2006 .

[18]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[19]  W. K. Anderson,et al.  Shape sensitivity analysis for the compressible Navier–Stokes equations via discontinuous Galerkin methods , 2012 .

[20]  Venkat Venkatakrishnan,et al.  Higher Order Schemes for the Compressible Navier-Stokes Equations , 2003 .

[21]  Anthony Skjellum,et al.  Using MPI - portable parallel programming with the message-parsing interface , 1994 .

[22]  Dimitri J. Mavriplis,et al.  An hp-Adaptive Discontinuous Galerkin Solver for Aerodynamic flows on Mixed-Element Meshes , 2011 .

[23]  Timothy J. Barth,et al.  Numerical Methods for Gasdynamic Systems on Unstructured Meshes , 1997, Theory and Numerics for Conservation Laws.

[24]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: II. Beyond SUPG , 1986 .

[25]  Steven M. Klausmeyer,et al.  Comparative Results From a CFD Challenge Over a 2D Three-Element High-Lift Airfoil , 1997 .

[26]  David L. Darmofal,et al.  A triangular cut-cell adaptive method for high-order discretizations of the compressible Navier-Stokes equations , 2007, J. Comput. Phys..

[27]  S. Rebay,et al.  Numerical evaluation of two discontinuous Galerkin methods for the compressible Navier–Stokes equations , 2002 .

[28]  Dimitri J. Mavriplis,et al.  A Parallel hp-Multigrid Solver for Three-Dimensional Discontinuous Galerkin Discretizations of the Euler Equations , 2007 .

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

[30]  Arnold Afb,et al.  Comparison of SU/PG and DG Finite-Element Techniques for the Compressible Navier-Stokes Equations on Anisotropic Unstructured Meshes , 2013 .

[31]  Jaime Peraire,et al.  Navier-Stokes Solution Using Hybridizable Discontinuous Galerkin methods , 2011 .

[32]  D. Bonhaus,et al.  A higher order accurate finite element method for viscous compressible flows , 1999 .

[33]  Dimitri J. Mavriplis,et al.  Efficient Solution Techniques for Discontinuous Galerkin Discretizations of the Navier-Stokes Equations on Hybrid Anisotropic Meshes , 2010 .

[34]  Ralf Hartmann,et al.  Symmetric Interior Penalty DG Methods for the Compressible Navier-Stokes Equations II: Goal--Oriented A Posteriori Error Estimation , 2005 .

[35]  Miguel R. Visbal,et al.  An implicit discontinuous Galerkin method for the unsteady compressible Navier–Stokes equations , 2009 .

[36]  Dimitri J. Mavriplis,et al.  Multigrid algorithms for high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2009, J. Comput. Phys..

[37]  David L. Darmofal,et al.  p-Multigrid solution of high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2005 .

[38]  Li Wang,et al.  Implicit solution of the unsteady Euler equations for high-order accurate discontinuous Galerkin discretizations , 2006, J. Comput. Phys..

[39]  J. V. D. Vegt,et al.  Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics , 2002 .

[40]  David L. Darmofal,et al.  An Optimization Framework for Anisotropic Simplex Mesh Adaptation: Application to Aerodynamic Flows , 2012 .

[41]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[42]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[43]  D. Mavriplis,et al.  Robust Computation of Turbulent Flows Using a Discontinuous Galerkin Method , 2012 .

[44]  Zhi J. Wang,et al.  Adaptive high-order methods in computational fluid dynamics , 2011 .

[45]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[46]  J. Remacle,et al.  Shock detection and limiting with discontinuous Galerkin methods for hyperbolic conservation laws , 2004 .

[47]  D. Darmofal,et al.  Review of Output-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics , 2011 .

[48]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations , 1991 .

[49]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .