UFR 2-05 Test Case: Difference between revisions
David.Fowler (talk | contribs) No edit summary |
David.Fowler (talk | contribs) mNo edit summary |
||
Line 24: | Line 24: | ||
<center>From Large Eddy Simulation of the Flow Around a High Lift Aerofoil by Simon Dahlström, PhD Thesis, Chalmers University Sweden (2003) p9. Reproduced with permission. </center> | <center>From Large Eddy Simulation of the Flow Around a High Lift Aerofoil by Simon Dahlström, PhD Thesis, Chalmers University Sweden (2003) p9. Reproduced with permission. </center> | ||
The case had its origin in a French co-operative research project Opération Décrochage driven by Aerospatiale and including ONERA, CERT and IMFT as partners. A complete experimental database was compiled (see below) covering transition, boundary layer development and the separated region. It was noted even at that time that algebraic and k-epsilon models failed to predict the separation [U2-05ref.htm#ref6 [6]]. | The case had its origin in a French co-operative research project Opération Décrochage driven by Aerospatiale and including ONERA, CERT and IMFT as partners. A complete experimental database was compiled (see below) covering transition, boundary layer development and the separated region. It was noted even at that time that algebraic and k-epsilon models failed to predict the separation [U2-05ref.htm#ref6 [6]]. | ||
Line 34: | Line 33: | ||
<center>[[Image:U2-05d32_files_image004.jpg]] Figure 2. Lift and drag coefficients for Re<sub>c</sub><nowiki>=2.0x10</nowiki><sup>6 </sup>(circles); Re<sub>c</sub><nowiki>=2.07x10</nowiki><sup>6</sup>(squares);</center> | <center>[[Image:U2-05d32_files_image004.jpg]] Figure 2. Lift and drag coefficients for Re<sub>c</sub><nowiki>=2.0x10</nowiki><sup>6 </sup>(circles); Re<sub>c</sub><nowiki>=2.07x10</nowiki><sup>6</sup>(squares);</center> | ||
<center>Re<sub>c</sub | <center>Re<sub>c</sub>=3.13x10</nowiki><sup>6</sup> (crosses). The vertical dotted line is at α<nowiki>=13.3°.</center> | ||
<center>From Large Eddy Simulation of the Flow Around a High Lift Aerofoil by Simon Dahlström, </center> | <center>From Large Eddy Simulation of the Flow Around a High Lift Aerofoil by Simon Dahlström, </center> | ||
Line 40: | Line 39: | ||
<center>PhD Thesis, Chalmers University Sweden (2003) p10. Reproduced with permission.</center> | <center>PhD Thesis, Chalmers University Sweden (2003) p10. Reproduced with permission.</center> | ||
The results obtained in the F2 tunnel comprise pressure coefficient and skin friction distributions and LDV measurements of velocity and Reynolds stress profiles in the boundary layer and wake. From the latter, | The results obtained in the F2 tunnel comprise pressure coefficient and skin friction distributions and LDV measurements of velocity and Reynolds stress profiles in the boundary layer and wake. From the latter, δ, δ* and θ distributions were obtained. The laminar bubble leading to transition at 12% chord was identified using oil-flow visualisation. The flow was confirmed to be approximately two-dimensional up to about 13°, but at 13.3° some three-dimensionality may have been beginning to appear, leading to lower C<sub>L</sub> and higher C<sub>D</sub> compared to measurements in the F1 tunnel. | ||
== CFD Methods == | == CFD Methods == | ||
Line 83: | Line 82: | ||
BAe | BAe | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
Baldwin-Lomax, k- | Baldwin-Lomax, k-ε / Wolfshtein | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Line 92: | Line 91: | ||
CERFACS | CERFACS | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
k- | k-ε / Chien, Algebraic Stress Model | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Line 110: | Line 109: | ||
SAAB | SAAB | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
Baldwin-Lomax, k- | Baldwin-Lomax, k-ε / Wolfshtein | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Line 155: | Line 154: | ||
BAe | BAe | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
Baldwin-Lomax, k- | Baldwin-Lomax, k-ε / Chien, Bae k-t | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
RANSMB | RANSMB | ||
Line 182: | Line 181: | ||
Helsinki UT | Helsinki UT | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
k-w, Menter k-w, Shear-Stress-Transport, Cebeci-Smith, k- | k-w, Menter k-w, Shear-Stress-Transport, Cebeci-Smith, k-ε / Chien | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
FINFLO | FINFLO | ||
Line 191: | Line 190: | ||
IMFT | IMFT | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
k- | k-ε / Chien, Jones-L, Lam-Bremhorst, L-Sharma, modified, two-layer | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Tzabiras et al(1986), | Tzabiras et al(1986), | ||
Line 220: | Line 219: | ||
Parallax | Parallax | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
k- | k-ε / Chien | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Line 229: | Line 228: | ||
SAAB | SAAB | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
k- | k-ε / Shear-Stress-Transport | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Saab code NS2D | Saab code NS2D | ||
Line 238: | Line 237: | ||
TU Berlin | TU Berlin | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
Baldwin-Lomax, Cebeci-Smith, Johnson-King, k- | Baldwin-Lomax, Cebeci-Smith, Johnson-King, k-ε / LLYap, k-w Menter, k-w Wilcox, k-w (own) | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Line 265: | Line 264: | ||
UMIST | UMIST | ||
| style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | | style="width: 144.95pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="193" valign="top" | | ||
Lien-Leschziner, Non-Linear k- | Lien-Leschziner, Non-Linear k-ε Shih+Wolfshtein, Gibson-L/Kraft-L, Gibson-L/Shir | ||
| style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | | style="width: 99.15pt; border-top: none; border-left: none; border-bottom: solid windowtext 1.0pt; border-right: solid windowtext 1.0pt; padding: 0cm 5.4pt 0cm 5.4pt" width="132" valign="top" | | ||
Revision as of 18:43, 10 March 2009
Flow around airfoils (and blades) A-airfoil (Ma=0.15,
Re/m=2x106)
Underlying Flow Regime 2-05 © copyright ERCOFTAC 2004
Test Case
Brief description of the study test case
The present case is one of several measured and studied computationally. The chord Reynolds number is Rec=2.07x106 and the Mach number M=0.15, close enough to zero to be considered effectively incompressible but also high enough for compressible CFD methods to be used and the results compared with incompressible methods. The [C/geometry.dat geometry.dat] is purely two-dimensional. The case of primary interest, corresponding closely to the maximum lift for the foil, is at an angle of attack of α=13.3º. These conditions are roughly similar to the take-off and landing stall conditions for a small aeroplane like a Cessna Cardinale, though the Reynolds number for the real aircraft would be somewhat higher at about three million.
The case had its origin in a French co-operative research project Opération Décrochage driven by Aerospatiale and including ONERA, CERT and IMFT as partners. A complete experimental database was compiled (see below) covering transition, boundary layer development and the separated region. It was noted even at that time that algebraic and k-epsilon models failed to predict the separation [U2-05ref.htm#ref6 [6]].
Test Case Experiments
The Aerospatiale A airfoil was designed in 1986. The UFR corresponds precisely to one case measured in the ONERA F2 tunnel at M=0.15 and Rec=2.1x106 measured by [U2-05ref.htm#ref7 Gleyzes [7,8]]. The range of angles of attack included α=13.3º, at which the CL was measured to be 1.515. There were also measurements carried out in [U2-05ref.htm#ref9 the F1 tunnel [9,10]] where the incidence of 13.1º was measured at Rec=2.07x106, and measurements at Rec= 3.13x106 were gathered.
Aside from fixing (or correctly predicting) the position of transition (dealt with below) there appears to be little clear unambiguous evidence of superiority of any particular model. Second- moment closures have been recommended as superior, but since simpler models can give convincing prediction in some circumstances, perhaps with suitable tuning of parameters, this conclusion cannot be commended as best practice.
EUROVAL |
turbulence model |
code employed |
discretisation scheme(s) and solver(s) |
BAe |
Baldwin-Lomax, k-ε / Wolfshtein |
||
CERFACS |
k-ε / Chien, Algebraic Stress Model |
||
NLR |
Baldwin-Lomax, BL-Goldberg |
||
SAAB |
Baldwin-Lomax, k-ε / Wolfshtein |
||
SAAB/HUT |
Cebeci-Smith |
||
TUD |
Cebeci-Smith |
||
ECARP |
|||
BAe |
Baldwin-Lomax, k-ε / Chien, Bae k-t |
RANSMB |
a finite volume structured multiblock solver, an explicit four-stage Runge-Kutta scheme |
DLR |
Johnson-King, Algebraic Stress Model:Yakhot-Orszag |
DLR ''ViB Code'' |
a multi-grid Navier-Stokes solver, using a cell-centered finite-volume Runge-Kutta method |
Dornier |
Baldwin-Lomax -Granville |
Dornier IKARUS, and a boundary layer code (Horton, 1991) |
the Dornier N-S methods,and a boundary layer method |
Helsinki UT |
k-w, Menter k-w, Shear-Stress-Transport, Cebeci-Smith, k-ε / Chien |
FINFLO |
a 2D/3D finite volume N-S solver, an implicit pseudo-time integration |
IMFT |
k-ε / Chien, Jones-L, Lam-Bremhorst, L-Sharma, modified, two-layer |
Tzabiras et al(1986), Braza(1986)/Jin&Braza(1994) |
based on the principles of predictor-corrector pressure-schemes, either the SIMPLE algorithm or the predictor-corrector pressure schemes |
KTH |
Baldwin-Lomax, k-t (with Speziale & y+) |
NSMB (Navier-Stokes Multi Block) |
structured multi block grids, the finite volume method with a central difference, and an explicit four stage Runge-Kutta scheme |
ONERA |
k (half-eqn) and k-uv (2 half-eqn) using adaptive meshing |
VIS07 |
Viscous-Inviscid-Solvers (VIS) |
Parallax |
k-ε / Chien |
centered differencing for diffusive fluxes, Bounded Skew Hybrid Differencing Scheme for convective fluxes, PISO(Pressure Implicit Split Operator) algorithm | |
SAAB |
k-ε / Shear-Stress-Transport |
Saab code NS2D |
a cell-centered finite volume approximation in space, an explicit five-step Runge-Kutta scheme in time |
TU Berlin |
Baldwin-Lomax, Cebeci-Smith, Johnson-King, k-ε / LLYap, k-w Menter, k-w Wilcox, k-w (own) |
Interactive Boundary-Layer (IBL), consisting of a variable direct/inverse boundary-layer procedure | |
TU Chemnitz |
Baldwin-Lomax, Cebeci-Smith, k-e Wolfshtein (with and without Johnstone blending) |
MUFLO |
a cell-centred, finite-volume method, an explicit time-stepping Runge-Kutta algorithm |
TU Denmark |
Baldwin-Lomax, Baldwin-Barth |
a central difference expressions for the Laplacians and either simple upwinding or a so-called plus/minus upwinding for the convective term, the ADI procedure of Peaceman-Rachford. | |
UMIST |
Lien-Leschziner, Non-Linear k-ε Shih+Wolfshtein, Gibson-L/Kraft-L, Gibson-L/Shir |
a general non-orthogonal, block-structured finite-volume scheme, advection is approximated by the UMIST scheme. |
A great variety of wall functions and wall-region treatments has also been employed in the previous studies, with mixed results, on the basis of which we are not in a position to
recommend one specific approach.
Boundary Conditions
Based on practices used in the ECARP project, we recommend the following for a compressible RANS calculation:
• At the airfoil boundary use adiabatic walls with a no-slip condition. Set the turbulence energy and normal derivative of dissipation to zero (or equivalent).
• At the far-field boundary, use absorbing / non-reflecting boundary conditions (1D Riemann invariants).
• At the inflow boundary, use Riemann invariants.
For an incompressible RANS calculation:
• Use no-slip conditions at the wall, with zero pressure gradient normal to the wall.
• Use second-order extrapolations at outflow boundary and elsewhere as required.
© copyright ERCOFTAC 2004
Contributors: Peter Voke - University of Surrey