In this paper, we introduce the discrete Maxwellian equilibrium distribution function for incompressible flow and force term into the two-stage third-order Discrete Unified Gas-Kinetic Scheme (DUGKS) for simulating low-speed turbulent flows. The Wall-Adapting Local Eddy-viscosity (WALE) and Vreman sub-grid models for Large-Eddy Simulations (LES) of turbulent flows are coupled within the present framework. Meanwhile, the implicit LES are also presented to verify the effect of LES models. A parallel implementation strategy for the present framework is developed, and three canonical wall-bounded turbulent flow cases are investigated, including the fully developed turbulent channel flow at a friction Reynolds number (Re) about 180, the turbulent plane Couette flow at a friction Re number about 93 and lid-driven cubical cavity flow at a Re number of 12000. The turbulence statistics, including mean velocity, the r.m.s. fluctuations velocity, Reynolds stress, etc. are computed by the present approach. Their predictions match precisely with each other, and they are both in reasonable agreement with the benchmark data of DNS. Especially, the predicted flow physics of three-dimensional lid-driven cavity flow are consistent with the description from abundant literature. The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.

中文翻译：

---

高阶离散统一气体动力学方案的壁面湍流大涡模拟

在本文中，我们将不可压缩流和力项的离散麦克斯韦平衡分布函数引入到用于模拟低速湍流的两级三阶离散统一气体动力学方案（DUGKS）中。用于湍流的大涡模拟（LES）的适应壁的局部涡粘性（WALE）和弗里曼子网格模型在本框架内耦合。同时，还提出了隐式LES以验证LES模型的效果。开发了本框架的并行实施策略，并研究了三种规范的有边界边界的湍流情况，包括在约180的摩擦雷诺数（Re）下完全开发的湍流通道流，Re值为93时的湍流平面Couette流和Re值为12000时由盖驱动的立方腔流。采用本方法可计算湍流统计数据，包括平均速度，均方根波动速度，雷诺应力等。 。它们的预测彼此精确匹配，并且都与DNS的基准数据合理地吻合。特别是，三维盖驱动腔流的预测流物理学与大量文献中的描述是一致的。目前的数值结果证明，目前基于二级三阶基于DUGKS的LES方法能够模拟不均匀的壁面湍流，并能在相对粗糙的网格上获得可靠的结果。雷诺应力等通过本方法计算。它们的预测彼此精确匹配，并且都与DNS的基准数据合理地吻合。特别是，三维盖驱动腔流的预测流物理学与大量文献中的描述是一致的。目前的数值结果证明，目前基于二级三阶基于DUGKS的LES方法能够模拟不均匀的壁面湍流，并能在相对粗糙的网格上获得可靠的结果。雷诺应力等通过本方法计算。它们的预测彼此精确匹配，并且都与DNS的基准数据合理地吻合。特别是，三维盖驱动腔流的预测流物理学与大量文献中的描述是一致的。目前的数值结果证明，目前基于二级三阶基于DUGKS的LES方法能够模拟不均匀的壁面湍流，并能在相对粗糙的网格上获得可靠的结果。三维盖驱动腔流的预测流物理学与大量文献中的描述是一致的。目前的数值结果证明，目前基于二级三阶基于DUGKS的LES方法能够模拟不均匀的壁面湍流，并能在相对粗糙的网格上获得可靠的结果。三维盖驱动腔流的预测流物理学与大量文献中的描述是一致的。目前的数值结果证明，目前基于二级三阶基于DUGKS的LES方法能够模拟不均匀的壁面湍流，并能在相对粗糙的网格上获得可靠的结果。