一种新型亚格子尺度模型建立方法与流程

本发明属于计算流体力学领域，具体涉及一种新型亚格子尺度模型建立方法。

背景技术：

计算流体力学(cfd)是流体力学、计算数学、计算机的交叉学科，采用计算机对流体动力学方程进行模拟，获得流体运动的力、热、频率等信息，从而为相关工业设计提供数据支持。随着计算机技术的发展，计算流体力学在航空航天、交通运输、化工、机械、能源等领域正发挥着越来越重要的作用。

流体流动分为层流和湍流两种状态。实际的流动基本都是湍流，或至少包含部分湍流。准确预测流体运动的关键是湍流模拟技术。当前的湍流模拟技术包括：雷诺平均方程(rans)方法，大涡模拟(les)方法和直接数值模拟(dns)方法。其中，雷诺平均方法对计算资源需求较少，但准度较低。直接数值模拟方法准度最高，但计算开销极大，目前主要局限于简单的学术问题。大涡模拟方法，以及部分区域采用雷诺平均的les/rans混合方法，能够以目前的计算能力大大提升复杂湍流的模拟准度，正在迅速渗透到各个设计部门，解决了许多之前难以处理的复杂湍流问题。

湍流大涡模拟方法的基本思想是对流场中大于网格尺度的流动(“大涡”)采用控制方程直接计算求解，而对小于网格尺度的流动(“小涡”)采用亚格子尺度模型模拟。涡粘性模型通过构造湍流粘性系数来模拟亚格子尺度运动，是应用最为广泛的亚格子尺度模型。而涡粘性亚格子模型应用最广泛的是smagorinsky模型。该模型采用局部速度场的变形率张量模值作为构造涡粘性系数的特征物理量。由于变形率张量的分布与湍流结构分布的相关性很低(一方面，在很多没有湍流的区域也存在较大变形率张量模值，另一方面，变形率张量模值大的地方，湍流强度可能很弱)，使得现有模型在非湍流区域，或转捩区域产生过大的粘性耗散，从而导致湍流推迟出现或者不合理地消退。

技术实现要素：

针对现有技术中的上述不足，本发明提供的一种新型亚格子尺度模型建立方法解决了现有模型在非湍流区域和转捩区域产生过大的非物理耗散的问题。

为了达到上述发明目的，本发明采用的技术方案为：一种新型亚格子尺度模型建立方法，包括以下步骤：

s1、基于涡拉伸函数得出涡粘性系数定义；

s2、依据涡粘性系数得出涡粘性亚格子尺度模型；

s3、依据涡粘性系数得出动态涡粘性亚格子尺度模型。

进一步地：所述步骤s1中涡粘性系数的表达式为：

cvsv＝0.25

上式中，和c均为涡粘性系数，cvsv为模型常数，为滤波尺度或局部网格单元的尺度，为流场涡量矢量：为流场速度矢量，为速度场变形率张量，为涡拉伸函数，t表示矩阵转置。

进一步地：所述步骤s2中涡粘性亚格子尺度模型为：

cvsv＝0.25

上式中，为涡粘性湍流应力张量，为变形率张量的分量。

进一步地：所述步骤s3中动态涡粘性亚格子尺度模型为：

上式中，lij，mij为中间变量。

进一步地：所述mij的计算公式为：

上式中，ˉ为一次滤波，～为二次滤波，为二次滤波尺度；

其中，为涡矢量的分量，为变形率张量的分量。

进一步地：所述lij的计算公式为：

上式中，ui，uj(i,j＝1,2,3)为速度矢量的分量。

本发明的有益效果为：本发明所提出的基于涡拉伸函数的涡粘性系数定义，以及由此得到的涡粘性亚格子尺度模型和动态涡粘性亚格子尺度模型，比基于变形率张量的传统涡粘性亚格子尺度模型具有更明确的物理意义，能够更为准确地刻画湍流涡粘性的分布。

附图说明

图1为本发明流程图；

图2为本发明中涡粘性系数与smagorinsky模型得到的涡粘性系数比较示意图；

图3为本发明中涡粘性亚格子尺度模型和动态涡粘性亚格子尺度模型模拟得到的湍流能量衰减以及湍流能谱分布与试验结果的比较示意图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

如图1所示，一种新型亚格子尺度模型建立方法，包括以下步骤：

s1、基于涡拉伸函数得出涡粘性系数定义；

涡粘性系数的表达式为：

cvsv＝0.25(i)

上式中，ˉ为一次滤波，为涡粘性系数，cvsv为模型常数，为滤波尺度或局部网格单元的尺度，为滤波后流场涡量矢量：为流场速度矢量，为滤波后速度场变形率张量，为涡拉伸函数，t表示矩阵转置。

s2、依据涡粘性系数得出涡粘性亚格子尺度模型；

涡粘性亚格子尺度模型为：

上式中，为涡粘性湍流应力张量，为变形率张量的分量。

在壁面解析大涡模拟(wrles:wallresolvedles)计算过程中，式(ii)需要在近壁面对滤波尺度作修正。可采用经典的vandriest修正：

上式中，y⁺为基于粘性尺度的无量纲壁面法向距离。

s3、依据涡粘性系数得出动态涡粘性亚格子尺度模型。

动态涡粘性亚格子尺度模型为：

上式中，lij和mij为中间变量，ˉ为一次滤波，～为二次滤波，为二次滤波尺度，为涡矢量的分量，为变形率张量的分量。

图2为本发明给出的vsv模型(式(i)和式(ii))与smagorinsky模型得到的涡粘性系数比较。从图2中可以看出传统的smagorinsky模型在非湍流区域给出了相当大的涡粘性系数，这是由于smagorinsky模型所基于对变形率张量在这些区域存在较大的量值所致。而本发明所给出的涡粘性系数则较好地捕捉到了流场中湍流能量级串区域的存在。

图3为本发明给出的vsv模型(式(ii))和dvsv模型(式(iii))模拟得到的湍流能量衰减以及湍流能谱分布e(k)与试验结果的比较。从图3中可以看出本发明给出的亚格子模型能够准确捕捉湍流能量衰减和能谱演化。