有限体积法介绍.docxVIP

有限体积法 1有限体积法基本原理 上一章讲到的有限差分法将数值网格的节点上定义为计算节点， 并在网格节点上对微分 形式的流体基本方程进行离散，用网格节点上的物理量的代数方程作为原 PDE的近似。 在本章所要学习的有限体积法则采用了不同的离散形式。 首先，有限体积法离散的是积 分形式的流体力学基本方程： s V nds S n ds q d ? ⑴ 计算域用数值网格划分成若干小控制体。 和有限差分法不同的是， 有限体积法的网格定 义了控制体的边界，而不是计算节点。有限体积法的计算节点定义在小控制体内部。 一般有 限体积法的计算节点有两种定义方法， 一种是将网格节点定义在控制体的中心， 另一种方法 中，相邻两个控制体的计算节点到公共边界的距离相等。 第一种方法的优点在于用计算节点 第二种方法的好处是在控制体边界上的的值作为控制体上物理量的平均值具有二阶的精度; 第二种方法的好处是在控制体边界上的 中心差分格式具有较高的精度。 QQ Q Q O ■O C O O O O O O O Q 0 积分形式的守恒方程在小控制体和计算域上都是成立的。 为了获得每一个控制体上的代 数方程，面积分和体积分需要用求面积公式来近似。 2面积分的近似 采用结构化网格，在二维情况下，每一个控制体有 4个面，二维情况，每一个控制体有 6个表面。计算节点用大写字母表示， 控制体边界和节点用小写字母表示。 为了保