Observations Regarding Algorithms Required for Robust CFD Codes

Over the last three decades Computational Fluid Dynamics (CFD) has gradually joined the wind tunnel and flight test as a primary flow analysis tool for aerodynamic designers. CFD has had its most favorable impact on the aerodynamic design of the high-speed cruise configuration of a transport. This success has raised expectations among aerodynamicists that the applicability of CFD can be extended to the full flight envelope. However, the complex nature of the flows and geometries involved places substantially increased demands on the solution methodology and re- sources required. Currently most simulations involve Reynolds-Averaged Navier-Stokes (RANS) codes although Large Eddy Simulation (LES) and Detached Eddy Suimulation (DES) codes are occasionally used for component analysis or theoretical studies. Despite simplified underlying assumptions, current RANS turbulence models have been spectacularly successful for analyzing attached, transonic flows. Whether or not these same models are applicable to complex flows with smooth surface separation is an open question. A prerequisite for answering this question is ab- solute confidence that the CFD codes employed reliably solve the continuous equations involved. Too often, failure to agree with experiment is mistakenly ascribed to the turbulence model rather than inadequate numerics. Grid convergence in three dimensions is rarely achieved. Even residual convergence on a given grid is often inadequate. This paper discusses issues involved in residual and especially grid convergence.

[1]  D. P. Young,et al.  A locally refined rectangular grid finite element method: application to computational fluid dynamics and computational physics , 1990 .

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

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

[4]  Dimitri J. Mavriplis,et al.  Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes , 2003 .

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

[6]  Michael B. Bieterman,et al.  An adaptive grid method for analysis of 3D aircraft configurations , 1992 .

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

[8]  T. Hughes,et al.  MULTI-DIMENSIONAL UPWIND SCHEME WITH NO CROSSWIND DIFFUSION. , 1979 .

[9]  Michael B. Giles,et al.  Adjoint and defect error bounding and correction for functional estimates , 2003 .

[10]  L. Demkowicz One and two dimensional elliptic and Maxwell problems , 2006 .

[11]  James L. Thomas,et al.  Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio , 2013 .

[12]  N. Petrovskaya Discontinuous Weighted Least-Squares Approximation on Irregular Grids , 2008 .

[13]  Edward N. Tinoco,et al.  Summary of the Fourth AIAA CFD Drag Prediction Workshop , 2010 .

[14]  Edward N. Tinoco,et al.  Algorithm Issues and Challenges Associated with the Development of Robust CFD Codes , 2009 .

[15]  Forrester T. Johnson,et al.  THIRTY YEARS OF DEVELOPMENT AND APPLICATION OF CFD AT BOEING COMMERCIAL AIRPLANES, SEATTLE , 2005 .