An implicit, Navier-Stokes solution algorithm is presented for the computation of turbulent flow on unstructured grids. The inviscid fluxes are computed using an upwind algorithm and the solution is advanced in time using a backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with a point-implicit relaxation scheme. This methodology provides a viable and robust algorithm for computing turbulent flows on unstructured meshes.Results are shown for subsonic flow over a NACA 0012 airfoil and for transonic flow over a RAE 2822 airfoil exhibiting a strong upper-surface shock. In addition, results are shown for 3-element and 4-element airfoil configurations. For the calculations, two one-equation turbulence models are utilized. For the NACA 0012 airfoil, a pressure distribution and force data are compared with other computational results as well as with experiment. Comparisons of computed pressure distributions and velocity profiles with experimental data are shown for the RAE airfoil and for the 3-element configuration. For the 4-element case, comparisons of surface pressure distributions with experiment are made. In general, the agreement between the computations and the experiment is good.

SummaryAn implicit code for computing inviscid and viscous incompressible flows on unstructured grids is described. The foundation of the code is a backward Euler time discretization for which the linear system is approximately solved at each time step with either a point implicit method or a preconditioned Generalized Minimal Residual (GMRES) technique. For the GMRES calculations, several techniques are investigated for forming the matrix-vector product. Convergence acceleration is achieved through a multigrid scheme that uses non-nested coarse grids that are generated using a technique described in the present paper. Convergence characteristics are investigated and results are compared with an exact solution for the inviscid flow over a four-element airfoil. Viscous results, which are compared with experimental data, include the turbulent flow over a NACA 4412 airfoil, a three-element airfoil for which Mach number effects are investigated, and three-dimensional flow over a wing with a partial-span flap.

An aerodynamic design algorithm for turbulent flows using unstructured grids is described. The current approach uses adjoint (costate) variables to obtain derivatives of the cost function. The solution of the adjoint equations is obtained by using an implicit formulation in which the turbulence model is fully coupled with the flow equations when solving for the costate variables. The accuracy of the derivatives is demonstrated by comparison with finite-difference gradients and a few sample computations are shown. Recommendations on directions of further research into the Navier-Stokes design process are made.

SummaryAn implicit code for computing inviscid and viscous incompressible flows on unstructured grids is described. The foundation of the code is a backward Euler time discretization for which the linear system is approximately solved at each time step with either a point implicit method or a preconditioned Generalized Minimal Residual (GMRES) technique. For the GMRES calculations, several techniques are investigated for forming the matrix-vector product. Convergence acceleration is achieved through a multigrid scheme that uses non-nested coarse grids that are generated using a technique described in the present paper. Convergence characteristics are investigated and results are compared with an exact solution for the inviscid flow over a four-element airfoil. Viscous results, which are compared with experimental data, include the turbulent flow over a NACA 4412 airfoil, a three-element airfoil for which Mach number effects are investigated, and three-dimensional flow over a wing with a partial-span flap.

Abstract

A hierarchical multigrid algorithm for e cient steady solutions to the two-dimensional compressible Navier-Stokes equations is developed and demonstrated. The algorithm applies multigrid in twoways: a Full Approximation Scheme FAS for a nonlinear residual equation and a Correction Scheme CS for a linearized defect correction implicit equation. Multigrid analyses which include the e ect of boundary conditions in one direction are used to estimate the convergence rate of the algorithm for a model convection equation. Three alternating-lineimplicit algorithms are compared in terms of eciency. The analyses indicate that full multigrid e ciency is not attained in the general case; the number of cycles to attain convergence is dependent on the mesh density for high-frequency cross-stream variations. However, the dependence is reasonably small and fast convergence is eventually attained for any given frequency with either the FAS or the CS scheme alone. The paper summarizes numerical computations for which convergence has been attained to within truncation error in a few multigrid cycles for both inviscid and viscous ow simulations on highly stretched meshes. Head, Aerodynamic and Acoustic Methods Branch AAMB, AIAA Fellow. y Member AAMB, AIAA Senior Member. z Member AAMB, AIAA Associate Fellow.

No abstract

Committee Chairman: Bernard Grossman Aerospace and Ocean Engineering (ABSTRACT)The Streamline Upwind/Petrov-Galerkin (SU/PG) method is applied to higher-order finite-element discretizations of the Euler equations in one dimension and the NavierStokes equations in two dimensions. The unknown flow quantities are discretized on meshes of triangular elements using triangular Bezier patches. The nonlinear residual equations are solved using an approximate Newton method with a pseudotime term. The resulting linear system is solved using the Generalized Minimum Residual algorithm with block diagonal preconditioning.The exact solutions of Ringleb flow and Couette flow are used to quantitatively establish the spatial convergence rate of each discretization. Examples of inviscid flows including subsonic flow past a parabolic bump on a wall and subsonic and transonic flows past a NACA 0012 airfoil and laminar flows including flow past a a flat plate and flow past a NACA 0012 airfoil are included to qualitatively evaluate the accuracy of the discretizations. The scheme achieves higher order accuracy without modification. Based on the test cases presented, significant improvement of the solution can be expected using the higherorder schemes with little or no increase in computational requirements. The nonlinear system also converges at a higher rate as the order of accuracy is increased for the same number of degrees of freedom; however, the linear system becomes more difficult to solve.Several avenues of future research based on the results of the study are identified, including improvement of the SU/PG formulation, development of more general grid generation strategies for higher order elements, the addition of a turbulence model to extend the method to high Reynolds number flows, and extension of the method to three-dimensional flows. An appendix is included in which the method is applied to inviscid flows in three iii dimensions. The three-dimensional results are preliminary but consistent with the findings based on the two-dimensional scheme.iv