湍流的模拟分析.pptVIP

第一部分 湍流认识及N-S 方程 湍流的认识 从物理机理上说: 可以把湍流看成是各种不同尺度的涡旋叠合而成的，这些涡的大小及旋转方向分布是随机的。 大尺度的涡旋主要由流动的边界条件所决定，其尺寸可与流场的大小相比拟，是引起低频脉动的原因; 小尺度的涡主要是由粘性力所决定的，其尺寸可能只有流场尺度的千分之一，是引发高频脉动的原因。 湍流运动的构成: 大尺度的涡破裂后形成小尺度的涡，较小尺度的涡破裂后形成更小的涡。大尺度的涡从主流获得能量，通过涡间的转化将能量传给小尺度的涡。最后由于粘性的作用，小尺度的涡不断消失，机械能就转化 即耗散 为流体的热能。同是，由于边界的作用，扰动及速度梯度的影响，新的涡又不断产生。 湍流一个重要特点: 物理量脉动，非稳态N-S方程对湍流运动仍是适用的。 第二部分 传统湍流模型 传统模式理论 Reynolds平均法 在这类方程中，将非稳态N-S方程对时间作平均，即把湍流运动看成二个流动的叠加: 时间平均流动 瞬时脉动流动。 所得的时均的N-S方程中包含了脉动量乘积的时均值等未知量，称为Reynolds应力，它包括了六个未知量。显然方程的个数小于未知量的个数。要让方程封闭，必须作出假设。 大涡模拟 LES 按湍流的机理，湍流的脉动及混合主要由大尺度的涡造成。大尺度的涡高度非线性，其相互作用把能量传给小尺度的涡。小尺度的涡几乎是各向同性的，它们起到能量耗散的作用。 导致了只用非稳态N-S方程模拟大涡，不直接计算小涡，而将小涡对大涡的影响通过近似的模型来考虑，这种影响称为亚格子Reynolds应力。 大涡模拟对计算机的要求仍比较高，但比直接模拟低得多。 直接模拟 DNS 用三维非稳态N-S方程对湍流进行直接数值计算。要对高度复杂的湍流进行直接的计算，必须采用很小的时间和空间步长，才能分析出详细的空间结构及变化剧烈的时间特性。故计算量很大，对计算机的要求很高。 模型评价 零方程-- Boussinesq涡粘模型 湍流粘性系数法 基于Boussinesq1887年的假设，它将湍流脉动所造成的附加应力 Reynolds应力 同层流运动应力那样与时均的应变率关联起来. ?T为湍流粘性系数,是标量且为常数; 这一假设并无物理基础，且采用各向同性的湍流动力粘度来计算湍流应力，难于考虑旋转流动和表面曲率变化的影响，但以此为基础的湍流模型目前在工程计算却应用最为广泛。 零方程– Prandtal混合长度理论 单方程模型 单方程模型缺陷 双方程模型 双方程模型 近壁区的处理 以上模型均是针对湍流发展非常充分的湍流流动来建立的，是针对高Re数的湍流计算模型，适用于离开壁面一定距离的湍流区域。这里的Re数是以湍流脉动动能的平方根作为速度 又称湍流Re数 计算的，是分子扩散造成的动力粘性。 而在Re数比较低的区域，湍流发展不充分，湍流的脉动影响可能不如分子粘性大，在贴近壁面的底层内，流动可能处于层流状态。 必须采用特殊的处理，一般有二种解决方法， 1 壁面函数法 2 低Re数的k- ?模型。 壁面函数法的处理 壁面函数法的基本思想是： 对于湍流核心区的流动使用k-?模型求解; 而在壁面区不进行求解，直接使用半经验公式将壁面上的物理量与湍流核心区内的求解变量联系起来。它需要把第一个节点布置在对数律层，对第一个节点的值由公式确定。 不需要对壁面内的流动进行求解，可直接得到与壁面相邻控制体积的节点变量。各种改进的壁面函数法越来越准确的模拟壁面的相关特性。 低Re数的k-?模型的处理 为了使k-?模型能够计算到壁面，出现了各种低Re数的k-?模型，它实际上是对k-?模型的控制方程的各个项作出相应的修改，以体现壁面附近流动的各种真实特征。 如: 控制方程中的扩散项包括湍流扩散和分子扩散二部分; k的控制方程中壁面脉动动能的耗散是各向异性。 Reynolds应力模型 RSM和ASM 上述方程湍流模型都假定湍流粘性系数是各向同性的;采用了湍流粘性的假设,用有效粘性系数和平均速度梯度的乘积来模拟雷诺应力。这些模型难于反映旋转流动及流动方向表面曲率变化的影响，有必要对湍流脉动应力直接建立微分方程求解。 在应力方程模型中，对二个脉动值乘积的时均值方程直接求解，而对三个脉动值乘积的时均值，采用模拟方式计算，不少研究者认为这是目前最有发展前途的湍流模型，这就是Reynolds应力方程模型 RSM 。 为了减轻RSM的计算工作量，将Reynolds应力用代数方程式而不是用微分方程来求解，用代数方程去近似的模拟微分方程，这就是代数应力方程模型 ASM 。 雷诺应力模型 湍 流 大 涡 模 拟 简 介 二维大涡模拟亚网格模型 第三部分 湍流直接数值模拟 数值模拟大致步骤： （1）建立基本守衡方程 由流体力学、传热学、燃烧学热等离子体