Discontinuous Galerkin Turbulent Flow Simulations of NASA Turbulence Model Validation Cases and High Lift Prediction Workshop Test Case DLR-F11

A high-order Discontinuous Galerkin (DG) solver is used to simulate three turbulent flow simulations. Two of these simulations come from the NASA turbulence modeling resource and include turbulent flow over a hemisphere-cylinder and turbulent flow over a 3D bump in a channel. The third simulation is the DLR-F11 from the second high-lift prediction workshop. The flow is simulated by solving the Reynolds-Averaged Navier-Stokes equations closed by the negative Spalart-Allmaras turbulence model. For these simulations, lift, drag, pitching moment, pressure, and skin friction coefficients are provided for multiple grids and discretization orders and compared against other simulation results from well known finite-volume solvers. The simulations give very similar results to these benchmark solvers, pointing towards fully mesh resolved simulations and providing verification evidence of correct and consistent implementation of these discretizations. Results obtained using the high-order DG discretizations show higher accuracy using fewer degrees of freedom compared to the finite-volume discretizations. Also, it is shown for the DG simulations that p-refinement converges quicker to the mesh resolved solutions compared to h-refinement. However, there are still issues with the efficiency of DG and the time to solution is still an order of magnitude slower compared to finite-volumes methods such as NSU3D.

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

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

[3]  Kazuyoshi Takayama,et al.  An artificially upstream flux vector splitting scheme for the Euler equations , 2003 .

[4]  C. Rumsey,et al.  Grid-Adapted FUN3D Computations for the Second High Lift Prediction Workshop , 2015 .

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

[6]  Boris Diskin,et al.  Grid Convergence for Turbulent Flows(Invited) , 2015 .

[7]  W. K. Anderson,et al.  An implicit upwind algorithm for computing turbulent flows on unstructured grids , 1994 .

[8]  Dimitri J. Mavriplis,et al.  High-Order Discontinuous Galerkin Mesh Resolved Turbulent Flow Simulations of a NACA 0012 Airfoil (Invited) , 2015 .

[9]  T. Hsieh,et al.  An Investigation of Separated Flow about a Hemisphere-Cylinder at 0- to 19-Deg Incidence in the Mach Number Range from 0.6 to 1.5 , 1977 .

[10]  Forrester T. Johnson,et al.  Modi cations and Clari cations for the Implementation of the Spalart-Allmaras Turbulence Model , 2011 .

[11]  W. K. Anderson,et al.  Finite-Element Solutions for Turbulent Flow over the NACA 0012 Airfoil (Invited) , 2015 .

[12]  R. Hartmann,et al.  Adaptive discontinuous Galerkin finite element methods for the compressible Euler equations , 2002 .

[13]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

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

[15]  Dimitri J. Mavriplis,et al.  3D Mixed Element Discontinuous Galerkin with Shock Capturing , 2013 .

[16]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[17]  Omar D. Lopez,et al.  Detached Eddy Simulation of the DLR-F11 wing/body Configuration as a Contribution to the 2nd AIAA CFD High Lift Prediction Workshop , 2014 .

[18]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[19]  K. Fidkowski,et al.  Constrained pseudo‐transient continuation , 2015 .

[20]  W. K. Anderson,et al.  High-Order Finite-Element Method and Dynamic Adaptation for Two-Dimensional Laminar and Turbulent Navier-Stokes , 2014 .

[21]  W. K. Anderson,et al.  An Adjoint-Based hp-Adaptive Petrov-Galerkin Method for Turbulent Flows , 2015 .

[22]  Neal T. Frink,et al.  Tetrahedral Unstructured Navier-Stokes Method for Turbulent Flows , 1998 .

[23]  David L. Darmofal,et al.  Application of a Higher-order Adaptive Method to RANS Test Cases (Invited) , 2015 .

[24]  Ralf Hartmann,et al.  An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier-Stokes equations , 2008, J. Comput. Phys..

[25]  David L. Darmofal,et al.  An adaptive simplex cut-cell method for discontinuous Galerkin discretizations of the Navier-Stokes equations , 2007 .

[26]  Krzysztof J. Fidkowski,et al.  High-Order Output-Based Adaptive Simulations of Turbulent Flow in Two Dimensions , 2016 .

[27]  Ralf Rudnik,et al.  EUROLIFT Test Case Description for the 2nd High Lift Prediction Workshop , 2012 .

[28]  Dimitri J. Mavriplis,et al.  NSU3D Results for the Second AIAA High-Lift Prediction Workshop , 2015 .