5.2 Engquist-Majdal吸收边界条件 966 Trefethen-Halpern近似展开 >算子方程中的根式√1一s2用泰勒级数展开是近似方案之一 V1-s2也可以用有理函数逼近(Trefethen-Halpern提出):、 1-s2=P( q,(s) (m,n)型有理逼近:p,m(s)和qn(S)分别是区间[-1,1]上的m阶和n阶多项式 >(2,0)型有理逼近 D V1-32=p,+p,52;s D,c:E=D,-2- LU=0→ aU Poau oU Oxot c0-Pac =0 02 >po和p2的值由所选的插值方法确定,在[-1,1]上最佳地逼近V1-s2 Pade'逼近、最小二乘逼近、Chebyshevi逼近、 >二阶泰勒级数近似:p0=+1、p2=-0.5的(2,0)型逼近,在x=0边界的正入射方向准确地吸收外向波 >改变系数po和p2,将改变准确吸收外向波的角度 99 5.2 Engquist-Majda吸收边界条件 Trefethen-Halpern近似展开 算子方程中的根式 𝟏 − 𝒔 𝟐 用泰勒级数展开是近似方案之一 𝟏 − 𝒔 𝟐 也可以用有理函数逼近(Trefethen-Halpern提出): (m, n) 型有理逼近: 和 分别是区间[-1, 1]上的m阶和n阶多项式 (2, 0)型有理逼近 ; ; p0和p2的值由所选的插值方法确定,在[-1, 1]上最佳地逼近 𝟏 − 𝒔 𝟐 Pade’逼近、最小二乘逼近、Chebyshev逼近、…… 二阶泰勒级数近似:p0=+1、p2=-0.5的(2,0)型逼近,在x=0边界的正入射方向准确地吸收外向波 改变系数p0和p2,将改变准确吸收外向波的角度 1 2 s r s p s q s m n ( ) ( ) ( ) pm (s) q s n ( ) 1 2 0 2 2 s p p s y t D s D c 2 1 t x D L D s c L U 0 2 2 2 0 2 2 2 0 U U U p p c x t c t y