Transient Thermal and Structural Analysis of the Liquid Rocket Turbopump Turbine

Author(s):  
S. H. Yoon ◽  
S. M. Jeon ◽  
J. Kim

A turbopump unit is a key component of the liquid rocket engine assembly. Generally operated by an impulse-type turbine with supersonic hot gas for a gas generator cycle engine, it pressurizes liquid oxygen and fuel to attain high specific impulse. The turbine is exposed to hot gas while rotating at a high speed and experiences severe thermal and mechanical loads during the operation. Its structural integrity must be assured through thermal and structural analyses of the turbine. To investigate non-steady heat transfer characteristics from the hot gas to the turbine blade, the heat transfer of supersonic hot gas flow to the complicated turbine blade geometry is solved through a three-dimensional CFD analysis. The distribution of gas temperature, pressure and velocity on the blade surface are obtained by the CFD analysis with appropriate convection heat transfer coefficients to be used for the structural analysis of the turbine rotor. The thermal stress due to material expansion and mechanical stress due to centrifugal loads are considered in the analysis. The three-dimensional transient stress distributions of the turbine blades, shroud, and disk are obtained by a finite element analysis. Finally, the overall results are examined to evaluate the strength safety of the turbine.

2001 ◽  
Author(s):  
Thomas J. Martin ◽  
George S. Dulikravich

Abstract This paper presents the theoretical methodology, conceptual demonstration, and validation of a fully automated computer program for the inverse design and optimization of internal convectively cooled three-dimensional axial gas turbine blades. A parametric computer model of the three-dimensional internal cooling network was developed, including the automatic generation of computational grids. A boundary element computer program was written to solve the steady-state non-linear heat conduction equation inside the internally cooled and thermal barrier-coated turbine blade. A finite element algorithm was written to model an arbitrary network of internal coolant passages for the calculation of the internal heat transfer coefficients, pressure losses, local flow rates, the effects of centrifugal pumping, heating of the coolant fluid, and forced convection from the thermal model of the solid to the coolant fluid. The heat conduction and internal flow analyses were iteratively coupled to account for the heat balance between the blade and the coolant fluid. The computer-automated design and optimization system was demonstrated on the second high-pressure turbine blade row of the Pratt & Whitney F100 engine. The internal cooling configuration and local heat transfer enhancements (boundary layer trip strips and pin fins) inside the three-dimensional blade were optimized for maximum cooling effectiveness and durability against corrosion and thermo-mechanical fatigue.


2017 ◽  
Vol 11 (10) ◽  
pp. 189 ◽  
Author(s):  
Sutrisno Sutrisno ◽  
Deendarlianto Deendarlianto ◽  
Indarto Indarto ◽  
Sigit Iswahyudi ◽  
Muhammad Agung Bramantya ◽  
...  

The research on three dimensional (3-D) wind turbine blades has been introduced (Sutrisno, Prajitno, Purnomo, & B.W. Setyawan, 2016). In the current experiment, the 3-D wind turbine blades would be fitted with helicopter-like blade tips and additional fins to the blade hubs to demonstrate some laminarizing features. It was found that additional helicopter-like blade tip to the turbine blade creates strong laminar flows over the surface of the blade tips. Supplementary, finned hub, fitted to the blade body, creates rolled-up vortex flows, weakens the blade stall growth development, especially for blades at high-speed wind. A proposed mathematical form of modified lifting line model has been developed to pursue further 3-d blade development study of 3-d wind turbine blade. Rolled up vortex effects, developed by finned of the base hub, has been acknowledged could demolish the turbulent region, as well as laminarize the stall domain to intensify the induced wind turbine blade lift.


2016 ◽  
Vol 138 (11) ◽  
Author(s):  
Matteo Giovannini ◽  
Michele Marconcini ◽  
Filippo Rubechini ◽  
Andrea Arnone ◽  
Francesco Bertini

The present activity was carried out in the framework of the Clean Sky European Research Project ITURB (optimal high-lift turbine blade aeromechanical design), aimed at designing and validating a turbine blade for a geared open-rotor engine. A cold-flow, large-scale, low-speed (LS) rig was built in order to investigate and validate new design criteria, providing reliable and detailed results while containing costs. This paper presents the design of an LS stage and describes a general procedure that allows to scale three-dimensional (3D) blades for LS testing. The design of the stator row was aimed at matching the test-rig inlet conditions and at providing the proper inlet flow field to the blade row. The rotor row was redesigned in order to match the performance of the high-speed (HS) configuration, compensating for both the compressibility effects and different turbine flow paths. The proposed scaling procedure is based on the matching of the 3D blade loading distribution between the real engine environment and the LS facility one, which leads to a comparable behavior of the boundary layer and hence to comparable profile losses. To this end, the datum blade is parameterized, and a neural-network-based methodology is exploited to guide an optimization process based on 3D Reynolds-averaged Navier–Stokes (RANS) computations. The LS stage performance was investigated over a range of Reynolds numbers characteristic of modern low-pressure turbines (LPTs) by using a multi-equation, transition-sensitive, turbulence model. Some comparisons with experimental data available within the project finally proved the effectiveness of the proposed scaling procedure.


