ALBERT

Description: Coupled Aerodynamic-Thermal-Structural (CATS) Analysis is a focused effort within the Numerical Propulsion System Simulation (NPSS) program to streamline multidisciplinary analysis of aeropropulsion components and assemblies. Multidisciplinary analysis of axial-flow compressor performance has been selected for the initial focus of this project. CATS will permit more accurate compressor system analysis by enabling users to include thermal and mechanical effects as an integral part of the aerodynamic analysis of the compressor primary flowpath. Thus, critical details, such as the variation of blade tip clearances and the deformation of the flowpath geometry, can be more accurately modeled and included in the aerodynamic analyses. The benefits of this coupled analysis capability are (1) performance and stall line predictions are improved by the inclusion of tip clearances and hot geometries, (2) design alternatives can be readily analyzed, and (3) higher fidelity analysis by researchers in various disciplines is possible. The goals for this project are a 10-percent improvement in stall margin predictions and a 2:1 speed-up in multidisciplinary analysis times. Working cooperatively with Pratt & Whitney, the Lewis CATS team defined the engineering processes and identified the software products necessary for streamlining these processes. The basic approach is to integrate the aerodynamic, thermal, and structural computational analyses by using data management and Non-Uniform Rational B-Splines (NURBS) based data mapping. Five software products have been defined for this task: (1) a primary flowpath data mapper, (2) a two-dimensional data mapper, (3) a database interface, (4) a blade structural pre- and post-processor, and (5) a computational fluid dynamics code for aerothermal analysis of the drum rotor. Thus far (1) a cooperative agreement has been established with Pratt & Whitney, (2) a Primary Flowpath Data Mapper has been prototyped and delivered to General Electric Aircraft Engines and Pratt & Whitney for evaluation, (3) a collaborative effort has been initiated with the National Institute of Standards and Testing to develop a Standard Data Access Interface, and (4) a blade tip clearance capability has been implemented into the Structural Airfoil Blade Engineering Routine (SABER) program. We plan to continue to develop the data mappers and data management tools. As progress is made, additional efforts will be made to apply these tools to propulsion system applications.

Description: Gaseous oxygen and hydrogen propellants used in a special engine energy cycle called Full-Flow Staged Combustion are believed to significantly increase the lifetime of a rocket engine's pumps. The cycle can also reduce the operating temperatures of the engine. Improving the lifetime of the hardware reduces its overall maintenance and operations costs, and is critical to reducing costs for the joint NASA/industry Reusable Launch Vehicle (RLV). The work in this project will demonstrate the performance and lifetime of one-element and many-element combustors with gaseous O2/H2 injectors. This work supporting the RLV program is a cooperative venture of the NASA Lewis Research Center, the NASA Marshall Space Flight Center, Rocketdyne, and the Pennsylvania State University. Information about gas-gas rocket injector performance with O2/H2 is very limited. Because of this paucity of data, new testing is needed to improve the knowledge base for testing and designing new injectors for the RLV and to improve computer models that predict the combusting gas flows of new injector designs. Therefore, detailed observations and measurements of the combusting flow from many-element injectors in a rocket engine are being sought. These observations and measurements will be done with three different tools: schlieren photography, ultraviolet imaging, and Raman spectroscopy. The schlieren system will take photos of the density differences in combusting flow, the ultraviolet movies will determine the location of the hydroxyl (OH) radical in the combustion flow, and the Raman spectroscopic measurements will provide the combustion temperature and amount of water (H2O), hydrogen (H2), and oxygen (O2) in the combustor. Marshall is providing overall program management, design and computational fluid dynamics (CFD) analyses, as well as funding for the work at Penn State. An existing, windowed combustor and several injectors will be provided by Rocketdyne--two injectors for the initial screening tests and one with an optimized design based on the best design found in the screening tests. Lewis will provide a nozzle and several injectors for the screening test program. The configuration of the injectors will be based on a design chosen by all the participants, and their elements will be based on the coaxial and impinging flow. Lewis also will provide the instrumentation for the flow-field measurements: schlieren, ultraviolet imaging, and Raman spectroscopy. In addition, thermocouples will measure heat flow on the injector face. Other traditional measurements of rocket performance will be made as well: chamber pressure, mass flow of each propellant, purge flow, and the barrier cooling gas flow. Penn State will conduct single-element testing with the injector elements from both the Rocketdyne and the jointly designed injectors. A wide variety of traditional and nontraditional injector designs will be tested in this program. The results will be valuable in computational fluid dynamics code validation and overall rocket combustion efficiency measurements. Correlations between combustion efficiency, laser measurements of species, and ultraviolet and visible light photography will also be made. Thus far, several different single-element injectors have been tested at Penn State and Lewis. The experimental setup of a rocket engine with a viewing window is shown. The combusting flow is shown. The results are helping engineers design the many element injectors.

Description: The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic- source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in at-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.

Description: Results are presented of an experimental program that investigated the use of a secondary air stream to control the amount of flow through a convergent-divergent nozzle. These static tests utilized high pressure, ambient temperature air that was injected at the throat of the nozzle through an annular slot. Multiple injection slot sizes and injection angles were tested. The introduction of secondary flow was made in an opposing direction to the primary flow and the resulting flow field caused the primary stream to react as though the physical throat size had been reduced. The percentage reduction in primary flow rate was generally about twice the injected flow rate. The most effective throttling was achieved by injecting through the smallest slot in an orientation most nearly opposed to the approaching primary flow. Thrust edliciency, as measured by changes in nozzle thrust coefficient, was highest at high nozzle pressure ratios, NPR. The static test results agreed with predictions obtained prior from PABSD, a fully viscous computational fluid dynamics program. Since use of such an injection system on gas turbine engine exhaust nozzles would be primarily at high NPRs, it was concluded that fluidic control holds promise for reducing nozzle weight and complexity on future systems.

Description: The component integration of a class of hypersonic high-lift configurations known as waveriders into hypersonic cruise vehicles was evaluated. A wind-tunnel model was developed which integrates realistic vehicle components with two waverider shapes, referred to as the straight-wing and cranked-wing shapes. Both shapes were conical-flow-derived waveriders for a design Mach number of 4.0. Experimental data and limited computational fluid dynamics (CFD) predictions were obtained over a Mach number range of 1.6 to 4.63 at a Reynolds number of 2.0 x 10(exp 6) per foot. The CFD predictions and flow visualization data confirmed the shock attachment characteristics of the baseline waverider shapes and illustrated the waverider flow-field properties. Experimental data showed that no significant performance degradations, in terms of maximum lift-to-drag ratios, occur at off-design Mach numbers for the waverider shapes and the integrated configurations. A comparison of the fully-integrated waverider vehicles to the baseline shapes showed that the performance was significantly degraded when all of the components were added to the waveriders, with the most significant degradation resulting from aftbody closure and the addition of control surfaces. Both fully-integrated configurations were longitudinally unstable over the Mach number range studied with the selected center of gravity location and for unpowered conditions. The cranked-wing configuration provided better lateral-directional stability characteristics than the straight-wing configuration.

Description: A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flow fields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three-dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a non-catalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the non-catalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated in both the radial and streamwise directions. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three-dimensional case in comparison to an all thin-layer Navier-Stokes solution.

