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 initial effort was concentrated on developing the quasi-analytical approach for two-dimensional transonic flow. To keep the problem computationally efficient and straightforward, only the two-dimensional flow was considered and the problem was modeled using the transonic small perturbation equation.
A new algorithm for an iterative computation of solutions of Laplace's or Poisson's equations in two dimensions, using Green's second identity, is presented. This algorithm converges strongly and geometrically and can be applied to curved, irregular, or moving boundaries with nonlinear and/or discontinuous boundary conditions. It has been implemented in Pascal on a number of micro- and minicomputers and applied to several geometries. Cases with known analytic solutions have been tested. Convergence to within 0.1 percent to 0.01 percent of the theoretical values are obtained in a few minutes on a microcomputer.
The acoustical near field of a counterrotating propeller configuration of unducted fan (UDF) type is presently predicted by means of a numerical methodology employing the transonic flow field generated by the three-dimensional numerical flow solver designated SSTAGE. This numerical approach simulates the time-averaged, three-dimensional inviscid flow field of the UDF by means of a four-stage Runge-Kutta integration scheme which marches the equations forward in time toward an asymptotic limit. The near field pressure time history is obtained directly from the flow field solution; by converting these histories with a Fourier series representation, specific data on the acoustic signal frequency spectra can be obtained.
The upper surface boundary layer on a transport wing model was extensively surveyed with miniature yaw probes at a subsonic and a transonic cruise condition. Additional data were obtained at a second transonic test condition, for which a separated region was present at mid-semispan, aft of mid-chord. Significant variation in flow direction with distance from the surface was observed near the trailing edge except at the wing root and tip. The data collected at the transonic cruise condition show boundary layer growth associated with shock wave/boundary layer interaction, followed by recovery of the boundary layer downstream of the shock. Measurements of fluctuating surface pressure and wingtip acceleration were also obtained. The influence of flow field unsteadiness on the boundary layer data is discussed. Comparisons among the data and predictions from a variety of computational methods are presented. The computed predictions are in reasonable agreement with the experimental data in the outboard regions where 3-D effects are moderate and adverse pressure gradients are mild. In the more highly loaded mid-span region near the trailing edge, displacement thickness growth was significantly underpredicted, except when unrealistically severe adverse pressure gradients associated with inviscid calculations were used to perform boundary layer calculations.
The purpose of this research effort was to begin the study of the application of hp-version finite elements to the numerical solution of optimal control problems. Under NAG-939, the hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element. One possible drawback is the increased computational effort within each element required in implementing hp-version finite elements. We are trying to determine whether this computational effort is sufficiently offset by the reduction in the number of time elements used and improved Newton-Raphson convergence so as to be useful in solving optimal control problems in real time. Because certain of the element interior unknowns can be eliminated at the element level by solving a small set of nonlinear algebraic equations in which the nodal values are taken as given, the scheme may turn out to be especially powerful in a parallel computing environment. A different processor could be assigned to each element. The number of processors, strictly speaking, is not required to be any larger than the number of sub-regions which are free of discontinuities of any kind.
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.
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.
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.