创作编号：GB8878185555334563BT9125XW
创作者： 凤呜大王*
大多数飞行器都是在高Re数下飞行，表面的流态是湍流。为了准确地确定湍流流态下的摩阻、热流，湍流成为一个重要而困难的研究课题。 （一）DNS
目前处理湍流数值计算问题有三种方法，第一种方法即所谓直接数值模拟方法（DNS方法），直接求解湍流运动的N-S方程，得到湍流的瞬时流场，即各种尺度的随机运动，可以获得湍流的全部信息。随着现代计算机的发展和先进的数值方法的研究，DNS方法已经成为解决湍流的一种实际的方法。但由于计算机条件的约束，目前只能限于一些低Re数的简单流动，不能用于工程应用。目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数(Re～200)和非常简单的流动外形，如平板边界层、完全发展的槽道流，以及后台阶流动等。用直接数值模拟方法处理工程中的复杂流动问题，即使是当前最先进的计算机也还差三个量级。透析液
（二）LES
另一种方法称做大涡模拟方法（LES方法）。这是一种折衷的方法，即对湍流脉动部分直接地模拟，将N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。小涡对大涡的影响会出现在大涡方程中，再通过建立模型（亚格子尺度模型）来模拟小涡的影响。由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件，难以出普遍的湍流模型来描述具有不同的边界特征的大涡结构，宜做直接模拟。相反地，小尺度涡对边界条件不存在直接依赖关系，而且一般具有各向同性性质。所以亚格子模型具有更大的普适性，比较容易构造，这是它比雷诺平均方法要优越的地方。自从1970年Deardorff第一次给出具有工程意义的LES计算以来，LES方法已经成为计算湍流的最强有力的工具之一，应用的方向也在逐步扩展，但是仍然受计算机条件等的限制，使之成为解决大量工程问题的成熟方法仍有很长的路要走。
（三）RANS
目前能够用于工程计算的方法就是模式理论。所谓湍流模式理论，就是依据湍流的理论知识、实验数据或直接数值模拟结果，对Reynolds应力做出各种假设，即假设各种经验的和 龙摄天下摄影团
半经验的本构关系，从而使湍流的平均Reynolds方程封闭。随着计算流体力学的发展，湍流模式理论也有了很大的进步，有了非常丰硕的成果。从对模式处理的出发点不同，可以将湍流模式理论分类成两大类：一类称为二阶矩封闭模式，另一类称涡粘性封闭模式。
（1）雷诺应力模式
所谓二阶矩封闭模式，是从高中英语写作教学Reynolds应力满足的方程出发，将方程右端未知的项（生成项，扩散项，耗散项等）用平均流动的物理量和湍流的特征尺度表示出来。典型的平均流动的变量是平均速度和平均温度的空间导数。这种模式理论，由于保留了Reynolds应力所满足的方程，如果模拟的好，可以较好地反映Reynolds应力随空间和时间的变化规律，因而可以较好地反映湍流运动规律。因此，二阶矩模式是一种较高级的模式，但是，由于保留了Reynolds应力的方程，加上平均运动的方程整个方程组总计15个方程，是一个庞大的方程组，应用这样一个庞大的方程组来解决实际工程问题，计算量很大，这就极大地限制了二阶矩模式在工程问题中的应用。
（2）涡粘性模式
天津市财政地税网
在工程湍流问题中得到广泛应用的模式是涡粘性模式。这是由Boussinesq仿照分子粘性的思路提出的，即设Reynolds应力为，
（）
这里是湍动能，称为涡粘性系数，这是最早提出的基准涡粘性模式，即假设雷诺应力与平均速度应变率成线性关系，当平均速度应变率确定后，六个雷诺应力只需要通过确定一个涡粘性系数就可完全确定，且涡粘性系数各向同性，可以通过附加的湍流量来模化，比如湍动能k，耗散率，比耗散率ω以及其它湍流量，，，根据引入的湍流量的不同，可以得到不同的涡粘性模式，比如常见的，k-w模式，以及后来不断得到发展的非接触测量，q-w，k-l等模式，涡粘性系数可以分别表示为
，，，，
为了使控制方程封闭，引入多少个附加的湍流量，就要同时求解多少个附加的微分方程，根据求解的附加的微分方程的数目，一般可将涡粘性模式划分为三类：零方程模式，半方程模型，一方程模式，两方程模式。
1）零方程模式
所谓零方程模式是试图直接用平均流动物理量模化，而不引入任何湍流量（如等）。例如，Prandttl的混合长理论就是一种零方程模式：
（5.7）
式中称为混合长。
在零方程模式的框架下，得到最为广泛应用的是Baldwin-Lomax模式[22]。该模式是对湍流边界层的内层和外层采用不同的混合长假设。这是因为靠近壁面处，湍流脉动受到很大的抑制，含能涡的尺度减小很多，因此长度尺度减小很多；另一方面，在边界层外缘，湍流呈间歇状，质量、动量和能量的输运能力大大下降，即湍流的扩散能力减小。这样，应用混合长理论来确定涡粘性系数在这两个不同的区域应该有不同的形式。Baldwin-Lomax模式的具体数学描述如下。
（5.8）
这里是的离壁面最小距离值。
对于内层，即，有
（5.9）
是涡量，是长度尺度
（5.10）
其中k=0.4是Karman常数，A+是模化常数，是无量纲法向距离：
而是摩擦速度，其含义为
此处下标表示壁面。
对于外层，即，有在母亲心里流浪
（5.11）
其中
是下列函数的最大值：
而是达到最大值的位置。是所谓的Klebanoff间歇函数：
是平均速度分布中最大值和最小值之差。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议