W. Kyle Anderson (Texas A&M Class of ’81) was the professor for my last course for my PhD. He is a brilliant man; unfortunately, I didn’t do well in this course. But optimization in general was crucial for my PhD disssertation. 
W. Kyle Anderson and V. Venkatakrishnan
NASA Langley Research Center Contract NASI19480
January 1997
A continuous adjoint approach for obtaining sensitivity derivatives on unstructured grids is developed and analysed. The derivation of the costate equations is presented, and a secondorder accurate discretization method is described. The relationship between the continuous formulation and a discrete formulation is explored for inviscid, as well as for viscous flow. Several limitations in a strict adherence to the continuous approach are uncovered, and an approach that circumvents these difficulties is presented. The issue of grid sensitivities, which do not arise naturally in the continuous formulation, is investigated and is observed to be of importance when dealing with geometric singularities. A method is described for modifying inviscid and viscous meshes during the design cycle to accommodate changes in the surface shape. The accuracy of the sensitivity derivatives is established by comparing with finitedifference gradients and several design examples are presented. 

Liou MengSing
NASA Technical Memorandum 87329
May 1986
The present paper shows an efficient numerical procedure for solving a set of nonlinear partial differential equations, specifically the steady Euler equations. Solutions of the equations were obtained by Newton’s linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that first and secondorder derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new Iterates and the second derivatives Indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. Firstorder upwind discretization In conjunction with the StegerWarming fluxvector splitting is employed on the implicit operator and a diagonal dominant matrix is resulted. However the explicit operator is represented by first and secondorder upwind differencings, using both StegerWarming’s and van Leer’s splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one and twodimensional flows, we show detailed study as to the phlyctenae, accuracy, and convergence of the present method. 

P.T.Williams
Oak Ridge National Laboratory ORNL/TM12389
September 1993
As tile field of computational fluid dynamics (CFD) continues to mature, algorithms are required to exploit the most recent advances in approximation theory, numerical mathematics, computing architectures, and hardware. Meeting this requirement is particularly challenging in incompressible fluid mechanics, where primitivevariable CFD formulations that are robust, while also accurate and efficient in three dimensions, remain an elusive goal. This dissertation asserts that one key to accomplishing this goal is recognition of the dual role assumed by the pressure, i.e., a mechanism for instantaneously embracing conservation of mass and a force in tile mechanical balance law for conservation of momentum. Proving this assertion has motivated the development of a new, primitivevariable, incompressible, CFD algorithm called tile Continuity Constraint Method (CCM). The theoretical basis fbr the CCM consists ofa finiteelement spatial semidiscretization of a Galerkin weak statement, equalorder interpolation for all statevariables, a Oimplicit timeintegration scheme, and a quasiNewton iterative procedure extended by a Taylor Weak Statement (TWS) formulation for dispersion error control. Original contributions to algorithmic theory include: (a) formulation of the unsteady evolution of the divergence error, (b) investigation of the role of nonsmoothness in the discretized continuityconstraint function, (c) development of a uniformly H^{t} Galerkin weak statement for the Reynoldsaveraged NavierStokes pressure Poisson equation, (d) derivation of physically and numerically wellposed boundary conditions, and (e) investigation of sparse data structures and iterative methods for solving the matrix algebra statements generated by the algorithm. In contrast to tile general family of”pressurerelaxation” incompressible CFD algorithms, the CCM does not use the pressure as merely a mathematical device to constrain the velocity distribution to conserve mass. Rather, the mathematically smooth and physicallymotivated genuine pressure is an underlying replacement for tile nonsmooth continuityconstraint function to control inherent dispersiveerror mechanisms. Dominated by this dispersiveerror mode, the nonsmoothness of the discrete continuityconstraint function is proven to play a critical role in its ability to remove the divergence error in tile discrete velocity distribution. The genuine pressure is calculated by the diagnostic pressure Poisson equation, evaluated using the verified solenoidal velocity field. This new separation of tasks also produces a genuinely clear view of the totally distinct boundary conditions required for the continuityconstraint function and genuine pressure. A broad range of 3dimensional verification, benchmark, and validation lest problems, as computed by the code CFDI,.PHI3D, completes this dissertation. 
This is but a sample portion of the book released by one of the authors. To order the entirety of this excellent book, click here. 
Harvard Lomax, Thomas H. Pulliam and David W. Zingg
SpringerVerlag
2003
This book is intended as a textbook for a first course in computational fluid dynamics and will be of interest to researchers and practitioners as well. It emphasizes fundamental concepts in developing, analysing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. The linear convection and diffusion equations are used to illustrate concepts throughout. The chosen approach, in which the partial differential equations are reduced to ordinary differential equations, and finally to difference equations, gives the book its distinctiveness and provides a sound basis for a deep understanding of the fundamental concepts in computational fluid dynamics.
Webmaster’s Note: this is an outstanding book, gets to the core of the subject in a clear and concise manner. 

