计算流体力学绪论

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