一个基于双重尺度展开的重正化群雷诺应力二阶矩模型.docVIP

一个基于双重尺度展开的重正化群雷诺应力二阶矩模型利用尺度参数将脉动场与平均场联系起来，并利用YOrszag的湍流重正化群理论波数脉动分量性能更可靠的预报性。 关键词：湍流模式，雷诺应力，二阶矩模式，重正化群，双重尺度xhwang@ustc.edu.cn 电话：0551-3607224 1、引言 近年来，湍流二阶矩封闭模式取得了很大的发展，如Shih和Lumley提出了二阶非线性压力-应变率快速项模型，Launder、Reece和Rodi提出了著名的LRR模型，符松等提出了FLT二阶矩模型成功地描述了湍流的非线性变形等。二阶矩模式中含有丰富的物理信息，。但是，湍流模式都是通过简单的量纲分析和数学物理性质的合理性决定，模式常数也基本通过实验数据的因此，从理论出发推导湍流二阶矩封闭模型对于认识湍流和完善不同特性湍流的封闭模式具有重要的现实意义。 Yoshizawa在Kraichnan的直接相互作用理论基础上提出了双重尺度展开法——TSDIA。Yoshizawa考虑到湍流场中平均场和脉动场的尺度不一样，因此引入双重尺度，并假设平均场为大尺度缓变运动和脉动场为小尺度快变运动，利用Kraichnan的直接相互作用理论对湍流方程进行封闭。Yakhot和Orszag较系统地利用重正化群方法分析湍流场，建立了湍流重正化群K-RNG K-）模型RNG K-模型已被Fluent等计算流体软件。Rubinstein和Barton利用重正化群理论对雷诺应力的本构关系进行分析，并对雷诺应力二阶矩模式进行模拟封闭，但是他们的中存在着一些不自洽的地方，。YO理论非线性雷诺应力模型中没有，这YO理论本身的局限性，YO理论中并没有很好的区分平均场脉动场，脉动速度各向同性。Yoshizawa的双重尺度展开区分的平均场脉动场，利用尺度参数将脉动场与平均场联系起来，并利用YO的湍流重正化群理论波数脉动分量，对雷诺应力方程中各未知复杂项进行模拟，计算各常数，从湍流二阶矩封闭模式的控制方程满足 (1) Yoshizawa认为由于湍流平均场和脉动场具有不同的空间和时间尺度，故需引入双重尺度：，用以区分湍流平均场和脉动场 [4,5]： (2) 其中大写的和表示平均场的空间和时间尺度，小写的和表示脉动场的空间和时间尺度，。并将脉动量围绕其各向同性部分展开如下： (3) 这里分别表示脉动速度和脉动压力的各向同性部分，则分别表示脉动速度和脉动压力场由于剪切平均场导致的各向异性部分。把公式(2)和(3)代入方程(1)中,可得各阶脉动速度分量的偏微分方程。为简化起见，我们忽略脉动量关于大尺度量和的导数项，为了消除对流项的卷带作用，引入Taylor假设，即作Galilean变换：，整理后的各阶脉动速度分量为 (4) (5) 为使湍流各向同性部分处于统计定常状态，需在方程(4)右边引入随机力，该随机力可以看作大尺度通过NS方程中非线性项对小尺度的有效作用。假设随机力满足Gauss分布，并具有两点关联： (6) 利用YO重正化群方法对方程(4)在Fourier空间进行分析，可得[6,10] 　　(7) 其中，，，，，为重正化微扰展开参数，不失普遍性令。结合方程(6)和(7)，随机脉动速度各向同性部分的两点关联近似满足 (8) 脉动速度场的各向异性部分满足方程： 　 (9) 　(10) 本文将从湍流雷诺应力输运方程出发，利用方程(8)-(10) 建立雷诺应力二阶矩封闭模型。 湍流雷诺应力输运方程为 (11) 雷诺应力二阶矩模型就是对方程中复杂未知项建立模型。我们先将各未知项利用Yoshizawa建议的双尺度作展开，然后对各展开项利用YO湍流重正化群理论进行计算模拟，并利用尺度关系将平均场与脉动场联系起来，由于具体数学计算推导过程较复杂，在本文中我们只直接给出推导结果。对于附加对流扩散项，模拟结果为： (12) 压力梯度与速度关联项可以被模拟成 (13) 对于耗散项的模拟，我们可以把它与雷诺应力中各向异性部分联系起来，并得到 (1４) 把方程(12)-(14)代入雷诺应力输运方程(11)中可得： (15) 方程(15)中: （16） 其中,, 方程(16)中理论计算得到的模型参数如下 ,，，，　(17) 3、数值算例比较分析 我们对新的RNG二阶矩模型利用两个简单的算例分析其有效性，选取充分发展的定常二维旋转槽道流动验证湍流模型对于旋转湍流的预报性；选取二维后台阶流动验证湍流模型对于湍流二次流的预报性。这两个算例虽然简单，但都具有很强的工程应用背景，其中旋转槽道湍流流动在涡轮机械装置中很常见，而后台阶流动则代表了燃烧工程、水利工程中很普遍的突扩流动。下面给出了新的RNG二阶矩模型，Gibson－Launder雷诺应力微分模型(GL模型)和