计算流体动力学基础.pptx
流体动力学仿真,电子科技大学机械电子工程学院,第二部分：计算流体动力学,1.教学内容计算流体动力学基础知识计算区域及控制方程的离散化代数方程组的求解算法流场仿真的经典算法SIMPLE算法边界条件处理常用软件介绍,第二部分：计算流体动力学,2.教学目标初步掌握计算流体动力学的基本原理和基本方法。,3.参考书目计算流体动力学分析，王福军有限体积法基础，李人宪数值传热学，陶文铨,第七讲 计算流体动力学基础,Lecture 7 Fundamental of CFD,本讲内容,7-1 流体力学知识回顾7-2 计算流体动力学概述What?计算流体动力学是什么Why?为什么研究计算流体动力学How?计算流体动力学研究的一般方法7-3 计算流体动力学的发展与应用,7-1 流体力学知识回顾,一、基本概念,理想流体与粘性流体牛顿流体与非牛顿流体可压流体与不可压流体定常与非定常流动层流与湍流对流与扩散,7-1 流体动力学知识回顾,二、基本方程,连续性方程,连续性方程是质量守恒定律在流体运动中的体现，是一切流体运动必须遵循的普遍原则。,或用另一种散度符号写为另一种形式：,7-1 流体动力学知识回顾,二、基本方程,连续性方程,直角坐标系下，连续方程式可写为：,三元流动：,二元流动：,一元流动：,7-1 流体动力学知识回顾,二、基本方程,连续性方程,若为定常流动，连续方程式可简化为：,三元流动：,二元流动：,一元流动：,或：,7-1 流体动力学知识回顾,二、基本方程,连续性方程,若为不可压流动，连续方程式可简化为：,三元流动：,二元流动：,一元流动：,7-1 流体动力学知识回顾,二、基本方程,N-S方程,或：,N-S方程是动量守恒定律在流体运动中的体现，是粘性流体遵循的原则；当不考虑粘性时，方程退化为欧拉方程,7-1 流体动力学知识回顾,二、基本方程,N-S方程,若为定常流动，方程式可简化为：,若为不可压流动，方程式可简化为：,7-1 流体动力学知识回顾,二、基本方程,传热方程传热方程是能量守恒定律在流体运动中的体现传质方程传质方程是物质的量守恒在流体运动中的体现,7-1 流体动力学知识回顾,二、基本方程,方程的通用形式,以上流体动力学方程可写为统一形式：,通用变量，可代表1,u,v,w,T等求解变量 广义扩散系数S广义源项,基于通用形式的方程可以实现以同样的程序求解不同的变量，区别仅在于与S的具体表达式。,7-1 流体动力学知识回顾,二、基本方程,方程的通用形式,以上流体动力学方程可写为统一形式：,(2),(3),(4),(1),(1)：瞬态项流场的非定常性对物理量的影响(2)：对流项流体的宏观运动对物理量的影响(3)：扩散项流体中微观粒子运动对物理量的影响(4)：源项其他所有因素(外力等)对物理量的影响,7-2 计算流体动力学概述,一、What什么是计算流体动力学？,Computational Fluid Dynamics(CFD)：是应用计算机和离散化的数值方法对流体动力学问题进行数值模拟和分析的学科；借助计算机，在流动基本方程控制下对流动的数值模拟；是流体动力学仿真的基本原理和实现手段。,7-2 计算流体动力学概述,一、What什么是计算流体动力学？,时间域及空间域上连续物理量的场,有限个离散点上的变量值的集合,描述离散点上场变量之间关系的代数方程组,场变量的近似值,时域及空间域离散化,流动规则离散化,解方程,基本思想,7-2 计算流体动力学概述,一、What什么是计算流体动力学？,关键：离散化,7-2 计算流体动力学概述,二、Why为什么研究计算流体动力学？,实验往往受到模型尺寸、流场扰动、人身安全和测量精度的限制，有时可能很难通过试验方法得到结果。实验还会遇到经费投入、人力和物力的巨大耗费及周期长等许多困难。,理论或实验研究方法在实验应用中有诸多限制：,流动问题一般是非线性的，理论研究往往需要对计算对象进行抽象和简化，才有可能得出理论解。到目前为止，只有极少数流动能给出解析结果。,7-2 计算流体动力学概述,二、Why为什么研究计算流体动力学？,通过计算流体动力学的数值模拟，可以实现：得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布。,7-2 计算流体动力学概述,7-2 计算流体动力学概述,Velocity magnitude(0-6 m/s)on a dinosaur,7-2 计算流体动力学概述,Velocity vectors around a dinosaur,7-2 计算流体动力学概述,二、Why为什么研究计算流体动力学？,通过计算流体动力学的数值模拟，可以实现：得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布。得到物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量，如旋转式流体机械的转矩、水力损失和效率等。根据计算结果，进行结构优化设计等。,计算流体动力学研究特点：,CFD方法并不追求获得非线性流动问题的解析解，而是应用数值方法找出满足工程需要的近似解它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较,7-2 计算流体动力学概述,二、Why为什么研究计算流体动力学？,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,Step 1.给出物理模型（Physical model/description)Step 2.借助基本原理给出数学模型（Mathematical model),总体步骤：,物理模型是指把实际的问题，通过相关的物理定律概括和抽象出来并满足实际情况的物理表征。,数学模型是指对物理模型的数学描写，是CFD方法的基础和出发点。,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,完整的数学模型应不仅包括控制方程，还应包括初始条件和边界条件。,初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况。对于瞬态问题，必须给定初始条件。对于稳态问题，不需要初始条件。边界条件是在求解区域的边界上所求解的变量或其导数随地点和时间的变化规律。对于任何问题，都需要给定边界条件。,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,物理模型：后台阶突扩层流流动与换热问题（二维、稳态、不可压缩、常物性、无内热源）,例：后台阶突扩流动,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,连续方程：,动量方程：,能量方程：,数学模型：,例：后台阶突扩流动,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,(1)进口(速度u、温度T)分布要给定,(2)中心线,(3)壁面,(4)出口 数值处理方法,边界条件:,例：后台阶突扩流动,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,Step 3.离散化,经过四十多年的发展，CFD出现了多种数值解法。这些方法之间的主要区别在于对控制方程的离散方式。根据离散的原理不同，CFD大体上可分为三个分支：,有限差分法(Finite Different Method,FDM)有限元法(Finite EIement Method,FEM)有限体积法(Finite Volume Method,FVM),7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,有限差分法(FDM)有限差分法是应用最早、最经典的CFD方法有限差分法是一种直接将微分问题变为代数问题的近似数值解法,将求解域划分为差分网格，用有限个网格节点代替连续的求解域将偏微分方程的导数用差商代替特点：经典、成熟，数学理论基础明确,原理：,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,有限元法(FEM)有限元法在CFD中应用较晚，但具有较大发展潜力,将求解区域分成若干个小的单元（element）设定待求变量在单元上的分布函数采用了变分计算中选择逼近函数对区域进行积分,原理：,适应性强，适用于复杂的求解区域一度有取代有限差分法的趋势程序技巧要求高数学基础不如有限差分法明确,特点：,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,有限体积法(FVM)有限体积法是目前计算流体动力学中应用最广的一种方法,原理：将计算区域划分为一系列控制体积，将待解微分方程对每一个控制体积积分得出离散方程。,关键：在导出离散方程过程中，需要对界面上的被求函数本身及其导数的分布作出某种形式的假定。,特点：用有限体积法导出的离散方程可以保证具有守恒特性，而且离散方程系数物理意义明确，计算量相对较小。,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,其他方法,Particle-in-cell PICMarker-and-Cell MAC有限分析法 Finite Analytic Method边界元法和混合元法 BEM数值积分变换法 NITM 格子Boltzmann方法 LBM直接模拟 Monte Carlo方法分子动力学模拟,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,Step 4.求解,无论离散方法为FDM、FEM或FVM，处理的结果总是和单元或节点数量相对应的代数方程组方程组的求解方法包括直接法（如Cramer法则、高斯消元法）以及迭代法（如Jacob迭代）直接法只适用于网格数量很少的情况，一般较少应用，而实用一般为考虑方程组系数的带状稀疏性的特定迭代法（如TDMA)商业软件中求解过程一般在后台进行,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,Step 5.结果显示,线值图：在二维或三维空间上，将横坐标取为空间长度或时间历程，将纵坐标取为某一物理量，然后用光滑曲线或曲面在坐标系内绘制出某一物理量沿空间或时间的变化情况。,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,矢量图：直接给出二维或三维空间里矢量（如速度）的方向及大小，一般用不同颜色和长度的箭头表示速度矢量。矢量图能形象地显示流动特征。,Step 5.结果显示,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,等值线图(contour)：同一线上的物理量(如压强)相等。,振荡机翼附近的压力等值分布,Step 5.结果显示,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,云图：使用渲染的方式，将流场某个截面上的物理量(如压力或温度)用连续变化的颜色块表示其分布。,振荡机翼附近的压力云图,Step 5.结果显示,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,流线,Step 5.结果显示,7-2 计算流体动力学概述,三、HOW计算流体动力学研究的一般方法,建立控制方程,确立初始条件及边界条件,划分计算网格，生成计算节点,建立离散方程,给定求解参数，进行求解,解收敛否,显示和输出计算结果,否,时域及空间域离散化,控制方程离散化,初始和边界条件离散化,离散初始条件和边界条件,数学模型,离散化,物理模型,是,7-3 计算流体动力学的发展与应用,上世纪50年代，奠定了一定理论基础1965年Harlow（美国）发表“流体力学的计算机实验”，模拟卡门涡街。Macagno（法国）发表“水力学模拟的某些新方向”，计算突然扩大管流。20世纪70-80年代计算机技术的发展为计算流体力学的应运而生提供了条件。近代一大批高新技术产业对计算流体力学提出了新的要求，同时也为计算流体力学的发展提供了新的机遇。实践与理论的不断互动，形成计算流体力学的新热点、新动力，从而推动计算流体力学不断向前发展。,7-3 计算流体动力学的发展与应用,1970 年代，飞机设计主要依赖风洞实验 YF-17研制，风洞实验13,500小时1980年代，CFD逐渐发展，部分取代实验 YF-23，风洞实验5,500小时，CFD计算15,000小时,7-3 计算流体动力学的发展与应用,90年代，CFD 在飞机设计中发挥了主力作用 波音777，CFD占主角 2000 之后，CFD 取代了大部分风洞实验 波音787：全机风洞实验仅3次,7-3 计算流体动力学的发展与应用,Where is CFD used?AerospaceAutomotiveBiomedicalChemical ProcessingHVACHydraulicsMarineOil&GasPower GenerationSports,F18 Store Separation,Temperature and natural convection currents in the eye following laser heating.,Aerospace,Automotive,Biomedical,7-3 计算流体动力学的发展与应用,Polymerization reactor vessel-prediction of flow separation and residence time effects.,HVAC,Chemical Processing,Hydraulics,Where is CFD used?AerospaceAutomotiveBiomedicalChemical ProcessingHVACHydraulicsMarineOil&GasPower GenerationSports,Streamlines for workstation ventilation,7-3 计算流体动力学的发展与应用,Flow of lubricating mud over drill bit,Flow around cooling towers,Oil&Gas,Power Generation,Where is CFD used?AerospaceAutomotiveBiomedicalChemical ProcessingHVACHydraulicsMarineSportsOil&GasPower Generation,Marine,Sports,7-2 计算流体动力学概述,算例1：阀门内部流动数值模拟,开度10,7-2 计算流体动力学概述,算例1：阀门内部流动数值模拟(速度分布),开度50,开度100,开度10,7-2 计算流体动力学概述,开度50,开度100,算例1：阀门内部流动数值模拟(压力分布),开度10,开度50,开度100,7-2 计算流体动力学概述,算例1：阀门内部流动数值模拟(压力分布),7-2 计算流体动力学概述,算例2：导弹飞行器(1),大攻角下的分离流线,算例2：导弹飞行器(2),7-2 计算流体动力学概述,算例2：导弹飞行器(3),7-2 计算流体动力学概述,算例3：质子电解质膜(PEM)燃料电池,7-2 计算流体动力学概述,算列4：离心叶轮（1）,7-2 计算流体动力学概述,叶轮通道内的压力分布,93叶高处相对速度矢量线,7-2 计算流体动力学概述,