Description: The Advanced Subsonic Combustion Rig (ASCR), a unique, state-of-the-art facility for conducting combustion research, is located at the NASA Lewis Research Center in Cleveland, Ohio. The ASCR, which was nearing completion at the close of 1995, will be capable of simulating the very high pressure and high temperature conditions that are expected to exist in future, advanced subsonic gas turbine (jet) engines. Future environmental regulations will require much cleaner burning (more environmentally friendly) aircraft engines. The ASCR is critical to the development of these cleaner engines. It will allow NASA and U.S. aircraft engine industry researchers to identify and test promising clean-burning gas turbine engine combustion concepts under the pressure and temperature conditions that are expected for those future engines. Combustion processes will be investigated for a variety of next-generation aircraft engine sizes, including engines for large, long-range aircraft (with typical trip lengths of about 3000 mi) and for regional aircraft (with typical trip lengths of about 400 mi). The ASCR design was conceived and initiated in 1993, and fabrication and construction of the rig, including the buildup of an advanced control room, took place throughout 1994 and 1995. In early 1996, the ASCR will be operational for obtaining research data. The ASCR is an intricate part of the NASA Advanced Subsonic Technology Propulsion Program, which is aimed at developing technologies critical to the next generation of gas turbine engines. This effort is in collaboration with the U.S. aircraft gas turbine engine industry. A goal of the Advanced Subsonic Technology Propulsion Program is to develop combustion concepts and technologies that will result in gas turbine engines that produce 50 percent less nitrous oxide (NO_x) pollutants than current engines do. This facility is unique in its capability to simulate advanced subsonic engine pressure, temperature, and air flow rate conditions. Specifically, it will provide operating temperatures up to 3000 F and pressures up to 60 atm. Under these conditions, researchers will obtain detailed combustion temperatures, pressures, and flow velocities as well as the chemical compositions of the combustion exhaust. Researchers also will be able to obtain data by using nonintrusive laser diagnostic techniques. The ASCR facility will be used to test fundamental combustion configurations (flametubes) for detailed study of combustion processes, to test sectors of gas turbine combustors to study the process in configurations more like those of aircraft engines, and in some cases to test full annular combustors.

Description: The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Description: The traditional approach to the design of engine inlets for commercial transport aircraft is a tedious process that ends with a less-than-optimum design. With the advent of high-speed computers and the availability of more accurate and reliable computational fluid dynamics (CFD) solvers, numerical optimization processes can effectively be used to design an aerodynamic inlet lip that enhances engine performance. The designers' experience at Boeing Corporation showed that for a peak Mach number on the inlet surface beyond some upper limit, the performance of the engine degrades excessively. Thus, our objective was to optimize efficiency (minimize the peak Mach number) at maximum cruise without compromising performance at other operating conditions. Using a CFD code NPARC, the NASA Lewis Research Center, in collaboration with Boeing, developed an integrated procedure at Lewis to find the optimum shape of a subsonic inlet lip and a numerical optimization code, ADS. We used a GRAPE-based three-dimensional grid generator to help automate the optimization procedure. The inlet lip shape at the crown and the keel was described as a superellipse, and the superellipse exponents and radii ratios were considered as design variables. Three operating conditions: cruise, takeoff, and rolling takeoff, were considered in this study. Three-dimensional Euler computations were carried out to obtain the flow field. At the initial design, the peak Mach numbers for maximum cruise, takeoff, and rolling takeoff conditions were 0.88, 1.772, and 1.61, respectively. The acceptable upper limits on the takeoff and rolling takeoff Mach numbers were 1.55 and 1.45. Since the initial design provided by Boeing was found to be optimum with respect to the maximum cruise condition, the sum of the peak Mach numbers at takeoff and rolling takeoff were minimized in the current study while the maximum cruise Mach number was constrained to be close to that at the existing design. With this objective, the optimum design satisfied the upper limits at takeoff and rolling takeoff while retaining the desirable cruise performance. Further studies are being conducted to include static and cross-wind operating conditions in the design optimization procedure. This work was carried out in collaboration with Dr. E.S. Reddy of NYMA, Inc.

Description: The Reynolds stress budget in a full developed turbulent channel flow for three Reynolds numbers (Re = 180,395,590) are used to investigate the near wall scaling of various turbulence quantities. We find that as the Reynolds number increases, the extent of the region where the production of the kinetic energy is equal to the dissipation increases. At the highest Reynolds number the region of equilibrium extends from y+ - 120 to y+ = 240. As the Reynolds number increases, we find that wall scaling collapses the budgets for the streamwise fluctuating component, but the budgets for the other two components show Reynolds number dependency.

Description: Features of chemically reacting separated hypersonic flows are identified and issues concerning their analysis and simulation are discussed. Emphasis is placed on flows of high temperature dissociating and ionizing air and current methods for studying and characterizing these flows, including separation, are reviewed. The aeroassist orbital transfer vehicle and its flight trajectory are used for illustration. Thermochemical nonequilibrium phenomena are emphasized and extension of continuum analysis to the high altitude slip-flow regime is considered.

Description: A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Description: The talk will discuss current work in developing and applying the INS2D and INS3D flow solvers to problems in high-lift aerodynamics, primarily multi-element airfoils. These flow solvers solve the incompressible Navier-Stokes equations using the method of artificial compressibility. High-lift system flowfields are perhaps the most challenging aeronautical configurations for CFD. Difficulties with convergence for multi-element configurations using fine-grid overset meshes has led to research in using Krylov-space iterative matrix solvers. In particular, a Generalized Minimum Residual (GMRES) solver has found to be a dramatic improvement over a previously used Gauss-Seidal line-relaxation solver.

Description: The status of the three-dimensional parabolized Navier-Stokes solver UPS is described. The UPS code, initiated at NASA Ames Research Center in 1986, continues to develop and evolve through application to supersonic and hypersonic flow fields. Hypersonic applications have motivated enhancement of the physical modeling capabilities of the code, specifically real gas modeling, boundary conditions, and turbulence and transition modeling. The UPS code has also been modified to enhance robustness and efficiency in order to be practically used in concert with an optimization code for supersonic transport design. These developments are briefly described along with some relevant results for generic test problems obtained during verification of the enhancements. Included developments and results have previously been published and widely disseminated domestically.

Description: The development of a two-dimensional viscous incompressible flow generated by an off center thin oscillating bd on top of a cavity is studied computationally as a prototype of vortex generators. The lid is placed asymmetrically over the cavity so that the gap size is different on either side of the cavity. An adaptive numerical scheme, based on high resolution viscous vortex methods, is used to integrate the vorticity/velocity formulation of the Navier-Stokes equations with the no-slip boun.lary condition enforced on the lid and cavity walls. Depending on the a amplitude and frequency of the oscillation as well as the the gap size, vorticity is ejected in the fluid above the cavity either from the large and/or the small gap. The results of the computations complement ongoing experimental work.

Description: Numerous two-and three-dimensional computational simulations were performed for the inlet associated with the combustor model for the hypersonic propulsion experiment in the NASA Ames 16-Inch Shock Tunnel. The inlet was designed to produce a combustor-inlet flow that is nearly two-dimensional and of sufficient mass flow rate for large scale combustor testing. The three-dimensional simulations demonstrated that the inlet design met all the design objectives and that the inlet produced a very nearly two-dimensional combustor inflow profile. Numerous two-dimensional simulations were performed with various levels of approximations such as in the choice of chemical and physical models, as well as numerical approximations. Parametric studies were conducted to better understand and to characterize the inlet flow. Results from the two-and three-dimensional simulations were used to predict the mass flux entering the combustor and a mass flux correlation as a function of facility stagnation pressure was developed. Surface heat flux and pressure measurements were compared with the computed results and good agreement was found. The computational simulations helped determine the inlet low characteristics in the high enthalpy environment, the important parameters that affect the combustor-inlet flow, and the sensitivity of the inlet flow to various modeling assumptions.