Author(s):  
B. Woerz ◽  
Y. Mick ◽  
E. Findeisen ◽  
P. Jeschke ◽  
M. Rabs

This paper presents different numerical methods to predict the thermal load of a convection cooled gas turbine blade under realistic operating temperature conditions. The subject of the investigation is a gas turbine rotor blade which is equipped with a state-of-the-art convection cooling system. Firstly, two FEM based methods are introduced. One method, referred to as FEM1D method, uses empirical correlations from the open literature to obtain one dimensional heat transfer coefficients along one flow line inside the cooling channels while in the hot gas path a three dimensional CFD simulation is used. The second method (FEM2D) uses three dimensional CFD simulations to obtain two dimensional heat transfer coefficient distributions for both, the inner cooling channels and the hot gas path. The results from both numerical methods are compared with each other and are validated with experimental data, quantifying also their accuracy limits. In total this paper gives an evaluation of two different FEM methods to predict temperature distribution in convection cooled gas turbines. Their accuracy, numerical cost and limitations are evaluated. It turns out that the temperature profiles gained by both methods are generally in good agreement with the experiments. However, while causing higher numerical costs the more detailed FEM2D method achieves more accurate results.


Author(s):  
Vijay K. Garg ◽  
Ali A. Ameri

A three-dimensional Navier-Stokes code has been used to compute the heat transfer coefficient on two film-cooled turbine blades, namely the VKI rotor with six rows of cooling holes including three rows on the shower head, and the C3X vane with nine rows of holes including five rows on the shower head. Predictions of heat transfer coefficient at the blade surface using three two-equation turbulence models, specifically, Coakley’s q-ω model, Chien’s k-ε model and Wilcox’s k-ω model with Menter’s modifications, have been compared with the experimental data of Camci and Arts (1990) for the VKI rotor, and of Hylton et al. (1988) for the C3X vane along with predictions using the Baldwin-Lomax (B-L) model taken from Garg and Gaugler (1995). It is found that for the cases considered here the two-equation models predict the blade heat transfer somewhat better than the B-L model except immediately downstream of the film-cooling holes on the suction surface of the VKI rotor, and over most of the suction surface of the C3X vane. However, all two-equation models require 40% more computer core than the B-L model for solution, and while the q-ω and k-ε models need 40% more computer time than the B-L model, the k-ω model requires at least 65% more time due to slower rate of convergence. It is found that the heat transfer coefficient exhibits a strong spanwise as well as streamwise variation for both blades and all turbulence models.


Author(s):  
Alka Gupta ◽  
Abdulrahman Alsultan ◽  
R. S. Amano ◽  
Sourabh Kumar ◽  
Andrew D. Welsh

Energy is the heart of today’s civilization and the demand seems to be increasing with our growing population. Alternative energy solutions are the future of energy, whereas the fossil-based fuels are finite and deemed to become extinct. The design of the wind turbine blade is the main governing factor that affects power generation from the wind turbine. Different airfoils, angle of twist and blade dimensions are the parameters that control the efficiency of the wind turbine. This study is aimed at investigating the aerodynamic performance of the wind turbine blade. In the present paper, we discuss innovative blade designs using the NACA 4412 airfoil, comparing them with a straight swept blade. The wake region was measured in the lab with a straight blade. All the results with different designs of blades were compared for their performance. A complete three-dimensional computational analysis was carried out to compare the power generation in each case for different wind speeds. It was found from the numerical analysis that the slotted blade yielded the most power generation among the other blade designs.


Author(s):  
Shane Colón ◽  
Mark Ricklick ◽  
Doug Nagy ◽  
Amy Lafleur

