5.2 Engquist-.Majda吸收边界条件 966 Mur差分格式 近似的数值的吸收边界条件 aU 1aU > x=0边界,一阶近似:x c Ot 0 >展开点:x=△x/2处、t=(n+0.5)△t时刻 >中点差分格式+中点平均近似 "m:m( >Mur的一阶吸收边界条件 )(1)(-U] 7
7 5.2 Engquist-Majda吸收边界条件 Mur差分格式 近似的数值的吸收边界条件 x=0边界, 一阶近似: 展开点:x=x/2处、t=(n+0.5)t 时刻 中点差分格式+中点平均近似 ; ; Mur的一阶吸收边界条件 1 0 U U x c t 1 1 2 2 1 1 , 2 2 1 (1) (0) n n U x n t U U x x 1 1 1 , 1 1 1 1 2 2 2 2 n n U x n t U U c t c t U m U m U m n n n 1 2 1 1 2 ( ) U m U m U m n n n 1 2 1 2 1 1 1 0 1 1 0 n n n n c t x U U U U c t x U j 0, x 1 , 2 U j y U j 1,
5.2 Engquist-.Majdal吸收边界条件 966 Mur差分格式 0(0j+1) 0L+1) 近似的数值的吸收边界条件 a'U 18U caU x=0边界,二阶近似:xtca222 =0 >展开点:辅助网格点(0.5△x,j△y)处、t=(n+0.5)△t时刻 >中点差分格式+中点平均近似 U(j-1) >Mur的二阶吸收边界条件(5.46) =-w+,点+w+24,+肌) 0A-2吧,+w5+w-2+wr) (c△t)2△x 十 ■在角点处并不适用,因为其中要用到的某些网格点的数据无法知道,这些网格点位于网格区域之外 >三维情形的Mu吸收边界条件(5.49),类似推导 8
8 5.2 Engquist-Majda吸收边界条件 Mur差分格式 近似的数值的吸收边界条件 x=0边界, 二阶近似: 展开点:辅助网格点(0.5x,jy) 处、t=(n+0.5)t 时刻 中点差分格式+中点平均近似 Mur的二阶吸收边界条件(5.46) 在角点处并不适用,因为其中要用到的某些网格点的数据无法知道,这些网格点位于网格区域之外 三维情形的Mur吸收边界条件(5.49),类似推导 2 2 2 2 2 1 0 2 U U c U x t c t y 1 1 1 1 0, 1, 1, 0, 0, 1, 2 2 0, 1 0, 0, 1 1, 1 1, 1, 1 2 ( ) 2 2 2( ) ( ) n n n n n n j j j j j j n n n n n n j j j j j j c t x x W W W W W W c t x c t x c t x W W W W W W y c t x