Transonic Navier-Stokes wing solution using a zonal approach. Part 1: Solution methodology and code validation
Abstract : A fast diagonalized Beam Warming algorithm is coupled with a zonal approach to solve the three dimensional Euler/Navier Stokes equations. The computer code, called Transonic Navier Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin layer Navier Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method coded. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000 point mesh for a NACA 0012 wing. A three order of magnitude drop in the L2 norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.
Status of computational fluid dynamics in the United States
CFD-related progress in U.S. aerospace industries and research institutions is evaluated with respect to methods employed, their applications, and the computer technologies employed in their implementation. Goals for subsonic CFD are primarily aimed at greater fuel efficiency; those of supersonic CFD involve the achievement of high sustained cruise efficiency. Transatmospheric/hypersonic vehicles are noted to have recently become important concerns for CFD efforts. Attention is given to aspects of discretization, Euler and Navier-Stokes general purpose codes, zonal equation methods, internal and external flows, and the impact of supercomputers and their networks in advancing the state-of-the-art.
Noniterative three-dimensional grid generation using parabolic partial differential equations
A new algorithm for generating three-dimensional grids has been developed and implemented which numerically solves a parabolic partial differential equation (PDE). The solution procedure marches outward in two coordinate directions, and requires inversion of a scalar tridiagonal system in the third. Source terms have been introduced to control the spacing and angle of grid lines near the grid boundaries, and to control the outer boundary point distribution. The method has been found to generate grids about 100 times faster than comparable grids generated via solution of elliptic PDEs, and produces smooth grids for finite-difference flow calculations.