Abstract Thermal barrier coatings (TBC) found on turbine blades are a key element in the performance and reliability of modern gas turbines. TBC reduces the heat transfer into turbine blades by introducing an additional surface thermal resistance; consequently allowing for higher gas temperatures. During the service life of the blades, the TBC surface may be damaged due to manufacturing imperfections, handling damage, service spalling, or service impact damage, producing chips in the coating. While an increase in aerofoil temperature is expected, it is unknown to what degree the blade will be affected and what parameters of the chip shape affect this result. During routine inspections, the severity of the chipping will often fall to the discretion of the inspecting engineer. Without a quantitative understanding of the flow and heat transfer around these chips, there is potential for premature removal or possible blade failure if left to operate. The goal of this preliminary study is to identify the major driving parameters that lead to the increase in metal temperature when TBC is damaged, such that more quantitative estimates of blade life and refurbishing needs can be made. A two-dimensional computational Conjugate Heat Transfer model was developed; fully resolving the hot gas path and TBC, bond-coat, and super alloy solids. Representative convective conditions were applied to the cold side to emulate the characteristics of a cooled turbine blade. The hot gas path properties included an inlet temperature of 1600 K with varying Mach numbers of 0.30, 0.59, and 0.80 and Reynolds number of 5.1×105, 7.0×105, and 9.0×105 as referenced from the leading edge of the model. The cold side was given a coolant temperature of 750 K and a heat transfer coefficient of 1500 W/m2*K. The assigned thermal conductivities of the TBC, bond-coat, and metal alloys were 0.7 W/m*K, 7.0 W/m*K, and 11.0 W/m*K, respectively, and layer thicknesses of 0.50 mm, 0.25 mm, and 1.50 mm, respectively. A flat plate model without the presence of the chip was first evaluated to provide a basis of validation by comparison to existing correlations. Comparing heat transfer coefficients, the flat plate model matched within uncertainty to the Chilton-Colburn analogy. In addition, flat plate results captured the boundary layer thickness when compared with Prandtl’s 1/7th power-law. A chip was then introduced into the model, varying the chip width and the edge geometry. The most sensitive driving parameters were identified to be the chip width and Mach number. In cases where the chip width reached 16 times the TBC thickness, temperatures increased by almost 30% when compared to the undamaged equivalents. Additionally, increasing the Mach number of the incoming flow also increased metal temperatures. While the Reynolds number based on the leading edge of the model was deemed negligible, the Reynolds number based on the chip width was found to have a noticeable impact on the blade temperature. In conclusion, this study found that chip edge geometry was a negligible factor, while the Mach number, chip width, and Reynolds number based on the chip width had a significant effect on the total metal temperature.


Author(s):  
Karthik Krishnaswamy ◽  
◽  
Srikanth Salyan ◽  

The performance of a gas turbine during the service life can be enhanced by cooling the turbine blades efficiently. The objective of this study is to achieve high thermohydraulic performance (THP) inside a cooling passage of a turbine blade having aspect ratio (AR) 1:5 by using discrete W and V-shaped ribs. Hydraulic diameter (Dh) of the cooling passage is 50 mm. Ribs are positioned facing downstream with angle-of-attack (α) of 30° and 45° for discrete W-ribs and discerte V-ribs respectively. The rib profiles with rib height to hydraulic diameter ratio (e/Dh) or blockage ratio 0.06 and pitch (P) 36 mm are tested for Reynolds number (Re) range 30000-75000. Analysis reveals that, area averaged Nusselt numbers of the rib profiles are comparable, with maximum difference of 6% at Re 30000, which is within the limits of uncertainty. Variation of local heat transfer coefficients along the stream exhibited a saw tooth profile, with discrete W-ribs exhibiting higher variations. Along spanwise direction, discrete V-ribs showed larger variations. Maximum variation in local heat transfer coefficients is estimated to be 25%. For experimented Re range, friction loss for discrete W-ribs is higher than discrete-V ribs. Rib profiles exhibited superior heat transfer capabilities. The best Nu/Nuo achieved for discrete Vribs is 3.4 and discrete W-ribs is 3.6. In view of superior heat transfer capabilities, ribs can be deployed in cooling passages near the leading edge, where the temperatures are very high. The best THPo achieved is 3.2 for discrete V-ribs and 3 for discrete W-ribs at Re 30000. The ribs can also enhance the power-toweight ratio as they can produce high thermohydraulic performances for low blockage ratios.


2000 ◽  
Author(s):  
M. Singh ◽  
P. K. Panigrahi ◽  
G. Biswas

Abstract A numerical study of rib augmented cooling of turbine blades is reported in this paper. The time-dependent velocity field around a pair of symmetrically placed ribs on the walls of a three-dimensional rectangular channel was studied by use of a modified version of Marker-And-Cell algorithm to solve the unsteady incompressible Navier-Stokes and energy equations. The flow structures are presented with the help of instantaneous velocity vector and vorticity fields, FFT and time averaged and rms values of components of velocity. The spanwise averaged Nusselt number is found to increase at the locations of reattachment. The numerical results are compared with available numerical and experimental results. The presence of ribs leads to complex flow fields with regions of flow separation before and after the ribs. Each interruption in the flow field due to the surface mounted rib enables the velocity distribution to be more homogeneous and a new boundary layer starts developing downstream of the rib. The heat transfer is primarily enhanced due to the decrease in the thermal resistance owing to the thinner boundary layers on the interrupted surfaces. Another reason for heat transfer enhancement can be attributed to the mixing induced by large-scale structures present downstream of the separation point.


Sign in / Sign up

Export Citation Format

Share Document