Lecture 4 湍流数值模拟方法.pptVIP

大涡模拟 Large Eddy Simulation (LES) 大尺度各向异性- RANS无法准确模拟 小尺度各向同性- RANS可以准确模拟 是否能利用RANS来专门模拟小尺度+网格较大，步长较长的’DNS’来模拟大尺度？大涡模拟(LES)类似于这种方法。 用大网格来直接求解NS方程，方程好比被一个滤波器滤了波 此为求解的结果， 而真解为f 截取自Lesieur, Metais Comet Large eddy simulations of turbulence Cambridge press 脉动动能在各波数上的分布 如何处理被大网格滤掉的小尺度结构 Smagorinsky延伸了Prandtl混合长度理论 LES eddy viscosity 正比 subgrid-scale 特称长度 Δx 和 subgrid-scale 特征速度： 利用Smagorinsky涡流粘性可估算subgrid尺度雷诺应力, 此模型可展现Kolmogorov的k?5/3消散能量规律 这里 LES特点 抓大不放小 非常有利，有力的工具 是最近,可预见未来流体力学研究和应用的热点 近来又出现了VLES, DES等在LES上发展而来的工具 LES会最终取代RANS吗？ 不会。Journal of Fluids Engineering -V127, 5, pp. 831-839 (Will RANS survive LES? I think yes, at least for the next few decades.) 计算能力在增强，LES,DNS应用越来越多 但是RANS也在进步，应用范围越来越广 综合LES/RANS的计算方法会越来越多 湍流的几种数值模拟方法 2007年12月13日 内容 简要介绍CFD求解过程 NS方程的直接数值模拟(DNS) 求解雷诺平均后的NS方程(RANS) Prandtl 混合长度模型， 0方程 k-ε 模型， 2方程 雷诺应力模型 ，5-7方程 大涡模拟(LES) CFD求解过程 定义求解区域 划分网格 方程离散化 设置边界条件，速度或者压强等 循环直至收敛 若模拟非稳态场，进入下一时间，继续 定义域,划分网格 Le et al. 1996 Kaiktsis and Monkewitz 2003 背向台阶 Backward facing step 纳维 斯托克式方程Navier Stokes Equations 方程： 把动量方程在每个网格上离散， 如果三维计算，每个网格形成关于u1,u2,u3 (或者u,v,w)的3个方程。 *本讲所列方程均为不可压缩，忽略重力 利用泊松方程求压强Possion Equation 对动量方程的每项求散度 divergence 再利用连续性方程，可得泊松方程 div CFD求解 计算过程 给每个网格设定初值，设定边界值 For i=1 to N %从时间to开始，每步增加Δt While u,v,w 还在变化 For j= 1 to M %第1个到第M个网格 Compute u Compute v Compute w Compute p End End End 对网格大小和时间步长的要求 准确捕捉到流场所有特性 网格尺寸要小到最小涡旋长度尺度的若干分之一 时间步长要小到最小涡旋时间尺度的若干分之一 能够达到这种要求的模拟被称作湍流的直接数值模拟(Direct Numerical Simulation, DNS) DNS 一定是三维模拟 非常‘昂贵’，小尺度～1/ReL3/4 “the number of grid points and the cost required increase roughly with Re3” (Rodi 2006) 主要用在科研，来分析湍流的物理特性一般Re=2000~6000 Le et al, 1996, 背向台阶的直接模拟 Re=5100 grid 512x192x64 克雷机，每步10秒 共需24天计算时间 DNS 最近一两年，欧美几个研究机构开始计算一些工程中的流场 例如卡尔斯鲁厄大学W. Rodi小组，1千万网格，Re～105 200-400个CPU Rodi 2006 燃气轮机压缩机叶片的大尺度结构，Re=51800 雷诺平均后的方程 Reynolds Averaged Navier Stokes (RANS) 很多实际问题中只关心流场的平均量 利用雷诺分解 代入动量方程，对每项求平均，得 需要建立雷诺应力 与‘已知量’的关系，才能得到解 (close momentum equation) 涡流粘度 Eddy viscosity or tu