404
Views
0
CrossRef citations to date
0
Altmetric
Research Article

Errors and uncertainties in CFD validation for non-equilibrium turbulent boundary layer flows at high Reynolds numbers

, , , , , , , , & show all
Article: 2360195 | Received 16 Jan 2024, Accepted 21 May 2024, Published online: 30 May 2024

Figures & data

Figure 1. Left: Virginia Tech 2D smooth wall flow (2D VT) (from [Citation11]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.). Right: DLR/UniBw 2D smooth wall flow (reprinted with permission from [Citation16]).

Figure 1. Left: Virginia Tech 2D smooth wall flow (2D VT) (from [Citation11]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.). Right: DLR/UniBw 2D smooth wall flow (reprinted with permission from [Citation16]).

Figure 2. Left: 3D BeVERLI Hill flow (from [Citation17]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.). Right: Body-of-Revolution in turbulent flow.

Figure 2. Left: 3D BeVERLI Hill flow (from [Citation17]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.). Right: Body-of-Revolution in turbulent flow.

Figure 3. Left: Grid metrics using arithmetic average and root-mean-square cell areas for the grid levels for the 2D VT flow. Right: Grid levels 1, 3, 5 showing hierarchical refinement (zoom of the mesh for the 2D VT flow)(Figures from [Citation11]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.).

Figure 3. Left: Grid metrics using arithmetic average and root-mean-square cell areas for the grid levels for the 2D VT flow. Right: Grid levels 1, 3, 5 showing hierarchical refinement (zoom of the mesh for the 2D VT flow)(Figures from [Citation11]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.).

Figure 4. Convergence study for 2D VT flow. Left: Mesh convergence for Cf using the SST model using the CFD solver Kestrel. Middle and right: Code to code comparison for Cf and mean velocity profiles u+ vs. y+ on grid level 3(Figures from [Citation11]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.).

Figure 4. Convergence study for 2D VT flow. Left: Mesh convergence for Cf using the SST model using the CFD solver Kestrel. Middle and right: Code to code comparison for Cf and mean velocity profiles u+ vs. y+ on grid level 3(Figures from [Citation11]; reprinted by permission of the American Institute of Aeronautics and Astronautics, Inc.).

Figure 5. Illustration of the impact of an adjusted computational set-up compared to the original set-up for the 2D DLR/UniBw flow for the displacement thickness δ along the contour (left) and for the mean velocity profile in the region of strong adverse pressure gradient at x=9.944m (right).

Figure 5. Illustration of the impact of an adjusted computational set-up compared to the original set-up for the 2D DLR/UniBw flow for the displacement thickness δ∗ along the contour (left) and for the mean velocity profile in the region of strong adverse pressure gradient at x=9.944m (right).

Figure 6. Cp and Cf along the centerline in the streamwise x1-direction (left) and the centerspan in the spanwise x3-direction (right) from the 45 yaw case BeVERLI Hill at ReH650,000 on grid levels GL1 to GL4 (from [Citation13]; reprinted by permission of the American Society of Mechanical Engineers ASME). Note that x1x and x3z in Figure  (left). The dotted line shows the hill geometry.

Figure 6. Cp and Cf along the centerline in the streamwise x1-direction (left) and the centerspan in the spanwise x3-direction (right) from the 45∘ yaw case BeVERLI Hill at ReH≈650,000 on grid levels GL1 to GL4 (from [Citation13]; reprinted by permission of the American Society of Mechanical Engineers ASME). Note that x1≡x and x3≡z in Figure 2 (left). The dotted line shows the hill geometry.

Figure 7. Uncertainty due to the discretisation error (DE) in Cf along the centerline in x1-direction (top) and the centerspan in x3-direction (bottom) of the 45 yaw case BeVERLI Hill at ReH650,000 for SA-model (left) and SST-model (right) (from [Citation13]; reprinted by permission of the American Society of Mechanical Engineers ASME). Note that x1x and x3z in Figure  (left). The dotted line shows the hill geometry.

Figure 7. Uncertainty due to the discretisation error (DE) in Cf along the centerline in x1-direction (top) and the centerspan in x3-direction (bottom) of the 45∘ yaw case BeVERLI Hill at ReH≈650,000 for SA-model (left) and SST-model (right) (from [Citation13]; reprinted by permission of the American Society of Mechanical Engineers ASME). Note that x1≡x and x3≡z in Figure 2 (left). The dotted line shows the hill geometry.

Figure 8. Left: Uncertainty bars of the discretisation error for Cf for the SST model for two different families of meshes (baseline meshes by VT and meshes with refinement around the airfoil by MARIN/IST (MI)). Right: Estimate of the error of the turbulence model δtm for Cf with reference model SST 2003 for model-scale Re=2×106 and for full-scale Re=109(Figures reprinted from [Citation34]).

Figure 8. Left: Uncertainty bars of the discretisation error for Cf for the SST model for two different families of meshes (baseline meshes by VT and meshes with refinement around the airfoil by MARIN/IST (MI)). Right: Estimate of the error of the turbulence model δtm for Cf with reference model SST 2003 for model-scale Re=2×106 and for full-scale Re=109(Figures reprinted from [Citation34]).

Figure 9. Mean velocity profiles in viscous units for the 2D VT flow at model-scale Re=2×106 and at full scale Re=109 at x=3.5m for the airfoil incidence angles α=10 (left), α=0 (middle), and α=12 (right)(Figures reprinted from [Citation34]).

Figure 9. Mean velocity profiles in viscous units for the 2D VT flow at model-scale Re=2×106 and at full scale Re=109 at x=3.5m for the airfoil incidence angles α=−10∘ (left), α=0∘ (middle), and α=12∘ (right)(Figures reprinted from [Citation34]).