湍流直接数值模拟

第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
混淆误差产生原因： 当用伪谱方法计算时，谱空分量用，它等于：
利用三角级数公式
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
在上式中，p、q、k的取值范围都是{一N/2，N/2-1}，所 以取值不等于零的波数组合情况只有两种：p+q-k＝0，或 p+q-k＝ N，于是有
Kolmogorov耗散尺度
动速度均方根值)，将以上关系代入式，可得：

1/4 3


，而 ~ u '3 / l （ u ' 是脉
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率
那么三维总网格数N：
＝10４
这是一个天文数字的估计 ，假设 Rel ＝104，就要求网格数为109， 考虑到计算的流动变量数，需要约1010字长的计算机内存。直接数值 模拟实际工程湍流运动时，对网格分辨率的要求更高。 应当指出选定最小的网格长度还和数值方法有关。谱方法的数 值精度最高，差分法的精度和差分格式有关。
边界条件的提法： (3)渐近条件 对于湍流边界层或其他薄湍流切变层，在远离薄层和物面的渐近区 域，速度场趋近于无旋的均匀场，因此对于不可压缩流体可以采用如下：
刚盖假定：数值方法只能计算有限域内的流动，渐近条件只能采用 近似形式，一种方法是在离开薄层或物体横向一定距离的平面上设 置“虚拟边界”，在虚拟边界y＝H满足以下条件，称为刚盖假定：
式中，第一项是乘积w＝uv在谱空间中的投影，第二项是 在伪谱运算中产生的误差，故称混淆误差。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
混淆误差的消除--3/2规则 将函数在傅里叶空间中展开的系数进行延拓到原来的3/2倍， 延拓后函数的傅里叶系数有M＝3N/2项。延拓按以下的规则赋 值：
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.3 初始条件和边界条件
边界条件的提法： (3)渐近条件 另一种更好的方法是先做一个指数变换，将无限域变到有限域， 例如，令：
然后，在有限域里数值求解Navier－Stokes方程。如果y=0是 固壁，则在指数变换时，在y＝0附近自动加密网格，而在 方向 则是均匀网格。在(x，y，z)坐标系里，原渐近边界条件可写作：
(1)加权余量法
L u f (u)
L表示微分算子，f (u) 是已知函数。将未知函数 u 用一组完备的 线性独立函数族 展开： k
k 0,1
当展开式只取有限项时，上式是原函数 u 的近似。把uN代入原 来的微分方程，将产生误差，并称之为残差或余量，用RN表示：
第6章 湍流直接数值模拟
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率
以均匀各向同性湍流为例，假定各向同性湍流的含能尺度或积分 尺度为 l ，Kolmogorov耗散尺度等于 。为了足够准确地计算湍 流的大尺度运动，立方体的长度L必须大于含能尺度 l ，另一方 面，为了保证准确模拟湍流小尺度运动，网格长度 必须小于耗 散尺度 。因此，一维网格数至少应满足以下不等式：
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.3 平面槽道湍流的直接数值模拟
平面槽道湍流是典型的有固壁的简单湍流运动，如图6-8所示。流体 在平行平板之间流动，假定平板在流向和展向都是无限长。平均定常的 槽道湍流的流向平均压强梯度是常数，展向的平均压强梯度等于零。计 算过程中保持平均流量不变，因此以平均速度为特征长度的流动雷诺数 也保持常数。槽道宽度等于2H，计算域的长度和宽度根据流动雷诺数确 定，原则是计算域应当包含足够多的近壁结构。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
再根据-M/2≤p≤M/2-1和-M/2≤q≤M/2-1(M＝3N/2)的取值限制,满足 上式的p和q只有以下四种可能取值范围：
在以上p，q取值范围里，按照延拓公式有 U p Vq 0 ，因此按3/2延 拓规则,式(6．19)的第二项等于零，就是说，式(6．19)的最后结果 是：


