第二章 计算力学的基本知识.docVIP

第二章 计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。数学的发展计算机的不断进步以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了计算流体力学。从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。.2 计算流体力学的定义 计算流体动力学（Computational Fluid Dynamics ,简称CFD）是建立在经典流体力学与数值计算方法基础上的新型独立的学科，通过计算机数值计算和图像显示的方法，在时间和空间上定量描述流场的数值解，从而达到对物理问题研究的目的。它兼有理论性和实践性的双重特点，建立了理论和方法，为现代科学中许多复杂流动和传热问题提供了有效的计算技术。 计算流体动力学（CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。它的基本思想是：把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的几何来代替，通过一定的原则和方式建立起来的关于这些离散点上场变量之间关系的代数方程组，然后代数方程组获得场变量的近似值[5]。 CFD方法和传统的理论分析方法、实验测量方法组成了研究流体流动问题的完整体系。 理论分析方法的优点在于所得结果具有普遍性，各种影响因素清晰可见，是指导实验研究和验证数值计算方法的理论基础，但是它往往要求对计算进行抽象和简化，才可能得出理论解。对于非线性情况，只有少数流动才能给出解析结果。 实验测量方法所得到的实验结果真实可信，它是理论分析和数值方法的基础，其重要性不容低估。然而，实验往往受到模型尺寸、流场流动、人身安全和测量精度的限制，有时可能很难通过试验的方法得到满意的结果。 而CFD方法恰好克服了前面两种方法的弱点，在计算机上实现一个特定的计算，就好像在计算机上做一个物理实验。例如，机翼的绕流，通过计算机并将其结果在屏幕上显示，就可以看到流场的各种细节：如激波的运动、强度，涡的生成与传播，流动的分离、表面的压力分布、受力大小及其随时间的变化等。数值模拟可以形象地再现流动情景，与做实验没有什么区别。 2.1.3 计算流体力学的计算步骤 采用CFD的方法对流体流动进行数值模拟，通常包括如下步骤： （1）建立反映工程问题或物理问题本质的数学模型。具体的说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数学模型的出发点。没有正确完善的数学模型，数值模拟就没有任何意义。流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程，以及这些方程相应的定解条件。 （2）寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法，如有限差分法、有限元法、有限体积法等。这里的计算方法不仅包括微分方程的离散化方法及求解条件，还包括体坐标的建立，边界条件的处理等。这些内容可以说是CFD的核心。 （3）编制程序和进行计算。这部分工作包括计算网格划分、初始条件和边界条件的输入，控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较的复杂，比如Navier-Stokes方程就是一个十分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正是从这个意义上讲，数值模拟又叫数值实验。 （4）显示实验的结果，计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要的意义。 2.1.4 计算流体力学的局限性 虽然CFD具有许多的优点，但是也存在一定的局限性。首先，数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合在计算机上进行计算的离散的数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差；第二，它不像物理模型试验一开始就能给出流动现象并定性的描述，往往需要由原体观测或物理模型试验提供某