Numerical Modeling of Freestream Turbulence Decay Using Different Commercial Computational Fluid Dynamics Codes

This work models the spatial decay of freestream turbulence using three different commercial computational fluid dynamics (CFD) codes: Fluent, star-ccm+, and cfx. The two-equation shear stress transport k–ω (SST-k–ω) steady Reynolds-averaged-Navier–Stokes (RANS) model was used, within each of these three different commercial codes, and the modeling variations were analyzed. Comparison of the results from the SST-k–ω model with experiments and large eddy simulation (LES) (carried out using star-ccm+) were also made, which reveal that all the commercial CFD codes demonstrate either a higher or slower rate of spatial turbulent kinetic energy (TKE) decay. Attempts were then made to unify the resultant modeling approach between these three CFD tools, by careful manipulation of the inlet boundary conditions and subsequent fine-tuning of the SST-k–ω model constant (β∞∗). The results obtained not only displayed uniformity among the three CFD codes but also demonstrated a much better agreement to the experiments and the LES results. Thereafter, the optimized model coefficient (β∞∗) was integrated with the three-equation k–kl–ω transition model to examine its applicability in modeling a turbulent boundary layer flow over a flat plate with low incoming turbulence. The results showed good agreement with the theoretical boundary layer correlations, with correct prediction of the transition location. The findings from this study can be used as a suitable modeling method to accurately model the effects of freestream turbulence on bluff-body and boundary layer flows.