POD-based wall boundary conditions for the numerical simulation of turbulent channel flows

Abstract

Simulation of turbulent wall-bounded flows requires a high spatial resolution in the wall region, which limits the range of Reynolds numbers which can be effectively reached. In previous work, we proposed proper orthogonal decomposition (POD) based wall boundary conditions to bypass the simulation of the inner wall region. Tests were carried out for direct numerical simulation at a low Reynolds number Reτ = 180. The boundary condition is based on the POD spatial eigenfunctions which are determined a priori in the full channel. It consists of a three-component velocity field on the plane y+ = 50 which is reconstructed at each instant from a combination of selected eigenfunctions. The coefficients of the combination are estimated from the simulation in the reduced domain using the threshold-based reconstruction method described in Podvin et al. The study is now extended to large-eddy simulation at higher Reynolds numbers Reτ = 295 and Reτ = 590. Two versions of the reconstruction method are considered. In the first version, both the phases and the moduli of the coefficients are allowed to vary. In the second version, only the phases are adjusted. We find that the latter method is associated with improved statistics and is relatively robust with respect to the reconstruction threshold. However, it is sensitive to the details of the numerical simulation, unlike the former method, which is associated with less accurate statistics and is more dependent on the reconstruction threshold.