第6章 湍流直接数值模拟
6.1 湍流数值模拟的方法
三者对比与联系： （1）直接数值模拟要求模拟所有尺度的湍流脉动；雷 诺平均方法网格尺度应当大于脉动的积分尺度；大涡数 值模拟的网格分辨率介于DNS和RANS之间。
（2）三种湍流数值模拟方法给出的信息量有很大差别。
（3）必须根据需要来选择数值模拟的方法。
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.3 平面槽道湍流的直接数值模拟
控制方程：
边界条件：槽道的上下壁面应是无滑移条件
结合实际可将速度和脉动压强展开如下：
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.3 平面槽道湍流的直接数值模拟
网格坐标为：
理论上，直接数值模拟要求网格最小长度达到Kolmogorov尺 度的量级，因此在给定算例的雷诺数后，应当估计需要的网格数， 以满足分辨率的要求。 首先，确定槽道湍流的雷诺数，Um是槽道截面的平均流速，H是 槽道的半宽度。然后估算Kolmogorov的耗散尺度和时间尺度
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
再选择另一组完备的线性独立函数族 k 作为权函数．要求余量 的加权积分等于零：
如果微分算子L是线性的，则最后求解的是线性代数方程组； 如果微分算子L是非线性的，则最后求解的是非线性代数方程组。 求出代数方程的解，就得到微分方程的近似解。 ①伽辽金方法(Galerkin法) 加权余量法三种形式 ②Tau方法 ③配置点法
以延拓的谱系数做伪谱运算，得到类似的乘积公式
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理
上面式中p,q,k的取值范围都是{-M/2,M/2-1}，所以，当时，由 延拓公式(6．18a)~(6．18d)可得式(6.19)的第一项,该项恰好是准确的 卷积，即：
由于M＝3N/2，式(6.18)中的第二项中p+q＝k 3N/2，当 N / 2 k N / 2 1 时，p+q的取值范围是：
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
6.2.1 湍流直接数值模拟的空间分辨率 6.2.2 湍流直接数值模拟的时间分辨率 6.2.3 初始条件和边界条件
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理
湍流是多尺度的不规则运动，湍流直接数值模拟和层 流运动的数值计算有很大区别：
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理 谱方法的优点是精度高，计算速度快 。但是适应复杂边界的 试探函数十分难找，所以对于复杂边界的湍流问题，特别是对于 流场中存在间断的情况，只能采用差分离散方法。
(2)伪谱法和混淆误差 混淆误差概念： 对于非线性方程，有限项谱展开的非线性项会产生附加的 误差，这种误差在谱方法中称做混淆误差(aliasingerror)。
a. 近似非定常出口条件：
(5)出口条 件 b. 嵌边区(或强粘性区)方法
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理 6.3.2 格栅湍流的直接数值模拟 6.3.3 平面槽道湍流的直接数值模拟
第6章 湍流直接数值模拟
6.3 湍流直接数值模拟的谱方法
6.3.1 谱方法的基本原理 谱方法是一种加权余量的数值计算方法，它把微分方程离散化 为代数方程组，其基本原理如下。设有微分方程：
数值模拟结果：各向同性湍流中 的涡结构算例
强迫各向同性湍流直接数值模 拟的结果表明，在各向同性湍流中 强涡是以细长涡管形式出现，如图 6-5所示。涡管的直径是Kolmogorov 尺度，涡管的平均长度是积分尺度， 涡管的强度(环绕涡管的环量)随雷诺 数的增加而增大。图中管状结构是 涡量的等值面，虽然它们并非真正 的涡管，但由于强涡量集中在很细 的管状结构中，可以推断它们接近 于当地的涡管。
第六章 湍流直接数值模拟
第6章 湍流直接数值模拟
6.1 湍流数值模拟的方法 6.2 湍流数值模拟的基本原理 6.3 湍流直接数值模拟的谱方法 6.4 湍流直接数值模拟的差分法
第6章 湍流直接数值模拟
6.1 湍流数值模拟的方法
直接数值模拟（DNS）： 从完全精确的流动控制方程出发，对所有尺度的湍流运动进行 数值模拟，这种最精细的数值模拟称为直接数值模拟。 雷诺平均数值模拟（RANS）： 从雷诺平均方程出发，在这一层次上的数值模拟称为雷诺平均 数值模拟。 大涡数值模拟（LES）： 介于DNS和RANS之间的数值模拟方法称为大涡数值模拟。思想 是：大尺度脉动(或大尺度湍涡)用数值模拟方法计算，只将小尺度 脉动对大尺度运动的作用做模型假设。LES的理论依据是小尺度脉 动有局部平衡的性质，很可能存在某种局部普适的统计规律，如局 部各向同性或局部相似性等。
第一，由于湍流脉动具有宽带的波数谱和频谱，因此湍流直接 数值模拟要求有很高的时间和空间分辨率。
第二，为了求得湍流统计特性，需要足够多的样本流动；如果 湍流是时间平稳态，就要有足够长的时间序列，通常在充分发展的 湍流中，需要105以上的时间积分步。 由于这些特殊要求，需要有内存大、速度快的计算机才能实现 湍流直接数值模 拟。
第6章 湍流直接数值模拟
6.2 湍流直接数值模拟的基本原理