Simulations of turbulent rotating flows using a subfilter scale stress model derived from the partially integrated transport modeling method

Abstract

The partially integrated transport modeling (PITM) method [B. Chaouat and R. Schiestel, “A new partially integrated transport model for subgrid-scale stresses and dissipation rate for turbulent developing flows,” Phys. Fluids 17, 065106 (2005)10.1063/1.1928607; R. Schiestel and A. Dejoan, “Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations,” Theor. Comput. Fluid Dyn. 18, 443 (2005)10.1007/s00162-004-0155-z; B. Chaouat and R. Schiestel, “From single-scale turbulence models to multiple-scale and subgridscale models by Fourier transform,” Theor. Comput. Fluid Dyn. 21, 201 (2007)10.1007/s00162-007-0044-3; B. Chaouat and R. Schiestel, “Progress in subgrid-scale transport modelling for continuous hybrid non-zonal RANS/LES simulations,” Int. J. Heat Fluid Flow 30, 602 (2009)10.1016/j.ijheatfluidflow.2009.02.021] viewed as a continuous approach for hybrid RANS/LES (Reynolds averaged Navier-Stoke equations/large eddy simulations) simulations with seamless coupling between RANS and LES regions is used to derive a subfilter scale stress model in the framework of second-moment closure applicable in a rotating frame of reference. This present subfilter scale model is based on the transport equations for the subfilter stresses and the dissipation rate and appears well appropriate for simulating unsteady flows on relatively coarse grids or flows with strong departure from spectral equilibrium because the cutoff wave number can be located almost anywhere inside the spectrum energy. According to the spectral theory developed in the wave number space [B. Chaouat and R. Schiestel, “From single-scale turbulence models to multiple-scale and subgrid-scale models by Fourier transform,” Theor. Comput. Fluid Dyn. 21, 201 (2007)10.1007/s00162-007-0044-3], the coefficients used in this model are no longer constants but they are some analytical functions of a dimensionless parameter controlling the spectral distribution of turbulence. The pressure-strain correlation term encompassed in this model is inspired from the nonlinear SSG model [C. G. Speziale, S. Sarkar, and T. B. Gatski, “Modelling the pressure-strain correlation of turbulence: an invariant dynamical systems approach,” J. Fluid Mech. 227, 245 (1991)10.1017/S0022112091000101] developed initially for homogeneous rotating flows in RANS methodology. It is modeled in system rotation using the principle of objectivity. Its modeling is especially extended in a low Reynolds number version for handling non-homogeneous wall flows. The present subfilter scale stress model is then used for simulating large scales of rotating turbulent flows on coarse and medium grids at moderate, medium, and high rotation rates. It is also applied to perform a simulation on a refined grid at the highest rotation rate. As a result, it is found that the PITM simulations reproduce fairly well the mean features of rotating channel flows allowing a drastic reduction of the computational cost in comparison with the one required for performing highly resolved LES. Overall, the mean velocities and turbulent stresses are found to be in good agreement with the data of highly resolved LES [E. Lamballais, O. Metais, and M. Lesieur, “Spectral-dynamic model for large-eddy simulations of turbulent rotating flow,” Theor. Comput. Fluid Dyn. 12, 149 (1998)]. The anisotropy character of the flow resulting from the rotation effects is also well reproduced in accordance with the reference data. Moreover, the PITM2 simulations performed on the medium grid predict qualitatively well the three-dimensional flow structures as well as the longitudinal roll cells which appear in the anticyclonic wall-region of the rotating flows. As expected, the PITM3 simulation performed on the refined grid reverts to highly resolved LES. The present model based on a rational formulation appears to be an interesting candidate for tackling a large variety of engineering flows subjected to rotation.