Description: The Burnett equations have been shown to potentially violate the second law of thermodynamics. The objective of this investigation is to correlate the numerical problems experienced by the Burnett equations to the negative production of entropy. The equations have had a long history of numerical instability to small wavelength disturbances. Recently, Zhong corrected the instability problem and made solutions attainable for one dimensional shock waves and hypersonic blunt bodies. Difficulties still exist when attempting to solve hypersonic flat plate boundary layers and blunt body wake flows, however. Numerical experiments will include one-dimensional shock waves, quasi-one dimensional nozzles, and expanding Prandlt-Meyer flows and specifically examine the entropy production for these cases.

Description: Computations have been performed on a three-dimensional inlet associated with the NASA Ames combustor model for the hypersonic propulsion experiment in the 16-inch shock tunnel. The 3-dimensional inlet was designed to have the combustor inlet flow nearly two-dimensional and of sufficient mass flow necessary for combustion. The 16-inch shock tunnel experiment is a short duration test with test time of the order of milliseconds. The flow through the inlet is in chemical non-equilibrium. Two test entries have been completed and limited experimental results for the inlet region of the combustor-model are available. A number of CFD simulations, with various levels of simplifications such as 2-D simulations, 3-D simulations with and without chemical reactions, simulations with and without turbulent conditions, etc., have been performed. These simulations have helped determine the model inlet flow characteristics and the important factors that affect the combustor inlet flow and the sensitivity of the flow field to these simplifications. In the proposed paper, CFD modeling of the hypersonic inlet, results from the simulations and comparison with available experimental results will be presented.

Description: A dynamic model is developed to describe the evolution of a vapor bubble growing at a nucleation site on a superheated surface under arbitrary gravity. The bubble is separated from the surface by a thin microlayer and grows due to the evaporation from the microlayer interface. The average thickness of the microlayer increases as the bubble expands along the surface if the evaporation rate is lower than some critical value. The corresponding threshold value of the surface temperature has to be associated with the burn-out crisis. Two main reasons make for bubble separation, which are the buoyancy force and a force caused by the vapor momentum that comes to the bubble with vapor molecules. The latter force is somewhat diminished if condensation takes place at the upper bubble surface in subcooled liquids. The action of the said forces is opposed by inertia of the additional mass of liquid as the bubble center rises above the surface and by inertia of liquid being expelled by the growing bubble in radial directions. An extra pressure force arises due to the liquid inflow into the microlayer with a finite velocity. The last force helps in holding the bubble close to the surface during an initial stage of bubble evolution. Two limiting regimes with distinctly different properties can be singled out, depending on which of the forces that favor bubble detachment dominates. Under conditions of moderately reduced gravity, the situation is much the same as in normal gravity, although the bubble detachment volume increases as gravity diminishes. In microgravity, the buoyancy force is negligible. Then the bubble is capable of staying near the surface for a long time, with intensive evaporation from the microlayer. It suggests a drastic change in the physical mechanism of heat removal as gravity falls below a certain sufficiently low level. Inferences of the model and conclusions pertaining to effects caused on heat transfer processes by changes in bubble hydrodynamics induced by gravity are discussed in connection with experimental evidence, both available in current and in as yet unpublished literature.

Description: When liquid is expelled by a vapor bubble growing at a nucleation site on a superheated surface, a thin microlayer underneath the bubble is left behind. It is evaporated from the free microlayer surface that provides for bubble growth. The average thickness of the microlayer determining the evaporation rate increases with time if the latter does not exceed a threshold value associated with the burn-out crisis. The bubble is described as a spherical segment with its flattened part adjoining the microlayer. This introduces two independent variables - the radius of the spherical part of the bubble surface and the polar angle that defines the relative area of the flattened part. They are to be found out from a set of two strongly nonlinear equations resulting from mass and momentum conservation laws. The first one depends on both microlayer thickness and nonmonotonously changing bubble base area. The second involves two major factors favoring bubble detachment - the buoyancy and a force due to the initial momentum of vapor input into the bubble. The former force depends on gravity whereas the latter one does not. It is why the limiting regimes of bubble evolution that correspond to normal or moderately reduced gravity and to microgravity feature drastically different properties. In the first case, the buoyancy dominates and the bubble evolves in such a manner as to become a full sphere at a moment that can be viewed as that of detachment. The detachment volume grows as gravity decreases. In the second case, the buoyancy is negligible and the bubble stays near the surface, while its volume continues to increase for a sufficiently long time. The findings are discussed in connection with experimental data obtained under different gravity conditions, some unpublished experiments being included. They help to understand why the pool boiling heat transfer coefficient frequently increases as gravity falls down and eventually vanishes.

Description: The fractional step and the pseudocompressibility methods for the solution of the incompressible Navier-Stokes equations are outlined. The fractional step method is based on finite-volume formulation and uses the pressure and the volume fluxes across the faces of each cell as dependent variables. The momentum equations are solved implicitly and the Poisson equation for the pressure is solved by using the multigrid method. The pseudocompressibility approach uses an implicit-higher-order-upwind differencing scheme for the convective terms together with the Gauss-Seidel line relaxation method. The dependent variables in the pseudocompressibility approach are the pressure and the cartesian velocity components in unstaggered mesh orientation. The 90-degree square duct flow, the wing-tip vortex wake flow and unsteady turbulent flows over an oscillating NACA 0015 airfoil are computed using both the fractional step and the pseudocompressibility methods. The results obtained from two different schemes are compared against experimental measurements.

Description: The accuracy and efficiency of two types of subiterations in both explicit and implicit Navier-Stokes codes are explored for unsteady laminar circular-cylinder flow and unsteady turbulent flow over an 18-percent-thick circular-arc (biconvex) airfoil. Grid and time-step studies are used to assess the numerical accuracy of the methods. Nonsubiterative time-stepping schemes and schemes with physical time subiterations are subject to time-step limitations in practice that are removed by pseudo time sub-iterations. Computations for the circular-arc airfoil indicate that a one-equation turbulence model predicts the unsteady separated flow better than an algebraic turbulence model; also, the hysteresis with Mach number of the self-excited unsteadiness due to shock and boundary-layer separation is well predicted.

Description: Some new guidelines on the usage of implicit linear multistep methods (LMMs) as time-dependent approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) are explored. The commonly used implicit LMMs in CFD belong to the class of superstable time discretizations. It can be shown that the nonlinear asymptotic behavior in terms of bifurcation diagrams and basins of attractions of these schemes can provide an improved range of initial data and time step over the linearized stability limit.

