3西南交通大学计算流体力学.pptVIP

目录 第一章 绪论 第二章 扩散问题的有限体积法 第三章 对流扩散问题的有限体积法 第四章 差分格式问题 第五章 压力--速度耦合问题的有限体积法 第六章 有限体积法离散方程的解法 第七章 非稳态流动问题的有限体积法 第八章 边界条件处理 第三章 对流扩散问题的有限体积法 3-1 一维稳态对流扩散问题的FVM计算格式 3-2 多维稳态对流扩散问题的有限体积法求解 * * 计算流体力学电子教案 西南交通大学力学与工程学院 结构工程系 喻勇 2010-3-14 本章内容 3-1 一维稳态对流扩散问题的FVM计算格式 3-2 多维稳态对流扩散问题的FVM求解 复习：高斯公式(奥氏公式） 或： 通用变量方程 瞬态扩散方程? 稳态扩散方程? 瞬态对流扩散方程? 稳态对流扩散方程? 压力速度耦合方程? 非定常项 对流项 扩散项 源项 由通用变量方程得稳态对流扩散方程为： 将上式按张量运算法则展开得： 由上式得无源项一维稳态对流扩散控制微分方程： 上式中,u为φ在x方向的流动速度。φ为通用变量，可为温度、速度等变量,Г为扩散系数或粘性系数。 此外，流动还应满足连续性方程： 在控制体上对控制微分方程和连续性方程积分，分别得： 再由高斯公式，积分得： 注意：须对控制体界面上的物理量进行处理 对控制体界面上的物理量作近似处理(线性插值)： 将以上各式代入(3-5) 和(3-6)式，并令Aw =Ae=A, F=ρu , D=Г/δx ,得: 对第一式进行整理 ,得: 写成标准离散方程: 式中： 注： 例3.1设某场变量经过对流扩散过程从一维区域的x=0点输运到x=L点，如图。输运控制方程为 其中流体密度ρ =1kg/m3， L=1m,扩散系数Г =0.1kg/(m.s) 。本问题 的理论解为： 求：（1)当流速u=0.1m/s,离散成5个控制体时，φ在区域内的分布; （2)当流速u=2.5m/s,离散成5个控制体时，φ在区域内的分布; （3)当流速u=2.5m/s,离散成20个控制体时，φ在区域内的分布. 解析解的推导，以及excel和matlab的图形表达（fplot） fplot(1-(exp(10*2.5*x)-1)/(exp(10*2.5)-1),[0,1]) fplot([1-(exp(10*0.1*x)-1)/(exp(10*0.1)-1),1-(exp(10*2.5*x)-1)/(exp(10*2.5)-1)],[0,1]) 解：本题中控制微分方程在控制体上的积分方程已在面前推出， 即（3-5）式。对于中间控制体，近似处理后（3-5）式可写成离 散方程(3-9)式的标准形式。对于边界控制体，需从（3-5）式出发推 导出相应的离散方程。 式中： 提示：F=ρu , D=Г/δx 下面求最西端处控制体的离散方程，积分方程为(3-5)积 由于： 代入（3-5）式得： (a)式即为最西端控制体的离散方程（未经整理） 同理可得最东端处控制体的离散方程： 由于： 代入（3-5）式得： (b)式即为最东端控制体的离散方程（未经整理） 求：（1)当流速u=0.1m/s,离散成5个控制体时，φ在区域内的分布; （2)当流速u=2.5m/s,离散成5个控制体时，φ在区域内的分布; （3)当流速u=2.5m/s,离散成20个控制体时，φ在区域内的分布. 式中： 对中间控制体： 西： 东： 由以上三式可得到全部离散方程。 (1)当流速u=0.1m/s,离散成5个控制体时，φ在区域内的分布 （2)当流速u=2.5m/s,离散成5个控制体时，φ在区域内的分布 数值计算结果很不理想！ （3)当流速u=2.5m/s,离散成20个控制体时，φ在区域内的分布 计算精度与F/D的值有关，二者关系将在第四章介绍。 由通用变量方程得稳态对流扩散方程为： 将上式按张量运算法则展开得： 由上式得二维稳态对流扩散控制微分方程： 上式中,u、v为φ在x、y方向的流动速度。φ为通用变量，Г为扩散系数或粘性系数。 一、二维对流扩散问题计算格式 *