UFR 2-12 Evaluation

From KBwiki
Jump to navigation Jump to search

Turbulent Flow Past Two-Body Configurations

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Flows Around Bodies

Underlying Flow Regime 2-12

Evaluation

Comparison of CFD Calculations with Experiments

This section is organized as follows. First (Section 6.1), results of some sensitivity studies are presented and briefly discussed. These include evaluation of such effects as span-size of the domain, compressibility, time sample used for computing the mean flow and turbulent statistics, and numerical dissipation of the method used. Then, in Section 6.2, a comparison with the experimental data is shown for the main body of simulations carried out within the ATAAC project with the use of the physical and computational problem setups outlined in Section 5.

RESULTS OF SENSITIVITY STUDIES

Effect of span size of domain

As mentioned in Section 4, the aspect ratio of the CT configuration Lz/ D in the BART facility is equal to 12.4. Strictly speaking this demands carrying out simulations exactly at this value of Lz/ D and imposing no-slip boundary conditions on the floor and ceiling of the test section (see Figure 2). However such simulations would be very expensive. Considering this and, also, recommendations of the BANC-I Workshop based on simulations at different Lz/ D with periodic boundary conditions in the spanwise directions, most of the simulations in the ATAAC project were performed at Lz/ D = 3 assuming spanwise periodicity. In order to get an idea on how strong the effect of such a simplification could be, NTS conducted a series of simulations at different Lz/ D. Some results of these simulations are presented below [1]


Figure 4 compares flow visualisations from the SA DDES carried out in the "mandatory" (Lz/ D = 3) and the widest of the considered domains (Lz/ D = 16) in the form of instantaneous isosurface of the magnitude of the second eigenvalue of the velocity gradient tensor or "swirl" quantity, λ2. The figure is reassuring in the sense that it visibly displays that the narrow-domain simulation resolves not only fine-grained turbulent eddies but also large, nearly coherent, structures and exhibits all the complex flow features observed in the visualization of the wide-domain simulation, except for the initial region of the free shear-layer separated from the upstream cylinder, where a noticeable difference between the two flow-visualizations is observed.


UFR2-12 figure4a.png UFR2-12 figure4b.png
Figure 4: Isosurface of λ2 = 4.0(U0 /D ) from incompressible SA DDES at Lz = 3D and 16D.


As a result, sensitivity of predictions of the major characteristics of the flow in the wake of the downstream cylinder to the value of Lz/D turns out to be marginal (see Figure 5 and Table 4). At the same time, as seen in Figure 6, the flow features directly related to the details of the flow past the upstream cylinder (its boundary layers separation and shear-layers roll-up) vary with Lz/D variation rather significantly. Other than that, Figures 5, 6 suggest that the effect of Lz/D within different turbulence modelling approaches is different and is stronger pronounced for IDDES than for DDES. These findings should be kept in mind when analyzing agreement with the experiment of the simulations carried out at Lz/D = 3 with the use of different approaches to turbulence representation presented in the next section.


UFR2-12 figure5a.png UFR2-12 figure5b.png UFR2-12 figure5c.png
Figure5: Effect of Lz/D on centreline distributions of mean velocity and 2D kinetic energy in the wake of the downstream cylinder (left and middle frames respectively) and on rms of pressure coefficient on the downstream cylinder (right) from incompressible SA DDES and SA IDDES of NTS.


Table 4: Effect of modelling approach and span-size of computational domain on the global shedding frequency.
Shedding frequency, Hz
DDES, Lz = 3D 188
DDES, Lz = 16D 188
IIDDES, Lz = 3D 192
IIDDES, Lz = 16D 192


UFR2-12 figure6a.png UFR2-12 figure6b.png UFR2-12 figure6c.png
Figure6: Effect of Lz/D on centreline distributions of mean velocity and 2D kinetic energy in the gap between cylinders (left and middle frames respectively) and on rms of pressure coefficient on the upstream cylinder (right) from incompressible SA DDES and IDDES of NTS.



  1. The simulations at Lz/ D = 16 were conducted with the use of resources of the Argonne Leadership Computing Facility at Argonne National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under contract DE-AC02-06CH11357.

Compressibility effects

Considering that TUB simulations were carried out under assumption of incompressible flow and BTU and DLR performed compressible simulation at the experimental value of the Mach number (M = 0.128), it was important to find out how strongly this could affect the obtained results. In order to gain this knowledge NTS has carried out IDDES simulations both in the framework of incompressible and compressible problem statements. Results of these two simulations for the quantities which were found to be most sensitive to Lz/D (see Figure 6) are shown in Figure 7. They suggest that the role of the effects of compressibility in the considered flow is negligible.


UFR2-12 figure7a.png UFR2-12 figure7b.png UFR2-12 figure7c.png
Figure7: Comparison of results of SA IDDES at Lz/D = 3 carried out in incompressible and compressible (M = 0.128) modes. Left and middle frames: centreline distributions of mean velocity and 2D kinetic energy in the gap between cylinders; right frame: rms of pressure coefficient on the upstream cylinder.


Effect of time sample

This effect has been shown to be very strong for nominally 2D bodies with massive separation in many previous DES and LES studies. This is not surprising, since such flows typically have large scale coherent vortices in a vortex street pattern overlaid with finer random turbulent fluctuations at higher frequencies and random modulation and intermittency at frequencies lower than the vortex shedding one. This is true for the TC flow as well (see Figure 4), which dictates a need of rather long time samples for getting statistically representative mean flow characteristics and especially turbulence statistics. In order to exclude or at least minimize the effect of insufficient time sample when comparing results of different simulations, NTS has carried out a time-sample sensitivity study of the SA DDES of the TC flow. Its outcome is presented in Figure 8 and in Table 5. One can see that the time sample of about 150 convective time units (D / U0) is sufficient to obtain a reliable statistics not only for the mean drag but also for rms of the pressure coefficient. Exactly this value was recommended as a minimum one within the ATAAC project.


