3. 声学基础与波动方程

声学仿真的第一步不是打开 COMSOL 建几何，而是把“声音”翻译成一组可计算的物理量。声压、声速、声阻抗、波动方程、Helmholtz 方程，这五个概念构成了整个声学模块的地基。理解它们之间的推导关系，才能明白 COMSOL 中那些看似零散的设置项——为什么频域研究要选 Helmholtz 接口、为什么时域要选瞬态压力声学、为什么材料里要填密度和声速——究竟在约束什么。

声压：把空气的疏密变化写成一个标量场

声波本质是介质中微小的压强扰动。设介质处于静压 p₀（空气中约 101325 Pa），声波经过时总压强变为 p₀ + p，其中 p 就是声压，单位 Pa。人耳可听阈约 20 μPa，痛阈约 20 Pa，跨越六个数量级，因此工程上常用声压级 SPL = 20·log₁₀(p/p_ref) 表示，参考声压 p_ref = 20 μPa。例如 p = 0.2 Pa 对应 SPL ≈ 80 dB，属于嘈杂街道水平。

在 COMSOL 中，声压是压力声学接口的因变量，默认名 p。所有边界条件、源项、后处理都围绕这个标量场展开。这一点很关键：声压是标量，而结构位移是矢量，两者耦合时正是通过法向分量建立联系。

声速与声阻抗：介质的两个“身份证”

声速 c 描述扰动传播的快慢，由介质可压缩性与惯性共同决定。理想气体中 c = √(γRT)，空气中常温约 343 m/s；水中约 1500 m/s；钢中约 5900 m/s。声阻抗定义为 Z = p/u，其中 u 是质点速度，对平面波有 Z = ρc，ρ 为介质密度。空气的 ρc ≈ 1.21 × 343 ≈ 415 Pa·s/m，水约 1.5 × 10⁶ Pa·s/m，两者相差约 3600 倍。

这个巨大差异直接决定了声学工程的核心难题：空气与水的界面几乎全反射。在 COMSOL 中，声阻抗并不需要单独输入，它由材料节点的密度 ρ 和声速 c 自动组合而成。因此材料设置里填的 ρ 和 c 不是装饰，而是决定反射、透射、共振频率的关键参数。下表给出几个常用介质的典型值，可作为建模时的参考起点。

从波动方程到 Helmholtz 方程

把连续性方程、动量方程和状态方程线性化，消去质点速度，得到齐次波动方程：

∇²p − (1/c²)·∂²p/∂t² = 0

它描述声压随时间和空间的变化。若声源以单一角频率 ω 简谐振动，可令 p(x,t) = p(x)·e^{jωt}，代入后时间导数变为 jω，波动方程退化为 Helmholtz 方程：

∇²p + k²p = 0,   k = ω/c = 2πf/c

其中 k 是波数，单位 rad/m。Helmholtz 方程只含空间变量，把时域问题转成了频域问题。这就是 COMSOL 中“压力声学，频域”接口的数学内核。频域求解一次即得整个频率点的稳态复声压，适合扫频、传递函数、共振分析；时域则用“压力声学，瞬态”接口直接积分波动方程，适合脉冲、非线性、宽频瞬态过程。

COMSOL 中的对应设置

在“压力声学，频域”接口下，材料节点需指定流体密度和声速，软件据此计算波数 k。研究节点选择“频域”，输入频率列表，例如 range(100,50,2000) 表示 100 Hz 到 2000 Hz、步长 50 Hz。边界条件中，“声硬”对应法向速度为零（刚性壁面），“声软”对应声压为零（压力释放），“阻抗”边界则直接输入法向阻抗 Z，用于模拟吸声衬里或辐射负载。

若切换到“压力声学，瞬态”，研究节点需设置时间步长与总时长。时间步长应满足 Courant 条件：c·Δt/Δx ≤ 1，即每步声波传播距离不超过一个网格单元。例如空气中 c = 343 m/s，网格尺寸 0.01 m，则 Δt ≤ 2.9×10⁻⁵ s。COMSOL 可自动估计，但手动检查能避免数值色散。

下面这张图把本节的核心关系串起来：从介质参数 ρ、c 出发，经波动方程与简谐假设，分别通向时域与频域两条求解路径，并标出各自在 COMSOL 中的接口、研究类型和关键设置。

频域与时域的物理含义对照

频域求解的是稳态响应：声源持续以频率 f 振动，系统经过足够长时间后进入周期稳态，此时声压幅值和相位不再随时间变化。COMSOL 频域研究输出的是复数声压，其模为幅值、辐角为相位。这让我们能直接读取某点的传递函数、计算吸声系数、绘制指向性图。时域求解的是瞬态过程：从初始条件出发，逐步推进，能看到波前传播、反射叠加、驻波形成的完整时间历程。

一个常见的教学示例是：空气中 1 kHz 平面波垂直入射到 0.05 m 厚的吸声层，背衬刚性壁。频域扫描 100–2000 Hz，可观察到吸声系数在四分之一波长匹配频率附近出现峰值；若改用瞬态，给一个高斯脉冲，则能看到脉冲进入吸声层、部分反射、部分透射，最终在刚性壁处再次反射回来的时序。两种方法结果一致，但呈现的信息维度不同。

把本节内容落到操作上，可以形成一条清晰的检查链：先确认介质 ρ 和 c 是否与工况一致；再根据关心的物理过程选择频域或瞬态接口；然后检查边界条件是否与实际情况匹配；最后核对频率列表或时间步长是否覆盖目标范围且满足精度要求。这四步做完，波动方程才真正被“翻译”成了 COMSOL 里可求解的模型。