声传播算子计算方法及其改进.pdfVIP

声传播算子计算方法及其改进+ 卢晶徐柏龄 (南京大学声学研究所) AcousticalWaye andIts PropagatorImprovement Lu Jing，XuBoling Institute (The ofAcoustics，NanjingUniversity) l，引言 ． 近年来，随着数值计算技术的发展和计算机运算能力的提高，声波方程的时域解越来越受到人们的关 注[1】，但大多数计算方法由于受到差分方法计算精度的影响，在计算域边界情况比较复杂时，计算的准确 度和效率都受到很大的局限；有限元或边界元方法虽然计算精度较高，但其运算量太大，无法适用于可变 参数声场的研究。量子化学家和物理学家在解决时变薛定谔方程问题时提出的量子波动算子在研究不同量 子系统中瞬态电子波动方程的演化中是非常精确高效的工具。这一方法已经被成功的应用于量子化学和原 o-”量级)。2000年，澳 子物理中【2】，无论波传输的时间多长，其准确度都接近于机器精度(绝对误差在1 太利亚的潘杰教授将该方法引入到声场的时域计算中，提出了声传播算子方法。声传播算子方法[3]将频域 计算的谱方法和时域计算的切比雪夫展式相结合，是一种高效准确的时域计算方法。该方法有一个显著的 优点；只需要知道空间声学参数(包括声速，密度等)的分布，不需要人为构造边界条件，这样便使得计 算与实际情况匹配的更好，并能及时地追踪参数变化对声场的影响。 2．声传播算子的模型及其改进 空气，水以及固体中的声波波包随时间的演化特性由以下偏微分方程决定[4]： a ． ÷e(x，f)=-HdP(x，f) (1) 凹 公式中的x代表空间坐标系，巾是状态矢量，n是空间微分算子。以一维的声传播情况为例，则如果 m=[p，V】2代表一维管中的声压p(而r)和质点速度V(工，f)，那么 0 见e拿 ox H= 土； o p，ox 这里co是空气中的声速，岛是空气密度。 对公式(1)时间积分便可以获得声传播算子表达式 O(x，t)=-e-ftft％’o@，to) (3) 根据这一表达式，任意时刻的状态矢量中(工，f)都可以通过对初始状态中(工，to)加声传播算子P一“‘“‘’得到。 ‘自然科学基金基金号60340420325，博士点基金号20040284038 281 声传播算子需要以有限多项式展开的近似式完成对初始状态的操作，对时变薛定谔方程的研究表明切比雪 夫多项式展开法有其独到的优势[5]。 声传播算子的切比雪夫多项式展开方式如下； e-l’(t-to)=e-RI：I=艺％(R)瓦m) (4) n=O 雪夫多项式。进一步的研究[6]还表明采用修正切比雪夫多项式展开： e妒哪。=∑成(R)g(由