Description: The current study computationally examines one of the principle three-dimensional features of the flow over a high-lift system, the flow associated with a flap edge. Structured, overset grids were used in conjunction with an incompressible Navier-Stokes solver to compute the flow over a two-element high-lift configuration. The computations were run in a fully turbulent mode using the one-equation Baldwin-Barth model. Specific interest was given to the details of the flow in the vicinity of the flap edge, so the geometry was simplified to isolate this region. The geometry consisted of an unswept wing, which spanned a wind tunnel test section, equipped with a single element flap. Two flap configurations were computed; a full-span and a half-span Fowler flap. The chord based Reynolds number was 3.7 million for all cases. The results for the full-span flap agreed with two-dimensional experimental results and verified the method. Grid topologies and related issues for the half-span flap geometry are discussed. Results of the half-span flap case are presented with emphasis on the flow features associated with the flap edge.

Description: This lecture will discuss the computational tools currently available for high-lift multi-element airfoil analysis. It will present an overview of a number of different numerical approaches, their current capabilities, short-comings, and computational costs. The lecture will be limited to viscous methods, including inviscid/boundary layer coupling methods, and incompressible and compressible Reynolds-averaged Navier-Stokes methods. Both structured and unstructured grid generation approaches will be presented. Two different structured grid procedures are outlined, one which uses multi-block patched grids, the other uses overset chimera grids. Turbulence and transition modeling will be discussed.

Description: An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Description: Scientific visualization serves the dual purpose of exploration and exposition of the results of numerical simulations of fluid flow. Along with the basic visualization process which transforms source data into images, there are four additional components to a complete visualization system: Source Data Processing, User Interface and Control, Presentation, and Information Management. The requirements imposed by the desired mode of operation (i.e. real-time, interactive, or batch) and the source data have their effect on each of these visualization system components. The special requirements imposed by the wide variety and size of the source data provided by the numerical simulation of fluid flow presents an enormous challenge to the visualization system designer. We describe the visualization system components including specific visualization techniques and how the mode of operation and source data requirements effect the construction of computational fluid dynamics visualization systems.

Description: Homogeneous turbulent flows with various combinations of strain-rate, rotation rate and coriolis force capture some important aspects of more complex flows with streamline curvature and rotation. Presently, a situation is considered in which as a box of turbulence rotates, strain axes rotate with it. This is to be contrasted with the elliptic streamline flow in which the strain axes are fixed in an inertial frame. The elliptic flow is known to exhibit (inviscid) growth of turbulent energy and one might expect even more rapid growth with the strain-axes following the box. Instead, it is found that the sign of the Reynolds shear stress is reversed leading to a negative production term for turbulent energy. Partial understanding of the phenomenon is obtained from a consideration of the rotation of inertial waves relative to the strain axes as well as the "pressure-less" RDT argument put forward by Cambon etal. [J. Fluid Mech, 278, 175]. Some comparisons with the predictions of second-order closure models will be presented.

Description: The pacing items to reach the highly desirable goal of obtaining computationally accurate flow simulation of a wingtip vortex include; super-computer development, solver accuracy, grid generation and turbulence modeling. In these four areas, many advances have been made but the fact remains that most wing computations are, at best, five percent accurate (in drag coefficient, for example). This level of accuracy has been sufficient for many purposes such as airfoil design, rudimentary wing design, and some forms of optimization. However, this accuracy level will not allow commercial aircraft designers to extract the remaining few percent of efficiency theoretically possible for conventional aircraft configurations. Thus further research is needed, particularly in the areas of solver development and turbulence modeling, to advance the state of the art of viscous computational techniques as applied to problems in aerodynamics. During the course of this study, a substantial amount of measured and computed results have been acquired. In this paper, only a small selection of experimental and computational results will be presented. This paper will outline and discuss a simulation procedure for a wingtip vortex flow analysis using the method of artificial compressibility to solve the three-dimensional, incompressible, Navier-Stokes equations (INS3D-UP).

Description: A new technology has been developed for remote visualization that provides remote, 3D, high resolution, dynamic, interactive viewing of scientific data (such as fluid dynamics simulations or measurements). Based on this technology, some World Wide Web sites on the Internet are providing fluid dynamics data for educational or testing purposes. This technology is also being used for remote collaboration in joint university, industry, and NASA projects in computational fluid dynamics and wind tunnel testing. Previously, remote visualization of dynamic data was done using video format (transmitting pixel information) such as video conferencing or MPEG movies on the Internet. The concept for this new technology is to send the raw data (e.g., grids, vectors, and scalars) along with viewing scripts over the Internet and have the pixels generated by a visualization tool running on the viewer's local workstation. The visualization tool that is currently used is FAST (Flow Analysis Software Toolkit).

Description: Aeroelasticity that involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of aircraft which can experience complex flow/structure interactions. Supersonic transports can experience vortex induced aeroelastic oscillations whereas subsonic transports can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Using these high fidelity methods, design quantities such as structural stresses can be directly computed. Computations using these high fidelity equations require state-of-the-art large computational resources. This paper will describe the transition from potential equations to Euler/Navier-Stokes equations for fluids and from modal equations to finite element equations for structures. Complexities associated with grids for moving control surfaces will be discussed. All contents of the paper will be limited to those already presented at domestic conferences within U.S. such as those sponsored by AIAA. It will not include any materials related to HPCCP computational work.

Description: The influence or rotation on the spectral energy transfer of homogeneous turbulence is investigated in this paper. Given the fact that linear dynamics, e.g. the inertial waves regime tackled in an RDT (Rapid Distortion Theory) fashion, cannot Affect st homogeneous isotropic turbulent flow, the study of nonlinear dynamics is of prime importance in the case of rotating flows. Previous theoretical (including both weakly nonlinear and EDQNM theories), experimental and DNS (Direct Numerical Simulation) results are gathered here and compared in order to give a self-consistent picture of the nonlinear effects of rotation on tile turbulence. The inhibition of the energy cascade, which is linked to a reduction of the dissipation rate, is shown to be related to a damping due to rotation of the energy transfer. A model for this effect is quantified by a model equation for the derivative-skewness factor, which only involves a micro-Rossby number Ro(sup omega) = omega'/(2(OMEGA))-ratio of rms vorticity and background vorticity as the relevant rotation parameter, in accordance with DNS and EDQNM results fit addition, anisotropy is shown also to develop through nonlinear interactions modified by rotation, in an intermediate range of Rossby numbers (Ro(omega) = (omega)' and Ro(omega)w greater than 1), which is characterized by a marco-Rossby number Ro(sup L) less than 1 and Ro(omega) greater than 1 which is characterized by a macro-Rossby number based on an integral lengthscale L and the micro-Rossby number previously defined. This anisotropy is mainly an angular drain of spectral energy which tends to concentrate energy in tile wave-plane normal to the rotation axis, which is exactly both the slow and the two-dimensional manifold. In Addition, a polarization of the energy distribution in this slow 2D manifold enhances horizontal (normal to the rotation axis) velocity components, and underlies the anisotropic structure of the integral lengthscales. Finally is demonstrated the ability of a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) model to predict the underlying spectral transfer structure and all the subsequent developments of classic anisotropy indicators in physical space, when compared to recent LES results. Even if the applications mainly concern developed strong turbulence, a particular emphasis is given to the strong formal analogy of this EDQNM2 model with recent weakly nonlinear approaches to wave-turbulence.

