计算流体力学读书报告.pdf

《计算流体力学》读书笔记 计算流体动力学(Computational Fluid Dynamics ，简称CFD)是通过计算机数 值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分 析。其基本思想为：把原来在时间域及空间域上连续的物理量的场，如速度场和 压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方 式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组 获得场变量的近似值。 CFD 可以看作是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒 方程)控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问 题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等) 的分布，以 及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还 可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。 此外，与CAD 联合，还可进行结构优化设计等。 1. 计算流体力学的目的、内容和特点 计算流体力学的目的是用数值方法求解基本方程（连续性方程、N-S 方程、 Euler 方程），内容是其本身的理论方法以及应用。 计算流体动学的特点： ① 流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和 边界条件复杂，很难求得解析解，而用CFD 方法则有可能找出满足工程 需要的数值解。 ② 可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方 程中各项有效性和敏感性试验，从而进行方案比较。 ③ 它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给 出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实 条件和实验中只能接近而无法达到的理想条件。 ④ 数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、 适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提 供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的 计算误差。 ⑤ 它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要 由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模 型进行验证。 ⑥ 程序的编制及资料的收集、整理与正确利用，在很大程度上依赖于经验 与技巧。 ⑦ 因数值处理方法等原因有可能导致计算结果的不真实，例如产生数值粘 性和频散等伪物理效应。 ⑧ CFD 因涉及大量数值计算，因此，常需要较高的计算机软硬件配置。 理论分析成本最低、结果最理想、影响因素表达清楚；缺点：局限与非常简 单的问题。数值方法成本较低：数值实验、适用范围宽；缺点：可靠性差，表达 困难。实验测量可靠、成本高。将三种方法有机结合，互为补充，必然会取得相 得益彰的效果。 经过四十多年的发展，CFD 出现了多种数值解法。这些方法之间的主要区 别在于对控制方程的离散方式。根据离散的原理不同，CFD 大体上可分为三个 分支：有限差分法(Finite Different Method，FDM) ； 有限元法(Finite EIement Method ，FEM) ；有限体积法(Finite Volume Method，FVM) 。 有限差分法是应用最早、最经典的CFD 方法，它将求解域划分为差分网格， 用有限个网格节点代替连续的求解域，然后将偏微分方程的导数用差商代替，推 导出含有离散点上有限个未知数的差分方程组。求出差分万程组的解，就是微分 方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题的近似数值 解法。这种方法发展较早，比较成熟，较多地用于求解双曲型和抛物型问题。在 此基础上发展起来的方法有PIC(Particle-in-cell)法、MAC(Marker-and-Cell)法，以 及南美籍华人学者陈景广提出的有限分析法(Finite Analytic Method)等。 有限元法是20 世纪80 年代开始应用的—种数值解法，它吸收了有限差分法 中离散处理的内核，又采用了变分计算中选择逼近函数对区域进行积分的合理方 法。有限元法因求解速度较有限差分法和有限体积法慢，因此应用不是特别广泛。 在有限元法的基础上，英国C A ．BBrebbia 等提出了边界元