基于高阶数值方法的超声声场模拟与分析.docVIP

  基于高阶数值方法的超声声场模拟与分析# 马友能1,2，余锦华1,2，汪源源1,2** （1. 复旦大学专用集成电路与系统国家重点实验室，上海 200433； 5 10 15 20 25 30 35 2. 复旦大学电子工程系，上海 200433） 摘要：超声声场计算及模拟在超声的运用中具有重要地位。在利用数值方法进行声场模拟和 分析时，对空间和时间离散化引起的误差将随时间积累。当计算和仿真远场或仿真计算时 间较长时，低阶数值方法将产生较大的数值色散。必须进一步细化空间与时间步长，从而 导致存储空间和计算时间的显著增加。高阶数值方法是减少数值色散的有效途径。本文将 两种高阶数值计算方法（高阶 FDTD 算法和傅立叶伪谱法（FPSTD））用于声场仿真和计 算，详细讨论了其在超声仿真中的实现方法，实现了相应的高阶完全匹配层（PML）。通过 仿真实验，分析和研究了高阶算法的性能。结果显示：高阶数值方法在计算和模拟声场方 面十分有效，尤其适用于远场声场计算和仿真，可以节约运算时间和存储空间。在二维仿 真中发现，PML 的吸收效果十分理想。本文为高阶数值方法在超声仿真和分析中的运用提 供有效指导。 关键词：高阶时域有限差分算法 ; 傅里叶时域伪谱算法; 完全匹配层; 声场模拟与分析 中图分类号：TB553 Simulation and Analysis of Ultrasonic Sound Field with High Order Numerical Methods MA Youneng1,2, YU Jinhua1,2, WANG Yuanyuan1 (1. National Key Laboratory of ASIC and System of Fudan University,Shanghai 200433; 2. Department of Electronic Engineering,Fudan University,Shanghai 200433) Abstract: Simulation and calculation of ultrasonic sound field plays an important role in the utility of ultrasound. When numerical methods are used for the analysis and simulation, numerical dispersion due to the discretization of time and space increases over time. Low order numerical methods, which are subject to severe numerical dispersion, are not feasible to simulate far-field ultrasound or ultrasound propagation of long duration, unless refinement of space and time step is carried out which can result in dramatic increase in storage and calculation time. The higher order numeric method is an effective way to get rid of this problem. In this paper, two kinds of high order numeric methods (high order FDTD and FPSTD) are applied to the analysis and simulation of ultrasound field. Their application in ultrasound simulation are discussed in detail. A series of numerical experiments are carried out .The results show that higher order numerical methods could save storage and calculation time and maintain high accuracy at the same time. In the 2D simulation, PML obtains desired absorbing effect. This article could provide valuable gui