UFR 4-05 Test Case: Difference between revisions
m (Dave.Ellacott moved page SilverP:UFR 4-05 Test Case to UFR 4-05 Test Case over redirect) |
|||
| (16 intermediate revisions by the same user not shown) | |||
| Line 17: | Line 17: | ||
The DLR annular test cascade is a subsonic, low aspect ratio turbine stator. The stator is scaled by the factor of 2.348 and has constant hub and shroud radii. | The DLR annular test cascade is a subsonic, low aspect ratio turbine stator. The stator is scaled by the factor of 2.348 and has constant hub and shroud radii. | ||
Figure 2: shows a meridional view of the geometry. | Figure 2: shows a meridional view of the geometry. | ||
| Line 28: | Line 27: | ||
Figure 3: shows a 3D view of the blade geometry. The 3D blade coordinates can be obtained at the Institut für Antriebstechnik (Germany) (see AGARD-AR-355 for more details). | Figure 3: shows a 3D view of the blade geometry. The 3D blade coordinates can be obtained at the Institut für Antriebstechnik (Germany) (see AGARD-AR-355 for more details). | ||
Figure 3: 3D view of the geometry | <center>[[Image:U4-05d32_files_image005.jpg]]</center> | ||
<center>'''Figure 3:''' 3D view of the geometry</center> | |||
The following table gives the main dimensions of the stator. | The following table gives the main dimensions of the stator. | ||
| Line 69: | Line 70: | ||
<center>20.5°</center> | <center>20.5°</center> | ||
|} | |} | ||
== Test Case Experiments == | == Test Case Experiments == | ||
| Line 74: | Line 76: | ||
The test facility is shown in Figure 4:. It is an open loop continuously operating facility supplied with cold air. | The test facility is shown in Figure 4:. It is an open loop continuously operating facility supplied with cold air. | ||
Figure 4: Test installation (From AGARD-AR-355. « CFD Validation for Propulsion System Components», Publisher: Advisory Group for Aerospace Research & Development). | <center>[[Image:U4-05d32_files_image007.gif]]</center> | ||
<center>'''Figure 4:''' Test installation (From AGARD-AR-355. « CFD Validation for Propulsion System Components», Publisher: Advisory Group for Aerospace Research & Development).</center> | |||
The experimental conditions are summarized in the following table. | The experimental conditions are summarized in the following table. | ||
| Line 135: | Line 139: | ||
|} | |} | ||
The measurements were carried out on 4 planes at constant axial coordinate, one in front of the stator, 2 inside the blade passage and one downstream. These | The measurements were performed with a 5-hole probe downstream of the stator and with the three dimensional Laser-Two-Focus (3D-L2F) measuring technique inside the passage and downstream. Additionally, the upstream boundary layers were traversed using a ‘fish-mouth’ Pitot probe. | ||
The measurements were carried out on 4 planes at constant axial coordinate, one in front of the stator, 2 inside the blade passage and one downstream. These planes were located, respectively at −1.02, 0.8, 0.93 and 1.4 midspan axial chord from the leading edge. Only results on the last measuring plane are discussed here. | |||
The estimated uncertainty of the 5-hole probe data is: | The estimated uncertainty of the 5-hole probe data is: | ||
| Line 155: | Line 160: | ||
Two series of simulations on the DLR turbine stator are used in the present report. Both encompass a high number of simulations with various turbulence models. The first series came from the AGARD Advisory Report 355 whereas the second one came from new calculations performed with the FINE EURANUS/Turbo flow solver. These two series represent a total number of 29 simulations performed on this test case. | Two series of simulations on the DLR turbine stator are used in the present report. Both encompass a high number of simulations with various turbulence models. The first series came from the AGARD Advisory Report 355 whereas the second one came from new calculations performed with the FINE EURANUS/Turbo flow solver. These two series represent a total number of 29 simulations performed on this test case. | ||
The AGARD study | ===The AGARD study=== | ||
The aim of the working group involved in the AGARD Advisory Report 355 was to assess the ability of CFD to predict several of the aerodynamic phenomena, which dominate the performance of turbomachines. | The aim of the working group involved in the AGARD Advisory Report 355 was to assess the ability of CFD to predict several of the aerodynamic phenomena, which dominate the performance of turbomachines. | ||
| Line 161: | Line 166: | ||
Numerical results of ten different flow solvers were presented in the AGARD report. Here follows a non-exhaustive list of these codes: TRANSCode, CANARI, ADPAC, RVC3D, TIP3D, TRACE-S, HAH3D, and FLOS3D. 19 simulations were performed on this test case. Among them, 8 were using an algebraic/mixing length turbulence closure model, whereas 11 calculations were using a turbulent transport model. No common grid was imposed and some participants have performed calculation with the same model on different grids (for instance TRACE-S with a ''k-ε'' turbulence model was used both with a CH and a OHHH grid). | Numerical results of ten different flow solvers were presented in the AGARD report. Here follows a non-exhaustive list of these codes: TRANSCode, CANARI, ADPAC, RVC3D, TIP3D, TRACE-S, HAH3D, and FLOS3D. 19 simulations were performed on this test case. Among them, 8 were using an algebraic/mixing length turbulence closure model, whereas 11 calculations were using a turbulent transport model. No common grid was imposed and some participants have performed calculation with the same model on different grids (for instance TRACE-S with a ''k-ε'' turbulence model was used both with a CH and a OHHH grid). | ||
It should be stressed that those calculations were performed in the early | It should be stressed that those calculations were performed in the early 90’s and that those simulations are representative of the state of the art in CFD at this time. | ||
The inlet boundary condition is based on the data measured on the first plane. A profile of total pressure is imposed together with constant total temperature and angles. The total temperature is set to 306.6 K and the velocity direction is put to axial. Figure 5: shows the total pressure profile at the inlet. It must be stressed that there is uncertainties on the flow angle and velocity distributions which are were not measured at the inlet and which derived from the total pressure profile and the extrapolated pressure for the numerical experiments. At the downstream boundary the hub pressure if imposed together with a radial equilibrium. | The inlet boundary condition is based on the data measured on the first plane. A profile of total pressure is imposed together with constant total temperature and angles. The total temperature is set to 306.6 K and the velocity direction is put to axial. Figure 5: shows the total pressure profile at the inlet. It must be stressed that there is uncertainties on the flow angle and velocity distributions which are were not measured at the inlet and which derived from the total pressure profile and the extrapolated pressure for the numerical experiments. At the downstream boundary the hub pressure if imposed together with a radial equilibrium. | ||
<center>[[Image:U4-05d32_files_image008.gif]]</center> | |||
New calculations | |||
<center>'''Figure 5:''' Inlet Total Pressure Profile (From AGARD-AR-355. « CFD Validation for Propulsion System Components», Publisher: Advisory Group for Aerospace Research & Development).</center> | |||
===New calculations=== | |||
Ten new calculations were performed with the Code FINE(tm)/Turbo from NUMECA Int. | Ten new calculations were performed with the Code FINE(tm)/Turbo from NUMECA Int. | ||
| Line 176: | Line 183: | ||
An H-O-H non-matching grid, generated with IGG/Autogrid, is used. Two grids were generated with the same number of nodes but with different cell clustering at the solid walls: 1×10<sup>-6</sup> and 1×10<sup>-5</sup> meter corresponding to y<sup>+</sup> = 0.6-0.7 and y<sup>+</sup> = 6-7, respectively. There are three blocks in the H-O-H grid, related to the inlet block, principal block, outlet block respectively, with 29×81×21 29×81×221 and 33×81×49 grid points (azimuthal × spanwise × axial). Figure 6: shows the mesh at midspan and a meridional view. | An H-O-H non-matching grid, generated with IGG/Autogrid, is used. Two grids were generated with the same number of nodes but with different cell clustering at the solid walls: 1×10<sup>-6</sup> and 1×10<sup>-5</sup> meter corresponding to y<sup>+</sup> = 0.6-0.7 and y<sup>+</sup> = 6-7, respectively. There are three blocks in the H-O-H grid, related to the inlet block, principal block, outlet block respectively, with 29×81×21 29×81×221 and 33×81×49 grid points (azimuthal × spanwise × axial). Figure 6: shows the mesh at midspan and a meridional view. | ||
Figure 6: Mesh | <center>[[Image:U4-05d32_files_image010.gif]]</center> | ||
<center>'''Figure 6:''' Mesh</center> | |||
Different turbulence models were applied: the Baldwin-Lomax model and four k-epsilon models: the low Reynolds Chien model, the low Reynolds Yang-Shih model, the low Reynolds Launder-Sharma model and the standard model (based on the Yang-Shih model together with an extended wall functions derived from DNS curve fitting in order to work for values of y<sup>+</sup> greater than 20). For some turbulence models, two sets of boundary conditions for the turbulence quantities were used. The boundary condition for the turbulent kinetic energy ''k'' is set to match the turbulence intensity of the measurements (4.4%) to a constant value of 11.2 m<sup>2</sup>/s<sup>2</sup>. The turbulence dissipation ε is set to impose a turbulence ratio [[Image:U4-05d32_files_image012.gif]] equal to 5 and 50. The corresponding values of ε are 232000 m<sup>2</sup>/s<sup>3</sup> and 23200 m<sup>2</sup>/s<sup>3</sup>, respectively. | Different turbulence models were applied: the Baldwin-Lomax model and four k-epsilon models: the low Reynolds Chien model, the low Reynolds Yang-Shih model, the low Reynolds Launder-Sharma model and the standard model (based on the Yang-Shih model together with an extended wall functions derived from DNS curve fitting in order to work for values of y<sup>+</sup> greater than 20). For some turbulence models, two sets of boundary conditions for the turbulence quantities were used. The boundary condition for the turbulent kinetic energy ''k'' is set to match the turbulence intensity of the measurements (4.4%) to a constant value of 11.2 m<sup>2</sup>/s<sup>2</sup>. The turbulence dissipation ε is set to impose a turbulence ratio [[Image:U4-05d32_files_image012.gif]] equal to 5 and 50. The corresponding values of ε are 232000 m<sup>2</sup>/s<sup>3</sup> and 23200 m<sup>2</sup>/s<sup>3</sup>, respectively. | ||
Latest revision as of 14:09, 12 February 2017
Curved passage flow
Underlying Flow Regime 4-05 © copyright ERCOFTAC 2004
Test Case
Brief description of the study test case
The DLR annular test cascade is a subsonic, low aspect ratio turbine stator. The stator is scaled by the factor of 2.348 and has constant hub and shroud radii. Figure 2: shows a meridional view of the geometry.
Figure 3: shows a 3D view of the blade geometry. The 3D blade coordinates can be obtained at the Institut für Antriebstechnik (Germany) (see AGARD-AR-355 for more details).
The following table gives the main dimensions of the stator.
|
Blade Geometry | |
|
Number of blades |
|
|
Chord length at hub, CH |
|
|
Chord length at midspan, CMS |
|
|
Chord length at tip, CT |
|
|
Axial chord length at midspan, CAX, MS |
|
|
Outlet flow angle reference to tangential direction, α 3 |
|
Test Case Experiments
The test facility is shown in Figure 4:. It is an open loop continuously operating facility supplied with cold air.
The experimental conditions are summarized in the following table.
|
Experimental conditions | |
|
Mass flow rate, m0 |
|
|
Total pressure, Pto,MS |
|
|
Total temperature, Tt0 |
|
|
Inlet flow angle, α0 (circumf.) |
|
|
Inlet flow angle, β0(Radial) |
|
|
Inlet turbulence intensity level, Tu0 |
|
|
Inlet Mach number, Ma0 |
|
|
Mean outlet Mach number, Ma3 |
|
|
Static pressure at hub, Pst,3 |
|
|
Reynolds number, Re |
|
The measurements were performed with a 5-hole probe downstream of the stator and with the three dimensional Laser-Two-Focus (3D-L2F) measuring technique inside the passage and downstream. Additionally, the upstream boundary layers were traversed using a ‘fish-mouth’ Pitot probe.
The measurements were carried out on 4 planes at constant axial coordinate, one in front of the stator, 2 inside the blade passage and one downstream. These planes were located, respectively at −1.02, 0.8, 0.93 and 1.4 midspan axial chord from the leading edge. Only results on the last measuring plane are discussed here.
The estimated uncertainty of the 5-hole probe data is:
- Flow angles: ±0.5 deg.
- Pressures (static and total): ±0.1 mbar
The uncertainty for the 3D Two-Focus system is:
- Velocity: ±1.5 m/s
- Flow angle: ± 0.1 to 0.2 deg.
CFD Methods
Two series of simulations on the DLR turbine stator are used in the present report. Both encompass a high number of simulations with various turbulence models. The first series came from the AGARD Advisory Report 355 whereas the second one came from new calculations performed with the FINE EURANUS/Turbo flow solver. These two series represent a total number of 29 simulations performed on this test case.
The AGARD study
The aim of the working group involved in the AGARD Advisory Report 355 was to assess the ability of CFD to predict several of the aerodynamic phenomena, which dominate the performance of turbomachines.
Numerical results of ten different flow solvers were presented in the AGARD report. Here follows a non-exhaustive list of these codes: TRANSCode, CANARI, ADPAC, RVC3D, TIP3D, TRACE-S, HAH3D, and FLOS3D. 19 simulations were performed on this test case. Among them, 8 were using an algebraic/mixing length turbulence closure model, whereas 11 calculations were using a turbulent transport model. No common grid was imposed and some participants have performed calculation with the same model on different grids (for instance TRACE-S with a k-ε turbulence model was used both with a CH and a OHHH grid).
It should be stressed that those calculations were performed in the early 90’s and that those simulations are representative of the state of the art in CFD at this time.
The inlet boundary condition is based on the data measured on the first plane. A profile of total pressure is imposed together with constant total temperature and angles. The total temperature is set to 306.6 K and the velocity direction is put to axial. Figure 5: shows the total pressure profile at the inlet. It must be stressed that there is uncertainties on the flow angle and velocity distributions which are were not measured at the inlet and which derived from the total pressure profile and the extrapolated pressure for the numerical experiments. At the downstream boundary the hub pressure if imposed together with a radial equilibrium.
New calculations
Ten new calculations were performed with the Code FINE(tm)/Turbo from NUMECA Int. flow solver (http://www.numeca.com). In this code, the Reynolds-averaged Navier-Stokes equations are solved with finite volume method. Central-difference is employed for the spatial discretization. A Four stage explicit Runge-Kutta scheme is applied for the temporal discretization. The multi-grid method is used combined with local time step and implicit residual smoothing to speed up the convergence.
An H-O-H non-matching grid, generated with IGG/Autogrid, is used. Two grids were generated with the same number of nodes but with different cell clustering at the solid walls: 1×10-6 and 1×10-5 meter corresponding to y+ = 0.6-0.7 and y+ = 6-7, respectively. There are three blocks in the H-O-H grid, related to the inlet block, principal block, outlet block respectively, with 29×81×21 29×81×221 and 33×81×49 grid points (azimuthal × spanwise × axial). Figure 6: shows the mesh at midspan and a meridional view.
Different turbulence models were applied: the Baldwin-Lomax model and four k-epsilon models: the low Reynolds Chien model, the low Reynolds Yang-Shih model, the low Reynolds Launder-Sharma model and the standard model (based on the Yang-Shih model together with an extended wall functions derived from DNS curve fitting in order to work for values of y+ greater than 20). For some turbulence models, two sets of boundary conditions for the turbulence quantities were used. The boundary condition for the turbulent kinetic energy k is set to match the turbulence intensity of the measurements (4.4%) to a constant value of 11.2 m2/s2. The turbulence dissipation ε is set to impose a turbulence ratio 
equal to 5 and 50. The corresponding values of ε are 232000 m2/s3 and 23200 m2/s3, respectively.
Inlet boundary conditions are imposed as for the AGARD study. A radial equilibrium boundary condition is imposed at the outlet. The mean measured outlet pressure is imposed at midspan.
© copyright ERCOFTAC 2004
Contributors: Nouredine Hakimi - NUMECA International