7037-各种湍流模型详细推导
- 格式:pdf
- 大小:639.68 KB
- 文档页数:34
第六章湍流模型湍流模型湍流运动中动量与能量交换主要受大尺度涡的影响湍流的基本方程无论湍流运动多么复杂,非稳态的连续方程和Navier-stokes 方程对于瞬时运动仍然是使用的。
对不可压流动:=01+=-+(grad )1+=-+(grad )1+=-+(grad )u p u v u t x v p v v t y w p w v w t zρρρ∇∂∂∇∇∂∂∂∂∇∇∂∂∂∂∇∇∂∂u u u u ()(v )()一、“雷诺平均”模式(RANS)根据湍流统计平均理论,湍流的速度、压强都可以分解为平均量和脉动量'i i iu u u=+p p p '=+其中,,i u p 为系综统计平均量,任意变量ф的时间平均值定义为:1()t ttt dt t φφ+∆=∆⎰,i u p ''为脉动量一、“雷诺平均”模式(RANS)对N-S 方程做系综平均()0i iu x ∂=∂遵循求导和系综平均可交换的原则,上式的线性项可直接写出:i iu u t t∂∂=∂∂21()i i i j i j i j ju u pu u f t x x x x νρ∂∂∂∂+=-++∂∂∂∂∂一、“雷诺平均”模式(RANS)对非线性对流项()()(()())()()i j i j j i i j i j i j i j i j j j j j i j i j ju u u u u u u u u u u u u u u u x x x x u u u u x ∂∂∂∂''''''==++=+++∂∂∂∂∂''=+∂将以上方程代入N-S 方程的系综平均中:'2'''''''2'''''''2=01+=-+(grad )+[---]1+=-+(grad )[---]1+=-+(grad )[---]u p u u v u w u v u t x x y z v p u v v v w v v v t y x y z w p u w v w w w v w t zx y z ρρρ∇∂∂∂∂∂∇∇∂∂∂∂∂∂∂∂∂∂∇∇+∂∂∂∂∂∂∂∂∂∂∇∇+∂∂∂∂∂u u u u ()()()()0i iu x ∂=∂21()()i i i j i j iji j j i u p u u u v u u f t x x x x x ρ∂∂∂∂∂''+=-+-+∂∂∂∂∂∂()ij i j R u u ρ''=-为雷诺应力项一、“雷诺平均”模式(RANS)()0i iu t x ρρ∂∂+=∂∂()1()[()]i i i j i j i ji j j u p u u u u u s t x x x x ρρμρρ∂∂∂∂∂''+=-+-+∂∂∂∂∂()()[()]j i j j i ju u s t x x x φρφρφρφ∂∂∂∂''+=Γ-+∂∂∂∂RANS方程和原N-S方程在形式上很相似,只是多了雷诺应力项(6个)。
湍流模型概述范文湍流是液体或气体流动中的一种复杂的现象,它包含着三个主要特征:不规则性、不可预测性和多尺度性。
湍流的产生是由于流体内部存在多个尺度的涡旋,它们之间相互作用并且不断地改变尺度和形状。
湍流模型是用来解释和描述湍流现象的一种数学方法。
本文将对湍流模型进行概述。
湍流模型的目的是通过对流场中各个参数的统计平均来描述湍流的性质。
根据湍流模型的复杂程度和适用范围的不同,可以将湍流模型分为三个等级:经验模型、半经验模型和基于数值模拟的模型。
经验模型是最早发展的湍流模型,它基于观察和实验结果,将各个参数之间的关系表示为一些经验公式。
这种模型的优点是简单易用,适用于几乎所有的湍流问题。
然而,由于经验模型只是基于经验规律,对于复杂的湍流现象并不准确,所以它的适用范围有限。
半经验模型是在经验模型的基础上发展而来的一种湍流模型。
它基于一些经验关系,并结合一些理论和模型来提高预测的准确性。
半经验模型通常可以包括一些常见的湍流模型,如k-ε模型、k-ω模型等。
这些模型通过引入一些方程和参数,来描述湍流的速度、压力和湍动能等参数之间的动态平衡。
半经验模型在计算流体力学领域得到了广泛的应用,能够较准确地预测湍流的性质。
基于数值模拟的模型是通过计算流体力学(CFD)方法来模拟湍流现象的模型。
这种模型基于流体的基本方程和湍流模型,通过数值方法进行求解并得到流场的数值解。
基于数值模拟的模型具有较高的计算精度和更大的适用范围,能够模拟各种复杂的湍流流动现象。
然而,基于数值模拟的模型需要较大的计算资源和时间,并且对于湍流模型的选择和设定需要一定的经验。
湍流模型的开发和改进是一个长期而具有挑战性的研究领域。
目前,湍流模型的研究主要集中在发展更准确、更适用于特定流动条件的模型。
研究人员通过理论推导、实验验证和数值模拟等方法,不断改进湍流模型的参数设定和方程形式,以提高湍流模型的准确性和适用性。
总之,湍流模型是描述和解释湍流现象的数学方法。
湍流模型方程
湍流模型方程是用来描述湍流流动的数学方程。
其中最经典的湍流模型方程是雷诺平均纳维-斯托克斯方程,也称为RANS方程。
雷诺平均纳维-斯托克斯方程是对流体流动进行平均处理后得到的方程,可以描述湍流的运动规律。
其方程形式如下:
∂(ρu_i)/∂t + ∂(ρu_iu_j)/∂x_j = - ∂p/∂x_i + ∂(τ_ij)/∂x_j + ρg_i + F_i
其中,ρ是流体的密度,u_i是速度分量,t是时间,x_i是空间坐标,p是压力,τ_ij是应力张量,g_i是重力分量,F_i是外力分量。
这个方程描述了流体的连续性、动量守恒和能量守恒。
湍流模型方程还包括了湍流模型,用来描述湍流的统计性质。
最常用的湍流模型是k-ε模型,它基于湍流运动的能量和湍流耗散率进行描述。
k-ε模型的方程如下:
∂(ρk)/∂t + ∂(ρku_i)/∂x_i = ∂(μ+μ_t)∂x_j ∂u_i/∂x_j - ρε + ρg_i + F_i
∂(ρε)/∂t + ∂(ρεu_i)/∂x_i = C_1εk/μ (∂(μ+μ_t)∂x_i ∂u_i/∂x_j) - C_2ρε^2/k + ρg_iu_i + F_i
其中,k是湍流能量,ε是湍流耗散率,μ是动力粘度,μ_t是湍流粘度,C_1和C_2是经验常数。
这个模型方程描述了湍流能量和湍流耗散率的传输过程,可以用来计算湍流流动的各种统计量。
湍流模型理论§3.1 引言自然界中的实际流动绝大部分是三维的湍流流动,如河流,血液流动等。
湍流是流体粘性运动最复杂的形式,湍流流动的核心特征是其在物理上近乎于无穷多的尺度和数学上强烈的非线性,这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。
回顾计算流体力学的发展,特别是活跃的80年代,不仅提出和发展了一大批高精度、高分辨率的计算格式,从主控方程看相当成功地解决了Euler方程的数值模拟,可以说Euler方程数值模拟方法的精度已接近于它有效使用范围的极限;同时还发展了一大批有效的网格生成技术及相应的软件,具体实现了工程计算所需要的复杂外形的计算网格;且随着计算机的发展,无论从计算时间还是从计算费用考虑,Euler方程都已能适用于各种实践所需。
在此基础上,80年代还进行了求解可压缩雷诺平均方程及其三维定态粘流流动的模拟。
90年代又开始一个非定常粘流流场模拟的新局面,这里所说的粘流流场具有高雷诺数、非定常、不稳定、剧烈分离流动的特点,显然需要继续探求更高精度的计算方法和更实用可靠的网格生成技术。
但更为重要的关键性的决策将是,研究湍流机理,建立相应的模式,并进行适当的模拟仍是解决湍流问题的重要途径。
要反映湍流流场的真实情况,目前数值模拟主要有三种方法:1.平均N-S方程的求解,2.大涡模拟(LES),3.直接数值模拟(DNS)。
但是由于叶轮机械内部结构的复杂性以及目前计算机运算速度较慢,大涡模拟和直接数值模拟还很少用于叶轮机械内部湍流场的计算,更多的是通过求解平均N-S方程来进行数值模拟。
因为平均N-S方程的不封闭性,人们引入了湍流模型来封闭方程组,所以模拟结果的好坏很大程度上取决于湍流模型的准确度。
自70年代以来,湍流模型的研究发展迅速,建立了一系列的零方程、一方程、两方程模型和二阶矩模型,已经能够十分成功的模拟边界层和剪切层流动。
但是,对于复杂的工业流动,比如航空发动机中的压气机动静叶相互干扰问题,大曲率绕流,激波与边界层相互干扰,流动分离,高速旋转以及其他一些原因,常常会改变湍流的结构,使那些能够预测简单流动的湍流模型失效,所以完善现有湍流模型和寻找新的湍流模型在实际工作中显得尤为重要。
湍流模型公式湍流模型公式是研究复杂流场的理论模型,由美国科学家亨利莱茨史密斯(Henry L. Smith)在1941年提出,它可以有效地用于模拟和研究湍流性质。
湍流模型公式也被称为史密斯-史密斯-贝尔(SSB)湍流模型,它是要全面描述湍流的最常用的模型之一。
湍流模型公式的核心是以下一维湍流方程:u_t+(uU)_x=νu_xx+ f(x,t) (1)其中u(x,t)为流体速度,U(x,t)为水流平均速度,ν为流的粘性系数,f(x,t)为外力。
另外,有必要提到的是,湍流模型公式还包括一维流体剪切应力方程和一维湍流能量方程:τ_t+(uτ)_x =ντ_xx+συw (2)E_t+(uE)_x =νE_xx+σu^2w (3) 其中τ(x,t)为流体剪切应力,E(x,t)为湍流能量,σ为流的粘性系数,u(x,t)为流体速度,U(x,t)为水流平均速度,υw为对流湍流损失(eddy viscosity)。
湍流模型公式最引人注目的特点是它准确地描述了湍流动力学的复杂性,因此它可以用来模拟流动场中的湍流运动,以提出解决实际技术问题的方案。
因此,湍流模型公式已广泛应用于工程和科学领域,如水资源开发、船舶控制、海洋工程、航空设计、气动学等。
湍流模型公式的求解用到了数值方法,即有限差分法。
它可以将湍流动力学的三维湍流问题简化为一维湍流问题,然后使用有限求差分法(FD)或有限元法(FEM)对参数进行有效求解。
其主要过程如下:1. 使用有限差分法或有限元法对定义的流动模型进行离散;2.用叙述方法或迭代方法求解湍流模型公式;3. 使用反演方法或最优控制方法优化参数;4.析和模拟湍流的空间和时间变化及其特性。
湍流模型公式的应用范围越来越广泛,它对于实际工程外面的许多流体动力学问题具有重要的理论意义和工程意义,这些问题包括流体运动的抗阻,水动力学的结构优化设计,湍流动态特性,空气动力学的优化控制,结构调节系数的估计,气动学的分析等等。
湍流模型介绍因为湍流现象是高度复杂的,所以至今还没有一种方法能够全面、准确地对所有流动问题中的湍流现象进行模拟。
在涉及湍流的计算中,都要对湍流模型的模拟能力以及计算所需系统资源进行综合考虑后,再选择合适的湍流模型进行模拟。
FLUENT 中采用的湍流模拟方法包括Spalart-Allmaras模型、standard(标准)k −ε模型、RNG(重整化群)k −ε模型、Realizable(现实)k −ε模型、v2 −f 模型、RSM(Reynolds Stress Model,雷诺应力模型)模型和LES(Large Eddy Simulation,大涡模拟)方法。
7.2.1 雷诺平均与大涡模拟的对比因为直接求解NS 方程非常困难,所以通常用两种办法对湍流进行模拟,即对NS 方程进行雷诺平均和滤波处理。
这两种方法都会增加新的未知量,因此需要相应增加控制方程的数量,以便保证未知数的数量与方程数量相同,达到封闭方程组的目的。
雷诺平均NS 方程是流场平均变量的控制方程,其相关的模拟理论被称为湍流模式理论。
湍流模式理论假定湍流中的流场变量由一个时均量和一个脉动量组成,以此观点处理NS 方程可以得出雷诺平均NS 方程(简称RNS 方程)。
在引入Boussinesq 假设,即认为湍流雷诺应力与应变成正比之后,湍流计算就归结为对雷诺应力与应变之间的比例系数(即湍流粘性系数)的计算。
根据计算中使用的变量数目和方程数目的不同,湍流模式理论中所包含的湍流模型又被分为二方程模型、一方程模型和零方程模型(代数模型)等大类。
FLUENT 中使用的三种k −ε模型、Spalart-Allmaras 模型、k −ω模型及雷诺应力模型RSM)等都属于湍流模式理论。
大涡模拟(LES)方法是通过滤波处理计算湍流的,其主要思想是大涡结构(又称拟序结构)受流场影响较大,小涡则可以认为是各向同性的,因而可以将大涡计算与小涡计算分开处理,并用统一的模型计算小涡。
由于航发燃烧室中的流动特性极其复杂,要想提高数值计算的预测能力,必须要慎重选择湍流模型。
用四种不同的湍流模型对带双径向旋流杯的下游流场进行数值模拟,将计算结果与实验结果作对比,比较各湍流模型的原理和物理基础,优劣,并分析流场速度分布和回流区特性。
涉及的湍流模型:标准k-ε湍流模型(SKE)1标准k-ε湍流模型有较高的稳定性,经济性和计算精度,应用广泛,适合高雷诺数湍流,但不适合旋流等各向异性较强的流动。
2简单的湍流模型是两个方程的模型,需要解两个变量,即速度和长度。
在fluent中,标准k-ε湍流模型自从被Launderand Spalding 提出之后,就变成流场计算中的主要工具。
其在工业上被普遍应用,其计算收敛性和准确性都非常符合工程计算的要求。
3但其也有某些限制,如ε方程包含不能在壁面计算的项,因此必须使用壁面函数。
另外,其预测强分离流,包含大曲率的流动和强压力梯度流动的结果较弱。
它是个半经验的公式,是从实验现象中总结出来的。
动能输运方程是通过精确的方程推导得到,耗散率方程是通过物理推理,数学上模拟相似原型方程得到的。
应用范围:该模型假设流动为完全湍流,分子粘性的影响可以忽略,此标准κ-ε模型只适合完全湍流的流动过程模拟。
可实现的k-ε模型是才出现的,比起标准k-ε模型来有两个主要的不同点:·可实现的k-ε模型为湍流粘性增加了一个公式。
·为耗散率增加了新的传输方程,这个方程来源于一个为层流速度波动而作的精确方程。
术语“realizable”,意味着模型要确保在雷诺压力中要有数学约束,湍流的连续性。
应用范围:可实现的k-ε模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。
而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。
可实现的k-ε模型和RNG k-ε模型都显现出比标准k-ε模型在强流线弯曲、漩涡和旋转有更好的表现。
由于带旋流修正的k-ε模型是新出现的模型,所以还没有确凿的证据表明它比RNGk-ε模型有更好的表现。