Engineering topics
Publications and source records attributed to Sondak, D. L..
Contents include the following: Motivation. Code details. Test cases: aardvark (2D), phantom (3D). Algorithm basics. New feature I, II, and III. UTRC LSRR 1 - 1/2 stage turbine. Supersonic turbine. UTRC LSRR single-stage turbine. LOX elbow duct. Computational grid. High-pressure fuel pump inducer. Development plan. Summary. Appendix - Background.
The turbines used in rocket-engine applications are often partial-admission turbines, meaning that the flow enters the rotor over only a portion of the annulus. These turbines have been traditionally analyzed, however, assuming full-admission characteristics. This assumption enables the simulation of only a portion of the 360-degree annulus, with periodic boundary conditions applied in the circumferential direction. While this traditional approach to the simulating the flow in partial-admission turbines significantly reduces the computational requirements, the accuracy of the solutions has rarely been evaluated. In the current investigation, both full- and partial-admission three dimensional unsteady Navier-Stokes simulations were performed for a partial-admission turbine designed and tested at NASA Marshall Space Flight Center. The results indicate that the partial-admission nature of the turbine must be included in simulations to properly predict the performance and flow unsteadiness of the turbine.
A k-epsilon turbulence model suitable for compressible flow, including the new wall function formulation, has been incorporated into an existing compressible Reynolds-averaged Navier-Stokes code, F3D. The low Reynolds number k-epsilon model of Chien (1982) was added for comparison with the present method. A number of features were added to the F3D code including improved far-field boundary conditions and viscous terms in the streamwise direction. A series of computations of increasing complexity was run to test the effectiveness of the new formulation. Flow over a flat plate was computed by using both orthogonal and nonorthogonal grids, and the friction coefficients and velocity profiles compared with a semi-empirical equation. Flow over a body of revolution at zero angle of attack was then computed to test the method's ability to handle flow over a curved surface. Friction coefficients and velocity profiles were compared to test data. All models gave good results on a relatively fine grid, but only the wall function formulation was effective with coarser grids. Finally, in order to demonstrate the method's ability to handle complex flow fields, separated flow over a prolate spheroid at angle of attack was computed, and results were compared to test data. The results were also compared to a k-epsilon model by Kim and Patel (1991), in which one equation model patched in at the wall was employed. Both models gave reasonable solutions, but improvement is required for accurate prediction of friction coefficients in the separated regions.