Engineering topics
Publications and source records attributed to Kutler, Paul.
Charon is a software toolkit that enables engineers to develop high-performing message-passing programs in a convenient and piecemeal fashion. Emphasis is on rapid program development and prototyping. In this report a detailed description of the functional design of the toolkit is presented. It is illustrated by the stepwise parallelization of two representative code examples.
Schools can obtain scientific information over the information superhighway. However, information suppliers use formats that permit access and analysis by the "least common denominator" tools for access and analysis. The result: most sources of dynamic representations of science are in the format of flat movies. We can shorten the time to get the "scientific expeditions" into schools and provide a unifying focus to vendors and information suppliers by establishing a target and goals for the "least common denominator" for tools to be used to access and analyze information over the information superhighway.
One or the key objectives of the Applied Research Branch in the Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Allies Research Center is the accelerated introduction of highly parallel machines into a full operational environment. In this report we discuss the performance results obtained from the implementation of some computational fluid dynamics (CFD) applications on the Connection Machine CM-2 and the Intel iPSC/860. We summarize some of the experiences made so far with the parallel testbed machines at the NAS Applied Research Branch. Then we discuss the long term computational requirements for accomplishing some of the grand challenge problems in computational aerosciences. We argue that only massively parallel machines will be able to meet these grand challenge requirements, and we outline the computer science and algorithm research challenges ahead.
An overview of the computational effort to analyze forebody tangential slot blowing for use on the F/A-18 aircraft is presented. Tangential slot blowing generates side force and yawing moment which may be used to control the aircraft flying at high angle of attack. Computations using the isolated forebody are obtained at full-scale wind tunnel test conditions for direct comparison with available experimental data. The effects of jet exit conditions, jet length, and jet location are also studied using the isolated forebody. In addition, these computations are used to predict the effect of slot blowing at transonic maneuvering flight conditions. The effects of over- and under-blowing on force and moment production are analyzed. Non-time-accurate solutions are obtained to determine the steady-state side force and yawing moments generated by tangential slot blowing. Time-accurate solutions are obtained to study the force onset time lag of tangential slot blowing. The effect of blowing on the burst point location are then analyzed by obtaining computations using the aircraft geometry, which includes the wing, empennage, and faired-over inlets. The effect of blowing on the buffet loads on the vertical tails are analyzed using time-accurate computations. Comparison with available experimental data from full-scale wind tunnel and sub-scale wind tunnel tests are made. This computational analysis compliments the experimental results and provides a detailed understanding of the effects of tangential slot blowing on the flow field about the F/A-18.
A high speed phase locked interferometry system has been designed and developed for real-time measurements of the dynamic stall flow over a pitching airfoil. Point diffraction interferograms of incipient flow separation over a sinusoidally oscillating airfoil have been obtained at rates of up to 20 KHz and for free stream Mach numbers of 0.3 and 0.45. The images were recorded on ASA 125 and ASA 400 film using a drum camera. Special electronic timing and synchronizing circuits were developed to trigger the laser light source from the camera, and to initiate acquisition of the interferogram sequence from any desired phase angle of oscillation. The airfoil instantaneous angle of attack data provided by an optical encoder was recorded via a FIFO and in EPROM into a microcomputer. The interferograms have been analyzed using software developed in-house to get quantitative flow density and pressure distributions.
A hyperbolic grid generation method is used to generate prismatic grids and an approach using prismatic grids to solve the Euler equations is presented. The theory of the stability and feasibility of the hyperbolic grid generation method is presented. The hyperbolic grid generation method of Steger et al for structured grids is applied to a three dimensional triangularized surface definition to generate a grid that is unstructured on each successive layer. The grid, however, retains structure in the body-normal direction and has a computational cell shaped like a triangular prism. In order to take advantage of the structure in the normal direction, a finite-volume scheme that treats the unknowns along the normal direction implicitly is introduced and the flow over a sphere is simulated.
MIT flapping foil experiment was used as a validation case to evaluate the current incompressible Navier-Stokes approach with overlapped grid schemes. Steady-state calculations were carried out for overlapped and patched grids. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.
The existing pseudocompressibility method for the system of incompressible Navier-Stokes equations is extended to heat transfer problems by including the energy equation. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. Both forced and natural convection problems are examined. Numerical results from turbulent reattaching flow behind a backward-facing step will be compared against experimental measurements for the forced convection case. The validity of Boussinesq approximation to simplify the buoyancy force term will be investigated. The natural convective flow structure generated by heat transfer in a vertical rectangular cavity will be studied. The numerical results will be compared by experimental measurements by Morrison and Tran.
Computational fluid dynamics (CFD) is beginning to play a major role in the aircraft industry of the United States because of the realization that CFD can be a new and effective design tool and thus could provide a company with a competitive advantage. It is also playing a significant role in research institutions, both governmental and academic, as a tool for researching new fluid physics, as well as supplementing and complementing experimental testing. In this presentation, some of the progress made to date in CFD at NASA Ames will be reviewed. The presentation addresses the status of CFD in terms of methods, examples of CFD solutions, and computer technology. In addition, the role CFD will play in supporting the revolutionary goals set forth by the Aeronautical Policy Review Committee established by the Office of Science and Technology Policy is noted. The need for validated CFD tools is also briefly discussed.
Study of turbulent flows in rotating reference frames has proven to be one of the more challenging areas of turbulence research. The large number of theoretical, experimental, and computational studies performed over the years have demonstrated that the effect of solid-body rotation on turbulent flows is subtle and remains exceedingly difficult to predict. Because of the complexities associated with non-homogeneous turbulence, it is worthwhile to examine the effect of steady system rotation on the evolution of an initially isotropic turbulent flow. The assumption of statistical homogeneity considerably simplifies analysis and computation; calculation of homogeneous turbulence is further motivated since it possesses the essential physics found in more complex rotating flows. The principal objectives of the present study have therefore been to increase our fundamental understanding of turbulent flows in rotating reference frames through an examination of the asymptotic state of homogeneous rotating turbulence; particularly as to the existence of an asymptotic state which is self similar. Knowledge of an asymptotic similarity state permits prediction of the ultimate statistical evolution of the flow without requiring detailed knowledge of the complex, and not well understood, non-linear transfer processes. Aside from examination of possible similarity states in rotating turbulence, of further interest in this study has been an examination of the degree to which solid-body rotation induces a two-dimensional state in an initially isotropic flow.
One of the goals of NASA's High Alpha Technology Program is to provide flight-validated design methods for the high-angle-of-attack regime. This is an integrated effort utilizing computational simulations, wind tunnel experiments, and flight tests using the F-18 High Alpha Research Vehicle (HARV). The dominant physics of the aircraft flows in the high alpha regime changes as the angle of attack is increased. At moderate angle of attack the flow is characterized by boundary layer separation and the formation of tight vortices. As the angle of attack is increased, these vortices break down producing unsteady wakes. With further increase in angle of attack, the, vortex breakdown moves progressively upstream until the entire flowfield becomes dominated by the unsteady wake. Previous computational work has demonstrated the ability to simulate flows about the F-18 HARV in the medium-to-high angle of attack range, where the flowfield is characterized by the vortex formation and subsequent breakdown. This paper extends the previous computations to include conditions of 45 degree angle of attack where the flowfield becomes dominated by the unsteady wake shed from the Leading Edge Extension (LEX), and regions of laminar and transitional flow appear on the fuselage forebody. A more complete surface geometry is utilized, which includes the features of the engine nacelle, inlet diffuser, and the boundary layer diverter duct. A volume grid sensitivity study was also performed to extend the accuracy of the results, most notably in the prediction of the LEX vortex breakdown position. This paper includes comparisons of computational results with both in-flight surface pressure measurements, and flow visualizations of the surface and off-surface particle trajectories.
In recent years significant advances have been made for parallel computers in both hardware and software. Now parallel computers have become viable tools in computational mechanics. Many application codes developed on conventional computers have been modified to benefit from parallel computers. Significant speedups in some areas have been achieved by parallel computations. For single-discipline use of both fluid dynamics and structural dynamics, computations have been made on wing-body configurations using parallel computers. However, only a limited amount of work has been completed in combining these two disciplines for multidisciplinary applications. The prime reason is the increased level of complication associated with a multidisciplinary approach. In this work, procedures to compute aeroelasticity on parallel computers using direct coupling of fluid and structural equations will be investigated for wing-body configurations. The parallel computer selected for computations is an Intel iPSC/860 computer which is a distributed-memory, multiple-instruction, multiple data (MIMD) computer with 128 processors. In this study, the computational efficiency issues of parallel integration of both fluid and structural equations will be investigated in detail. The fluid and structural domains will be modeled using finite-difference and finite-element approaches, respectively. Results from the parallel computer will be compared with those from the conventional computers using a single processor. This study will provide an efficient computational tool for the aeroelastic analysis of wing-body structures on MIMD type parallel computers.
The performance of a workstation cluster used for the solution of the Reynolds-averaged Navier-Stokes equations is compared with a conventional vector supercomputer architecture. The application simulation of the steady flowfield about a transonic transport was computed using an implicit diagonal scheme in an overset mesh framework. Static load balancing was used, while coarse grain decomposition was achieved by solution of a grid zone per processor. Price/performance ratios are estimated for several scenarios in which such clusters may be utilized.
Fluid dynamics of turbomachines are complicated by inherently three dimensional structures such as endwall boundary layers, hub corner separation bubbles and tip-leakage flows. In addition, the relative motion between rotors and stators causes unsteady aerodynamic interactions to occur between blade rows. It is necessary to understand the aerodynamics associated with these interactions in order to design turbomachines that are both light and compact as well as reliable and efficient. An unsteady, three-dimensional, thin-layer, Navier-Stokes zonal algorithm is used to investigate the unsteady aerodynamics of multi-stage turbines and compressors. Relative motion between rotors and stators is made possible by the use of systems of patched and overlaid grids. Time-averaged pressures and pressure envelopes have been computed for several two- and three-dimensional single- and multi-stage configurations. Flow visualizations and computed results are in good agreement with experimental data.
The ability to control the extent of laminar flow on swept wings at supersonic speeds may be a critical element in developing the enabling technology for a High Speed Civil Transport (HSCT). Laminar boundary layers are less resistive to forward flight than their turbulent counterparts, thus the farther downstream that transition from laminar to turbulent flow in the wing boundary layer is extended can be of significant economic impact. Due to the complex processes involved experimental studies of boundary layer stability and transition are needed, and these are performed in "quiet" wind tunnels capable of simulating the low-disturbance environment of free flight. At Ames, a wind tunnel has been built to operate at flow conditions which match those of the HSCT laminar flow flight demonstration 'aircraft, the F-16XL, i.e. at a Mach number of 1.6 and a Reynolds number range of 1 to 3 million per foot. This will allow detailed studies of the attachment line and crossflow on the leading edge area of the highly swept wing. Also, use of suction as a means of control of transition due to crossflow and attachment line instabilities can be studied. Topics covered include: test operating conditions required; design requirements to efficiently make use of the existing infrastructure; development of an injector drive system using a small pilot facility; plenum chamber design; use of computational tools for tunnel and model design; and early operational results.
It is well known that slender bodies of revolution will develop an asymmetric, unsteady flow pattern in experimental tests, if the angle of incidence to the oncoming stream is above a critical value. It has been suggested that the origin of these asymmetric flows may stem from geometric imperfections of the model being tested, or from disturbances in the oncoming stream. In numerical simulations, it is possible to generate bodies of revolution which are perfectly symmetric about their longitudinal axis, and to impose uniform flow conditions which are free from disturbances. The current work presents numerical simulations of the flow about an ogive-cylinder configuration at 40 and 60 degree angle of incidence. These simulations. were performed using numerical algorithms which are also symmetric about the lateral plane of the cylinder body. The flowfields at 40 degree angle of attack were seen to remain symmetric to the round-off accuracy of the computer. At 60 degree angle of attack, a lateral force coefficient developed of O(1) which progressed to an alternate vortex shedding in time. The nature of this lateral force generation and vortex shedding was dependent on the choice of numerical algorithm. The origin of the asymmetries observed in the 60 degree angle of attack computations were traced to round-off errors in the implicit block-matrix inverter. A means of inverting the implicit operator matrices, which maintains the symmetry of the overall numerical algorithm was implemented.
Steady and unsteady flows for propulsion systems are efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard kappa - epsilon model equations. The resulting computer code is applied to the flow analysis inside an advanced rocket pump impeller in steadily rotating reference frames. Numerical results are compared with experimental measurements. The effects of exit and shroud cavities with the leak-age flow are investigated. Time-accurate incompressible Navier-Stokes formulation with the overlapped grid scheme capability was evaluated by using MIT flapping foil experiment. The grid dependency, turbulence model effects, and the effect of order of differencing were investigated. Numerical results were compared against experimental data. The resulting procedure were applied to unsteady flapping foil calculations. Two upstream NACA 0025 foils perform high frequency synchronized motion and generate unsteady flow conditions to the downstream larger stationary foil. Comparison between unsteady experimental data and numerical results from two different moving boundary procedures will be presented.
The dynamics of a class of TVD schemes for model hyperbolic and parabolic equations is studied numerically using a highly parallel supercomputer (CM-5). The objective is to utilize the highly parallel property of the CM-5 to reveal the reliable time step and entropy parameter ranges, and the degree of compressible flux limiters to avoid slow convergence and the production of nonphysical numerical solutions. We choose to study the nonlinear stability property of TVD schemes numerically since it is otherwise not amenable analytically.