计算流体力学电子教案ppt课件.ppt

《计算流体力学电子教案ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、计算流体力学电子教案目录第一章绪论第二章扩散问题的有限体积法第三章对流扩散问题的有限体积法第四章差分格式问题第五章压力--速度耦合问题的有限体积法第六章有限体积法离散方程的解法第七章非稳态流动问题的有限体积法第八章边界条件处理第二章扩散问题的有限体积法2-1一维稳态扩散问题的FVM计算格式2-2多维稳态扩散问题的FVM求解预备知识：高斯公式(奥氏公式）或：写出一维条件下的奥氏公式通用变量方程瞬态扩散方程稳态扩散方程瞬态对流扩散方程稳态对流扩散方程压力速度耦合方程非定常项对流项扩散项源项2-1一维稳态扩散问题的FVM计算

2、格式由通用变量方程得稳态扩散方程为：2-1-1一维稳态扩散方程将上式按张量运算法则展开得：由上式得一维条件下的稳态扩散方程：上式中,为通用变量，可为温度、速度等变量；为扩散系数或粘性系数，S为源项。2-1-2求解一维稳态扩散问题的步骤第一步生成离散网格第二步由控制方程(积分形式)形成离散方程组第三步求解方程组第一步：生成离散网格控制体的划分（先划分控制体后定节点,节点在控制体中心）相关的尺寸定义（约定：大写字母代表节点，小写字母代表边界。）第二步：由控制方程(积分形式)形成离散方程组一维稳态扩散控制方程为：将此控制方程在某控制体

3、上积分:则由奥氏公式或高斯散度定理有：式中，控制体的体积为ΔV，全部表面积为A，源项在控制体中的平均值为S上式有明确的物理意义：场变量的净增扩散量（即自西侧界面流入的扩散流量减去东侧界面流出的扩散流量）等于源项产生的扩散流量。积分方程中的下标e、w表示控制体的界面（不是节点处），意味着我们需要知道扩散系数___和场变量的梯度在控制体东西边界上的值。这些值可由节点处的值插值得到。若采用线性插值(近似处理)，对于均匀网格有：同理，有：于是，通过界面的扩散流量为接下来处理源项，源项可能为常数，也可能为场变量的函数，对其进行线性化处理，

4、得：将以上三式代入积分后的控制方程（即下式）中将上式按场变量的节点值进行整理，得：令得离散方程：对于每一个节点（控制体）都可建立一个离散方程，所有节点的离散方程构成一个方程组。由上式形成的方程组是三元一次的线性方程组，该方程的特点是具有三条对角线，故称为三对角线性方程。目前可暂用matlab中Ab语句求解（高斯消元法）。下面用两个例题说明有限体积法如何求一维稳态扩散问题。第三步：解方程组例2.1用有限体积法求解无热源一维稳态导热问题图示绝热棒长0.5m，截面积A=10-2m2，左右端温度保持为TA=100ºC，TB=500

5、ºC。棒材料导热系数k=1000W/(m·K)。求绝热棒在稳定状态下的温度分布。解：本问题的控制微分方程为可将此式与（2-1）式比较,可采用三步求解方法。第一步：生成离散网格（先控制体后节点）,生成5个单元（本问题有解析解）第二步：构造离散方程对求解域中的2、3、4节点应用离散方程(2-8)因故有：式中：由此得到3个方程对于2号控制体对于3号控制体对于4号控制体5个未知数，3个方程。可见不引入边界条件是没法求解的对求解域中的边界节点1、5的离散方程需作特殊处理。方法仍然是对微分方程在边界控制体内积分。微分方程为：上式在左边界控制

6、体上积分，得：即在上述过程中有一假定：认为A点的温度梯度dT/dx与A点和1点的温度线性相关将（2-12）式按节点温度整理得：将上式与（2-8）式对照可知，边界条件可以转化成源项进入控制容积积分方程。即左边界的离散方程可以写成：式中同理可对右边界控制体进行处理，得式中式中根据以上过程可以得到左右边界控制体的离散方程：右端控制体左端控制体TA=100ºC，TB=500ºC第三步解线性方程组本问题的解析解为：T=800x+100例2.2用有限体积法求解有内热源一维稳态导热问题图示厚度为L=2cm的无限大平板，导热系数k=0.5W

7、/(m·K)，板内有均匀内热源q=1000kW/m3,表面温度A、B分别保持为TA=100ºC，TB=200ºC。求板内x向的温度分布。解：由于板在y、z方向为无限大，因此可作为一维问题处理，即只考虑x方向。相对于无源问题，控制方程中增加了源项。即第一步：生成离散网格（先控制体后节点）,生成5个单元第二步：构造离散方程方法一：可以直接套用公式(2-8)，但边界节点需特殊处理。式中注：方法二：也可通过对控制方程的积分推导出离散方程，同例2.1的过程，以下用方法求离散方程。由于有厚度L=2cm，导热系数k=0.5W/(m·K)，板内

8、均匀内热源q=1000kW/m3,表面温度A、B分别保持为TA=100ºC，TB=200ºC。求板内x向的温度分布。由得控制体2、3、4、1、5的离散方程为（打一4个字母的英文单词）2-2多维稳态扩散问题的FVM求解2-2-1二维稳态扩散问题的有限体积法第一步：