第15章 计算声场的有限差分法 理论声学 教学课件.ppt

理论声学(2) 第十五章 第十五章 计算声场的有限差分法 理想的模型，比较简单，解析的结果，深入的分析，深刻地揭示物理本质 实际问题常常比较复杂，研究的定量化，近代飞速发展的计算机,数值模拟 计算声学 数值计算不能直接得到一般的规律和物理意义 复杂问题的计算结果的准确性需要认真分析 三者结合 数值计算的方法非常多 纯数值计算方法，如有限差分、有限元 半解析的方法，包括变分法和微扰法 有限差分法，直观 二维问题 三维问题根据二维问题推广得到 差分算法 离散化 导数用采样值的线性组合来近似 泰勒展开 差分格式的精度 差分格式的长度越大，精度越高 实际计算中适当选取差分近似的阶数 交错网格差分算法 中心差分格式求导数时没有用到本身的函数值 为了提高效率，导数的采样位置和函数值的采样位置不同 交错网格差分格式 交错网格差分近似是 交错网格方法的精度和效率更高，稳定性更好 声场问题 一阶偏微分方程 时间二阶差分迭代格式 初始条件和边界条件 反复利用差分方程逐步迭代，求出各个时刻的声场 计算中采用当地的数值，易于处理非均匀问题 柱面波被平面散射的声场 正方形中左上角的线源激发声场 点源在立方体房间中激发的声场的二维模拟 稳态振动和共振频率的计算 线性代数方程组 共振频率 有限差分计算的性能 收敛性，步长趋于零时数字计算结果趋向于对应的微分方程的解 精确度，截断误差的大小，与差分近似的阶数相似 稳定性，各种误差在计算过程中的变化和影响 有限差分算法的稳定性 差分算法是一个迭代过程 迭代表示成矩阵的乘法 模式在这个迭代过程中的稳定性 对平面波的复振幅迭代 每个时间步长只要迭代一次 迭代的结果由 确定 考虑在什么条件下迭代结果是等比级数 特征值问题 为了保证迭代结果不发散，公比的绝对值不能大于1 稳定性要求 两个方向的步长一样 正弦函数换成1 三维问题稳定性条件 步长选择 奈奎斯特定理 实际的步长还要小 空间步长对时间步长的限制 差分计算的物理类比 数字频散 计算带来的声速的误差问题 稳态平面波 如果空间步长非常小，正弦函数用宗量代替，上式趋于微分方程组的结果 相速度 时间步长使相速度随频率增大，相对误差 空间的步长使相速度随频率减小，相对误差 和一个波长里的节点数的平方成反比 随传播方向变化 高阶差分的频散，相速度的相对误差和高次幂成比例 这种频散不是介质的性质，而是计算方法引起的，网格频散，数字频散 脉冲波 脉冲波是频率不同的正弦波组成的 脉冲波在传播中不同频率的成分传播速度不一样，使脉冲波的波形畸变 有限元方法 线性剖分 插值 声波方程 声波方程 边界条件 弱形式的声波方程 检验函数 如果 得到变分形式的声波方程 另一种弱形式 有限元方程 剖分和插值 选取一组检验函数，使弱形式成立，得到近似解 加权余量法 如果检验函数选为各节点的基函数，方法中出现的矩阵常常是对称的，称为伽列金方法。 瞬态振动问题, 声压在时间上离散化 稳态受迫振动问题 高阶有限元算法 多项式，更好地逼近待求的函数 一般方向四边形 谱元法 切比雪夫多项式 有限元法和谱方法结合 快速收敛特性 在一定精度要求下节点数最小 二维谱元法的节点 辐射和散射问题的边界元方法 无限大均匀介质 赫尔姆兹惠更斯积分 * 交错网格差分近似 频散关系 *