A Hybrid Large Eddy Simulation Algorithm Based on the Implicit Domain Decomposition

The resolution of small near-wall eddies encountered in high-Reynolds number flows using large eddy simulation (LES) requires very fine meshes that may be computationally prohibitive. As a result, the use of wall-modeled LES as an alternative is becoming more popular. In this paper, the near-wall domain decomposition (NDD) approach that was originally developed for Reynolds-averaged Navier–Stokes simulations (RANSs) is extended to the hybrid RANS/LES zonal decomposition. The algorithm is implemented in two stages. First, the solution is computed everywhere with LES on a coarse grid using a new non-local slip boundary condition for the instantaneous velocity at the wall. The solution is then recomputed in the near-wall region with RANS. The slip boundary conditions used in the first stage guarantee that the composite solution is smooth at the inner/outer region interface. Another advantage of the model is that the turbulent viscosity in the inner region is computed based on the corresponding RANS velocity. This shows improvement over those hybrid models that have only one velocity field in the whole domain obtained from LES. The model is realized in the open source code OpenFOAM with different approximations of turbulent viscosity and is applied to the planar channel flow at frictional Reynolds numbers of Re τ = 950 , 2000, and 4200. Mean streamwise velocity and Reynolds stress intensities are predicted reasonably well in comparison to the solutions obtained with unresolved LES and available DNS benchmarks. No additional forcing at the interface is required, while the log–layer mismatch is essentially reduced in all cases.