Description: ElVis is a prototype system with allows for the interactive visualization of unsteady fluid flow. The increasing computational power applied to fluid dynamics simulations presents the enormous challenge to the visualization system designer to apply a wide range of technologies to the analysis process with ever increasing demands on performance. Visualization of the results of unsteady fluid flow simulations presents the challenge of exploring very large and complex data sets. Since exploration is a trial and error process, it is of utmost importance that the time required to execute a trial (i.e., create a visualization) be at a minimum in order to provide real time interaction.

Description: Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

Description: Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Description: We consider an isolated bubble boiling regime in which vapour bubbles are intermittently produced at a prearranged set of nucleation site on an upward facing overheated wall plane. In this boiling regime, the bubbles depart from the wall and move as separate entities. Except in the matter of rise velocity, the bubbles do not interfere and are independent of one another. However, the rise velocity is dependent on bubble volume concentration in the bulk. Heat transfer properties specific to this regime cannot be described without bubble detachment size, and we apply our previously developed dynamic theory of vapour bubble growth and detachment to determine this size. Bubble growth is presumed to be thermally controlled. Two limiting cases of bubble evolution are considered: the one in which buoyancy prevails in promoting bubble detachment and the one in which surface tension prevails. We prove termination of the isolated regime of pool nucleate boiling to result from one of the four possible causes, depending on relevant parameters values. The first cause consists in the fact that the upward flow of rising bubbles hampers the downward liquid flow, and under certain conditions, prevents the liquid from coming to the wall in an amount that would be sufficient to compensate for vapour removal from the wall. The second cause is due to the lateral coalescence of growing bubbles that are attached to their corresponding nucleation sites, with ensuing generation of larger bubbles and extended vapour patches near the wall. The other two causes involve longitudinal coalescence either 1) immediately in the wall vicinity, accompanied by the establishment of the multiple bubble boiling regime, or 2) in the bulk, with the formation of vapour columns. The longitudinal coalescence in the bulk is shown to be the most important cause. The critical wall temperature and the heat flux density associated with isolated bubble regime termination are found to be functions of the physical and operating parameters and are discussed in detail.

Description: Flow uniformity in a high enthalpy facility is investigated. The scramjet research facility is composed of a rectangular combustor duct connected to a 100 MW electric arc air heater. The Mach 3.3 flow is accelerated through a two-dimensional contoured nozzle. Instream measurements were made with water-cooled Pitot probes and stagnation point heat flux gages at stream enthalpy levels ranging from 4 to 7 Mj/kg. Flow surveys were made on the flow centerline and off centerline in order to measure the three dimensional uniformity of the flow in the rectangular duct. Measurements indicated that although the flow in the aspect ratio 6:1 duct was relatively uniform on the centerline, three dimensional viscous effects were apparent near the corners. Flow through the nozzle and constant area duct was modeled computationally using a two dimensional, Navier-Stokes, reacting gas code. The computations predict that the flow in the test section is in vibrational equilibrium. The computed and measured Pitot pressure and heat flux profiles are in reasonable agreement with the experimental data.

Description: One concept for high lift aerodynamics is the trapped-vortex. Recently it has been suggested that several spoilers located at different chordwise locations on an airfoil could be deployed to form several cavities in the chordwise direction. This may provide a means of increasing upper surface camber and thereby increase wing lift. It is envisioned that the cavities would be located on an airfoil in its pressure recovery region. Investigations of the effect of one or more cavities on the adjacent boundary layer in an adverse boundary layer were not found in a literature survey. Since flow separation can be caused by flow over cavities, boundary layer measurements were made in the vicinity of cavities which were located in an adverse pressure gradient to study the associated flow in an effort to identify any benefits. The experimental data from the present investigation indicate a minimal impact on the boundary layer profiles is caused by the presence of cavities in either a constant pressure freestream or an adverse pressure gradient. The data from the experimental investigation has been obtained and analyzed. The computational investigation using an incompressible Navier-Stokes code show that reasonable results are being obtained when they are compared with the experimental data.

Description: Large eddy simulation is used to investigate the development of anisotropies and the evolution towards a quasi two-dimensional state in rotating homogeneous turbulence at high Reynolds number. The present study demonstrates the existence of two transitions in the development of anisotropy. The first transition marks the onset of anisotropy and occurs when a macro-Rossby number (based on a longitudinal integral lengthscale) has decreased to near unity while the second transition occurs when a micro-Rossby number (defined in this work as the ratio of the rms fluctuating vorticity to background vorticity) has decreased to unity. The anisotropy marked by the first transition corresponds to a reduction in dimensionality while the second transition corresponds to a polarization of the flow, i.e., relative dominance of the velocity components in the plane normal to the rotation axis. Polarization is reflected by emergence of anisotropy measures based on the two-dimensional component of the turbulence. Investigation of the vorticity structure shows that the second transition is also characterized by an increasing tendency for alignment between the fluctuating vorticity vector and the background angular velocity vector with a preference for corrotative vorticity.

Description: Heat transfer rates specific to nucleate pool boiling under various conditions are determined by the dynamics of vapour bubbles that are originated and grow at nucleation sites of a superheated surface. A new dynamic theory of these bubbles has been recently developed on the basis of the thermodynamics of irreversible processes. In contrast to other existing models based on empirically postulated equations for bubble growth and motion, this theory does not contain unwarrantable assumptions, and both the equations are rigorously derived within the framework of a unified approach. The conclusions of the theory are drastically different from those of the conventional models. The bubbles are shown to detach themselves under combined action of buoyancy and a surface tension force that is proven to add to buoyancy in bubble detachment, but not the other way round as is commonly presumed. The theory ensures a sound understanding of a number of so far unexplained phenomena, such as effect caused by gravity level and surface tension on the bubble growth rate and dependence of the bubble characteristics at detachment on the liquid thermophysical parameters and relevant temperature differences. The theoretical predictions are shown to be in a satisfactory qualitative and quantitative agreement with observations. When being applied to heat transfer at nucleate pool boiling, this bubble dynamic theory offers an opportunity to considerably improve the main formulae that are generally used to correlate experimental findings and to design boiling heat removal in various industrial applications. Moreover, the theory makes possible to pose and study a great deal of new problems of essential impact in practice. Two such problems are considered in detail. One problem concerns the development of a principally novel physical model for the first crisis of boiling. This model allows for evaluating critical boiling heat fluxes under various conditions, and in particular at different gravity levels, with a good agreement with experimental evidence. The other problem bears upon equilibrium shapes of a detached bubble near a heated surface in exceedingly low gravity. In low gravity or in weightlessness, the bubble can remain in the close vicinity of the surface for a long time, and its shape is greatly affected by the Marangoni effect due to both temperature and possible surfactant concentration being nonuniform along the interface. The bubble performs at these conditions like a heat pipe, with evaporation at the bubble lower boundary and condensation at its upper boundary, and ultimately ensures a substantial increase in heat removal as compared with that in normal gravity. Some other problems relevant to nucleate pool and forced convection boiling heat transfer are also discussed.

