Araya Guillermo, Lagares Christian
High Performance Computing and Visualization Laboratory, Department of Mechanical Engineering, University of Puerto Rico, Mayaguez 00681, Puerto Rico.
Entropy (Basel). 2022 Apr 15;24(4):555. doi: 10.3390/e24040555.
We employ numerically implicit subgrid-scale modeling provided by the well-known streamlined upwind/Petrov-Galerkin stabilization for the finite element discretization of advection-diffusion problems in a Large Eddy Simulation (LES) approach. Whereas its original purpose was to provide sufficient algorithmic dissipation for a stable and convergent numerical method, more recently, it has been utilized as a subgrid-scale (SGS) model to account for the effect of small scales, unresolvable by the discretization. The freestream Mach number is 2.5, and direct comparison with a DNS database from our research group, as well as with experiments from the literature of adiabatic supersonic spatially turbulent boundary layers, is performed. Turbulent inflow conditions are generated via our dynamic rescaling-recycling approach, recently extended to high-speed flows. Focus is given to the assessment of the resolved Reynolds stresses. In addition, flow visualization is performed to obtain a much better insight into the physics of the flow. A weak compressibility effect is observed on thermal turbulent structures based on two-point correlations (IC vs. supersonic). The Reynolds analogy (u' vs. t') approximately holds for the supersonic regime, but to a lesser extent than previously observed in incompressible (IC) turbulent boundary layers, where temperature was assumed as a passive scalar. A much longer power law behavior of the mean streamwise velocity is computed in the outer region when compared to the log law at Mach 2.5. Implicit LES has shown very good performance in Mach 2.5 adiabatic flat plates in terms of the mean flow (i.e., Cf and UVD+). iLES significantly overpredicts the peak values of u', and consequently Reynolds shear stress peaks, in the buffer layer. However, excellent agreement between the turbulence intensities and Reynolds shear stresses is accomplished in the outer region by the present iLES with respect to the external DNS database at similar Reynolds numbers.
我们采用了由著名的流线迎风/彼得罗夫-伽辽金稳定化方法提供的数值隐式亚网格尺度建模,用于大涡模拟(LES)方法中平流扩散问题的有限元离散化。尽管其最初目的是为稳定且收敛的数值方法提供足够的算法耗散,但最近它已被用作亚网格尺度(SGS)模型,以考虑离散化无法解析的小尺度效应。自由流马赫数为2.5,并与我们研究小组的直接数值模拟(DNS)数据库以及绝热超声速空间湍流边界层文献中的实验进行了直接比较。通过我们最近扩展到高速流动的动态重缩放-循环方法生成湍流流入条件。重点是对解析雷诺应力的评估。此外,进行流动可视化以更好地洞察流动物理。基于两点相关性(不可压缩流与超声速流),在热湍流结构上观察到了弱压缩性效应。雷诺类比(u' 与 t')在超声速区域大致成立,但程度小于在不可压缩(IC)湍流边界层中先前观察到的情况,在不可压缩湍流边界层中温度被视为被动标量。与马赫数为2.5时的对数律相比,在外区计算出的平均流向速度的幂律行为要长得多。隐式LES在马赫数为2.5的绝热平板的平均流(即摩擦系数Cf和UVD+)方面表现出非常好的性能。隐式LES在缓冲层中显著高估了u' 的峰值,进而高估了雷诺剪切应力峰值。然而,在相似雷诺数下,当前的隐式LES在外区的湍流强度和雷诺剪切应力方面与外部DNS数据库达成了极好的一致性。