CFD2015-第1讲-基本方程.pptx

2015年春季中国科学院大学课程计算流体力学 李新亮mailto:lixl@lnm.imech.ac.cnlixl@imech.ac.cnTel:学所主楼214参考数目： 傅德薰等：《计算空气动力学》 任玉新等： 《计算流体力学基础》；阎超：《计算流体力学方法及应用》， J. Blazek: Computational Fluid Dynamics: Principles and Applications E. F. Toro: Riemann Solvers and numerical methods for fluid dynamics 第一讲 流体力学基本方程计算流体力学(CFD) 的概念及意义 流体力学的基本方程 偏微分方程组的类型 重点: 流体力学基本概念：连续介质假设，流动描述方法 N-S方程及其无量纲化（熟记）； 双曲型方程性质；Copyright by Li Xinliang§ 1.1 绪论计算流体力学: Computational Fluid Dynamics简称CFD计算流体力学是通过数值方法求解流体力学控制方程，得到流场的离散的定量描述，并以此预测流体运动规律的学科CFD: 通过离散求解流动方程得到流动信息Copyright by Li Xinliang 理论解（解析解）精确解： Poiseuille解， Blasius解， Plantdl 湍流边界层解流动控制方程渐进解、近似解： Stokes解数值解方程复杂（非线性偏微方程组）， 解析解很难获得差分法、 有限体积法、边界元法、谱（元）方法、 粒子方法 ……借助计算机来实现数值求解在计算机产生之前，数值方法已然产生 Copyright by Li Xinliang计算流体力学(CFD): 在航空航天领域得到广泛应用● 1970 年代， 飞机设计主要依赖风洞实验 YF-17研制，风洞实验13,500小时● 1980年代，CFD逐渐发展， 部分取代实验 YF-23，风洞实验5,500小时，CFD计算15,000机时YF17YF23Copyright by Li XinliangYF17● 90年代， CFD 在飞机设计中发挥了主力作用 波音777， CFD占主角 ●2000 之后， CFD 取代了大部分风洞实验 波音787：全机风洞实验仅3次● 航天领域，CFD发挥着实验无法取代的作用 实验难点：复现高空高速流动条件 波音787波音777Copyright by Li XinliangCFD 面临的挑战及主要任务：多尺度复杂流动的数学模型化; 湍流的计算模型； 转捩的预测模型； 燃烧及化学反应模型； 噪声模型……可处理间断及多尺度流场的高分辨率、强鲁棒性、高效数值方法； 高精度激波捕捉法； 间断有限元法; …… 可处理复杂外形、易用性强的算法； 复杂外形—— 网格生成工作量大 多块分区算法； 无网格法； 粒子算法；Copyright by Li Xinliang传统计算方法： 有限差分法， 有限体积法 ， 有限元法， 谱方法（谱元法）等；最近发展的方法:基于粒子的算法（格子-Boltzmann, BGK），无网格 优点缺点适用范围有限差分法简单成熟，可构造高精度格式处理复杂网格不够灵活相对简单外形的高精度计算有限体积法守恒性好，可处理复杂网格不易提高精度（二阶以上方法复杂）复杂外形的工程计算有限元法基于变分原理，守恒性好对于复杂方程处理困难多用于固体力学等间断有限元法(DG)精度高、守恒性好、易于处理复杂网格计算量大；捕捉激波（限制器）难度大复杂外形的高精度计算谱方法精度高外形、边界条件简单简单外形的高精度计算粒子类方法算法简单，可处理复杂外形精度不易提高复杂外形的工程计算Copyright by Li Xinliang课程安排 流体力学基本方程 双曲型方程组及其间断（Riemann）解 差分法 （1）： 数学基础及Fourier分析方法 差分法 （2）： 高精度激波捕捉格式 差分法 （3）： 通量分裂技术 有限体积法（1） 有限体积法（2） 间断有限元方法 代数方程组的求解 不可压方程的数值方法 网格生成技术 湍流与转捩（1） 湍流与转捩（2） 燃烧及化学反应流动初步； 并行计算编程初步 （MPI Part1) 并行计算编程初步 (MPI part2, OpenMP)Copyright by Li Xinliang§ 1.2 流体力学基本方程组1. 基本概念流体连续地充满整个空间连续介质假设流体质点：微观充分大，宏观充分小 举例说明流体密度定义平均密度： 控制体内流动的总质量/控制体体积体积为V的控制体控制体太大，有宏观波动控制体太小，有微观波动微观充分大，宏观充分小控制体内的平均密度随体积变化规律