D. Bohn, T. Uno, Takeshi Yoshida, Christian Betcher, Jan Frohnheiser, Kristof Weidtmann
One common approach for anti-erosion measures in low pressure steam turbines is to equip a hollow stator vane with slots on the airfoil surface in order to remove the water film by suction and consequently reduce the amount of secondary droplets. The purpose of this paper is to build an understanding of the predominant effects in fluid-film interaction and to examine the suitability of modern numerical methods for the design process of such slots. The performance of a suction slot in terms of collection rate and air leakage is investigated numerically in a flatplate setup with upstream injection of water. In order to model the relevant phenomena (film transport, edge stripping of droplets, transport of droplets in the surrounding fluid, wall impingement of droplets) an unsteady Eulerian-Lagrangian simulation setup is applied. The accuracy of the numerical approach is assessed by comparison with experimental measurements. The comparison of four cases with the measured data demonstrates that the chosen simulation approach is able to predict the main features of film flow and interaction with the surrounding fluid. The collection rate as well as fluid film properties show the same qualitative dependency from water mass flow rate and air velocity.
{"title":"Numerical and Experimental Study of Droplet-Film-Interaction for Low Pressure Steam Turbine Erosion Protection Applications","authors":"D. Bohn, T. Uno, Takeshi Yoshida, Christian Betcher, Jan Frohnheiser, Kristof Weidtmann","doi":"10.33737/gpps20-tc-65","DOIUrl":"https://doi.org/10.33737/gpps20-tc-65","url":null,"abstract":"One common approach for anti-erosion measures in low pressure steam turbines is to equip a hollow stator vane with slots on the airfoil surface in order to remove the water film by suction and consequently reduce the amount of secondary droplets. The purpose of this paper is to build an understanding of the predominant effects in fluid-film interaction and to examine the suitability of modern numerical methods for the design process of such slots. The performance of a suction slot in terms of collection rate and air leakage is investigated numerically in a flatplate setup with upstream injection of water. In order to model the relevant phenomena (film transport, edge stripping of droplets, transport of droplets in the surrounding fluid, wall impingement of droplets) an unsteady Eulerian-Lagrangian simulation setup is\u0000applied. The accuracy of the numerical approach is assessed by comparison with experimental\u0000measurements. The comparison of four cases with the measured data demonstrates that the chosen simulation approach is\u0000able to predict the main features of film flow and interaction with the surrounding fluid. The collection\u0000rate as well as fluid film properties show the same qualitative dependency from water mass flow rate and air velocity.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45347811","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The flow field in a one-stage axial flow turbine with 30 stator and 62 rotor blades including the wheel space is investigated by large-eddy simulation (LES). The Navier-Stokes equations are solved using a massively parallel finite-volume solver based on a Cartesian mesh with immersed boundaries. The strict conservation of mass, momentum, and energy is ensured by an efficient cut-cell/level-set ansatz, where a separate level-set solver describes the motion of the rotor. Both solvers use individual subsets of a shared Cartesian mesh, which they can adapt independently. The focus of the analysis is on the flow field inside the rotor stator cavity between the stator and rotor disks. Two cooling gas mass flow rates are investigated for the same rim seal geometry. First, the time averaged flow field for both simulations is compared, followed by a detailed investigation of the unsteady flow field. The results for the cooling effectiveness are compared to experimental data. Both cases show good agreement with experimental data. It is shown that for the lower cooling gas mass flux several of the wheel space’s acoustic waves are excited. This is not observed for the higher cooling gas mass flux. The excited waves lead to stable, i.e., bounded, fluctuations inside the wheel space and result in a significantly higher hot gas ingestion.
{"title":"Large-Eddy Simulations of Rim Seal Flow in a One-Stage Axial Turbine","authors":"Thomas Hösgen, M. Meinke, W. Schröder","doi":"10.33737/gpps20-tc-104","DOIUrl":"https://doi.org/10.33737/gpps20-tc-104","url":null,"abstract":"The flow field in a one-stage axial flow turbine with 30 stator and 62 rotor blades including the wheel space is investigated by large-eddy simulation (LES). The Navier-Stokes equations are solved using a massively parallel finite-volume solver based on a Cartesian mesh with immersed boundaries. The strict conservation of mass, momentum, and energy is ensured by an efficient cut-cell/level-set ansatz, where a separate level-set solver describes the motion of the rotor. Both solvers use individual subsets of a shared Cartesian mesh, which they can adapt independently. The focus of the analysis is on the flow field inside the rotor stator cavity between the stator and rotor disks. Two cooling gas mass flow rates are investigated for the same rim seal geometry. First, the time averaged flow field for both simulations is compared, followed by a detailed investigation of the unsteady flow field. The results for the cooling effectiveness are compared to experimental data. Both cases show good agreement with experimental data. It is shown that for the lower cooling gas mass flux several of the wheel space’s acoustic waves are excited. This is not observed for the higher cooling gas mass flux. The excited waves lead to stable, i.e., bounded, fluctuations inside the wheel space and result in a significantly higher hot gas ingestion.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41633623","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper presents Laser-Doppler Anemometry (LDA) measurements obtained from the Sussex Multiple Cavity test facility. This facility comprises a number of heated disc cavities with a cool bore flow and is intended to emulate the secondary air system flow in an H.P compressor. Measurements were made of the axial and tangential components of velocity over the respective range of Rossby, Rotational and Axial Reynolds numbers, (Ro, Reθ andRez),0.32<Ro<1.28,Reθ=7.1×105, 1.2×104<Rez<4.8×104 and for the values of the buoyancy parameter (βΔT) :0.50<βΔT<0.58. The frequency spectra analysis of the tangential ve
{"title":"Experimental and Computational Investigation of Flow Structure of Buoyancy Induced Flow in Heated Rotating Cavities","authors":"S. M. Fazeli, V. Kanjirakkad, C. Long","doi":"10.33737/gpps20-tc-37","DOIUrl":"https://doi.org/10.33737/gpps20-tc-37","url":null,"abstract":"This paper presents Laser-Doppler Anemometry (LDA) measurements obtained from the Sussex Multiple Cavity test facility. This facility comprises a number of heated disc cavities with a cool bore flow and is intended to emulate the secondary air system flow in an H.P compressor. Measurements were made of the axial and tangential components of velocity over the respective range of Rossby, Rotational and Axial Reynolds numbers, (Ro, <inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mtext>R</mml:mtext><mml:msub><mml:mtext>e</mml:mtext><mml:mi>θ</mml:mi></mml:msub></mml:math></inline-formula> and<inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mrow><mml:mspace width=\"0.25em\"/><mml:mi mathvariant=\"normal\">R</mml:mi></mml:mrow><mml:msub><mml:mtext>e</mml:mtext><mml:mi>z</mml:mi></mml:msub></mml:math></inline-formula>),<inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mrow><mml:mspace width=\"0.25em\"/></mml:mrow><mml:mn>0.32</mml:mn><mml:mo><</mml:mo><mml:mtext>Ro</mml:mtext><mml:mo><</mml:mo><mml:mn>1.28</mml:mn></mml:math></inline-formula>,<inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mrow><mml:mspace width=\"0.25em\"/><mml:mi mathvariant=\"normal\">R</mml:mi></mml:mrow><mml:msub><mml:mtext>e</mml:mtext><mml:mi>θ</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>7.1</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>5</mml:mn></mml:msup></mml:math></inline-formula>, <inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mn>1.2</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>4</mml:mn></mml:msup><mml:mo><</mml:mo><mml:mrow><mml:mspace width=\"0.25em\"/><mml:mi mathvariant=\"normal\">R</mml:mi></mml:mrow><mml:msub><mml:mtext>e</mml:mtext><mml:mi>z</mml:mi></mml:msub><mml:mo><</mml:mo><mml:mn>4.8</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>4</mml:mn></mml:msup></mml:math></inline-formula> and for the values of the buoyancy parameter <inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mo stretchy=\"false\">(</mml:mo><mml:mrow><mml:mi>β</mml:mi><mml:mrow><mml:mi mathvariant=\"normal\">Δ</mml:mi></mml:mrow><mml:mtext>T</mml:mtext></mml:mrow><mml:mo stretchy=\"false\">)</mml:mo></mml:math></inline-formula> :<inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mrow><mml:mspace width=\"0.25em\"/></mml:mrow><mml:mn>0.50</mml:mn><mml:mo><</mml:mo><mml:mrow><mml:mspace width=\"0.25em\"/><mml:mi>β</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant=\"normal\">Δ</mml:mi></mml:mrow><mml:mtext>T</mml:mtext><mml:mo><</mml:mo><mml:mn>0.58</mml:mn></mml:math></inline-formula>. The frequency spectra analysis of the tangential ve","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49229513","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The stall margin of tip-critical axial compressors can be improved by using circumferential casing grooves. From previous studies, in the literature, the stall margin improvement due to the casing grooves can be attributed to the reduction of the near casing blockage. The pressure rise across the compressor as the compressor is throttled intensifies the tip leakage flow. This results in a stronger tip leakage vortex that is thought to be the main source of the blockage. In this paper, the near casing blockage due to the tip region aerodynamics in a low-speed axial compressor rotor is numerically studied and quantified using a mass flow-based blockage parameter. The peak blockage location at the last stable operating point for a rotor with smooth casing is found to be at about 10% of the tip chord aft of the tip leading edge. Based on this information, an optimised single casing groove design that minimises the peak blockage is found using a surrogate-based optimisation approach. The implementation of the optimised groove is shown to produce a stall margin improvement of about 5%.
{"title":"Stall margin improvement in a low-speed axial compressor rotor using a blockage-optimised single circumferential casing groove","authors":"A. Mustaffa, V. Kanjirakkad","doi":"10.33737/gpps20-tc-89","DOIUrl":"https://doi.org/10.33737/gpps20-tc-89","url":null,"abstract":"The stall margin of tip-critical axial compressors can be improved by using circumferential casing grooves. From previous studies, in the literature, the stall margin improvement due to the casing grooves can be attributed to the reduction of the near casing blockage. The pressure rise across the compressor as the compressor is throttled intensifies the tip leakage flow. This results in a stronger tip leakage vortex that is thought to be the main source of the blockage. In this paper, the near casing blockage due to the tip region aerodynamics in a low-speed axial compressor rotor is numerically studied and quantified using a mass flow-based blockage parameter. The peak blockage location at the last stable operating point for a rotor with smooth casing is found to be at about 10% of the tip chord aft of the tip leading edge. Based on this information, an optimised single casing groove design that minimises the peak blockage is found using a surrogate-based optimisation approach. The implementation of the optimised groove is shown to produce a stall margin improvement of about 5%.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44965647","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Henrik Hoffmann, Lukas Stuhldreier, R. V. Rennings, P. Jeschke
This paper presents a numerical investigation of the effects of compressing various gases, for example, carbon dioxide (CO2) and methane (CH4), on an eight-stage axial air compressor. Several adaptation methods are applied to achieve a similar operating point as for air. Theoretically, the operating point depends on Mach number, flow angles, Reynolds number and isentropic exponent. Numerical results show mismatch effects which arise in the parameters using a non-adapted geometry. A rematching procedure is described, including deduced speed adjustments, in order to achieve Mach number equality at compressor inlet. Only shroud modifications are performed to rematch the flow angles of the air simulation. Although Reynolds and Mach number are kept constant at compressor inlet, an inevitable deviation in downstream flow causes mismatches in efficiency and pressure ratio. Both analytical and numerical methodologies show that the scale of shroud adjustments, as well as the size of mismatch in Mach and Reynolds number, can be correlated to the isentropic gas exponent. In summary, the main impact on gas behavior in an axial air compressor is attributable to the change in isentropic exponent. Derivations of shroud adaptation and analyses of inevitable aerodynamic mismatch are therefore developed depending on the isentropic exponent.
{"title":"Influence of Different Gases on the Design Point of an Industrial Axial Compressor and Deduced Aerodynamical Rematching Methodology","authors":"Henrik Hoffmann, Lukas Stuhldreier, R. V. Rennings, P. Jeschke","doi":"10.33737/gpps20-tc-84","DOIUrl":"https://doi.org/10.33737/gpps20-tc-84","url":null,"abstract":"This paper presents a numerical investigation of the effects of compressing various gases, for example, carbon dioxide (CO2) and methane (CH4), on an eight-stage axial air compressor. Several adaptation methods are applied to achieve a similar operating point as for air.\u0000Theoretically, the operating point depends on Mach number, flow angles, Reynolds number and isentropic exponent. Numerical results show mismatch effects which arise in the parameters using a non-adapted geometry. A rematching procedure is described, including deduced speed adjustments, in order to achieve Mach number equality at compressor inlet. Only shroud modifications are performed to rematch the flow angles of the air simulation. Although Reynolds and Mach number are kept constant at compressor inlet, an inevitable deviation in downstream flow causes mismatches in efficiency and pressure ratio. Both analytical and numerical methodologies show that the scale of shroud adjustments, as well as the size of mismatch in Mach and Reynolds number, can be correlated to the isentropic gas exponent.\u0000In summary, the main impact on gas behavior in an axial air compressor is attributable to the change in isentropic exponent. Derivations of shroud adaptation and analyses of inevitable aerodynamic mismatch are therefore developed depending on the isentropic exponent.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44866076","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Vahid Iranidokht, A. Kalfas, R. Abhari, Shigeki Senoo, Kazuhiro Momma
This paper presents an experimental investigation on the impact of different design and operational variations on the instabilities induced at the hub cavity outlet of a turbine. The experiments were conducted at the “LISA” test facility at ETH Zurich. The axial gap at the 2nd stage hub cavity exit was varied, and also three different flow deflectors were implemented at the cavity exit to control the cavity modes (CMs). Furthermore, the turbine pressure ratio was altered to mimic the off-design condition and study the sensitivity of the CMs to this parameter. Measurements were performed using pneumatic, and Fast Response Aerodynamic Probes (FRAP) at stator and rotor exit. In addition, unsteady pressure transducers were installed at the cavity exit wall to measure the characteristic parameters of the CMs. For the small axial gap, distinct and strong CMs were generated, which actively interacted with stator and rotor hub flow structures. Increasing the gap damped the fluctuations; however, a broader range of frequencies was amplified. The flow deflectors successfully suppressed the CMs by manipulating the shear layer velocity profile and blocking the growing instabilities. Eventually, the increase in the turbine pressure ratio strengthened the CMs and vice versa.
{"title":"Sensitivity analysis on the impact of geometrical and operational variations on turbine hub cavity modes and practical methods to control them","authors":"Vahid Iranidokht, A. Kalfas, R. Abhari, Shigeki Senoo, Kazuhiro Momma","doi":"10.33737/gpps20-tc-143","DOIUrl":"https://doi.org/10.33737/gpps20-tc-143","url":null,"abstract":"This paper presents an experimental investigation on the impact of different design and operational variations on the instabilities induced at the hub cavity outlet of a turbine. The experiments were conducted at the “LISA” test facility at ETH Zurich. The axial gap at the 2nd stage hub cavity exit was varied, and also three different flow deflectors were implemented at the cavity exit to control the cavity modes (CMs). Furthermore, the turbine pressure ratio was altered to mimic the off-design condition and study the sensitivity of the CMs to this parameter. Measurements were performed using pneumatic, and Fast Response Aerodynamic Probes (FRAP) at stator and rotor exit. In addition, unsteady pressure transducers were installed at the cavity exit wall to measure the characteristic parameters of the CMs.\u0000For the small axial gap, distinct and strong CMs were generated, which actively interacted with stator and rotor hub flow structures. Increasing the gap damped the fluctuations; however, a broader range of frequencies was amplified. The flow deflectors successfully suppressed the CMs by manipulating the shear layer velocity profile and blocking the growing instabilities. Eventually, the increase in the turbine pressure ratio strengthened the CMs and vice versa.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45816839","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Taofei Chen, Bijie Yang, Miles C. Robertson, R. Martinez-Botas
Real-gas effects have a significant impact on compressible turbulent flows of dense gases, especially when flow properties are in proximity of the saturation line and/or the thermodynamic critical point. Understanding of these effects is key for the analysis and improvement of performance for many industrial components, including expanders and heat exchangers in organic Rankine cycle systems. This work analyzes the real-gas effect on the turbulent boundary layer of fully developed channel flow of two organic gases, R1233zd(E) and MDM - two candidate working fluids for ORC systems. Compressible direct numerical simulations (DNS) with real-gas equations of state are used in this research. Three cases are set up for each organic vapour, representing thermodynamic states far from, close to and inside the supercritical region, and these cases refer to weak, normal and strong real-gas effect in each fluid. The results within this work show that the real-gas effect can significantly influence the profile of averaged thermodynamic properties, relative to an air baseline case. This effect has a reverse impact on the distribution of averaged temperature and density. As the real-gas effect gets stronger, the averaged centre-to-wall temperature ratio decreases but the density drop increases. In a strong real-gas effect case, the dynamic viscosity at the channel center point can be lower than at channel wall. This phenomenon can not be found in a perfect gas flow. The real-gas effect increases the normal Reynolds stress in the wall-normal direction by 7–20% and in the spanwise direction by 10–21%, which is caused by its impact on the viscosity profile. It also increases the Reynolds shear stress by 5–8%. The real-gas effect increases the turbulence kinetic energy dissipation in the viscous sublayer and buffer sublayer (y∗<30) but not in the outer layer. The turbulent viscosity hypthesis is checked in these two fluids, and the result shows that the standard two-function RANS model (k−ϵ and k−ω) with a constant Cμ=0.09 is still suitable in the outer layer
{"title":"Direct Numerical Simulation of real-gas effects within turbulent boundary layers for fully-developed channel flows","authors":"Taofei Chen, Bijie Yang, Miles C. Robertson, R. Martinez-Botas","doi":"10.33737/gpps20-tc-68","DOIUrl":"https://doi.org/10.33737/gpps20-tc-68","url":null,"abstract":"Real-gas effects have a significant impact on compressible turbulent flows of dense gases, especially when flow properties are in proximity of the saturation line and/or the thermodynamic critical point. Understanding of these effects is key for the analysis and improvement of performance for many industrial components, including expanders and heat exchangers in organic Rankine cycle systems.\u0000\u0000This work analyzes the real-gas effect on the turbulent boundary layer of fully developed channel flow of two organic gases, R1233zd(E) and MDM - two candidate working fluids for ORC systems. Compressible direct numerical simulations (DNS) with real-gas equations of state are used in this research. Three cases are set up for each organic vapour, representing thermodynamic states far from, close to and inside the supercritical region, and these cases refer to weak, normal and strong real-gas effect in each fluid.\u0000\u0000The results within this work show that the real-gas effect can significantly influence the profile of averaged thermodynamic properties, relative to an air baseline case. This effect has a reverse impact on the distribution of averaged temperature and density. As the real-gas effect gets stronger, the averaged centre-to-wall temperature ratio decreases but the density drop increases. In a strong real-gas effect case, the dynamic viscosity at the channel center point can be lower than at channel wall. This phenomenon can not be found in a perfect gas flow.\u0000\u0000The real-gas effect increases the normal Reynolds stress in the wall-normal direction by 7–20% and in the spanwise direction by 10–21%, which is caused by its impact on the viscosity profile. It also increases the Reynolds shear stress by 5–8%. The real-gas effect increases the turbulence kinetic energy dissipation in the viscous sublayer and buffer sublayer <inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mo stretchy=\"false\">(</mml:mo><mml:msup><mml:mi>y</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mo><</mml:mo><mml:mn>30</mml:mn><mml:mo stretchy=\"false\">)</mml:mo></mml:math></inline-formula> but not in the outer layer. The turbulent viscosity hypthesis is checked in these two fluids, and the result shows that the standard two-function RANS model (<inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>ϵ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>ω</mml:mi></mml:math></inline-formula>) with a constant <inline-formula><mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\" display=\"inline\" overflow=\"scroll\"><mml:msub><mml:mi>C</mml:mi><mml:mi>μ</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>0.09</mml:mn></mml:math></inline-formula> is still suitable in the outer layer <inline-formula><mml:math xm","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45465153","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
With the imminent digitalisation of the manufacturing processes of gas turbine components, a large volume of geometric data of as-manufactured parts is being generated. This geometric data can be used in aerodynamic simulations to predict component performance. Both the cost and accuracy of these simulations increase with their fidelity. To efficiently exploit Digital Twin technology, one must therefore understand how realistic the aerodynamic simulations need to be to give useful performance predictions. This paper considers this issue for a sample of scrapped high-pressure turbine rotor blades from a civil aero engine. The measured geometric data was used to build aerodynamic models of varying degrees of realism, ranging from quasi-three-dimensional blade sections for an Euler solver to three-dimensional, multi-passage and multi-stage Reynolds-Averaged-Navier-Stokes models. The flow near the tip of these shrouded blades is sensitive to manufacturing variability and can switch between two quasi-stable horseshoe vortex modes. In general, capacity and exit flow angle can be adequately predicted by three-dimensional, single-passage calculations: averaging single-passage calculations gives a good prediction of the multi-passage behaviour. For efficiency and stage loading, the approach of averaging single-passage calculations is less accurate as the multi-passage behaviour requires an accurate prediction of the horseshoe vortex modes.
{"title":"The required aerodynamic simulation fidelity to usefully support a gas turbine Digital Twin for manufacturing","authors":"Wen Yao Lee, W. Dawes, J. Coull","doi":"10.33737/gpps20-tc-35","DOIUrl":"https://doi.org/10.33737/gpps20-tc-35","url":null,"abstract":"With the imminent digitalisation of the manufacturing processes of gas turbine components, a large volume of geometric data of as-manufactured parts is being generated. This geometric data can be used in aerodynamic simulations to predict component performance. Both the cost and accuracy of these simulations increase with their fidelity. To efficiently exploit Digital Twin technology, one must therefore understand how realistic the aerodynamic simulations need to be to give useful performance predictions. \u0000This paper considers this issue for a sample of scrapped high-pressure turbine rotor blades from a civil aero engine. The measured geometric data was used to build aerodynamic models of varying degrees of realism, ranging from quasi-three-dimensional blade sections for an Euler solver to three-dimensional, multi-passage and multi-stage Reynolds-Averaged-Navier-Stokes models. The flow near the tip of these shrouded blades is sensitive to manufacturing variability and can switch between two quasi-stable horseshoe vortex modes. In general, capacity and exit flow angle can be adequately predicted by three-dimensional, single-passage calculations: averaging single-passage calculations gives a good prediction of the multi-passage behaviour. For efficiency and stage loading, the approach of averaging single-passage calculations is less accurate as the multi-passage behaviour requires an accurate prediction of the horseshoe vortex modes.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47889400","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
An analytical model to describe the flutter onset of straight-through labyrinth seals has been numerically validated using a frequency domain linearized Navier-Stokes solver. A comprehensive set of simulations has been conducted to assess the stability criterion of the analytical model originally derived by Corral and Vega (2018), “Conceptual Flutter Analysis of Labyrinth Seals Using Analytical Models - Part I: Theoretical Support,” ASME J. Turbomach., 140 (12), pp. 121006. The accuracy of the model has been assessed by using a simplified geometry consisting of a two-fin straight-through labyrinth seal with identical gaps. The effective gaps and the kinetic energy carried over are retained and their effects on stability are evaluated. It turns out that is important to inform the model with the correct values of both parameters to allow a proper comparison with the numerical simulations. Moreover, the non-isentropic perturbations included in the formulations are observed in the simulations at relatively low frequencies whose characteristic time is of the same order as the discharge time of the seal. This effect is responsible for the bending of the stability limit in the 0th ND stability map obtained both in the model and the simulations. It turns out that the analytical model can predict accurately the stability of the seal in a wide range of pressure ratios, vibration mode-shapes, and frequencies provided that this is informed with the fluid dynamic gaps and the energy carried over to the downstream fin from a steady RANS simulation. The numerical calculations show for the first time that the model can be used to predict accurately not only the trends of the work-per-cycle of the seal but also quantitative results.
{"title":"Numerical Validation of an Analytical Seal Flutter Model","authors":"M. Greco, R. Corral","doi":"10.33737/gpps20-tc-142","DOIUrl":"https://doi.org/10.33737/gpps20-tc-142","url":null,"abstract":"An analytical model to describe the flutter onset of straight-through labyrinth seals has been numerically validated using a frequency domain linearized Navier-Stokes solver. A comprehensive set of simulations has been conducted to assess the stability criterion of the analytical model originally derived by Corral and Vega (2018), “Conceptual Flutter Analysis of Labyrinth Seals Using Analytical Models - Part I: Theoretical Support,” ASME J. Turbomach., 140 (12), pp. 121006. The accuracy of the model has been assessed by using a simplified geometry consisting of a two-fin straight-through labyrinth seal with identical gaps. The effective gaps and the kinetic energy carried over are retained and their effects on stability are evaluated. It turns out that is important to inform the model with the correct values of both parameters to allow a proper comparison with the numerical simulations. Moreover, the non-isentropic perturbations included in the formulations are observed in the simulations at relatively low frequencies whose characteristic time is of the same order as the discharge time of the seal. This effect is responsible for the bending of the stability limit in the 0th ND stability map obtained both in the model and the simulations. It turns out that the analytical model can predict accurately the stability of the seal in a wide range of pressure ratios, vibration mode-shapes, and frequencies provided that this is informed with the fluid dynamic gaps and the energy carried over to the downstream fin from a steady RANS simulation. The numerical calculations show for the first time that the model can be used to predict accurately not only the trends of the work-per-cycle of the seal but also quantitative results.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42567989","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The generation mechanism of a diffuser stall in a centrifugal compressor with a vaneless diffuser was investigated by experimental and computational analyses. The diffuser stall generated as the mass flow rate decreased. The diffuser stall cell rotated at 25-30 % of the impeller rotational speed, with diffuser stall fluctuations observed at 180° from the cutoff. The diffuser stall fluctuation magnitude gradually increased near the cutoff. According to the CFD analysis, the mass flow fluctuations at the diffuser exit showed a low mass flow region, rotating at approximately 25% of the impeller rotational speed. They began at 180° from the cutoff and developed as this region approached the cutoff. Therefore, the diffuser stall could be simulated by CFD analysis. First, the diffuser stall cell originated at 180° from the cutoff by interaction with boundary separation and impeller discharge vortex. Then, the diffuser stall cell further developed by boundary separation accumulation and the induced low velocity area The low velocity region formed a blockage across the diffuser passage span. The diffuser stall cell expanded due to boundary separation caused by a positive flow angle. Finally, the diffuser stall cell vanished when it passed the cutoff, because mass flow recovery occurred.
{"title":"Generation Mechanism of Diffuser Stall in a Centrifugal Compressor with Vaneless Diffuser","authors":"N. Fujisawa, K. Tajima, H. Miida, Y. Ohta","doi":"10.33737/gpps20-tc-59","DOIUrl":"https://doi.org/10.33737/gpps20-tc-59","url":null,"abstract":"The generation mechanism of a diffuser stall in a centrifugal compressor with a vaneless diffuser was investigated by experimental and computational analyses. The diffuser stall generated as the mass flow rate decreased. The diffuser stall cell rotated at 25-30 % of the impeller rotational speed, with diffuser stall fluctuations observed at 180° from the cutoff. The diffuser stall fluctuation magnitude gradually increased near the cutoff. According to the CFD analysis, the mass flow fluctuations at the diffuser exit showed a low mass flow region, rotating at approximately 25% of the impeller rotational speed. They began at 180° from the cutoff and developed as this region approached the cutoff. Therefore, the diffuser stall could be simulated by CFD analysis. First, the diffuser stall cell originated at 180° from the cutoff by interaction with boundary separation and impeller discharge vortex. Then, the diffuser stall cell further developed by boundary separation accumulation and the induced low velocity area The low velocity region formed a blockage across the diffuser passage span. The diffuser stall cell expanded due to boundary separation caused by a positive flow angle. Finally, the diffuser stall cell vanished when it passed the cutoff, because mass flow recovery occurred.","PeriodicalId":53002,"journal":{"name":"Journal of the Global Power and Propulsion Society","volume":null,"pages":null},"PeriodicalIF":0.9,"publicationDate":"2020-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43524286","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：