基本方程的数值解.pptVIP

基本方程组的数值求解;一、引言;二、积分区域与微分方程的离散化 ;2．微分方程的离散化;;三、交错网格;;四、差分格式;1）物理上的真实性;2）迭代求解的收敛性;3）解的精确性;2．对流项和扩散项的差分;定常一维的流动和扩散过程，其控制方程为 （8） 在图2所示的控制容积上积分得 （9） 另一方面，方程（8）有精确解，通解为 （10） 代入W-P段两端的边界条件：当 时， ；当 时， 可得该段的解为 其中 为w截面的佩克莱特数。 ;将上式代入式（9）的后一项，可得 （11） 其中 为w截面的参数， 。通解中代入P-E段两端的边界条件，同样??得 （12） 这里 ， 。 ;将式（11）和式（12）代入式（9），并利用连续条件 可得 其中 ， ， （13）;几种对流项差分格式;隐式中心差分格式;（15）式中各系数分别为 ;显式中心差分格式 ;2）迎风差分格式;将精确解给出的系数，以及中心差分和迎风差分的系数 [分别为式（16）和式（19）] 同时画出在图4中 从图中可看出，当时，中心差分格式比迎风差分格式更接近于精确值。但当数增大或减小时，中心差分很快就远离精确值 迎风差分格式在数较小时，精度不如中心差分格式高，但它对精确值的偏离不随数或流动雷诺数改变，适合解高雷诺数的流动。; 图4 差分方程系数的比较 根据图4，可推荐混合差分、指数差分和乘方差分等几种差分格式;3）混合差分格式;4）指数差分格式;5）乘方差分格式;当 时，此格式和混合差式相同 系数可用如下的简洁形式表示：