Difference between revisions of "UFR 405 Test Case"
m (Dave.Ellacott moved page SilverP:UFR 405 Test Case to UFR 405 Test Case over redirect) 

(9 intermediate revisions by the same user not shown)  
Line 70:  Line 70:  
<center>20.5°</center>  <center>20.5°</center>  
}  }  
+  
== Test Case Experiments ==  == Test Case Experiments ==  
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.  
+  
<center>[[Image:U405d32_files_image007.gif]]</center>  <center>[[Image:U405d32_files_image007.gif]]</center>  
−  Figure 4: Test installation (From AGARDAR355. « CFD Validation for Propulsion System Components», Publisher: Advisory Group for Aerospace Research & Development).  +  <center>'''Figure 4:''' Test installation (From AGARDAR355. « 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 136:  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 5hole probe downstream of the stator and with the three dimensional LaserTwoFocus (3DL2F) measuring technique inside the passage and downstream. Additionally, the upstream boundary layers were traversed using a ‘fishmouth’ 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 5hole probe data is:  The estimated uncertainty of the 5hole probe data is:  
Line 156:  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 162:  Line 166:  
Numerical results of ten different flow solvers were presented in the AGARD report. Here follows a nonexhaustive list of these codes: TRANSCode, CANARI, ADPAC, RVC3D, TIP3D, TRACES, 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 TRACES 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 nonexhaustive list of these codes: TRANSCode, CANARI, ADPAC, RVC3D, TIP3D, TRACES, 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 TRACES 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.  
−  
−  Figure 5: Inlet Total Pressure Profile (From AGARDAR355. « CFD Validation for Propulsion System Components», Publisher: Advisory Group for Aerospace Research & Development).  +  <center>[[Image:U405d32_files_image008.gif]]</center> 
+  
+  
+  <center>'''Figure 5:''' Inlet Total Pressure Profile (From AGARDAR355. « CFD Validation for Propulsion System Components», Publisher: Advisory Group for Aerospace Research & Development).</center>  
−  New calculations  +  ===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 177:  Line 183:  
An HOH nonmatching 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.60.7 and y<sup>+</sup> = 67, respectively. There are three blocks in the HOH 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 HOH nonmatching 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.60.7 and y<sup>+</sup> = 67, respectively. There are three blocks in the HOH 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:U405d32_files_image010.gif]]</center> 
+  
+  <center>'''Figure 6:''' Mesh</center>  
+  
Different turbulence models were applied: the BaldwinLomax model and four kepsilon models: the low Reynolds Chien model, the low Reynolds YangShih model, the low Reynolds LaunderSharma model and the standard model (based on the YangShih 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:U405d32_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 BaldwinLomax model and four kepsilon models: the low Reynolds Chien model, the low Reynolds YangShih model, the low Reynolds LaunderSharma model and the standard model (based on the YangShih 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:U405d32_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
Contents
Curved passage flow
Underlying Flow Regime 405 © 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 AGARDAR355 for more details).
The following table gives the main dimensions of the stator.
Blade Geometry  
Number of blades 

Chord length at hub, C_{H} 

Chord length at midspan, C_{MS} 

Chord length at tip, C_{T} 

Axial chord length at midspan, C_{AX, 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, m_{0} 

Total pressure, P_{to,MS} 

Total temperature, T_{t0} 

Inlet flow angle, α_{0 }(circumf.) 

Inlet flow angle, β_{0}(Radial) 

Inlet turbulence intensity level, Tu_{0} 

Inlet Mach number, Ma_{0} 

Mean outlet Mach number, Ma_{3} 

Static pressure at hub, P_{st,3} 

Reynolds number, Re 

The measurements were performed with a 5hole probe downstream of the stator and with the three dimensional LaserTwoFocus (3DL2F) measuring technique inside the passage and downstream. Additionally, the upstream boundary layers were traversed using a ‘fishmouth’ 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 5hole probe data is:
 Flow angles: ±0.5 deg.
 Pressures (static and total): ±0.1 mbar
The uncertainty for the 3D TwoFocus 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 nonexhaustive list of these codes: TRANSCode, CANARI, ADPAC, RVC3D, TIP3D, TRACES, 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 TRACES 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 Reynoldsaveraged NavierStokes equations are solved with finite volume method. Centraldifference is employed for the spatial discretization. A Four stage explicit RungeKutta scheme is applied for the temporal discretization. The multigrid method is used combined with local time step and implicit residual smoothing to speed up the convergence.
An HOH nonmatching 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.60.7 and y^{+} = 67, respectively. There are three blocks in the HOH 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 BaldwinLomax model and four kepsilon models: the low Reynolds Chien model, the low Reynolds YangShih model, the low Reynolds LaunderSharma model and the standard model (based on the YangShih 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 m^{2}/s^{2}. The turbulence dissipation ε is set to impose a turbulence ratio equal to 5 and 50. The corresponding values of ε are 232000 m^{2}/s^{3} and 23200 m^{2}/s^{3}, 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