Flow in a model of the Space Shuttle Main Engine main injector bowl
Explore the source record for details and available documents.
Thompson, B. E.↗
Engineering topics
Mcdonald, H.
Publications and source records attributed to Mcdonald, H..
Velocity and pressure characteristics of a model SSME high pressure fuel turbopump
Under the present effort an experiment rig has been constructed, an instrumentation package developed and a series of mean and rms velocity and pressure measurements made in a turbopump which modelled the first stage of the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump. The rig was designed so as to allow initial experiments with a single configuration consisting of a bell-mouth inlet, a flight impeller, a vaneless diffuser and a volute. Allowance was made for components such as inlet guide vanes, exit guide vanes, downstream pumps, etc. to be added in future experiments. This flexibility will provide a clear baseline set of experiments and allow evaluation in later experiments of the effect of adding specific components upon the pump performance properties. The rotational speed of the impeller was varied between 4260 and 7680 rpm which covered the range of scaled SSME rotation speeds when due allowance is made for the differing stagnation temperature, model to full scale. The results at the inlet obtained with rotational speeds of 4260, 6084 and 7680 rpm showed that the axial velocity at the bell-mouth inlet remained roughly constant at 2.2 of the bulk velocity at the exit of the turbopump near the center of the inlet, but it decreased rapidly with increasing radius at all three speeds. Reverse flow occurred at a radius greater than 0.9 R for all three speeds and the maximum negative velocity reduced from 1.3 of the bulk velocity at the exit of the turbopump at 4260 rpm to 0.35 at 7680 rpm, suggesting that operating at a speed closer to the design condition of 8700 rpm improved the inlet characteristics. The reverse flow caused positive prerotation at the impeller inlet which was negligibly small near the center but reached 0.7 of the impeller speed at the outer annulus. The results in the diffuser and the volute obtained at 7680 rpm show that the hub and shroud walls of the diffuser were characterized by regions of transient reverse flow with negative revolution-averaged velocity of 8 percent of the maximum forward revolution-averaged velocity at the center of the diffuser passage near the shroud wall.
Tse, D. G-N.↗
Reciprocal interactions of hairpin-shaped vortices and a boundary layer
The present effort simulates the flow dynamics resulting from hairpin-shaped vortices in a boundary layer through the solution of the time-dependent, three-dimensional, compressible Navier-Stokes equations. Two simulations have been carried out. In the first case, the initial condition contains only one imposed hairpin-shaped vortex. In the second case, there are two incipient vortices separated by a short distance in the streamwise direction. The results add to the detailed understanding of the dynamic role of hairpin vortices in the entrainment, production and continual replenishment of turbulent structure, as well as the process of the interaction, intertwining and eventual conglomeration of vortices in boundary layers.
Liu, N.-S.↗
Navier-Stokes analysis of solid propellant rocket motor internal flows
A multidimensional implicit Navier-Stokes analysis that uses numerical solution of the ensemble-averaged Navier-Stokes equations in a nonorthogonal, body-fitted, cylindrical coordinate system has been applied to the simulation of the steady mean flow in solid propellant rocket motor chambers. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. The code was validated by comparing computed results with the experimental data obtained in cylindrical-port cold-flow tests. The agreement between the computed and experimentally measured mean axial velocities is excellent. The axial location of transition to turbulent flow predicted by the two-equation (k-epsilon) turbulence model used in the computations also agrees well with the experimental data. Computations performed to simulate the axisymmetric flowfield in the vicinity of the aft field joint in the Space Shuttle solid rocket motor using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot.
Sabnis, J. S.↗
Solution of the incompressible Navier-Stokes equations using artificial compressibility methods
A modified artificial compressibility formulation is derived from a preconditioned low Mach number compressible formulation, and these two methods are compared on a two-dimensional laminar leading edge flow using a LBI/ADI solution algorithm. The two methods are essentially equivalent with appropriate preconditioning parameters and have the same convergence rates and efficiency, giving 4 orders of residual reduction in about 75 iterations with a vectorized CRAY-XMP run time of 20 seconds for 3000 grid points.
Briley, W. R.↗
Flow in a model of the SSME main injector bowl
Mean-flow and turbulence quantities obtained in a 32-percent acrylic model of the SSME main injector bowl are presented. Flow visualization obtained with laser sheets and conventional photography shows the basic patterns and identifies regions of recirculating flow. Measurements of mean-flow and turbulence normal stresses are obtained with laser velocimetry to confirm these results and to quantify flow characteristics in the inlet ducts, entrance region, racetrack, LOX-post bundle, and exit nozzle. Measured velocity distributions show the upstream effect of the injector bowl in the transfer ducts, the plenum effect of the racetrack and LOX-post bundle, and the nonuniform distribution of flow into the main exit nozzle.
Thompson, B. E.↗
Two- and three-dimensional turbine blade row flow field simulations
Work performed in the numerical simulation of turbine passage flows via a Navier-Stokes approach is discussed. Both laminar and turbulent simulations in both two and three dimensions are discussed. An outline of the approach, background, and an overview of the results are given.
Buggeln, R. C.↗
A combined Eulerian-Lagrangian analysis for computation of two-phase flows
A combined Eulerian-Lagrangian analysis, which combines a linearized block implicit Navier-Stokes analysis for the continuous phase with a Lagrangian analysis for the motion of the droplet phase to simulate evaporating two-phase flows, has been developed. A unique aspect of this analysis is that the Lagrangian equations for the droplet motion have been transformed into the Eulerian computational space using coordinate transformation resulting in better computational efficiency. Use of the present implicit procedure for the continuous phase makes it possible to efficiently use a locally highly refined mesh. The coupling of the continuous phase and the droplet phase analyses is such that it allows the use of convergence acceleration techniques for steady-state problems, as well as making it possible to simulate transient flows. A rapidly evaporating two-phase flow caused by spray of liquid nitrogen in gaseous nitrogen has been used as a test problem to demonstrate the computer code using this analysis.
Sabnis, J. S.↗
Navier-Stokes analysis of two- and three-dimensional flow field in solid rocket motors with segment joints
A multidimensional implicit Navier-Stokes analysis which uses numerical solution of ensemble-averaged Navier-Stokes equations in a nonorthogonal bodyfitted cylindrical-polar coordinate system has been applied to simulation of the internal flow field in solid-propellant rocket motor chambers with segment joints. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. Computations performed to simulate the axisymmetric flow field in the vicinity of the aft field joint in the Space Shuttle SRB using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot. Calculations were also performed for two cases involving asymmetric three-dimensional flow in the vicinity of the aft field joint in the SRB using 721,525 grid points to estimate circumferential velocities and pressure gradients at the joint.
Sabnis, J. S.↗
Computation of multi-dimensional viscous supersonic jet flow
A new method has been developed for two- and three-dimensional computations of viscous supersonic flows with embedded subsonic regions adjacent to solid boundaries. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases relevant to internal supersonic flow is presented and compared with data. Comparison between data and computation are in general excellent thus indicating that the computational technique has great promise as a tool for calculating supersonic flow with embedded subsonic regions. Finally, a User's Manual is presented for the computer code used to perform the calculations.
Kim, Y. N.↗
Computation of multi-dimensional viscous supersonic flow
A method has been developed for two- and three-dimensional computations of viscous supersonic jet flows interacting with an external flow. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases associated with supersonic jet flow is presented and compared with other calculations for axisymmetric cases. Demonstration calculations indicate that the computational technique has great promise as a tool for calculating a wide range of supersonic flow problems including jet flow. Finally, a User's Manual is presented for the computer code used to perform the calculations.
Buggeln, R. C.↗
Turbine stator flow field simulations
The increased capability and accessibility of modern computers, coupled with increasingly sophisticated and accurate numerical and physical modeling, has led to a marked impact of numerical simulations upon current turbine design and research problems. The turbine section represents a considerable challenge as it contains significant regions of complex three-dimensional flow, including both aerodynamic and heat transfer phenomena. The focus of the present effort is the development of an efficient and accurate three-dimensional Navier-Stokes calculation procedure for application to the turbine stator and rotor problems. In particular, an effective procedure is sought which: (1) adequately represents the flow physics, (2) allows for sufficient resolution in regions of small length scale, and (3) has sufficiently good convergence properties so as to allow use on a regular basis.
Buggeln, R. C.↗
Numerical analysis of some supersonic viscous flows related to inlet and nozzle systems
A numerical method originally developed for three-dimensional supersonic inlet flow calculations is extended and applied to the study of two-dimensional and three-dimensional flows associated with arbitrary propulsion systems. The method is based on the forward spatial marching solution of a reduced form of the three-dimensional steady Navier-Stokes equations in which streamwise pressure gradients are retained in both the subsonic and supersonic regions. The present paper briefly describes the analysis and then shows three applications. In the first application, a wall transpiration study has been performed for the two-dimensional shock wave/turbulent boundary layer interaction flow field with application to an inlet configuration. The second application treats a rectangular high speed inlet with a swept sideplate including the effects of sideplate spillage. Finally, the method is utilized to analyze the interaction of an under-expanded supersonic jet with an ambient flow. Computed results are examined and compared with available experimental measurements. It is demonstrated that the present numerical method is capable of numerically simulating complex two- and three-dimensional flows relevant to hypersonic propulsion systems in a manner which both shows good agreement with data when such data is available, and which shows the complex flow features in the absence of data.
Kim, Y.-N.↗
Influence of large-scale motion on turbulent transport for confined coaxial jets. Volume 2: Navier-Stokes calculations of swirling and nonswirling confined coaxial jets
The existence of large scale coherent structures in turbulent shear flows has been well documented. Discrepancies between experimental and computational data suggest a necessity to understand the roles they play in mass and momentum transport. Using conditional sampling and averaging on coincident two-component velocity and concentration velocity experimental data for swirling and nonswirling coaxial jets, triggers for identifying the structures were examined. Concentration fluctuation was found to be an adequate trigger or indicator for the concentration-velocity data, but no suitable detector was located for the two-component velocity data. The large scale structures are found in the region where the largest discrepancies exist between model and experiment. The traditional gradient transport model does not fit in this region as a result of these structures. The large scale motion was found to be responsible for a large percentage of the axial mass transport. The large scale structures were found to convect downstream at approximately the mean velocity of the overall flow in the axial direction. The radial mean velocity of the structures was found to be substantially greater than that of the overall flow.
Weinberg, B. C.↗
Numerical and experimental investigation of nonswirling and swirling confined jets
An investigation of the influence of large scale structures on the flow development for coaxial jets with sudden expansion (with and without swirl) is presented. Both an experimental study and numerical predictions were performed for a configuration corresponding to that considered by Johnson and Bennett and Roback and Johnson. The effects of large scale structures on the swirling and nonswirling flows are documented, in particular their influence on turbulence modeling and the numerical simulation. The ensemble-averaged, time-dependent Navier-Stokes equations are solved by an LBI procedure to predict the turbulent flow field. Effects of artificial dissipation and placement of the upstream boundary in the numerical computation are also discussed.
Brondum, D. C.↗
Turbine Vane External Heat Transfer. Volume 2. Numerical Solutions of the Navier-stokes Equations for Two- and Three-dimensional Turbine Cascades with Heat Transfer
The application of the time-dependent ensemble-averaged Navier-Stokes equations to transonic turbine cascade flow fields was examined. In particular, efforts focused on an assessment of the procedure in conjunction with a suitable turbulence model to calculate steady turbine flow fields using an O-type coordinate system. Three cascade configurations were considered. Comparisons were made between the predicted and measured surface pressures and heat transfer distributions wherever available. In general, the pressure predictions were in good agreement with the data. Heat transfer calculations also showed good agreement when an empirical transition model was used. However, further work in the development of laminar-turbulent transitional models is indicated. The calculations showed most of the known features associated with turbine cascade flow fields. These results indicate the ability of the Navier-Stokes analysis to predict, in reasonable amounts of computation time, the surface pressure distribution, heat transfer rates, and viscous flow development for turbine cascades operating at realistic conditions.
Yang, R. J.↗
Calculations of two and three-dimensional transonic cascade flow fields using the Navier-Stokes equations
The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations have been used to study the turbulent flow field in two and three-dimensional turbine cascades. The viscous regions of the flow were resolved and non-slip boundary conditions were utilized on solid surfaces. The calculations were performed in a constructive 'O'-type grid which allows representation of the blade rounded trailing edge. Converged solutions were obtained in relatively few time steps (about 80-150) and comparisons for both surface pressure and heat transfer showed good agreement with data. The three-dimensional turbine cascade calculation showed many of the expected flow-field features.
Weinberg, B. C.↗
A solution procedure for two- and three-dimensional unsteady viscous flows
An efficient computational procedure for solving three-dimensional unsteady turbulent flows is described. The consistently split Linearized Block Implicit (LBI) scheme is used in conjunction with the QR operator scheme to solve an approximate form of the Navier-Stokes equations in generalized nonorthogonal coordinates employing physical velocity components. As a demonstration calculation the turbulent oscillating flow over a flat plate corresponding to the experiment of Karlsson is considered in both two and three dimensions. New inflow boundary conditions are proposed which yield physically plausible solutions near the upstream boundary. The results obtained agree both qualitatively and quantitatively with Karlsson's data and shed new light on the controversy concerning the interpretation of the skin friction phase angle as a function of reduced frequency.
Weinberg, B. C.↗