计算流体力学_对流方程_有限差分法_Lax格式_迎风格式_FTCS格式.docxVIP

有限差分法数值求解一维伯格斯方程 作者：潭花林 1.引言 本文利用有限差分法计算了一维伯格斯方程的初边值问题。采用 FTCS格式，并深入讨 论了它的相容性、收敛性与稳定。有限差分法在计算流体力学、数值传热学中都有众多的应 用，而且可以用于高维情形。所有问题都是采用 matlab编程计算。本文只是一个简单的一 维问题的算例。 关键词：计算流体力学，有限差分法，一维对流方程 题目 用计算机求对流方程的初值问题 u u 0 t x 1 x 0 1 u u x 0 18,t 0 2 t 0 x 0 的数值解(由于对流方程的计算结果只依赖与上游， 只需要给出上有的边界条件就可以了) (1)分别用C格式，Lax格式，FTCS格 式在 t 1 t 2两种情况下计算。 x 2 x (2)计算范围为 18 x 18，取x 1,t 0， 计算 80个时间步长。 (3)写出计算报告，内容为 (1)计算课题 (II )计算框图 (III)计算程序 (IV)计算结果，t 0,10,20,40时的， u x图 (V)体会 计算原理 迎风格式 内点采用如下差分格式 t ui ,n 1 ui,n ui,n ui 1,n x 2 i 37,1 n 80 初值为 Ui,i u Xi,0 X 18 i 1 x x 1,1 i 37 边界条件为 u1,n u2, n 稳定性： 差分格式的稳定性：误差方程与差分方程相同 i ,n 1 设误差为i ) i,n x Ikxi n Cne ， i,n 则 i 1,n Ikx Ikxi t Ikxi Ikx cn 岸 Cne cne cne x t ’ Ik x cn 1 1 1 e Cn x 放大因子 G 1 —1 e Ik x x 所以 2 2 t 2 t Gx 1 Gy x x 为使 G 1 应有 t 0 1 对于本问题，初值和边界条件并不影响稳定性和收敛性问题。 根据Lax等价定理，对于适定的处置问题，只要差分格式相容，稳定的就是收敛的。 Lax格式 内点采用如下差分格式 1 t i,n 1 2 i 1,n i 1,n 2 x i 1,n i 1,n 2 i 36,1 n 80 初值为 Ui,1 u Xi,0 x 18 i 1 x x 1,1 i 37 边界条件为 Wn U2,n U36,n U37」 稳定性：差分格式的稳定性:误差方程与差分方程相同i ,n 11,ni 1,nti,n i 1,n2 x设误差为i,n cneIkXi放大因子所以 稳定性： 差分格式的稳定性: 误差方程与差分方程相同 i ,n 1 1,n i 1,n t i,n i 1,n 2 x 设误差为 i,n cne IkXi 放大因子 所以 Gx2 为使G cos k cos k x Gy t 应有 △sin k x cn x -sin k x x 1x对于本问题，初值和边界条件并不影响稳定性和收敛性问题。根据Lax 1 x 对于本问题，初值和边界条件并不影响稳定性和收敛性问题。 根据Lax等价定理，对于适定的处置问题，只要差分格式相容，稳定的就是收敛的。 FTCS格式 内点采用如下差分格式 t i,n 1 i ,n i 1,n 2 x 2 i 36,1 n 80 初值为 U,1 u知0 x 18 i 1 x x 1,1 i 37 i 1,n 边界条件为 U1,n U2,n U36,n U37,n 稳定性： 差分格式的稳定性：误差方程与差分方程相同 i ,n 1i 1,ni 1,n i ,n 1 i 1,n i 1,n 设误差为i,n CneIkXi，则 cn 1 1 .t— v c 1 Olli X 放大因子 G 1 I t .. sin k X X 所以 2 2 t G 1 G 1 —sin k x 1 X 所以该格式是完全不稳定的。 计算结果及其讨论 t 1 —-时，结果是收敛的 x 2 所示，这是由于还t 0时刻，三种格式的处理是相同的，处理结果如 错误！未找到引用源。 所示，这是由于还 00没有开始数值计算。 00 —-Ljh格或 F兀沖氏 OS ■■ 04 02 -15 -1010Ifi -15 -10 10 Ifi 图1 t 10时刻，计算结果如 错误!未找到引用源。 所示，Lax格式和迎风格式的结果几乎是相同 的，而FTCS格式却出现了很大的数值震荡。 i= i= 计亘络卑 图2 t 20时刻，计算结果如 错误!未找到引用源。所示，Lax格式和迎风格式的结果几乎是相同 的，并且结果已经接近稳态。而 FTCS格式却数值震荡更大了，在几百的量级。 匸扯时訂的计 图3 t 40时刻，三种格式的处理是相同的，处理结果如 错误！未找到引用源。 所示，Lax格式和 迎风格式的结果几乎是相同的，并且结果已经接近稳态。而 FTCS格式