二维圆柱绕流的数值模拟.docxVIP

二维圆柱绕流的数值模拟 0 大试现象的应用 一个世纪以来，圆曲流问题一直是经典水流问题之一，也是许多理论分析、数值模拟和实验研究的对象。尽管研究历经时间如此之长,但这种最简单的流动现象物理本质的理解依然不完整。闵强利计算了2种雷诺数Re=200,800情况下的圆柱绕流,发现在Re=200时产生的卡门涡街非常规则,而随着雷诺数的增加,Re=800时,涡街出现不规则现象,并表现出三维效应。文献研究了不同截面形状在低雷诺数下的绕流场,并重点研究了阻力、升力的周期性变化特性。 随着计算机硬件的不断发展,大涡模拟作为一种新的湍流模型相比原来的湍流模型在求解非稳态流动方面具有巨大的潜力。大涡模拟方法的基本思想是用瞬时的N-S方程直接模拟计算湍流中的大尺度涡,而小尺度涡对大涡的影响则通过建立近似的模型来考虑,这种影响模型称为亚格子尺度模型。王汉青等人介绍了大涡模拟的理论进展和发展趋势, 描述了当前大涡模拟在工程中的具体应用。指出大涡模拟在模拟计算从层流到湍流转换、非定常湍流和高速湍流方面具有其他湍流模型无可比拟的优势。傅慧萍以潜艇模型SUBOFF为研究对象,采用大涡模拟(LES)方法求解了流动的非定常解。通过与试验值以及采用RNGk-ε湍流模型得到的定常结果比较,验证了大涡模拟方法的有效性。 本文采用大涡模拟方法求解低雷诺数下的圆柱绕流场,研究卡门涡街的周期性特性。对一个完整的涡街脱落周期进行详细讨论,揭示了不同时刻旋涡的生成、发展、脱落和演化过程。通过研究圆柱表面的压力变化说明产生周期性脱落涡对的原因,以及产生周期性变化的升、阻力的原因。 1 亚品格尺度模型 自然界中的流动现象基本可以通过基于连续性假设的纳维叶-斯托克斯方程(N-S方程)来描述。在应用于大涡模拟方法时,需要对N-S方程进行一定过滤处理。 ?ˉui?t+?ˉuiˉuj?xj=-1ρ?ˉp?xi+v?2ˉui?xj?xj-?τij?xj?uˉi?t+?uˉiuˉj?xj=?1ρ?pˉ?xi+v?2uˉi?xj?xj??τij?xj, (1) ?ˉui?xi=0?i?{1,2,3}?uˉi?xi=0?i?{1,2,3}。 (2) 对于亚格子尺度模型,采用标准Smagorinsky模型的方法。 τij-13τkkδij=-2vtˉSij=-2(CSˉΔ)2∥S∥ˉSijτij?13τkkδij=?2vtSˉˉij=?2(CSΔˉˉˉ)2∥S∥Sˉˉij。 (3) 式中:S为应变速率张量,S=12(?uj?xi+?ui?xj)S=12(?uj?xi+?ui?xj),‖S‖为进行滤过后的应变速率张量,∥S∥=√2ˉSijˉSij∥S∥=2SˉˉijSˉˉij???????√;vt为涡粘度,vt=l2‖S‖。本文计算所采用的网格均为四边形结构网格,滤过范围可取为 ˉΔ=2vol13Δˉˉˉ=2vol13。 (4) 式中,vol是计算单元的体积。 亚格子尺度可表示为 l=CsˉΔ(1-exp(-y+?25)3)1/2l=CsΔˉˉˉ(1?exp(?y+/25)3)1/2。 (5) 式中,Smagorinsky常量Cs取为0.255。 2 计算雷诺数和流动 本文所采用的计算模型为直径D=0.04 m的圆柱。来流速度为U=0.01 m/s,计算的雷诺数Re=400。计算域为方形,沿着来流方向长度为7.5D,尾流去流段长度为40D,两侧宽度为7.5D。流动在选定的计算域内能够充分发展,圆柱周围采用边界层网格,能够捕捉圆柱周围的细微流动。图1为计算域示意图,图2为圆柱周围的网格划分。 3 计算与分析 3.1 实验结果与实验结果的比较 表1为计算所得的斯特劳哈尔数(St)以及阻力系数(Cd)与文献中实验结果的对比。 计算结果与实验结果吻合良好,说明采用本文的数值计算方法是准确而可靠的。 3.2 u3000学习前后建设内自由振动的分布 图3～图8为圆柱绕流非稳态流动中一个典型的涡街脱落周期中不同时刻圆柱周围的压力场分布。反映了一个典型涡街脱落周期内漩涡的初生、形成、脱落及在尾流中不断发展的整个过程。从图3可看出,t=0时刻,上一周期脱落形成的涡在圆柱尾流中逐渐发展,圆柱的右下方约45°的位置形成一片明显的负压区,旋涡在负压区内逐渐生成。图4中t=1/5T时刻,圆柱右下方的漩涡完全生成,沿着尾流方向移动,并有逐步脱落的趋势。而图5中t=2/5T时刻,圆柱右下方形成的漩涡完全脱落并沿着尾流方向发展。图6中t=3/5T时刻,之前由圆柱右下方脱落的漩涡在尾流中逐渐发展,同时在圆柱的右上方约45°的位置形成一片负压区,旋涡在负压区内逐渐生成。图7中t=4/5T时刻圆柱右上方的漩涡完全生成,沿着尾流方向移动,并有逐渐脱落的趋势。最后圆柱右上