哈工大高等流体力学介绍.docxVIP

高等流体力学研究报告考核科目:计算流体力学学生所在院（系）：能源科学与工程学院学生所在学科: 学生姓名: 学号: 湍流的数值模拟综述摘要：通过查找相关书籍和阅读大量相关文献，本文对湍流的形成条件以及在FLUENT中常用的湍流模型进行了分析总结，并对部分模型的机理和使用条件进行了详细说明，此外，本文还综述了当下发展较快的K-ε双方程模型和湍流燃烧大涡模拟(LES)的发展现状，浅析了其发展方向。关键词：湍流；FLUENT；数值模拟序言湍流是自然界广泛存在的流动现象，在多数工程问题中流体的流动往往处于湍流状态。如大气、海洋环境的流动，飞行器与船舰的绕流，叶轮机械、化学反应器、核反应器中的流体运动都是湍流。判断湍流的标准如下：对于外流而言，流动为湍流的判断标准一般为：外掠平板Rex500000。外掠障碍物Red20000。对于内流而言，Red2300；其中，，L=x、d、dh等。对于自然对流而言，，为瑞利数，Pr=为普朗特数。湍流流动的核心特征是其在物理上近乎于无穷多的尺度和数学上的强烈的非线性，这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。相比之下，计算机模拟要更具有可操作性，但是需要我们继续探求高精度的计算方法和实用可靠的网格生成技术，在研究湍流的同时，建立相应的模式，从而进行适当的模拟。Navier-Stokes方程为描述流体运行的控制方程。运用三维非稳态Navier-Stokes方程对流动进行直接数值计算时，必须采用很小的时间与空间步长才能分辨出湍流中详细的空间结构及变化剧烈的时间特征，而这对内存空间和CPU运行速度的要求非常高，除少数拥有超级计算机的研究者以外，Navier-Stokes方程的直接数值计算还无法运用于工程计算。因此，需要引入某种处理方法，对Navier-Stokes方程进行简化，使之可以适用于工程计算，这就是接下来即将介绍的大涡模拟方法和雷诺平均模拟方法。大涡模拟法。湍流的脉动与混合主要由大尺度的涡运动造成。大尺度的涡从主流中获取能量，通过相互作用，大尺度Navier-Stokes尺度的涡，在此过程中能量从大尺度涡传递到小尺度的涡，小尺度涡的主要作用是耗散能量，分析不同尺度涡的特点，发先大尺度的涡是各向异性的，且随流动的情形而异；而小尺度的涡几乎是各向同性的，且不同流动情形中小尺度涡有许多共性。根据对流动涡运动的认识，就产生了大涡模拟法。大涡模拟法的本质为用非稳态的Navier-Stokes方程来直接模拟大尺度涡运动，而通过模型来模拟小尺度涡运动，即大涡直接求解，小涡用模型。小涡对大涡的影响通过近似的模型来考虑，这称为亚格子雷诺应力( Subgrid Reynolds)。最早的亚格子模型为Smagorinsky模型，以后很多人对此进行了改进。大涡模拟对内存和CPU运行速度的要求虽然仍然很高，但远低于直接模拟方法对计算机资源的要求，在工作站和PC上都可以进行大涡模拟。雷诺平均模拟法。上述介绍的大涡模法对大尺度的涡仍然进行了直接求解，计算量仍然较大，而应用雷诺平均法，将Navier-Stokes方程对时间进行平均，引入模型进行求解，计算量更小，更适合应用于工程数值计算。雷诺平均模拟法需要引入模型的原因为：将Navier-Stokes方程对时间进行平均后，得出的关于时均物理量的控制方程中包含含了未知量，因此得到的方程的个数小于未知量的个数，无法直接求解。要使方程组封闭，使方程组可解，必须使方程组中方程的个数等于未知数的个数，这就要通过假设，创建模型来实现。1 FLUENT提供的湍流模型湍流模型如下图1所示，其中从零方程模拟到直接模拟，包含的机理越来越多，每次迭代计算量也增加。FLUENT提供的模型包括单方程模型和双方程模型。图1湍流模型框图下面将介绍基础的单方程模型、模型、模型、RSM模型和LES模型。1.1单方程模型在介绍一方程（单方程）模型前，首先给出湍流时均的连续性方程与雷诺方程： (1) (2) (3)单方程模型在湍流时均的连续性方程与雷诺方程的基础之上，再建立了一个湍动能k的输运方程，将湍动粘度μ表示成k的函数，从而使方程组封闭。湍动能k的输运方程为：(4)从左至右，方程(4)中各项依次为瞬态项、对流项、扩散项、产生项、耗散项。由普朗特表达式有： (5)式中，、CD、Cμ为经验常数，=1，Cμ=0.09，CD=0.08~0.38。一方程模型中如何确定长度尺寸l仍然是不易解决的问题，因此很难得到推广应用。1.2 K-ε模型模型为双方程模型，是在一方程模型的基础之上，再引入了一个关于湍动耗散率的方程，该模型是目前使用最为广泛的湍流模型，又分为标准模型、重组化群模型和可实