A finite element algorithm for computational fluid dynamics
Previously cited in issue 16, p. 2810, Accession no. A81-37558
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.
Previously cited in issue 16, p. 2810, Accession no. A81-37558
Previously cited in issue 16, p. 2809, Accession no. A81-37528
Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.
A three-dimensional parabolic Navier-Stokes code, PEPSIG, was used to analyze the flow in the subsonic diffuser section of a typical modern inlet design. The effect of curvature of the diffuser centerline and transitioning cross sections was studied to determine the primary cause of flow distortion in the duct. Total pressure values at the engine compressor face are reported.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Techniques for displaying two- and three-dimensional flowfield solutions are described. Several methods of illustrating flow structure are addressed including particle tracing, simulated oil flow, and shock finding. These are incorporated into an interactive graphics program for CFD flowfields, called PLOT3D. Emphasis is made on the difficulty in visualizing three-dimensional flow features, and the importance of color, fast 3D image manipulation, and dynamic movie play-back in displaying such flows. The need for advanced algorithms to identify shock waves, vortices, and separation lines is pointed out. It is likely that the supercomputer will be needed for this process because of the size of 3D and/or unsteady CFD databases.
A three-dimensional parabolic Navier-Stokes code, PEPSIG, was used to analyze the flow in the subsonic diffuser section of a typical modern inlet design. The effect of curvature of the diffuser centerline and transitioning cross sections was studied to determine the primary cause of flow distortion in the duct. Total pressure values at the engine compressor face are reported.
The field of helicopter rotor aerodynamics is basically the study of unsteady aerodynamic flows in a rotating and translating coordinate system. Current trends in this field are briefly reviewed by examining recent advances in lifting-surface theory, wake modeling, panel methods, and finite-difference models. Examples are used to illustrate selected current methods and some indications of promising future directions are highlighted.
The introduction of the supercomputer and recent advances in both Reynolds averaged, and large eddy simulation fluid flow approximation techniques to the Navier-Stokes equations, have created a robust environment for the exploration of problems of interest to the Navy in general, and the Naval Underwater Systems Center in particular. The nature of problems that are of interest, and the type of resources needed for their solution are addressed. The goal is to achieve a good engineering solution to the fluid-structure interaction problem. It is appropriate to indicate that a paper by D. Champman played a major role in developing the interest in the approach discussed.
At the present time the efforts on this project are organized into two areas: the development of an approximate stagnation point solution and approximate flowfield studies which can be used to develop and investigate shock jump, electron temperature, radiation, vibration-dissociation coupling, and chemistry models. Progress in each area is discussed.
The Jacobi matrix technique was applied to the direct calculation of inviscid supersonic flow about two dimensional airfoils of varying thickness, angle of attack and camber and axisymmetric bodies of varying thickness and taper; and the design (inverse) calculation of inviscid supersonic flow past airfoils described by a given family of pressure distributions and axisymmetric bodies described by a given family of pressure distributions. The method was also applied to subsonic potential flow about two dimensional airfoils by modifying Jameson's FLO36. Results are discussed.
An approach to analyzing CFD knowledge-based systems is proposed which is based, in part, on the concept of knowledge-level analysis. Consideration is given to the expert cooling fan design system, the PAN AIR knowledge system, grid adaptation, and expert zonal grid generation. These AI/CFD systems demonstrate that current AI technology can be successfully applied to well-formulated problems that are solved by means of classification or selection of preenumerated solutions.
A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented. The methods are demonstrated in fully vectorized codes for a CRAY type architecture. The emphasis will be on the Euler equations in generalized coordinates.
This paper highlights some recent accomplishments regarding CFD numerical algorithm constructions for generation of discrete approximate solutions to classes of Reynolds-averaged Navier-Stokes equations. Following an overview of turbulent closure modeling, and development of appropriate conservation law systems, a Taylor weak-statement semi-discrete approximate solution algorithm is developed. Various forms for completion to the final linear algebra statement are cited, as are a range of candidate numerical linear algebra solution procedures. This development sequence emphasizes the key building blocks of a CFD RNS algorithm, including solution trial and test spaces, integration procedure and added numerical stability mechanisms. A range of numerical results are discussed focusing on key topics guiding future research directions.
Explore the source record for details and available documents.