An improved anisotropy-resolving subgrid-scale model for flows in laminar-turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, kSGS, containing a production term proportional to the square root of kSGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes kSGS to self-reproduce. Consequently, the laminar-turbulent transition region predicted by this model depends on the inflow or initial condition of kSGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the kSGS transport equation. In the MTS model, the SGS turbulence energy, kes, estimated by filtering the instantaneous flow field is used. Since the kes approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional kSGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.