Search NASASearch

Engineering topics

Chima, Rodrick V.

Publications and source records attributed to Chima, Rodrick V..

At least 37 records · Page 2

Preconditioning for Numerical Simulation of Low Mach Number Three-Dimensional Viscous Turbomachinery Flows

A preconditioning scheme has been implemented into a three-dimensional viscous computational fluid dynamics code for turbomachine blade rows. The preconditioning allows the code, originally developed for simulating compressible flow fields, to be applied to nearly-incompressible, low Mach number flows. A brief description is given of the compressible Navier-Stokes equations for a rotating coordinate system, along with the preconditioning method employed. Details about the conservative formulation of artificial dissipation are provided, and different artificial dissipation schemes are discussed and compared. The preconditioned code was applied to a well-documented case involving the NASA large low-speed centrifugal compressor for which detailed experimental data are available for comparison. Performance and flow field data are compared for the near-design operating point of the compressor, with generally good agreement between computation and experiment. Further, significant differences between computational results for the different numerical implementations, revealing different levels of solution accuracy, are discussed.

Tweedt, Daniel L.

SWIFT: Multistage Turbomachinery Analysis Code Developed

A new turbomachinery analysis code called SWIFT has been developed at the NASA Lewis Research Center. SWIFT solves thin-layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and an explicit finite-difference scheme. Preconditioning allows the code to be used for all speed ranges, from incompressible to supersonic flows. Multiblock capability allows three types of grids to be patched together to simulate many types of turbomachinery geometries, including hub and tip clearances. In addition, code users can analyze multistage turbomachinery by using a steady averaging-plane approach. This approach uses linearized characteristic boundary conditions to pass information accurately between the stages.

Chima, Rodrick V.

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.

Rapid Numerical Simulation of Viscous Axisymmetric Flow Fields

A two-dimensional Navier-Stokes code has been developed for rapid numerical simulation of axisymmetric flow fields, including flow fields with an azimuthal velocity component. The azimuthal-invariant Navier-Stokes equations in a cylindrical coordinate system are mapped to a general body-fitted coordinate system, with the streamwise viscous terms then neglected by applying the thin-layer approximation. Turbulence effects are modeled using an algebraic model, typically the Baldwin-Lomax turbulence model, although a modified Cebeci-Smith model can also be used. The equations are discretized using central finite differences and solved using a multistage Runge-Kutta algorithm with a spatially varying time step and implicit residual smoothing. Results are presented for calculations of supersonic flow over a waisted body-of-revolution, transonic flow through a normal shock wave in a straight circular duct of constant cross sectional area, swirling supersonic (inviscid) flow through a strong shock in a straight radial duct, and swirling subsonic flow in an annular-to-circular diffuser duct. Comparisons between computed and experimental results are in fair to good agreement, demonstrating that the viscous code can be a useful tool for practical engineering design and analysis work.

Tweedt, Daniel L.

A k-Omega Turbulence Model for Quasi-Three-Dimensional Turbomachinery Flows

A two-equation k-omega turbulence model has been developed and applied to a quasi-three-dimensional viscous analysis code for blade-to-blade flows in turbomachinery. the code includes the effects of rotation, radius change, and variable stream sheet thickness. The flow equations are given and the explicit runge-Kutta solution scheme is described. the k-omega model equations are also given and the upwind implicit approximate-factorization solution scheme is described. Three cases were calculated: transitional flow over a flat plate, a transonic compressor rotor, and transonic turbine vane with heat transfer. Results were compared to theory, experimental data, and to results using the Baldwin-Lomax turbulence model. The two models compared reasonably well with the data and surprisingly well with each other. Although the k-omega model behaves well numerically and simulates effects of transition, freestream turbulence, and wall roughness, it was not decisively better than the Baldwin-Lomax model for the cases considered here.

Chima, Rodrick V.

The effect of adding roughness and thickness to a transonic axial compressor rotor

