计算流体力学绪论.ppt

课程安排1.教学内容第一讲 绪论第二讲 计算区域及控制方程的离散化第三讲 差分格式第四讲 原始变量法第五讲 湍流模型第六讲 Gambit的使用第七讲 Fluent的使用2.教学课时选修课：40课时 包括：讲解和实践课程安排3.参考书n自带的帮助文件n计算流体动力学分析 王福军 清华大学出版社nFluent流体工程仿真实例与应用 韩占忠等 北京理工大学出版社 1-1 1-1 概述概述n什么是计算流体力学？借助计算机，在流动基本方程控制下对流动的数值模拟。nCFD的基本思想：把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。1-1 1-1 概述概述 CFD可以看做是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。1-2 1-2 计算流体力学在流体力学中的地位计算流体力学在流体力学中的地位n实验流体力学n理论流体力学：连续场模型、分析解法n计算流体力学：离散场模型、数值解法 单纯实验测试 单纯理论分析 计算流体力学1-2 1-2 计算流体力学在流体力学中的地位计算流体力学在流体力学中的地位 实验测量方法实验测量方法:所得到的实验结果真实可信，它是理论分析和数值方法的基础。局限性：局限性：(1)实验往往受到模型尺寸、流场扰动、人身安全和测量精度的限制，有时可能很难通过试验方法得到结果。(2)实验还会遇到经费投入、人力和物力的巨大耗费及周期长等许多困难。Important!1-2 1-2 计算流体力学在流体力学中的地位计算流体力学在流体力学中的地位 理论分析方法理论分析方法:所得结果具有普遍性，各种影响因素清晰可见，是指导实验研究和验证新的数值计算方法的理论基础。局限性：局限性：它往往要求对计算对象进行抽象和简化，才有可能得出理论解。对于非线性情况，只有少数流动才能给出解析结果。CFD方法克服了前面两种方法的弱点，在计算机上实现个特定的计算，就好像在计算机上做一次物理实验。计算流体动力学的特点计算流体动力学的特点p流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解p可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较p它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。p数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差。p它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证。计算流体动力学的特点计算流体动力学的特点p程序的编制及资料的收集、整理与正确利用，在很大程度上依赖于经验与技巧。p因数值处理方法等原因有可能导致计算结果的不真实，例如产生数值粘性和频散等伪物理效应。pCFD因涉及大量数值计算，因此，常需要较高的计算机软硬件配置。计算流体动力学的特点计算流体动力学的特点n理论分析成本最低n结果最理想n影响因素表达清楚n缺点：局限与非常简单的问题,n数值方法n成本较低：数值实验n适用范围宽n缺点：可靠性差，表达困难n实验测量n可靠n成本高 将三种方法有机结合，互为补充，必然会取得相得益彰的效果各种方法比较各种方法比较n给出物理模型（给出物理模型（Physical model/description)Physical model/description)n借助基本原理借助基本原理/定律给出定律给出数学模型数学模型数学模型数学模型（Mathematical Mathematical model)model)n质量守恒（Mass Conservation）n能量守恒（Energy Conservation）n动量守恒（Momentum Conservation）n傅立叶定律（Fouriers heat conduction law)n菲克定律（Ficks mass diffusion law）n牛顿内摩擦定律（Newtons friction law）n。CFDCFD：总体步骤：总体步骤 比如，我们研究管道内的流体流动，抽象出来一个直管，和粘性流体模型，或者我们认为管道内的液体是没有粘性的，使用一个直管和无粘流体模型.还有，我们根据热传导定律，认为固体的热流率是温度梯度的线形函数，相应的傅立叶定律就是导热问题的物理模型。因此，不难理解物理模型是对实际问题的抽象概念，对实际问题的一种描述方式，这种抽象包括了实际问题的几何模型，时间尺度，以及相应的物理规律。把实际的问题，通过相关的物理定律概括和抽象出来并满足实际情况的物理表征。物理模型 对物理模型的数学描写。数学模型 比如N-S方程就是对粘性流体动力学的一种数学描写，值得注意的是，数学模型对物理模型的描写也要通过抽象，简化的过程。1-3 流动与传热问题的数学描述p控制方程式（不可压缩流体）连续性方程：动量方程：能量方程：通用微分方程：非稳态项 对流项 扩散项 源项 广义变量 广义扩散系数1 u v T w 连续性方程 N-S方程 能量方程 不同方程中 代表的物理量p 单值性条件l初始条件：l边界条件：（1）第一类边界条件：边界上 给定（2）第二类边界条件：边界上 给定（3）第三类边界条件：边界上 与 关系给定给定 边界条件的通用形式B=0 第一类边界条件A=0 第二类边界条件 第三类边界条件对于初始条件和边界条件的处理，直接影响计算结果的精度 举例后台阶突扩流动物理模型：后台阶突扩层流流动与换热问题（二维、稳态、不可压缩、常物性、无内热源）连续方程：动量方程：能量方程：举例后台阶突扩流动数学模型：举例后台阶突扩流动(1)进口 (速度u、温度T)分布要给定(2)中心线(3)壁面(4)出口 数值处理方法边界条件:初始条件与边界条件是控制方程有确定解的前提，初始条件与边界条件是控制方程有确定解的前提，控制方程与相应的初始条件、边界条件的组合控制方程与相应的初始条件、边界条件的组合构成对一个物理过程完整的数学描述。构成对一个物理过程完整的数学描述。1-4 流动与传热问题数值求解的基本思路n初始条件：n基于连续介质模型(1)用有限个离散点上的值代表连续体 区域离散(2)建立离散点上物理量的代数方程式 方程的离散(3)求解代数方程物理问题数值求解的基本过程物理问题数值求解的基本过程1-5 常用数值方法 经过四十多年的发展，CFD出现了多种数值解法。这些方法之间的主要区别在于对控制方程的离散方式。根据离散的原理不同，CFD大体上可分为三个分支：p有限差分法(Finite Different Method，FDM)p有限元法(Finite EIement Method，FEM)p有限体积法(Finite Volume Method，FVM)有限差分法是应用最早、最经典的CFD方法，它将求求解域划分为差分网格，用有限个网格节点代替连续的求解解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的导数用差商代替，推导出含有离域，然后将偏微分方程的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。散点上有限个未知数的差分方程组。求出差分万程组的解，就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题的近似数值解法。这种方法发展较早，比较成熟，较多地用于求解双曲型和抛物型问题。在此基础上发展起来的方法有PIC(Particle-in-cell)法、MAC(Marker-and-Cell)法，以及南美籍华人学者陈景广提出的有限分析法(Finite Analytic Method)等.有限差分法有限元法 有限元法是20世纪80年代开始应用的种数值解法，它吸收了有限差分法中离散处理的内核吸收了有限差分法中离散处理的内核，又采用了变分计算中选择逼近函数对区域进行积分变分计算中选择逼近函数对区域进行积分的合理方法。有限元法因求解速度较有限差分法和有限体积法慢，因此应用不是特别广泛。在有限元法的基础上，英国C ABBrebbia等提出了边界元法边界元法和混合元法混合元法等方法。有限体积法是将计算区域划分为一系列控制体积，将待解微分方程对每一个控制体积积分得出离散方程。有限体积法的关键关键是在导出离散方程过程中，需要对界面上的对界面上的被求函数本身及其导数的分布作出某种形式的假定被求函数本身及其导数的分布作出某种形式的假定。用有限体积法导出的离散方程可以保证具有守恒特性，而且离散方程系数物理意义明确，计算量相对较小。有限体积法 1980年，S.V.Patanker在其专著Numericacl Heat Transfer and Fluid Flow中对有限体积法作了全面的阐述。此后，该方法得到了广泛应用，是目前CFD应用最广的一种方法。当然，对这种方法的研究和扩展也在不断进行，如PChow提出了适用于任意多边形非结构网格的扩展有限体积法。有限体积法离散方法分类小结n有限差分法（Finite difference method）n用差商与代替导数n经典、成熟n数学理论基础明确n主导方法n有限容积法(Finite volume method)n控制容积法（Control volume method）n基本上属于有限差分法的范畴n有限元法（Finite element method)n将求解区域分成若干个小的单（element）n设定待求变量在单元上的分布函数n适应性强，适用于复杂的求解区域n一度有取代有限差分法的趋势n程序技巧要求高n数学基础不如有限差分法明确Particle-in-cell PICMarker-and-Cell MAC有限分析法 Finite Analytic Method边界元法和混合元法 BEM数值积分变换法 NITM 格子Boltzmann方法 LBM直接模拟Monte Carlo方法分子动力学模拟还有其它一些方法：1-5 流动与传热问题控制方程的基本类型及其对数值解的影响n从数学角度从数学角度0 双曲型方程椭圆型问题离散方程：联立求解（整场求解）抛物型问题离散方程：步进求解。单向坐标单向坐标扰动沿该坐标轴仅向一个方向传递，或该坐标轴上任一点处的物理量仅受其一侧条件的影响（如非稳态：t，边界层：x）。双向坐标双向坐标扰动可以向两个方向传递，坐标轴上一点物理量之值受其两侧条件的影响（非边界层：x，y；边界层：y）。定义定义：如果微分方程在任意有限大区域V内积分所得表达式都表示该区域的守恒定律，该方程为守恒型方程，否则为非守恒型方程。(1)对流项表示成散度形式。(2)对有限大小的区域只有守恒型方程才能保证物理守恒定律成立。(3)要使离散方程具有守恒特性，必须从守恒型方程出发。n从物理角度从物理角度守恒型与非守恒型有限容积法中将所计算的区域划分成有限容积法中将所计算的区域划分成一系列控制容积，每个控制容积都有一系列控制容积，每个控制容积都有一个节点作代表。通过将守恒型的控制一个节点作代表。通过将守恒型的控制方程对控制容积做积分来导出离散方程。方程对控制容积做积分来导出离散方程。问题：1）计算流体力学的基本任务是什么？3）CFD的基本工作思路是什么？2)什么叫控制方程？给出不可压流体对流换热问题的控制方程。Fluent工程实例汽车空调Fluent工程实例Fluent工程实例Fluent工程实例呼吸过程