Steady-state laminar flow solutions for NACA 0012 airfoil

Abstract In this paper we consider the solution of the compressible Navier–Stokes equations for a class of laminar airfoil flows. The principal objective of this paper is to demonstrate that members of this class of laminar flows have steady-state solutions. These laminar airfoil flow cases are often used to evaluate accuracy, stability and convergence of numerical solution algorithms for the Navier–Stokes equations. In recent years such flows have also been used as test cases for high-order numerical schemes. While generally consistent steady-state solutions have been obtained for these flows using higher order schemes, a number of results have been published with various solutions, including unsteady ones. We demonstrate with two different numerical methods and a range of meshes with a maximum density that exceeds 8 × 10 6 grid points that steady-state solutions are obtained. Furthermore, numerical evidence is presented that even when solving the equations with an unsteady algorithm, one obtains steady-state solutions.

[1]  Eli Turkel,et al.  On Central-Difference and Upwind Schemes , 1992 .

[2]  Vít Dolejší,et al.  A design of residual error estimates for a high order BDF‐DGFE method applied to compressible flows , 2013 .

[3]  A. Jameson,et al.  A multigrid method for the Navier Stokes equations , 1986 .

[4]  R. C. Swanson,et al.  A Study of Multigrid Preconditioners Using Eigensystem Analysis , 2005 .

[5]  Eli Turkel,et al.  Artificial dissipation and central difference schemes for the Euler and Navier-Stokes equations , 1987 .

[6]  Eli Turkel,et al.  Convergence acceleration of Runge-Kutta schemes for solving the Navier-Stokes equations , 2007, J. Comput. Phys..

[7]  Dimitri J. Mavriplis,et al.  Multigrid solution of the Navier-Stokes equations on triangular meshes , 1989 .

[8]  A. Jameson.  Solution of the Euler equations for two dimensional transonic flow by a multigrid method , 1983 .

[9]  R. Swanson,et al.  Multistage Schemes With Multigrid for Euler and Navier-Stokes Equations , 1997 .

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

[11]  Arthur Rizzi,et al.  Computer-aided analysis of the convergence to steady state of discrete approximations to the euler equations , 1985 .

[12]  Stefan Langer.  Computer Aided Analysis of Preconditioned Multistage Runge-Kutta Methods Applied to Solve the Compressible Reynolds Averaged Navier-Stokes Equations , 2016 .

[13]  Eli Turkel.  Improving the accuracy of central difference schemes , 1989 .

[14]  Norbert Kroll,et al.  The DLR Flow Solver TAU - Status and Recent Algorithmic Developments , 2014 .

[15]  Eli Turkel,et al.  Analysis of a RK/Implicit Smoother for Multigrid , 2011 .

[16]  T. Gerhold,et al.  Calculation of Complex Three-Dimensional Configurations Employing the DLR-tau-Code , 1997 .

[17]  P. Laasonen.  ON THE DISCRETIZATION ERROR OF THE DIRICHLET PROBLEM IN A PLANE REGION WITH CORNERS. , 1967 .

[18]  V. Venkatakrishnan.  Newton solution of inviscid and viscous problems , 1988 .

[19]  Stefan Langer,et al.  Investigation and application of point implicit Runge–Kutta methods to inviscid flow problems , 2012 .

[20]  Antony Jameson,et al.  Multigrid solution of the Navier-Stokes equations for flow over wings , 1988 .

[21]  E. Turkel,et al.  A multistage time-stepping scheme for the Navier-Stokes equations , 1985 .

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

[23]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[24]  R. C. Swanson,et al.  An efficient solver for the RANS equations and a one-equation turbulence model , 2011 .

[25]  A computational study of vortex–airfoil interaction using high-order finite difference methods , 2010 .

[26]  Stefan Langer,et al.  Agglomeration multigrid methods with implicit Runge-Kutta smoothers applied to aerodynamic simulations on unstructured grids , 2014, J. Comput. Phys..

[27]  Norbert Kroll,et al.  The ADIGMA Project , 2010 .

[28]  A. Jameson,et al.  Lower-upper Symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations , 1988 .

[29]  Cord-Christian Rossow,et al.  Efficient computation of compressible and incompressible flows , 2007, J. Comput. Phys..

[30]  Bram van Leer,et al.  Design of Optimally Smoothing Multi-Stage Schemes for the Euler Equations , 1989 .

[31]  Cord-Christian Rossow,et al.  Convergence Acceleration for Solving the Compressible Navier-Stokes Equations , 2006 .

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

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

[34]  Implicit multigrid algorithm for the Navier-Stokes equations , 1991 .

[35]  M. D. Salas,et al.  Far-field boundary conditions for transonic lifting solutions to the Euler equations , 1986 .

[36]  D. Venditti,et al.  Anisotropic grid adaptation for functional outputs: application to two-dimensional viscous flows , 2003 .

[37]  Norbert Kroll.  ADIGMA - A European Initiative on the Development of Adaptive Higher-Order Variational Methods for Aerospace Applications : Results of a collaborative research project funded by the European Union, 2006-2009 , 2010 .

[38]  Jacques Periaux,et al.  Presentation of Problems and Discussion of Results , 1987 .

[39]  R. Abgrall,et al.  High order preserving residual distribution schemes for the laminar and turbulent Navier Stokes on arbitrary grids , 2013 .