直接数值模拟和大涡模拟后台阶湍流流动.pdfVIP

直接数值模拟与大涡模拟后台阶湍流流动 摘 要 本文编制了一套二维湍流的计算机程序，可用于在交错网格系统下对后台 阶湍流流动的流场进行直接数值模拟(DNS)与大涡模拟(LES) 。在DNS 中，采用 改进的MAC 方法求解非定常N-S 方程，时间项离散采用二阶格式；在LES 中， 采用 k 方程亚网格尺度模型(KSGS)描述亚网格湍流粘性，SIMPLE 算法和混合 差分格式对控制方程进行求解。并将两种数值模拟方法(改进的MAC 与KSGS) 的计算结果与实验数据进行对比，研究结果表明，两种方法都能得到令人满意 的计算结果，但DNS 结果要比LES 更加符合实验数据，其精度也高于LES ，同 时直接数值模拟的计算工作量也要比大涡模拟大得多。 本文在DNS 中还分别采用二阶显式混合差分格式和四阶紧致格式离散动量 方程；采用隐式五点中心差分和九点中心差分格式离散压力Poisson 方程，并数 值分析不同差分格式对DNS 计算结果和误差的影响，研究表明不同差分格式对 后台阶湍流流动和收敛性都有较大的影响。数值研究不同的进口速度分布对 DNS 所得的后台阶湍流流场的影响，与实验数据比较表明，采用正切分布所得 的DNS 结果与实验值符合要比正弦分布更好些。 关键词：直接数值模拟，大涡模拟，后台阶湍流流动 iv 南京航空航天大学硕士论文 Abstract A general computer code has been developed for direct numerical simulation(DNS) and large eddy simulation(LES) of turbulent flow over a backward-facing step in a two-dimensional duct. For DNS，the modified Marker-and Cell(MAC) method is adopted for solving time-dependent viscous incompressible Navier-Stokes equations on a staggered grid system and the time derivative is discretized using the second-order accurate implicit formulation. For LES，the k-equation sub-grid scale model was used to simulate the sub-grid eddy viscosity and SIMPLE algorithm and the hybrid difference scheme are applied to large eddy simulation for turbulent flows. Predictions of recirculation region length and axial velocity profiles are in fairly agreement with experimental data. It reveals that both two numerical methods can be used to predict transient turbulent flows and the accuracy of DNS is higher than LES， but computation time of DNS is more than LES. For DNS，the explicit finite difference methods based on second-order accurate hybrid and fourth-order accurate compact difference formulae are used to solve momentum equations. A Poisson equation for pressure corre