交通大学课件 计算流体力学.ppt

交通大学课件计算流体力学;目录;第七章非稳态流动问题的有限体积法;7-1非稳态流动问题的守恒方程;利用奥氏公式将对流项和扩散项的体积积分转换为对控制体的表面积分，

并将第一项的时间积分与体积积分的顺序对调，得：;7-2非稳态扩散问题的离散方程;上式左端，令;为计算（7-8）式右端对时间的积分，需要知道积分式中各个变量（主要为TE、TP、TW）随时间的变化规律，但这一变化规律并不知道。通常的处理方法是利用t时刻的温度和(t+Δt)时刻的温度的加权温度作为这一时间段的平均温度，即（以P点以例）：;积分时间域中的三种温度变化规律（以P点以例）;式中：;二、显式计算格式;Δt为时间步长，可知，上式对时间步长作出了相当严格的限制：网格越密,即Δx

越小，上式要求的Δt就越小，计算工作量就越大。故显示格式不适合于一般情况下

的非稳态问题。它只适合简单扩散问题的计算。;若k为常数，并采用均匀网格，则上述条件可写成：;式中：;例7.1一无限大薄板，初始状态为均匀温度200℃，在某一时刻t=0，板东侧面温度

突然降到0℃，西侧面保持绝热。板厚L=2cm，热传导系数k=10W/m/K,

ρc=10×106J/m3/K。

求:(1)采用有限体积法的显式格式，选用一合理的时间步长，计算t=40s，80s，120s

时板的温度分布.

(2)采用显式格式的极限稳定时间重新计算t=40s的过程，并比较显式格式、隐式格

式和解析解的计算结果。

解：本问题的控制微分方程为;将计算区域离散成5个控制体，每个控制体长0.004m，如图。;显式格式的极限稳定时间步长由（7-15）式决定;在（7-14）中，源项S为0。代入数值，得2~4节点的离散方程：;对于边界控制体5，当t0时东侧界面温度TB=0，故对通式（a)进行修改：;解控制体1~5组成的离散方程组，得到第一时间步的结果。将此结果作为第二

时间步中的初始温度，即可算出第二步的结果，如此推进下去，直到收敛。;以上迭代过程可由excel完成;(2)采用显式格式的极限稳定时间重新计算t=40s的过程，并比较显式格式、隐式格式和解析解的计算结果。;若采用全隐格式计算本题（例7.2）仍取Δt=2s，得如下结果;下图全隐格式（例7.2），Δt=8s，t=40s的结果，可知

全???格式稳定性好，计算精度高于显示格式，所以在高

维问题成为推荐算法;五、高维非稳态扩散问题全隐格式;7-3非稳态对流扩散问题的离散方程;二、QUICK格式非稳态对流扩散问题;若本题为无源、稳态问题，则控制微分方程为：;上式除第一项外，采用全隐格式，有：;即：;对于控制体1，由于：;对于控制体2，由于：;对于中间控制体有：;;本题通过matlab编程求解;7-4非稳态压力-速度耦合问题求解过程;本章作业：7-1，7-2，7-3

下周上机作业：7-4