论文部分内容阅读
基于湍流模式的数值模拟在解决工程实际中湍流问题时发挥了重要的作用。然而这些常用的湍流模型主要是通过量纲分析和数学物理性质的合理性来实现的,模型中的湍流常数利用典型的湍流实验数据来标定。随着湍流模型在工程中的应用日益广泛,湍流模型的合理性、准确性面临着考验。因此从理论上来建立湍流模型,对人们研究和认识湍流有着重要的意义。基于湍流与临界现象的相似性,湍流研究工作者将用于分析相变临界现象的重正化群(Renormalization group)方法应用到湍流研究中来。Forster,Nelson和Stephen(FNS)最早将重正化群方法用来分析Navier-Stokes方程。而后,Yakhot和Orszag(YO)较系统地利用重正化群方法分析了各向同性湍流,并取得了一系列成果。Rubinstein和Barton等开创性地利用Yakhot和Orszag发展的湍流重正化群理论对Reynolds应力剪切模型进行了推导,然而他们的工作中存着一些不合理的地方。本文将利用湍流重正化群理论对一些常用的剪切湍流模型进行分析,主要包括Reynolds应力非线性涡粘性模型和二阶矩封闭模型的分析,并从理论上求得包含在这些模型中的湍流常数。在Rubinstein和Barton等工作的基础上,本文首先直接利用重正化群方法对Reynolds应力的涡粘性模型进行了分析,从理论上得到了一个二阶非线性涡粘性模型;计算上,本文通过积分波数的变换,采用单一的值y=3,即ε=4,回避了Rubinstein和Barton以及Yakhot等工作中的同一参数采用不同取值的不自洽问题。理论分析所得的非线性涡粘性模型中不仅含有平均应变的平方项(▽U)~2,还含有平均速度的二阶导数项U▽~2U,这一项可以看成Reynolds应力“松弛”项,能更好的反映湍流运动的非线性本质,但是,该模型非线性项中少了平均应变率张量与平均涡量张量的相互作用项,这是该理论模型的一个缺陷。本文还对Reynolds应力微分方程进行理论分析,对其中的惯性项、速度与压力梯度关联项以及耗散项予以模拟,并从理论上计算了各湍流常数,得到了完整的二阶矩封闭模型。结果表明重正化群二阶矩模型和著名的LRR模型很接近。需要指出的是,本文通过对Reynolds应力微分方程中耗散项的非线性部分模拟,得到二阶矩模型中的各向同性回归项。工程实际中Reynolds数为有限值的剪切湍流,其大涡运动在空间和时间上存在结构,不再满足自相似性。因此,直接将Yakhot—Orszag湍流重正化群理论推广到剪切湍流场分析时,其理论基础不完善。剪切湍流场中,湍流平均场是缓慢变化的,脉动场是快速变化的,它们具有不同的时间和空间尺度。因此,本文引入Yoshizawa采用的双重尺度展开,区分脉动场(快变量)和平均场(缓变量),并利用尺度比例常数对脉动场围绕其各向同性部分进行展开,得到脉动中的各向异性部分和各向同性的迭代关系,同时将剪切应变率引入脉动场各向异性分量中。在此基础上,利用Yakhot-Orszag湍流重正化群理论中的模式耦合原理,消除脉动场中的各向同性部分高波数分量,对Reynolds应力非线性涡粘性模型和二阶矩模型进行分析。引入双重尺度展开后,新的非线性模型比直接利用重正化群方法分析所得非线性模型有所改进,模型中含有平均应变率张量与平均涡量张量的相互作用项。新模型满足Speziale建议的合理性条件,不含有平均涡量的平方项。通过与其他模型的系数比较,本文所得理论模型与Speziale建议的非线性模型很接近;对于Reynolds应力二阶矩模型,主要通过对Reynolds应力输运方程中的三元速度关联项、脉动与压力梯度关联项以及耗散项的模拟对其进行封闭。本文同样是从耗散项的高阶项(各向异性项)中导出二阶矩模型中的Reynolds应力各向同性回归项。相比于直接利用湍流重正化群理论分析的剪切模型,引入双重尺度展开后,分析过程大大简化。两个典型湍流算例的数值结果表明,基于双重尺度展开所得的剪切模型比标准K-ε模型均有一定程度的改进,其中新的重正化群二阶矩模型具有相当于Gibson-Launder模型的精度。