动力前言课..docxVIP

湍流数值模拟方法简介摘要：湍流是自然界普遍存在的一种流体流动现象，对湍流的准确数值预测是当前基础学科中迫切需要解决的空气动力学前沿问题之一。目前湍流的数值模拟方法有：直接数值模拟（DNS），雷诺平均数值模拟（RANS）以及大涡模拟（LES）三种类型。本文分别介绍了这三种模拟方法的特点，并对这三种方法在计算量及精度等方面进行了对比。最后对湍流模拟方法的未来发展趋势做出了展望。关键词：直接数值模拟；雷诺平均数值模拟；大涡模拟；湍流模拟发展趋势0 序言湍流是自然界普遍存在的一种流体流动状态，湍流的研究是目前许多自然科学学科及工程部门中一个十分重要的基础问题。自从1883年Reynolds发现这一现象以来，一直是无数学者热门的研究课题。从物理结构来说，湍流是由大小不同的涡旋叠合而成的流动，涡旋的大小和旋转的方向分布都是随机的。大尺度涡旋主要由流动的边界所决定，其尺寸与流场的大小可比拟，是引起低频脉动的原因；小尺度涡旋主要是由于黏性力引起，其尺度小至流场千分之一的量级，造成高频脉动。大涡破裂后形成小涡，小涡破裂形成更小的涡旋，因而涡旋大小尺寸可在一个相当大的范围内连续变化。由于流体的黏性作用，小尺度涡旋不断消失,其机械能变成了热能；同时由于边界的作用、扰动及速度梯度等原因，不断有新的涡旋产生，因而构成了湍流运动。由于大小涡旋的随机运动造成了湍流物理量的随机脉动[1]。自从20世纪70年代以来，随着电子计算机的迅速发展，数值模拟成为研究湍流的有效方法。根据计算机的条件和研究湍流的不同目的，湍流数值模拟的精细程度有不同的层次。为了对湍流物理性质进行深入的了解，需要最精细的数值计算，这时，必须从完全精确的流动控制方程出发，对所有尺度的湍流运动进行数值模拟，这种最精细的数值模拟称为直接数值模拟，简写为DNS（Direetly NumericalSimulation）。实用上，只需要预测湍流的平均速度场、平均标量场和平均作用力时，可以从雷诺平均方程出发，在这一层次上的数值模拟称为雷诺平均数值模拟，简称RANS（ Reynolds-Averaged Navier-Stokes equations）。雷诺平均方程是不封闭的，必须引入雷诺应力的封闭模型才可解出平均。雷诺应力的主要贡献来自大尺度脉动，而大尺度脉动的性质和流动的边界条件密切相关，因此雷诺应力的封闭模式不可能是普适的，就是说，不存在对一切复杂流动都适用的统一封闭模式。介于DNS和RANS之间的数值模拟方法称为大涡数值模拟，简称LES（Large Eddy Simulation）。大涡数值模拟的思想是：大尺度脉动（或大尺度湍涡）用数值模拟方法计算，只将小尺度脉动对大尺度运动的作用做模型假设。LES的理论依据是小尺度脉动有局部平衡的性质，很可能存在某种局部普适的统计规律，如局部各向同性或局部相似性等。因而，小尺度脉动对大尺度运动的统计规律可能是普适的[2]。本文将分别探讨直接数值模拟（DNS）、雷诺平均数值模拟(RANS)以及大涡模拟(LES)各自的特点，并做湍流数值模拟的未来发展趋势做出一定的预测。1 直接数值模拟自20世纪70年代以来，湍流直接数值模拟取得显著的成果。Orzag和Patterson[3]于1972最早用直接数值模拟计算了各向同性湍流，当时网格数只有323,相应的雷诺数。虽然网格分辨率不够，计算结果不是很精确，但表明了谱方法可以用来对三维湍流大尺度结构的计算，对湍流高级数值模拟中最常用的数值方法--谱方法的发展作出了重大的贡献。Rogallo[4]发展了orzag-patterson算法，计算了各向同性湍流，验证了平均剪应力，无旋应力及其旋转对湍流产生的影响，并将结果同理论和实验数据进行对比，验证了湍流模型。从20世纪80年代初开始，Standfod大学的Reynolds、Ferziger和Moin[5]也在直接数值模拟方面做过大量的工作，如Moin等对弯曲槽道进行了模拟，1996年Moin[6]等又进行过对平板边界层分离的研究。1.1 湍流直接模拟的基本要求湍流流动是由许多大小不同尺度的涡旋组成，这些涡旋进行不规则的随机运动，为了能数值模拟各种大小不同尺度涡旋结构，对直接模拟有以下要求：1）空间分辨率高，网格尺寸小为了能准确描述小尺度的湍流特征，必须把研究体系至少分成相当于Kolmogorov微尺度大小网格，可认为与湍流动能秏散率和流体运动粘性有关，用量纲分析可得：（1.1）根据小尺度涡旋能量局部平衡假定，即认为大涡旋向小涡旋提供的能量与小涡旋本身消耗的能量达到平衡，可得，其中为脉动速度均方根值，为大旋涡尺度，它的最大尺寸可相当于流场宽度尺寸，因此一维网格数至少应为：（1.2）式中为流动雷诺数，。对于三维网格数至少应有：。如果为 104，则要求网格总数为109，对于目前计算机