第10章 湍流模型.pdf
《第10章 湍流模型.pdf》由会员分享,可在线阅读,更多相关《第10章 湍流模型.pdf(55页珍藏版)》请在咨信网上搜索。
1、第十章湍流模型本章主要介绍Fluent所使用的各种湍流模型及使用方法。各小节的具体内容是:10.1 简介10.2 选择湍流模型10.3 Spalart-Allmaras 模型10.4 标准、RNG和k-e相关模型10.5 标准和SSTk-a模型10.6 雷诺应力模型10.7 大涡模拟模型LEN10.8 边界层湍流的近壁处理10.9 湍流仿真模型的网格划分10.10 湍流模型的问题提出10.11 湍流模型问题的解决方法10.12 湍流模型的后处理10.1 简介湍流出现在速度变动的地方。这种波动使得流体介质之间相互交换动量、能量和浓度变 化,而且引起了数量的波动。由于这种波动是小尺度且是高频率的,
2、所以在实际工程计算中 直接模拟的话对计算机的要求会很高。实际上瞬时控制方程可能在时间上、空间上是均匀的,或者可以人为的改变尺度,这样修改后的方程耗费较少的计算机。但是,修改后的方程可能 包含有我们所不知的变量,湍流模型需要用已知变量来确定这些变量。FLUENT提供了以下湍流模型:,Spalart-Allmaras 模型,k-e模型标准k-e模型Renormalization-group(RNG)k-e 模型带旋流修正k-e模型 k-s模型标准k-3模型压力修正k-。模型 雷诺兹压力模型 大漩涡模拟模型10.2 选择一个湍流模型不幸的是没有一个湍流模型对于所有的问题是通用的。选择模型时主要依靠以
3、下几点:流体是否可压、建立特殊的可行的问题、精度的要求、计算机的能力、时间的限制。为了选 择最好的模型,你需要了解不同条件的适用范围和限制这一章的目的是给出在FLUENT中湍流模型的总的情况。我们将讨论单个模型对cpu 和内存的要求。同时陈述一下一种模型对那些特定问题最适用,给出一般的指导方针以便对 于你需要的给出湍流模型。10.2.1 雷诺平均逼近vsLES在复杂形体的高雷诺数湍流中要求得精确的N-S方程的有关时间的解在近期内不太可 能实现。两种可选择的方法用于把N-S方程不直接用于小尺度的模拟:雷诺平均和过滤。两种方法都介绍了控制方程的附加条件,这些条件用于使模型封闭(封闭意味着有足够的方
4、 程来解所有的未知数。)对于所有尺度的湍流模型,雷诺平均N-S方程(The Reynolds-averaged Navier-Stokes(RANS)只是传输平均的数量。找到一种可行的平均流动变量可以大大的减少计算机的工 作量。如果平均流动是稳态的,那么控制方程就不必包含时间分量,并且稳态状态解决方法 会更加有效。甚至在暂态过程中计算也是有利的,因为时间步长在平均流动中取决于全局的 非稳态。雷诺平均逼近主要用于实际工程计算中,还有使用的模型比如Spalart-Allmaras,k-e系列,k-3系列和RSM。LES提供了一种方式,让依靠时间尺度模拟的大边界计算问题可以利用一系列的过滤方 程。对
5、于解确切的N-S方程,过滤是一种必要的方法,用于改变比过滤法尺度小的边界通 常用于网格大小。和雷诺平均一样,过滤法加入了未知的变量,必须模拟出来以便方程能够 封闭。必须强调的是LES应用于工业的流产模拟还处于起步阶段。回顾近期的出版物,典型的 方法已经用于简单的几何形体。这主要是因为解决含有能量的湍流漩涡需要大量的计算机资 源。很多成功的LES模型已经用于高度空间的离散化,而且花了很多精力来解决尺度比惯性 附属区域大的方面。在中间流中用LES降低精度的方法没有很多的资料。另外,用LES解 决平板问题还需要进一步的证实。作为一个一般性的介绍,在这里推荐一般的湍流模型用雷诺平均对于实际的计算是十分
6、 有用的。在10.7中将会详细介绍的LES逼近,对你十分有用,如果你的计算机能力很强大 或者有意更新你的计算机的话。这一章余下的部分将会介绍选择雷诺平均逼近模型。10.2.2 雷诺平均在雷诺平均中,在瞬态N-S方程中要求的变量已经分解为时均常量和变量。以速度为 例:%=%+4.(10.2-1)这里耳和孙时时均速度和波动分量。相似的,像压力和其它的标量a.(10.22)这里,表示一个标量如压力,动能,或粒子浓度。用这种形式的表达式把流动的变量放入连续性方程和动量方程并且取一段一段时间的 平均,这样可以写成一下的形式:,十o(1023)d d万)+西(PW)=一患+5偿+需一施整+台一丽)(102
7、4)方程1023和1024称为雷诺平均N-S方程。它和瞬态雷诺方程又相同的形式,速度和其 它的变量表示成为了其时均形式。由于湍流造成的附加的条件现在表现出来了。这些雷诺压 力一鬲可,必须被模拟出来以便使方程10.2-4封闭。对于变密度的流体,方程10.2-3和10.2-4认为是Favre平均N-S方程,速度表示为了 平均值。这样,方程10.2-3和10.2-4可以应用于变密度的流体。10.2.3 Boussinesq逼近VS雷诺压力转化模型对于湍流模型,雷诺平均逼近要求在方程10.2-4的雷诺压力可以被精确的模拟。一般 的方法利用Boussinesq假设把雷诺压力和平均速度梯度联系起来:Bou
8、ssinesq假设使用在Spalart-Allmaras模型、k-e模型和k-。模型中。这种逼近方法好处是对计 算机的要求不高。在Spalart-Allmaras模型中只有一个额外的方程要解。k-e模型和k-s模型中又 两个方程要解。Boussinesq假设的不足之处是假设3是个等方性标量,这是不严格的。可选的逼近,在RSM中,是用来解决在方程中的雷诺压力张量。另外要加一个方程。这就意味着在二维流场中要加五个方程,而在三维方程中要加七个方程。在很多情况下基于Boussinesq假设的模型很好用,而且计算量并不是很大。但是RSM 模型对于对层流有主要影响的各向异性湍流的状况十分适用。10.2.4
9、 The Spalart-Allmaras 模型对于解决动力漩涡粘性,Spalart-Allmaras模型是相对简单的方程。它包含了一组新的 方程,在这些方程里不必要去计算和剪应力层厚度相关的长度尺度。Spalart-Allmaras模型 是设计用于航空领域的,主要是墙壁束缚流动,而且已经显示出和好的效果。在透平机械中 的应用也愈加广泛。在原始形式中Spalart-Allmaras模型对于低雷诺数模型是十分有效的,要求边界层中粘 性影响的区域被适当的解决。在FLUENT中,Spalart-Allmaras模型用在网格划分的不是很 好时。这将是最好的选择,当精确的计算在湍流中并不是十分需要时。再
10、有,在模型中近壁 的变量梯度比在k-e模型和k-3模型中的要小的多。这也许可以使模型对于数值的误差变得 不敏感。想知道数值误差的具体情况请看5.1.2。需要注意的是Spalart-Allmaras模型是一种新出现的模型,现在不能断定它适用于所有 的复杂的工程流体。例如,不能依靠它去预测均匀衰退,各向同性湍流。还有要注意的是,单方程的模型经常因为对长度的不敏感而受到批评,例如当流动墙壁束缚变为自由剪切流。10.2.5 标准k-醴型最简单的完整湍流模型是两个方程的模型,要解两个变量,速度和长度尺度。在FLUENT 中,标准k-e模型自从被Launder and Spalding提出之后,就变成工程
11、流场计算中主要的工具了。适用范围广、经济、合理的精度,这就是为什么它在工业流场和热交换模拟中有如此广泛的 应用了。它是个半经验的公式,是从实验现象中总结出来的。由于人们已经知道了k-e模型适用的范围,因此人们对它加以改造,出现了RNGk-e模型 和带旋流修正k-e模型 10.2.6 RNGk-醴型RNGk-e模型来源于严格的统计技术。它和标准k-e模型很相似,但是有以下改进:RNG模型在e方程中加了一个条件,有效的改善了精度。考虑到了湍流漩涡,提高了在这方面的精度。RNG理论为湍流Prandtl数提供了一个解析公式,然而标准k-e模型使用的是用户提供的 常数。然而标准k-e模型是一种高雷诺数的
12、模型,RNG理论提供了一个考虑低雷诺数流动粘性 的解析公式。这些公式的效用依靠正确的对待近壁区域这些特点使得RNG k-e模型比标准k-e模型在更广泛的流动中有更高的可信度和精度。10.2.7 带旋流修正的k-醴型带旋流修正的k-e模型是近期才出现的,比起标准k-e模型来有两个主要的不同点。带旋流修正的k-e模型为湍流粘性增加了一个公式。为耗散率增加了新的传输方程,这个方程来源于一个为层流速度波动而作的精确方程 术语realizable,意味着模型要确保在雷诺压力中要有数学约束,湍流的连续性。带旋流修正的k-e模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。而且它对于旋转流动、强逆
13、压梯度的边界层流动、流动分离和二次流有很好的表现。带旋流修正的k-e模型和RNG k-e模型都显现出比标准k-e模型在强流线弯曲、漩涡和 旋转有更好的表现。由于带旋流修正的k-e模型是新出现的模型,所以现在还没有确凿的证 据表明它比RNG k-e模型有更好的表现。但是最初的研究表明带旋流修正的k-e模型在所有 k-e模型中流动分离和复杂二次流有很好的作用。带旋流修正的k-e模型的一个不足是在主要计算旋转和静态流动区域时不能提供自然 的湍流粘度。这是因为带旋流修正的k-e模型在定义湍流粘度时考虑了平均旋度的影响。这 种额外的旋转影响已经在单一旋转参考系中得到证实,而且表现要好于标准k-e模型。由
14、于 这些修改,把它应用于多重参考系统中需要注意。10.2.8 标准k-s模型标准卜。模型是基于Wilcoxk-s模型,它是为考虑低雷诺数、可压缩性和剪切流传播而 修改的。Wilcoxk模型预测了自由剪切流传播速率,像尾流、混合流动、平板绕流、圆柱 绕流和放射状喷射,因而可以应用于墙壁束缚流动和自由剪切流动。标准k-e模型的一个变 形是SSTk-s模型,它在FLUENT中也是可用的,将在1029中介绍它。10.2.9 剪切压力传输(SST)k-3模型SSTk-s模型由Menter发展,以便使得在广泛的领域中可以独立于k-e模型,使得在近壁 自由流中k-3模型有广泛的应用范围和精度。为了达到此目的
15、,k-e模型变成了k-s公式。SST k-3模型和标准k-。模型相似,但有以下改进:SSTky模型和k-e模型的变形增长于混合功能和双模型加在一起。混合功能是为近壁 区域设计的,这个区域对标准上。模型有效,还有自由表面,这对k-e模型的变形有效。SSTk-s模型合并了来源于。方程中的交叉扩散。湍流粘度考虑到了湍流剪应力的传波。模型常量不同这些改进使得SST k模型比标准卜模型在在广泛的流动领域中有更高的精度和可信 度。10.2.10 雷诺压力模型(RSM)在FLUENT中RSM是最精细制作的模型。放弃等方性边界速度假设,RSM使得雷诺平均 N-S方程封闭,解决了关于方程中的雷诺压力,还有耗散速
16、率。这意味这在二维流动中加入 了四个方程,而在三维流动中加入了七个方程。由于RSM比单方程和双方程模型更加严格的考虑了流线型弯曲、漩涡、旋转和张力快速 变化,它对于复杂流动有更高的精度预测的潜力。但是这种预测仅仅限于与雷诺压力有关的 方程。压力张力和耗散速率被认为是使RSM模型预测精度降低的主要因素。RSM模型并不总是因为比简单模型好而花费更多的计算机资源。但是要考虑雷诺压力的 各向异性时,必须用RSM模型。例如飓风流动、燃烧室高速旋转流、管道中二次流。10.2.11 计算成效:cpu时间和解决方案从计算的角度看Spalart-Allmaras模型在FLUENT中是最经济的湍流模型,虽然只有
17、一种方程可以解。由于要解额外的方程,标准k-e模型比Spalart-Allmaras模型耗费更多的计算 机资源。带旋流修正的k-e模型比标准k-e模型稍微多一点。由于控制方程中额外的功能和非 线性,RNGk-e模型比标准k-e模型多消耗1015%的CPU时间。就像k-e模型,卜。模型也是两 个方程的模型,所以计算时间相同。比较一下k-e模型和k-模型,RSM模型因为考虑了雷诺压力而需要更多的CPU时间。然 而高效的程序大大的节约了CPU时间。RSM模型比k-e模型和k-s模型要多耗费5060%的CPU 时间,还有1520%的内存。除了时间,湍流模型的选择也影响FLUENT的计算。比如标准k-e
18、模型是专为轻微的扩散 设计的,然而RNGk-e模型是为高张力引起的湍流粘度降低而设计的。这就是RNG模型的缺点。同样的,RSM模型需要比k-e模型和k-s模型更多的时间因为它要联合雷诺压力和层流。10.3 Spalart-Allmaras 模型在湍流模型中利用Boussinesq逼近,中心问题是怎样计算漩涡粘度。这个模型被Spalart andAllmaras提出,用来解决因湍流动粘滞率而修改的数量方程。10.3.1 Spa/ar/maras膜型的偏微方程Spa/a”/加aras模型的变量中v是湍流动粘滞率除了近壁区域,方程是:K+S/(10.3-1)这里Gv是湍流粘度生成的,Yv是被湍流粘度
19、消去,发生在近壁区域。S一是用户定义的。注 意到湍流动能在Spalart-Allmaras没有被计算,但估计雷诺压力时没有被考虑。10.3.2 湍流粘度的建模湍流粘度W由以下公式计算:%=pifvi(10.3-2)储由下式:3fvl=3;10.3-3并且=(10.3-4)V10.3.3 湍流生产的建模Gv由下式Gy=CblPSi)(10.3-5)v zS=Sfv2 10.3-6fv2=1-,(10.3-7)1+xfvlCbi和k是常数,d是离墙的距离,S是变形张量。在FLUENT中,S由下式给出:S三河而(10.3-8)这里dj是层流旋转张量,由下式定义:(10.3-9)当模型给出时,我们最感
20、兴趣的是墙壁束缚流动中s表达式的修正,湍流漩涡只发生在近壁。但是,我们知道要把湍流产生的平均应变考虑进去,并且按照建议改变模型。这种修改包括旋度和应变,在s中定义:S=|5j|+Cprod min(10.3-10)Cprod=2。三,20*5 三,在平均应变率中Sij定义为:(10.3-11)包括旋度和应变张量减少了漩涡粘度从而减少了漩涡粘度本身。这样的例子可以在漩 涡流动中找到。旋度和应变张量更多正确的考虑湍流旋度。一般的方法是预测漩涡粘度的产 生并且预测漩涡粘度本身。你可以选择模型,在Viscous Model面板。10.3.4 湍流耗散的建模耗散的模型是:Yy=CwlPfw(10.3-1
21、2)tw=91+死心z;6 ii9 十 w3g=+CW2(r6-r(10.3-13)(10.3-14)vSK2d2(10.3-15)Cwi、Cw2和Cw3是常量,S由方程10.3-6给出。注意到考虑大平均应力而修改的S也会影响用S去计算厂。10.3.5 模型常量模型常量包括C72CV1,。小和k,下面是它们的值:9Qi=0.1335,Cb2=0.622.(jy=-?Cv=7.1 Jcwl=Ui+J_.Cw2=0.3,Cw3=2.0,K=0.4187 M(Jy10.3.6 墙壁边界条件在墙壁上,修改后的湍流动粘度,V,被认为是0。当网格划分的较好可以解决层状亚 层,壁面剪应力可以由下面的关系式得
22、出:U=叫一了 4(10.3-1G)如果网格太粗糙不足以解决,那么就假设工InEUT K(pUrVk 4(10.3-17)这里u是平行于壁面的速度,Ur是切速度,y是离墙壁的距离,k是von Karman常量E=9.793。10.3.6 热对流和质量转移模型在FLUENT中,湍流热交换使用的是对湍流动能交换的雷诺分析,能量方程如下:袅E)+岛k(pE+p)=/G+器)器+如(%任+Sh(10.3-18)k是导热系数,E是总能,Tq)ef是偏应力张量:dUj dUi dxi+a叼2 dm 铲F%T(ij)ef考虑到了由于粘性而产生的热,并且总是联合方程中。它在不能单个中解出,但是可 以在粘性模型
23、面板中找到。默认的湍流Prandtl数是0.85,你可以在粘性模型面板中改变它。湍流物质交换可以按照相似的方法,Schmidt数是0.7,可以在粘性模型面板中改变它。标量的墙壁边界条件可以类似于动量,可以用墙壁法则。10.4 标准、RNG和带旋流修正k-e模型这一章讲述标准、RNG和带旋流修正k-e模型这三种模型有相似的形式,有k方程和e方程,它们主要的不同点是:计算湍流粘性的方法 湍流Prandtl数由k和e方程的湍流扩散决定 在e方程中湍流的产生和消失每个模型计算湍流粘性的方法和模型的常数不一样。但从本质上它们在其它方面是一样的。10.4.1标准k-e模型标准k-e模型是个半经验公式,主要
24、是基于湍流动能和扩散率。k方程是个精确方程,e方程是个由经验公式导出的方程。k-e模型假定流场完全是湍流,分子之间的粘性可以忽略。标准k-e模型因而只对完 全是湍流的流场有效。标准k-e模型的方程湍流动能方程k,和扩散方程e:(pfc)十-(pA-Ui)=7 十J 十 GA 十 Gb +S*(10.4-1)and晶+总=就(+知蜀十%(Gk+QG)-二十 S6(10.4-2)K方程中Gk表示由层流速度梯度而产生的湍流动能,计算方法在10.4.4中有介绍。伍是由浮力 产生的湍流动能,10.4.5中有介绍,YM由于在可压缩湍流中,过渡的扩散产生的波动,10.4.6 中有介绍,Cl,C2,C3,是常
25、量,Ok和。e是k方程和e方程的湍流Prandtl数,Sk和Se是用户定 义的。湍流速度模型湍流速度ut由下式确定4之3P。一(10.4-3)G是常量模型常量Ce=1.44,C2c=1.92,C世=0.09,。卜=1.0.cre=1.3这些常量是从试验中得来的,包括空气、水的基本湍流。他们已经发现了怎样很好的处 理墙壁束缚和自由剪切流。虽然这些常量对于大多数情况是适用的,你还是可以在粘性模型面板中来改变它们。10.4.2 RNG k-e 模型RNG k-e模型是从暂态N-S方程中推出的,使用了一种叫“renormalization group”的 数学方法。解析性是由它直接从标准k-e模型变来
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第10章 湍流模型 10 湍流 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【曲****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【曲****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。