The performance deterioration of a high speed axial compressor rotor due to surface roughness and airfoil thickness variations is reported. A 0.025 mm (0.001 in.) thick rough coating with a surface finish of 2.54-3.18 RMS microns (100-125 RMS microinches) is applied to the pressure and suction surface of the rotor blades. Coating both surfaces increases the leading edge thickness by 10% at the hub and 20% at the tip. Application of this coating results in a loss in efficiency of 6 points and a 9% reduction in the pressure ratio across the rotor at an operating condition near the design point. To separate the effect of thickness and roughness, a smooth coating of equal thickness is also applied to the blade. The smooth coating surface finish is 0.254-0.508 RMS microns (10-20 RMS microinches), compared to the bare metal blade surface finish of 0.508 RMS microns (20 RMS microinches). The smooth coating results in approximately half of the performance deterioration found from the rough coating. Both coatings are then applied to different portions of the blade surface to determine which portions of the airfoil are most sensitive to thickness/roughness variations. Aerodynamic performance measurements are presented for a number of coating configurations at 60%, 80%, and 100% of design speed. The results indicate that thickness/roughness over the first 10% of blade chord accounts for virtually all of the observed performance degradation for the smooth coating, compared to about 70% of the observed performance degradation for the rough coating. The performance deterioration is investigated in more detail at design speed using laser anemometer measurements as well as predictions generated by a quasi-3D Navier-Stokes flow solver which includes a surface roughness model. Measurements and analysis are performed on the baseline blade and the full-coverage smooth and rough coatings. The results indicate that coating the blade causes a thickening of the blade boundary layers. The interaction between the rotor passage shock and the thickened suction surface boundary layer then results in an increase in blockage which reduces the diffusion level in the rear half of the blade passage, thus reducing the aerodynamic performance of the rotor.

Suder, Kenneth L.

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.

An algebraic turbulence model for turbomachinery

This paper presents a description and verification of RVC3D (rotor viscous code 3-D) which provides a Euler or Navier-Stokes analysis for steady three dimensional flows in turbomachinery. A motivation for this analysis is the calculation of turbine endwall heat transfer. Features of the turbulence model code include thin-layer formulation, Baldwin-Lomax or Cebeci-Smith turbulence models, node-centered finite difference formulation, and explicit four-stage Runge-Kutta time marching scheme. Results for flat plate, annular turbine cascade, turbine endwall heat transfer, and supersonic compressor blade test cases are presented.

Chima, Rodrick V.

TCGRID: A three dimensional C-grid generator for turbomachinery

A fast 3-D grid code for turbomachinery was developed. The code, TCGRID (Turbomachinery C-GRID), can generate either C or H type grids for fairly arbitrary axial or radical turbomachinery geometries. The code also has limited blocked grid capability and can generate an axisymmetric H type grid upstream of the blade row or an O type grid within the tip clearance region. Hub and tip geometries are input as a simple list of pairs. All geometric data is handled using parametric splines so that geometries that turn 90 degrees can be handled without difficulty. Blade input is in standard MERIDL or Lewis compressor design code format. TCGRID adds leading and trailing edge circles to MERIDL geometries and intersects the blade with the hub and tip if necessary using a novel intersection algorithm. The procedure used to generate the grid is given. Output is in PLOT3D format, which can also be read by the RVC3D (Rotor Viscous Code 3-D) Navier-Stokes code for turbomachinery. Intermediate 2-D or 3-D grids useful for debug and other purposes can also be output using a convenient output flag. A grid generated figure is given.

Chima, Rodrick V.

Viscous three-dimensional calculations of transonic fan performance

A 3-D flow analysis code was used to compute the design speed operating line of a transonic fan rotor, and the results were compared with experimental data. The code is an explicit finite difference code with an algebraic turbulence model. The transonic fan, called Rotor 67, was tested experimentally at NASA Lewis conventional aerodynamic probes and with user anemometry and was included as one of the AGARD test cases for the computation of internal flows. The experimental data are described. Maps of total pressure ratio and adiabatic efficiency vs mass flow were computed and are compared with the experimental maps, with good agreement. Detailed comparisons between calculations and experiment are made at two operating points, one near peak efficiency and the other near stall. Blade-to-blade contour plots are used to show the shock structure. Comparisons of circumferentially integrated flow quantities downstream of the rotor show spanwise distributions of several aerodynamic parameters. Calculated Mach number distributions are compared with laser anemometer data within the blade row and the wake to quantify the accuracy of the calculations. Particle traces are used to show the nature of secondary flow.

