Search NASASearch

Engineering topics

Anderson, W. K.

Publications and source records attributed to Anderson, W. K..

First-Order Model Management With Variable-Fidelity Physics Applied to Multi-Element Airfoil Optimization

First-order approximation and model management is a methodology for a systematic use of variable-fidelity models or approximations in optimization. The intent of model management is to attain convergence to high-fidelity solutions with minimal expense in high-fidelity computations. The savings in terms of computationally intensive evaluations depends on the ability of the available lower-fidelity model or a suite of models to predict the improvement trends for the high-fidelity problem, Variable-fidelity models can be represented by data-fitting approximations, variable-resolution models. variable-convergence models. or variable physical fidelity models. The present work considers the use of variable-fidelity physics models. We demonstrate the performance of model management on an aerodynamic optimization of a multi-element airfoil designed to operate in the transonic regime. Reynolds-averaged Navier-Stokes equations represent the high-fidelity model, while the Euler equations represent the low-fidelity model. An unstructured mesh-based analysis code FUN2D evaluates functions and sensitivity derivatives for both models. Model management for the present demonstration problem yields fivefold savings in terms of high-fidelity evaluations compared to optimization done with high-fidelity computations alone.

Alexandrov, N. M.

An O(Nm(sup 2)) Plane Solver for the Compressible Navier-Stokes Equations

A hierarchical multigrid algorithm for efficient steady solutions to the two-dimensional compressible Navier-Stokes equations is developed and demonstrated. The algorithm applies multigrid in two ways: 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 effect of boundary conditions in one direction are used to estimate the convergence rate of the algorithm for a model convection equation. Three alternating-line- implicit algorithms are compared in terms of efficiency. The analyses indicate that full multigrid efficiency 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.

Thomas, J. L.

The Use of Heavy Gas for Increased Reynolds Numbers in Transonic Wind Tunnels

The use of a high molecular weight test gas to increase the Reynolds number range of transonic wind tunnels is explored. Modifications to a small transonic wind tunnel are described and the real gas properties of the example heavy gas (sulfur hexafluoride) are discussed. Sulfur hexafluoride is shown to increase the test Reynolds number by a factor of more than 2 over air at the same Mach number. Experimental and computational pressure distributions on an advanced supercritical airfoil configuration at Mach 0.7 in both sulfur hexafluoride and nitrogen are presented. Transonic similarity theory is shown to be partially successful in transforming the heavy gas results to equivalent nitrogen (air) results, provided the correct definition of gamma is used.

Anders, J. B.

Navier-Stokes computations and experimental comparisons for multielement airfoil configurations

A two-dimensional unstructured Navier-Stokes code is utilized for computing the flow around multielement configurations. Comparisons are shown for a landing configuration with an advanced-technology flap as well as for a takeoff configuration. For each, a grid convergence study is conducted to assess inaccuracies caused by inadequate grid resolution. Although adequate resolution is obtained for determining the pressure distributions, further refinement is needed to sufficiently resolve the velocity profiles at high angles of attack. For the advanced flap configuration, comparisons of pressure distributions and lift are made with experimental data. Here, two flap riggings and two Reynolds numbers are considered. In general, the trends caused by variations in these quantities are well predicted by the computations, although the angle of attack for maximum lift is overpredicted. For the takeoff configuration, comparisons between pressures, velocity profiles, and Reynolds stresses are made at a single angle of attack with generally good agreement.

Anderson, W. K.

Grid convergence for adaptive methods

The ability of adaptive methods to obtain accurate results is examined using two different Euler solvers for a near-sonic flow containing several important flow features. It is shown that the accuracy obtained can be greatly affected by the lack of resolution of smooth portions of the flow caused by adapting only to the more prevalent flow features such as discontinuities. In particular, common methods of adaptation can lead to results in which shocks are well resolved but whose locations are highly inaccurate due to the lack of resolution of the smoother regions. An explanation for this behavior is given and a correction is proposed.

Warren, Gary P.

Navier-Stokes computations of vortical flows over low aspect ratio wings

