Effective programming of finite element methods for computational fluid dynamics on supercomputers
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The hierarchy of codes; boundary layers and turbulence models; uniqueness and artificial viscosity; efficiency of solution algorithms; data display and analysis; and cost of experimental validation are briefly examined.
Approximations applicable to the radiating, reacting, and conducting stagnation region of a hypervelocity vehicle were incorporated into a method for rapidly obtaining approximate solutions. This solution utilizes a coordinate system based upon the origin of the radiative losses and includes in a phenomenologically correct manner the effects of chemical and thermal nonequilibrium, and nonequilibrium, nongray radiative transfer. Results were presented which demonstrate the usefulness of the method and indicate which radiation parameters require further study and definition. Excellent comparisons were obtained with published results for the Fire2 data. An axisymmetric nonequilibrium inverse method was modified and extended and used to investigate and compare various vibration dissociation chemistry coupling models and radiative heat transfer approximations. The similarities, differences, and consequences of using these models in the Aero-assist Orbital Transfer Vehicles flight regime will be discussed.
Results from the Viscous Transonic Airfoil Workshop are compared with each other and with experimental data. Test cases used include attached and separated transonic flows for the NACA 0012 airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical method used vary widely and include: 16 Navier-Stokes methods, 2 Euler boundary layer methods, and 5 potential boundary layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.
The growing application of computational aerodynamics to nonlinear rotorcraft problems is outlined, with particular emphasis on the development of new methods based on the Euler and thin-layer Navier-Stokes equations. Rotor airfoil characteristics can now be calculated accurately over a wide range of transonic flow conditions. However, unsteady 3-D viscous codes remain in the research stage, and a numerical simulation of the complete flow field about a helicopter in forward flight is not now feasible. Nevertheless, impressive progress is being made in preparation for future supercomputers that will enable meaningful calculations to be made for arbitrary rotorcraft configurations.
An analytic method is described which uses static experimental data to predict the separated flow effect on rigid and elastic vehicle dynamics. Key parameters in the analytic relationship between steady and nonsteady aerodynamics are: the time lag occurring before a change of flow conditions can affect the separation-induced aerodynamic loads; the accelerated flow effect (i.e., the pressure gradient lag relative to the static aerodynamic characteristics); and the moving wall effect (i.e., the effect of the nonsteady boundary condition at the vehicle surface). Using the existing experimental data base, an analytic theory is formulated that can predict the separation-induced unsteady aerodynamics if the static characteristics are known from theory or experiment. Reference is made to increased-maneuverability advanced aircraft and to Space Shuttle Orbiter aerodynamics.
The mathematical properties of the pressure-velocity and vorticity-velocity formulations of the equations of viscous flow are compared. It is shown that a vorticity-velocity formulation exists which has the interesting property that noninertial effects only enter into the problem through the implementation of initial and boundary conditions. This valuable characteristic, along with other advantages of the vorticity-velocity approach, are discussed in detail.
Explore the source record for details and available documents.
The concept of a viscous shock-layer and several related versions of continuum theories/methods are examined for their adequacy as a viable framework to study flow physics and aerothermodynamics of relevance to sustained hypersonic flights. Considering the flat plate at angle of attack, or the wedge, as a generic example for the major aerodynamic component of a hypersonic vehicle, the relative importance of the molecular-transport effects behind the shock (in the form of the 'shock slip') and the wall-slip effects are studied. In the flow regime where the shock-transition-zone thickness remains small compared to the shock radius of curvature, a quasi-one-dimensional shock structure under the Burnett/thirteen-moment approximation, as well as particulate/collisional models, can be consistently developed. The fully viscous version of the shock-layer model is shown to provide the crucial boundary condition downstream the shock in this case. The gas-kinetic basis of the continuum description for the flow behind the bow shock, and certain features affecting the non-equilibrium flow chemistry, are also discussed.
The effect of periodic variations in the magnitude and direction of gravitational acceleration on buoyancy-driven convection in the melt during microgravity growth of low Prandtl number semiconductor materials in Brigdman-Stockbarger configuration is investigated through the numerical solution of the transient Navier-Stokes and energy equations. Results are presented over a wide range of Rayleigh numbers and frequency spectrum of periodic g-level variations corresponding to conditions present on earth and on board spacecrafts. The frequency response of convection in the melt during growth is found to be similar to that of a first order dynamic system with a cutoff frequency determined by the diffusion of momentum in the system. Results indicate that during growth of small diameter Germanium crystals, convection in the melt is insensitive to the measured g-gitter on board the shuttle.
The growing application of computational aerodynamics to nonlinear rotorcraft problems is outlined, with particular emphasis on the development of new methods based on the Euler and thin-layer Navier-Stokes equations. Rotor airfoil characteristics can now be calculated accurately over a wide range of transonic flow conditions. However, unsteady 3-D viscous codes remain in the research stage, and a numerical simulation of the complete flow field about a helicopter in forward flight is not now feasible. Nevertheless, impressive progress is being made in preparation for future supercomputers that will enable meaningful calculations to be made for arbitrary rotorcraft configurations.
The primary objective was the development of nonequilibrium radiation and chemistry models suitable for engineering applications associated with the flow fields about aeroassisted orbital transfer vehicles (AOTVs), the aero-assisted flight experiment vehicle (AFE), and other vehicles operating at superorbital velocities and very high attitudes.
Patched-grid algorithms for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver were investigated. Conservative and non-conservative approaches for performing zonal interpolations were implemented. The latter approach yields the most flexible technique in that it can handle both patched and overlaid grids. Results for a two-dimensional blunt body problem show that either approach yield accurate steady-state shock locations and jump conditions. In addition, calculations of the turbulent flow through a hypersonic inlet on a three-zone grid show that the numerical prediction is in good agreement with the experimental results. Through the use of a generalized coordinate transformation at the zonal interface between two or more blocks, the algorithm can be applied to highly stretched viscous grids and to arbitrarily-shaped zonal boundaries. Applications were made to the F-18 aircraft at subsonic, high-alpha conditions, in support of the NASA High-Alpha Research Program. The calculations were compared to ground-based and flight test experiments and were used as a guide to understanding the ground-based tests, which are laminar and transitional, and their relationship to flight. Calculations about a complete reconnaissance aircraft were also performed in order to further demonstrate the capability of the patched-grid algorithm.
The results of an investigation into the feasibility of using the MPP for direct and large eddy simulations of the Navier-Stokes equations is presented. A major part of this study was devoted to the implementation of two of the standard numerical algorithms for CFD. These implementations were not run on the Massively Parallel Processor (MPP) since the machine delivered to NASA Goddard does not have sufficient capacity. Instead, a detailed implementation plan was designed and from these were derived estimates of the time and space requirements of the algorithms on a suitably configured MPP. In addition, other issues related to the practical implementation of these algorithms on an MPP-like architecture were considered; namely, adaptive grid generation, zonal boundary conditions, the table lookup problem, and the software interface. Performance estimates show that the architectural components of the MPP, the Staging Memory and the Array Unit, appear to be well suited to the numerical algorithms of CFD. This combined with the prospect of building a faster and larger MMP-like machine holds the promise of achieving sustained gigaflop rates that are required for the numerical simulations in CFD.
Results from the Viscous Transonic Airfoil Workshop held in January 1987, are compared with each other and with experimental data. Test cases used include attached and separated transonic flows for the NACA 0012 airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical methods used vary widely and include: 16 Navier-Stokes methods, 2 Euler/boundary-layer methods, and 5 potential/boundary-layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The implementation and performance of a finite-difference algorithm for the compressible Navier-Stokes equations in two or three dimensions on the Connection Machine are described. This machine is a single-instruction multiple-data machine with up to 65536 physical processors. The implicit portion of the algorithm is of particular interest. Running times and megadrop rates are given for two- and three-dimensional problems. Included are comparisons with the standard codes on a Cray X-MP/48.