H. C. Yee, R. F. Warming,and A. Harten
NASA Technical Memorandum 84342
March 1983
We examine the application of a new implicit unconditionally stable highresolution TVD scheme to steadystate calculations. It is a member of a oneparameter family of explicit and implicit secondorder accurate schemes developed by Harten for the computation of weak solutions of onedimensional hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a fairly rapid convergence rate, but also generates a highly resolved approximation to the steadystate solution. A detailed implementation of the implicit scheme for the one and twodimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and twodimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme. 
Edwin Sereno Holdredge taught Mechanical Engineering at Texas A&M University for many years. His use of dimensionless analysis, illustrated in this work, became a hallmark of his teaching.
A native of Lenoir City, TN and graduate of the University of Tennessee at Knoxville, he had a dry sense of humour that frequently eluded many of his students. Concerning this report, he stated in class that “The Army told me they didn’t know anything about fluid flow around buildings.” He paused and then said, “They were right, they didn’t.”

Edwin S. Holdredge and Bob H. Reed, Texas Engineering Experiment Station
August 1956
This report describes the principal work performed for the period 15 June 195431 August 1956 by the Texas Engineering Experiment Station under contract between the Texas Agricultural and Mechanical College System as Contractor and the Department of the Army. The referenced contracts cover investigations to determine air flow pattern and infiltration characteristics of militarytype buildings by means of scale models thereof. The following items were studied:

Air flow patterns around model buildings at wind speeds ranging from approximately one to twenty mph.

Amount of infiltration of particulate substances such as smokes into model buildings at wind speeds ranging from one to twenty mph.

Ths minimum pressure required to prevent infiltration of particulate substances into the model buildings at wind speeds of approximately one to twenty mph and gusts up to approximately twenty mph.


Thomas H. Pulliam
NASA Ames Research Centre
1985
Implicit finite difference schemes for solving two dimensional and three dimensional Euler and NavierStokes equations will be addressed. The methods are demonstrated in fully vectorized codes for a CRAY type architecture. We shall concentrate on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated nontime accurate steady state schemes. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented. Examples for 2D inviscid and viscous calculations (e.g. airfoils with a deflected spoiler, circulation control airfoils and unsteady buffeting) and also 3D viscous flow are included. 

H.C. Yee
NASA Technical Memorandum 89464
May 1987
The development of numerical methods for hyperbolic conservation laws has been a rapidly growing area for the last ten years. Many of the fundamental concepts and stateoftheart developments can only be found in meeting proceedings or internal reports. This review paper attempts to give an overview and a unified formulation of a class of shockcapturing methods. Special emphasis will be on the construction of the basic nonlinear scalar secondorder schemes and the methods of extending these nonlinear scalar schemes to nonlinear systems via the exact Riemann solver, approximate Riemann solvers, and fluxvector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. The performance of some of these schemes is illustrated by numerical examples for one, two and threedimensional gasdynamics problems. 