Description: The discipline of computational fluid dynamics (CFD) is at a crossroad. Most of the significant advances related to computational methods have taken place. The emphasis is now shifting from methods to results. Significant efforts are made in applying CFD to solve design problems. The value of CFD results in design depends on the credibility of computed results for the intended use. The process of establishing credibility requires a standard so that there is a consistency and uniformity in this process and in the interpretation of its outcome. The key element for establishing the credibility is the quantification of uncertainty. This paper presents salient features of a proposed standard and a procedure for determining the uncertainty. A customer of CFD products - computer codes and computed results - expects the following: A computer code in terms of its logic, numerics, and fluid dynamics and the results generated by this code are in compliance with specified requirements. This expectation is fulfilling by verification and validation of these requirements. The verification process assesses whether the problem is solved correctly and the validation process determines whether the right problem is solved. Standards for these processes are recommended. There is always some uncertainty, even if one uses validated models and verified computed results. The value of this uncertainty is important in the design process. This value is obtained by conducting a sensitivity-uncertainty analysis. Sensitivity analysis is generally defined as the procedure for determining the sensitivities of output parameters to input parameters. This analysis is a necessary step in the uncertainty analysis, and the results of this analysis highlight which computed quantities and integrated quantities in computations need to be determined accurately and which quantities do not require such attention. Uncertainty analysis is generally defined as the analysis of the effect of the uncertainties involved in all stages of a process on the final responses. There are two approaches for conducting the uncertainty analysis: experimental and computational. These analyses and approaches are briefly described.

Description: The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Description: Heat exchange to the oscillating flows in regenerative coolers generates entropy. These flows are characterized by oscillating mass flows and oscillating temperatures. Heat is transferred between the flow and heat exchangers and regenerators. In the former case, there is a steady temperature difference between the flow and the heat exchangers. In the latter case, there is no mean temperature difference. In this paper a mathematical model of the entropy generated is developed for both cases. Estimates of the entropy generated by this process are given for oscillating flows in heat exchangers and in regenerators. The practical significance of this entropy is also discussed.

Description: In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Description: A model of the dynamic regime of gas injection through a submerged plate orifice into an ideally wetting liquid is developed in the circumstance when successively detached bubbles may be regarded as independent objects. Two major factors favor bubble detachment: the buoyancy force and a force due to the momentum inflow into the bubble with injected gas. In normal and moderately reduced gravity, the first factor dominates. At relatively low flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole period of its evolution till detachment. If the flow rate exceeds a critical value, another stage of bubble evolution occurs in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical gaseous stem that lengthens as the bubble rises above the plate. The bubble continues to grow until the gas supply through the stem is completely cut off, after which back flow of gas into the stem from the bubble becomes possible. In microgravity, the second factor prevails, and the latter stage is always present irrespective of the flow rate. However, the gas envelope rises and the stem lengthens very slowly. This difference in the underlying physical mechanisms provides for key properties of bubble growth and detachment being drastically different in appreciable and sufficiently reduced gravity. The frequency of bubble formation slightly decreases with and the detachment volume is almost proportional to the gas flow rate in the first case, in accordance with familiar relations. In the second case, the first variable is proportional to the flow rate whereas the second one is independent of it. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Description: Computational fluid dynamics has been developed to a level where it has become an Indispensable part of aerospace research and design. Technology developed foe aerospace applications am also be utilized for the benefit of human health. For example, a flange-to-flange rocket engine fuel-pump simulation includes the rotating and non-rotating components: the flow straighteners, the impeller, and diffusers A Ventricular Assist Device developed by NASA Johnson Space Center and Baylor College of Medicine has a design similar to a rocket engine fuel pump in that it also consists of a flow straightener, an impeller, and a diffuser. Accurate and detailed knowledge of the flowfield obtained by incompressible flow calculations can be greatly beneficial to designers in their effort to reduce the cost and improve the reliability of these devices. In addition to the geometric complexities, a variety of flow phenomena are encountered in biofluids Then include turbulent boundary layer separation, wakes, transition, tip vortex resolution, three-dimensional effects, and Reynolds number effects. In order to increase the role of Computational Fluid Dynamics (CFD) in the design process the CFD analysis tools must be evaluated and validated so that designers gain Confidence in their use. The incompressible flow solver, INS3D, has been applied to flow inside of a liquid rocket engine turbopump components and extensively validated. This paper details how the computational flow simulation capability developed for liquid rocket engine pump component analysis has bean applied to the Left Ventricular Assist Device being developed jointly by NASA JSC and Baylor College of Medicine.

Description: Physical phenomena of fluid flows are often analyzed in computational fluid dynamics. With the progress in computer systems, complex 3D time-varying flow simulations are becoming feasible. There are many existing techniques for visualizing fluid flows; however, most of them do not consider the time-dependent nature of the flow. Hence, they may not provide adequate insights to the flow phenomena. New techniques for visualizing time-varying fluid flows are presented.

Description: The 16-Inch Shock Tunnel at NASA Ames Research Center is a unique test facility used for hypersonic propulsion testing. To provide information necessary to understand the hypersonic testing of the combustor model, computational simulations of the facility nozzle were performed and results are compared with available experimental data, namely static pressure along the nozzle walls and pitot pressure at the exit of the nozzle section. Both quasi-one-dimensional and axisymmetric approaches were used to study the numerous modeling issues involved. The facility nozzle flow was examined for three hypersonic test conditions, and the computational results are presented in detail. The effects of variations in reservoir conditions, boundary layer growth, and parameters of numerical modeling are explored.

Description: This paper is largely based on the paper titled "Guide to Credible Computational Fluid Dynamics Simulations (Invited)" (AIAA Paper 95-2225). The significance of computational fluid dynamics (CFD) simulations depends solely on their credibility. A customer of CFD products -- simulations and software -- expects that these products are credible for the intended use. Simulation model verification and validation are critical in establishing the credibility of simulations and in certifying simulation software. The only thing that matters in establishing the credibility is uncertainty, not veracity or validity. The sensitivity-uncertainty analysis is the key to the establishment of this credibility. Assessing the credibility of complex simulation results poses a significant challenge. Terminology, concepts, framework, and guidelines are presented for addressing this challenge. Verification assesses whether the problem is solved correctly and estimates the level of computational accuracy of simulations; validation determines whether the right problem is solved and assesses the level of validity of the simulation model by estimating the degree to which simulations accurately represent reality. These concepts and the related guidelines are discussed with examples from high-speed flow fields.

Description: Three direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of various irrotational plane strains have been generated. A pseudospectral numerical method with up to 26 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2000. Three different plane strains (of the same magnitude) are imposed. In the first two simulations the strain is in a plane normal to the streamwise wake direction (the two cases having strain of opposite sign); in the third the wake is compressed in the streamwise direction and stretched in the inhomogeneous cross-stream direction. The two flows that are stretched in the cross-stream direction experience an exponential increase of Re; flow visualization indicates many small-scale vortices with little or no organized large-scale structure. In the flow that is compressed in the cross-stream direction Re decays exponentially and the layer appears to be relaminarizing. The evolution of several turbulence statistics in each of these flows is examined.

Description: Numerical simulations of complex 3D time-dependent (unsteady) flows are becoming increasingly feasible because of the progress in computing systems. Unfortunately, many existing flow visualization systems were developed for time-independent (steady) solutions and do not adequately depict solutions from unsteady flow simulations. Furthermore, most systems only handle one time step of the solutions individually and do not consider the time-dependent nature of the solutions. For example, instantaneous streamlines are computed by tracking the particles using one time step of the solution. However, for streaklines and timelines, particles need to be tracked through all time steps. Streaklines can reveal quite different information about the flow than those revealed by instantaneous streamlines. Comparisons of instantaneous streamlines with dynamic streaklines are shown. For a complex 3D flow simulation, it is common to generate a grid system with several millions of grid points and to have tens of thousands of time steps. The disk requirement for storing the flow data can easily be tens of gigabytes. Visualizing solutions of this magnitude is a challenging problem with today's computer hardware technology. Even interactive visualization of one time step of the flow data can be a problem for some existing flow visualization systems because of the size of the grid. Current approaches for visualizing complex 3D time-dependent CFD solutions are described. The flow visualization system developed at NASA Ames Research Center to compute time-dependent particle traces from unsteady CFD solutions is described. The system computes particle traces (streaklines) by integrating through the time steps. This system has been used by several NASA scientists to visualize their CFD time-dependent solutions. The flow visualization capabilities of this system are described, and visualization results are shown.