UFR2-12 figure8a.png UFR2-12 figure8b.png UFR2-12 figure8c.png
Figure8: Effect of time sample on running time average of the drag coefficient (left frame) and rms of pressure coefficient on the surface of upstream (middle frame) and downstream (right frame) cylinders of tandem at Lz/D = 3 from incompressible SA IDDES


Table 5: Effect of time sample on time average of integral of TC configuration
Time, Convective time units (D/U0) CD upstream cylinder CD downstream cylinder
100 0.501 0.410
200 0.510 0.403
300 0.505 0.404


Effect of numerical dissipation

In order to evaluate this effect, NTS has carried compressible SA IDDES of the flow at Lz = 3D on the mandatory grid with the use of different approximations of the inviscid fluxes, Finv, available in the NTS code. All these approximations are based on the weighted upwind-biased (Fupw ) and central (Fctr ) schemes

Finv = (1 - σupw</sub> )  Fctr + σupw Fupw

with the solution-dependent empiric weight-function of the upwind scheme, σupw, computed as

σupw = max {σmax tanh(ACH1), σmin }

.

The function is designed so that σupw ∈ [σminσmax ] and is close to σmin in the LES region of DES and close to σmax in the RANS, irrotational, and departure regions [14]. The series of simulations performed included three simulations, all carried out for compressible flow in free air rather than in the wind tunnel. The first and the second ones had σmin = 0 and combined the 4th order central scheme with 5th and 3rd order upwind-biased schemes, respectively, whereas the third one had σmin = 0.2 and combined 4th order central scheme with 3rd order upwind-biased scheme. Thus the first scheme has the minimum and the third scheme the maximum numerical dissipation with the second scheme somewhere in between [1]. Some results from this series are presented in Figures 9, 10.


UFR2-12 figure9a.png UFR2-12 figure9b.png UFR2-12 figure9c.png
UFR2-12 figure9d.png UFR2-12 figure9e.png
Figure9: Effect of increase of numerical dissipation on snapshots of spanwise vortricity (upper row) and Power Spectral Density (PSD) of SA-IDDES streamwise velocity at a hot point in the wake of the downstream cylinder (lower row).


UFR2-12 figure10a.png UFR2-12 figure10b.png UFR2-12 figure10c.png
Figure10: Effect of increase of numerical dissipation on centreline distributions of SA-IDDES mean velocity and 2D TKE in the gap between cylinders (left and middle frames) and on rms of pressure coefficient on the upstream cylinder (right frame).


They clearly show that increase of numerical dissipation causes a visible delay of the roll-up of the shear layer separated from the upstream cylinder and general damping of resolved turbulent structures (Figure 9). This, in turn, leads to a tangible variation not only of the unsteady pressure but also the mean flow parameters (Figure 10), which turns out to be comparable with their variation caused by the variation of the span size of the computational domain (see Figure 6 above). Based on these results, it can be concluded that using as low dissipative scheme as possible is essential for simulation of the considered and similar flows. Also, the fact that even relatively mild increase of the upwinding results in a considerable variation of the mean flow predictions suggests that the mandatory grid is not "fine-enough", which should be kept in mind when assessing agreement of results of simulations with the use of different turbulence modelling approaches with experiment presented in Section 2.1.2 below.




  1. Note that all the NTS results presented in the section 2.1.2 below are obtained with the use of the first, less dissipative, scheme.

COMPARISON WITH EXPERIMENT

In this section we present a comparison of predictions based on the simulations summarized in Table 3 above with experimental data [1]


UFR2-12 figure11a.png Experiment UFR2-12 figure11b.png
UFR2-12 figure11c.png UFR2-12 figure11d.png
UFR2-12 figure11e.png UFR2-12 figure11f.png
UFR2-12 figure11g.png UFR2-12 figure11h.png
UFR2-12 figure11i.png
Figure 11: Comparison of predicted and experimental instantaneous vorticity fields


Figure 11 presents flow visualizations from the PIV and from all the simulations in the form of instantaneous contours of the spanwise vorticity component. One can see that all the simulations predict qualitatively similar and generally plausible (consistent with the experiment) turbulent structures both in the gap between the cylinders and in the wake of the downstream cylinder. This suggests that the modelling approaches used, in principle, are capable of representing the key physics of the flow and that the numerical methods applied in all the codes possess essential features needed for turbulence- resolving simulations. At the same time, the comparison reveals considerable difference between the simulations which is most pronounced in the initial region of the free shear layers separated from the upstream cylinder, the peculiarities being associated with both turbulence models and codes.


We now move to a quantitative comparison with the experiment of the computed mean flow characteristics, particularly, distributions of mean pressure coefficient Cp over the surface of both cylinders and centreline distributions of the mean streamwise velocity in the gap between the cylinders and in the wake of the downstream cylinder.


UFR2-12 figure12a.png UFR2-12 figure12b.png
UFR2-12 figure12c.png UFR2-12 figure12d.png
Figure 12: Mean pressure coefficient distributions on the upstream (left column) and downstream (right column) cylinders.




  1. For NTS, results of the incompressible simulations are shown obtained at Lz = 3 with the use of hybrid, 5/4, scheme at σmin = 0.


Contributed by: A. Garbaruk, M. Shur and M. Strelets — New Technologies and Services LLC (NTS) and St.-Petersburg State Polytechnic University

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2024