流动湍流及其控制的数值模拟研究.pptVIP

----------------------------------------------------------- 流动湍流及其控制的数值模拟 唐国宁 2005．11．30 内容提要 流动湍流研究的简要回顾 湍流控制的一般方法 直接数值模拟的常用方法 流动湍流控制的数值模拟结果 (一)流动湍流研究的简要回顾 层流 流体质点的运动轨迹光滑而有规则 湍流 是一种有结构、不规则、多尺度的随机流动 最早由雷诺圆管实验, 人们认识到湍流的随机性 1894年,O.Reynolds 首次将湍流流动分解为平均量与脉动量之和 ，并由此导出湍流的统计方程 ． 雷诺平均方程 脉动方程 封闭湍流统计方程的方法叫湍流模式 代数涡黏模式、 涡黏模式（包括线性和非线性），代数应力模式和雷诺应力模式（或称二价矩模式） 20世纪初，一些杰出的流体力学家相继对涡黏系数提出各种流体力学模型，如Taylor(1921年)的涡模型，Prandtl(1925年)的混合长度模型和Vonkaman(1930年)相似模型等. 考虑连续介质不规则运动的特点－－运动的多尺度性。导致产生描述多尺度的谱概念和谱方法，并最终产生了Kolmogorov的局部各向同性的通用谱（即-5/3谱） 近50年来，大量的物理实验和数值模拟使得人们对湍流本质的认识越来越深入，人们注意到湍流并不是完全随机的，而是有某种相干结构，湍流是一种有结构的不规则多尺度流动． 实验和数值模拟表明：无论是简单还是复杂湍流，都存在一定的大涡结构。湍流中大涡拟序结构对湍流生成和发展有主宰作用，因此抑制或消除大涡结构可能抑制整体的湍流强度，甚至使流动层流化，这是近代湍流和降噪的主要思想 湍流的有结构、不规则、多尺度流动导致了研究湍流的困难． 最好的湍流模型就是N-S方程 湍流的直接数值计算也存在局限性。为了能模拟最小尺度的涡，计算网格的分辨率要足够高,而计算区域的尺寸应足以容纳最大的涡，因此三维计算域的网格点数具有 　 ，例如当Re=5000时,三维计算域的网格点数至少有 2.1* 10**8, 由于计算机内存的限制，目前直接数值模拟只能研究低雷诺数湍流Re5000. 从20世纪以来，科学家虽然对湍流的不规则性进行了深入的理论探讨，以Kolmogorov为代表的莫斯科学派和以Taylor、batchelor为代表的剑桥学派对湍流有杰出的贡献，但是认真考察一百多年的湍流文献，人们发现推动湍流研究的理论屈指可数。可以称做奠基性的理论只有Kolmogorov的局部各向同性湍流理论。 （二）湍流控制的一般方法 　在流体力学领域，对流体流动的控制就是要使壁湍流和无剪切流动发生有益的变化，例如使转捩延迟或提前；抑制或增强湍流；防止或促进分离等。产生有用的结果是：使阻力减少、升力增加、混合增强和抑制流动引起的噪声。迄今人们提出了许多湍流控制方法，有主动控制和被动控制两大类，主要控制方法有： 壁面的抽吸或注入。 添加聚合物、表面话化剂、纤维、微型水泡。 用磁场控制流场。 电磁力控制流场。 化学反应改变流场温度或流体密度。 改变流场几何参数．如在壁面开小尺度的纵向糟或肋条，使壁面振荡等。 多孔壁。 在壁面中嵌入物体。 （三）直接数值模拟常用方法 湍流的数值模拟方法有两种： 大涡模拟 直接数值模拟 　. 大涡模拟的基本思想是把包括脉动在内的湍流的瞬时运动通过某种滤波方法分解为大尺度运动和小尺度运动两部分。 大尺度运动（大涡）强烈依赖于初始条件，其形态与强度因平均流动的不同而不同。 大涡与平均运动之间有强烈的相互作用，它直接由平均运动提供能量，反过来又对平均运动有强烈的影响，大部分的质量、动量和能量的输运是由大涡运动引起。 大尺度量要通过数值求解运动微分方程计算出来． 根据经典的湍流统计理论，认为小尺度运动（小涡）是完全随机的，与平均运动和流场的边界形状几乎没有关系，因此对小尺度运动用建立模型进行模拟。 小尺度对大尺度运动的影响在运动方程中表现为类似于雷诺应力项，称为亚格子雷诺应力. 湍流的直接数值模拟就是数值求解完整的、三维非定常的Navier-Stokes方程组，计算包括脉动运动在内的湍流所有瞬时流动量在三维流场中的时间演化． 最早的直接数值模拟工作开始于上个世纪70年代。 到目前为止，在现有的计算能力限制下，直接数值模拟还只能计算中等雷诺数以下且具有简单几何边界条件的湍流流动，网格可达到512＊512＊512。在这种情况下，数值模拟完全可以代替物理实验。 对于湍流的直接数值模拟，数值方法主要有有限元法、有限体积法、高阶有限差分法和谱方法等。 对于直槽道湍流的直接模拟，由于其边界简单，通常采用谱方法 谱方法具有以下优点：1）精度高．2）具有准确的