An upwind-biased finite-volume algorithm is applied to the low-speed flow over a low aspect ratio delta wing from zero to forty degrees angle of attack. The differencing is second-order accurate spatially, and a multigrid algorithm is used to promote convergence to the steady state. The results compare well with the detailed experiments of Hummel (1983) and others for a Re(L) of 0.95 x 10 to the 6th. The predicted maximum lift coefficient of 1.10 at thirty-five degrees angle of attack agrees closely with the measured maximum lift of 1.06 at thirty-three degrees. At forty degrees angle of attack, a bubble type of vortex breakdown is evident in the computations, extending from 0.6 of the root chord to just downstream of the trailing edge.

Thomas, J. L.

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.

The effects of various implicit operators on a flux vector splitting method

Three different implicit operators in a numerical method for solving the two-dimensional steady Euler equations using flux vector splitting are investigated. These include the implementation of the true Jacobian matrices in the implicit part, the use of their splitting approximate form, and a modification of the implicit part of the scheme that uses the scalar form of the implicit part based on the spectral radii of the split Jacobian matrices. The three versions of the basic algorithm are compared with the results given by two common numerical methods using several two-dimensional test cases. Special attention is paid to the quality of results as well as computational efficiency and convergence properties.

Von Lavante, E.

Multigrid acceleration of the flux split Euler equations

Multigrid acceleration is applied to a flux-split algorithm for solving the Euler equations in two and three dimensions. The basic algorithm is an implicit spatially-split approximate factorization method. The stability of the scheme in comparison to other factorization is examined. Results are presented for two-dimensional airfoil flows and three-dimensional wing flows which demonstrate substantially improved convergence with the multigrid algorithm. An asymptotic spectral radius of 0.89 and 0.93 is attained for a 97 x 17 x 17 wing solution at subcritical and supercritical conditions, respectively.

Anderson, W. K.

A comparison of finite volume flux vector splittings for the Euler equations

A comparison is made between the computational results of the Steger-Warming (1981) and van Leer (1982) flux splitting methods, which have been applied in generalized coordinates to quasi-one-dimensional transonic flow in a nozzle and two-dimensional subsonic, transonic, and supersonic flow over airfoils. The latter splitting method leads to higher convergence rates and a sharper representation of shocks in the transition region. The second-order accurate, one-sided-difference model is extended to a third-order, upwind-biased model with only small additional computational effort.

Anderson, W. K.

A combined direct/inverse three-dimensional transonic wing design method for vector computers

A three-dimensional transonic-wing design algorithm for vector computers is developed, and the results of sample computations are presented graphically. The method incorporates the direct/inverse scheme of Carlson (1975), a Cartesian grid system with boundary conditions applied at a mean plane, and a potential-flow solver based on the conservative form of the full potential equation and using the ZEBRA II vectorizable solution algorithm of South et al. (1980). The accuracy and consistency of the method with regard to direct and inverse analysis and trailing-edge closure are verified in the test computations.

Weed, R. A.

Application of thin-layer Navier-Stokes equations near maximum lift

The flowfield about a NACA 0012 airfoil at a Mach number of 0.3 and Reynolds number of 1 million is computed through an angle of attack range, up to 18 deg, corresponding to conditions up to and beyond the maximum lift coefficient. Results obtained using the compressible thin-layer Navier-Stokes equations are presented as well as results from the compressible Euler equations with and without a viscous coupling procedure. The applicability of each code is assessed and many thin-layer Navier-Stokes benchmark solutions are obtained which can be used for comparison with other codes intended for use at high angles of attack. Reasonable agreement of the Navier-Stokes code with experiment and the viscous-inviscid interaction code is obtained at moderate angles of attack. An unsteady solution is obtained with the thin-layer Navier-Stokes code at the highest angle of attack considered. The maximum lift coefficient is overpredicted, however, in comparison to experimental data, which is attributed to the presence of a laminar separation bubble near the leading edge not modeled in the computations. Two comparisons with experimental data are also presented at a higher Mach number.

Anderson, W. K.