基于k-两方程模式的湍流绕流流场数值模拟.docxVIP

基于k-两方程模式的湍流绕流流场数值模拟 1 大攻角特性分析 现代战争对战术导弹的射程和机动性提出了很高的要求。为了适应这种要求,先进战术导弹多采用小升力面甚至无翼布局,以减小巡航段的零升阻力,在需要过载时提高飞行攻角,充分利用大攻角时弹体的非线性涡升力输出所需的过载。飞行攻角增大到一定程度,由于三维边界层的分离,细长体背风面将形成分离涡,卷起形成特定的空间涡结构。背风面分离流动的特性对飞行器的气动性能有着重要的影响,是影响飞行器大攻角稳定性和操纵性的关键因素,给飞行器气动外形和控制系统的设计带来许多新问题。因此尽可能准确地预测前体大攻角特性具有很大的工程应用价值。 分离涡是三维边界层分离形成的,这一物理本质决定了对细长体背风面分离涡的准确模拟必须考虑粘性。此外,真实飞行器飞行雷诺数很高,可达到107～109量级,其绕流绝大部分呈湍流流态。因此细长体背风面分离流动模拟的准确度在很大程度上又与所采用的湍流模式有关。近几年来,国内外学者在湍流模式研究中发展了一系列的非线性涡粘性模式,对研究湍流的近壁特性和流场曲率影响效果显著。为探讨不同湍流模式对细长体分离流的模拟能力,作者曾选用不同湍流模式,针对尖头细长旋成体,研究过超音速中等大攻角下的背风面分离流。为了深入分析细长体绕流的物理特征以及湍流模式的适用性,本文对前阶段工作做了进一步的深化和细化,选用Shih根据可实现性条件提出的k-ε模式对不同速度下旋成体背风面分离流进行了数值模拟,从背风面分离涡的强度和位置、物面压力分布、集中力和力矩等方面与已有试验结果进行了对比,并分析了分离涡内的湍流特性。 2 数值模拟方法 2.1 u3000u3000e+ 对密度和压力采用Renolds平均,速度、能量和温度采用Favre平均,结合Boussinesque涡粘性假设,可压缩N-S方程可以写为 ?ρ?t+?ρuj?xj=0?ρui?t+?ρujui?xj=-?p?xi+?σij?xj?ρE?t+?ρujΗ?xj=??xj(uiσij-qj)}(1)?ρ?t+?ρuj?xj=0?ρui?t+?ρujui?xj=??p?xi+?σij?xj?ρE?t+?ρujH?xj=??xj(uiσij?qj)?????????????(1) 其中 σij=(μ+μt)(2Sij-23Sllδij)-23ρkδijqj=-γ(μΡr+μtΡrt)?e?xjE=e+12uiui+kΗ=h+12uiui+k}(2)σij=(μ+μt)(2Sij?23Sllδij)?23ρkδijqj=?γ(μPr+μtPrt)?e?xjE=e+12uiui+kH=h+12uiui+k???????????????(2) 式中e=CvT,h=CpT,p=(γ-1)ρ e。 用有限体积法对方程(1)进行空间离散,多步Runge-Kutta显式方法进行时间推进。 2.2 可实现性条件是构造通用湍流模式的一个极端重要和有效的约束条件,要求模拟后的量不应当产生物理上不可能的值,如负的正应力或湍动能,关联系数大于1等。Shih等人根据可实现性条件提出的k-ε模式为 ??t(ρk)+??xi(ρkui)=??xi[(μ+μtσk)?k?xi]+Ρk-ρε-YΜ(3)??t(ρk)+??xi(ρkui)=??xi[(μ+μtσk)?k?xi]+Pk?ρε?YM(3) ??t(ρε)+??xi(ρεui)=??xi[(μ+μtσε)?ε?xi]+ρC1Sε-ρC2ε2k+√νε(4)??t(ρε)+??xi(ρεui)=??xi[(μ+μtσε)?ε?xi]+ρC1Sε?ρC2ε2k+νε√(4) 式中C2是常数,σk和σε是对应于k和ε的Prandtl数,Pk和YM分别是湍动能生成项和压缩性影响项。湍流粘性定义为 μt=ρCμk2εμt=ρCμk2ε 其中 Ρk=-ρˉu′iu′jSij?YΜ=2ρεΜ2tC1=max(0.43?ηη+5)?η=SkεC2=1.9?σk=1.0?σε=1.2Cμ=1A0+ASU(*)k/ε 2.3 k-模式下的k-模式 定义湍流雷诺数Rey Rey=ρ√ky/μ 根据Rey的大小将壁面附近的计算域分为两层:Rey200时认为是完全湍流区,采用两方程k-ε模式;而当Rey200时,采用Wolfstein提出的近壁处理方法,只求解k方程,分别用关系式 μt=ρCμ√klμ?ε=k32/lε 计算湍流粘性系数和湍流耗散率ε,其中的长度尺度取为 lμ=C1y[1-e(-ReyAμ)]?lε=C1y[1-e(-ReyAε)] 式中的常数取为 C1=κC-3/4μ,Au=70,Aε=2C1 3 流场的对称面 文献提供了一个尖拱形旋成体的试验结果。试验模型是13倍长细比的尖拱形旋成体,头部长细比等于3,柱段长细比等于10。