计算流体力学概论.pptVIP

计算流体力学(CFD)概论 金晗辉 湍 流 大 涡 模 拟 简 介 * * 主 要 内 容 历史、发展 数学模型 计算方法 计算软件 1933年，英国科学家Thom用手摇计算机完成了对外掠圆柱流动的数值计算；此时，鉴于计算机本身的局限性，实验流体力学更具优势。（风洞、水槽等） 1965-1974 萌芽时期（交错网格、抛物型流动P-S方法、椭圆形流动SIMPLE等方法提出） 1975-1984 开始工业应用（GENMIX/ PHOENICS/ FLUENT/ STARCD/ CFX等软件开发面世） 1985－ 进一步发展（前、后处理软件/并行计算机/ DNSLES提出等） 计算流体力学的发展概况 历史、发展 数学模型 计算方法 计算软件 主 要 内 容 湍流数学模型的建立 基于N-S方程的湍流模式理论 基于分子运动学的湍流统计理论 直接模拟、湍流模式理论、大涡模拟、离散涡方法 格子器 N-S（Navier-Stokes)方程简介 连续性方程 动量方程 理论基础：连续介质假设、流体微团假设 浓淡燃烧器流场示意图 双方程模型 非线性 模型 多尺度 模型 RNG 模型 Reynolds应力模型（RSM） 代数应力模型（ASM） FLT模型 SSG模型 湍流模式理论以Reynolds时均运动方程和脉动运动方程为基础，依靠理论与经验的接合，引进一系列模型假设，从而建立一组描写湍流平均量的方程组。 湍流模式理论简介 平板尾迹流场示意图 气相流场全场涡结构 湍流模式理论局限性 对经验数据的依赖性; 将脉动运动的全部细节一律抹平从而丢失大量重要信息; 目前各种模型，都只能适用于解决一种或者几种特定的湍流运动。 湍流直接模拟（DNS）简介 计算机发展 数值算法发展 直接模拟 （DNS)技术的应用 出现大型并行计算机 Petaflops （1015）级 有限差分 有限元 谱方法 小波变换 自适应网格 并行计算技术 方程本身是精确的，不含任何认为假设和经验常数，仅有的误差只是由数值方法引入的误差 。 计算包括脉动运动在内的湍流所有瞬时流动量在三维流场中的时间演变; 不用任何湍流模型，直接数值求解完整的三维非定常的N-S方程组; 湍流旋涡结构包括大尺度涡和小尺度涡 大尺度涡 小尺度涡 湍 流 流 场 涡 结 构 图 直接模拟计算量太大，很难计算工程实际高雷诺数湍流流场。 为 什 么 要大 涡 模 拟？ 流 场 大尺度涡 小尺度涡 决定湍流流场的基本形态和性质； 流场质量、能量的主要携带者； 高度各向异性,无法建立统一模型。 由大涡非线性作用产生； 流场能量的主要耗散者； 近似各向同性,可以考虑建立统一模型。 小尺度涡对大涡的影响用模型进行模拟 大涡模拟思想 对大尺度涡进行直接模拟 离散涡方法 湍流模式理论所得到的流场变量都是时均值，人们无法从中获知流体脉动运动细节，而时变的离散涡方法则不存在这个问题。离散涡方法以大量的离散旋涡代替流动空间中连续分布的涡量场。 通过计算这些涡元的相互作用并计及涡元与固体的相互作用，根据毕奥－萨瓦定律（Biot-savart），最终得到大尺度涡结构的流场全貌。 离散涡方法优、缺点 优点： 对处理高Reynolds数，大尺度分离流有较好的效果。 与有限差分及有限元方法相比的一个重要特点是计算过程不需要网格，即Lagrangian描述而非Euler描述流体的运动过程，物理量并非在网格点之间进行传输，而是随涡元在流场中的位置变更而发生对流扩散，因此用于求解大Reynolds数流动时能避免网格方法存在的困难。 缺点： 尽管Lagrangian描述可以避免数值误差，但这就使得人们对平滑流动中产生的误差无法理解； 对粘性影响的处理也会带来较大的麻烦；以及数值稳定性、初生涡的放置等方面也都存在着缺陷。 湍流数学模型的建立 基于N-S方程的湍流模式理论 基于分子运动学的湍流统计理论 直接模拟、湍流模式理论、大涡模拟、离散涡方法 格子器 格子气与格子Boltzmann理论 理论基础： 一种观点认为，湍流的脉动运动无论在空间上还是在时间上都是很剧烈的随机运动，要想得到任何有意义的结果，都必须如同统计物理学中研究气体分子运动一样，应用统计方法。将统计方法与流体力学相结合即形成湍流统计理论。 特点：流体的离散性 将流体及其存在的时间空间完全离散，把流体想象成由大量只有质量没有体积的流体粒子组成，所有的这些粒子同步的随着离散的整时间步长，依据给定的简单碰撞规则在网格节点上相互作用（碰撞），沿网格线在节点间之间运动，碰撞规则遵循质量、动量和能量守恒定律。