Description: Previous investigations have shown that a compressible reacting mixing layer can develop two peaks in the mean density weighted vorticity profile. Linear stability analyses show that at these peaks two distinct 'outer' instability modes appear in addition to the more common central mode, which exists unaccompanied in incompressible nonreacting flows. The present study parametrically analyzes the effects of compressibility, heat release, stoichiometry, and density ratio on the amplification rate and obliquity of each stability mode. The mean profiles used in the spatial stability calculation are generated by self-similar solutions of the compressible boundary layer equations combined with the assumption of infinitely fast chemistry. It is shown that the influence of stoichiometry and density ratio on the peaks of the density weighted vorticity profile determines which modes will dominate. Of particular interest are the conditions where two modes are equally amplified, causing the mixing layer to develop into a 'colayer' structure.

Description: We have developed a data parallel multizone compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the "chimera" approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. The design choices can be summarized as: 1. finite differences on structured grids; 2. implicit time-stepping with either distributed solves or data motion and local solves; 3. sequential stepping through multiple zones with interzone data transfer via a distributed data structure. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran (HPF). One interesting feature is the issue of turbulence modeling, where the architecture of a parallel machine makes the use of an algebraic turbulence model awkward, whereas models based on transport equations are more natural. We will present some performance figures for the code on the CM-5, and consider the issues involved in transitioning the code to HPF for portability to other parallel platforms.

Description: A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Description: The previously developed dynamic theory of growth and detachment of vapour bubbles under conditions of nucleate pool boiling is applied to study motion and deformation of a bubble evolving at a single nucleation site. The bubble growth is presumed to be thermally controlled, and two components of heat transfer to the bubble are accounted of: the one from the bulk of surrounding liquid and the one due to heat conduction across a liquid microlayer formed underneath the bubble. Bubble evolution is governed by the buoyancy and an effective surface tension force, both the forces making the bubble centre of mass move away from the wall and, thus, assisting its detachment. Buoyancy-controlled and surface-tension-controlled regimes are considered separately in a meticulous way. The duration of the whole process of bubble evolution till detachment, the rate of growth, and the bubble departure size are found as functions of time and physical and operating parameters. Some repeatedly observed phenomena, such as an influence of gravity on the growth rate, are explained. Inferences of the model agree qualitatively with available experimental evidence, and conclusions pertaining to the dependence on gravity of the bubble radius at detachment and the whole time of the bubble development when being attached to the wall are confirmed quantitatively.

Description: We consider the behaviour of a vapour bubble formed at a nucleation site on a heated horizontal wall. There is no forced convection of an ambient liquid, and the bubble is presumably separated from the wall by a thin liquid microlayer. The energy conservation law results in a variational equation for the mechanical energy of the whole system consisting of the bubble and liquid. It leads to a set of two strongly nonlinear equations which govern bubble expansion and motion of its centre of mass. A supplementary equation to find out the vapour temperature follows from consideration of heat transfer to the bubble, both from the bulk of surrounding liquid and through the microlayer. The average thickness of the microlayer is shown to increase monotonously with time as the bubble meniscus spreads along the wall. Bubble expansion is driven by the pressure head between vapour inside and liquid far away from the bubble, with due allowance for surface tension and gravity effects. It is resisted by inertia of liquid being placed into motion as the bubble grows. The inertia originates also a force that presses the bubble to the wall. This force is counteracted by the buoyancy and an effective surface tension force that tends to transform the bubble into a sphere. The analysis brings about quite a new formulation of the familiar problem of bubble growth and detachment under conditions of nucleate pool boiling.

Description: New exact solutions of the Navier-Stokes equations are obtained for the unbounded and bounded oscillatory and impulsive tangential edgewise motion of touching half-infinite plates in their own plane. In contrast to Stokes classical solutions for the harmonic and impulsive motion of an infinite plane wall, where the solutions are separable or have a simple similarity form, the present solutions have a two-dimensional structure in the near region of the contact between the half-infinite plates. Nevertheless, it is possible to obtain relatively simple closed-form solutions for the flow field in each case by defining new variables which greatly simplify the r- and theta- dependence of the solutions in the vicinity of the contact region. These solutions for flow in a half-infinite space are then extended to bounded flows in a channel using an image superposition technique. The impulsive motion has application to the motion near geophysical faults, whereas the oscillatory motion has arisen in the design of a novel oscillating half-plate flow chamber for examining the effect of fluid shear stress on cultured cell monolayers.

Description: The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed. With the use of the parabolized stability equations (PSE) theory for the numerical computations, either asymptotic, Dirichlet, Neumann or mixed boundary conditions are imposed at various distances from the wall. The results indicate that asymptotic and mixed boundary conditions yield the most accurate mean-flow distortion and unsteady instability modes in comparison with the results obtained with either Dirichlet or Neumann conditions.

Description: Laminar fuel-air counterflow diffusion flames (CFDFs) were studied using axisymmetric convergent-nozzle and straight-tube opposed jet burners (OJBs). The subject diagnostics were used to probe a systematic set of H2/N2-air CFDFs over wide ranges of fuel input (22 to 100% Ha), and input axial strain rate (130 to 1700 Us) just upstream of the airside edge, for both plug-flow and parabolic input velocity profiles. Laser Doppler Velocimetry (LDV) was applied along the centerline of seeded air flows from a convergent nozzle OJB (7.2 mm i.d.), and Particle Imaging Velocimetry (PIV) was applied on the entire airside of both nozzle and tube OJBs (7 and 5 mm i.d.) to characterize global velocity structure. Data are compared to numerical results from a one-dimensional (1-D) CFDF code based on a stream function solution for a potential flow input boundary condition. Axial strain rate inputs at the airside edge of nozzle-OJB flows, using LDV and PIV, were consistent with 1-D impingement theory, and supported earlier diagnostic studies. The LDV results also characterized a heat-release hump. Radial strain rates in the flame substantially exceeded 1-D numerical predictions. Whereas the 1-D model closely predicted the max I min axial velocity ratio in the hot layer, it overpredicted its thickness. The results also support previously measured effects of plug-flow and parabolic input strain rates on CFDF extinction limits. Finally, the submillimeter-scale LDV and PIV diagnostics were tested under severe conditions, which reinforced their use with subcentimeter OJB tools to assess effects of aerodynamic strain, and fueVair composition, on laminar CFDF properties, including extinction.

Description: This talk summarizes three applications of the overset grid method for CFD using some level of automated grid generation, flow solution and post-processing. These applications are 2D high-lift airfoil analysis (INS2D code), turbomachinery applications (ROTOR2/3 codes), and subsonic transport wing/body configurations (OVERFLOW code). These examples provide a forum for discussing the advantages and disadvantages of overset gridding for use in an automated CFD process. The goals and benefits of the automation incorporated in each application will be described, as well as the shortcomings of the approaches.

