極坐標(biāo)下Helmholtz方程的高階緊致差分方法研究.pdf

1、Helmholtz方程是一類重要的物理方程，如電磁場中的波導(dǎo)問題、噪聲的控制、薄膜振動(dòng)問題等都是由Helmholtz方程控制的.近年來，許多學(xué)者運(yùn)用有限體積法、有限元法、有限差分法等數(shù)值方法在數(shù)值求解Helmholtz方程方面做了大量工作，方法的求解精度可達(dá)四階、六階或更高階.當(dāng)波數(shù)為分段常數(shù)時(shí)，采用有限差分方法求解Helmholtz方程，求解精度常無法達(dá)到預(yù)期精度，方法甚至不可用.對(duì)這一類情況，浸入界面方法能有效的處理。此方法主要是利

2、用跳躍條件將界面一側(cè)的網(wǎng)格點(diǎn)過渡到界面另一側(cè)，對(duì)格式進(jìn)行修正，最終達(dá)到格式的預(yù)期精度。
　　本文利用有限差分方法，結(jié)合浸入界面方法的思想對(duì)極坐標(biāo)系下帶有不連續(xù)波數(shù)的Helmholtz方程進(jìn)行求解.首先，利用文獻(xiàn)中提出的用于求解波數(shù)連續(xù)情況下一維Helmholtz方程的四階緊致差分格式，結(jié)合浸入界面方法和Taylor級(jí)數(shù)展開，對(duì)波數(shù)是分段常數(shù)的一維Helmholtz方程進(jìn)行處理。通過在已建立的格式上添加修正項(xiàng)，利用跳躍條件將界面附近

3、的網(wǎng)格點(diǎn)過渡到界面的另一側(cè)，由此確定修正項(xiàng)的系數(shù)，最終建立波數(shù)不連續(xù)時(shí)一維Helmholtz方程四階緊致差分格式.
　　其次，對(duì)于二維Helmholtz方程，波數(shù)連續(xù)時(shí)將其寫成兩個(gè)Laplacian算子的形式分別進(jìn)行離散，將兩個(gè)算子的四階緊致差分格式聯(lián)合在一起就得到了二維Helmholtz方程的四階緊致差分格式.對(duì)于波數(shù)不連續(xù)的情況，和一維處理方法類似，利用跳躍條件將界面附近的網(wǎng)格點(diǎn)過渡到界面的另一側(cè)，由此確定修正項(xiàng)的系數(shù)，從而可