Chima, Rodrick V.

Viscous three-dimensional calculations of transonic fan performance

A 3-D flow analysis code was used to compute the design speed operating line of a transonic fan rotor, and the results were compared with experimental data. The code is an explicit finite difference code with an algebraic turbulence model. The transonic fan, called rotor 67, was tested experimentally at NASA-Lewis with conventional aerodynamic probes and with user anemometry and was included as one of the AGARD test cases for the computation of internal flows. The experimental data are described. Maps of total pressure ratio and adiabatic efficiency versus mass flow were computed and are compared with the experimental maps, with good agreement. Detailed comparisons between calculations and experiment are made at two operating points, one near peak efficiency and the other near stall. Blade-to-blade contour plots are used to show the shock structure. Comparisons of circumferentially integrated flow quantities downstream of the rotor show spanwise distributions of several aerodynamic parameters. Calculated Mach number distributions are compared with laser anemometer data within the blade row and the wake to quantify the accuracy of the calculations. Particle traces are used to show the nature of secondary flow.

Chima, Rodrick V.

Analysis of three-dimensional viscous flow in a supersonic throughflow fan

A 3-D Navier-Stokes code was developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system rotating about the x-axis, and then mapped to a general body-fitted coordinate system. Streamwise viscous terms are neglected using the thin layer assumption, and turbulence effects are modelled using the Baldwin-Lomax turbulence model. The equations are discretized using finite differences on stacked C-type grids and are solved using a multistage Runge-Kutta algorithm with a spatially varying time step and implicit residual smoothing. Calculations were made of the flow around a supersonic throughflow fan blade. The fan was designed as a key component in a supersonic cruise engine. The 3-D calculations were done on a 129x29x33 grid and took 50 minutes of cpu time. Comparisons with the quasi-3-D results show minor differences in loading due to 3-D effects. Particle traces show nearly 2-D flows near the pressure surface, but large secondary flows within the suction surface boundary layer. The horseshoe vortex ahead of the leading edge is clearly seen.

Chima, Rodrick V.

An unconditionally stable Runge-Kutta method for unsteady rotor-stator interaction

A quasi-three-dimensional analysis has been developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stress-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using an explicit four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing is used to increase the stability limit of the time-accurate computations. The scheme is described, and stability and accuracy analyses are given.

Chima, Rodrick V.

An unconditionally stable Runge-Kutta method for unsteady flows

A quasi-three-dimensional analysis was developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream surface thickness. The Baldwin-Lomax eddy viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing was employed to accelerate the solution of the time accurate computations. The scheme is described and accuracy analyses are given. Results are shown for a supersonic through-flow fan designed for NASA Lewis. The rotor:stator blade ratio was taken as 1:1. Results are also shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Here the blade ratio is 2:3. Implicit residual smoothing was used to increase the time step limit of the unsmoothed scheme by a factor of six with negligible differences in the unsteady results. It is felt that the implicitly smoothed Runge-Kutta scheme is easily competitive with implicit schemes for unsteady flows while retaining the simplicity of an explicit scheme.

Jorgenson, Philip C. E.

An explicit Runge-Kutta method for unsteady rotor/stator interaction

A quasi-three-dimensional rotor/stator analysis has been developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using a four-stage Runge-Kutta scheme with a constant timestep. Results are shown for the first stage of the Space Shuttle Main Engine high pressure fuel turbopump. Euler and Navier-Stokes results are compared on the scaled single- and multi-passage machine. The method is relatively fast and the quasi-three-dimensional formulation is applicable to a wide range of turbomachinery geometries.

Jorgenson, Philip C. E.