Description: A numerical study was performed to investigate the accuracy and validity of cold-gas simulation of actual hot scramjet exhaust within a Mach 10 free stream over a representative single-stage-to-orbit airbreathing configuration. In particular, exhausts of various noncombusting chemistry models were studied to characterize their effects on the vehicle aftbody performance and the plume flow field definition. Two approximations of the hot scramjet combustion products were utilized to determine the requirement for expensive, multi-species numerical modeling, and to establish a baseline for the validation of cold-gas simulation. Cold-gas simulation at Mach 10 is shown to be a viable technique using an appropriate thermally perfect gas mixture for reproducing hot scramjet exhaust effects.

Description: The objective of this work is to reshape a nacelle to achieve a specified nacelle pressure distribution. The nacelle may be either isolated or installed on an airplane. There are no restrictions on the attitude (toe, incidence, and roll) and position of the nacelle. The design algorithm is coupled to two different multi-block 3-D Navier Stokes flow solvers. The coupling between design and analysis is automated to the point where the design proceeds with minimal user input.

Description: The Advisory Group for Aerospace Research and Development (AGARD) has requested my participation in the lecture series entitled Parallel Computing in Computational Fluid Dynamics to be held at the von Karman Institute in Brussels, Belgium on May 15-19, 1995. In addition, a request has been made from the US Coordinator for AGARD at the Pentagon for NASA Ames to hold a repetition of the lecture series on October 16-20, 1995. I have been asked to be a local coordinator for the Ames event. All AGARD lecture series events have attendance limited to NATO allied countries. A brief of the lecture series is provided in the attached enclosure. Specifically, I have been asked to give two lectures of approximately 75 minutes each on the subject of parallel solution techniques for the fluid flow equations on unstructured meshes. The title of my lectures is "Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes" (Parts I-II). The contents of these lectures will be largely review in nature and will draw upon previously published work in this area. Topics of my lectures will include: (1) Mesh partitioning algorithms. Recursive techniques based on coordinate bisection, Cuthill-McKee level structures, and spectral bisection. (2) Newton's method for large scale CFD problems. Size and complexity estimates for Newton's method, modifications for insuring global convergence. (3) Techniques for constructing the Jacobian matrix. Analytic and numerical techniques for Jacobian matrix-vector products, constructing the transposed matrix, extensions to optimization and homotopy theories. (4) Iterative solution algorithms. Practical experience with GIVIRES and BICG-STAB matrix solvers. (5) Parallel matrix preconditioning. Incomplete Lower-Upper (ILU) factorization, domain-decomposed ILU, approximate Schur complement strategies.

Description: An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Lagrange multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Description: The level of nitric oxide contamination in the test gas of the NASA Langley Research Center Arc-Heated Scramjet Test Facility and the effect of the contamination on scramjet test engine performance were investigated analytically. The study was conducted for standard facility conditions corresponding to Mach 6, 7, and 8 flight simulations. The analytically determined levels of nitric oxide produced in the facility are compared with experimentally measured levels. Results of the analysis indicate that nitric oxide levels range from one to three mole percent, which corroborates the measured levels. A three-stream combustor code with finite rate chemistry was used to investigate how nitric oxide affects scramjet performance in terms of combustor pressure rise, heat release, and thrust performance. Results indicate minimal effects on engine performance for the test conditions of this investigation.

Description: Direct numerical simulation was used to study the formation and growth of a hairpin vortex in a flat-plate boundary layer and its later development into a young turbulent spot. Fluid injection through a slit in the wall triggered the initial vortex. The legs of the vortex were stretched into a hairpin shape as it traveled downstream. Multiple hairpin vortex heads developed between the stretched legs. New vortices formed beneath the streamwise-elongated vortex legs. The continued development of additional vortices resulted in the formation of a traveling region of highly disturbed ow with an arrowhead shape similar to that of a turbulent spot.

Description: The multiblock three-dimensional Navier-Stokes method PAB3D was utilized by the Component Integration Branch (formerly Propulsion Aerodynamics Branch) at the NASA-Langley Research Center in an international study sponsored by AGARD Working Group #17 for the assessment of the state-of-the-art of propulsion-airframe integration testing techniques and CFD prediction technologies. Three test geometries from ONERA involving fundamental flow physics and four geometries from NASA-LaRC involving realistic flow interactions of wing, body, tail, and jet plumes were chosen by the Working Group. An overview of results on four (1 ONERA and 3 LaRC) of the seven test cases is presented. External static pressures, integrated pressure drag and total drag were calculated for the Langley test cases and jet plume velocity profiles and turbulent viscous stresses were calculated for the ONERA test case. Only selected data from these calculations are presented in this paper. The complete data sets calculated by the participants will be presented in an AGARD summary report. Predicted surface static pressures compared favorably with experimental data for the Langley geometries. Predicted afterbody drag compared well with experiment. Predicted nozzle drag was typically low due to over-compression of the flow near the trailing edge. Total drag was typically high. Predicted jet plume quantities on the ONERA case compared generally well with data.

Description: The fluid management system comprises a mixing/recirculation system including an external recirculation pump for receiving fluid from a zero gravity storage system and returning an output flow of the fluid to the storage system. An internal axial spray injection system is provided for receiving a portion of the output flow from the recirculation pump. The spray injection system thermally de-stratifies liquid and gaseous cryogenic fluid stored in the storage system.

Description: The Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Ames Research Center has begun the development of a Computational Fluid Dynamics (CFD) data set archive on the World Wide Web (WWW) at URL http://www.nas.nasa.gov/NAS/DataSets/. Data sets are integrated with related information such as research papers, metadata, visualizations, etc. In this paper, four classes of users are identified and discussed: students, visualization developers, CFD practitioners, and management. Bandwidth and security issues are briefly reviewed and the status of the archive as of May 1995 is examined. Routine network distribution of data sets is likely to have profound implications for the conduct of science. The exact nature of these changes is subject to speculation, but the ability for anyone to examine the data, in addition to the investigator's analysis, may well play an important role in the future.

Description: Detailed measurements of aerodynamic heating rates in the wake of a Mars-Pathfinder configuration model have been made. Heating data were obtained in a conventional wind tunnel, the NASA LaRC 31" Mach 10 Air Tunnel, and in a high-enthalpy impulse facility, the NASA HYPULSE expansion tube, in which air and CO2 were employed as test gases. The enthalpy levels were 0.7 MJ/kg in the Mach 10 Tunnel, 12 MJ/kg at Mach 9.8 for HYPULSE CO2 tests and 14 MJ/kg at Mach 7.9 for HYPULSE air tests. Wake heating rates were also measured on three similar parametric configurations, and forebody heating measurements were made in order to facilitate CFD comparisons. The ratio of peak wake heating to forebody stagnation point heating in the Mach 10 Tunnel varied from 7% to 15% depending on the freestream Reynolds number. In HYPULSE, the ratio was ~5% for both air and CO 2. It was observed that an increase in the ratio of forebody corner radius to nose radius resulted in a decrease in peak wake heating, and moved the peak closer to the base of the forebody. The wake flow establishment process in HYPULSE was studied, and a method was developed to determine when the wake has become fully established.