二维对流方程在 A、 B、 C 三种差分格式.pdfVIP

一、 实验目的 1、 用数值方法计算二维对流方程在A、B、C 三种差分格式下的解以及用时间分裂法计算C 格式。 2、 计算范围为在−16 ≤x ≤ 16及−16 ≤y ≤ 16，取Δx = Δy = 1，计算至少90 个的时间步长。 二、实验原理及步骤 u u u    0 t x y 1,当4  x  4,4  y  4   u x, y,0   0, 其他 t t 1、编程用迎风C 格式、A 格式在 , =0.5, 1, 2 三种情况下计算。 x y 2、计算范围为在16  x 16,16  y 16 ，取x 1, y 1,t  0, 计算至90 个的时间步长。 t t 3、选择另外两种格式，讨论在 , =0.5, 1, 2 下的差分格式的稳定性。 x y 4、分别