Development of an Interpolated RANS-LES solver for wall-bounded turbulent fluid flows

Abstract

Scale-resolving hybrid RANS-LES models are being increasingly used to numerically simulate wall-bounded turbulent flows since they are computationally efficient compared to LES (Large Eddy Simulation) while also being more accurate compared to RANS (Reynolds-Averaged Navier–Stokes) models. Existing hybrid RANS-LES models typically either suffer from modeled stress depletion or they can introduce a sharp jump in the mean Reynolds stress across the RANS-LES zone. We present a novel interpolated RANS-LES (IRL) solver implemented in OpenFOAM, in which the RANS equations (for eddy-viscosity) and LES equations (for filtered velocity) are evolved simultaneously on the same computational grid. Based on the distance from the no-slip wall, the flow domain is demarcated into a near-wall, free-stream, and a “hybrid” region, sandwiched between the near-wall and free-stream regions. Spalart Allmaras (S-A) model has been used to evolve the RANS eddy viscosity, while the Wall Adapting Local Eddy Viscosity (WALE) model is used for evolving the filtered LES velocity. Using a novel scheme based on solving a partial differential equation, the turbulent eddy-viscosity is interpolated between the RANS eddy viscosity in the near wall region and the effective eddy viscosity from the resolved LES fluctuations in the free-stream region. The mean subgrid stress in the LES momentum equation is then corrected in the near-wall and hybrid regions using the interpolated turbulent eddy-viscosity. IRL’s grid resolution requirements are the same as DES’s (Detached Eddy Simulation). Grid convergence and sensitivity studies have been carried out for two canonical flow geometries. For turbulent channel flow simulations, the IRL solver does not display the typical mean stress depletion observed in DES while also predicting resolved Reynolds stresses accurately in the free-stream region. For simulations of flow over backward facing step, the IRL solver can predict the friction factor and coefficient of pressure better than that DES, especially after the reattachment point.