计算流体力学常用数值方法简介.pdfVIP

2004年第3期 广 东 造 船 5 计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 — L —-L 一 日JI 吾 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的，这些 基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运 动的控制方程，揭示流体运动的物理规律，研究流体运动的时一空物理特征，这样的学科称 为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科，它涉及流体力学理论、 计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世 纪60年代中后期逐步发展起来的，大致经历了四个发展阶段：无粘性线性、无粘性非线性、 雷诺平均的N—s方程以及完全的N—s方程。随着计算机技术、网络技术、计算方法和后 处理技术的迅速发展，利用计算流体力学解决流动问题的能力越来越高，现在许多复杂的流 动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展，计算流体力学己经成为一种有力的数值实验与设计手段，在许多 工业领域如航天航空 、汽车、船舶等部门解决了大量的工程设计实际问题，其中在航天航空 领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外 形，确定其气动载荷，从而有效地提高了设计效率，减少了风洞试验次数，大大地降低了设计 成本。此外，计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以 及工业中的化学反应等各个领域，显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加，计算流体力学也己 经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线，包括计算机技术、 计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向 二个方向发展：一方面是研究流动非定常稳定性以及湍流流动机理，开展高精度、高分辩率 的计算方法和并行算法等的流动机理与算法研究；另一方面是将计算流体力学直接应用于 模拟各种实际流动，解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种，其数学原理各不相同，但有二点是所有方法都具备的，即 离散化和代数化。总的来说其基本思想是：将原来连续的求解区域划分成网格或单元子区 6 广 东 造 船 2004年第3期 域，在其中设置有限个离散点(称为节点)，将求解区域中的连续函数离散为这些节点上的 函数值；通过某种数学原理，将作为控制方程的偏微分方程转化为联系节点上待求函数值之 间关系的代数方程(离散方程)，求解所建立起来的代表方程以获得求解函数的节点值。不 同的数值方法，其主要区别在于求解区域的离散方式和控制方程的离散方式上。在流体力 学数值方法中，应用比较广泛的是有限差分法、有限元法、边界元法、有限体积法和有限分析 法，现简述如下。 l 有限差分法 这是最早采用的数值方法，它是将求解区域划分为矩形或正交曲线网格，在网格线交点 (即节点)上，将控制方程中的每一个微商用差商来代替，从而将连续函数的微分方程离散 为网格节点上定义的差分方程，每个方程中包含了本节点及其附近一些节点上的待求函数 值，通过求解这些代数方程就可获得所需的数值解。 有限差分法的优点是它建立在经典的数学逼近理论的基础上，容易为人们理解和接受； 有限差分法的主要缺点是对于复杂流体区域的边界形状处理不方便，处理得不好将影响计 算精度。 2 有限元法 有限元法的基本原理是把适定的微分问题的解域进行离散化，将其剖分成相连结又互 不重叠的具有一定规则几何形状的有限个子区域 (如：在二维问题中可以划分为三角形或 四边形；在三维问题中可以划分为四面体或六面体等)，这些子区域称之为单元，单元之间 以节点相联结。函数值被定义在节点上，在单元中选择基函数(又称插值函