求解声波动方程的隐格式有限体积法  

An implicit finite volume method for the acoustic wave equation

作　　者：宣领宽[1] 龚京风[1] 张文平[2] 袁兴军[1] 

机构地区：[1]中国舰船研究设计中心,武汉430064 [2]哈尔滨工程大学动力与能源工程学院,哈尔滨150001

出　　处：《哈尔滨工业大学学报》2016年第7期145-149,共5页Journal of Harbin Institute of Technology

基　　金：国家自然科学基金(51509232)

摘　　要：为研究声传播问题,提出一种声波动方程的隐格式有限体积法,该方法将格点型有限体积法与Newmark格式相结合.模拟平面波的传播过程,对比分析隐格式有限体积法和文献中显格式有限体积法的精度、稳定性及计算消耗等方面的性能.数值结果表明:当λ/Δx≥10时,两种算法均能得到满足精度要求的解;采用无条件稳定的隐格式算法,当满足ω0Δt≤0.3时,预测声压的相对峰值误差<1%;当采用相同时间、空间步长时,隐格式算法精度高于显格式算法;隐格式算法对吸收边界的处理精度高于显格式算法,但对全反射边界的处理精度低于显格式算法;两种算法内存消耗比较接近,显格式算法的CPU耗时较少.An implicit finite volume method （ FVM） with Newmark scheme is proposed for the sound propagation problem, and it is used to simulate the propagation of the plane wave. By comparing to the results from the FVM with that of explicit central difference algorithm, the numerical error, stability and computational consumption are analyzed. It is shown that both methods can acquire accurate numerical solutions when λ/Δx≥10. The relative peak error of the implicit scheme is less than 1% when ω_0Δt≤0.3. With the same time step and spatial step, the numerical result of the implicit scheme agrees better with the exact result. The implicit method treats the absorbing boundary condition more accurate than the explicit method, but it is reverse for the total reflecting boundary condition. The two methods consume similar memory and the explicit method consumes less CPU time.

关 键 词：有限体积法 差分格式 Newmark格式 声波动方程 平面波 

分 类 号：O422[理学—声学]